Locally authored papers of the past 5 days

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this "back-reaction" of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

Internal feedback from massive stars and active galactic nuclei (AGN) play a key role in galaxy evolution, but external environmental effects can also strongly influence galaxies. We investigate the impact of environment on galaxy evolution, and its dependence on baryonic physics implementation, using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) spanning a wide range of stellar and AGN feedback implementations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. We show that satellite galaxies are significantly affected by the environment in all simulation models, with their gas fraction and star formation rate (SFR) suppressed in overdense regions compared to similar mass satellites in underdense environments at $z=0$. Central galaxies are less sensitive to environment but tend to show lower gas fraction and SFR in overdense regions at low stellar mass, transitioning to higher gas fraction and SFR for massive galaxies in higher-density environments. Halo baryon fraction ($f_{\rm B}$) and circumgalactic medium mass fraction ($f_{\rm CGM}$) at $z=0$ show clear environmental effects. In SIMBA, low-mass haloes in overdense regions have systematically lower $f_{\rm B}$ and $f_{\rm CGM}$ at fixed halo mass, while Swift-EAGLE haloes in overdense regions have systematically higher $f_{\rm B}$ and $f_{\rm CGM}$ across the full halo mass range, and IllustrisTNG and ASTRID show opposite trends at the low and high mass ends. Environmental effects can flip at higher redshift, with SFR and $f_{\rm B}$ increasing with local density in low-mass haloes before quenching at an increasing overdensity threshold. Our results demonstrate that the impact of environment on galaxy evolution depends significantly on galaxy formation model, and higher-density environments can either suppress or enhance star formation depending on galaxy mass and cosmic epoch.

Using deep medium-resolution JWST rest-optical spectra of a sample of typical star-forming galaxies (Lyman break galaxies and Lyman-$\alpha$ emitters) from the LyC22 survey at $z \sim 3$, we determined the nebular abundances of N, O, and Ne relative to H for a subsample of 25 objects with the direct method, based on auroral [OIII]4363 line detections. Our measurements increases the number of accurate N/O determinations at $z \sim 2-4$ using a homogeneous approach. We found a mean value of $\log({\rm N/O})=-1.29^{+0.25}_{-0.21} $ over a metallicity range 12+log(O/H)=7.5 to 8.44. The observed N/O ratio and scatter are indistinguishable from that observed in low-z galaxies and HII regions over the same metallicity range, showing thus no redshift evolution of N/O for typical galaxies over a significant fraction of cosmic time. We also show that typical $z \sim 3$ galaxies show a similar offset in the BPT diagram as galaxies from the low-z Lyman Continuum Survey (LzLCS), when compared to the average of SDSS galaxies, and show that this offset is not due to enhanced nitrogen abundances. Our results establish a basis for future studies of the evolution of N and O at higher redshifts.

The cycling of carbon between its ionized, atomic, and molecular phases shapes the chemical compositions and physical conditions of the interstellar medium (ISM). However, ground-based studies of the full carbon cycle have been limited by atmospheric absorption. Dome~A, the most promising site for submillimeter astronomy, has long resisted successful submillimeter astronomical observations. Using the 60~cm Antarctic Terahertz Explorer, we present the first successful CO ($4-3$) and [CI] ($^3P_1 - ^3P_0$) mapping observations of two archetypal triggered massive star-formation regions at Dome~A. These data, together with archival [CII], provide the first complete characterization of all three carbon phases in these environments. We find elevated C$^{0}$/CO abundance ratios in high-extinction regions, plausibly driven by deep penetration of intense radiation fields from massive stars into a clumpy ISM. These findings mark a major milestone for submillimeter astronomy at Dome~A and offer valuable insights into the impact of massive star feedback on the surrounding ISM.

