Binary Star Systems Research Articles

How the QCD trace anomaly behaves at the core of twin stars?

We investigate the behavior of the dense and cold (normalized) quantum chromodynamics (QCD) trace anomaly, Δ, in the interior of twin neutron stars (obtained from several sets of equations of state in agreement with modern compact-star and multimessenger data) satisfying static and dynamic stability conditions. We scan the formed twin-star parameter space in order to look for effects caused by the presence of a strong first-order phase transition connecting hadron and quark phases by means of a Maxwell construction. We found robustly that Δ suffers an abrupt decrease around the transition point, even reaching large negative values (Δ≃−0.35), in marked contrast to current studies pointing out a smooth behavior with Δ≳0 at all densities. Besides, we characterize the behavior of the recently defined conformal factor, dc, in the four categories of twin stars for which we perform comparisons with theoretical constraints put, e.g., by Bayesian studies sometimes adjusted to agree with perturbative QCD at high densities. These analyses allow us to hypothesize modifications in the strong QCD coupling in dense nuclear matter with a strong thermodynamic discontinuity. Published by the American Physical Society 2024

Two Novel Hot Jupiter Formation Pathways: How White Dwarf Kicks Shape the Hot Jupiter Population

Abstract The origin of Hot Jupiters (HJs) is disputed between a variety of in situ and ex situ formation scenarios. One of the early proposed ex situ scenarios was the Eccentric Kozai–Lidov (EKL) mechanism combined with tidal circularization, which can produce HJs with the aid of a stellar or planetary companion. However, observations have revealed a lack of stellar companions to HJs, which challenges the importance of the binary star-driven-EKL-plus-tides scenario. In this work, we explore so far unaccounted-for stellar evolution effects on HJ formation, in particular the effect of white dwarf (WD) formation. Gaia observations have revealed that WDs often undergo a kick during formation, which can alter a binary’s orbital configuration or even unbind it. Based on this WD kick, in this Letter, we propose and explore two novel HJ formation pathways: (1) HJs that are presently orbiting single stars but were initially formed in a binary that was later unbound by a WD kick; (2) binaries that survive the WD kick can trigger enhanced EKL oscillations and lead to second-generation HJ formation. We demonstrate that the majority of seemingly single HJs could have formed in binary star systems. As such, HJ formation in binaries via the EKL mechanism could be one of the dominant HJ formation pathways, and our results highlight that unaccounted-for stellar evolution effects, like WD formation, can obscure the actual origin of observed exoplanet populations.

Characterization of the Blue Straggler Star Populations, through Statistical, Photometric, and Spectral Energy Distribution Analysis, in the Old Open Cluster: NGC 2243

We present a statistical, photometric, and spectral energy distribution (SED) analysis of the poorly studied old open cluster: NGC 2243, to characterize its blue straggler star (BSS) population. We applied ensemble-based unsupervised machine learning methods to estimate the membership probabilities using Gaia Data Release 3 (DR3) astrometric data. NGC 2243 is an open cluster that is 3.67 Gyr old with a metallicity of −0.375 dex, situated at a distance of 3.65 kpc. By analyzing the position of cluster members on the color–magnitude diagram using MIST isochrones, we have identified 12 potential BSSs in NGC 2243. We fitted the radial surface density profile and investigated the dynamical state and mass segregation effect of the cluster. It is found that the BSSs are significantly concentrated within the central region. We used data from Swift/UVOT, Gaia DR3, Pan-STARRS1 DR2, 2MASS, and WISE to fit the SEDs of the 12 identified BSSs using VOSA. We estimated the masses of the BSSs from the Hertzsprung–Russell diagram and found that they ranged from 1.25 to 2.22 M ☉. Consequently, we concluded that the BSSs likely gained 0.11–1.08 M ☉ through the mass transfer or merger channels. We discovered a hot companion associated with one BSS candidate, which has a temperature of 19,000 K, a luminosity of 0.55 L ☉, and a radius of 0.065 R ☉. The hot companion is probably a white dwarf, with its mass estimated to be approximately 0.18–0.20 M ☉ and an age of 186 Myr, suggesting it is a post-mass-transfer (Case A or Case B) system.

