Statistical structure and physical interpretation of the fillout factor distribution in contact binary stars
Abstract
Contact binary systems exhibit complex evolutionary behavior governed by mass transfer, angular momentum loss (AML), and Roche geometry. The fillout factor (f) is widely used as a descriptor of the degree of overcontact, yet its physical and statistical role remains nonuniformly defined across the literature. We used a compiled observational sample of W UMa-type contact binaries to construct a data-driven classification framework based on f and its relation to fundamental system parameters. Unsupervised clustering applied to the one-dimensional f distribution reveals a robust three-class structure corresponding to shallow-, medium-, and deep-contact configurations. The statistically determined boundaries are located at f ≃ 0.257 and f ≃ 0.561, with a stable silhouette score of S ≃ 0.613. The stability of this classification is confirmed through perturbation tests in f and parameter-conditioned samples. A random forest regression model shows that 60.7\% of the variance in f can be explained by the considered physical and geometric system parameters, with geometric parameters providing the largest predictive contribution. Among the considered parameters, the mass ratio emerges as the dominant predictor. The identified subclasses therefore appear to reflect a continuous geometric sequence rather than discrete physical states, governed by the progressive evolution of the Roche equipotential configuration under mass transfer and AML. Notably, the second boundary separating medium- and deep-contact systems shows lower stability than the first boundary, suggesting that the transition toward deep contact may be more sensitive to the physical composition of the analyzed sample.
Related
- The physical and chemical structure of TMC-1: The spatial distribution and excitation conditions of carbon, oxygen, and sulfur bearing species -- Depletion of CO and CS
- The physical and chemical structure of TMC-1. The spatial distribution and excitation conditions of carbon, oxygen, and sulfur bearing species: Depletion of CO and CS
- Probing Triton's Space Environment and Internal Structure: An Integrated Detection-and-Interpretation Framework
- Effects of Physical Cooling on the Structure of Circumbinary Disks
- A model for the image structure and spectral energy distribution of relativistic jets
BibTeX
@misc{poro2026statistical,
title = {{Statistical structure and physical interpretation of the fillout factor distribution in contact binary stars}},
author = {A. Poro and K. Li and E. Paki and F. Alicavus and N. Alan},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.07287},
eprint = {2609.07287},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.07287},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Statistical structure and physical interpretation of the fillout factor distribution in contact binary stars
AU - A. Poro
AU - K. Li
AU - E. Paki
AU - F. Alicavus
AU - N. Alan
PY - 2026
DA - 2026/09/07
PB - arXiv
DO - 10.48550/arxiv.2609.07287
UR - https://arxiv.org/abs/2609.07287
AB - Contact binary systems exhibit complex evolutionary behavior governed by mass transfer, angular momentum loss (AML), and Roche geometry. The fillout factor (f) is widely used as a descriptor of the degree of overcontact, yet its physical and statistical role remains nonuniformly defined across the literature. We used a compiled observational sample of W UMa-type contact binaries to construct a data-driven classification framework based on f and its relation to fundamental system parameters. Unsupervised clustering applied to the one-dimensional f distribution reveals a robust three-class structure corresponding to shallow-, medium-, and deep-contact configurations. The statistically determined boundaries are located at f ≃ 0.257 and f ≃ 0.561, with a stable silhouette score of S ≃ 0.613. The stability of this classification is confirmed through perturbation tests in f and parameter-conditioned samples. A random forest regression model shows that 60.7\% of the variance in f can be explained by the considered physical and geometric system parameters, with geometric parameters providing the largest predictive contribution. Among the considered parameters, the mass ratio emerges as the dominant predictor. The identified subclasses therefore appear to reflect a continuous geometric sequence rather than discrete physical states, governed by the progressive evolution of the Roche equipotential configuration under mass transfer and AML. Notably, the second boundary separating medium- and deep-contact systems shows lower stability than the first boundary, suggesting that the transition toward deep contact may be more sensitive to the physical composition of the analyzed sample.
ER -