Magnetars are the most strongly magnetized compact objects known in the Universe and are regarded as one of the primary engines powering a variety of enigmatic, high-energy transients. However, our understanding of magnetars remains highly limited, constrained by observational sample size and radiative variability. XTE~J1810$-$197, which re-entered a radio-active phase in 2018, is one of only six known radio-pulsating magnetars. Leveraging the distinctive capability for simultaneous dual-frequency observations, we utilized the Shanghai Tianma Radio Telescope (TMRT) to monitor this magnetar continuously at both 2.25 and 8.60~GHz, capturing its entire evolution from radio activation to quenching. This enabled precise characterization of the evolution in its integrated profile, spin frequency, flux density, and spectral index ($\alpha$, defined by $S \propto f^{\alpha}$). The first time derivative of its spin frequency $\dot{\nu}$ passed through four distinct phases -- rapid decrease, violent oscillation, steady decline, and stable recovery -- before returning to its pre-outburst value concomitant with the cessation of radio emission. Remarkably, both the amplitudes and the characteristic time-scales of these $\dot{\nu}$ variations match those observed during the previous outburst that began in 2003, providing the first demonstration that post-outburst rotational evolution and radiative behavior in a magnetar are repeatable. A twisted-magnetosphere model can qualitatively account for this repeatability as well as for the progressive narrowing and abrupt disappearance of the radio pulse radiation, thereby receiving strong observational support.

Small-scale discrepancies in the standard Lambda cold dark matter paradigm have motivated the exploration of alternative dark matter (DM) models, such as warm dark matter (WDM). We investigate the constraining power of galaxy scaling relations on cosmological, astrophysical, and WDM parameters through a joint analysis of hydrodynamic simulations and observational data. Our study is based on the DREAMS project and combines large-volume uniform-box simulations with high-resolution Milky Way zoom-in runs in a $\Lambda$WDM cosmology. To ensure consistency between the different simulation sets, we apply calibrations to account for resolution effects, allowing us to exploit the complementary strengths of the two suites. We compare simulated relations, including stellar size, DM mass and fraction within the stellar half-mass radius, and the total-to-stellar mass ratio, with two complementary galaxy samples: the SPARC catalog of nearby spirals and the LVDB catalog of dwarf galaxies in the Local Volume. Using a bootstrap-based fitting procedure, we show that key cosmological parameters ($\Omega_m$, $\sigma_8$) and supernova feedback strength can be recovered with good accuracy, particularly from the uniform-box simulations. While the WDM particle mass remains unconstrained, the zoom-in simulations reveal subtle WDM-induced trends at low stellar masses in both the DM mass and total-to-stellar mass ratio. We also find that the galaxy stellar mass function exhibits a measurable dependence on the WDM particle mass below log10(M_*/Msun) <~ 8, which appears separable from the impact of feedback, suggesting it as a promising complementary probe. Our results highlight the importance of combining multi-resolution simulations with diverse observational datasets to jointly constrain baryonic processes and DM properties.

Flares are usually observed during the afterglow phase of Gamma-Ray Bursts (GRBs) in soft X-ray, optical and radio bands, but rarely in gamma-ray band. Despite the extraordinary brightness, GECAM-C has accurately measured both the bright prompt emission and flare emission of GRB 221009A without instrumental effects, offering a good opportunity to study the relation between them. In this work, we present a comprehensive analysis of flare emission of GRB 221009A, which is composed of a series of flares. Among them, we identify an exceptionally bright flare with a record-breaking isotropic energy $E_{\rm iso} = 1.82 \times 10^{53}$ erg of GRB flares. It exhibits the highest peak energy ever detected in GRB flares, $E_{\rm peak} \sim 300$ keV, making it a genuine gamma-ray flare. It also shows rapid rise and decay timescales, significantly shorter than those of typical X-ray flares observed in soft X-ray or optical band, but comparable to those observed in prompt emissions. Despite these exceptional properties, the flare shares several common properties with typical GRB flares. We note that this is the first observation of a GRB flare in the keV-MeV band with sufficiently high temporal resolution and high statistics, which bridges the last gap between prompt emission and flare.

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1\sigma$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