Realistic binary neutron star initial data with Elliptica

Abstract This work introduces the Elliptica pseudo-spectral code for generating initial data of binary neutron star systems. Building upon the recent Elliptica code update, we can now construct initial data using not only piecewise polytropic equations of state, but also tabulated equations of state for these binary systems. Furthermore, the code allows us to endow neutron stars within the binary system with spins. These spins can have a magnitude close to the mass shedding limit and can point in any direction.

2019odp -- A massive oxygen-rich Type Ib supernova

Stripped envelope (SE) supernovae are explosions of stars that have somehow lost most of their outer envelopes. We present the discovery and analyse the observations of the supernova 2019odp (a.k.a. ZTF19abqwtfu) covering epochs within days of the explosion to late nebular phases at 360\,d post-explosion. Our observations include an extensive set of photometric observations and low- to medium-resolution spectroscopic observations, both covering the complete observable time range. We analysed the data using analytic models for the recombination cooling emission of the early excess emission and the diffusion of the peak light curve. We expanded on existing methods to derive oxygen mass estimates from nebular phase spectroscopy, and briefly discuss progenitor models based on this analysis. Our spectroscopic observations confirm the presence of He in the supernova ejecta and we thus (re)classify SN 2019odp as a supernova . From the pseudo-bolometric light curve, we estimate a high ejecta mass of $M_ ej 4 - 7 M_ The high ejecta mass, large nebular O I Ca II $ line flux ratio ($1.2 - 1.9$), and an oxygen mass above $ 0.5\, M_ point towards a progenitor with a pre-explosion mass higher than $18\,M_ Whereas a majority of analysed SE supernovae in the literature seem to have low ejecta masses, indicating stripping in a binary star system, SN 2019odp instead has parameters that are consistent with an origin in a single massive star. The compact nature of the progenitor ($ 10\,R_ suggests that a Wolf-Rayet star is the progenitor.

The First Catalog of Candidate White Dwarf–Main-sequence Binaries in Open Star Clusters: A New Window into Common Envelope Evolution

Close binary systems are the progenitors to both Type Ia supernovae and the compact object mergers that can be detected via gravitational waves. To achieve a binary with a small radial separation, it is believed that the system likely undergoes common envelope (CE) evolution. Despite its importance, CE evolution may be one of the largest uncertainties in binary evolution due to a combination of computational challenges and a lack of observed benchmarks where both the post-CE and pre-CE conditions are known. Identifying post-CE systems in star clusters can partially circumvent this second issue by providing an independent age constraint on the system. For the first time, we conduct a systematic search for white dwarf and main-sequence binary systems in 299 Milky Way open star clusters. Coupling Gaia DR3 photometry and kinematics with multiband photometry from Pan-STARRS1 and the Two Micron All Sky Survey, we apply a machine learning-based approach and find 52 high-probability candidates in 38 open clusters. For a subset of our systems, we present follow-up spectroscopy from the Gemini and Lick Observatories and archival light curves from the Transiting Exoplanet Survey Satellite, Kepler/K2, and the Zwicky Transient Facility. Examples of M dwarfs with hot companions are spectroscopically observed, along with regular system variability. While the kinematics of our candidates are consistent with their host clusters, some systems have spatial positions offset relative to their hosts, potentially indicative of natal kicks. Ultimately, this catalog is a first step to obtaining a set of observational benchmarks to better link post-CE systems to their pre-CE progenitors.

Quantitative spectroscopy of multiple OB stars

Context. The majority of massive stars are located in binary or multiple star systems. Compared to single stars, these objects pose additional challenges to quantitative analyses based on model atmospheres. In particular, little information is currently available on the chemical composition of such systems. Aims. The members of the quadruple star system HD 37061, which excites the H II region Messier 43 in Orion, are fully characterised. Accurate and precise abundances for all elements with lines traceable in the optical spectrum are derived for the first time. Methods. A hybrid non-local thermodynamic equilibrium (non-LTE) approach, using line-blanketed hydrostatic model atmospheres computed with the ATLAS12 code in combination with non-LTE line-formation calculations with DETAIL and SURFACE, was employed. A high-resolution composite spectrum was analysed for the atmospheric parameters and elemental abundances of the individual stars. Fundamental stellar parameters were derived based on stellar evolution tracks, and the interstellar reddening was characterised. Results. We determined the fundamental parameters and chemical abundances for three stars in the HD 37061 system. The fourth and faintest star in the system shows no distinct spectral features, as a result of its fast rotation. However, this star has noticeable effects on the continuum. The derived element abundances and determined ages of the individual stars are consistent with each other, and the abundances coincide with the cosmic abundance standard. We find an excellent agreement between our spectroscopic distance and the Gaia Data Release 3 parallax distance.

The Possible Factors that Influence the Stability of Binary Star System

Binary star systems are a crucial subject in astrophysics, as they provide insights into stellar formation, evolution, and dynamics. Despite extensive research, there remains a gap in understanding the factors influencing the longterm stability of such systems. This paper discusses the factors that may affect the stability of binary star systems by analyzing and exploring the physical laws and formulas of the operation of binary star systems and using the operation data of specific binary star systems as evidence. It can be concluded that the factors affecting the stability of a binary star system may be as follows: 1. Mass ratio of two stars, 2. Stellar spacing, 3. Orbital eccentricity, 4. Influence of stellar evolution stage, 5. The influence of external bodies. Future research should further explore the combined effects of these factors in more complex systems to better understand binary star stability in various astrophysical environments. Additionally, understanding how these factors interact in multi-body systems could provide deeper insights into the evolution of complex stellar systems.

A Statistical, Photometric, and SED Analysis to Characterize the BSS Populations in Old Open Cluster: Berkeley 39

ABSTRACT We present a statistical, photometric, and spectral energy distribution (SED) analysis to characterize the blue straggler stars (BSSs) populations in the Galactic old open cluster Berkeley 39. Berkeley 39 is a 6.16 Gyr old open cluster located at a distance of 3.99 Kpc.Gaia DR3 astrometry data have been used to estimate the membership probabilities using ensemble-based unsupervised machine learning techniques. We identified 21 BSS candidates on the colour–magnitude diagram, with 19 of them being detected in the Swift/UVOT UVW2 filter. We analysed the radial surface density profile and examined the cluster dynamical states and mass segregation effect. The SEDs of 19 BSSs are constructed using multiwavelength data covering UV to IR wavelengths. A single-component SED is fitted successfully for 14 BSS candidates. We discovered hot companions in five BSS candidates. These hot companions have temperatures of approximately 14 000 to 23 000 K, radii ranging from 0.04 to 0.13 R$_{\odot }$, and luminosities ranging from 0.16 to 2.91 L$_{\odot }$. Among these, three are most likely extremely low-mass white dwarfs (WDs) with masses around 0.17 to 0.18 M$_{\odot }$, and two are low-mass WDs with masses around 0.18 to 0.39 M$_{\odot }$. This confirms that they are post-mass transfer (Case A or Case B) systems. We also investigated the variable characteristics of BSSs by analysing their light curves using data from TESS. Our analysis confirms that two BSSs identified as eclipsing binaries in Gaia DR3 are indeed eclipsing binaries. Additionally, one of the two eclipsing binary BSSs shows evidence of having hot companions, as indicated by the multiwavelength SEDs.

Five new eclipsing binaries with low-mass companions

Precise space-based photometry from the Transiting Exoplanet Survey Satellite results in a huge number of exoplanetary candidates. However, the masses of these objects are unknown and must be determined by ground-based spectroscopic follow-up observations, frequently revealing the companions to be low-mass stars rather than exoplanets. We present the first orbital and stellar parameter solutions for five such eclipsing binary-star systems using radial-velocity follow-up measurements together with spectral-energy-distribution solutions. TOI-416 and TOI-1143 are totally eclipsing F+M star systems with well-determined secondary masses, radii, and temperatures. TOI-416 is a circular system with an F6 primary and a secondary with a mass of M2 = 0.131(8) M⊙. TOI-1143 consists of an F6 primary with an M2 = 0.142(3) M⊙ secondary on an eccentric orbit with a third companion. With respect to the other systems, TOI-1153 shows ellipsoidal variations, TOI-1615 contains a pulsating primary, and TOI-1788 has a spotted primary, while all have moderate mass ratios of 0.2–0.4. However, these systems are in a grazing configuration, which limits their full description. The parameters of TOI-416B and TOI-1143B are suitable for the calibration of the radius-mass relation for dwarf stars.

Exotic Stable Branches with Efficient TOV Sequences

Modern inference schemes for the neutron star (NS) equation of state (EoS) require large numbers of stellar models constructed with different EoS, and these stellar models must capture all the behavior of stable stars. I introduce termination conditions for sequences of stellar models for cold nonrotating NSs that can identify all stable stellar configurations up to arbitrarily large central pressures along with an efficient algorithm to build accurate interpolators for macroscopic properties. I explore the behavior of stars with both high- and low-central pressures. Interestingly, I find that EoSs with monotonically increasing sound-speed can produce multiple stable branches (twin stars) and that large phase transitions at high densities can produce stable branches at nearly any mass scale, including sub-solar masses while still supporting stars with M > 2 M ⊙. I conclude with some speculation about the astrophysical implications of this behavior.

Many body gravity and the galaxy rotation curves

A novel theory was proposed earlier to model systems with thermal gradients, based on the postulate that the spatial and temporal variation in temperature can be recast as a variation in the metric. Combining the variation in the metric due to the thermal variations and gravity, leads to the concept of thermal gravity in a 5-D space-time-temperature setting. When the 5-D Einstein field equations are projected on to a 4-D space, they result in additional terms in the field equations. This may lead to unique phenomena such as the spontaneous symmetry breaking of scalar particles in the presence of a strong gravitational field. This theory, originally conceived in a quantum mechanical framework, is now adapted to explain the galaxy rotation curves. A galaxy is not in a state of thermal equilibrium. A parameter called the “degree of thermalization” is introduced to model partially thermalized systems. The generalization of thermal gravity to partially thermalized systems, leads to the theory of many-body gravity. The theory of many-body gravity is now shown to be able to explain the rotation curves of the Milky Way and the M31 (Andromeda) galaxies, to a fair extent. The radial acceleration relation (RAR) for 63 galaxies, with their galactic masses spanning three orders of magnitude, has been replicated. Finally, the wide binary star (WBS) system is touched upon.

Revising Properties of Planet–Host Binary Systems. IV. The Radius Distribution of Small Planets in Binary Star Systems Is Dependent on Stellar Separation* *Based on observations obtained with the Hobby–Eberly Telescope, which is a joint project of the University of Texas at Austin, the Pennsylvania State University, Ludwig-Maximilians-Universität München, and Georg-August-Universität Göttingen.

Small planets (R p ≤ 4 R ⊕) are divided into rocky super-Earths and gaseous sub-Neptunes separated by a radius gap, but the mechanisms that produce these distinct planet populations remain unclear. Binary stars are the only main-sequence systems with an observable record of the protoplanetary disk lifetime and mass reservoir, and the demographics of planets in binaries may provide insights into planet formation and evolution. To investigate the radius distribution of planets in binary star systems, we observed 207 binary systems hosting 283 confirmed and candidate transiting planets detected by the Kepler mission, then recharacterized the planets while accounting for the observational biases introduced by the secondary star. We found that the population of planets in close binaries (ρ ≤ 100 au) is significantly different from the planet population in wider binaries (ρ > 300 au) or single stars. In contrast to planets around single stars, planets in close binaries appear to have a unimodal radius distribution with a peak near the expected super-Earth peak of R p ∼ 1.3 R ⊕ and a suppressed population of sub-Neptunes. We conclude that we are observing the direct impact of a reduced disk lifetime, smaller mass reservoir, and possible altered distribution of solids reducing the sub-Neptune formation efficiency. Our results demonstrate the power of binary stars as a laboratory for exploring planet formation and as a controlled experiment of the impact of varied initial conditions on mature planet populations.

Exploring pathways to forming twin stars

Exploring pathways to forming twin stars

Short-period post-common envelope binaries with Balmer emission from SDSS and LAMOST based on ZTF photometric data

ABSTRACT We present here 55 short-period post-common envelope binaries (PCEBs) containing a hot white dwarf (WD) and a low-mass main sequence (MS). Based on the photometric data from Zwicky Transient Facility survey data Release 19 (ZTF DR19), the light curves are analysed for about 200 WDMS binaries with emission line(s) identified from the Sloan Digital Sky Survey (SDSS) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) spectra, in which 55 WDMS binaries are found to exhibit variability in their luminosities with a short period and are thus short-period binaries (i.e. PCEBs). In addition, it is found that the orbital periods of these PCEBs locate in a range from 2.2643 to 81.1526 h. However, only six short-period PCEBs are newly discovered and the orbital periods of 19 PCEBs are improved in this work. Meanwhile, it is found that three objects are newly discovered eclipsing PCEBs, and a object (i.e. SDSS J1541) might be the short-period PCEB with a late M-type star or a brown dwarf companion based on the analysis of its spectral energy distribution. At last, the mechanism(s) being responsible for the emission features in the spectra of these PCEBs are discussed, the emission features arising in their optical spectra might be caused by the stellar activity or an irradiated component owing to a hot WD companion because most of them contain a WD with an effective temperature higher than $\sim$10 000 K.

