# On the Metamorphosis of the Bailey Diagram for RR Lyrae Stars

Guillaume Bono, V. Braga, J. Crestani, M. Fabrizio, C. Sneden, M. Marconi, G. Preston, J. Mullen, C. Gilligan, G. Fiorentino, et al.

► To cite this version:

Guillaume Bono, V. Braga, J. Crestani, M. Fabrizio, C. Sneden, et al. On the Metamorphosis of the Bailey Diagram for RR Lyrae Stars. The Astrophysical Journal Letters, 2020, 896 (1), pp.L15. 10.3847/2041-8213/ab9538. hal-03047116

HAL Id: hal-03047116

https://hal.science/hal-03047116v1

Submitted on 1 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# On the Metamorphosis of the Bailey Diagram for RR Lyrae Stars

G. Bono<sup>1,2</sup>①, V. F. Braga<sup>2,3</sup>①, J. Crestani<sup>1,2,4</sup>, M. Fabrizio<sup>2,3</sup>①, C. Sneden<sup>5</sup>①, M. Marconi<sup>6</sup>①, G. W. Preston<sup>7</sup>, J. P. Mullen<sup>8</sup>, C. K. Gilligan<sup>9</sup>, G. Fiorentino<sup>2</sup>①, A. Pietrinferni<sup>10</sup>, G. Altavilla<sup>2,3</sup>, R. Buonanno<sup>10</sup>, B. Chaboyer<sup>9</sup>, R. da Silva<sup>2,3</sup>, M. Dall'Ora<sup>6</sup>, S. Degl'Innocenti<sup>11,12</sup>, E. Di Carlo<sup>10</sup>, I. Ferraro<sup>2</sup>, E. K. Grebel<sup>13</sup>, G. Iannicola<sup>2</sup>, L. Inno<sup>6,14</sup>, V. Kovtyukh<sup>15,16</sup>, A. Kunder<sup>17</sup>①, B. Lemasle<sup>13</sup>, M. Marengo<sup>8</sup>①, S. Marinoni<sup>2,3</sup>, P. M. Marrese<sup>2,3</sup>, C. E. Martínez-Vázquez<sup>18</sup>①, N. Matsunaga<sup>19</sup>, M. Monelli<sup>20</sup>①, J. Neeley<sup>21</sup>①, M. Nonino<sup>22</sup>, P. G. Prada Moroni<sup>11,12</sup>, Z. Prudil<sup>13</sup>, P. B. Stetson<sup>23</sup>①, F. Thévenin<sup>24</sup>, E. Tognelli<sup>12</sup>, E. Valenti<sup>25</sup>①, and A. R. Walker<sup>18</sup>①

## Abstract

We collected over 6000 high-resolution spectra of four dozen field RR Lyrae (RRL) variables pulsating either in the fundamental (39 RRab) or in the first overtone (9 RRc) mode. We measured radial velocities (RVs) of four strong metallic and four Balmer lines along the entire pulsational cycle and derived RV amplitudes with accuracies better than 1–2 km s⁻<fnref n="1" />. The new amplitudes were combined with literature data for 23 RRab and 3 RRc stars (total sample of 74 RRLs), which allowed us to investigate the variation of the Bailey diagram (photometric amplitude versus period) when moving from optical to mid-infrared bands and to recast the Bailey diagram in terms of RV amplitudes. We found that RV amplitudes for RRab are minimally affected by nonlinear phenomena (shocks) and multiperiodicity (Blazhko effect). The RV slope (logP–A(V<sub>r</sub>)) when compared with the visual slope (logP–A(V)) is shallower, and the dispersion, at fixed period, decreases by a factor of two. We constructed homogeneous sets of horizontal branch evolutionary models and nonlinear, convective pulsation models of RRLs to constrain the impact of evolutionary effects on their pulsation properties. Evolution causes, on the Bailey diagram based on RV amplitudes, a modest variation in pulsation period and a large dispersion in amplitude. The broad dispersion in period of the Bailey diagram is mainly caused by variation in RRL intrinsic parameters (stellar mass, chemical composition). Empirical evidence indicates that RV amplitudes are an optimal diagnostic for tracing the mean effective temperature across the RRab instability strip.

*Unified Astronomy Thesaurus concepts: Horizontal branch stars (746); Stellar oscillations (1617); RR Lyrae variable stars (1410); Radial velocity (1332)*

## 1. Introduction

The RR Lyrae (RRLs) are a very popular group of radial variables, because they can be easily identified thanks to the coupling between the pulsation periods (0.25 ≲ P ≲ 1.0 day) and the shapes of their light curves (da Costa et al. 2010; Monelli et al. 2017).

RRLs pulsate in three different flavors: fundamental (RRab), first overtone (RRc), and mixed mode (pulsating simultaneously in the fundamental and in the first overtone). The

Bailey diagram (Bailey et al. 1919)—luminosity amplitude versus logarithmic period—is a fundamental plane to identify the RRL pulsation mode. In this plane RRc variables cover the short-period, low-amplitude domain, and the shape of their amplitude-period relationship displays either a “bell” or a “hairpin” shape (Bono et al. 2011; Kunder et al. 2013; Fiorentino et al. 2017). The RRab variables cover the long-period, large-amplitude area and display a well-defined decrease in luminosity amplitude as period increases.

<footnote n="1"><sup>1</sup> Dipartimento di Fisica, Università di Roma Tor Vergata, via della Ricerca Scientifica 1, I-00133 Roma, Italy</footnote>

<footnote n="2"><sup>2</sup> INAF-Osservatorio Astronomico di Roma, via Frascati 33, I-00078 Monte Porzio Catone, Italy</footnote>

<footnote n="3"><sup>3</sup> Space Science Data Center, via del Politecnico snc, I-00133 Roma, Italy</footnote>

<footnote n="4"><sup>4</sup> Departamento de Astronomia, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves 6500, Porto Alegre 91501-970, Brazil</footnote>

<footnote n="5"><sup>5</sup> Department of Astronomy and McDonald Observatory, The University of Texas, Austin, TX 78712, USA</footnote>

<footnote n="6"><sup>6</sup> INAF-Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, I-80131 Napoli, Italy</footnote>

<footnote n="7"><sup>7</sup> The Observatories of the Carnegie Institution for Science, 813 Santa Barbara Street, Pasadena, CA 91101, USA</footnote>

<footnote n="8"><sup>8</sup> Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA</footnote>

<footnote n="9"><sup>9</sup> Department of Physics and Astronomy, Dartmouth College, Hanover, NH 03755, USA</footnote>

<footnote n="10"><sup>10</sup> INAF-Osservatorio Astronomico d'Abruzzo, Via Mentore Maggini snc, Loc. Collurania, I-64100 Teramo, Italy</footnote>

<footnote n="11"><sup>11</sup> INFN, Sezione di Pisa, Largo Pontecorvo 3, I-56127 Pisa, Italy</footnote>

<footnote n="12"><sup>12</sup> Dipartimento di Fisica "Enrico Fermi", Università di Pisa, Largo Pontecorvo 3, I-56127 Pisa, Italy</footnote>

<footnote n="13"><sup>13</sup> Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstr. 12-14, D-69120 Heidelberg, Germany</footnote>