Multiwavelength observations, from radio to X-rays, have revealed the presence of multiphase high-velocity gas near the center of the Milky Way likely associated with powerful galactic outflows. This region offers a unique laboratory to study the physics of feedback and the nature of multiphase winds in detail. To this end, we have developed physically motivated semi-analytical models of a multiphase outflow consisting of a hot gas phase ($T \gg 10^6$ K) that embeds colder clouds ($T \sim 5000$ K). Our models include the gravitational potential of the Milky Way; the drag force exerted by the hot phase onto the cold clouds; and the exchange of mass, momentum, and energy between gas phases. Using Bayesian inference, we compared the predictions of our models with observations of a population of HI high-velocity clouds detected up to $\sim$1.5 kpc above the Galactic plane near the Galactic center. We find that a class of supernova-driven winds launched by star formation in the central molecular zone can successfully reproduce the observed velocities, spatial distribution, and masses of the clouds. In our two-phase models, the mass and energy loading factors of both phases are consistent with recent theoretical expectations. The cold clouds are accelerated by the hot wind via ram pressure drag and via accretion of high-velocity material, resulting from the turbulent mixing and subsequent cooling. However, this interaction also leads to gradual cloud disruption, with smaller clouds losing over 70\% of their initial mass by the time they reach $\sim$2 kpc.

Stellar-mass black holes (BHs) can form from the near-complete collapse of massive stars, causing them to abruptly disappear. The star M31-2014-DS1 in the Andromeda galaxy was reported to exhibit such a disappearance between 2014 and 2022, with properties consistent with the failed explosion of a $\approx 12 - 13$ M$_\odot$ yellow supergiant leading to the formation of a $\approx 5$ M$_\odot$ BH. We present mid-infrared (MIR) observations of the remnant obtained with the James Webb Space Telescope (JWST) and X-ray observations from the Chandra X-ray Observatory in 2024. The JWST MIRI/NIRSpec data reveal an extremely red source, showing strong blueshifted absorption from molecular gas (CO, CO$_2$, H$_2$O, SO$_2$) and deep silicate dust features. Modeling the dust continuum confirms continued bolometric fading of the central source to $\log(L/L_\odot)\approx3.88$ ($\approx7-8$% of the progenitor luminosity), surrounded by a dust shell spanning $\approx40-200$ au. Modeling of the molecular gas indicates $\sim 0.1$ M$_\odot$ of gas expanding at $\approx 100$ km s$^{-1}$ near the inner edge of the dust shell. No X-ray source is detected down to a luminosity limit of $L_X\lesssim1.5\times10^{35}$ erg s$^{-1}$. We show that the panchromatic observations are explained by (i) a low-energy ($\approx10^{46}$ erg) ejection of the outer H-rich progenitor envelope and (ii) a fading central BH powered by inefficient ($\sim0.1$% in mass) accretion of loosely bound fallback material. The analysis robustly establishes the bolometric fading of M31-2014-DS1 and provides the first cohesive insights into BH formation via low-energy explosions and long-term fallback.

Stars that exhibit prominent emission lines in their spectra are referred to as emission-line stars, encompassing a wide range of stellar types and indicative of intriguing physical properties. The Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) has released millions of spectra from its Medium-Resolution Survey (MRS). A small fraction of these spectra exhibit emission lines, yet they remain undiscovered and unanalyzed due to being buried in the vast dataset. We have developed a method based on Derivative Spectroscopy (DS), which provides a novel approach for detecting and identifying emission-line stars by extracting signals from complex backgrounds and estimating spectral line profiles. Applying this method to the $\mathrm{H\alpha}$ spectral line profiles from the LAMOST-MRS Data Release 7 (DR7), we compiled a catalog of emission-line stars using the second- and third-derivative spectra for automated peak detection. This approach also facilitates the classification of $\mathrm{H\alpha}$ emission-line morphologies through a simplified scheme. The catalog comprises 56\,649 spectra with relatively prominent $\mathrm{H\alpha}$ emission lines from 45\,206 unique stars, with each emission-line component accompanied by approximate estimates of its wavelength, amplitude, and width. All $\mathrm{H\alpha}$ spectral lines were classified into three morphological classes under a unified classification scheme: single emission peak (83.0 \%), double emission peaks (5.6 \%), and P Cygni-type profiles (11.5 \%), which encompass both P Cygni and inverse P Cygni features. Through cross-referencing with SIMBAD, 39\,497 stars represent new emission-line sources discovered in our research.

