Showing 1–50 of 57 results for author: Spada, F

Revisiting the trajectory of the interstellar object 'Oumuamua: preference for a radially directed non-gravitational acceleration?

Authors: Federico Spada

Abstract: I present a re-analysis of the available observational constraints on the trajectory of 'Oumuamua, the first confirmed interstellar object discovered in the solar system. 'Oumuamua passed through the inner solar system on a hyperbolic (i.e., unbound) trajectory. Its discovery occurred after perihelion passage, and near the time of its closest approach to Earth. After being observable for approxima… ▽ More I present a re-analysis of the available observational constraints on the trajectory of 'Oumuamua, the first confirmed interstellar object discovered in the solar system. 'Oumuamua passed through the inner solar system on a hyperbolic (i.e., unbound) trajectory. Its discovery occurred after perihelion passage, and near the time of its closest approach to Earth. After being observable for approximately four months, the object became too faint and was lost at a heliocentric distance of around 3 au. Intriguingly, analysis of the trajectory of 'Oumuamua revealed that a dynamical model including only gravitational accelerations does not provide a satisfactory fit of the data, and a non-gravitational term must be included. The detected non-gravitational acceleration is compatible with either solar radiation pressure or recoil due to outgassing. It has, however, proved challenging to reconcile either interpretation with the existing quantitative models of such effects without postulating unusual physical properties for 'Oumuamua (such as extremely low density and/or unusual geometry, non-standard chemistry). My analysis independently confirms the detection of the non-gravitational acceleration. After comparing several possible parametrizations for this effects, I find a strong preference for a radially directed non-gravitational acceleration, pointing away from the Sun, and a moderate preference for a power-law scaling with the heliocentric distance, with an exponent between 1 and 2. These results provide valuable constraints on the physical mechanism behind the effect; a conclusive identification, however, is probably not possible on the basis of dynamical arguments alone. △ Less

Submitted 14 April, 2023; originally announced April 2023.

Comments: MATLAB code will be shared upon reasonable request to the author. Comments are welcome!

Asteroseismology of evolved stars to constrain the internal transport of angular momentum. VI. Testing a parametric formulation for the azimuthal magneto-rotational instability

Authors: F. D. Moyano, P. Eggenberger, B. Mosser, F. Spada

Abstract: Asteroseismic measurements of the internal rotation rate in evolved stars pointed out to a lack of angular momentum (AM) transport in stellar evolution models. Several physical processes in addition to hydrodynamical ones were proposed as candidates for the missing mechanism. Nonetheless, no current candidate can satisfy all the constraints provided by asteroseismology. We revisit the role of a ca… ▽ More Asteroseismic measurements of the internal rotation rate in evolved stars pointed out to a lack of angular momentum (AM) transport in stellar evolution models. Several physical processes in addition to hydrodynamical ones were proposed as candidates for the missing mechanism. Nonetheless, no current candidate can satisfy all the constraints provided by asteroseismology. We revisit the role of a candidate process whose efficiency scales with the contrast between the rotation rate of the core and the surface which was proposed to be related to the azimuthal magneto-rotational instability (AMRI) by Spada et al. We compute stellar evolution models of low- and intermediate-mass stars with the parametric formulation of AM transport proposed by Spada et al. until the end of the core-helium burning for low- and intermediate-mass stars and compare our results to the latest asteroseismic constraints available in the post main sequence phase. Both hydrogen-shell burning stars in the red giant branch and core-helium burning stars of low- and intermediate-mass in the mass range $1 M_{\odot} \lesssim M \lesssim 2.5 M_{\odot}$ can be simultaneously reproduced by this kind of parametrisation. Given current constraints from asteroseismology, the core rotation rate of post-main sequence stars seems to be well explained by a process whose efficiency is regulated by the internal degree of differential rotation in radiative zones. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 10 pages, 10 figures, 1 appendix

Journal ref: A&A 673, A110 (2023)

KIC 7955301: a hierarchical triple system with eclipse timing variations and an oscillating red giant

Authors: P. Gaulme, T. Borkovits, T. Appourchaux, K. Pavlovski, F. Spada, C. Gehan, J. Ong, A. Miglio, A. Tkachenko, B. Mosser, M. Vrard, M. Benbakoura, S. D. Chojnowski, J. Perkins, A. Hedlund, J. Jackiewicz

Abstract: KIC 7955301 is a hierarchical triple system with eclipse timing and depth variations discovered by the Kepler mission. It is composed of a non-eclipsing primary star at the bottom of the red giant branch on a 209-day orbit with a K/G-type main-sequence inner eclipsing binary, orbiting in 15.3 days. This system was noted for the large amplitude of its eclipse timing variations (4 hours), and the cl… ▽ More KIC 7955301 is a hierarchical triple system with eclipse timing and depth variations discovered by the Kepler mission. It is composed of a non-eclipsing primary star at the bottom of the red giant branch on a 209-day orbit with a K/G-type main-sequence inner eclipsing binary, orbiting in 15.3 days. This system was noted for the large amplitude of its eclipse timing variations (4 hours), and the clear solar-like oscillations of the red-giant component, including p-modes of degree up to l=3 and mixed l=1 modes. The system is a single-lined spectroscopic triple. We perform a dynamical model by combining the Kepler photometric data, eclipse timing variations, and radial-velocity data obtained at Apache Point (ARCES) and Haute Provence (SOPHIE) observatories. The dynamical mass of the red-giant is determined with a 2% precision at 1.30 (+0.03,-0.02) solar mass. We perform asteroseismic modeling based on the global seismic parameters and on the individual frequencies. Both methods lead to a mass of the red giant that matches the dynamical mass within the uncertainties. Asteroseismology also reveals the rotation rate of the core (15 days), the envelope (150 days), and the inclination (75 deg) of the red giant. Three different approaches lead to an age between 3.3 and 5.8 Gyr, which highlights the difficulty of determining stellar ages despite the exceptional wealth of available information. On short timescales, the inner binary exhibits eclipses with varying depths during a 7.3-year long interval, and no eclipses during the consecutive 11.9 years. This is why Kepler could detect its eclipses, TESS will not, and the future ESA PLATO mission should. Over the long term, the system owes its evolution to the evolution of its individual components. It could end its current smooth evolution by merging by the end of the red giant or the asymptotic giant branch of the primary star. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Accepted in A&A, 25 pages, 17 figures, 7 tables

Journal ref: A&A 668, A173 (2022)

Stellar evolution models with entropy-calibrated mixing-length parameter: application to red giants

Authors: Federico Spada, Pierre Demarque, Friedrich Kupka

Abstract: We present evolutionary models for solar-like stars with an improved treatment of convection that results in a more accurate estimate of the radius and effective temperature. This is achieved by improving the calibration of the mixing-length parameter, which sets the length scale in the 1D convection model implemented in the stellar evolution code. Our calibration relies on the results of 2D and 3… ▽ More We present evolutionary models for solar-like stars with an improved treatment of convection that results in a more accurate estimate of the radius and effective temperature. This is achieved by improving the calibration of the mixing-length parameter, which sets the length scale in the 1D convection model implemented in the stellar evolution code. Our calibration relies on the results of 2D and 3D radiation hydrodynamics simulations of convection to specify the value of the adiabatic specific entropy at the bottom of the convective envelope in stars as a function of their effective temperature, surface gravity and metallicity. For the first time, this calibration is fully integrated within the flow of a stellar evolution code, with the mixing-length parameter being continuously updated at run-time. This approach replaces the more common, but questionable, procedure of calibrating the length scale parameter on the Sun, and then applying the solar-calibrated value in modeling other stars, regardless of their mass, composition and evolutionary status. The internal consistency of our current implementation makes it suitable for application to evolved stars, in particular to red giants. We show that the entropy calibrated models yield a revised position of the red giant branch that is in better agreement with observational constraints than that of standard models. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: 12 pages, 10 figures; accepted for publication in MNRAS

Spectroscopic and seismic analysis of red giants in eclipsing binaries discovered by Kepler

Authors: M. Benbakoura, P. Gaulme, J. McKeever, S. Sekaran, P. G. Beck, F. Spada, J. Jackiewicz, S. Mathis, S. Mathur, A. Tkachenko, R. A. García

Abstract: Eclipsing binaries (EBs) are unique benchmarks for stellar evolution. On the one hand, detached EBs hosting at least one star with detectable solar-like oscillations constitute ideal test objects to calibrate asteroseismic measurements. On the other hand, the oscillations and surface activity of stars that belong to EBs offer unique information about the evolution of binary systems. This paper bui… ▽ More Eclipsing binaries (EBs) are unique benchmarks for stellar evolution. On the one hand, detached EBs hosting at least one star with detectable solar-like oscillations constitute ideal test objects to calibrate asteroseismic measurements. On the other hand, the oscillations and surface activity of stars that belong to EBs offer unique information about the evolution of binary systems. This paper builds upon previous works dedicated to red giant stars (RG) in EBs -- 20 known systems so far -- discovered by the NASA Kepler mission. Here we report the discovery of 16 RGs in EBs also from the Kepler data. This new sample includes three SB2-EBs with oscillations and six close systems where the RG display a clear surface activity and complete oscillation suppression. Based on dedicated high-resolution spectroscopic observations (Apache Point Observatory, Observatoire de Haute Provence), we focus on three main aspects. From the extended sample of 14 SB2-EBs, we first confirm that the simple application of the asteroseismic scaling relations to RGs overestimates masses and radii of RGs, by about 15% and 5%. This bias can be reduced by employing either new asteroseismic reference values for RGs, or model-based corrections of the asteroseismic parameters. Secondly, we confirm that close binarity leads to a high level of photometric modulation (up to 10%), and a suppression of solar-like oscillations. In particular, we show that it reduces the lifetime of radial modes by a factor of up to 10. Thirdly, we use our 16 new systems to complement previous observational studies that aimed at constraining tidal dissipation in interacting binaries. In particular, we identify systems with circular orbits despite relatively young ages, which suggests exploring complementary tidal dissipation mechanisms in the future. Finally, we report the measurements of mass, radius, and age of three M-dwarf companion stars. △ Less