Signatures of hadron-quark phase transition through the r-mode instability in twin stars

Signatures of hadron-quark phase transition through the r-mode instability in twin stars

UOCS XIV: Study of the Open Cluster NGC 2627 Using UVIT/AstroSat

We study the intermediate-age open cluster NGC 2627, located at a distance of ∼2 kpc, using UVIT/AstroSat and other archival data. Using a machine learning-based algorithm, ML-MOC, on the Gaia DR3 data, we identify 422 cluster members, including four blue straggler stars (BSSs), one yellow straggler star (YSS), one blue lurker (BL), one red clump (RC) star, and two binary candidates with detection in both UVIT/F148W and UVIT/F169M filters. We characterise them using multiwavelength spectral energy distributions (SEDs). Out of the above nine sources, one BSS, the BL, and one binary candidate have a source nearby; hence, we did not fit their SEDs. Of the remaining six sources, we successfully fit two with single-component SEDs and four with binary-component SEDs. The binary-component SED-based parameters indicate that the hot companions of BSSs, the YSS, the RC star, and the binary candidate are extremely low-mass white dwarfs, confirming that at least four out of nine stars (44%) are formed via the mass transfer channel. We fit King’s profile function to the high-probability (p > 0.8) cluster members and estimate the cluster core radius (r C ) to be 3.84′ and the tidal radius (r t ) to be 36.85′. We find that the equal-mass binaries are most concentrated towards the cluster center, followed by the single massive stars, and single low-mass stars. The BSS population of the cluster is also found to be located within a radius r ∼ 10 × r C from the cluster center, suggesting the dynamical evolution of the cluster.