We present a rotational velocity (vsini) survey of 32 stellar/substellar objects and giant planets using Keck/KPIC high-resolution spectroscopy, including 6 giant planets (2-7 M$_\mathrm{Jup}$) and 25 substellar/stellar companions (12-88 M$_\mathrm{Jup}$). Adding companions with spin measurements from the literature, we construct a curated spin sample for 43 benchmark stellar/substellar companions and giant planets and 54 free-floating brown dwarfs and planetary mass objects. We compare their spins, parameterized as fractional breakup velocities at 10 Myr, assuming constant angular momentum evolution. We find the first clear evidence that giant planets exhibit distinct spins versus low-mass brown dwarf companions (10 to 40 M$_\mathrm{Jup}$) at 4-4.5 $\sigma$ significance assuming inclinations aligned with their orbits, while under randomly oriented inclinations the significance is at 1.6-2.1 $\sigma$. Our findings hold when considering various assumptions about planets, and the mass ratio below 0.8% gives a clean cut for rotation between giant planets and brown dwarf companions. The higher fractional breakup velocities of planets can be interpreted as less angular momentum loss through circumplanetary disk braking during the planet formation phase. Brown dwarf companions exhibit evidence of slower rotation compared to isolated brown dwarfs, while planets and planetary mass objects show similar spins. Finally, our analysis of specific angular momentum versus age of 221 stellar/substellar objects below 0.1 MSun with spin measurements in the literature indicates that the substellar objects of 5-40 M$_\mathrm{Jup}$ retain much higher angular momenta compared to stellar and substellar objects of 40-100 M$_\mathrm{Jup}$ after 10 Myr, when their initial angular momenta were set.

We present a search for ultralight axion dark matter coupled to electron spins using a levitated ferromagnetic torsional oscillator (FMTO). This platform directly measures axion-induced torques on a macroscopic spin-polarized body, combining large spin density with strong mechanical isolation to probe magnetic fluctuations below 10 Hz while suppressing gradient-field noise. In a first implementation, the experiment yielded 18000 s of analyzable data at room temperature under high vacuum with optical readout and triple-layer magnetic shielding. A likelihood-based statistical framework, incorporating stochastic fluctuations in the axion-field amplitude, was used to evaluate the data. No excess consistent with an axion-induced pseudo-magnetic field was observed near 2e-14 eV. To account for possible shielding-induced signal attenuation, we quantify its effect and report both the uncorrected (g_aee < 1e-7) and attenuation-corrected (g_aee < 6e-5) 90% CL limits on the axion-electron coupling. Looking ahead, improvements guided by both noise-budget analysis and shielding-attenuation considerations, including optimized levitation geometry, cryogenic operation, and superconducting shielding, are expected to boost sensitivity by multiple orders of magnitude.

Abell 3571 is a nearby, X-ray bright galaxy cluster located in the Shapley Supercluster. Although it appears morphologically relaxed in X-ray images, multiwavelength observations reveal subtle indications of residual dynamical activity, likely associated with past merger events. Using wide-field ($1^{\circ} \times 1^{\circ}$) data from the Einstein Probe Follow-up X-ray Telescope (EP-FXT), we extend measurements of the cluster's properties beyond its $R_{500}$ radius. We detect surface-brightness excesses on both the northern and southern sides, consistent with previous XMM-Newton results. The temperature, pressure, and entropy in the northern excess region are lower than the average values, whereas those on the southern side are slightly higher. However, we find no evidence for cold fronts or shocks. These features can be interpreted as sloshing motions triggered by an off-center minor merger. Our findings suggest that, despite its symmetric appearance, A3571 is still recovering from a minor merger and is currently in a post-merger phase. This work also demonstrates the excellent capability of EP-FXT for probing the outskirts of galaxy clusters.

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