Submitted 26 January, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 36 pages, 17 figures, radial velocity tables, Accepted in A&A

Journal ref: A&A 648, A113 (2021)

Tunneling Enhancement of the Gas-Phase CH + CO2 Reaction at Low Temperature

Authors: Dianailys Nunez-Reyes, Kevin M. Hickson, Jean-Christophe Loison, Rene F. K. Spada, Rafael M. Vichietti, Francisco B. C. Machado, Roberto L. A. Haiduke

Abstract: The rates of numerous activated reactions between neutral species increase at low temperatures through quantum mechanical tunneling of light hydrogen atoms. Although tunneling processes involving molecules or heavy atoms are well known in the condensed phase, analogous gas-phase processes have never been demonstrated experimentally. Here, we studied the activated CH + CO2 -> HCO + CO reaction in a… ▽ More The rates of numerous activated reactions between neutral species increase at low temperatures through quantum mechanical tunneling of light hydrogen atoms. Although tunneling processes involving molecules or heavy atoms are well known in the condensed phase, analogous gas-phase processes have never been demonstrated experimentally. Here, we studied the activated CH + CO2 -> HCO + CO reaction in a supersonic flow reactor, measuring rate constants that increase rapidly below 100 K. Mechanistically, tunneling is shown to occur by CH insertion into the C-O bond, with rate calculations accurately reproducing the experimental values. To exclude the possibility of H-atom tunneling, CD was used in additional experiments and calculations. Surprisingly, the equivalent CD + CO2 reaction accelerates at low temperature as zero point energy effects remove the barrier to product formation. In conclusion, heavy-particle tunneling effects might be responsible for the observed reactivity increase at lower temperatures for the CH + CO2 reaction, while the equivalent effect for the CD + CO2 reaction results instead from a submerged barrier with respect to reactants. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: Accepted for publication in the Journal of Physical Chemistry A. 30 pages including supplementary materials

Seismic evidence for near solid-body rotation in two Kepler subgiants and implications for angular momentum transport

Authors: S. Deheuvels, J. Ballot, P. Eggenberger, F. Spada, A. Noll, J. W. den Hartogh

Abstract: Asteroseismic measurements of the internal rotation of subgiants and red giants all show the need for invoking a more efficient transport of angular momentum than theoretically predicted. Constraints on the core rotation rate are available starting from the base of the red giant branch (RGB) and we are still lacking information on the internal rotation of less evolved subgiants. We identified two… ▽ More Asteroseismic measurements of the internal rotation of subgiants and red giants all show the need for invoking a more efficient transport of angular momentum than theoretically predicted. Constraints on the core rotation rate are available starting from the base of the red giant branch (RGB) and we are still lacking information on the internal rotation of less evolved subgiants. We identified two young Kepler subgiants, KIC8524425 and KIC5955122, whose mixed modes are clearly split by rotation. Using the full Kepler data set, we extracted the mode frequencies and rotational splittings for the two stars using a Bayesian approach. We then performed a detailed seismic modeling of both targets and used the rotational kernels to invert their internal rotation profiles. We found that both stars are rotating nearly as solid bodies, with core-envelope contrasts of $Ω_{\rm g}/Ω_{\rm p}=0.68\pm0.47$ for KIC8524425 and $0.72\pm0.37$ for KIC5955122. This result shows that the internal transport of angular momentum has to occur faster than the timescale at which differential rotation is forced in these stars (between 300 Myr and 600 Myr). By modeling the additional transport of angular momentum as a diffusive process with a constant viscosity $ν_{\rm add}$, we found that values of $ν_{\rm add}>5\times10^4$~cm$^2$.s$^{-1}$ are required to account for the internal rotation of KIC8524425, and $ν_{\rm add}>1.5\times10^5$~cm$^2$.s$^{-1}$ for KIC5955122. These values are lower than or comparable to the efficiency of the core-envelope coupling during the main sequence, as given by the surface rotation of stars in open clusters. On the other hand, they are higher than the viscosity needed to reproduce the rotation of subgiants near the base of the RGB. Our results yield further evidence that the efficiency of the internal redistribution of angular momentum decreases during the subgiant phase. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: Accepted in A&A, 16 pages, 11 figures

Journal ref: A&A 641, A117 (2020)

A Knee-Point in the Rotation-Activity Scaling of Late-type Stars with a Connection to Dynamo Transitions

Authors: Jyri J. Lehtinen, Maarit J. Käpylä, Nigul Olspert, Federico Spada

Abstract: The magnetic activity of late-type stars is correlated with their rotation rates. Up to a certain limit, stars with smaller Rossby numbers, defined as the rotation period divided by the convective turnover time, have higher activity. A more detailed look at this rotation-activity relation reveals that, rather than being a simple power law relation, the activity scaling has a shallower slope for th… ▽ More The magnetic activity of late-type stars is correlated with their rotation rates. Up to a certain limit, stars with smaller Rossby numbers, defined as the rotation period divided by the convective turnover time, have higher activity. A more detailed look at this rotation-activity relation reveals that, rather than being a simple power law relation, the activity scaling has a shallower slope for the low-Rossby stars than for the high-Rossby ones. We find that, for the chromospheric CaII H&K activity, this scaling relation is well modelled by a broken two-piece power law. Furthermore, the knee-point of the relation coincides with the axisymmetry to non-axisymmetry transition seen in both the spot activity and surface magnetic field configuration of active stars. We interpret this knee-point as a dynamo transition between dominating axi- and non-axisymmetric dynamo regimes with a different dependence on rotation and discuss this hypothesis in the light of current numerical dynamo models. △ Less

Submitted 15 December, 2021; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 10 pages, 5 figures, 2 tables

Report number: NORDITA-2020-052

Journal ref: ApJ 910 110 (2021)

Active red giants: close binaries versus single rapid rotators

Authors: Patrick Gaulme, Jason Jackiewicz, Federico Spada, Drew Chojnowski, Benoit Mosser, Jean McKeever, Anne Hedlund, Mathieu Vrard, Mansour Benbakoura, Cilia Damiani

Abstract: The objective of this work is to determine what fraction of red-giant (RG) stars shows photometric rotational modulation, and understand its origin. One of the underlying questions is the role of close binarity in this population, standing upon the fact that RGs in short-period binary systems (<150 days or so) have been observed to display strong rotational modulation. We select a sample of about… ▽ More The objective of this work is to determine what fraction of red-giant (RG) stars shows photometric rotational modulation, and understand its origin. One of the underlying questions is the role of close binarity in this population, standing upon the fact that RGs in short-period binary systems (<150 days or so) have been observed to display strong rotational modulation. We select a sample of about 4500 relatively bright RGs observed by Kepler, and show that 370 of them (8%) display rotational modulation. Almost all have oscillation amplitudes below the median of the sample, while 30 of them are not oscillating at all. Of the 85 of these RGs with rotational modulation chosen for follow-up radial-velocity observation and analysis, 34 show clear evidence of spectroscopic binarity. Surprisingly, 26 of the 30 non-oscillators are in this group of binaries. To the contrary, about 85% of the active RGs with detectable oscillations are not part of close binaries. With the help of stellar masses and evolutionary states computed from the oscillation properties, it appears that low-mass red-giant branch stars tend to be magnetically inactive, while intermediate-mass ones tend to be highly active. The opposite trends are true for helium-core burning (red clump) stars, whereby the lower-mass clump stars are comparatively more active and the higher-mass ones less so. In other words, we find that low-mass red-giant branch stars gain angular momentum as they evolve to clump stars, while higher-mass ones lose angular momentum. The trend observed with low-mass stars leads to possible scenarios of planet engulfment or other merging events during the shell-burning phase. Regarding intermediate-mass stars, the rotation periods are long with respect to theoretical expectations reported in the literature, which reinforces the existence of an unidentified sink of angular momentum after the main sequence. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: Accepted in Astronomy & Astrophysics (18 pages, 18 figures, 3 tables)

Journal ref: A&A 639, A63 (2020)

Common dynamo scaling in slowly rotating young and evolved stars

Authors: Jyri J. Lehtinen, Federico Spada, Maarit J. Käpylä, Nigul Olspert, Petri J. Käpylä

Abstract: One interpretation of the activity and magnetism of late-type stars is that these both intensify with decreasing Rossby number up to a saturation level, suggesting that stellar dynamos depend on both rotation and convective turbulence. Some studies have claimed, however, that rotation alone suffices to parametrise this scaling adequately. Here, we tackle the question of the relevance of turbulence… ▽ More One interpretation of the activity and magnetism of late-type stars is that these both intensify with decreasing Rossby number up to a saturation level, suggesting that stellar dynamos depend on both rotation and convective turbulence. Some studies have claimed, however, that rotation alone suffices to parametrise this scaling adequately. Here, we tackle the question of the relevance of turbulence to stellar dynamos by including evolved, post main sequence stars in the analysis of the rotation-activity relation. These stars rotate very slowly compared with main sequence stars, but exhibit similar activity levels. We show that the two evolutionary stages fall together in the rotation-activity diagram and form a single sequence in the unsaturated regime in relation only to Rossby numbers derived from stellar models, confirming earlier preliminary results that relied on a more simplistic parametrisation of the convective turnover time. This mirrors recent results of fully convective M dwarfs, which likewise fall on the same rotation-activity sequence as partially convective solar-type stars. Our results demonstrate that turbulence plays a crucial role in driving stellar dynamos and suggest that there is a common turbulence-related dynamo mechanism explaining the magnetic activity of all late-type stars. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: This is the initial submission version of the paper "Common dynamo scaling in slowly rotating young and evolved stars" (Lehtinen et al. 2020). The final version accepted in Nature Astronomy can be found at https://www.nature.com/articles/s41550-020-1039-x