<footnote n="14"><sup>14</sup> Università Parthenope di Napoli, Science and Technology Department, CDN IC4, I-80143 Naples, Italy</footnote>

<footnote n="15"><sup>15</sup> Astronomical Observatory, Odessa National University, Shevchenko Park, UA-65014 Odessa, Ukraine</footnote>

<footnote n="16"><sup>16</sup> Newton Institute of Chile, Odessa branch, Shevchenko Park, UA-65014 Odessa, Ukraine</footnote>

<footnote n="17"><sup>17</sup> Saint Martin's University, 5000 Abbey Way SE, Lacey, WA 98503, USA</footnote>

<footnote n="18"><sup>18</sup> Cerro Tololo Inter-American Observatory, NSF's National Optical-Infrared Astronomy Research Laboratory, Casilla 603, La Serena, Chile
<sup>19</sup> Department of Astronomy, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan</footnote>

<footnote n="20"><sup>20</sup> Instituto de Astrofísica de Canarias, Calle Via Lactea s/n, E-38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="21"><sup>21</sup> Department of Physics, Florida Atlantic University, 777 Glades Road, Boca Raton, FL 33431 USA</footnote>

<footnote n="22"><sup>22</sup> INAF-Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, I-34143 Trieste, Italy</footnote>

<footnote n="23"><sup>23</sup> Herzberg Astronomy and Astrophysics, National Research Council, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada</footnote>

<footnote n="24"><sup>24</sup> Université de Nice Sophia-antipolis, CNRS, Observatoire de la Côte d'Azur, Laboratoire Lagrange, BP 4229, F-06304 Nice, France</footnote>

<footnote n="25"><sup>25</sup> European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei Munchen, Germany</footnote>

<footnote>*Received 2020 April 21; revised 2020 May 20; accepted 2020 May 21; published 2020 June 11*</footnote>
The Bailey diagram has also been the crossroad of a long-
standing astrophysical problem, the so-called Oosterhoff dichot-
omy, identified over 80 yr ago (Oosterhoff 1939). RRLs in
globular clusters (GCs) can be split according to their mean
periods in two groups: the Oosterhoff type I (OoI), with
⟨P<sub>ab</sub>⟩ ~ 0.56, ⟨P<sub>c</sub>⟩ ~ 0.31 days, and the Oosterhoff type II
(OoII), with ⟨P<sub>ab</sub>⟩ ~ 0.66, ⟨P<sub>c</sub>⟩ ~ 0.36 days. Subsequent inves-
tigations identified a similar dichotomic distribution for the RRc,
albeit with a difference of only ~0.05 days. Spectroscopic studies
demonstrated that OoI GCs are more metal-rich and have a
smaller ratio of RRc to the total number of RRLs than OoII GCs.
More recently, Fabrizio et al. (2019) found a linear correlation
between mean period and metallicity, thus, suggesting that this
dichotomy is a natural consequence of the lack of Galactic,
metal-intermediate GCs rich in RRLs.

RRL luminosity amplitudes steadily decrease when shifting
to longer wavelengths. However, the amplitude attains an
almost constant value at wavelengths longer than the *K* band.
This behavior is caused by the fact that luminosity amplitudes
in the optical bands are mainly driven by temperature
variations, while in the near-infrared (NIR)/mid-infrared
(MIR), they are mainly driven by radius variations (Madore
et al. 2013; Neeley et al. 2017).

The Bailey diagram relies on two observables that are
independent of distance and reddening. The RRL distribution
in the Bailey diagram depends on their intrinsic parameters
(stellar mass, luminosity, effective temperature) and chemical
composition, but it is also affected by degeneracies caused by
three different phenomena. (i) Secondary modulations. Empiri-
cal evidence suggests that a fraction of RRLs are affected by a
secondary modulation (Blazhko effect, Blažko 1907) with a
timescale ranging from tens of days to a few years. This
fraction in large samples (in the Galactic bulge) ranges from
~6% for RRc to ~40% for RRab (Prudil & Skarka 2017;
Netzel et al. 2018). This means that the position of an RRL in
the Bailey diagram can change by several tenths of a magnitude
along the Blazhko cycle. The Blazhko effect partially explains
why the RRLs display, at fixed pulsation period and chemical
composition, a well-defined intrinsic spread in luminosity
amplitude. (ii) Evolutionary phase. RRLs are low-mass helium-
burning stars, and during this evolutionary phase they evolve
off the zero-age horizontal branch (ZAHB). In their evolution
they experience changes both in effective temperature and in
surface gravity. These changes cause a variation both in
pulsation period and in luminosity amplitude. We still lack
empirical evidence concerning the variation of the position in
the Bailey diagram directly caused by evolutionary effects.
(iii) Nonlinear phenomena. Theory and observations indicate
that RRLs experience violent shocks approaching maximum
compression (Preston et al. 2019). Shocks are strongly affected
by the variation of the efficiency of convective motions (Bono
& Stellingwerf 1994). This is the reason why RRab amplitudes
display a steady decrease close to the red edge of the instability
strip (IS; Bono et al. 1997; Braga et al. 2016).

The main aim of this investigation is to study the physical
mechanisms affecting the distribution of RRLs in the Bailey
diagram.

## 2. Spectroscopic and Photometric Data Sets

The largest and most homogeneous subsample of high-resolution (HR) spectra was collected with the Las Campanas Observatory du Pont echelle spectrograph. It includes 6206 HR

(R ≡ λ/δλ = 35,000) spectra of 192 RRLs. Their typical
signal-to-noise ratio (S/N), in the optical regime (λ = 5100 Å),
is ~40.

These spectra were augmented with 272 HR (R = 34,540–107,200, S/N ~ 20; 68 RRLs) UVES@VLT spectra; 41 medium-resolution (R = 4,300–18,000, S/N ~ 44; 18 RRLs) X-Shooter@VLT spectra; 166 HR (R = 115,000, S/N ~ 10; 5 RRLs) HARPS@3.6 m spectra, and 55 HR (R = 48,000, S/N ~ 4; 3 RRLs) FEROS@2.2 m (ESO, Chile) spectra.

Additionally we included HR spectra from SES@STELLA (15, R = 55,000, S/N ~ 35; 1 RRL), HRS@SALT (82, R ~ 40,000, S/N ~ 50; 70 RRLs), and HARPS-N@TNG (10, R = 115,000, S/N ~ 40; 4 RRLs).

In total we collected 6572 spectra for 316 RRLs. Among them, 48 RRLs have a good/optimal coverage of the pulsation cycle (6226 spectra). We label these variables as RV calibrators. These data were complemented with RV curves for 26 field and cluster RRLs available in the literature from Baade-Wesselink (BW) analyses (Storm et al. 1994; Bono et al. 2003). We ended up with a sample of 74 (62 RRab, 12 RRc) RV calibrators (see Table 1).

The periods from the literature are usually precise enough to
provide good phasing of the RV curves. However, we derived
the pulsation period from RV measurements for 2 RRab and 1
RRc. We adopted Gaia *G*-band amplitudes (Gaia Collaboration
et al. 2018; Marrese et al. 2019) transformed into the *V* band by
using the equations provided by Evans et al. (2018). We also
used *V*-band photometry from ASAS (Pojmanski 2002),
ASAS-SN (Shappee et al. 2014; Jayasinghe et al. 2018), and
from Catalina (Drake et al. 2013a, 2013b; Torrealba et al.
2015).

