Source author record

F. Spada

F. Spada appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2019arXiv

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

(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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

PSR J1906+0722: An Elusive Gamma-ray Pulsar

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.

preprint2015arXiv

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

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).

preprint2015arXiv

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

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\%$.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

The angular momentum transport by unstable toroidal magnetic fields

We demonstrate with a nonlinear MHD code that angular momentum can be transported due to the magnetic instability of toroidal fields under the influence of differential rotation, and that the resulting effective viscosity may be high enough to explain the almost rigid-body rotation observed in radiative stellar cores. Only stationary current-free fields and only those combinations of rotation rates and magnetic field amplitudes which provide maximal numerical values of the viscosity are considered. We find that the dimensionless ratio of the effective over molecular viscosity, $ν_T/ν$;, linearly grows with the Reynolds number of the rotating fluid multiplied with the square-root of the magnetic Prandtl number - which is of order unity for the considered red sub-giant KIC 7341231. For the considered interval of magnetic Reynolds numbers - which is restricted by numerical constraints of the nonlinear MHD code - there is a remarkable influence of the magnetic Prandtl number on the relative importance of the contributions of the Reynolds stress and the Maxwell stress to the total viscosity, which is magnetically dominated only for Pm $\gtrsim$ 0.5. We also find that the magnetized plasma behaves as a non-Newtonian fluid, i.e. the resulting effective viscosity depends on the shear in the rotation law. The decay time of the differential rotation thus depends on its shear and becomes longer and longer during the spin-down of a stellar core.

preprint2014arXiv

Upgrading the Solar-Stellar Connection: News about activity in Cool Stars

In this splinter session, ten speakers presented results on solar and stellar activity and how the two fields are connected. This was followed by a lively discussion and supplemented by short, one-minute highlight talks. The talks presented new theoretical and observational results on mass accretion on the Sun, the activity rate of flare stars, the evolution of the stellar magnetic field on time scales of a single cycle and over the lifetime of a star, and two different approaches to model the radial-velocity jitter in cool stars that is due to the granulation on the surface. Talks and discussion showed how much the interpretation of stellar activity data relies on the sun and how the large number of objects available in stellar studies can extend the parameter range of activity models.

preprint2013arXiv

The radius discrepancy in low mass stars: single vs. binaries

A long-standing issue in the theory of low mass stars is the discrepancy between predicted and observed radii and effective temperatures. In spite of the increasing availability of very precise radius determinations from eclipsing binaries and interferometric measurements of radii of single stars, there is no unanimous consensus on the extent (or even the existence) of the discrepancy and on its connection with other stellar properties (e.g. metallicity, magnetic activity). We investigate the radius discrepancy phenomenon using the best data currently available (accuracy about 5%). We have constructed a grid of stellar models covering the entire range of low mass stars (0.1-1.25 M_sun) and various choices of the metallicity and of the mixing length parameter α. We used an improved version of the Yale Rotational stellar Evolution Code (YREC), implementing surface boundary conditions based on the most up-to-date PHOENIX atmosphere models. Our models are in good agreement with others in the literature and improve and extend the low mass end of the Yale-Yonsei isochrones. Our calculations include rotation-related quantities, such as moments of inertia and convective turnover time scales, useful in studies of magnetic activity and rotational evolution of solar-like stars. Consistently with previous works, we find that both binaries and single stars have radii inflated by about 3% with respect to the theoretical models; among binaries, the components of short orbital period systems are found to be the most deviant. We conclude that both binaries and single stars are comparably affected by the radius discrepancy phenomenon.

preprint2011arXiv

Modelling the rotational evolution of solar-like stars: the rotational coupling timescale

We investigate the rotational evolution of solar-like stars with a focus on the internal angular momentum transport processes. The double zone model, in which the star's radiative core and convective envelope are assumed to rotate as solid bodies, is used to test simple relationships between the core-envelope coupling timescale, tau_c, and rotational properties, like the envelope angular velocity or the differential rotation at the core-envelope interface. The trial relationships are tested by fitting the model parameters to available observations via a Monte Carlo Markov Chain method. The synthetic distributions are tested for compatibility with their observational counterparts by means of the standard Kolmogorov-Smirnov (KS) test. A power-law dependence of tau_c on the inner differential rotation leads to a more satisfactory agreement with observations than a two-value prescription for tau_c, which would imply a dichotomy between the initially slow (P_rot > 3 d) and fast (P_rot < 3 d) rotators. However, we find it impossible to reconcile the high fraction of fast rotators in alpha Per with the rotation period distributions in stellar systems at earlier and later evolutionary stages. This could be explained by local environmental effects (e.g. early removal of circumstellar discs due to UV radiation and winds from nearby high-mass stars) or by observational biases. The low KS probability that the synthetic and observed distributions are not incompatible, found in some cases, may be due to over-simplified assumptions of the double zone model, but the large relative uncertainties in the age determination of very young clusters and associations are expected to play a relevant role. Other possible limitations and uncertainties are discussed.

preprint2009arXiv

A semi-analytic approach to angular momentum transport in stellar radiative interiors

We address the problem of angular momentum transport in stellar radiative interiors with a novel semi-analytic spectral technique, using an eigenfunction series expansion, that can be used to derive benchmark solutions in hydromagnetic regimes with very high Reynolds number (10^7 - 10^8). The error arising from the truncation of the series is evaluated analytically. The main simplifying assumptions are the neglect of meridional circulation and of non-axisymmetric magnetic fields. The advantages of our approach are shown by applying it to a spin-down model for a 1 M_sun main-sequence star. The evolution of the coupling between core and envelope is investigated for different values of the viscosity and different geometries and values of the poloidal field. We confirm that a viscosity enhancement by 10^4 with respect to the molecular value is required to attain a rigid rotation in the core of the Sun within its present age. We suggest that a quadrupolar poloidal field may explain the short coupling time-scale needed to model the observed rotational evolution of fast rotators on the ZAMS, while a dipolar geometry is indicated in the case of slow rotators. Our novel semi-analytic spectral method provides a conceptually simple and rigorous treatment of a classic MHD problem and allows us to explore the influence of various parameters on the rotational history of radiative interiors.