Issue 
A&A
Volume 607, November 2017



Article Number  A33  
Number of page(s)  16  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201731230  
Published online  03 November 2017 
Pulsating lowmass white dwarfs in the frame of new evolutionary sequences
V. Asteroseismology of ELMV white dwarf stars
^{1} Grupo de Evolución Estelar y Pulsaciones, Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque s/n, 1900 La Plata, Argentina
email: lcalcaferro@fcaglp.unlp.edu.ar; acorsico@fcaglp.unlp.edu.ar; althaus@fcaglp.unlp.edu.ar
^{2} Instituto de Astrofísica La Plata, CONICETUNLP, Paseo del Bosque s/n, 1900 La Plata, Argentina
Received: 23 May 2017
Accepted: 19 July 2017
Context. Many pulsating lowmass white dwarf stars have been detected in the past years in the field of our Galaxy. Some of them exhibit multiperiodic brightness variation, therefore it is possible to probe their interiors through asteroseismology.
Aims. We present a detailed asteroseismological study of all the known lowmass variable white dwarf stars based on a complete set of fully evolutionary models that are representative of lowmass Hecore white dwarf stars.
Methods. We employed adiabatic radial and nonradial pulsation periods for lowmass white dwarf models with stellar masses ranging from 0.1554 to 0.4352 M_{⊙} that were derived by simulating the nonconservative evolution of a binary system consisting of an initially 1 M_{⊙} zeroage mainsequence (ZAMS) star and a 1.4 M_{⊙} neutron star companion. We estimated the mean period spacing for the stars under study (where this was possible), and then we constrained the stellar mass by comparing the observed period spacing with the average of the computed period spacings for our grid of models. We also employed the individual observed periods of every known pulsating lowmass white dwarf star to search for a representative seismological model.
Results. We found that even though the stars under analysis exhibit few periods and the period fits show multiplicity of solutions, it is possible to find seismological models whose mass and effective temperature are in agreement with the values given by spectroscopy for most of the cases. Unfortunately, we were not able to constrain the stellar masses by employing the observed period spacing because, in general, only few periods are exhibited by these stars. In the two cases where we were able to extract the period spacing from the set of observed periods, this method led to stellar mass values that were substantially higher than expected for this type of stars.
Conclusions. The results presented in this work show the need for further photometric searches, on the one hand, and that some improvements of the theoretical models are required on the other hand in order to place the asteroseismological results on a firmer ground.
Key words: stars: evolution / stars: interiors / stars: oscillations / white dwarfs
© ESO, 2017
1. Introduction
White dwarf (WD) stars are the last stage in the life of the majority of stars (Winget & Kepler 2008; Fontaine & Brassard 2008; Althaus et al. 2010). Most WDs have envelopes rich in H, and they define the spectral class DA WD, whose distribution peaks at 0.59 M_{⊙}. This class also shows a peak at low mass: M_{⋆}/M_{⊙} ≲ 0.45. These stars are thought to be the result of strong massloss episodes in interactive binary systems before the He flash during the red giant branch phase of lowmass stars (Althaus et al. 2013; Istrate et al. 2016b, for recent studies). At variance with average WDs with C and O cores, they are expected to contain He cores, since He burning is avoided. Specifically, this interactive binary evolutionary scenario is thought to be the most plausible origin for the socalled extremely lowmass (ELM) WDs, which have masses below ~ 0.18−0.20 M_{⊙}.
In the past years, numerous lowmass WDs, including ELM WDs, have been discovered via the ELM survey and the SPY and WASP surveys (see Koester et al. 2009; Brown et al. 2010, 2012, 2013; Maxted et al. 2011; Kilic et al. 2011, 2012, 2015; Gianninas et al. 2014a, 2015). The detection of pulsation g modes (gravity modes) in some of them (Hermes et al. 2012, 2013b,a; Kilic et al. 2015; Bell et al. 2015, 2017) has given rise to a new class of variable WDs, the ELMVs. These pulsating lowmass WDs provide us an exceptional chance for probing the interiors of these stars and possibly to test their formation scenarios by employing the tools of asteroseismology. Since g modes in ELMVs are restricted mainly to the core regions (Steinfadt et al. 2010; Córsico et al. 2012b; Córsico & Althaus 2014a), we would be able to constrain their core chemical structure. Furthermore, as shown by stability computations (Córsico et al. 2012b; Van Grootel et al. 2013; Córsico & Althaus 2016), a combination of the κ−γ mechanism (Unno et al. 1989) and the “convective driving” mechanism (Brickhill 1991), both acting at the Hionization region, excite longperiod g modes, in agreement with observations. Moreover, some unstable shortperiod g modes could be driven by the ε mechanism as a result of stable H burning (Córsico & Althaus 2014b).
In addition to ELM stars, several objects are considered as their precursors, the socalled preELMs. These stars exhibit metals in their atmospheres (e.g., Gianninas et al. 2014b; Hermes et al. 2014; Istrate et al. 2016b). Interestingly enough, pulsations in a number of objects have been detected in the past years (Maxted et al. 2013, 2014; Zhang et al. 2016; Gianninas et al. 2016; Corti et al. 2016). Evolutionary models that consider only element diffusion cannot explain these properties (e.g., Córsico et al. 2016; Istrate et al. 2016a) and might be an indication that the missing physics might also affect the evolution of the objects on the cooling track (e.g., the thickness of the H envelope). Moreover, there are indications that a preELMV WD will later be observed as an ELMV (Fontaine et al. 2017).
The definition of an ELM WD is still under debate. In the context of the ELM survey, an ELM WD is defined as an object with surface gravity of 5 ≲ log (g) ≲ 7 and effective temperature in the range of 8000 ≲ T_{eff} ≲ 22 000 (e.g., Brown et al. 2010; Kilic et al. 2011; Brown et al. 2016). In addition, an ELM WD should be part of a tight binary system and therefore show shortperiod or highamplitude velocity variability (e.g., Brown et al. 2017). Córsico & Althaus (2014a) suggested defining an ELM WD as a WD resulting from a progenitor star that did not undergo Hshell flashes as the pulsational properties are quite different from those of systems that experienced flashes. However, this mass limit depends on the metallicity of the progenitor stars (Istrate et al. 2016b).
White dwarf asteroseismology has already proven to be a very useful technique for peering into the interior of several pulsating WDs, and it has been applied by employing two different methods: one method considers stellar models with parametrized chemical composition profiles, and the other involves fully evolutionary models characterized by a consistent chemical structure. The former has the advantage of allowing a full exploration of the parameter space (the total mass, the mass of the H envelope, and the chemical composition of the core, among others) to find an optimal asteroseismological model. Examples of this approach are the pioneer works by Bradley (1998, 2001). More recent works using this avenue are from Pech et al. (2006), Pech & Vauclair (2006), BischoffKim et al. (2008), Castanheira & Kepler (2008, 2009), Paparó et al. (2013), Bognár et al. (2016) and the recent developments of the core parameterization by Giammichele et al. (2016, 2017a,b). The second approach, developed at La Plata Observatory, is different but complementary, as it employs fully evolutionary models that are the result of the complete evolution of the progenitor stars, from the zeroage main sequence (ZAMS) until the WD phase. Examples of the application of this method to GW Virginis stars (pulsating PG1159 stars) are the works by Córsico et al. (2007a,b, 2008, 2009), Kepler et al. (2014) and Calcaferro et al. (2016). The method has also been applied in DBV WDs (Herich atmosphere) by Córsico et al. (2012a) and Bognár et al. (2014). In ZZ Ceti stars, this approach has been successfully employed by Kepler et al. (2012) and Romero et al. (2012, 2013). In particular, this method has the added advantage that the chemical structure of the background models is consistent with the preWD evolution.
Here, we apply our asteroseismological approach for the first time to the complete set of the known ELMVs, whose spectroscopic parameters we describe in Table 1, and whose periods are listed in Tables 2 to 10, along with the corresponding frequencies and amplitudes. SDSS J184037.78+642312.3 (hereafter J1840) is the ELMV that was discovered first (Hermes et al. 2012, 2013a). SDSS J111215.82+111745.0 (hereafter J1112) was reported by Hermes et al. (2013b). This case is particularly interesting because this star shows seven periods, two of which are very short, probably associated with p modes or radial modes. SDSS J151826.68+065813.2 (hereafter J1518) is the hottest ELMV so far, according to Hermes et al. (2013b). This star shows seven independent periods, which makes a more detailed asteroseismological analysis possible. SDSS J161431.28+191219.4 (hereafter J1614) is an ELMV, according to Hermes et al. (2013a). SDSS J222859.93+362359.6 (hereafter J2228) is the coolest EMLV known to date, according to Hermes et al. (2013a). This star exhibits only three independent periods in the range ~ [3255−6235] s, so these periods seem to be approximately in the asymptotic regime (see Córsico & Althaus 2014a). In particular, the period 6234.9 s is the longest period ever measured in a pulsating WD star. PSR J1738+0333 is a millisecond pulsar that has an ELMV companion (which we call J1738 for short), according to Kilic et al. (2015). This case is particularly interesting because it is the only binary system with a millisecond pulsar and a pulsating WD. SDSS J161831.69+385415.15 (hereafter J1618) is an ELMV, according to Bell et al. (2015). Finally, SDSS J1735+2134 (hereafter J1735) and SDSS J2139+2227 (hereafter J2139) are two recently detected ELMVs, according to Bell et al. (2017). In particular, the former has very long periods, which seem to be in the asymptotic regime of nonradial g modes (see Córsico & Althaus 2014a). It is worth mentioning how the spectroscopic masses of the ELM WDs are determined. Brown et al. (2017) showed that for the same metallicity of the progenitor stars there is a 15% difference in the mass of ELM WDs for the same log (g) and T_{eff} parameters, using either the Althaus et al. (2013) evolutionary tracks or the Istrate et al. (2016b) evolutionary tracks. We can adopt this difference as the true uncertainty in the spectroscopic determination of the masses of ELM WDs. In Fig. 1 we show the location of the different families of pulsating WDs, including all the known ELMV stars (red circles). The total number of ELMVs rises to nine because there is a high probability that the star discovered by Bell et al. (2017), SDSS J1355+1956, is a δ Scuti pulsator, as claimed by these authors. In this sense, it is important to stress here that some of the stars under analysis in this work may not be pulsating ELM WDs. The analysis reported by Brown et al. (2017) suggests that there are only four pulsating ELM WDs: J1840, J1112, J1518, and J1738. In addition, as discussed by Bell et al. (2015, 2017), the stars J1618, J1735, and J2139 may not be ELM WD stars. A measurement of the rate of period change for these stars could help to shed light on this question (Calcaferro et al. 2017). Despite this, given the exploratory nature of this work, we consider that these stars are genuine ELMVs and they are included in our analysis.
Stellar parameters (derived using 1D and 3D model atmospheres) and observed pulsation properties of all the known ELMV WD stars.
The five independent frequencies in the data of J1840 from Hermes et al. (2012).
The seven independent frequencies in the data of J1112 from Hermes et al. (2013b).
The seven independent frequencies in the data of J1518 from Hermes et al. (2013b).
In this paper we report a further step in the study of lowmass WD stars by performing an asteroseismological analysis of all the known ELMVs. This is the fifth work of a series dedicated to these stars. The first paper (Córsico & Althaus 2014a) was focused on the adiabatic properties of these stars; the second paper (Córsico & Althaus 2016) was dedicated to the nonadiabatic pulsation stability features of these stars. The third work (Córsico et al. 2016) was aimed at studying the pulsation properties of the preELMV WDs. The fourth paper (Calcaferro et al. 2017) was focused on studying the theoretical temporal rates of the period change of ELMV and preELMV stars. In this work, we follow the asteroseismological approach that employs fully evolutionary models resulting from the complete evolution of the progenitor stars. The employment of fully evolutionary models is a crucial requirement because some models (particularly those with the lowest mass) are characterized by strong Hnuclear burning that depends sensitively on the thickness of the H envelope, a quantity that results from the previous evolution. We employ the adiabatic radial (ℓ = 0) and nonradial (ℓ = 1,2) p and gmode pulsation periods computed in Córsico & Althaus (2014a) on lowmass Hecore WD evolutionary models with stellar masses ranging from 0.1554 to 0.4352 M_{⊙}, extracted from the computations of Althaus et al. (2013), which take the binary evolution of the progenitor stars into account.
The paper is organized as follows. A brief summary of the stellar models and the pulsational code employed is provided in Sect. 2. In Sect. 3 we describe the asteroseismological analyses we carried out. Next, in Sect. 3.1, we try to determine (when possible) the observed period spacing for the target stars, and then in Sect. 3.2 we constrain the stellar mass by comparing the observed period spacing with the average of the computed period spacings. In Sect. 3.3 we search for the bestfit asteroseismological model by comparing the individual periods from each ELMV star with theoretical periods from our grid of models. Finally, in Sect. 4 we summarize the main findings of this work.
The two independent frequencies in the data of J1614 from Hermes et al. (2013a).
The three independent frequencies in the data of J2228 from Hermes et al. (2013a).
The three independent frequencies in the data of J1738 from Kilic et al. (2015).
The three independent frequencies in the data of J1618 from Bell et al. (2015).
Fig. 1 Location of the known ELMVs (red circles) along with the other classes of pulsating WD stars (dots of different colors) in the log T_{eff}−log g plane. The nature of the three stars emphasized with squares surrounding the light green circles is currently unclear (they could be preELMV stars as well as SX Phe and/or δ Scuti stars). In parenthesis we include the number of known members of each class. Two postVLTP evolutionary tracks for Hdeficient WDs and two evolutionary tracks for lowmass Hecore WDs are plotted for reference. Dashed lines indicate the theoretical blue edge for the different classes of pulsating WDs. 
The four independent frequencies in the data of J1735 from Bell et al. (2017).
The three independent frequencies in the data of J2139 from Bell et al. (2017).
2. Evolutionary models and pulsational code
In this work, we have employed the fully evolutionary models of lowmass Hecore WDs generated with the LPCODE stellar evolution code. This code computes in detail the complete evolutionary stages that lead to the WD formation, allowing the study of the WD evolution consistently with the predictions of the evolutionary history of progenitors. Details of the LPCODE can be found in Althaus et al. (2005, 2009, 2013, 2015) and references therein. Here, we briefly mention the ingredients employed that are relevant for our analysis of lowmass Hecore WD (see Althaus et al. 2013, for details). The standard mixing length theory (MLT) for convection in the ML2 prescription is used (see Tassoul et al. 1990, for its definition), but adiabatic periods do not sensitively depend on the specific version of the MLT convection theory that is employed (Bradley 1998). We assumed the metallicity of the progenitor stars to be Z = 0.01. We considered the radiative opacities for arbitrary metallicity in the range of 0 to 0.1 from the OPAL project (Iglesias & Rogers 1996). Conductive opacities are those of Cassisi et al. (2007). For the mainsequence evolution, we considered the equation of state from OPAL for H and Herich compositions. We also adopted from Itoh et al. (1996) the neutrino emission rates for pair, photo, and bremsstrahlung processes, and for plasma processes, we included the treatment of Haft et al. (1994). For the WD regime we have employed an updated version of the equation of state of Magni & Mazzitelli (1979). The nuclear network takes into account 16 elements and 34 thermonuclear reaction rates for ppchains, CNO bicycle, He burning, and C ignition. We also considered timedependent diffusion due to gravitational settling and chemical and thermal diffusion of nuclear species following the multicomponent gas treatment of Burgers (1969). We have computed abundance changes according to element diffusion, nuclear reactions, and convective mixing, a treatment that represents a very significant aspect in evaluating the importance of residual nuclear burning during the cooling stage of lowmass WDs.
The pulsation analysis was carried out for radial (ℓ = 0) and nonradial (ℓ = 1,2) p and g modes, on the basis of the set of adiabatic and nonadiabatic pulsation periods presented in Córsico & Althaus (2014a, 2016), computed employing the adiabatic and nonadiabatic versions of the LPPUL pulsation code (Córsico & Althaus 2006; Córsico et al. 2006). The adiabatic version of the LPPUL pulsation code is coupled to the LPCODE evolutionary code and is based on a general NewtonRaphson technique that solves the fourthorder (secondorder) set of real equations and boundary conditions governing linear, adiabatic, nonradial (radial) stellar pulsations following the dimensionless formulation of Dziembowski (1971; see also, Unno et al. 1989). On the other hand, the nonradial (radial) nonadiabatic version of the LPPUL pulsation code solves the sixthorder (fourthorder) complex system of linearized equations and boundary conditions as given by Unno et al. (1989; see also Saio et al. 1983). Our nonadiabatic computations rely on the frozenconvection approximation, in which the perturbation of the convective flux is neglected. To compute the BruntVäisälä frequency (N), we follow the socalled Ledoux modified treatment (Tassoul et al. 1990; Brassard et al. 1991).
Realistic configurations of the evolutionary sequences for lowmass Hecore WD stars were derived by Althaus et al. (2013) by mimicking the binary evolution of progenitor stars. Binary evolution was assumed to be fully nonconservative, and the losses of angular momentum due to mass loss, gravitational wave radiation, and magnetic braking were considered. All of the Hecore WD initial models were derived from evolutionary calculations for binary systems consisting of an evolving mainsequence lowmass component (donor star) of initially 1 M_{⊙} and a 1.4 M_{⊙} neutron star companion as the other component. A total of 14 initial Hecore WD models with stellar masses of 0.1554, 0.1612, 0.1650, 0.1706, 0.1762, 0.1805, 0.1863, 0.1917, 0.2019, 0.2389, 0.2707, 0.3205, 0.3624, and 0.4352 M_{⊙} were computed for initial orbital periods at the beginning of the Roche lobe phase in the range of 0.9 to 300 d. The evolution of these models was computed down to the range of luminosities of cool WDs, including the stages of multiple thermonuclear CNO flashes at the beginning of the cooling branch.
3. Asteroseismological analysis
Asteroseismology has so far been applied to infer the fundamental parameters of numerous pulsating WD stars. Specifically, by comparing the observed frequencies (or periods) of pulsating WDs and appropriate theoretical models, it has been possible to infer details about their origin, internal structure, and evolution. The larger the number of frequencies detected in a given pulsating WD, the more information can be inferred, such as gravity, effective temperature, stellar mass, and also the internal chemical stratification, the rate of rotation, the existence of magnetic fields, the cooling timescale, and the core composition, to name a few. For instance, the works of Bradley (1998), Romero et al. (2012) and Giammichele et al. (2016, 2017a,b) have proven that asteroseismology is a powerful technique for exploring the interior of WDs.
In the next subsections, we describe the asteroseismological methods we employed.
3.1. Searching for a constant period spacing
In the asymptotic limit of highradial order k, nonradial g modes with the same harmonic degree ℓ are expected to be equally spaced in period (Tassoul 1980): (1)where N is the BruntVäisälä frequency. In principle, the asymptotic period spacing or the average of the computed period spacings calculated from a grid of models (with different masses and effective temperatures) can be compared with the mean period spacing exhibited by a pulsating WD star, and then a value of the stellar mass can be inferred. The initial step to do so is to obtain (if it exists) a mean period spacing underlying the observed periodicities. We searched for a constant period spacing in the data of the target stars by using the KolmogorovSmirnov (KS; see Kawaler 1988), the inverse variance (IV; see O’Donoghue 1994) and the Fourier transform (FT; see Handler et al. 1997) significance tests. In the KS test, the quantity Q is defined as the probability that the observed periods are randomly distributed. Thus, any uniform – or at least systematically nonrandom – period spacing present in the period spectrum of the star under analysis will appear as a minimum in Q. In the IV test, a maximum of the inverse variance will indicate a constant period spacing. Finally, in the FT test, we calculate the Fourier transform of a Dirac comb function (created from a set of observed periods), and then we plot the square of the amplitude of the resulting function in terms of the inverse of the frequency. Once again, a maximum in the square of the amplitude will indicate a constant period spacing.
Fig. 2 IV (upper panel), KS (middle panel), and FT significance (bottom panel) tests applied to the period spectrum of J1518 to search for a constant period spacing. The periods used here are those indicated in Table 4. 
Tables 2 to 10 show that the number of periods exhibited by all the known ELMV WDs varies from 2 to 7. In particular, because of the few periods exhibited by J1614, J2228, J1738, J1618 and J2139, it is not possible to search for a constant period spacing in these cases. For J1840, J1112, J1518 and J1735, however, we were able to carry out this procedure. Unfortunately, for J1840 and J1112, we were unable to estimate any unambiguous constant period spacing, and furthermore, there was no agreement between the significance tests. The reason might be that the periods exhibited by these stars are not fully in the asymptotic regime and/or there are not as many periods as required to determine a mean period spacing. However, for J1518 and J1735, as shown in Figs. 2 and 3, respectively, we found clear indication of a constant period spacing for the three independent significance tests for both stars. For J1518, it lies at roughly ΔΠ ~ 44 s, although there is also another possible value at ~ 22.2 s, both for the three significance tests. However, the latter is too short and probably represents the harmonic of the main period spacing (). In addition, a period separation of ~ 22 s is not likely to be the asymptotic period spacing because the resulting stellar mass would be prohibitively high (Córsico & Althaus 2014a). We therefore assume that the period spacing associated with J1518 is ΔΠ^{O} ~ 44 s. In the case of J1735, there is a possible value for the period spacing at ~ 26 s, but once again, as we mentioned before, this value is too low and we discarded it. Two other possibilities are at ~ 47 s and ~ 59 s, the latter being a more expectable value for the period spacing, according to the asymptotic predictions. We therefore adopt ΔΠ^{O} ~ 59 s as the period spacing for this star.
3.2. Determination of the stellar mass of J1518 and J1735 from the observed period spacing
In this section, we aim to estimate the masses of J1518 and J1735 by comparing the average of the computed period spacings () for our grid of models with the observed period spacing () determined in the previous section for each star. We caution that this approach is problematic: the period spacing in this type of stars could also be sensitive to the thickness of the outer H envelope in addition to the stellar mass (Tassoul et al. 1990; Fontaine & Brassard 2008). We defer a full exploration of this dependence to a future publication.
Fig. 4 Average of the computed dipole (ℓ = 1) period spacings () assessed in the range of the periods observed in J1518, corresponding to each WD model sequence considered in this work, in terms of the logarithm of the effective temperature. Numbers along each curve denote the stellar mass (in solar units). The observed period spacing derived for J1518 is depicted with a horizontal solid line. We also indicate the T_{eff} (vertical solid line) in the 3D model, together with its uncertainties (vertical dashed lines). 
The average of the computed period spacings is assessed as , where the “forward” period spacing is defined as ΔΠ_{k} = Π_{k + 1}−Π_{k} (k being the radial order), and n is the number of theoretical periods within the range of the periods observed in the target star. For J1518, Π_{k} ∈ [1330,3900] s, while for J1735, Π_{k} ∈ [3350,5000] s.
In Fig. 4 we show the run of the average of the computed period spacings (ℓ = 1) for J1518 in terms of the logarithm of the effective temperature for our ELM WD evolutionary sequences, along with the observed period spacing for J1518. The figure shows that the lower the values of , the greater the stellar mass. In this case, it clearly shows that such a low value (~ 44 s) of the observed period spacing would lead to a stellar mass greater than 0.4352 M_{⊙}, which is higher than expected for this type of stars. Hence, this analysis does not seem to lead to a proper value of the mass for J1518. The reason might be that this star is not pulsating in the asymptotic regime (see Córsico & Althaus 2014a).
In Fig. 5 we show the run of the average of the computed period spacings (ℓ = 1) for J1735 in terms of the logarithm of the effective temperature for our ELM WD evolutionary sequences, along with the observed period spacing for J1735. Once again, the value we obtain for the stellar mass is higher than expected (≳0.43 M_{⊙}), even though this star may be pulsating in the asymptotic regime. Alternatively, if the value we have obtained for the observed period spacing is real, in the sense that it can be associated with the behavior of high radial order g modes, then it would indicate that this star has a mass somewhat higher than 0.4352 M_{⊙}, and that their spectroscopic parameters (T_{eff},log g) could be incorrect.
In the next section we follow another approach to estimate the stellar mass and other features of all the known ELMVs through the search of theoretical models that best fit the individual observed periods. The advantage of this procedure is that once we have chosen a model, we have access to information of the star that is otherwise very difficult (if not impossible) to obtain by any other method.
3.3. Constraints from the individual observed periods: searching for the bestfit model
In this approach we search for a pulsation model that best matches the individual pulsation periods of a given star under study. The goodness of the match between the theoretical pulsation periods () and the observed individual periods () is measured by means of a merit function defined as (2)where m is the number of observed periods. The ELM model that shows the lowest value of χ^{2}, if exists, is adopted as the bestfit model. We assess the function χ^{2} = χ^{2}(M_{⋆},T_{eff}) for stellar masses of 0.1554, 0.1612, 0.1650, 0.1706, 0.1762, 0.1805, 0.1863, 0.1917, 0.2019, 0.2389, 0.2707, 0.3205, 0.3624, and 0.4352 M_{⊙}. For the effective temperature we also cover a wide range: 13 000 ≳ T_{eff} ≳ 6000 K.
Fig. 6 Inverse of the quality function of the period fit considering ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel) versus T_{eff} for J1840 (see text for details). The vertical strips depict the spectroscopic T_{eff} (solid vertical lines) and their uncertainties (dashed vertical lines) for the 1D (dark gray lines) and 3D models (red lines). 
We have carried out asteroseismological fits for all the known ELMV WD stars. This is the first time that this procedure is used for this type of stars. We started our analysis assuming that all of the observed periods correspond to g modes associated with ℓ = 1, and considering the set of observed periods, , of each star to compute the quality function given by Eq. (2). Next, we considered the case in which all of the observed periods correspond to g modes associated with ℓ = 2, and finally, we considered the case of a mix of g modes associated with ℓ = 1 and ℓ = 2. In the case of J1112, we performed a more detailed analysis. For this star we worked with two different sets of observed periods. On the one hand, the five longest periods, for which we carried out the analysis previously mentioned. On the other hand, we adopted the whole set of periods (seven) and explored the possibility that they correspond to a mix of g and p modes with ℓ = 1, and we also considered the case in which the observed periods correspond to radial (ℓ = 0) and p and g modes (ℓ = 1,2).
Figures 6 to 15 show the quantity (χ^{2})^{1} in terms of the effective temperature for different stellar masses for each known ELMV, taking the corresponding set of observed periods into account. We also include the effective temperatures and their uncertainties for the 1D (dark gray vertical lines) and 3D model atmosphere (red vertical lines) determinations. As mentioned before, the goodness of the match between the theoretical and the observed periods is measured by the value of χ^{2}: the better the period match, the lower the value of χ^{2} – in the figures, the greater the value of (χ^{2})^{1}. In some cases, there is a single maximum, and we adopt that model as the asteroseismological solution for that star. Sometimes, there are multiple possible solutions, and we need to employ some external constraint in order to choose one. Generally, the constraint is the uncertainty in the effective temperature, given by the spectroscopy. In some cases, when there are still multiple possible solutions, we choose a model with a mass as close as possible to the mass given by the spectroscopic determinations. It is important to mention that as it is more likely to observe ℓ = 1 than ℓ = 2 modes (because geometric cancellation effects become stronger with higher values of ℓ; see Dziembowski 1977), we typically chose, when possible, the asteroseismological solutions that fit observed periods with a larger number of ℓ = 1 modes.
Fig. 8 Same as Fig. 7, but for the case of the seven observed periods of J1112. In the upper panel, the case of radial (ℓ = 0) and p and g modes (ℓ = 1,2) is shown, while in the bottom panel the case of p and g modes with ℓ = 1 is displayed. 
Fig. 16 Location of the nine analyzed stars, according to the spectroscopic parameters from the 3D modelatmosphere (black circles), and the corresponding asteroseismological models adopted for each star (white circles), along with our evolutionary tracks of lowmass Hecore WDs in the log T_{eff}−log g plane. The grayshaded region bounded by the thin blue line corresponds to the instability domain of ℓ = 1g modes according to the nonadiabatic computations of Córsico & Althaus (2016). 
3.3.1. The case of J1840
In Fig. 6 we show the match between the theoretical and the five observed periods of J1840, assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). The case ℓ = 2 is only depicted for the sake of completeness, since we do not expect that a pulsating star can exhibit all the periods associated with ℓ = 2 and none of them correspond to ℓ = 1 because of geometric cancellation arguments (see above). For the ℓ = 1 case, the upper panel shows that there is more than one solution. In particular, the absolute maximum (the best solution) lies at a much higher effective temperature than that allowed by both spectroscopic determinations. The secondbest solution lies within the range of allowed T_{eff} for the 1D model atmosphere (T_{eff} = 9390 ± 140 K), so that we may adopt this model because it represents a good period fit. It corresponds to a mass of M_{⋆} = 0.2389 M_{⊙} at T_{eff} ~ 9529 K. For the ℓ = 1,2 case, once again there is no a single solution, and the best period fit has a very high value of T_{eff}. In the ranges of allowed T_{eff}, there is no unambiguous asteroseismological model. However, we may adopt the solution given for the model with M_{⋆} = 0.1805 M_{⊙} at T_{eff} ~ 9007 K, which is the best fit in the range of allowed T_{eff} for the 3D model atmosphere (T_{eff} = 9120 ± 140 K).
Comparison of the observed and theoretical periods for J1840, corresponding to the asteroseismological model with M_{⋆} = 0.2389 M_{⊙} and T_{eff} = 9529 K (ℓ = 1).
In order to determine how good the agreement of theoretical and observed periods is, we can compare them by computing the absolute period differences  δΠ  =  Π^{O}−Π^{T} . The results for J1840 are shown in Table 11 for the case of ℓ = 1. Column 6 of Table 11 shows the value of the linear nonadiabatic growth rate (η), defined as η (≡−ℑ(σ) / ℜ(σ), where ℜ(σ) and ℑ(σ) are the real and the imaginary part, respectively, of the complex eigenfrequency σ, computed with the nonadiabatic version of the LPPUL pulsation code (Córsico et al. 2006; Córsico & Althaus 2016). A value of η> 0 (η< 0) implies an unstable (stable) mode (see column 6 of Table 11). For ℓ = 1,2 the results are shown in Table 12.
Considering that the period fit of the ℓ = 1,2 case is better because it has a higher value of (χ^{2})^{1} than the solution of the ℓ = 1 case, and the mass of this model (M_{⋆} = 0.1805 M_{⊙}) is in line with the spectroscopic masses determined for this star ( and ), we adopt the model with M_{⋆} = 0.1805 M_{⊙} and T_{eff} = 9007 K as the asteroseismological solution for J1840, which has a value of T_{eff} in agreement with the spectroscopy (for the 3D model), even though this model has more ℓ = 2 than ℓ = 1 modes. We note that most of the periods of this model correspond to pulsationally unstable modes.
3.3.2. The case of J1112
Because of the two short periods at ~ 108 s and ~ 134 s in the pulsation spectrum of J1112, which are probably associated with p or radial modes, we divided the analysis for this star into two parts: on the one hand, we considered the five longest observed periods, assuming that they are all associated with ℓ = 1, ℓ = 2 or a mix of ℓ = 1 and ℓ = 2g modes. On the other hand, we considered the whole set of periods (seven), considering first, that they are associated with a mix of ℓ = 1p and g modes^{1}, and second, with a mix of p and g modes with ℓ = 1 and ℓ = 2, and also radial modes (ℓ = 0). At this point, it is worth mentioning that the reality of these two short periods are not definitively confirmed, as claimed by Hermes et al. (2013b). It is therefore worthwhile performing these separate analyses.
In Fig. 7 we show the match between the theoretical and the five longest observed periods of J1112, assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel, shown for completeness), and ℓ = 1,2 (bottom panel). In the case of ℓ = 1, the absolute maximum is at a very high effective temperature, and there are also many other solutions for low values of T_{eff}. Within the range of allowed T_{eff} for the 3D model atmosphere analysis (T_{eff} = 9240 ± 140 K), there may be a solution for the model with M_{⋆} = 0.3205 M_{⊙} at T_{eff} ~ 9253 K. We compare the observed and the theoretical periods as we did for the previous star, and the results are displayed in Table 13. If we consider the case of ℓ = 1,2, there are multiple local maxima that are either too hot or too cold in comparison with the allowed values of T_{eff}. Nevertheless, there may be a possible solution within the range of allowed T_{eff} (for the 3D model) for the case of M_{⋆} = 0.2389 M_{⊙} at T_{eff} ~ 9300 K. In Table 14 we show the comparison between the observed and the theoretical periods for this model.
Next, we considered the case in which the whole set of observed periods (seven) corresponds to p and g modes with ℓ = 1, and also the case in which it corresponds to radial (ℓ = 0) and p and g modes (ℓ = 1,2). The results are shown in Fig. 8. For the case of the mix of p and g modes with ℓ = 1 (bottom panel), the absolute maximum lies at a high value of T_{eff} and there is no unambiguous solution in the allowed ranges of T_{eff}. However, there may be one possible solution for M_{⋆} = 0.1612 M_{⊙} that lies in the allowed range of T_{eff} for the 1D model (T_{eff} = 9590 ± 140 K). In a more complete analysis, considering the mix of radial (ℓ = 0) and p and g modes (ℓ = 1,2), we find that the absolute maximum is at the same model and it represents a better match because it has a higher value of (χ^{2})^{1}, as shown in the top panel of Fig. 8. It is the best period fit for this case, and corresponds to M_{⋆} = 0.1612 M_{⊙} at T_{eff} ~ 9709 K. Hence, this represents a very good solution for the case of the whole set of periods. Once again, we show the comparison between the observed and the theoretical periods in Table 15. This table shows that one of the short periods may be associated with a p mode (with ℓ = 2) and the other period a radial mode (ℓ = 0).
Considering all the results from this analysis, we may conclude that the best solution corresponds to the model with M_{⋆} = 0.1612 M_{⊙} and T_{eff} = 9709 K, which is quite in line with the spectroscopic masses determined for this star ( and ) and also in line with the T_{eff} given by the spectroscopy (for the 1D model). Unfortunately, our nonadiabatic computations (see Table 15) predict that all the modes of this possible solution are pulsationally stable, which forces us to discard this solution. When we consider this and ignore the two shortest periods of this star, we could adopt the solution that is found with a mass of M_{⋆} = 0.2389 M_{⊙} and T_{eff} = 9300 K (see Table 14), although this value of the stellar mass does not agree as well with the masses resulting from the spectroscopy. We note, however, that all the periods of this model are associated with pulsationally unstable modes. Finally, the fact that we are not able to find an asteroseismological model with unstable modes with periods that fit the seven periods observed in J1112 (including the shortest ones) could be indicating that the periods at ~ 108 s and ~ 134 s reported by Hermes et al. (2013b) are not real. This calls for the need of further photometric work on this star.
3.3.3. The case of J1518
In Fig. 9 we show the match between the theoretical periods and the seven observed periods of J1518 assuming they are associated with g modes for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). In the ℓ = 1 case, we see multiple local maxima. However, in the range of allowed T_{eff} given by the 3D model atmosphere calculations (T_{eff} = 9650 ± 140 K), there is a possible solution that may be chosen as a representative model for J1518. This corresponds to M_{⋆} = 0.3205 M_{⊙} at T_{eff} ~ 9625 K. In Table 16 we compare the theoretical and the observed periods for this model. For the case of ℓ = 1,2, the best period fit lies at a very low T_{eff}, and there are multiple local maxima within the ranges of allowed T_{eff}. However, there is a possible solution characterized by M_{⋆} = 0.2707 M_{⊙} and T_{eff} ~ 9789 K, which, although not so clear in the figure, is the best period fit that lies within the ranges of allowed T_{eff} (K and K). We present in Table 17 the comparison between the observed and the theoretical periods for this case.
Considering that the solution for the ℓ = 1,2 case implies a much better fit (a higher value of (χ^{2})^{1}) than the one for the ℓ = 1 case, we may adopt the model with M_{⋆} = 0.2707 M_{⊙} and T_{eff} = 9789 K as the asteroseismological solution for this star, which is in line with the T_{eff} given by the spectroscopy. Moreover, most of the periods of the adopted model are associated with pulsationally unstable modes. It is necessary to stress, however, that none of the solutions found are in good agreement with the masses resulting from the spectroscopic determinations ( and ).
3.3.4. The case of J1614
In Fig. 10 we depict the match between the theoretical and the two observed periods of J1614, assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). It is worth mentioning in advance that this period fit is not reliable because this star only shows two independent periods.
In the case of ℓ = 1, there is no unambiguous solution, and the best solutions are located beyond the ranges of allowed T_{eff}. However, we may choose the model with M_{⋆} = 0.1762 M_{⊙} and T_{eff} ~ 8862 K that lies within the ranges of allowed T_{eff} (K and K), and also has a mass value consistent with the spectroscopic determinations. In Table 18 we show the comparison between the observed and the theoretical periods for this model. In the case of ℓ = 1,2, the best fit is located at a high value of T_{eff}, but the secondbest fit lies within the range of allowed T_{eff} (for the 3D model). It is characterized by M_{⋆} = 0.3205 M_{⊙}, at T_{eff} ~ 8610 K. However, as can be seen in Table 19, the comparison between the observed and the theoretical periods shows that both periods are associated with ℓ = 2, which is not usually the case because, as already stated, it is more likely to observe periods associated with ℓ = 1 than ℓ = 2. Hence, the asteroseismological model we adopt corresponds to the case of ℓ = 1, with M_{⋆} = 0.1762 M_{⊙} and T_{eff} = 8862 K, with a mass in line with the spectroscopic determinations ( and ) and a T_{eff} in agreement with the spectroscopy. Finally, as can be seen from Table 18, both periods are associated with pulsationally unstable modes.
3.3.5. The case of J2228
In Fig. 11 we show the match between the theoretical and the three observed periods of J2228, assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). In the case of ℓ = 1, there are multiple possible solutions and the bestfit solution is located at a high value of T_{eff}. Within the ranges of allowed T_{eff} (K and K), there is a possible solution for the model with M_{⋆} = 0.1650 M_{⊙} at T_{eff} ~ 7828 K. In Table 20 we show the comparison between the observed and the theoretical periods for this case. In the case of ℓ = 1,2, the absolute maximum lies very close to the ranges of allowed T_{eff}. It corresponds to the model with M_{⋆} = 0.1554 M_{⊙} at T_{eff} ~ 7710 K. In Table 21 we display the comparison between the observed and the theoretical periods for this model.
Since the solution for the case of ℓ = 1,2 (M_{⋆} = 0.1554 M_{⊙} and T_{eff} = 7710 K) implies a much better period fit than the solution for the case of ℓ = 1, since it lies at a value of T_{eff} almost compatible with the values given by spectroscopy, and because its mass is in line with the spectroscopic determinations for the mass ( and ), we adopt this model as the asteroseismological solution for J2228. According to our nonadiabatic computations (Table 21), most of the periods of the adopted model are associated with pulsationally unstable modes.
3.3.6. The case of J1738
In Fig. 12 we depict the match between the theoretical and the three observed periods of J1738, assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). In the case of ℓ = 1, the best solution lies at a very high value of T_{eff}, but the secondbest solution lies within the range of allowed T_{eff} (for the 1D atmosphere model determination, T_{eff} = 9130 ± 140 K). This solution is characterized by M_{⋆} = 0.4352 M_{⊙} at T_{eff} ~ 9177 K, and the comparison between the observed and the theoretical periods is shown in Table 22. In the case of ℓ = 1,2, the absolute maximum is located at a higher effective temperature than the allowed by spectroscopy (T_{eff} = 9130 ± 140 K and T_{eff} = 8910 ± 150 K, 1D and 3D models, respectively), and there are many other solutions. However, the models with M_{⋆} = 0.3205 M_{⊙} at T_{eff} ~ 8922 K and M_{⋆} = 0.4352 M_{⊙} at T_{eff} ~ 9192 K are relatively good period fits that lie within the ranges of allowed T_{eff}. When we analyze the periodtoperiod fit in detail, we see that the latter (shown in Table 23) may be more realistic because more modes are associated with ℓ = 1 than ℓ = 2. The opposite is true with the former, so we may rather choose the solution with M_{⋆} = 0.4352 M_{⊙} although is not the best one.
From this analysis, since the values of (χ^{2})^{1} for the two possible solutions are not significantly different and the value of the mass is the same for both of them, but the solution for the ℓ = 1 case has more periods associated with pulsationally unstable modes, we conclude that this is the best asteroseismological solution, characterized by M_{⋆} = 0.4352 M_{⊙} and T_{eff} = 9177 K, which is in line with the T_{eff} given by the spectroscopy (for the 1D model atmosphere computations). However, when we compare the mass of this model with the masses from the spectroscopic determinations ( and ), we see that they do not agree well. In summary, we cannot find any agreement between the asteroseismological and the spectroscopic results for J1738.
3.3.7. The case of J1618
In Fig. 13 we show the match between the theoretical and the three observed periods of J1618, assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). In the case of ℓ = 1, there is no unambiguous solution. Within the range of allowed T_{eff} for the 1D model (T_{eff} = 9144 ± 120 K), there is a possible solution for the model with M_{⋆} = 0.4352 M_{⊙} at T_{eff} ~ 9136 K, and in the range of allowed T_{eff} for the 3D model (T_{eff} = 8965 ± 120 K) there is another possible solution for the model characterized by M_{⋆} = 0.2019 M_{⊙} at T_{eff} ~ 8863 K. The latter may be more suitable as a solution because its mass is in line with the spectroscopic determinations of the stellar mass (and in comparison, the period fit for the other solution is not significantly better). In Table 24 we show the comparison between the observed and the theoretical periods for the model with M_{⋆} = 0.2019 M_{⊙}. The panel for the case of a mix of modes with ℓ = 1,2 shows an absolute maximum for a higher value of T_{eff} than allowed, and does not show any unambiguous solution in the ranges of allowed T_{eff}. However, there is a possible solution characterized by M_{⋆} = 0.1706 M_{⊙} at T_{eff} ~ 9076 K because although it is not the best period fit in the ranges of allowed T_{eff}, it corresponds to modes associated both with ℓ = 1 and ℓ = 2 (and not only ℓ = 2, see Table 25), and the mass is also quite well in line with the mass of the spectroscopic determination.
Although the solution for the ℓ = 1 case (M_{⋆} = 0.2019 M_{⊙}) has a mass slightly closer to the masses from the spectroscopic determinations for J1618 ( and ) than the solution for the ℓ = 1,2 case (M_{⋆} = 0.1706 M_{⊙}), we conclude that the model with M_{⋆} = 0.1706 M_{⊙} and T_{eff} = 9076 K, is a more suitable solution when we consider that the latter is a better match between the observed and the theoretical periods and also that this value of the stellar mass is more realistic for this type of stars. Moroever, the T_{eff} is in line with the spectroscopy. Hence, this is the model we adopt for J1618. We note that all of the periods of the adopted model are associated with pulsationally unstable modes.
3.3.8. The case of J1735
In Fig. 14 we plot the match between the theoretical and the four observed periods of J1735 assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). In the ℓ = 1 case, there are multiple local maxima that have values of T_{eff} that are either too high or too low in comparison with the range of allowed T_{eff} (T_{eff} = 7940 ± 130 K). However, there is a possible solution within that range, corresponding to the model with M_{⋆} = 0.3624 M_{⊙} at T_{eff} ~ 7991 K. In Table 26 we show the comparison between the observed and the theoretical periods for the mentioned model. We note, however, that in this stellar model, which constitutes a possible seismological solution for J1735, all the modes are pulsationally stable. As for the case of ℓ = 1,2, the bestfit models have values of T_{eff} higher than allowed. Although there is a possible solution within the range of allowed T_{eff} for the model characterized by M_{⋆} = 0.1650 M_{⊙} at T_{eff} ~ 7963 K, when we compare the observed and the theoretical periods, we find that they are all associated with ℓ = 2. The model with M_{⋆} = 0.1612 M_{⊙} that lies at a slightly higher value of T_{eff} than allowed (~ 8075 K) may therefore be a good solution for this case (see Table 27).
Taking into consideration the spectroscopic determination for the mass of J1735, , and comparing the quality of the period fit of the asteroseismological results, we find that the model with M_{⋆} = 0.1612 M_{⊙} and T_{eff} = 8075 K is an appropriate solution (with the T_{eff} almost compatible with the spectroscopy), and this is the one we adopt. This model has most of its periods associated with pulsationally unstable modes.
3.3.9. The case of J2139
In Fig. 15 we depict the match between the theoretical and the three observed periods of J2139 assuming they are associated with g modes, for the cases of ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel). In the first case, the absolute maximum, located at T_{eff} ~ 8221 K for a model with M_{⋆} = 0.1863 M_{⊙}, is very close to the range of allowed T_{eff} (T_{eff} = 7990 ± 130 K). Then, it represents a good solution. In the case of a mix of modes with ℓ = 1,2, the absolute maximum corresponds to the same model, and the three periods are associated with ℓ = 1. There are other possible solutions, but they lie far from the range of allowed T_{eff}, and inside this range, the fits are poor. For the model with M_{⋆} = 0.1863 M_{⊙}, we show the comparison between the observed and the theoretical periods in Table 28.
Considering these results, we may adopt the mentioned model, with M_{⋆} = 0.1863 M_{⊙} and T_{eff} = 8221 K, which although it is not quite in line with the spectroscopic result for the mass () and the T_{eff}, is a very good period fit (with all the periods associated with ℓ = 1g modes), and also has all the periods associated with pulsationally unstable modes.
Stellar masses (in solar units) for all of the studied ELMV WD stars.
Main characteristics of the adopted asteroseismological model for every known ELMV WD.
4. Summary and conclusions
We have presented a detailed asteroseismological study of all the known pulsating ELM WD stars (ELMVs), considering the pulsation spectrum they exhibit and employing the set of evolutionary models of Althaus et al. (2013). This is the fifth paper in a series of works dedicated to pulsating lowmass Hecore WDs (including ELMV WDs). The present paper was devoted to performing the first asteroseismological analysis of all the known ELMV stars. For this purpose we employed some asteroseismological tools. One of them is based on the comparison between the observed period spacing of the star under analysis with the average of the period spacings computed on our grid of models. We therefore first tried to determine the observed period spacing for each target star through three independent significance tests. Because the stars under study exhibit few periods, we could only follow this approach for the cases of the stars showing four periods or more, that is, for J1840, J1112, J1518, and J1735. However, for the first two stars we could not find any unambiguous constant period spacing. In the case of J1518 and J1735, on the other hand, we found a clear indication of a constant period spacing at roughly 44 s and 59 s, respectively, from the three significance tests we applied. After comparing these values with the average of the computed period spacings for our grid of models, we found that the resulting stellar masses (greater than 0.4352 M_{⊙} in both cases) are higher than expected for this type of stars. In the case of J1518, it may be associated with the fact that this star is not pulsating in the asymptotic regime (Córsico & Althaus 2014a). The case of J1735 is more intriguing because this star seems to be in that regime.
Next, we searched for the bestfit model, that is to say, the theoretical model that provides the best match between the individual pulsation periods exhibited by the star and the theoretical pulsation periods. We assessed the function χ^{2} = χ^{2}(M_{⋆},T_{eff}) (given by Eq. (2) of Sect. 3.3) for our complete set of model sequences, covering a wide range in effective temperatures (13 000 ≳ T_{eff} ≳ 6000 K). Because of the multiplicity of solutions, we were forced to employ some external constraints (for instance, the uncertainty in the T_{eff}, given by spectroscopy). We assumed that all of the observed periods correspond to ℓ = 1g modes and considered them to compute the quality function for each target star. We also considered the (unlikely) case in which all of the observed periods correspond to ℓ = 2g modes. Finally, we considered the case of a mix of ℓ = 1 and ℓ = 2g modes. For the particular case of the star J1112, we performed two different analyses. Since the two shortest periods reported for this star are not confirmed (Hermes et al. 2013b), we first carried out a period fit applied to the subset of the five longest periods exhibited by this star considering they are associated with ℓ = 1, ℓ = 2 and a mix of ℓ = 1 and ℓ = 2 g modes. Second, for the whole set of seven periods, we explored two possibilities: that all of the observed periods correspond to a mix of g and p modes (ℓ = 1), and also the case in which the observed periods correspond to radial (ℓ = 0) and p and g modes (ℓ = 1,2). In Table 29 we show a compilation of the mass determinations for the ELMVs both from spectroscopic (other works) and periodfit results (this work). Considering the obtained results, we found that the seismological mass is in good agreement with the spectroscopic determinations for J1840 (in the case of a mix of ℓ = 1,2g modes), J1614 (for the case of ℓ = 1g modes), J2228 (for the case of ℓ = 1,2 g modes), J1618 (for the case of ℓ = 1,2g modes), and J1735 (for the case of ℓ = 1,2g modes). We consider that there is a good agreement between the seismological and spectroscopic mass when the difference is below the uncertainty of 15%, that is, the typical difference in the mass value derived from independent sets of evolutionary tracks. Then, we conclude that for most of the target stars, the adopted models from the asteroseismological analyses have masses that are in line with the spectroscopic results. At variance with this, for four stars (J1738, J1518, J1112 and J2139) we obtained a higher value of the stellar mass in comparison with the spectroscopic determinations. In particular, gathering together the mass determinations for J1518, we conclude that there is no agreement between the mass given by the spectroscopy ( and ), the mass obtained from the comparison between the observed period spacing, and the average of the computed period spacings (which is higher than 0.4352 M_{⊙}), and the mass from the adopted asteroseismological model (M_{⋆} = 0.2707 M_{⊙}). We also mention that although we were able to adopt a seismological model for J1735 whose mass (M_{⋆} = 0.1612 M_{⊙}) is in line with the spectroscopic determination (M_{⋆} = 0.142 ± 0.010 M_{⊙}), we failed to find such an agreement for the mass resulting from the comparison between the observed period spacing and the average of the computed period spacings (M_{⋆} ≳ 0.44 M_{⊙}). Reversing the argument, as this star seems to be in the asymptotic regime, if the value we have obtained for the period spacing of this star were in fact associated with high radial order g modes, it could indicate that the stellar mass is higher than the determined by the spectroscopic and the periodtoperiod fit analysis, although in that case the star could not be classified as an ELMV WD star. Finally, it is worth mentioning that in general, the pulsation periods corresponding to the asteroseismological models adopted in this work for the analyzed ELMV WDs are pulsationally unstable, according to our nonadiabatic computations. This agreement between the adiabatic and nonadiabatic predictions gives more relevance to our asteroseismological results.
From the results presented in this paper for all the known ELMVs, the power of this approach is evident once again because in most of the cases we were able to constrain the value of the stellar mass. Moreover, when a model has been adopted, we can access additional information, as can be seen in Table 30, which is another advantage of asteroseismology. Taking into account these results, four of the stars we analyzed (J1840, J1518, J1738, and J2139) are not strictly ELM WD according to our definition previously stated, that is, the progenitors of these stars might have experienced multiple flashes.
In Fig. 16 we show the location of the nine analyzed ELMVs (according to the 3D modelatmosphere parameters) and the corresponding values of T_{eff} and log g of the asteroseismological models adopted for each star, along with our evolutionary tracks of lowmass Hecore WDs and the instability domain of ℓ = 1g modes computed by Córsico & Althaus (2016). As we mentioned, for five stars we found good agreement between the seismological mass and the spectroscopic one, and for the remaining four stars, the agreement is poor. In addition, Fig. 16 demonstrates that for eight out of nine stars we analyzed, the asteroseismological models are more massive (i.e., they are characterized by higher gravities) in comparison with the spectroscopy results. This systematic trend is also found in the case (not shown) in which the T_{eff} and log (g) values derived from calculations of 1D atmospheres are adopted. This trend could be related, in part, to the fact that we did not consider lowmass Hecore WD models characterized by outer H envelopes thinner than those predicted by the complete binary evolutionary history of the progenitor stars. Alternatively, it could be an indication that the spectroscopic determinations of log g and T_{eff} in this class of stars were not correct.
In this paper, we have considered lowmass Hecore WDs coming from solar metallicity progenitors, typical of the population of the Galactic disk. The threshold in the stellar mass value below which CNO flashes on the early WD cooling branch are not expected to occur is ~ 0.18 M_{⊙}. If we had adopted progenitors with lower metallicities, representative of the population of the Galactic halo, the threshold mass limit should be higher (see Serenelli et al. 2002; Nelson et al. 2004; Istrate et al. 2016b). The H envelope of the lowmass WDs should also be thicker than those obtained in Althaus et al. (2013; e.g., Istrate et al. 2016b). This means that when we assume that some of the ELMVs studied in this work are objects of the Galactic halo, the asteroseismological analysis should be based on evolutionary models coming from lowmetallicity progenitors, and therefore the characteristics of the asteroseismological models for each analyzed star could be different to those obtained in this work.
We here considered lowmass Hecore WD models characterized by thick outer H envelopes, consistent with the previous evolution. We are well aware that there are strong uncertainties about the precise value of the thickness of this envelope. We cannot discard that WD models with H envelopes thinner than those characterizing our set of models could result from binary evolution computations that assume different angularmomentum loss prescriptions due to mass loss, different initial massratio, etc., than that adopted in Althaus et al. (2013; see the detailed works by Istrate 2015; Istrate et al. 2016b). Asteroseismological analyses considering lowmass Hecore WD models characterized by thinner outer H envelopes will be the core feature of a future work.
Acknowledgments
We wish to thank our anonymous referee for the constructive comments and suggestions that greatly improved the original version of the paper. Part of this work was supported by AGENCIA through the Programa de Modernización Tecnológica BID 1728/OCAR, and by the PIP 11220080100940 grant from CONICET. This research made use of NASA Astrophysics Data System.
References
 Althaus, L. G., Serenelli, A. M., Panei, J. A., et al. 2005, A&A, 435, 631 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Althaus, L. G., Panei, J. A., Romero, A. D., et al. 2009, A&A, 502, 207 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Althaus, L. G., Córsico, A. H., Isern, J., & GarcíaBerro, E. 2010, A&ARv, 18, 471 [NASA ADS] [CrossRef] [Google Scholar]
 Althaus, L. G., Miller Bertolami, M. M., & Córsico, A. H. 2013, A&A, 557, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Althaus, L. G., Camisassa, M. E., Miller Bertolami, M. M., Córsico, A. H., & GarcíaBerro, E. 2015, A&A, 576, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bell, K. J., Kepler, S. O., Montgomery, M. H., et al. 2015, in 19th European Workshop on White Dwarfs, eds. P. Dufour, P. Bergeron, & G. Fontaine, ASP Conf. Ser., 493, 217 [Google Scholar]
 Bell, K. J., Gianninas, A., Hermes, J. J., et al. 2017, ApJ, 835, 180 [NASA ADS] [CrossRef] [Google Scholar]
 BischoffKim, A., Montgomery, M. H., & Winget, D. E. 2008, ApJ, 675, 1505 [NASA ADS] [CrossRef] [Google Scholar]
 Bognár, Z., Paparó, M., Córsico, A. H., Kepler, S. O., & Győrffy, Á. 2014, A&A, 570, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bognár, Z., Paparó, M., Molnár, L., et al. 2016, MNRAS, 461, 4059 [NASA ADS] [CrossRef] [Google Scholar]
 Bradley, P. A. 1998, ApJS, 116, 307 [Google Scholar]
 Bradley, P. A. 2001, ApJ, 552, 326 [NASA ADS] [CrossRef] [Google Scholar]
 Brassard, P., Fontaine, G., Wesemael, F., Kawaler, S. D., & Tassoul, M. 1991, ApJ, 367, 601 [NASA ADS] [CrossRef] [Google Scholar]
 Brickhill, A. J. 1991, MNRAS, 251, 673 [NASA ADS] [Google Scholar]
 Brown, W. R., Kilic, M., Allende Prieto, C., & Kenyon, S. J. 2010, ApJ, 723, 1072 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, W. R., Kilic, M., Allende Prieto, C., & Kenyon, S. J. 2012, ApJ, 744, 142 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, W. R., Kilic, M., Allende Prieto, C., Gianninas, A., & Kenyon, S. J. 2013, ApJ, 769, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, W. R., Gianninas, A., Kilic, M., Kenyon, S. J., & Allende Prieto, C. 2016, ApJ, 818, 155 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, W. R., Kilic, M., & Gianninas, A. 2017, ApJ, 839, 23 [NASA ADS] [CrossRef] [Google Scholar]
 Burgers, J. M. 1969, Flow Equations for Composite Gases (New York: Academic Press) [Google Scholar]
 Calcaferro, L. M., Córsico, A. H., & Althaus, L. G. 2016, A&A, 589, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Calcaferro, L. M., Córsico, A. H., & Althaus, L. G. 2017, A&A, 600, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cassisi, S., Potekhin, A. Y., Pietrinferni, A., Catelan, M., & Salaris, M. 2007, ApJ, 661, 1094 [NASA ADS] [CrossRef] [Google Scholar]
 Castanheira, B. G., & Kepler, S. O. 2008, MNRAS, 385, 430 [Google Scholar]
 Castanheira, B. G., & Kepler, S. O. 2009, MNRAS, 396, 1709 [NASA ADS] [CrossRef] [Google Scholar]
 Córsico, A. H., & Althaus, L. G. 2006, A&A, 454, 863 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., & Althaus, L. G. 2014a, A&A, 569, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., & Althaus, L. G. 2014b, ApJ, 793, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Córsico, A. H., & Althaus, L. G. 2016, A&A, 585, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Althaus, L. G., & Miller Bertolami, M. M. 2006, A&A, 458, 259 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Althaus, L. G., Miller Bertolami, M. M., & Werner, K. 2007a, A&A, 461, 1095 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Miller Bertolami, M. M., Althaus, L. G., Vauclair, G., & Werner, K. 2007b, A&A, 475, 619 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Althaus, L. G., Kepler, S. O., Costa, J. E. S., & Miller Bertolami, M. M. 2008, A&A, 478, 869 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Althaus, L. G., Miller Bertolami, M. M., & GarcíaBerro, E. 2009, A&A, 499, 257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Althaus, L. G., Miller Bertolami, M. M., & BischoffKim, A. 2012a, A&A, 541, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Romero, A. D., Althaus, L. G., & Hermes, J. J. 2012b, A&A, 547, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Córsico, A. H., Althaus, L. G., Serenelli, A. M., et al. 2016, A&A, 588, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Corti, M. A., Kanaan, A., Córsico, A. H., et al. 2016, A&A, 587, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dziembowski, W. A. 1971, Acta Astron., 21, 289 [NASA ADS] [Google Scholar]
 Dziembowski, W. 1977, Acta Astron., 27, 203 [NASA ADS] [Google Scholar]
 Fontaine, G., & Brassard, P. 2008, PASP, 120, 1043 [NASA ADS] [CrossRef] [Google Scholar]
 Fontaine, G., Istrate, A., Gianninas, A., Brassard, P., & Van Grootel, V. 2017, in 20th European White Dwarf Workshop, eds. P.E. Tremblay, B. Gaensicke, & T. Marsh, ASP Conf. Ser., 509, 347 [Google Scholar]
 Giammichele, N., Fontaine, G., Brassard, P., & Charpinet, S. 2016, ApJS, 223, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Giammichele, N., Charpinet, S., Brassard, P., & Fontaine, G. 2017a, A&A, 598, A109 [Google Scholar]
 Giammichele, N., Charpinet, S., Fontaine, G., & Brassard, P. 2017b, ApJ, 834, 136 [NASA ADS] [CrossRef] [Google Scholar]
 Gianninas, A., Dufour, P., Kilic, M., et al. 2014a, ApJ, 794, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Gianninas, A., Hermes, J. J., Brown, W. R., et al. 2014b, ApJ, 781, 104 [NASA ADS] [CrossRef] [Google Scholar]
 Gianninas, A., Kilic, M., Brown, W. R., Canton, P., & Kenyon, S. J. 2015, ApJ, 812, 167 [NASA ADS] [CrossRef] [Google Scholar]
 Gianninas, A., Curd, B., Fontaine, G., Brown, W. R., & Kilic, M. 2016, ApJ, 822, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Haft, M., Raffelt, G., & Weiss, A. 1994, ApJ, 425, 222 [NASA ADS] [CrossRef] [Google Scholar]
 Handler, G., Pikall, H., O’Donoghue, D., et al. 1997, MNRAS, 286, 303 [NASA ADS] [CrossRef] [Google Scholar]
 Hermes, J. J., Montgomery, M. H., Winget, D. E., et al. 2012, ApJ, 750, L28 [Google Scholar]
 Hermes, J. J., Montgomery, M. H., Gianninas, A., et al. 2013a, MNRAS, 436, 3573 [NASA ADS] [CrossRef] [Google Scholar]
 Hermes, J. J., Montgomery, M. H., Winget, D. E., et al. 2013b, ApJ, 765, 102 [Google Scholar]
 Hermes, J. J., Gänsicke, B. T., Koester, D., et al. 2014, MNRAS, 444, 1674 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Istrate, A. G. 2015, in 19th European Workshop on White Dwarfs, eds. P. Dufour, P. Bergeron, & G. Fontaine, ASP Conf. Ser., 493, 487 [Google Scholar]
 Istrate, A. G., Fontaine, G., Gianninas, A., et al. 2016a, A&A, 595, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Istrate, A. G., Marchant, P., Tauris, T. M., et al. 2016b, A&A, 595, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Itoh, N., Hayashi, H., Nishikawa, A., & Kohyama, Y. 1996, ApJS, 102, 411 [NASA ADS] [CrossRef] [Google Scholar]
 Kawaler, S. D. 1988, in Advances in Helio and Asteroseismology, eds. J. ChristensenDalsgaard, & S. Frandsen, IAU Symp., 123, 329 [Google Scholar]
 Kepler, S. O., Pelisoli, I., Peçanha, V., et al. 2012, ApJ, 757, 177 [NASA ADS] [CrossRef] [Google Scholar]
 Kepler, S. O., Fraga, L., Winget, D. E., et al. 2014, MNRAS, 442, 2278 [NASA ADS] [CrossRef] [Google Scholar]
 Kilic, M., Brown, W. R., Allende Prieto, C., et al. 2011, ApJ, 727, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Kilic, M., Brown, W. R., Allende Prieto, C., et al. 2012, ApJ, 751, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Kilic, M., Hermes, J. J., Gianninas, A., & Brown, W. R. 2015, MNRAS, 446, L26 [NASA ADS] [CrossRef] [Google Scholar]
 Koester, D., Voss, B., Napiwotzki, R., et al. 2009, A&A, 505, 441 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Magni, G., & Mazzitelli, I. 1979, A&A, 72, 134 [NASA ADS] [Google Scholar]
 Maxted, P. F. L., Anderson, D. R., Burleigh, M. R., et al. 2011, MNRAS, 418, 1156 [NASA ADS] [CrossRef] [Google Scholar]
 Maxted, P. F. L., Serenelli, A. M., Miglio, A., et al. 2013, Nature, 498, 463 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Maxted, P. F. L., Serenelli, A. M., Marsh, T. R., et al. 2014, MNRAS, 444, 208 [NASA ADS] [CrossRef] [Google Scholar]
 Nelson, L. A., Dubeau, E., & MacCannell, K. A. 2004, ApJ, 616, 1124 [NASA ADS] [CrossRef] [Google Scholar]
 O’Donoghue, D. 1994, MNRAS, 270, 222 [NASA ADS] [CrossRef] [Google Scholar]
 Paparó, M., Bognár, Z., Plachy, E., Molnár, L., & Bradley, P. A. 2013, MNRAS, 432, 598 [NASA ADS] [CrossRef] [Google Scholar]
 Pech, D., & Vauclair, G. 2006, A&A, 453, 219 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pech, D., Vauclair, G., & Dolez, N. 2006, A&A, 446, 223 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Romero, A. D., Córsico, A. H., Althaus, L. G., et al. 2012, MNRAS, 420, 1462 [Google Scholar]
 Romero, A. D., Kepler, S. O., Córsico, A. H., Althaus, L. G., & Fraga, L. 2013, ApJ, 779, 58 [NASA ADS] [CrossRef] [Google Scholar]
 Saio, H., Winget, D. E., & Robinson, E. L. 1983, ApJ, 265, 982 [NASA ADS] [CrossRef] [Google Scholar]
 Serenelli, A. M., Althaus, L. G., Rohrmann, R. D., & Benvenuto, O. G. 2002, MNRAS, 337, 1091 [NASA ADS] [CrossRef] [Google Scholar]
 Steinfadt, J. D. R., Bildsten, L., & Arras, P. 2010, ApJ, 718, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, M., Fontaine, G., & Winget, D. E. 1990, ApJS, 72, 335 [Google Scholar]
 Tremblay, P.E., Gianninas, A., Kilic, M., et al. 2015, ApJ, 809, 148 [NASA ADS] [CrossRef] [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial oscillations of stars, 2nd edn. (Tokyo: University of Tokyo Press) [Google Scholar]
 Van Grootel, V., Fontaine, G., Brassard, P., & Dupret, M.A. 2013, ApJ, 762, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Winget, D. E., & Kepler, S. O. 2008, ARA&A, 46, 157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhang, X. B., Fu, J. N., Li, Y., Ren, A. B., & Luo, C. Q. 2016, ApJ, 821, L32 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Stellar parameters (derived using 1D and 3D model atmospheres) and observed pulsation properties of all the known ELMV WD stars.
The seven independent frequencies in the data of J1112 from Hermes et al. (2013b).
The seven independent frequencies in the data of J1518 from Hermes et al. (2013b).
The three independent frequencies in the data of J2228 from Hermes et al. (2013a).
Comparison of the observed and theoretical periods for J1840, corresponding to the asteroseismological model with M_{⋆} = 0.2389 M_{⊙} and T_{eff} = 9529 K (ℓ = 1).
Main characteristics of the adopted asteroseismological model for every known ELMV WD.
All Figures
Fig. 1 Location of the known ELMVs (red circles) along with the other classes of pulsating WD stars (dots of different colors) in the log T_{eff}−log g plane. The nature of the three stars emphasized with squares surrounding the light green circles is currently unclear (they could be preELMV stars as well as SX Phe and/or δ Scuti stars). In parenthesis we include the number of known members of each class. Two postVLTP evolutionary tracks for Hdeficient WDs and two evolutionary tracks for lowmass Hecore WDs are plotted for reference. Dashed lines indicate the theoretical blue edge for the different classes of pulsating WDs. 

In the text 
Fig. 2 IV (upper panel), KS (middle panel), and FT significance (bottom panel) tests applied to the period spectrum of J1518 to search for a constant period spacing. The periods used here are those indicated in Table 4. 

In the text 
Fig. 3 Same as Fig. 2, but for J1735. The periods used here are those indicated in Table 9. 

In the text 
Fig. 4 Average of the computed dipole (ℓ = 1) period spacings () assessed in the range of the periods observed in J1518, corresponding to each WD model sequence considered in this work, in terms of the logarithm of the effective temperature. Numbers along each curve denote the stellar mass (in solar units). The observed period spacing derived for J1518 is depicted with a horizontal solid line. We also indicate the T_{eff} (vertical solid line) in the 3D model, together with its uncertainties (vertical dashed lines). 

In the text 
Fig. 5 Same as Fig. 4, but for J1735. 

In the text 
Fig. 6 Inverse of the quality function of the period fit considering ℓ = 1 (top panel), ℓ = 2 (middle panel), and ℓ = 1,2 (bottom panel) versus T_{eff} for J1840 (see text for details). The vertical strips depict the spectroscopic T_{eff} (solid vertical lines) and their uncertainties (dashed vertical lines) for the 1D (dark gray lines) and 3D models (red lines). 

In the text 
Fig. 7 Same as Fig. 6, but for J1112 for the five longest periods. 

In the text 
Fig. 8 Same as Fig. 7, but for the case of the seven observed periods of J1112. In the upper panel, the case of radial (ℓ = 0) and p and g modes (ℓ = 1,2) is shown, while in the bottom panel the case of p and g modes with ℓ = 1 is displayed. 

In the text 
Fig. 9 Same as Fig. 6, but for the case of J1518. 

In the text 
Fig. 10 Same as Fig. 6, but for J1614. 

In the text 
Fig. 11 Same as Fig. 6, but for J2228. 

In the text 
Fig. 12 Same as Fig. 6, but for J1738. 

In the text 
Fig. 13 Same as Fig. 6, but for J1618. 

In the text 
Fig. 14 Same as Fig. 6, but for J1735. 

In the text 
Fig. 15 Same as Fig. 6, but for J2139. 

In the text 
Fig. 16 Location of the nine analyzed stars, according to the spectroscopic parameters from the 3D modelatmosphere (black circles), and the corresponding asteroseismological models adopted for each star (white circles), along with our evolutionary tracks of lowmass Hecore WDs in the log T_{eff}−log g plane. The grayshaded region bounded by the thin blue line corresponds to the instability domain of ℓ = 1g modes according to the nonadiabatic computations of Córsico & Althaus (2016). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.