The Aarhus red giants challenge II. Stellar oscillations in the red giant branch phase

Authors: J. Christensen-Dalsgaard, V. Silva Aguirre, S. Cassisi, M. Miller Bertolami, A. Serenelli, D. Stello, A. Weiss, G. Angelou, C. Jiang, Y. Lebreton, F. Spada, E. P. Bellinger, S. Deheuvels, R. M. Ouazzani, A. Pietrinferni, J. R. Mosumgaard, R. H. D. Townsend, T. Battich, D. Bossini, T. Constantino, P. Eggenberger, S. Hekker, A. Mazumdar, A. Miglio, K. B. Nielsen , et al. (1 additional authors not shown)

Abstract: Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated… ▽ More Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated using different stellar evolution codes, we investigate the extent to which the differences in the model calculation affect the model oscillation frequencies. Methods. For each of the models, which cover four different masses and different evolution stages on the red-giant branch, we computed full sets of low-degree oscillation frequencies using a single pulsation code and, from these frequencies, typical asteroseismic diagnostics. In addition, we carried out preliminary analyses to relate differences in the oscillation properties to the corresponding model differences. Results. In general, the differences in asteroseismic properties between the different models greatly exceed the observational precision of these properties, in particular for the nonradial modes whose mixed acoustic and gravity-wave character makes them sensitive to the structure of the deep stellar interior. In some cases, identifying these differences led to improvements in the final models presented here and in Paper I; here we illustrate particular examples of this. Conclusions. Further improvements in stellar modelling are required in order fully to utilise the observational accuracy to probe intrinsic limitations in the modelling. However, our analysis of the frequency differences and their relation to stellar internal properties provides a striking illustration of the potential of the mixed modes of red-giant stars for the diagnostics of stellar interiors. △ Less

Submitted 10 February, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: Astronomy and Astrophysics, in the press. Paper resubmitted to arXiv to correct author list in metadata

The Aarhus Red Giants Challenge I: Stellar structures in the red giant branch phase

Authors: V. Silva Aguirre, J. Christensen-Dalsgaard, S. Cassisi, M. Miller Bertolami, A. M. Serenelli, D. Stello, A. Weiss, G. Angelou, C. Jiang, Y. Lebreton, F. Spada, E. P. Bellinger, S. Deheuvels, R. M. Ouazzani, A. Pietrinferni, J. R. Mosumgaard, R. H. D. Townsend, T. Battich, D. Bossini, T. Constantino, P. Eggenberger, S. Hekker, A. Mazumdar, A. Miglio, K. B. Nielsen , et al. (1 additional authors not shown)

