ABSTRACT
Stellar magnetic activity produces time-varying distortions in the photospheric line profiles of solar-type stars. These lead to systematic errors in high-precision radial-velocity ...measurements, which limit efforts to discover and measure the masses of low-mass exoplanets with orbital periods of more than a few tens of days. We present a new data-driven method for separating Doppler shifts of dynamical origin from apparent velocity variations arising from variability-induced changes in the stellar spectrum. We show that the autocorrelation function (ACF) of the cross-correlation function used to measure radial velocities is effectively invariant to translation. By projecting the radial velocities on to a subspace labelled by the observation identifiers and spanned by the amplitude coefficients of the ACF’s principal components, we can isolate and subtract velocity perturbations caused by stellar magnetic activity. We test the method on a 5-yr time sequence of 853 daily 15-min observations of the solar spectrum from the HARPS-N instrument and solar-telescope feed on the 3.58-m Telescopio Nazionale Galileo. After removal of the activity signals, the heliocentric solar velocity residuals are found to be Gaussian and nearly uncorrelated. We inject synthetic low-mass planet signals with amplitude K = 40 cm s−1 into the solar observations at a wide range of orbital periods. Projection into the orthogonal complement of the ACF subspace isolates these signals effectively from solar activity signals. Their semi-amplitudes are recovered with a precision of ∼ 6.6 cm s−1, opening the door to Doppler detection and characterization of terrestrial-mass planets around well-observed, bright main-sequence stars across a wide range of orbital periods.
State-of-the-art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3 yr period during ...the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial-velocity Planet Searcher for the Northern hemisphere (HARPS-N) provides disk-integrated solar spectra, from which we extract RVs and log R HK ′ . The Solar Dynamics Observatory (SDO) provides disk-resolved images of magnetic activity. The Solar Radiation and Climate Experiment (SORCE) provides near-continuous solar photometry, analogous to a Kepler light curve. We verify that the SORCE photometry and HARPS-N log R HK ′ correlate strongly with the SDO-derived magnetic filling factor, while the HARPS-N RV variations do not. To explain this discrepancy, we test existing models of RV variations. We estimate the contributions of the suppression of convective blueshift and the rotational imbalance due to brightness inhomogeneities to the observed HARPS-N RVs. We investigate the time variation of these contributions over several rotation periods, and how these contributions depend on the area of active regions. We find that magnetic active regions smaller than 60 Mm2 do not significantly suppress convective blueshift. Our area-dependent model reduces the amplitude of activity-induced RV variations by a factor of two. The present study highlights the need to identify a proxy that correlates specifically with large, bright magnetic regions on the surfaces of exoplanet-hosting stars.
Abstract
The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets’ masses, obstructing detection of Earth twins. Since 2015 July, ...we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 yr of observation, interpreting them in the context of spatially resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise-limited radial-velocity precision better than 0.43 m s−1 per 5-min observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s−1. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s−1 and an autocorrelation half-life of 15 min. On longer time-scales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s−1 and is correlated with the line asymmetry, but leads it by 3 d. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.
Context.
The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such a high-cadence, long-baseline data set is crucial for understanding spurious ...radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub- m s
−1
systematics that needed to be corrected for.
Aims.
The goals of this manuscript are to (i) present a new data reduction software for HARPS-N, (ii) demonstrate the improvement brought by this new software during the first three years of the HARPS-N solar data set, and (iii) release all the obtained solar products, from extracted spectra to precise radial velocities.
Methods.
To correct for the instrumental systematics observed in the data reduced with the current version of the HARPS-N data reduction software (DRS version 3.7), we adapted the newly available ESPRESSO DRS (version 2.2.3) to HARPS-N and developed new optimised recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS.
Results.
The most significant improvement brought by the new DRS is a strong decrease in the day-to-day radial-velocity scatter, from 1.27 to 1.07 m s
−1
; this is thanks to a more robust method to derive wavelength solutions, but also to the use of calibrations closer in time. The newly derived solar radial-velocities are also better correlated with the chromospheric activity level of the Sun in the long term, with a Pearson correlation coefficient of 0.93 compared to 0.77 before, which is expected from our understanding of stellar signals. Finally, we also discuss how HARPS-N spectral ghosts contaminate the measurement of the calcium activity index, and we present an efficient technique to derive an index free of instrumental systematics.
Conclusions.
This paper presents a new data reduction software for HARPS-N and demonstrates its improvements, mainly in terms of radial-velocity precision, when applied to the first three years of the HARPS-N solar data set. Those newly reduced solar data, representing an unprecedented time series of 34 550 high-resolution spectra and precise radial velocities, are released alongside this paper. Those data are crucial to understand stellar activity signals in solar-type stars further and develop the mitigating techniques that will allow us to detect other Earths.
Context. Understanding stellar activity in solar-type stars is crucial for the physics of stellar atmospheres as well as for ongoing exoplanet programmes. Aims. We aim to test how well we understand ...stellar activity using our own star, the Sun, as a test case. Methods. We performed a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D1 D2, and He I D3) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. We made use of periodogram analyses to study solar rotation, and we used the pool variance technique to study the temporal evolution of active regions. The correlations between the different activity indicators as well as the correlations between activity indexes and the derived parameters from the cross-correlation technique are analysed. We also study the temporal evolution of these correlations and their possible relationship with indicators of inhomogeneities in the solar photosphere like sunspot number or radio flux values. Results. The value of the solar rotation period is found in all the activity indicators, with the only exception being Hδ. The derived values vary from 26.29 days (Hγ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several “bands” of periods around 26 and 30 days. They might be explained by the migration of active regions between the equator and a latitude of ∼30°, spot evolution, or a combination of both effects. A typical lifetime of active regions of approximately ten rotation periods is inferred from the pooled variance diagrams, which is in agreement with previous works. We find that Hα, Hβ, Hγ, Hϵ, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and Hα changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index and radial velocity correlation is also present in the data. Conclusions. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the Hα activity indicators and between the S index and radial velocity variations.
Abstract
Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, ...suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s
−1
. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass (∼0.5 m s
−1
) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus (∼0.1 m s
−1
) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets.
Purpose
Breast cancer is the most common malignancy in women in terms of incidence and mortality. Age is undoubtedly the biggest breast cancer risk factor. In this study we examined clinical, ...histological, and biological characteristics and mortality of breast cancer in elderly women along with their changes with advancing age.
Methods
We reviewed 63 original articles published between 2006 and 2016 concerning women over 70 years with breast cancer.
Results
Compared to patients 70–79 years, patients aged 80 and over had larger tumor size with fewer T1 (42.9% vs 57.7%, p < 0.01) and more T2 lesions (43.5% vs 33.0%, p < 0.01). Lymph nodes and distant metastases were more frequent, with more N + (49.5% vs 44.0%, p < 0.01) and more M1 (8.0% vs 5.9%, p < 0.01). Infiltrating mucinous carcinomas were more frequent (4.3% vs 3.7%, p < 0.01). Tumors had lower grades, with more grade 1 (23.2% vs 19.8%, p = 0.01) and fewer grade 3 (21.5% vs 25.5%, p < 0.01), and were more hormone-sensitive: PR was more often expressed (72.6% vs 67.3%, p < 0.01). Lympho-vascular invasion was less frequent in the 80 years and over (22.9% vs 29.7%, p = 0.01). Breast cancer-specific mortality was higher both at 5 years (25.8% vs 17.2%, p < 0.01) and 10 years (32.7% vs 26.6%, p < 0.01).
Conclusion
Clinico-pathological characteristics, increased incidence, and mortality associated with aging can be explained on one hand by biological changes of the breast such as increased estrogen sensitivity, epithelial cell alterations, immune senescence, and tumor microenvironment modifications. However, sociologic factors such as increased life expectancy, under-treatment, late diagnosis, and insufficient individual screening, are also involved.
Central pattern generators are small networks that contribute to generating animal locomotion. The models used to study gait generation and gait transition mechanisms often require biologically ...accurate neuron and synapse models, with high dimensionality and complex dynamics. Tuning the parameters of these models to elicit network dynamics compatible with gait features is not a trivial task, due to the impossibility of inferring a priori the effects of each parameter on the nonlinear system's emergent dynamics. In this paper we explore the use of global optimization strategies for parameter optimization in multigait central pattern generator (CPG) models with complex cell dynamics and minimal topology. We first consider an existing quadruped CPG model as a test bed for the objective function formulation, then proceed to optimize the parameters of a newly proposed multigait, interlimb hexapod CPG model. We successfully obtain hexapod gaits and prompt gait transitions by varying only control currents, while all CPG parameters, once optimized, are kept fixed. This mechanism of gait transitions is compatible with short-term synaptic plasticity.