Microwave cavities operated in the presence of a background magnetic field provide a promising avenue for detecting high-frequency gravitational waves (HFGWs). We demonstrate for the first time that the distinct antenna patterns of multiple electromagnetic modes within a single cavity enable localization and reconstruction of key properties of an incoming HFGW signal, including its polarization ratio and frequency drift rate. Using a 9-cell cavity commonly employed in particle accelerators as a representative example, we analyze the time-domain response of 18 nearly degenerate modes, which can be sequentially excited by a frequency-drifting signal. The sensitivity is further enhanced by the number of available modes, in close analogy to the scaling achieved by a network of independent detectors, enabling sensitivity to astrophysically plausible binary sources.

The origin of non-thermal power-law distribution functions ubiquitously observed in astrophysical/space (e.g., the solar wind) and laboratory kinetic plasmas, is not well understood. Another puzzling phenomenon is temperature inversion in the solar corona. These two issues are deeply connected. We develop a self-consistent quasilinear theory (QLT) for electromagnetically driven kinetic plasmas, deriving a Fokker-Planck equation for the simultaneous relaxation of multiple species, with (i) a drive diffusion coefficient for the heating of dressed particles directly by the drive and indirectly by waves, and (ii) Balescu-Lenard diffusion and drag coefficients for internal turbulence and Coulomb collisions. Both electron and ion distributions relax towards a universal attractor with a $v^{-5}$ $(E^{-2})$ tail, akin to a $\kappa = 1.5$ distribution, under a super-Debye (but sub-Larmor) drive with a steep power-spectrum. This is an outcome of Debye screening: large-scale fields accelerate the unscreened, fast particles but not the screened, slow ones. The universality may be broken by shallow power-spectra and incomplete relaxation. Collisions cannot decelerate suprathermal particles, rendering a high $v$ tail immune to Maxwellianization. Such a tail may be generated in the solar corona by chromospheric convection despite collisional losses. The suprathermal particles escape sun's gravity (velocity filtration), inverting the temperature profile and raising it to $10^6$ K. A proper analysis of velocity filtration with a $\kappa \approx 1.5-2$ distribution inspired by QLT provides a reasonable fit to the spectroscopic data of heavy ions and explains the abrupt temperature rise, a consequence of the divergence of pressure in the $\kappa \to 1.5$ limit.

Massive short-period binaries involving hot subdwarf stars (sdO/Bs) are rare but important to constraining pathways for binary star evolution. Moreover, some of the most promising candidate progenitor systems leading to Type Ia supernovae (SNe Ia) involve sdO/Bs. LAMOST J065816.72+094343.1 has been identified as such a candidate. To explore the nature and evolutionary future of LAMOST J065816.72+094343.1, we complemented archival spectroscopic data with additional time series spectra and high-resolution spectroscopy of the object. After combining these with photometric data, we determined the orbital parameters of the system and the mass of the companion. We solved the orbit of the system by analyzing 68 low- and medium-resolution spectra using state-of-the-art mixed local thermodynamic equilibrium (LTE) and non-LTE model atmospheres. Additionally, we gathered nine high-resolution spectra to determine atmospheric parameters and the projected rotational velocity of the sdOB. The inclination angle of the system was constrained assuming tidal synchronization of the sdOB, which was verified via analysis of the ellipsoidal variations in the TESS light curve. We determine LAMOSTJ065816.72+094343.1 to be a binary consisting of a massive $0.82 \pm 0.17 \mathrm{M}_{\odot}$ sdOB component with a $1.30^{+0.31}_{-0.26} \mathrm{M}_{\odot}$ unseen companion. Due to the companion's mass being very close to the Chandrasekhar mass limit and high for a white dwarf, it is unclear whether it is a white dwarf or a neutron star. We find the system to be in a close orbit, with a period of $P=0.31955193 \mathrm{d}$ and an inclination angle of $i = 49.6^{+5.2}_{-4.2} \mathrm{deg}$. While the exact nature of the companion remains unknown, we determine the system to either lead to a SN Ia or an intermediate mass binary pulsar, potentially after a phase as an intermediate-mass X-ray binary.