Abstract: (Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our c… ▽ More (Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at fixed radius along the red-giant branch for masses of 1.0 M$_\odot$, 1.5 M$_\odot$, 2.0 M$_\odot$, and 2.5 M$_\odot$, and compared the predicted stellar properties. Once models have been calibrated on the main sequence we find a residual spread in the predicted effective temperatures across all codes of ~20 K at solar radius and ~30-40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2-5% (increasing with stellar mass) which we track down to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ~0.1 mag in the optical V-band. We also compare the predicted [C/N] abundance ratio and found a spread of 0.1 dex or more for all considered masses. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and call for further investigations on these matters in the context of using properties of red giants as benchmarks for astrophysical studies. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 635, A164 (2020)

Testing the entropy calibration of the radii of cool stars: models of alpha Centauri A and B

Authors: F. Spada, P. Demarque

Abstract: We present models of alpha Centauri A and B implementing an entropy calibration of the mixing-length parameter alpha_MLT, recently developed and successfully applied to the Sun (Spada et al. 2018, ApJ, 869, 135). In this technique the value of alpha_MLT in the 1D stellar evolution code is calibrated to match the adiabatic specific entropy derived from 3D radiation-hydrodynamics simulations of stel… ▽ More We present models of alpha Centauri A and B implementing an entropy calibration of the mixing-length parameter alpha_MLT, recently developed and successfully applied to the Sun (Spada et al. 2018, ApJ, 869, 135). In this technique the value of alpha_MLT in the 1D stellar evolution code is calibrated to match the adiabatic specific entropy derived from 3D radiation-hydrodynamics simulations of stellar convective envelopes, whose effective temperature, surface gravity, and metallicity are selected consistently along the evolutionary track. The customary treatment of convection in stellar evolution models relies on a constant, solar-calibrated alpha_MLT. There is, however, mounting evidence that this procedure does not reproduce the observed radii of cool stars satisfactorily. For instance, modelling alpha Cen A and B requires an ad-hoc tuning of alpha_MLT to distinct, non-solar values. The entropy-calibrated models of alpha Cen A and B reproduce their observed radii within 1% (or better) without externally adjusted parameters. The fit is of comparable quality to that of models with freely adjusted alpha_MLT for alpha Cen B (within 1 sigma), while it is less satisfactory for alpha Cen A (within ~ 2.5 sigma). This level of accuracy is consistent with the intrinsic uncertainties of the method. Our results demonstrate the capability of the entropy calibration method to produce stellar models with radii accurate within 1%. This is especially relevant in characterising exoplanet-host stars and their planetary systems accurately. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Comments: 9 pages, 7 figures; accepted for publication in Monthly Notices of the Royal Astronomical Society

On the competing effect of wind braking and interior coupling in the rotational evolution of solar-like stars

Authors: F. Spada, A. C. Lanzafame

Abstract: Solar-like stars (M < 1.3 Msun) lose angular momentum through their magnetized winds. The resulting evolution of the surface rotation period, which can be directly measured photometrically, has the potential to provide an accurate indicator of stellar age, and is constrained by observations of rotation periods of coeval stars, such as members of Galactic open clusters. A prominent observational fe… ▽ More Solar-like stars (M < 1.3 Msun) lose angular momentum through their magnetized winds. The resulting evolution of the surface rotation period, which can be directly measured photometrically, has the potential to provide an accurate indicator of stellar age, and is constrained by observations of rotation periods of coeval stars, such as members of Galactic open clusters. A prominent observational feature of the mass-rotation period diagrams of open clusters is a sequence of relatively slower rotators. The formation and persistence of this slow rotators sequence across several billion years imply an approximately coherent spin-down of the stars that belong to it. In particular, the sequence is observed to evolve coherently towards longer periods in progressively older clusters. Recent observations of the 700-Myr old Praesepe and the 1-Gyr old NGC 6811 clusters, however, seem to contradict this general pattern. While the 1 Msun stars on the slow rotators sequence of the older NGC 6811 have longer periods than their counterparts in the younger Praesepe, as expected, the two sequences essentially merge at lower masses (< 0.8 Msun). In other words, low-mass stars seem to have not been spinning down at all in the intervening 300 Myr. Here we show that this behavior is a manifestation of the variable rotational coupling in solar-like stars. The resurfacing of angular momentum from the interior can temporarily compensate for the loss due to wind braking at the surface. The internal redistribution of angular momentum has a steep mass dependence; as a result, the re-coupling occurs at different ages for stars of different masses. Our model explains the stalled surface spin-down of low-mass stars between Praesepe and NGC 6811, and predicts that the same behavior should be observable at other ages in other mass ranges. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: 9 pages, 7 figures, 1 table. Submitted to Astronomy and Astrophysics

Asteroseismology of evolved stars to constrain the internal transport of angular momentum. I. Efficiency of transport during the subgiant phase

Authors: P. Eggenberger, S. Deheuvels, A. Miglio, S. Ekström, C. Georgy, G. Meynet, N. Lagarde, S. Salmon, G. Buldgen, J. Montalbán, F. Spada, J. Ballot

Abstract: Context: The observations of solar-like oscillations in evolved stars have brought important constraints on their internal rotation rates. To correctly reproduce these data, an efficient transport mechanism is needed in addition to meridional circulation and shear instability. Aims: We study the efficiency of the transport of angular momentum during the subgiant phase. Results: The precise asteros… ▽ More Context: The observations of solar-like oscillations in evolved stars have brought important constraints on their internal rotation rates. To correctly reproduce these data, an efficient transport mechanism is needed in addition to meridional circulation and shear instability. Aims: We study the efficiency of the transport of angular momentum during the subgiant phase. Results: The precise asteroseismic measurements of both core and surface rotation rates available for the six Kepler targets enable a precise determination of the efficiency of the transport of angular momentum needed for each of these subgiants. These results are found to be insensitive to all the uncertainties related to the modelling of rotational effects before the post-main sequence phase. An interesting exception in this context is the case of young subgiants (typical values of log(g) close to 4), because their rotational properties are sensitive to the degree of radial differential rotation on the main sequence. These young subgiants constitute therefore perfect targets to constrain the transport of angular momentum on the main sequence from asteroseismic observations of evolved stars. As for red giants, we find that the efficiency of the additional transport process increases with the mass of the star during the subgiant phase. However, the efficiency of this undetermined mechanism decreases with evolution during the subgiant phase, contrary to what is found for red giants. Consequently, a transport process with an efficiency that increases with the degree of radial differential rotation cannot account for the core rotation rates of subgiants, while it correctly reproduces the rotation rates of red giant stars. This suggests that the physical nature of the additional mechanism needed for the internal transport of angular momentum may be different in subgiant and red giant stars. △ Less

Submitted 12 December, 2018; originally announced December 2018.

Comments: 10 pages, 11 figures, accepted for publication in A&A

Improved calibration of the radii of cool stars based on 3D simulations of convection: implications for the solar model

Authors: F. Spada, P. Demarque, S. Basu, J. D. Tanner

Abstract: Main sequence, solar-like stars (M < 1.5 Msun) have outer convective envelopes that are sufficiently thick to affect significantly their overall structure. The radii of these stars, in particular, are sensitive to the details of inefficient, super-adiabatic convection occurring in their outermost layers. The standard treatment of convection in stellar evolution models, based on the Mixing-Length T… ▽ More Main sequence, solar-like stars (M < 1.5 Msun) have outer convective envelopes that are sufficiently thick to affect significantly their overall structure. The radii of these stars, in particular, are sensitive to the details of inefficient, super-adiabatic convection occurring in their outermost layers. The standard treatment of convection in stellar evolution models, based on the Mixing-Length Theory (MLT), provides only a very approximate description of convection in the super-adiabatic regime. Moreover, it contains a free parameter, alpha_MLT, whose standard calibration is based on the Sun, and is routinely applied to other stars ignoring the differences in their global parameters (e.g., effective temperature, gravity, chemical composition) and previous evolutionary history. In this paper, we present a calibration of alpha_MLT based on three-dimensional radiation-hydrodynamics (3D RHD) simulations of convection. The value of alpha_MLT is adjusted to match the specific entropy in the deep, adiabatic layers of the convective envelope to the corresponding value obtained from the 3D RHD simulations, as a function of the position of the star in the (log g, log T_eff) plane and its chemical composition. We have constructed a model of the present-day Sun using such entropy-based calibration. We find that its past luminosity evolution is not affected by the entropy calibration. The predicted solar radius, however, exceeds that of the standard model during the past several billion years, resulting in a lower surface temperature. This illustrative calculation also demonstrates the viability of the entropy approach for calibrating the radii of other late-type stars. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: 16 pages, 14 figures, accepted for publication in the Astrophysical Journal

Solar radius and luminosity variations induced by the internal dynamo magnetic fields

Authors: F. Spada, R. Arlt, M. Kueker, S. Sofia

Abstract: Although the occurrence of solar irradiance variations induced by magnetic surface features (e.g., sunspots, faculae, magnetic network) is generally accepted, the existence of intrinsic luminosity changes due to the internal magnetic fields is still controversial. This additional contribution is expected to be accompanied by radius variations, and to be potentially significant for the climate of t… ▽ More Although the occurrence of solar irradiance variations induced by magnetic surface features (e.g., sunspots, faculae, magnetic network) is generally accepted, the existence of intrinsic luminosity changes due to the internal magnetic fields is still controversial. This additional contribution is expected to be accompanied by radius variations, and to be potentially significant for the climate of the Earth. We aim to constrain theoretically the radius and luminosity variations of the Sun that are due to the effect of the variable magnetic fields in its interior associated with the dynamo cycle. We have extended a one-dimensional stellar evolution code to include several effects of the magnetic fields on the interior structure. We investigate different magnetic configurations, based on both observational constraints and on the output of state-of-the-art mean field dynamo models. We explore both step-like and simply periodic time dependences of the magnetic field peak strength. We find that the luminosity and radius variations are in anti-phase and in phase, respectively, with the magnetic field strength. For peak magnetic field strengths of the order of tens of kilogauss, luminosity variations ranging between 10^{-6} and 10^{-3} (in modulus) and radius variations between 10^{-6} and 10^{-5} are obtained. Modest but significant radius variations (up to 10^{-5} in relative terms) are obtained for magnetic fields of realistic strength and geometry, providing a potentially observable signature of the intrinsic variations. Establishing their existence in addition to the accepted surface effects would have very important implications for the understanding of solar-induced long-term trends on climate. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 18 pages, 7 figures; accepted for publication in Astronomische Nachrichten

Evidence of New Magnetic Transitions in Late-Type Dwarfs from Gaia DR2

Authors: Alessandro C. Lanzafame, Elisa Distefano, Sydney A. Barnes, Federico Spada

Abstract: The second Gaia data release contains the identification of 147 535 low-mass ($\le 1.4 M_{\odot}$) rotational modulation variable candidates on (or close to) the main sequence, together with their rotation period and modulation amplitude. The richness, the period and amplitude range, and the photometric precision of this sample make it possible to unveil, for the first time, signatures of differen… ▽ More The second Gaia data release contains the identification of 147 535 low-mass ($\le 1.4 M_{\odot}$) rotational modulation variable candidates on (or close to) the main sequence, together with their rotation period and modulation amplitude. The richness, the period and amplitude range, and the photometric precision of this sample make it possible to unveil, for the first time, signatures of different surface inhomogeneity regimes in the amplitude-period density diagram. The modulation amplitude distribution shows a clear bimodality, with an evident gap at periods $P \le 2$ d. The low amplitude branch, in turn, shows a period bimodality with a main clustering at periods $P \approx$ 5 - 10 d and a secondary clustering of ultra-fast rotators at $P \le 0.5$ d. The amplitude-period multimodality is correlated with the position in the period-absolute magnitude (or period-color) diagram, with the low- and high-amplitude stars occupying different preferential locations. Here we argue that such a multimodality represents a further evidence of the existence of different regimes of surface inhomogeneities in young and middle-age low-mass stars and we lay out possible scenarios for their evolution, which manifestly include rapid transitions from one regime to another. In particular, the data indicate that stars spinning up close to break-up velocity undergo a very rapid change in their surface inhomogeneities configuration, which is revealed here for the first time. The multimodality can be exploited to identify field stars of age $\sim$ 100 -- 600 Myr belonging to the slow-rotator low-amplitude sequence, for which age can be estimated from the rotation period via gyrochronology relationships. △ Less

Submitted 16 April, 2019; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: 15 pages, 6 figures, Accepted by ApJ

Effects of a subadiabatic layer on convection and dynamos in spherical wedge simulations

Authors: Petri J. Käpylä, Mariangela Viviani, Maarit J. Käpylä, Axel Brandenburg, Federico Spada

Abstract: We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We fi… ▽ More We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We find that the convective heat transport is strongly concentrated toward the equatorial and polar regions in the cases without a substantial radiative layer below the convection zone. The presence of a stable layer below the convection zone significantly reduces the anisotropy of radial enthalpy transport. Furthermore, the dynamo solutions are sensitive to subtle changes in the convection zone structure. We find that the kinetic helicity changes sign in the deeper parts of the convection zone at high latitudes in all runs. This region expands progressively toward the equator in runs with a thicker stably stratified layer. △ Less

Submitted 2 January, 2019; v1 submitted 15 March, 2018; originally announced March 2018.

Comments: 32 pages, 18 figures, accepted to GAFD Special issue on 'Recent Developments in Natural Dynamos'

Report number: NORDITA 2018-022

Journal ref: Geophys. Astrophys. Fluid Dyn. 113, 149-183 (2020)

PEPSI deep spectra. III. A chemical analysis of the ancient planet-host star Kepler-444

Authors: C. E. Mack III, K. G. Strassmeier, I. Ilyin, S. C. Schuler, F. Spada, S. A. Barnes

Abstract: We obtained an LBT/PEPSI spectrum with very high resolution and high signal-to-noise ratio (S/N) of the K0V host Kepler-444, which is known to host 5 sub-Earth size rocky planets. The spectrum has a resolution of R=250,000, a continuous wavelength coverage from 4230 to 9120A, and S/N between 150 and 550:1 (blue to red). We performed a detailed chemical analysis to determine the photospheric abunda… ▽ More We obtained an LBT/PEPSI spectrum with very high resolution and high signal-to-noise ratio (S/N) of the K0V host Kepler-444, which is known to host 5 sub-Earth size rocky planets. The spectrum has a resolution of R=250,000, a continuous wavelength coverage from 4230 to 9120A, and S/N between 150 and 550:1 (blue to red). We performed a detailed chemical analysis to determine the photospheric abundances of 18 chemical elements, in order to use the abundances to place constraints on the bulk composition of the five rocky planets. Our spectral analysis employs the equivalent width method for most of our spectral lines, but we used spectral synthesis to fit a small number of lines that require special care. In both cases, we derived our abundances using the MOOG spectral analysis package and Kurucz model atmospheres. We find no correlation between elemental abundance and condensation temperature among the refractory elements. In addition, using our spectroscopic stellar parameters and isochrone fitting, we find an age of 10+/-1.5 Gyr, which is consistent with the asteroseismic age of 11+/-1 Gyr. Finally, from the photospheric abundances of Mg, Si, and Fe, we estimate that the typical Fe-core mass fraction for the rocky planets in the Kepler-444 system is approximately 24 per cent. If our estimate of the Fe-core mass fraction is confirmed by more detailed modeling of the disk chemistry and simulations of planet formation and evolution in the Kepler-444 system, then this would suggest that rocky planets in more metal-poor and alpha-enhanced systems may tend to be less dense than their counterparts of comparable size in more metal-rich systems. △ Less

Submitted 19 December, 2017; originally announced December 2017.

Comments: in press, 11 pages, 3 figures, data available from pepsi.aip.de

Journal ref: A&A 612, A46 (2018)

Activity cycles in members of young loose stellar associations

Authors: E. Distefano, A. C. Lanzafame, A. Lanza, S. Messina, F. Spada

Abstract: Magnetic cycles have been detected in tens of solar-like stars. The relationship between the cycle properties and global stellar parameters is not fully understood yet. We searched for activity cycles in 90 solar-like stars with ages between 4 and 95 Myr aiming to investigate the properties of activity cycles in this age range. We measured the length $P_{ cyc}$ of a given cycle by analyzing th… ▽ More Magnetic cycles have been detected in tens of solar-like stars. The relationship between the cycle properties and global stellar parameters is not fully understood yet. We searched for activity cycles in 90 solar-like stars with ages between 4 and 95 Myr aiming to investigate the properties of activity cycles in this age range. We measured the length $P_{ cyc}$ of a given cycle by analyzing the long-term time-series of three activity indexes. For each star, we computed also the global magnetic activity index <IQR> that is proportional to the amplitude of the rotational modulation and is a proxy of the mean level of the surface magnetic activity. We detected activity cycles in 67 stars. Secondary cycles were also detected in 32 stars. The lack of correlation between $P_{ cyc}$ and $P_{ rot}$ suggest that these stars belong to the Transitional Branch and that the dynamo acting in these stars is different from the solar one. This statement is also supported by the analysis of the butterfly diagrams. We computed the Spearman correlation coefficient $r_{ S}$ between $P_{ cyc}$, <IQR> and different stellar parameters. We found that $P_{ cyc}$ is uncorrelated with all the investigated parameters. The <IQR> index is positively correlated with the convective turn-over time-scale, the magnetic diffusivity time-scale $τ_{ diff}$, and the dynamo number $D_{ N}$, whereas it is anti-correlated with the effective temperature $T_{ eff}$, the photometric shear $ΔΩ_{\rm phot}$ and the radius $R_{ C}$ at which the convective zone is located. We found that $P_{ cyc}$ is about constant and that <IQR> decreases with the stellare age in the range 4-95 Myr. We investigated the magnetic activity of AB Dor A by merging ASAS time-series with previous long-term photometric data. We estimated the length of the AB Dor A primary cycle as $P_{ cyc} = 16.78 \pm 2 \rm yr$. △ Less

Submitted 24 June, 2017; originally announced June 2017.

Comments: 19 pages , 15 figures, accepted

Journal ref: A&A 606, A58 (2017)

Cosmic-ray electron+positron spectrum from 7 GeV to 2 TeV with the Fermi Large Area Telescope

Authors: Fermi-LAT Collaboration, :, S. Abdollahi, M. Ackermann, M. Ajello, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. D. Bloom, R. Bonino, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, R. A. Cameron, R. Caputo, M. Caragiulo, D. Castro, E. Cavazzuti, C. Cecchi, A. Chekhtman, S. Ciprini, J. Cohen-Tanugi , et al. (76 additional authors not shown)

Abstract: We present a measurement of the cosmic-ray electron+positron spectrum between 7 GeV and 2 TeV performed with almost seven years of data collected with the Fermi Large Area Telescope. We find that the spectrum is well fit by a broken power law with a break energy at about 50 GeV. Above 50 GeV, the spectrum is well described by a single power law with a spectral index of… ▽ More We present a measurement of the cosmic-ray electron+positron spectrum between 7 GeV and 2 TeV performed with almost seven years of data collected with the Fermi Large Area Telescope. We find that the spectrum is well fit by a broken power law with a break energy at about 50 GeV. Above 50 GeV, the spectrum is well described by a single power law with a spectral index of $3.07 \pm 0.02 \; (\text{stat+syst}) \pm 0.04 \; (\text{energy measurement})$. An exponential cutoff lower than 1.8 TeV is excluded at 95\% CL. △ Less

Submitted 24 April, 2017; originally announced April 2017.

Journal ref: Phys. Rev. D 95, 082007 (2017)

Lithium abundance and 6Li/7Li ratio in the active giant HD123351 I. A comparative analysis of 3D and 1D NLTE line-profile fits

Authors: A. Mott, M. Steffen, E. Caffau, F. Spada, K. G. Strassmeier

Abstract: Current three-dimensional (3D) hydrodynamical model atmospheres together with NLTE spectrum synthesis, permit to derive reliable atomic and isotopic chemical abundances from high-resolution stellar spectra. Not much is known about the presence of the fragile 6Li isotope in evolved solar-metallicity RGB stars, not to mention its production in magnetically active targets like HD123351. From fits of… ▽ More Current three-dimensional (3D) hydrodynamical model atmospheres together with NLTE spectrum synthesis, permit to derive reliable atomic and isotopic chemical abundances from high-resolution stellar spectra. Not much is known about the presence of the fragile 6Li isotope in evolved solar-metallicity RGB stars, not to mention its production in magnetically active targets like HD123351. From fits of the observed CFHT spectrum with synthetic line profiles based on 1D and 3D model atmospheres, we seek to estimate the abundance of the 6Li isotope and to place constraints on its origin. We derive A(Li) and the 6Li/7Li isotopic ratio by fitting different synthetic spectra to the Li-line region of a high-resolution CFHT spectrum (R=120 000, S/R=400). The synthetic spectra are computed with four different line lists, using in parallel 3D hydrodynamical CO5BOLD and 1D LHD model atmospheres and treating the line formation of the lithium components in non-LTE (NLTE). We find A(Li)=1.69+/-0.11 dex and 6Li/7Li=8.0+/-4.4 % in 3D-NLTE, using the line list of Meléndez et al. (2012), updated with new atomic data for V I, which results in the best fit of the lithium line profile of HD123351. Two other line lists lead to similar results but with inferior fit qualities. Our 2-sigma detection of the 6Li isotope is the result of a careful statistical analysis and the visual inspection of each achieved fit. Since the presence of a significant amount of 6Li in the atmosphere of a cool evolved star is not expected in the framework of standard stellar evolution theory, non-standard, external lithium production mechanisms, possibly related to stellar activity or a recent accretion of rocky material, need to be invoked to explain the detection of 6Li in HD123351. △ Less

Submitted 21 April, 2017; originally announced April 2017.

Comments: 16 pages, 11 figures. Accepted for publication in A&A

Journal ref: A&A 604, A44 (2017)

The Sirius System and its Astrophysical Puzzles: Hubble Space Telescope and Ground-Based Astrometry

Authors: Howard E. Bond, Gail H. Schaefer, Ronald L. Gilliland, Jay B. Holberg, Brian D. Mason, Irving W. Lindenblad, Miranda Seitz-McLeese, W. David Arnett, Pierre Demarque, Federico Spada, Patrick A. Young, Martin A. Barstow, Matthew R. Burleigh, Donald Gudehus

Abstract: Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, brightest star in the sky, orbited in a 50.13-year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years, and nearly 2300 historica… ▽ More Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, brightest star in the sky, orbited in a 50.13-year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years, and nearly 2300 historical measurements dating back to the 19th century, we determine precise orbital elements for the visual binary. Combined with the parallax and the motion of the A component, these elements yield dynamical masses of 2.063+/-0.023 Msun and 1.018+/-0.011 Msun for Sirius A and B, respectively. Our precise HST astrometry rules out third bodies orbiting either star in the system, down to masses of ~15-25 Mjup. The location of Sirius B in the H-R diagram is in excellent agreement with theoretical cooling tracks for WDs of its dynamical mass, and implies a cooling age of ~126 Myr. The position of Sirius B in the mass-radius plane is also consistent with WD theory, assuming a carbon-oxygen core. Including the pre-WD evolutionary timescale of the assumed progenitor, the total age of Sirius B is about 228+/-10 Myr. We calculated evolutionary tracks for stars with the dynamical mass of Sirius A, using two independent codes. We find it necessary to assume a slightly sub-solar metallicity, of about 0.85 Zsun, to fit its location in the luminosity-radius plane. The age of Sirius A based on these models is about 237-247 Myr, with uncertainties of +/-15 Myr, consistent with that of the WD companion. We discuss astrophysical puzzles presented by the Sirius system, including the probability that the two stars must have interacted in the past, even though there is no direct evidence for this, and the orbital eccentricity remains high. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: Accepted by Astrophysical Journal

The Yale-Potsdam Stellar Isochrones (YaPSI)

Authors: F. Spada, P. Demarque, Y. -C. Kim, T. S. Boyajian, J. M. Brewer

Abstract: We introduce the Yale-Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei-Yale database, special emphasis is placed on the construction of accurate low-mass models (Mstar < 0.6 Msun), and in particular of their mass-luminosity and mass-radius relations, both crucial in characterizing exoplanet-hos… ▽ More We introduce the Yale-Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei-Yale database, special emphasis is placed on the construction of accurate low-mass models (Mstar < 0.6 Msun), and in particular of their mass-luminosity and mass-radius relations, both crucial in characterizing exoplanet-host stars and, in turn, their planetary systems. The YaPSI models cover the mass range 0.15 to 5.0 Msun, densely enough to permit detailed interpolation in mass, and the metallicity and helium abundance ranges [Fe/H] = -1.5 to +0.3, and Y = 0.25 to 0.37, specified independently of each other (i.e., no fixed Delta Y/Delta Z relation is assumed). The evolutionary tracks are calculated from the pre-main sequence up to the tip of the red giant branch. The isochrones, with ages between 1 Myr and 20 Gyr, provide UBVRI colors in the Johnson-Cousins system, and JHK colors in the homogeneized Bessell & Brett system, derived from two different semi-empirical Teff-color calibrations from the literature. We also provide utility codes, such as an isochrone interpolator in age, metallicity, and helium content, and an interface of the tracks with an open-source Monte Carlo Markov-Chain tool for the analysis of individual stars. Finally, we present comparisons of the YaPSI models with the best empirical mass- luminosity and mass-radius relations available to date, as well as isochrone fitting of well-studied ste △ Less

Submitted 11 March, 2017; originally announced March 2017.

Comments: 17 pages, 14 figures; accepted for publication in the Astrophysical Journal

Theoretical interpretation of Pass 8 Fermi-LAT $e^+ + e^-$ data

Authors: Mattia Di Mauro, Silvia Manconi, Andrea Vittino, Fiorenza Donato, Nicolao Fornengo, Luca Baldini, Raffaella Bonino, Nicolò Di Lalla, Luca Latronico, Simone Maldera, Alberto Manfreda, Michela Negro, Melissa Pesce-Rollins, Carmelo Sgrò, Francesca Spada

Abstract: The flux of positrons and electrons ($e^+ + e^-$) has been measured by the $Fermi$ Large Area Telescope (LAT) in the energy range between 7 GeV and 2 TeV. We discuss a number of interpretations of Pass 8 $Fermi$-LAT $e^+ + e^-$ spectrum, combining electron and positron emission from supernova remnants (SNRs) and pulsar wind nebulae (PWNe), or produced by the collision of cosmic rays with the inter… ▽ More The flux of positrons and electrons ($e^+ + e^-$) has been measured by the $Fermi$ Large Area Telescope (LAT) in the energy range between 7 GeV and 2 TeV. We discuss a number of interpretations of Pass 8 $Fermi$-LAT $e^+ + e^-$ spectrum, combining electron and positron emission from supernova remnants (SNRs) and pulsar wind nebulae (PWNe), or produced by the collision of cosmic rays with the interstellar medium. We find that the $Fermi$-LAT spectrum is compatible with the sum of electrons from a smooth SNR population, positrons from cataloged PWNe, and a secondary component. If we include in our analysis constraints from AMS-02 positron spectrum, we obtain a slightly worse fit to the $e^+ + e^-$ $Fermi$-LAT spectrum, depending on the propagation model. As an additional scenario, we replace the smooth SNR component within 0.7 kpc with the { individual sources} found in Green's catalog of Galactic SNRs. We find that separate consideration of far and near sources helps to reproduce the $e^+ + e^-$ $Fermi$-LAT spectrum. However, we show that the fit degrades when the radio constraints on the positron emission from Vela SNR (which is the main contributor at high energies) are taken into account. We find that a break in the power-law injection spectrum at about 100 GeV can also reproduce the measured $e^+ + e^-$ spectrum} and, among the cosmic-ray propagation models that we consider, no reasonable break of the power-law dependence of the diffusion coefficient can modify the electron flux enough to reproduce the observed shape. △ Less

Submitted 12 November, 2020; v1 submitted 1 March, 2017; originally announced March 2017.

Comments: 15 pages, 12 figures. The Astrophysical Journal, Volume 845, Number 2

Report number: TUM-HEP-1078/17

Fermi-LAT Observations of High-energy Behind-the-limb Solar Flares

Authors: M. Ackermann, A. Allafort, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, E. Charles, S. Ciprini, F. Costanza, S. Cutini, F. D'Ammando, F. de Palma, R. Desiante, S. W. Digel , et al. (64 additional authors not shown)

Abstract: We report on the Fermi-LAT detection of high-energy emission from the behind-the-limb (BTL) solar flares that occurred on 2013 October 11, and 2014 January 6 and September 1. The Fermi-LAT observations are associated with flares from active regions originating behind both the eastern and western limbs, as determined by STEREO. All three flares are associated with very fast coronal mass ejections (… ▽ More We report on the Fermi-LAT detection of high-energy emission from the behind-the-limb (BTL) solar flares that occurred on 2013 October 11, and 2014 January 6 and September 1. The Fermi-LAT observations are associated with flares from active regions originating behind both the eastern and western limbs, as determined by STEREO. All three flares are associated with very fast coronal mass ejections (CMEs) and strong solar energetic particle events. We present updated localizations of the >100 MeV photon emission, hard X-ray (HXR)and EUV images, and broadband spectra from 10 keV to 10 GeV, as well as microwave spectra. We also provide a comparison of the BTL flares detected by Fermi-LAT with three on-disk flares and present a study of some of the significant quantities of these flares as an attempt to better understand the acceleration mechanisms at work during these occulted flares. We interpret the HXR emission to be due to electron bremsstrahlung from a coronal thin-target loop top with the accelerated electron spectra steepening at semirelativistic energies. The >100 MeV gamma-rays are best described by a pion-decay model resulting from the interaction of protons (and other ions) in a thick-target photospheric source. The protons are believed to have been accelerated (to energies >10 GeV) in the CME environment and precipitate down to the photosphere from the downstream side of the CME shock and landed on the front side of the Sun, away from the original flare site and the HXR emission. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 13 pages, 14 figures published on ApJ

Journal ref: ApJ 835, 219, 1 February 2017

Search for extended sources in the Galactic Plane using 6 years of Fermi-Large Area Telescope Pass 8 data above 10 GeV

Authors: The Fermi LAT Collaboration, M. Ackermann, M. Ajello, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, R. A. Cameron, M. Caragiulo, P. A. Caraveo, D. Castro, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, C. C. Cheung , et al. (95 additional authors not shown)

Abstract: The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two… ▽ More The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain. △ Less

Submitted 11 April, 2018; v1 submitted 1 February, 2017; originally announced February 2017.

Comments: 33 pages, 22 figures & 3 tables. Published by The Astrophysical Journal. Available on the Fermi Science Support Center (FSSC) together with the 3FHL catalog

The second catalog of flaring gamma-ray sources from the Fermi All-sky Variability Analysis

Authors: S. Abdollahi, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, S. Buson, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, A. Chekhtman , et al. (102 additional authors not shown)

Abstract: We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data relea… ▽ More We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6σ$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2. △ Less

Submitted 12 September, 2017; v1 submitted 9 December, 2016; originally announced December 2016.

Comments: Accepted for publication by The Astrophysical Journal Supplement. Online analysis results available at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/FAVA/index.php . Corresponding authors: M. Giomi, R. Buehler, D. Kocevski, and M. Ajello

Gamma-ray blazar spectra with H.E.S.S. II mono analysis: the case of PKS 2155-304 and PG 1553+113

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, F. Brun , et al. (311 additional authors not shown)

Abstract: The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft sp… ▽ More The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 13 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 600, A89 (2017)

Evidence of radius inflation in stars approaching the slow-rotator sequence

Authors: Alessandro C. Lanzafame, Federico Spada, Elisa Distefano

Abstract: Average stellar radii in open clusters can be estimated from rotation periods and projected rotational velocities under the assumption of random orientation of the spin axis. Such estimates are independent of distance, interstellar absorption, and models, but their validity can be limited by missing data (truncation) or data that only represent upper/lower limits (censoring). We present a new stat… ▽ More Average stellar radii in open clusters can be estimated from rotation periods and projected rotational velocities under the assumption of random orientation of the spin axis. Such estimates are independent of distance, interstellar absorption, and models, but their validity can be limited by missing data (truncation) or data that only represent upper/lower limits (censoring). We present a new statistical analysis method to estimate average stellar radii in the presence of censoring and truncation. We use theoretical distribution functions of the projected stellar radius $R \sin i$ to define a likelihood function in the presence of censoring and truncation. Average stellar radii in magnitude bins are then obtained by a maximum likelihood parametric estimation procedure. This method is capable of recovering the average stellar radius within a few percent with as few as $\approx$ 10 measurements. Here it is applied for the first time to the dataset available for the Pleiades. We find an agreement better than $\approx$ 10 percent between the observed $R$ vs $M_K$ relationship and current standard stellar models for 1.2 $\ge M/M_{\odot} \ge$ 0.85 with no evident bias. Evidence of a systematic deviation at $2σ$ level are found for stars with 0.8 $\ge M/M_{\odot} \ge$ 0.6 approaching the slow-rotator sequence. Fast-rotators ($P$ < 2 d) agree with standard models within 15 percent with no systematic deviations in the whole 1.2 $\ge M/M_{\odot} \ge$ 0.5 range. The evidence found of a possible radius inflation just below the lower mass limit of the slow-rotator sequence indicates a possible connection with the transition from the fast to the slow-rotator sequence. △ Less

Submitted 21 September, 2016; v1 submitted 6 September, 2016; originally announced September 2016.

Comments: Accepted by Astronomy and Astrophysics, 11 pages, 6 figures

Journal ref: A&A 597, A63 (2017)

Sunspot positions, areas, and group tilt angles for 1611-1631 from observations by Christoph Scheiner

Authors: R. Arlt, V. Senthamizh Pavai, C. Schmiel, F. Spada

Abstract: Digital images of the observations printed in the books "Rosa Ursina sive solis" and "Prodromus pro sole mobili" by Christoph Scheiner as well as the drawings from Scheiner's letters to Marcus Welser are analysed in order to obtain information on positions and sizes of sunspots that appeared before the Maunder minimum. In most cases, the given orientation of the ecliptic is used to set up the heli… ▽ More Digital images of the observations printed in the books "Rosa Ursina sive solis" and "Prodromus pro sole mobili" by Christoph Scheiner as well as the drawings from Scheiner's letters to Marcus Welser are analysed in order to obtain information on positions and sizes of sunspots that appeared before the Maunder minimum. In most cases, the given orientation of the ecliptic is used to set up the heliographic coordinate system for the drawings. Positions and sizes are measured manually on the screen. Very early drawings have no indication of their orientation. A rotational matching using common spots of adjacent days is used in some cases, while in other cases, the assumption of images being aligned with a zenith-horizon coordinate system appeared to be the most probable. In total, 8167 sunspots were measured. A distribution of sunspot latitudes versus time (butterfly diagram) is obtained for Scheiner's observations. The observations of 1611 are very inaccurate, the drawings of 1612 have at least an indication of their orientation, while the remaining part of the spot positions from 1618-1631 have good to very good accuracy. We also computed 697 tilt angles of apparently bipolar sunspot groups observed in the period 1618-1631. We find that the average tilt angle of nearly 4 degrees is not significantly different from 20th-century values. △ Less

Submitted 25 August, 2016; originally announced August 2016.

Journal ref: A&A 595, A104 (2016)

Searching the Gamma-ray Sky for Counterparts to Gravitational Wave Sources: Fermi GBM and LAT Observations of LVT151012 and GW151226

Authors: J. L. Racusin, E. Burns, A. Goldstein, V. Connaughton, C. A. Wilson-Hodge, P. Jenke, L. Blackburn, M. S. Briggs, J. Broida, J. Camp, N. Christensen, C. M. Hui, T. Littenberg, P. Shawhan, L. Singer, J. Veitch, P. N. Bhat, W. Cleveland, G. Fitzpatrick, M. H. Gibby, A. von Kienlin, S. McBreen, B. Mailyan, C. A. Meegan, W. S. Paciesas , et al. (116 additional authors not shown)

Abstract: We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techn… ▽ More We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: 27 pages, 1 table, 11 figures, Submitted to ApJ

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the dif… ▽ More This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands. △ Less

Submitted 21 July, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: For the main Letter, see arXiv:1602.08492

Report number: LIGO-P1600137-v2

Journal ref: The Astrophysical Journal Supplement Series, 225:8 (15pp), 2016 July

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

Authors: M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, A. Chekhtman, J. Chiang, G. Chiaro, S. Ciprini , et al. (90 additional authors not shown)

Abstract: We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interac… ▽ More We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data. △ Less

Submitted 13 April, 2016; v1 submitted 12 April, 2016; originally announced April 2016.

Comments: Published by PRD, 16 pages, 11 figures, corresponding authors: F. Loparco and M. N. Mazziotta

Lower limit for differential rotation in members of young loose stellar associations

Authors: E. Distefano, A. C. Lanzafame, A. F. Lanza, S. Messina, F. Spada

Abstract: Surface differential rotation (SDR) plays a key role in dynamo models. SDR estimates are therefore essential for constraining theoretical models. We measure a lower limit to SDR in a sample of solar-like stars belonging to young associations with the aim of investigating how SDR depends on global stellar parameters in the age range (4-95 Myr). The rotation period of a solar-like star can be recove… ▽ More Surface differential rotation (SDR) plays a key role in dynamo models. SDR estimates are therefore essential for constraining theoretical models. We measure a lower limit to SDR in a sample of solar-like stars belonging to young associations with the aim of investigating how SDR depends on global stellar parameters in the age range (4-95 Myr). The rotation period of a solar-like star can be recovered by analyzing the flux modulation caused by dark spots and stellar rotation. The SDR and the latitude migration of dark-spots induce a modulation of the detected rotation period. We employ long-term photometry to measure the amplitude of such a modulation and to compute the quantity DeltaOmega_phot =2p/P_min -2pi/P_max that is a lower limit to SDR. We find that DeltaOmega_phot increases with the stellar effective temperature and with the global convective turn-over time-scale tau_c. We find that DeltaOmega_phot is proportional to Teff^2.18pm 0.65 in stars recently settled on the ZAMS. This power law is less steep than those found by previous authors, but closest to recent theoretical models. We find that DeltaOmega_phot steeply increases between 4 and 30 Myr and that itis almost constant between 30 and 95 Myr in a 1 M_sun star. We find also that the relative shear increases with the Rossby number Ro. Although our results are qualitatively in agreement with hydrodynamical mean-field models, our measurements are systematically higher than the values predicted by these models. The discrepancy between DeltaOmega_phot measurements and theoretical models is particularly large in stars with periods between 0.7 and 2 d. Such a discrepancy, together with the anomalous SDR measured by other authors for HD 171488 (rotating in 1.31 d), suggests that the rotation period could influence SDR more than predicted by the models. △ Less

Submitted 7 April, 2016; originally announced April 2016.

Comments: 23 pages, 15 figures, 5 tables,accepted by Astronomy and Astrophysics

Journal ref: A&A 591, A43 (2016)

Angular momentum transport efficiency in post-main sequence low-mass stars

Authors: F. Spada, M. Gellert, R. Arlt, S. Deheuvels

Abstract: Context. Using asteroseismic techniques, it has recently become possible to probe the internal rotation profile of low-mass (~1.1-1.5 Msun) subgiant and red giant stars. Under the assumption of local angular momentum conservation, the core contraction and envelope expansion occurring at the end of the main sequence would result in a much larger internal differential rotation than observed. This su… ▽ More Context. Using asteroseismic techniques, it has recently become possible to probe the internal rotation profile of low-mass (~1.1-1.5 Msun) subgiant and red giant stars. Under the assumption of local angular momentum conservation, the core contraction and envelope expansion occurring at the end of the main sequence would result in a much larger internal differential rotation than observed. This suggests that angular momentum redistribution must be taking place in the interior of these stars. Aims. We investigate the physical nature of the angular momentum redistribution mechanisms operating in stellar interiors by constraining the efficiency of post-main sequence rotational coupling. Methods. We model the rotational evolution of a 1.25 Msun star using the Yale Rotational stellar Evolution Code. Our models take into account the magnetic wind braking occurring at the surface of the star and the angular momentum transport in the interior, with an efficiency dependent on the degree of internal differential rotation. Results. We find that models including a dependence of the angular momentum transport efficiency on the radial rotational shear reproduce very well the observations. The best fit of the data is obtained with an angular momentum transport coefficient scaling with the ratio of the rotation rate of the radiative interior over that of the convective envelope of the star as a power law of exponent ~3. This scaling is consistent with the predictions of recent numerical simulations of the Azimuthal Magneto-Rotational Instability. Conclusions. We show that an angular momentum transport process whose efficiency varies during the stellar evolution through a dependence on the level of internal differential rotation is required to explain the observed post-main sequence rotational evolution of low-mass stars. △ Less

Submitted 3 March, 2016; originally announced March 2016.

Comments: 8 pages, 6 figures; accepted for publication in Astronomy & Astrophysics

Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared wit… ▽ More A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams. △ Less

Submitted 21 July, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: For Supplement, see https://arxiv.org/abs/1604.07864

Report number: LIGO-P1500227-v12

Journal ref: The Astrophysical Journal Letters, 826:L13 (8pp), 2016 July 20

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Authors: F. Acero, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, C. Cecchi , et al. (109 additional authors not shown)

Abstract: Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Ga… ▽ More Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center. △ Less

Submitted 23 February, 2016; originally announced February 2016.

Comments: Accepted for publication in ApJS

Journal ref: ApJS 223 26 (2016)

Very-high-energy gamma-rays from the Universe's middle age: detection of the z=0.940 blazar PKS 1441+25 with MAGIC

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, W. Bednarek, E. Bernardini, B. Biassuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, A. Chatterjee, R. Clavero, P. Colin, E. Colombo , et al. (229 additional authors not shown)

Abstract: The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies… ▽ More The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy. △ Less

Submitted 12 January, 2018; v1 submitted 14 December, 2015; originally announced December 2015.

Comments: Corresponding Authors: J. Becerra (josefa.becerra@nasa.gov), M. Nievas Rosillo (mnievas@ucm.es), M. Manganaro (manganaro@iac.es), F. Tavecchio (fabrizio.tavecchio@brera.inaf.it) Published in The Astrophysical Journal Letters (ApJL), 2015ApJ...815L..23A, DOI: 10.1088/2041-8205/815/2/L23

Journal ref: ApJL, 2015, 815, L23

The 1st Fermi Lat Supernova Remnant Catalog

Authors: Fabio Acero, Markus Ackermann, Marco Ajello, Luca Baldini, Jean Ballet, Guido Barbiellini, Denis Bastieri, Ronaldo Bellazzini, E. Bissaldi, Roger Blandford, E. D. Bloom, Raffaella Bonino, Eugenio Bottacini, J. Bregeon, Philippe Bruel, Rolf Buehler, S. Buson, G. A. Caliandro, Rob A. Cameron, R Caputo, Micaela Caragiulo, Patrizia A. Caraveo, Jean Marc Casandjian, Elisabetta Cavazzuti, Claudia Cecchi , et al. (134 additional authors not shown)

Abstract: To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245… ▽ More To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: Resubmitted to ApJS

Journal ref: ApJS 224 8 (2016)

A color-period diagram for the open cluster M 48 (NGC 2548), and its rotational age

Authors: Sydney A. Barnes, Jörg Weingrill, Thomas Granzer, Federico Spada, Klaus G. Strassmeier

Abstract: Rotation periods are increasingly being used to derive ages for cool single field stars. Such ages are based on an empirical understanding of how cool stars spin down, acquired by constructing color-period diagrams (CPDs) for a series of open clusters. Our main aims here are to construct a CPD for M 48, to compare this with other clusters of similar age to check for consistency, and to derive a ro… ▽ More Rotation periods are increasingly being used to derive ages for cool single field stars. Such ages are based on an empirical understanding of how cool stars spin down, acquired by constructing color-period diagrams (CPDs) for a series of open clusters. Our main aims here are to construct a CPD for M 48, to compare this with other clusters of similar age to check for consistency, and to derive a rotational age for M 48 using gyrochronology. We monitored M 48 photometrically for over 2 months with AIP's STELLA I 1.2 m telescope and the WiFSIP 4K imager in Tenerife. Light curves with 3 mmag precision for bright (V~14 mag) stars were produced and then analysed to provide rotation periods. A cluster CPD has then been constructed. We report 62 rotation periods for cool stars in M 48. The CPD displays a clear slow/I-sequence of rotating stars, similar to those seen in the 625 Myr-old Hyades and 590 Myr-old Praesepe clusters, and below both, confirming that M 48 is younger. A similar comparison with the 250 Myr-old M 34 cluster shows that M 48 is older and does not possess any fast/C-sequence G or early K stars like those in M 34, although relatively fast rotators do seem to be present among the late-K and M stars. A more detailed comparison of the CPD with rotational evolution models shows that the cluster stars have a mean age of 450 Myr, and its (rotating) stars can be individually dated to +-117 Myr (26%). Much of this uncertainty stems from intrinsic astrophysical spread in initial periods, and almost all stars are consistent with a single age of 450 Myr. The gyro-age of M 48 as a whole is 450+-50 Myr, in agreement with the previously determined isochrone age of 400+-100 Myr. △ Less

Submitted 2 November, 2015; originally announced November 2015.

Journal ref: A&A 2015, 583, A73

Search for extended gamma-ray emission from the Virgo galaxy cluster with Fermi-LAT

Authors: M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, R. Caputo, M. Caragiulo, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti , et al. (96 additional authors not shown)

Abstract: Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $γ$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for search… ▽ More Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $γ$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3°that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $γ$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle σv \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $γ$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$. △ Less

Submitted 30 September, 2015; originally announced October 2015.

Comments: 15 pages, 11 figures, 4 tables, accepted for publication in ApJ; corresponding authors: T. Jogler, S. Zimmer & A. Pinzke

Journal ref: ApJ 812 159 (2015)

Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113

Authors: The Fermi LAT collaboration, M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, R. Caputo, M. Caragiulo, P. A. Caraveo , et al. (117 additional authors not shown)

Abstract: We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated osci… ▽ More We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity. △ Less

Submitted 12 October, 2015; v1 submitted 7 September, 2015; originally announced September 2015.

Comments: 8 pages, 5 figures. Accepted to The Astrophysical Journal Letters. Corresponding authors: S. Ciprini (ASDC/INFN), S. Cutini (ASDC/INFN), S. Larsson (Stockholm Univ/KTH), A. Stamerra (INAF/SNS), D. J. Thompson (NASA GSFC)

PSR J1906+0722: An Elusive Gamma-ray Pulsar

Authors: C. J. Clark, H. J. Pletsch, J. Wu, L. Guillemot, M. Ackermann, B. Allen, A. de Angelis, C. Aulbert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, O. Bock, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo , et al. (95 additional authors not shown)

Abstract: We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT s… ▽ More We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($Δf / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: 7 pages, 4 figures, accepted for publication in Astrophysical Journal Letters

Journal ref: ApJ 809 L2 (2015)

Rotational evolution of slow-rotators sequence stars

Authors: Alessandro C. Lanzafame, Federico Spada

Abstract: The observed mass-age-rotation relationship in open clusters shows the progressive development of a slow-rotators sequence. The observed clustering on this sequence suggests that it corresponds to some equilibrium or asymptotic condition that still lacks a complete theoretical interpretation, crucial to our understanding of the stellar angular momentum evolution. We couple a rotational evolution m… ▽ More The observed mass-age-rotation relationship in open clusters shows the progressive development of a slow-rotators sequence. The observed clustering on this sequence suggests that it corresponds to some equilibrium or asymptotic condition that still lacks a complete theoretical interpretation, crucial to our understanding of the stellar angular momentum evolution. We couple a rotational evolution model, which takes into account internal differential rotation, with classical and new proposals for the wind braking law, and fit models to the data using a MCMC method. The description of the evolution of the slow-rotators sequence requires taking into account the transfer of angular momentum from the radiative core to the convective envelope; we find that, in the mass range 0.85-1.10 $M_{\odot}$, the core-envelope coupling time-scale for stars in the slow-rotators sequence scales as $M^{-7.28}$. Quasi-solid body rotation is achieved only after 1-2 Gyr, depending on stellar mass, which implies that observing small deviations from the Skumanich law ($P \propto \sqrt{t}$) would require period data of older open clusters than available to date. The observed evolution in the 0.1-2.5 Gyr age range and in the 0.85-1.10 $M_{\odot}$ mass range is best reproduced by assuming an empirical mass dependence of the wind angular momentum loss proportional to the convective turnover time-scale and to the stellar moment of inertia. Period isochrones based on our MCMC fit provide a tool for inferring stellar ages of solar-like main-sequence stars from their mass and rotation period largely independent from the wind braking model adopted. These effectively represent gyro-chronology relationships that take into account the physics of the two-zone model for the stellar angular momentum evolution. △ Less

Submitted 12 October, 2015; v1 submitted 17 June, 2015; originally announced June 2015.

Comments: 16 pages, 8 figures, 4 tables, accepted by A&A

Journal ref: A&A 584, A30 (2015)

Search for Early Gamma-ray Production in Supernovae Located in a Dense Circumstellar Medium with the Fermi LAT

Authors: M. Ackermann, I. Arcavi, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang , et al. (86 additional authors not shown)

Abstract: Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We se… ▽ More Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf). △ Less

Submitted 26 June, 2015; v1 submitted 4 June, 2015; originally announced June 2015.

Comments: Accepted for publication in ApJ. Corresponding author: A. Franckowiak (afrancko@slac.stanford.edu), updated author list and acknowledgements

Journal ref: ApJ, 807, 169 (2015)

Search for Gamma-Ray Emission from DES Dwarf Spheroidal Galaxy Candidates with Fermi-LAT Data

Authors: The Fermi-LAT Collaboration, The DES Collaboration, :, A. Drlica-Wagner, A. Albert, K. Bechtol, M. Wood, L. Strigari, M. Sanchez-Conde, L. Baldini, R. Essig, J. Cohen-Tanugi, B. Anderson, R. Bellazzini, E. D. Bloom, R. Caputo, C. Cecchi, E. Charles, J. Chiang, A. de Angelis, S. Funk, P. Fusco, F. Gargano, N. Giglietto, F. Giordano , et al. (102 additional authors not shown)

Abstract: Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of… ▽ More Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels. △ Less

Submitted 16 August, 2015; v1 submitted 9 March, 2015; originally announced March 2015.

Comments: 9 pages, 4 figures, 1 table. Updated to published version. Readers may be interested in the related work by Ackermann et al. arXiv:1503.02641 (The Fermi-LAT Collaboration), Bechtol, Drlica-Wagner, et al. arXiv:1503.02584 (The DES Collaboration), and Koposov, Belokurov, Torrealba, & Evans arXiv:1503.02079

Report number: FERMILAB-PUB-15-082-AE

Journal ref: ApJ 809 L4 (2015)

The Third Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope

Authors: M. Ackermann, M. Ajello, W. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, J. Gonzalez, R. Bellazzini, E. Bissaldi, R. Blandford, E. Bloom, R. Bonino, E. Bottacini, T. Brandt, J. Bregeon, R. Britto, P. Bruel, R. Buehler, S. Buson, G. Caliandro, R. Cameron, M. Caragiulo, P. Caraveo, J. Casandjian , et al. (118 additional authors not shown)

Abstract: The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10°), a 71% increase over the second catalog based o… ▽ More The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10°), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs. △ Less

Submitted 26 August, 2015; v1 submitted 24 January, 2015; originally announced January 2015.

Comments: Accepted by ApJ. Updated to published version

Stellar Diameters and Temperatures VI. High angular resolution measurements of the transiting exoplanet host stars HD 189733 and HD 209458 and implications for models of cool dwarfs

Authors: Tabetha Boyajian, Kaspar von Braun, Gregory A. Feiden, Daniel Huber, Sarbani Basu, Pierre Demarque, Debra A. Fischer, Gail Schaefer, Andrew W. Mann, Timothy R. White, Vicente Maestro, John Brewer, C. Brooke Lamell, Federico Spada, Mercedes López-Morales, Michael Ireland, Chris Farrington, Gerard T. van Belle, Stephen R. Kane, Jeremy Jones, Theo A. ten Brummelaar, David R. Ciardi, Harold A. McAlister, Stephen Ridgway, P. J. Goldfinger , et al. (2 additional authors not shown)

Abstract: We present direct radii measurements of the well-known transiting exoplanet host stars HD 189733 and HD 209458 using the CHARA Array interferometer. We find the limb-darkened angular diameters to be theta_LD = 0.3848 +/- 0.0055 and 0.2254 +/- 0.0072 milliarcsec for HD 189733 and HD 209458, respectively. HD 189733 and HD 209458 are currently the only two transiting exoplanet systems where detection… ▽ More We present direct radii measurements of the well-known transiting exoplanet host stars HD 189733 and HD 209458 using the CHARA Array interferometer. We find the limb-darkened angular diameters to be theta_LD = 0.3848 +/- 0.0055 and 0.2254 +/- 0.0072 milliarcsec for HD 189733 and HD 209458, respectively. HD 189733 and HD 209458 are currently the only two transiting exoplanet systems where detection of the respective planetary companion's orbital motion from high resolution spectroscopy has revealed absolute masses for both star and planet. We use our new measurements together with the orbital information from radial velocity and photometric time series data, Hipparcos distances, and newly measured bolometric fluxes to determine the stellar effective temperatures (T_eff = 4875 +/- 43, 6093 +/- 103 K), stellar linear radii (R_* = 0.805 +/- 0.016, 1.203 +/- 0.061 R_sun), mean stellar densities (rho_* = 1.62 +/- 0.11, 0.58 +/- 0.14 rho_sun), planetary radii (R_p = 1.216 +/- 0.024, 1.451 +/- 0.074 R_Jup), and mean planetary densities (rho_p = 0.605 +/- 0.029, 0.196 +/- 0.033 rho_Jup) for HD 189733 b and HD 209458 b, respectively. The stellar parameters for HD 209458, a F9 dwarf, are consistent with indirect estimates derived from spectroscopic and evolutionary modeling. However, we find that models are unable to reproduce the observational results for the K2 dwarf, HD 189733. We show that, for stellar evolutionary models to match the observed stellar properties of HD 189733, adjustments lowering the solar-calibrated mixing length parameter from 1.83 to 1.34 need to be employed. △ Less

Submitted 20 November, 2014; originally announced November 2014.