Exoplanet Aeronomy: A Case Study of WASP-69 b’s Variable Thermosphere

Aeronomy, the study of Earth’s upper atmosphere and its interaction with the local space environment, has long traced changes in the thermospheres of Earth and other solar system planets to solar variability in the X-ray and extreme-ultraviolet (collectively, XUV) bands. Extending comparative aeronomy to the short-period extrasolar planets may illuminate whether stellar XUV irradiation powers atmospheric outflows that change planetary radii on astronomical timescales. In recent years, near-IR transit spectroscopy of metastable Hei has been a prolific tracer of high-altitude planetary gas. We present a case study of exoplanet aeronomy using metastable Hei transit observations from Palomar Observatory's Wide Field InfraRed Camera and follow-up high-energy data from the Neil Gehrels Swift Observatory that were taken within 1 month of the WASP-69 system, a K-type main-sequence star with a well-studied hot Jupiter companion. Supplemented by archival data, we find that WASP-69's X-ray flux in 2023 was less than 50% of what was recorded in 2016 and that the metastable Hei absorption from WASP-69 b was lower in 2023 versus past epochs from 2017 to 2019. Via atmospheric modeling, we show that this time-variable metastable Hei signal is in the expected direction given the observed change in stellar XUV, possibly stemming from WASP-69's magnetic activity cycle. Our results underscore the ability of multiepoch, multiwavelength observations to paint a cohesive picture of the interaction between an exoplanet’s atmosphere and its host star.