The optical photometry was combined with MIR W1-band
photometry collected with NEOWISE (Mainzer et al. 2011).
The MIR light curves have very good phase coverage (~120
phase points), and good photometric precision (σ<sub>W1</sub> ≤
0.05 mag). Optical and MIR light curves were fitted with
analytical functions (splines, PLOESS, Fourier parameters)
following the same approach by Braga et al. (2018). The
photometric amplitudes are based on the analytical fits.

3. Radial Velocity Measurements

To measure RVs we chose three Fe I lines, from the multiplet 43, at λ = 4045.81, 4063.59, 4071.74 Å and the Sr II line at λ = 4077.71 Å. These metallic lines were strong enough to be easily identified in most spectra even at very low S/N or very high effective temperatures, enabling us to trace the whole pulsation cycle using individual spectra, i.e., without coadding spectra. We also used four Balmer lines (Hα, Hβ, Hγ, Hδ, λ = 6562.80, 4861.36, 4340.46, 4101.74 Å). The RVs were measured following the same approach employed by Fabrizio et al. (2011).

Figure 1 shows RV curves for the RRab, V Ind, and the RRc, MT Tel. The RV measurements plotted in this figure display the expected difference between metallic and Balmer lines (Preston 1959); the RV amplitudes for the latter ones are systematically larger (see also Figure 9 of Chadid et al. 2017 for montages of RV and photometric phase plots). The difference is caused by the physical conditions at which the different lines form: the smaller the optical depth, the larger the RV amplitude. The difference among Hα and metallic lines ranges from 50% for MT Tel to almost a factor of two for
<table><caption>Table 1
Pulsation Properties for the Current RRLs</caption><thead><tr><th>Gaia DR2 ID</th><th>Classical ID</th><th>Type<sup>a</sup></th><th>Period (days)</th><th>A(V) (mag)</th><th>A(W1) (mag)</th><th>A(Vr)<sub>re</sub> (km s<sup>-1</sup>)</th><th>Flag<sup>b</sup></th><th>A(Vr)<sub>Hα</sub> (km s<sup>-1</sup>)</th><th>A(Vr)<sub>Hβ</sub> (km s<sup>-1</sup>)</th><th>A(Vr)<sub>Hγ</sub> (km s<sup>-1</sup>)</th><th>A(Vr)<sub>Hδ</sub> (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>15489408711727488</td><td>X Ari</td><td>0</td><td>0.65117537</td><td>0.96</td><td>0.32</td><td>58.2 ± 2.4</td><td>0</td><td>108.2 ± 4.4</td><td>85.1 ± 3.3</td><td>76.4 ± 3.0</td><td>66.5 ± 2.9</td></tr><tr><td>234108363683247616</td><td>AR Per</td><td>0</td><td>0.42556048</td><td>0.88</td><td>0.27</td><td>65.7 ± 5.1</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>630421935431871232</td><td>RR Leo</td><td>0</td><td>0.45240000</td><td>1.25</td><td>0.32</td><td>63.3 ± 6.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>886793515494085248</td><td>RR Gem</td><td>0</td><td>0.39730000</td><td>1.17</td><td>0.28</td><td>70.0 ± 11.3</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1009665142487836032</td><td>TT Lyn</td><td>0</td><td>0.59744301</td><td>0.67</td><td>0.25</td><td>50.5 ± 3.1</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1058066262817534336</td><td>SU Dra</td><td>0</td><td>0.66041178</td><td>0.98</td><td>0.28</td><td>57.6 ± 3.5</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1191510003353849472</td><td>AN Ser</td><td>0</td><td>0.52207295</td><td>0.99</td><td>0.29</td><td>60.2 ± 3.8</td><td>0</td><td>94.5 ± 6.4</td><td>73.9 ± 5.2</td><td>70.5 ± 4.9</td><td>57.9 ± 3.8</td></tr><tr><td>1286188056265485952</td><td>RS Boo</td><td>0</td><td>0.37736549</td><td>1.04</td><td>0.25</td><td>67.1 ± 3.2</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1360405567883886720</td><td>CI<sup>*</sup> NGC 6341 SAW V1</td><td>0</td><td>0.70279828</td><td>0.87</td><td>0.26</td><td>61.2 ± 4.9</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1360408905076366208</td><td>CSS J171712.0+431225</td><td>0</td><td>0.63746100</td><td>0.76</td><td>0.26</td><td>61.4 ± 3.9</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1492230556717187456</td><td>TV Boo</td><td>1</td><td>0.31256000</td><td>0.63</td><td>0.12</td><td>35.0 ± 5.6</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1683444631037761024</td><td>SW Dra</td><td>0</td><td>0.56967009</td><td>0.94</td><td>0.28</td><td>62.6 ± 4.5</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>1760981190300823808</td><td>DX Del</td><td>0</td><td>0.47261773</td><td>0.70</td><td>0.24</td><td>54.0 ± 3.2</td><td>0</td><td>94.1 ± 6.1</td><td>71.3 ± 4.1</td><td>66.3 ± 4.1</td><td>58.3 ± 3.6</td></tr><tr><td>1793460115244988800</td><td>AV Peg</td><td>0</td><td>0.39038090</td><td>1.01</td><td>0.28</td><td>63.4 ± 2.8</td><td>0</td><td>92.6 ± 4.4</td><td>74.4 ± 3.4</td><td>66.2 ± 3.3</td><td>59.8 ± 3.0</td></tr><tr><td>2381771781829913984</td><td>DN Aqr</td><td>0</td><td>0.63376712</td><td>0.74</td><td>0.24</td><td>54.6 ± 2.1</td><td>0</td><td>104.3 ± 4.0</td><td>79.6 ± 3.4</td><td>79.6 ± 3.5</td><td>63.8 ± 3.0</td></tr><tr><td>2414817603803476864</td><td>UU Cet</td><td>0</td><td>0.60610163</td><td>0.60</td><td>0.23</td><td>53.0 ± 3.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>2558296724402139392</td><td>RR Cet</td><td>0</td><td>0.55302505</td><td>0.90</td><td>0.28</td><td>60.1 ± 3.1</td><td>0</td><td>106.2 ± 5.5</td><td>84.1 ± 4.3</td><td>81.2 ± 4.3</td><td>68.7 ± 3.6</td></tr><tr><td>2689556491246048896</td><td>SW Aqr</td><td>0</td><td>0.45930070</td><td>1.24</td><td>0.31</td><td>61.9 ± 1.8</td><td>0</td><td>105.1 ± 3.3</td><td>82.7 ± 2.6</td><td>79.5 ± 2.6</td><td>64.8 ± 2.3</td></tr><tr><td>2720896455287475584</td><td>DH Peg</td><td>1</td><td>0.25551624</td><td>0.50</td><td>0.10</td><td>25.4 ± 1.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>2857456207478683776</td><td>SW And</td><td>0</td><td>0.44230000</td><td>0.92</td><td>0.27</td><td>62.5 ± 4.6</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>2981136563930816128</td><td>RX Eri</td><td>0</td><td>0.58725159</td><td>0.87</td><td>0.28</td><td>60.3 ± 4.5</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>3479598373678136832</td><td>DT Hya</td><td>0</td><td>0.56798140</td><td>0.95</td><td>0.27</td><td>61.8 ± 3.1</td><td>1</td><td>109.9 ± 5.7</td><td>83.7 ± 4.3</td><td>78.1 ± 4.2</td><td>70.9 ± 3.7</td></tr><tr><td>3535368455297644928</td><td>ASAS J110522-2641.0</td><td>1</td><td>0.29445920</td><td>0.33</td><td>...</td><td>14.3 ± 1.7</td><td>0</td><td>16.4 ± 2.2</td><td>...</td><td>...</td><td>...</td></tr><tr><td>3546458301374134528</td><td>W Crt</td><td>0</td><td>0.41201190</td><td>1.36</td><td>0.32</td><td>66.1 ± 2.5</td><td>0</td><td>101.9 ± 4.4</td><td>83.9 ± 3.7</td><td>77.1 ± 3.3</td><td>66.6 ± 2.9</td></tr><tr><td>3604450388616968576</td><td>AM Vir</td><td>0</td><td>0.61510000</td><td>0.66</td><td>0.25</td><td>59.0 ± 2.5</td><td>0</td><td>98.1 ± 4.2</td><td>81.1 ± 3.6</td><td>77.3 ± 3.5</td><td>63.2 ± 2.8</td></tr><tr><td>3626569264033312896</td><td>AS Vir</td><td>0</td><td>0.55341000</td><td>0.63</td><td>0.22</td><td>59.1 ± 2.0</td><td>0</td><td>85.2 ± 2.9</td><td>77.2 ± 2.7</td><td>74.6 ± 2.8</td><td>65.7 ± 2.5</td></tr><tr><td>3652665558338018048</td><td>ST Vir</td><td>0</td><td>0.41080754</td><td>1.22</td><td>0.33</td><td>60.3 ± 3.4</td><td>0</td><td>97.9 ± 5.8</td><td>72.2 ± 4.4</td><td>78.5 ± 4.8</td><td>61.2 ± 4.0</td></tr><tr><td>3698725337375614464</td><td>UU Vir</td><td>0</td><td>0.47560267</td><td>1.10</td><td>0.30</td><td>67.8 ± 3.1</td><td>0</td><td>114.4 ± 7.2</td><td>91.7 ± 4.6</td><td>83.6 ± 4.3</td><td>74.9 ± 3.6</td></tr><tr><td>3797319369672686592</td><td>SS Leo</td><td>0</td><td>0.62632619</td><td>1.09</td><td>0.30</td><td>62.5 ± 2.0</td><td>0</td><td>97.2 ± 4.0</td><td>83.8 ± 3.0</td><td>80.9 ± 3.1</td><td>71.4 ± 2.8</td></tr><tr><td>3846786226007324160</td><td>T Sex</td><td>1</td><td>0.32468493</td><td>0.42</td><td>0.08</td><td>27.0 ± 2.3</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>3915998558830693888</td><td>ST Leo</td><td>0</td><td>0.47797595</td><td>1.23</td><td>0.31</td><td>64.0 ± 4.3</td><td>0</td><td>112.4 ± 7.6</td><td>82.9 ± 5.4</td><td>85.0 ± 6.1</td><td>65.3 ± 4.9</td></tr><tr><td>4022618712476736896</td><td>TU UMa</td><td>0</td><td>0.55690000</td><td>0.96</td><td>0.28</td><td>60.6 ± 4.3</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>4055098870077726976</td><td>V0494 Sco</td><td>0*</td><td>0.42727100</td><td>0.87</td><td>0.20</td><td>44.5 ± 1.7</td><td>0</td><td>68.6 ± 3.1</td><td>53.5 ± 2.5</td><td>53.7 ± 2.4</td><td>45.7 ± 2.4</td></tr><tr><td>4352084489819078784</td><td>V0445 Oph</td><td>0</td><td>0.39702600</td><td>0.77</td><td>0.26</td><td>57.4 ± 4.6</td><td>0</td><td>91.1 ± 7.1</td><td>68.6 ± 5.1</td><td>61.6 ± 5.1</td><td>58.8 ± 4.8</td></tr><tr><td>4417888542753226112</td><td>VY Ser</td><td>0</td><td>0.71410000</td><td>0.65</td><td>0.25</td><td>53.3 ± 4.3</td><td>0</td><td>90.6 ± 6.3</td><td>76.3 ± 11.2</td><td>74.2 ± 6.5</td><td>69.2 ± 6.0</td></tr><tr><td>4421571803630954752</td><td>CSS J151841.9+020232</td><td>0</td><td>0.54625080</td><td>0.65</td><td>0.26</td><td>51.8 ± 3.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>4421571803630958848</td><td>CI<sup>*</sup> NGC 5904 SAW V28</td><td>0</td><td>0.54383658</td><td>0.71</td><td>...</td><td>60.9 ± 4.7</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>4433070255716036864</td><td>ASAS J162158+0244.5</td><td>1</td><td>0.32362580</td><td>0.41</td><td>0.09</td><td>20.2 ± 2.3</td><td>0</td><td>29.2 ± 3.2</td><td>29.0 ± 3.8</td><td>...</td><td>...</td></tr><tr><td>4454183799545435008</td><td>AT Ser</td><td>0</td><td>0.74655408</td><td>0.89</td><td>0.26</td><td>56.8 ± 4.1</td><td>0</td><td>99.8 ± 7.8</td><td>77.2 ± 5.6</td><td>77.0 ± 5.7</td><td>70.0 ± 5.4</td></tr><tr><td>4467433017735910912</td><td>VX Scl</td><td>0</td><td>0.45518030</td><td>1.28</td><td>0.33</td><td>60.6 ± 3.9</td><td>0</td><td>105.9 ± 7.1</td><td>81.0 ± 5.3</td><td>74.7 ± 5.0</td><td>64.9 ± 4.5</td></tr><tr><td>4596935593202765184</td><td>TW Her</td><td>0</td><td>0.39959577</td><td>1.33</td><td>0.33</td><td>68.3 ± 4.7</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>4689637956899105792</td><td>CI<sup>*</sup> NGC 104 SAW V9</td><td>0</td><td>0.73720048</td><td>1.06</td><td>...</td><td>57.0 ± 5.3</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>4709830423483623808</td><td>W Tuc</td><td>0</td><td>0.64224028</td><td>1.14</td><td>0.31</td><td>63.7 ± 2.9</td><td>0</td><td>114.0 ± 5.8</td><td>79.9 ± 4.1</td><td>80.0 ± 4.0</td><td>66.2 ± 3.2</td></tr><tr><td>4860671839583430912</td><td>SX For</td><td>0</td><td>0.60534530</td><td>0.66</td><td>...</td><td>51.9 ± 2.9</td><td>0</td><td>93.5 ± 5.9</td><td>70.6 ± 4.3</td><td>66.6 ± 3.8</td><td>61.4 ± 3.7</td></tr><tr><td>4947090013255935616</td><td>CS Eri</td><td>1</td><td>0.31133020</td><td>0.52</td><td>...</td><td>28.9 ± 4.0</td><td>0</td><td>48.6 ± 3.9</td><td>34.7 ± 3.0</td><td>33.8 ± 2.9</td><td>28.2 ± 2.6</td></tr><tr><td>5022411786734718208</td><td>SV Scl</td><td>1</td><td>0.37735860</td><td>0.51</td><td>...</td><td>26.1 ± 2.1</td><td>0</td><td>42.2 ± 3.2</td><td>30.6 ± 2.6</td><td>29.4 ± 2.6</td><td>24.0 ± 2.6</td></tr><tr><td>5151789464548893184</td><td>RV Cet</td><td>0*</td><td>0.61802300</td><td>0.57</td><td>...</td><td>49.8 ± 2.3</td><td>0</td><td>89.3 ± 4.3</td><td>70.1 ± 3.3</td><td>65.6 ± 3.0</td><td>57.9 ± 2.7</td></tr></tbody></table>
<table><caption>**Table 1**
(Continued)</caption><thead><tr><th>Gaia DR2 ID</th><th>Classical ID</th><th>Type<sup>a</sup></th><th>Period (days)</th><th>A(V) (mag)</th><th>A(W1) (mag)</th><th>A(Vr)<sub>Fe</sub> (km s<sup>-1</sup>)</th><th>Flag<sup>b</sup></th><th>A(Vr)<sub>Hα</sub> (km s<sup>-1</sup>)</th><th>A(Vr)<sub>Hβ</sub> (km s<sup>-1</sup>)</th><th>A(Vr)<sub>Hγ</sub> (km s<sup>-1</sup>)</th><th>A(Vr)<sub>Hδ</sub> (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>5412243359495900928</td><td>CD Vel</td><td>0</td><td>0.57350788</td><td>0.86</td><td>0.27</td><td>52.5 ± 1.8</td><td>0</td><td>96.6 ± 4.1</td><td>74.0 ± 2.0</td><td>69.8 ± 2.8</td><td>60.7 ± 2.5</td></tr><tr><td>5461994297841116160</td><td>WY Ant</td><td>0</td><td>0.57434364</td><td>0.92</td><td>0.27</td><td>60.6 ± 2.9</td><td>0</td><td>106.5 ± 5.6</td><td>83.3 ± 4.4</td><td>77.1 ± 4.2</td><td>66.5 ± 3.7</td></tr><tr><td>5510293236607430656</td><td>HH Pup</td><td>0</td><td>0.39081190</td><td>1.22</td><td>0.32</td><td>67.9 ± 2.6</td><td>0</td><td>109.3 ± 4.3</td><td>81.4 ± 3.3</td><td>78.1 ± 3.4</td><td>68.9 ± 2.9</td></tr><tr><td>5707380936404638336</td><td>BB Pup</td><td>0</td><td>0.48055043</td><td>0.88</td><td>0.29</td><td>64.8 ± 4.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>5768557209320424320</td><td>UV Oct</td><td>0*</td><td>0.54257500</td><td>0.80</td><td>0.27</td><td>57.1 ± 1.1</td><td>0</td><td>90.3 ± 2.0</td><td>73.5 ± 1.6</td><td>73.0 ± 1.8</td><td>62.3 ± 1.7</td></tr><tr><td>5769986338215537280</td><td>RV Oct</td><td>0</td><td>0.57117800</td><td>1.05</td><td>0.30</td><td>65.1 ± 2.3</td><td>0</td><td>103.9 ± 3.6</td><td>82.7 ± 2.9</td><td>81.9 ± 3.0</td><td>70.1 ± 2.7</td></tr><tr><td>5773390391856998656</td><td>XZ Aps</td><td>0</td><td>0.58726739</td><td>1.15</td><td>0.29</td><td>63.0 ± 1.8</td><td>0</td><td>102.0 ± 3.2</td><td>81.7 ± 2.5</td><td>79.3 ± 2.7</td><td>65.1 ± 2.5</td></tr><tr><td>5806921716937210496</td><td>BS Aps</td><td>0</td><td>0.58256590</td><td>0.76</td><td>0.24</td><td>51.1 ± 1.7</td><td>0</td><td>85.2 ± 3.2</td><td>67.1 ± 2.5</td><td>62.0 ± 2.3</td><td>56.4 ± 2.3</td></tr><tr><td>5947570591534602240</td><td>S Ara</td><td>0*</td><td>0.45184650</td><td>0.96</td><td>0.29</td><td>60.6 ± 2.0</td><td>0</td><td>...</td><td>82.1 ± 3.0</td><td>77.9 ± 2.8</td><td>66.7 ± 2.5</td></tr><tr><td>6045464303643755008</td><td>Cl* NGC 6121 SAW V2</td><td>0</td><td>0.53568150</td><td>0.92</td><td>...</td><td>54.1 ± 6.5</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>6045485228725626752</td><td>Cl* NGC 6121 SAW V32</td><td>0</td><td>0.57909353</td><td>0.75</td><td>0.25</td><td>61.0 ± 6.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>6045491928874778624</td><td>Cl* NGC 6121 SAW V33</td><td>0</td><td>0.61484317</td><td>0.85</td><td>0.31</td><td>55.8 ± 8.3</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>6045502305516138368</td><td>Cl* NGC 6121 SAW V15</td><td>0</td><td>0.44379722</td><td>1.20</td><td>0.30</td><td>43.7 ± 6.0</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>6473212809741024256</td><td>ASAS J200431-5352.3</td><td>1</td><td>0.30022700</td><td>0.28</td><td>...</td><td>14.2 ± 1.9</td><td>0</td><td>16.1 ± 1.9</td><td>13.2 ± 2.2</td><td>15.4 ± 2.2</td><td>14.1 ± 2.3</td></tr><tr><td>6483680332235888896</td><td>V Ind</td><td>0</td><td>0.47959915</td><td>1.12</td><td>0.29</td><td>56.2 ± 1.2</td><td>0</td><td>101.7 ± 2.6</td><td>77.9 ± 2.0</td><td>73.8 ± 2.0</td><td>61.9 ± 1.7</td></tr><tr><td>6526559499016401408</td><td>RV Phe</td><td>0</td><td>0.59641862</td><td>0.67</td><td>0.25</td><td>54.6 ± 4.8</td><td>1</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>6541769554459131648</td><td>BO Gru</td><td>1</td><td>0.28272000</td><td>0.18</td><td>...</td><td>13.7 ± 1.5</td><td>0</td><td>16.0 ± 2.3</td><td>12.2 ± 2.3</td><td>15.6 ± 2.6</td><td>12.7 ± 2.6</td></tr><tr><td>6570585628216929408</td><td>TY Gru</td><td>0</td><td>0.57006515</td><td>0.91</td><td>0.27</td><td>59.2 ± 2.1</td><td>0</td><td>112.8 ± 4.7</td><td>82.5 ± 3.3</td><td>76.2 ± 3.2</td><td>64.6 ± 3.2</td></tr><tr><td>6625215584995450624</td><td>AE PsA</td><td>0</td><td>0.54674000</td><td>1.13</td><td>0.32</td><td>65.1 ± 6.5</td><td>0</td><td>102.3 ± 8.7</td><td>76.6 ± 7.0</td><td>79.1 ± 7.6</td><td>64.5 ± 6.2</td></tr><tr><td>6662886605712648832</td><td>MT Tel</td><td>1</td><td>0.31689945</td><td>0.55</td><td>0.13</td><td>22.4 ± 1.4</td><td>0</td><td>35.6 ± 1.9</td><td>27.0 ± 1.6</td><td>26.0 ± 2.3</td><td>21.7 ± 1.9</td></tr><tr><td>6680420204104678272</td><td>V1645 Sgr</td><td>0</td><td>0.55292000</td><td>0.76</td><td>0.24</td><td>50.9 ± 1.9</td><td>0</td><td>74.0 ± 2.8</td><td>66.2 ± 2.5</td><td>65.8 ± 2.7</td><td>56.5 ± 2.4</td></tr><tr><td>6701821205809488384</td><td>ASAS J181215-5206.9</td><td>0</td><td>0.83753980</td><td>0.52</td><td>0.20</td><td>42.6 ± 5.6</td><td>0</td><td>...</td><td>72.0 ± 9.5</td><td>61.3 ± 8.2</td><td>54.2 ± 7.4</td></tr><tr><td>6771307454464848768</td><td>V0440 Sgr</td><td>0</td><td>0.47750000</td><td>0.77</td><td>0.30</td><td>62.7 ± 5.1</td><td>0</td><td>113.2 ± 8.9</td><td>84.0 ± 6.5</td><td>76.5 ± 6.4</td><td>67.1 ± 5.2</td></tr><tr><td>6787617919184986496</td><td>Z Mic</td><td>0</td><td>0.58692775</td><td>0.64</td><td>0.22</td><td>53.2 ± 1.8</td><td>0</td><td>92.1 ± 3.7</td><td>74.3 ± 3.0</td><td>69.3 ± 2.8</td><td>61.7 ± 2.6</td></tr><tr><td>6856027093125912064</td><td>ASAS J203145-2158.7</td><td>1</td><td>0.31071060</td><td>0.37</td><td>0.08</td><td>22.0 ± 2.0</td><td>0</td><td>29.7 ± 2.8</td><td>24.3 ± 2.6</td><td>21.1 ± 2.5</td><td>24.6 ± 4.2</td></tr><tr><td>6883653108749373568</td><td>RV Cap</td><td>0*</td><td>0.44774990</td><td>0.97</td><td>...</td><td>49.9 ± 2.0</td><td>0</td><td>90.8 ± 3.6</td><td>72.7 ± 3.1</td><td>62.4 ± 2.7</td><td>59.2 ± 2.7</td></tr><tr><td>6884361748289023488</td><td>YZ Cap</td><td>1</td><td>0.27345290</td><td>0.49</td><td>0.10</td><td>26.3 ± 2.1</td><td>0</td><td>35.8 ± 3.0</td><td>28.3 ± 2.8</td><td>29.1 ± 2.9</td><td>23.2 ± 2.6</td></tr></tbody></table>

**Notes.**

<footnote><sup>a</sup> Pulsation mode: 0 = RRab; 1 = RRc. Candidate Blazhko RRLs are marked with an asterisk.</footnote>

<footnote><sup>b</sup> RV amplitudes: 0 = current investigation; 1 = BW sample.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Panels (a) and (b): radial velocity curves vs. phase for the RRab V Ind based on metallic lines Fe I (λ: 4045.81, 4063.59, 4071.74 Å) and Sr II (λ: 4077.71 Å), and on Balmer (λ: 6562.80, 4861.36, 4340.46, 4101.74 Å) lines. Measurements based on different spectrographs are marked with different colors. The cyan lines show the analytical fits. The RV amplitude ($A(V_r)$ km s⁻¹) and the pulsation period (days) are also labeled. Panels (c) and (d): same as the top, but for the RRc MT Tel.</figcaption></figure>

V Ind. In spite of the difference in amplitude, the secondary features (bumps, dips) appear to be largely co-phased.

We analytically fit the RV curves by using the PLOESS algorithm (Braga et al. 2018). The cyan lines in Figure 1 show these curves. Note that RVs for V Ind based on $H_γ$ show a well-defined “spike” around phases 0.70–0.76, near its point of minimum radius, as implied by the optical and MIR light curves of this star in the top panels of Figure 2.

This secondary spike feature was predicted by pulsation models (see Figures 3 and 16 in Bono & Stellingwerf 1994 and Figure 12 in Marconi et al. 2015) suggesting that it is associated with the formation and propagation of shocks. Note that the duration of the spike is ~25 minutes and we succeeded in tracing it only because the du Pont spectra were collected with exposure times of 400 ± 200 s. Moreover, exposure times for V Ind are ≈30% shorter than the mean, since it has a
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Top: optical (V, left) and MIR (W1, right) light curves for the RRab V Ind. The optical data were provided by ASAS-SN (Shappee et al. 2014; Jayasinghe et al. 2018), while the MIR ones by NEOWISE (Mainzer et al. 2011). The thin red lines show the analytical fits. Bottom: same as the top, but for the RRc MT Tel.</figcaption></figure>

brighter mean visual magnitude (V = 9.92 ± 0.05 mag). A similar feature is also present among other Balmer (mainly Hβ, Hδ) and metallic (mainly Sr) lines, but it is less pronounced.

Note that the absolute minimum in RV for V Ind takes place, due to the phase lag (Castor 1971), at the same phases of the absolute maximum in the light curves. However, the absolute minimum in luminosity is not at the same phases of the absolute maximum in RV. Indeed, the latter takes place around the phases in which the optical light curves show a well-defined bump. The difference appears to be caused by the interplay between the shock propagation and the increased efficiency of the convective transport in approaching the phases of minimum radius.

The RV curve of MT Tel is smoother than that of V Ind over the entire pulsation cycle. For this variable we have extensive measurements from both du Pont and HARPS spectrographs. Inspection of panels (c) and (d) of Figure 1 reveals that the RV amplitude from the HARPS data set is ~2–3 km s⁻¹ larger than that from du Pont. MT Tel is the only variable for which we found such a difference. We performed a number of tests that verified the reality of this offset. It cannot be caused by phasing problems, since the combined optical and NIR photometry allowed us very accurate phasing of the two data sets. The analytical fit displayed in Figure 1 is mainly following the du Pont measurements, because they have a more uniform sampling over the pulsation cycle. We suggest that the difference is intrinsic and possibly caused by the presence of a secondary modulation (Blazhko effect), even though the light curves do not show yet a clear modulation at fixed phase. Note that HARPS spectra were collected in a single night seven years before the du Pont spectra (over five consecutive nights) and the periodogram lacks relevant secondary peaks. This finding would suggest that accurate RV measurements might have a stronger sensitivity in detecting secondary modulations and in constraining the tomography of variable star atmospheres. In the next paper we will explore more aspects of these interesting issues.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Bailey diagram, amplitude vs. logarithmic period. Panels (a) and (b) display the luminosity amplitude in the visual and in the W1 band. The black lines display observed relations for OoI (solid), Oo-intermediate (dotted) and OoII (dashed) RRLs provided by Fabrizio et al. (2019). Panels (c) and (d) show RV amplitudes based either on metallic or on Balmer lines. The red symbols plotted in panel (c) display BW RRLs. Hydrogen RV amplitudes are plotted with different colors (panel (d)). Blazhko RRLs are marked with crosses. The error bars plotted in the bottom left corner display uncertainties in amplitudes.</figcaption></figure>

## 4. Results

We have investigated on a more quantitative basis the variation of the luminosity amplitudes when moving from optical to MIR, and to RV amplitudes and when moving from weak to strong lines. In Figure 3 we show the Bailey diagram for the RRLs for which we estimated both luminosity and velocity amplitudes. Panels (a) and (b) display quite clearly that the luminosity amplitudes in the W1 band are roughly a factor of 3 smaller than in the V band. This trend is expected because the luminosity variation in the former band is mainly driven by temperature variations, while in the latter it is mainly driven by radius variations. Amplitude dispersion at fixed pulsation period also is smaller in the infrared. It is of the order of 20% for the V band ($\sigma_V \sim 0.19$ mag), while for the W1 band it decreases to 10% ($\sigma_{W1} \sim 0.03$ mag). This outcome applies to both RRab and RRc variables and suggests that the dispersion is driven by phenomena affecting the temperature variation
along the pulsation cycle. Moreover, Blazhko RRLs (black crosses) in our sample display a variation in the visual band of ~0.4 mag and of ≈0.1 mag in the W1 band (Jurcsik et al. 2018).

Note that the current sample is quite representative of the RRL pulsation properties; indeed, the analytical relations for OoI, Oosterhoff-intermediate, and OoII RRLs (Fabrizio et al. 2019) provide a plausible fit of optical amplitudes.

Figures 3(c) and (d) display, for the first time, a Bailey diagram based on RV amplitudes. In panel (c) together with the current velocity amplitudes we also included the velocity amplitudes for field and cluster RRLs (red symbols) for which a BW solution is available in the literature. Interestingly, the relative dispersion in RV amplitudes is smaller than for visual amplitudes. Indeed, the RV amplitudes are purely tracing the radius variation along the pulsation cycle. The result is quite clear for RRab for which we find that the dispersion is of the order of 10% ($\sigma_{V_r} \sim 6$ km s⁻¹). Moreover, Blazhko RRLs show a variation in velocity amplitude significantly smaller (~20%) than in visual amplitude (~40%). Unfortunately, we cannot reach firm conclusions concerning RRc, since the sample is too small.

Figure 3(d) also shows an interesting trend. The velocity amplitudes for RRab variables based on Balmer lines attain an almost constant value when moving from the hot (short-period) to the cool (long-period) edge of the fundamental IS. The mean RV amplitudes decrease from $H_\alpha$ to $H_\beta$, but their values are constant over the entire period range. The dispersion is also modest and ranges from 5 to 11 km s⁻¹.

The RV amplitudes plotted in Figures 3(c) and (d) clearly show that the dispersion, at fixed period, is significantly larger than intrinsic errors; thus, suggesting that other physical mechanisms can contribute to the observed dispersion. To investigate this intrinsic property quantitatively, we constructed four HB models covering a broad range in chemical compositions ($Z = 0.0001, 0.001, 0.008, 0.02$; $Y = 0.245, 0.245, 0.256, 0.27$). The stellar masses ($M = 0.75, 0.65, 0.57, 0.55$ M☉) were selected to be centrally located inside the IS. The reason for this broad range in metallicity was driven by the relevant change in the evolutionary path inside the IS when moving from metal-poor to metal-rich stellar structures. We also constructed, following Marconi et al. (2015), a set of nonlinear, convective hydrodynamical models of RRL stars along the HB evolutionary models to investigate the variation of the pulsation properties in their off-ZAHB evolution.

The top panels of Figure 4 show small Hertzsprung-Russell diagram regions of the HB evolutionary models (black lines) together with the variable models that are pulsationally stable. The color and the size of the symbols change along the off-ZAHB evolution.

As expected, the evolution across the IS changes significantly when moving from the metal-poor to the metal-rich regime. Indeed, for the most metal-poor chemical composition the off-ZAHB evolution is redward, i.e., HB models become systematically brighter and cooler as a function of time. The evolutionary path becomes more complex at metal-intermediate regime, because the evolutionary path undergoes, in sequence, a redward, a blueward and eventually a redward evolution approaching its asymptotic giant branch track. This means that we are dealing with triple points, i.e., evolutionary phases covering the same range in effective temperature, but slightly different luminosities ($\Delta \log L/L_☉ \sim 0.07$, $Z = 0.0001$;

~0.09, $Z = 0.001$; ~0.11, $Z = 0.008$; ~0.05, $Z = 0.02$). This means a steady decrease in surface gravity and, in turn, a steady increase in pulsation period.

The change in the physical structure also implies variation in the pulsation properties. The middle and the bottom panels of Figure 4 display the comparison between observed and predicted amplitudes. This is the first time that evolutionary effects are investigated in detail along the off-ZAHB evolution. A glance at the data plotted in this figure shows that the models, when moving from metal-poor to metal-rich structures, attain, as expected, shorter periods and larger amplitudes. Interestingly enough, the most metal-poor models display a steady decrease in amplitude when moving from hotter to cooler effective temperatures. Note that the smaller the symbols, the faster the evolutionary phase. This means that the lack of long-period small amplitude RRLs is a mix between the probability of finding them and the sample size.

The comparison between theory and observations becomes even more interesting in the metal-intermediate and in the metal-rich regimes, because pulsation models associated with blueward and redward evolutionary phases attain very similar values in the Bailey diagram. This means that RRLs in their off-ZAHB evolution move back and forth in the Bailey diagram covering a modest range in period, but a large range in amplitude. This indicates that a significant fraction of the dispersion in the Bailey diagram is caused by the evolutionary status of RRLs. Theory and observations agree quite well over the entire metallicity range. Note that predicted amplitudes are affected by uncertainties in the mixing length parameter adopted to treat time-dependent convective transport in hydrodynamical models of radial variables (Bono & Stellingwerf 1994). However, we are using predicted amplitudes in the relative sense, i.e., we are interested in the variation of the pulsation amplitudes when moving from the blue to the red edge of the IS and not in their absolute value.

The comparison in the Bailey diagram for RV amplitudes shows similar results. Note that to compare theory and observations we adopted a projection factor, the ratio between pulsation velocity and RV, $p = 1.33$ (Jones et al. 1992). The comparison in this plane is even more compelling because RV amplitudes are marginally affected by nonlinear phenomena and by secondary modulations. Interestingly, predicted and observed velocity amplitudes display a quadratic variation when moving from the blue to the red edge of the IS, while the luminosity amplitudes show a more linear trend.

## 5. Conclusions

We investigated the properties of the Bailey diagram based on luminosity (optical/MIR) and on RV amplitudes. We collected homogeneous RV measurements for four dozen of field RRLs. The current measurements double the sample of RRLs with a good coverage of the entire pulsation cycle. We investigated the dependence of luminosity amplitudes on nonlinear phenomena, secondary modulations, and evolutionary effects. The impact of these effects is mitigated when moving from the visual to the MIR bands, because the variations along the pulsation cycle are mainly driven by radius than by temperature changes. This trend becomes even more compelling in dealing with RV amplitudes; indeed, we found that the dispersion, at fixed period, is systematically smaller when compared with luminosity amplitudes. Moreover,
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Top: Hertzsprung-Russell diagram for HB models (Pietrinferni et al. 2006) constructed by assuming different chemical compositions and stellar masses (see labeled values). The symbols mark the position of the pulsation models we constructed by using evolutionary prescriptions (chemical composition, stellar mass, luminosity, effective temperature). The light blue numbers mark the first three and the last three models in the sequence of pulsation models. The symbols are large and bluish close to the ZAHB and become smaller and greenish-reddish along the off-ZAHB evolution. Middle: Bailey diagram for visual amplitudes. The symbols are the same as in the top panels and they are connected with solid lines following the same color-coding. Bottom: Bailey diagram for RV amplitudes based on metallic lines. The red symbols display BW RRLs.</figcaption></figure>

Blazhko RRLs display smaller variations in velocity amplitudes than in luminosity amplitudes.

We constructed four different HB evolutionary models covering a broad range in metal contents (Z = 0.0001, 0.001, 0.008, 0.02; Y = 0.245, 0.245, 0.256, 0.27) and a sizable sample of nonlinear, convective hydrodynamical models of RRLs to constrain the impact of the off-ZAHB evolution on their pulsation properties. We found that evolutionary effects

take account of the vertical structure (dispersion in amplitude) of the Bailey diagram based on RV amplitudes. Thus, we suggest that the dispersion in period is mainly caused by variation in the intrinsic parameters of RRLs (stellar mass, chemical composition).

The current findings are suggesting that RV amplitudes are not only independent of distance and reddening, but also less prone to amplitude variations caused by nonlinear phenomena
and secondary modulations. This means that RV amplitudes are
an optimal diagnostic for the variation of the mean effective
temperature across the IS of RRab.

We thank an anonymous referee for constructive suggestions on
an early draft of this Letter. We acknowledge financial support
from INAF/ASI 2014-049-R.0 (SSDC), US NSF under grants
AST-1714534 (M.M., J.P.M.) and AST-1616040 (C.S.), and DFG
for SFB 881 (Project-ID 138713538; B.L., E.K.G.). Some of the
observations were obtained with SALT. This publication uses
NEOWISE data (JPL/Caltech, University of Arizona, NASA).

ORCID iDs

G. Bono https://orcid.org/0000-0002-4896-8841

V. F. Braga https://orcid.org/0000-0001-7511-2830

M. Fabrizio https://orcid.org/0000-0001-5829-111X

C. Sneden https://orcid.org/0000-0002-3456-5929

M. Marconi https://orcid.org/0000-0002-1330-2927

G. Fiorentino https://orcid.org/0000-0003-0376-6928

B. Chaboyer https://orcid.org/0000-0003-3096-4161

A. Kunder https://orcid.org/0000-0002-2808-1370

M. Marengo https://orcid.org/0000-0001-9910-9230

C. E. Martínez-Vázquez https://orcid.org/0000-0002-9144-7726

M. Monelli https://orcid.org/0000-0001-5292-6380

J. Neeley https://orcid.org/0000-0002-8894-836X

P. B. Stetson https://orcid.org/0000-0001-6074-6830

E. Valenti https://orcid.org/0000-0002-6092-7145

A. R. Walker https://orcid.org/0000-0002-7123-8943

References

Bailey, S. I., Leland, E. F., Woods, I. E., & Pickering, E. C. 1919, AnHar, 78, 195

Blažko, S. 1907, AN, 175, 325

Bono, G., Caputo, F., Castellani, V., & Marconi, M. 1997, A&AS, 121, 327

Bono, G., Caputo, F., Castellani, V., et al. 2003, MNRAS, 344, 1097

Bono, G., Dall'Ora, M., Caputo, F., et al. 2011, in RR Lyrae Stars, Metal-Poor Stars, and the Galaxy, Vol. 5, ed. A. McWilliam (Pasadena, CA: The Observatories of the Carnegie Institution of Washington), 1

Bono, G., & Stellingwerf, R. F. 1994, ApJS, 93, 233

Braga, V. F., Stetson, P. B., Bono, G., et al. 2016, AJ, 152, 170

Braga, V. F., Stetson, P. B., Bono, G., et al. 2018, AJ, 155, 137

Castor, J. I. 1971, ApJ, 166, 109

Chadid, M., Sneden, C., & Preston, G. W. 2017, ApJ, 835, 187

da Costa, G. S., Rejkuba, M., Jerjen, H., & Grebel, E. K. 2010, ApJL, 708, L121

Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013a, ApJ, 763, 32

Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013b, ApJ, 765, 154

Evans, D. W., Riello, M., de Angeli, F., et al. 2018, A&A, 616, A4

Fabrizio, M., Bono, G., Braga, V. F., et al. 2019, ApJ, 882, 169

Fabrizio, M., Nonino, M., Bono, G., et al. 2011, PASP, 123, 384

Fiorentino, G., Monelli, M., Stetson, P. B., et al. 2017, A&A, 599, A125

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1

Jayasinghe, T., Kochanek, C. S., Stanek, K. Z., et al. 2018, MNRAS, 477, 3145

Jones, R. V., Carney, B. W., Storm, J., & Latham, D. W. 1992, ApJ, 386, 646

Jurcsik, J., Hajdu, G., Dékány, I., et al. 2018, MNRAS, 475, 4208

Kunder, A., Stetson, P. B., Cassisi, S., et al. 2013, AJ, 146, 119

Madore, B. F., Hoffman, D., Freedman, W. L., et al. 2013, ApJ, 776, 135

Mainzer, A., Bauer, J., Grav, T., et al. 2011, ApJ, 731, 53

Marconi, M., Coppola, G., Bono, G., et al. 2015, ApJ, 808, 50

Marrese, P. M., Marinoni, S., Fabrizio, M., & Altavilla, G. 2019, A&A, 621, A144

Monelli, M., Fiorentino, G., Bernard, E. J., et al. 2017, ApJ, 842, 60

Neeley, J. R., Marengo, M., Bono, G., et al. 2017, ApJ, 841, 84

Netzel, H., Smolec, R., Soszyński, I., & Udalski, A. 2018, MNRAS, 480, 1229

Oosterhoff, P. T. 1939, Obs, 62, 104

Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2006, ApJ, 642, 797

Pojmanski, G. 2002, AcA, 52, 397

Preston, G. W. 1959, ApJ, 130, 507

Preston, G. W., Sneden, C., Chadid, M., Thompson, I. B., & Shectman, S. A. 2019, AJ, 157, 153

Prudil, Z., & Skarka, M. 2017, MNRAS, 466, 2602

Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48

Storm, J., Carney, B. W., & Latham, D. W. 1994, A&A, 290, 443

Torrealba, G., Catelan, M., Drake, A. J., et al. 2015, MNRAS, 446, 2251