-
Non-evolutionary effects on period change in Magellanic Cepheids I. New binary systems revealed from Light Travel Time Effect
Authors:
Rajeev Singh Rathour,
Gergely Hajdu,
Radosław Smolec,
Paulina Karczmarek,
Vincent Hocdé,
Oliwia Ziółkowska,
Igor Soszyński,
Andrzej Udalski
Abstract:
Period change studies give a window to probe into the evolution and dynamics of Cepheids. While evolutionary period changes have been well studied both observationally and theoretically, non-evolutionary period changes lack a systematic and quantitative description. The overall objective is to have a quantitative understanding of the full picture of non-evolutionary period changes in Cepheids, to…
▽ More
Period change studies give a window to probe into the evolution and dynamics of Cepheids. While evolutionary period changes have been well studied both observationally and theoretically, non-evolutionary period changes lack a systematic and quantitative description. The overall objective is to have a quantitative understanding of the full picture of non-evolutionary period changes in Cepheids, to develop a formalism to disentangle it from the secular evolutionary period change. In the first part of the series of works, we aim to conduct a systematic search for non-evolutionary period changes to search for Cepheids in likely binary configuration and quantify their incidence rates in the Magellanic Clouds. We collect more than decade-long time-series photometry from the publically available survey, Optical Gravitational Lensing Experiment (OGLE), with more than 7200 Cepheids collectively from the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). Our sample contains both fundamental-mode and first overtone-mode Cepheids. Then we calculate observed minus calculated ($O-C$) diagrams to reveal the light-travel time effect (LTTE). In our search, out of an overall sample of more than 7200 Cepheids, we found 52 candidate Cepheid binary systems in the LMC (30 fundamental and 22 first overtone-mode) and 145 in the SMC (85 fundamental and 60 first overtone-mode). The majority of the sample is characterized by orbital periods of 2000-4000\,d and eccentricities of 0.2-0.5. Moreover, we report two candidates in each galaxy with the Cepheid likely existing with a giant companion. The incidence rate ratio for SMC to LMC calculated from our sample is in agreement with binary Cepheid population synthesis predictions.
△ Less
Submitted 27 March, 2024; v1 submitted 20 March, 2024;
originally announced March 2024.
-
Pulsation modelling of the Cepheid Y Ophiuchi with RSP/MESA. Impact of the circumstellar envelope and a high projection factor on Baade-Wesselink method
Authors:
V. Hocdé,
R. Smolec,
P. Moskalik,
R. Singh Rathour,
O. Ziółkowska
Abstract:
Y~Ophiuchi (Y~Oph) is a classical Cepheid reported to be as dim as a Cepheid of about half its pulsation period, and exhibits a low radial velocity and light-curves amplitude. Our objective is to conduct hydrodynamical pulsation modeling of Y~Oph to derive its distance and provide physical insight to its low amplitude and luminosity, constrained by an extensive set of observations. We first perfor…
▽ More
Y~Ophiuchi (Y~Oph) is a classical Cepheid reported to be as dim as a Cepheid of about half its pulsation period, and exhibits a low radial velocity and light-curves amplitude. Our objective is to conduct hydrodynamical pulsation modeling of Y~Oph to derive its distance and provide physical insight to its low amplitude and luminosity, constrained by an extensive set of observations. We first perform a linear analysis on a grid of models using hydrodynamical pulsation code \texttt{MESA-RSP} in order to find the combinations of mass, metallicity, effective temperature and luminosity resulting in linear excitation of pulsations with period of about 17$\,$days. Then, for the best combinations of stellar parameters, we perform non-linear computations to obtain the full-amplitude pulsations of these models. Last, we compare the results to a complete set of observations along the pulsation cycle. We adjust simultaneously the distance, the color excess and circumstellar envelope (CSE) model to fit the light curves and the angular diameter. We find that all pulsation models at high effective temperature are in remarkable agreement with the observations along the pulsation cycle. This result suggests that the low amplitude of Y~Oph can be explained by its location close to the blue edge of the instability strip. We also find that a pulsational mass of about 7-8$\,\mathrm{M}_\odot$ is consistent with a non-canonical evolutionary model with moderate overshooting, PL relation and \textit{Gaia} parallax. However, a much lower mass below 5$\,$M$_\odot$ is required to match Baade-Wesselink (BW) distance measurements from the literature. We show that the combination of the impact of the CSE on the photometry together with a projection factor of about 1.5 explains the discrepant distance and luminosity values obtained from BW methods.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Precise Fourier parameters of Cepheid Radial Velocity Curves
Authors:
V. Hocdé,
P. Moskalik,
N. A. Gorynya,
R. Smolec,
R. Singh Rathour,
O. Ziółkowska
Abstract:
The primary goal of this paper is to derive precise Fourier parameters of the radial velocity (RV) curves for fundamental and first-overtone Galactic Cepheids. For each star, we carefully selected RV measurements available in the literature which yield the highest precision of Fourier parameters. We performed a Fourier decomposition of the RV curves. We subtracted RV modulations caused by binary m…
▽ More
The primary goal of this paper is to derive precise Fourier parameters of the radial velocity (RV) curves for fundamental and first-overtone Galactic Cepheids. For each star, we carefully selected RV measurements available in the literature which yield the highest precision of Fourier parameters. We performed a Fourier decomposition of the RV curves. We subtracted RV modulations caused by binary motion and have removed other residual trends. Finally, we have displayed and analyzed qualitatively the progressions found for Fourier parameters. We applied a standard identification of their pulsation mode based on their Fourier phase $φ_{21}$. Our final sample includes 178 fundamental-mode and 33 first overtone pulsators, as well as 7 additional Cepheids whose pulsation mode is uncertain or undetermined according to our criteria. For the fundamental-mode Cepheids, we improved the precision of Fourier parameters in comparison of previous results from the literature. We are able to firmly identify V495 Cyg as a new first-overtone Cepheid. We confirm first-overtone nature of several other stars. We also show that $α$ UMi should be definitely classified as a first-overtone pulsator. In 3 objects (VY Per, AQ Pup and QZ Nor) we found significant $γ$-velocity variations that we attribute to spectroscopic binarity. Finally, the analysis of the F-mode Fourier parameters up to 7th order reveals tight progression of Fourier phases. For $P<10\,$day we find a well defined upper limit for the Fourier amplitude $A_1$. The pulsation period coverage and the precision obtained, in particular for Fourier phase $φ_{21}$, will be useful for studying the dynamics of Cepheid pulsations with the help of hydrodynamical models. Further radial velocity measurements from modern high-resolution spectroscopic instruments will be important to improve these results.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
Metallicity estimation of MW, SMC and LMC classical Cepheids from the shape of the $V$- and $I$-band light curves
Authors:
V. Hocdé,
R. Smolec,
P. Moskalik,
O. Ziółkowska,
R. Singh Rathour
Abstract:
Estimating metallicity of classical Cepheids is of prime importance for studying metallicity effect on stellar evolution, chemical evolution of galaxies, and ultimately its impact on period-luminosity relation used in the extragalactic distance scale. We aim at establishing new empirical relations for estimating the iron content of classical Cepheids for short and long-periods based on Fourier par…
▽ More
Estimating metallicity of classical Cepheids is of prime importance for studying metallicity effect on stellar evolution, chemical evolution of galaxies, and ultimately its impact on period-luminosity relation used in the extragalactic distance scale. We aim at establishing new empirical relations for estimating the iron content of classical Cepheids for short and long-periods based on Fourier parameters from the $V$-band light curves. We calibrate new interrelations of Fourier parameters to convert $V$-band empirical relations into the $I$-band. Then we apply these relation in $V$ and $I$-bands to Cepheids from Milky Way (MW), Small and Large Magellanic Clouds (SMC and LMC) available in the literature. Last, we map the metallicity distribution in these galaxies for investigating potential application in galactic archeology. These empirical relations in $V$ and $I$ bands are able to derive the mean metallicity of a sample of MW, SMC and LMC Cepheids in agreement with literature values within 1$σ$. We also show that these relations are precise enough to reconstruct the radial metallicity gradients within the MW from OGLE data. The empirical relations in the $V$ and $I$ bands calibrated in this paper for short and long-period Cepheids provide a new useful tool to estimate the metallicity of Cepheids which are not accessible by spectroscopy. The calibration can be improved with further high-resolution spectroscopic observations of metal-poor Cepheids and homogeneous photometries in $V$ and $I$ bands.
△ Less
Submitted 10 February, 2023; v1 submitted 31 December, 2022;
originally announced January 2023.
-
First Overtone Cepheids of the OGLE Magellanic Cloud Collection -- beyond radial modes
Authors:
R. Smolec,
O. Ziółkowska,
M. Ochalik,
M. Śniegowska
Abstract:
We have analysed Optical Gravitational Lensing Experiment photometry for first overtone classical Cepheids in the Large and Small Magellanic Clouds in search for additional periodicities beyond radial modes. We have used standard consecutive prewhitening technique in some cases followed by time-dependent prewhitening. We report new candidates for double-mode radial pulsations. However, majority of…
▽ More
We have analysed Optical Gravitational Lensing Experiment photometry for first overtone classical Cepheids in the Large and Small Magellanic Clouds in search for additional periodicities beyond radial modes. We have used standard consecutive prewhitening technique in some cases followed by time-dependent prewhitening. We report new candidates for double-mode radial pulsations. However, majority of signals we have detected cannot be interpreted in terms of radial modes. We report 516 double-periodic stars with period ratios, $P_{\rm x}/P_{\rm 1O}$, in the range 0.60 and 0.65. We study the properties of this class and implications for model explaining these periodicities. We also report 28 stars in which additional variability is of longer period, below radial fundamental mode, with median period ratio, $P_{\rm 1O}/P_{\rm x}$, of 0.684. This class is an analogue of a class known in RR Lyrae stars. Hundreds of other signals were detected that cannot be attributed to radial modes or the above-mentioned classes. Statistical properties of these signals are analysed. We suggest that majority of these signals correspond to non-radial modes. In particular, a significant fraction can be attributed to non-radial modes of moderate degrees, tightly connected to a class with period ratios in between 0.60 and 0.65. In tens of stars, close to radial mode frequency, relatively large-amplitude and coherent signals are observed, that may represent yet another class. In 27 stars periodic modulation of pulsation was detected. Differences in additional frequency content between the two Clouds are discussed.
△ Less
Submitted 13 December, 2022;
originally announced December 2022.
-
Evolved eclipsing binary systems in the Galactic bulge: Precise physical and orbital parameters of OGLE-BLG-ECL-305487 and OGLE-BLG-ECL-116218
Authors:
K. Suchomska,
D. Graczyk,
C. Gałan,
O. Ziółkowska,
R. Smolec,
G. Pietrzyński,
W. Gieren,
S. Villanova,
M. Górski,
I. B. Thompson,
P. Wielgórski,
B. Zgirski,
P. Karczmarek,
B. Pilecki,
M. Taormina,
W. Narloch,
G. Hajdu,
M. Lewis,
M. Kałuszyński,
G. Rojas García
Abstract:
Our goal is to determine, with high accuracy, the physical and orbital parameters of two double-lined eclipsing binary systems, where the components are two giant stars. We also aim to study the evolutionary status of the binaries, to derive the distances towards them by using a surface brightness-colour relation, and to compare these measurements with the measurements presented by the Gaia missio…
▽ More
Our goal is to determine, with high accuracy, the physical and orbital parameters of two double-lined eclipsing binary systems, where the components are two giant stars. We also aim to study the evolutionary status of the binaries, to derive the distances towards them by using a surface brightness-colour relation, and to compare these measurements with the measurements presented by the Gaia mission. In order to measure the physical and orbital parameters of the systems, we analysed the light curves and radial-velocity curves with the Wilson-Devinney code. We used V band and I-band photometry from the OGLE catalogue and near-infrared photometry obtained with the New Technology Telescope (NTT) equipped with the SOFI instrument. The spectroscopic data were collected with the HARPS spectrograph mounted at the ESO 3.6m telescope and the MIKE spectrograph mounted at the 6.5m Clay telescope. We present the first analysis of this kind for two evolved eclipsing binary systems from the OGLE catalogue: OGLE-BLG-ECL-305487 and OGLE-BLG-ECL-116218. The masses of the components of OGLE-BLG-ECL-305487 are $M_1$ = 1.059 $\pm$ 0.019 and $M_2$ = 0.991 $\pm$ 0.018 $M_\odot$, and the radii are $R_1$ = 19.27 $\pm$ 0.28 and $R_2$ = 29.99 $\pm$ 0.24 R$_\odot$. For OGLE-BLG-ECL-116218, the masses are $M_1$= 0.969 $\pm$ 0.012 and $M_2$= 0.983 $\pm$ 0.012 $M_\odot$, while the radii are $R_1$= 16.73 $\pm$ 0.28 and $R_2$= 22.06 $\pm$ 0.26 $R_\odot$. The evolutionary status of the systems is discussed based on the PARSEC and MIST isochrones. The ages of the systems were established to be between 7.3-10.9 Gyr for OGLE-BLG-ECL-305487 and around 10 Gyr for OGLE-BLG-ECL-116218. We also measured the distances to the binaries. For OGLE-BLG-ECL-305487, $d$= 7.80 $\pm$ 0.18 (stat.) $\pm$ 0.19 (syst.) kpc and for OGLE-BLG-ECL-116218, $d$= 7.57 $\pm$ 0.28 (stat.) $\pm$ 0.19 (syst.) kpc.
△ Less
Submitted 29 October, 2022; v1 submitted 22 October, 2022;
originally announced October 2022.
-
Early recognition of Microlensing Events from Archival Photometry with Machine Learning Methods
Authors:
I. Gezer,
Ł. Wyrzykowski,
P. Zieliński,
G. Marton,
K. Kruszyńska,
K. A. Rybicki,
N. Ihanec,
M. Jabłońska,
O. Ziółkowska
Abstract:
Gravitational microlensing method is a powerful method to detect isolated black holes in the Milky Way. During a microlensing event brightness of the source increases and this feature is used by many photometric surveys to alert on potential events. A typical microlensing event shows a characteristic light curve, however, some outbursting variable stars may show similar light curves to microlensin…
▽ More
Gravitational microlensing method is a powerful method to detect isolated black holes in the Milky Way. During a microlensing event brightness of the source increases and this feature is used by many photometric surveys to alert on potential events. A typical microlensing event shows a characteristic light curve, however, some outbursting variable stars may show similar light curves to microlensing events especially when the cadence of observations is not dense enough. Our aim is to device a method for distinguishing candidates for microlensing events from any other types of alerts using solely their archival photometric multi-wavelength data. The most common contaminants in the microlensing event search are Classical Be-type stars, Young Stellar Objects and Asymptotic Giant Branch stars such as Miras. We build a training set using thousands of examples for the main classes of alerting stars combining optical to mid-infrared magnitudes from Gaia, 2MASS and AllWISE catalogues. We used supervised machine learning techniques to build models for classification of alerts. We verified our method on 120 microlensing events reported by Gaia Science Alerts which were studied spectroscopically and photometrically. With the use of only archival information at 90% probability threshold we correctly identified one-third of the microlensing events. We also run our classifier on positions of 368 Gaia alerts which were flagged as potential candidates for microlensing events. At the 90% probability threshold we classified 38 microlensing events and 29 other types of variables. The machine learning supported method we developed can be universally used for current and future alerting surveys in order to quickly assess the classification of galactic transients and help decide on further follow-up observations.
△ Less
Submitted 5 February, 2024; v1 submitted 28 January, 2022;
originally announced January 2022.
-
Single-lens mass measurement in the high-magnification microlensing event Gaia19bld located in the Galactic disc
Authors:
K. A. Rybicki,
Ł. Wyrzykowski,
E. Bachelet,
A. Cassan,
P. Zieliński,
A. Gould,
S. Calchi Novati,
J. C. Yee,
Y. -H. Ryu,
M. Gromadzki,
P. Mikołajczyk,
N. Ihanec,
K. Kruszyńska,
F. -J. Hambsch,
S. Zoła,
S. J. Fossey,
S. Awiphan,
N. Nakharutai,
F. Lewis,
F. Olivares E.,
S. Hodgkin,
A. Delgado,
E. Breedt,
D. L. Harrison,
M. vanLeeuwen
, et al. (44 additional authors not shown)
Abstract:
We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up ob…
▽ More
We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up observations, which allowed us to determine all parameters vital for the characterisation of the lens and the source in the microlensing event. Gaia19bld was discovered by the Gaia satellite and was subsequently intensively followed up with a network of ground-based observatories and the Spitzer Space Telescope. We collected multiple high-resolution spectra with Very Large Telescope (VLT)/X-Shooter to characterise the source star. The event was also observed with VLT Interferometer (VLTI)/PIONIER during the peak. Here we focus on the photometric observations and model the light curve composed of data from Gaia, Spitzer, and multiple optical, ground-based observatories. We find the best-fitting solution with parallax and finite source effects. We derived the limit on the luminosity of the lens based on the blended light model and spectroscopic distance. We compute the mass of the lens to be $1.13 \pm 0.03~M_{\odot}$ and derive its distance to be $5.52^{+0.35}_{-0.64}~\mathrm{kpc}$. The lens is likely a main sequence star, however its true nature has yet to be verified by future high-resolution observations. Our results are consistent with interferometric measurements of the angular Einstein radius, emphasising that interferometry can be a new channel for determining the masses of objects that would otherwise remain undetectable, including stellar-mass black holes.
△ Less
Submitted 2 December, 2021;
originally announced December 2021.
-
Lens parameters for Gaia18cbf -- a long gravitational microlensing event in the Galactic plane
Authors:
Katarzyna Kruszyńska,
Ł. Wyrzykowski,
K. A. Rybicki,
M. Maskoliūnas,
E. Bachelet,
N. Rattenbury,
P. Mróz,
P. Zieliński,
K. Howil,
Z. Kaczmarek,
S. T. Hodgkin,
N. Ihanec,
I. Gezer,
M. Gromadzki,
P. Mikołajczyk,
A. Stankevičiūtė,
V. Čepas,
E. Pakštienė,
K. Šiškauskaitė,
J. Zdanavičius,
V. Bozza,
M. Dominik,
R. Figuera Jaimes,
A. Fukui,
M. Hundertmark
, et al. (7 additional authors not shown)
Abstract:
Context: The timescale of a microlensing event scales as a square root of a lens mass. Therefore, long-lasting events are important candidates for massive lenses, including black holes.
Aims: Here we present the analysis of the Gaia18cbf microlensing event reported by the Gaia Science Alerts system. It exhibited a long timescale and features that are common for the annual microlensing parallax e…
▽ More
Context: The timescale of a microlensing event scales as a square root of a lens mass. Therefore, long-lasting events are important candidates for massive lenses, including black holes.
Aims: Here we present the analysis of the Gaia18cbf microlensing event reported by the Gaia Science Alerts system. It exhibited a long timescale and features that are common for the annual microlensing parallax effect. We deduce the parameters of the lens based on the derived best fitting model.
Methods: We used photometric data collected by the Gaia satellite as well as the follow-up data gathered by the ground-based observatories. We investigated the range of microlensing models and used them to derive the most probable mass and distance to the lens using a Galactic model as a prior. Using known mass-brightness relation we determined how likely it is that the lens is a main-sequence (MS) star.
Results: This event is one of the longest ever detected, with the Einstein timescale of $t_\mathrm{E}=491.41^{+128.31}_{-84.94}$ days for the best solution and $t_\mathrm{E}=453.74^{+178.69}_{-105.74}$ days for the second-best. Assuming Galaxy priors, this translates to the most probable lens mass of $M_\mathrm{L} = 2.65^{+5.09}_{-1.48} M_\odot$ and $M_\mathrm{L} = 1.71^{+3.78}_{-1.06} M_\odot$, respectively.
The limits on the blended light suggest that this event was most likely not caused by a MS star, but rather by a dark remnant of stellar evolution.
△ Less
Submitted 7 April, 2022; v1 submitted 16 November, 2021;
originally announced November 2021.
-
Gaia 18dvy: a new FUor in the Cygnus OB3 association
Authors:
E. Szegedi-Elek,
P. Ábrahám,
L. Wyrzykowski,
M. Kun,
A. Kóspál,
L. Chen,
G. Marton,
A. Moór,
Cs. Kiss,
A. Pál,
L. Szabados,
J. Varga,
E. Varga-Verebélyi,
C. Andreas,
E. Bachelet,
R. Bischoff,
A. Bódi,
E. Breedt,
U. Burgaz,
T. Butterley,
V. Čepas,
G. Damljanovic,
I. Gezer,
V. Godunova,
M. Gromadzki
, et al. (37 additional authors not shown)
Abstract:
We present optical-infrared photometric and spectroscopic observations of Gaia18dvy, located in the Cygnus OB3 association at a distance of 1.88 kpc. The object was noted by the Gaia alerts system when its lightcurve exhibited a $\gtrsim$4 mag rise in 2018-2019. The brightening was also observable at mid-infared wavelengths. The infrared colors of Gaia18dvy became bluer as the outburst progressed.…
▽ More
We present optical-infrared photometric and spectroscopic observations of Gaia18dvy, located in the Cygnus OB3 association at a distance of 1.88 kpc. The object was noted by the Gaia alerts system when its lightcurve exhibited a $\gtrsim$4 mag rise in 2018-2019. The brightening was also observable at mid-infared wavelengths. The infrared colors of Gaia18dvy became bluer as the outburst progressed. Its optical and near-infrared spectroscopic characteristics in the outburst phase are consistent with those of bona fide FU Orionis-type young eruptive stars. The progenitor of the outburst is probably a low-mass K-type star with an optical extinction of $\sim$3 mag. A radiative transfer modeling of the circumstellar structure, based on the quiescent spectral energy distribution, indicates a disk with a mass of $4{\times}10^{-3}\,M_{\odot}$. Our simple accretion disk modeling implies that the accretion rate had been exponentially increasing for more than 3 years until mid-2019, when it reached a peak value of $6.9 \times 10^{-6}\,M_{\odot}$yr$^{-1}$. In many respects, Gaia18dvy is similar to the FU Ori-type object HBC 722
△ Less
Submitted 16 June, 2020; v1 submitted 23 May, 2020;
originally announced May 2020.