Two dynamically discovered compact object candidate binary systems from LAMOST low-resolution survey

ABSTRACT We report two binary systems, LAMOST J035540 + 381550 (hereafter J035540) and LAMOST J035916 + 400732 (hereafter J035916), identified through the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS). Each of these two systems contains an M-type star orbiting with a invisible compact object candidate. Follow-up spectroscopic observations of Palomar 200-in. telescope (P200) enhance radial velocity measurements. We use radial velocities from LAMOST and P200, as well as light curves from Zwicky Transient Facility (ZTF) to constrain orbital parameters. The masses of the visible M-type stars are estimated by fitting the MIST isochrones and spectral energy distributions. The mass functions for the unseen companions are: $0.22\pm 0.01\,{\rm M}_{\odot }$ for J035540 and $0.16\pm 0.01\, {\rm M}_{\odot }$ for J035916. With the orbital and stellar parameters derived above and assuming the orbital inclination is 90° (edge-on), we find that the minimum masses of the invisible companions exceeds that of the visible stars. The single-lined feature and the dynamical evidence suggest the presence of compact objects. J035540’s ZTF light curve, modelled with phoebe, yields a compact object mass of $0.70^{+0.12}_{-0.05}\, {\rm M}_{\odot }$. For J035916, ellipsoidal modulation analysis constrains the light-curve amplitude, yielding a compact object mass range of $0.57-0.90\, {\rm M}_{\odot }$. The mass estimates indicate that both are likely white dwarfs. These findings underscore the efficiency of optical time-domain surveys and dynamical methods in identifying faint, massive white dwarfs, along with other compact objects in binaries.

The DBL Survey I: discovery of 34 double-lined double white dwarf binaries

ABSTRACT We present the first discoveries of the double-lined double white dwarf (DBL) survey that targets overluminous sources with respect to the canonical white dwarf cooling sequence according to a set of well-defined criteria. The primary goal of the DBL survey is to identify compact double white dwarf binary star systems from a unique spectral detection of both stars, which then enables a precise quantification of the atmospheric parameters and radial velocity variability of a system. Our search of 117 candidates that were randomly selected from a magnitude-limited sample of 399 yielded a 29 per cent detection efficiency with 34 systems exhibiting a double-lined signature. A further 38 systems show strong evidence of being single-lined or potentially DBL binaries and seven single-lined sources from the full observed sample are radial velocity variable. The 45 remaining candidates appear as a single WD with no companion or a non-DA white dwarf, bringing the efficiency of detecting binaries to 62 per cent. Atmospheric fitting of all double-lined systems reveals a large fraction that have two similar mass components that combine to a total mass of 1.0–1.3 $\mathrm{M}_\odot$ – a class of double white dwarf binaries that may undergo a sub-Chandrasekhar mass type Ia detonation or merge to form a massive O/Ne WD, although orbital periods are required to infer on which time-scales. One double-lined system located 49 pc away, WDJ181058.67+311940.94, is super-Chandrasekhar mass, making it the second such double white dwarf binary to be discovered.