Search

The recently discovered Pa 30 nebula, the putative type Iax supernova remnant associated with the historical supernova of 1181 AD, shows puzzling characteristics that make it unique among known supernova remnants. In particular, Pa 30 exhibits a complex morphology, with a unique radial and filamentary structure, and it hosts a hot stellar remnant at its center, which displays oxygen-dominated, ultra-fast winds. Because of the surviving stellar remnant and the lack of hydrogen and helium in its filaments, it has been suggested that Pa 30 is the product of a failed thermonuclear explosion in a near- or super-Chandrasekhar white dwarf, which created a sub-luminous transient, a rare sub-type of the Ia class of supernovae called type Iax. We here present a detailed study of the 3D structure and velocities of a full radial section of the remnant. The Integral Field Unit (IFU) observations, obtained with the new red channel of the Keck Cosmic Web Imager spectrograph, reveal that the ejecta are consistent with being ballistic, with velocities close to the free-expansion velocity. Additionally, we detect a large cavity inside the supernova remnant and a sharp inner edge to the filamentary structure, which coincides with the outer edge of a bright ring detected in infrared images. Finally, we detect a strong asymmetry in the amount of ejecta along the line of sight, which might hint to an asymmetric explosion. Our analysis provides strong confirmation that the explosion originated from SN 1181., Comment: 13 pages, 8 figures. Accepted for publication in ApJL

Over the past several decades, time-series photometry of CSPNe has yielded significant results including, but not limited to, discoveries of nearly 100 binary systems, insights into pulsations and winds in young white dwarfs, and studies of stars undergoing very late thermal pulses. We have undertaken a systematic study of optical photometric variability of cataloged CSPNe, using the epochal photometric data from the Zwicky Transient Facility (ZTF). By applying appropriate variability metrics, we arrive at a list of 94 significantly variable CSPNe. Based on the timescales of the light-curve activity, we classify the variables broadly into short- and long-timescale variables. In this first paper in this series, we focus on the former, which is the majority class comprising 83 objects. We infer periods for six sources for the first time, and recover several known periodic variables. Among the aperiodic sources, most exhibit a jitter around a median flux with a stable amplitude, and a few show outbursts. We draw attention to WeSb 1, which shows a different kind of variability: prominent deep and aperiodic dips, resembling transits from a dust/debris disk. We find strong evidence for a binary nature of WeSb 1 (possibly an A- to G-type companion). The compactness of the emission lines and inferred high electron densities make WeSb 1 a candidate for either an EGB 6-type planetary nucleus, or a symbiotic system inside an evolved planetary nebula, both of which are rare objects. To demonstrate further promise with ZTF, we report three additional newly identified periodic sources that do not appear in the list of highly variable sources. Finally, we also introduce a two-dimensional metric space defined by the von Neumann statistics and Pearson Skew and demonstrate its effectiveness in identifying unique variables of astrophysical interest, like WeSb 1., Comment: 19 pages + 8 pages appendix, 5 tables, 17 figures; Submitted to PASP; Comments are welcome!

This work focuses on the numerical solution of the dynamics of a poroelastic material in the frequency domain. We provide a detailed stability analysis based on the application of the Fredholm alternative in the continuous case, considering a total pressure formulation of the Biot's equations. In the discrete setting, we propose a stabilized equal order finite element method complemented by an additional pressure stabilization to enhance the robustness of the numerical scheme with respect to the fluid permeability. Utilizing the Fredholm alternative, we extend the well-posedness results to the discrete setting, obtaining theoretical optimal convergence for the case of linear finite elements. We present different numerical experiments to validate the proposed method. First, we consider model problems with known analytic solutions in two and three dimensions. As next, we show that the method is robust for a wide range of permeabilities, including the case of discontinuous coefficients. Lastly, we show the application for the simulation of brain elastography on a realistic brain geometry obtained from medical imaging.

We present volume-limited samples of cataclysmic variables (CVs) and AM CVn binaries jointly selected from SRG/eROSITA eRASS1 and \textit{Gaia} DR3 using an X-ray + optical color-color diagram (the ``X-ray Main Sequence"). This tool identifies all CV subtypes, including magnetic and low-accretion rate systems, in contrast to most previous surveys. We find 23 CVs, 3 of which are AM CVns, out to 150 pc in the Western Galactic Hemisphere. Our 150 pc sample is spectroscopically verified and complete down to $L_X = 1.3\times 10^{29} \;\textrm{erg s}^{-1}$ in the 0.2--2.3 keV band, and we also present CV candidates out to 300 pc and 1000 pc. We discovered two previously unknown systems in our 150 pc sample: the third nearest AM CVn and a magnetic period bouncer. We find the mean $L_X$ of CVs to be $\langle L_X \rangle \approx 4.6\times 10^{30} \;\textrm{erg s}^{-1}$, in contrast to previous surveys which yielded $\langle L_X \rangle \sim 10^{31}-10^{32} \;\textrm{erg s}^{-1}$. We construct X-ray luminosity functions that, for the first time, flatten out at $L_X\sim 10^{30} \; \textrm{erg s}^{-1}$. We find average number, mass, and luminosity densities of $\rho_\textrm{N, CV} = (3.7 \pm 0.7) \times 10^{-6} \textrm{pc}^{-3}$, $\rho_M = (5.0 \pm 1.0) \times 10^{-5} M_\odot^{-1}$, and $\rho_{L_X} = (2.3 \pm 0.4) \times 10^{26} \textrm{erg s}^{-1}M_\odot^{-1}$, respectively, in the solar neighborhood. Our uniform selection method also allows us to place meaningful estimates on the space density of AM CVns, $\rho_\textrm{N, AM CVn} = (5.5 \pm 3.7) \times 10^{-7} \textrm{pc}^{-3}$. Magnetic CVs and period bouncers make up $35\%$ and $25\%$ of our sample, respectively. This work, through a novel discovery technique, shows that the observed number densities of CVs and AM CVns, as well as the fraction of period bouncers, are still in tension with population synthesis estimates., Comment: Submitted to PASP, comments welcome

It has been suggested that giant planet occurrence peaks for stars with $M_{\ast}~\approx~3~M_{\odot}$ at a value a factor of four higher than observed for solar-mass stars. This population of giant planets predicted to frequently orbit main sequence B stars at $a~\approx~10$ AU is difficult to characterize observationally while fusion persists in their host stars. Fortunately, stars with $M_{\ast}~\gtrsim~3~M_{\odot}$ sustain fusion for only a few hundred million years. By the time those stars become massive, young white dwarfs, any giant planets present would still be luminous as a consequence of their recent formation. From an initial sample of 3268 Gaia-identified massive, young white dwarfs, we use homogeneous Spitzer Infrared Array Camera (IRAC) photometry to search for evidence of unresolved giant planets. For 30 systems, these IRAC data provide sensitivity to objects with $M~\lesssim~10~M_{\text{Jup}}$, and we identify one candidate with $M~\approx~4~M_{\text{Jup}}$ orbiting the white dwarf GALEX J071816.4+373139. Correcting for the possibility that some of the white dwarfs in our sample result from mergers, we find a giant planet occurrence $\eta_{\text{GP}}~=~0.11_{-0.07}^{+0.13}$ for stars with initial masses $M_{\ast}~\gtrsim~3~M_{\odot}$. Our occurrence inference is consistent with both the Doppler-inferred occurrence of giant planets orbiting $M_{\ast}~\approx~2~M_{\odot}$ giant stars and the theoretically predicted factor-of-four enhancement in the occurrence of giant planets orbiting $M_{\ast}~\approx~3~M_{\odot}$ stars relative to solar-mass stars. Future James Webb Space Telescope Near Infrared Camera observations of our sample would provide sensitivity to Saturn-mass planets and thereby a definitive estimate of the occurrence of giant planets orbiting stars with $M_{\ast}~\gtrsim~3~M_{\odot}$., Comment: 10 pages, 6 figures

Cataclysmic variables (CVs) are compact binary systems in which a white dwarf accretes matter from a Roche-lobe-filling companion star. In this study, we searched for new CVs in the Milky Way in the Chandra Source Catalog v2.0, cross-matched with Gaia Data Release 3 (DR3). We identified new CV candidates by combining X-ray and optical data in a color-color diagram called the ``X-ray Main Sequence". We used two different cuts in this diagram to compile pure and optically variable samples of CV candidates. We undertook optical spectroscopic follow-up observations with the Keck and Palomar Observatories to confirm the nature of these sources. We assembled a sample of 25,887 Galactic X-ray sources and found 14 new CV candidates. Seven objects show X-ray and/or optical variability. All sources show X-ray luminosity in the $\rm 10^{29}-10^{32}$ $\rm erg\ s^{-1}$ range, and their X-ray spectra can be approximated by a power-law model with photon indices in the $\rm \Gamma \sim 1-3$ range or an optically thin thermal emission model in the $\rm kT \sim 1-70$ keV range. We spectroscopically confirmed four CVs, discovering two new polars, one low accretion rate polar and a WZ~Sge-like low accretion rate CV. X-ray and optical properties of the other 9 objects suggest that they are also CVs (likely magnetic or dwarf novae), and one other object could be an eclipsing binary, but revealing their true nature requires further observations. These results show that a joint X-ray and optical analysis can be a powerful tool for finding new CVs in large X-ray and optical catalogs. X-ray observations such as those by Chandra are particularly efficient at discovering magnetic and low accretion rate CVs, which could be missed by purely optical surveys., Comment: 20 pages, 15 figures and 8 tables. Accepted for publication in Astronomy & Astrophysics

Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We use an existing \textsc{Arepo} hydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the 1 dimensional stellar evolution code MESA. We show that the supernova shock is not enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the $10^{5-6}$ year lifetimes inferred for observed hypervelocity runaways. Despite experiencing a shock from a supernova less than $\approx 0.02\,R_\odot$ away, our models do not experience significant interior heating, and all contract back to radii around $0.01\,R_\odot$ within about $10^4$\,years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution., Comment: Submitted to A\&A. 15 pages, 17 figures

We discuss how slip conditions for the Stokes equation can be handled using Nitsche method, for a stabilized finite element discretization. Emphasis is made on the interplay between stabilization and Nitsche terms. Well-posedness of the discrete problem and optimal convergence rates, in natural norm for the velocity and the pressure, are established, and illustrated with various numerical experiments. The proposed method fits naturally in the context of a finite element implementation while being accurate, and allows an increased flexibility in the choice of the finite element pairs., Comment: 16 pages, 10 figures

Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around $\approx 40-70\%$ of the Galactic CVs are post-period minimum systems referred to as "period bouncers", only a few dozen confirmed systems are known. We report the study and characterisation of a new eclipsing CV, SRGeJ041130.3+685350 (SRGeJ0411), discovered from a joint SRG/eROSITA and ZTF program. The optical spectrum of SRGeJ0411 shows prominent hydrogen and helium emission lines, typical for CVs. We obtained optical high-speed photometry to confirm the eclipse of SRGeJ0411 and determine the orbital period to be $P_\textrm{orb} \approx 97.530$ minutes. The spectral energy distribution suggests that the donor has an effective temperature of $\lesssim 1,800$ K. We constrain the donor mass with the period--density relationship for Roche-lobe-filling stars and find that $M_\textrm{donor} \lesssim 0.04\ M_\odot$. The binary parameters are consistent with evolutionary models for post-period minimum CVs, suggesting that SRGeJ0411 is a new period bouncer. The optical emission lines of SRGeJ0411 are single-peaked despite the system being eclipsing, which is typically only seen due to stream-fed accretion in polars. X-ray spectroscopy hints that the white dwarf in SRGeJ0411 could be magnetic, but verifying the magnetic nature of SRGeJ0411 requires further investigation. The lack of optical outbursts has made SRGeJ0411 elusive in previous surveys, and joint X-ray and optical surveys highlight the potential for discovering similar systems in the near future., Comment: 17 pages, 16 figures, accepted to MNRAS

The integration of experimental data into mathematical and computational models is crucial for enhancing their predictive power in real-world scenarios. However, the performance of data assimilation algorithms can be significantly degraded when measurements are corrupted by biased noise, altering the signal magnitude, or when the system dynamics lack smoothness, such as in the presence of fast oscillations or discontinuities. This paper focuses on variational state estimation using the so-called Parameterized Background Data Weak method, which relies on a parameterized background by a set of constraints, enabling state estimation by solving a minimization problem on a reduced-order background model, subject to constraints imposed by the input measurements. To address biased noise in observations, a modified formulation is proposed, incorporating a correction mechanism to handle rapid oscillations by treating them as slow-decaying modes based on a two-scale splitting of the classical reconstruction algorithm. The effectiveness of the proposed algorithms is demonstrated through various examples, including discontinuous signals and simulated Doppler ultrasound data.

We develop two new highly efficient estimators to measure the polarization (Stokes parameters) in experiments that constrain the position angle of individual photons such as scattering and gas-pixel-detector polarimeters, and analyse in detail a previously proposed estimator. All three of these estimators are at least fifty percent more efficient on typical datasets than the standard estimator used in the field. We present analytic estimates of the variance of these estimators and numerical experiments to verify these estimates. Two of the three estimators can be calculated quickly and directly through summations over the measurements of individual photons., Comment: 4 pages, accepted by the Open Journal of Astrophysics

In an accreting X-ray pulsar, a neutron star accretes matter from a stellar companion through an accretion disk. The high magnetic field of the rotating neutron star disrupts the inner edge of the disc, funneling the gas to flow onto the magnetic poles on its surface. Hercules X-1 is in many ways the prototypical X-ray pulsar; it shows persistent X-ray emission and it resides with its companion HZ Her, a two-solar-mass star, at about 7~kpc from Earth. Its emission varies on three distinct timescales: the neutron star rotates every 1.2~seconds, it is eclipsed by its companion each 1.7~days, and the system exhibits a superorbital period of 35~days which has remained remarkably stable since its discovery. Several lines of evidence point to the source of this variation as the precession of the accretion disc, the precession of the neutron star or both. Despite the many hints over the past fifty years, the precession of the neutron star itself has yet not been confirmed or refuted. We here present X-ray polarization measurements with the Imaging X-ray Polarimetry Explorer (IXPE) which probe the spin geometry of the neutron star. These observations provide direct evidence that the 35-day-period is set by the free precession of the neutron star crust, which has the important implication that its crust is somewhat asymmetric fractionally by a few parts per ten million. Furthermore, we find indications that the basic spin geometry of the neutron star is altered by torques on timescale of a few hundred days., Comment: 37 pages

We report on IXPE, NICER and XMM-Newton observations of the magnetar 1E 2259+586. We find that the source is significantly polarized at about or above 20% for all phases except for the secondary peak where it is more weakly polarized. The polarization degree is strongest during the primary minimum which is also the phase where an absorption feature has been identified previously (Pizzocaro et al. 2019). The polarization angle of the photons are consistent with a rotating vector model with a mode switch between the primary minimum and the rest of the rotation of the neutron star. We propose a scenario in which the emission at the source is weakly polarized (as in a condensed surface) and, as the radiation passes through a plasma arch, resonant cyclotron scattering off of protons produces the observed polarized radiation. This confirms the magnetar nature of the source with a surface field greater than about 1015 G, Comment: 13 pages, 13 figures, 5 tables

Despite their shared origin, members of globular clusters display star-to-star variations in composition. The observed pattern of element abundances is unique to these stellar environments, and cannot be fully explained by any proposed mechanism. It remains unclear whether stars form with chemical heterogeneity, or inherit it from interactions with other members. These scenarios may be differentiated by the dependence of chemical spread on stellar mass; however, obtaining a sufficiently large mass baseline requires abundance measurements on the lower main sequence that is too faint for spectroscopy even in the nearest globular clusters. We developed a stellar modelling method to obtain precise chemical abundances for stars near the end of the main sequence from multiband photometry, and applied it to the globular cluster 47 Tucanae. The computational efficiency is attained by matching chemical elements to the model components that are most sensitive to their abundance. We determined [O/Fe] for ~5000 members below the main sequence knee at the level of accuracy, comparable to the spectroscopic measurements of evolved members in literature. The inferred distribution disfavors stellar interactions as the origin of chemical spread; however, an accurate theory of accretion is required to draw a more definitive conclusion. We anticipate that future observations of 47 Tucanae with JWST will extend the mass baseline of our analysis into the substellar regime. Therefore, we present predicted color-magnitude diagrams and mass-magnitude relations for the brown dwarf members of 47 Tucanae.

We searched the Gaia DR3 database for ultramassive white dwarfs with kinematics consistent with having escaped the nearby Hyades open cluster, identifying three such candidates. Two of these candidates have masses estimated from Gaia photometry of approximately 1.1 solar masses; their status as products of single stellar evolution that have escaped the cluster was deemed too questionable for immediate follow-up analysis. The remaining candidate has an expected mass >1.3 solar masses, significantly reducing the probability of it being an interloper. Analysis of follow-up Gemini GMOS spectroscopy for this source reveals a non-magnetized hydrogen atmosphere white dwarf with a mass and age consistent with having formed from a single star. Assuming a single-stellar evolution formation channel, we estimate a 97.8% chance that the candidate is a true escapee from the Hyades. With a determined mass of 1.317 solar masses, this is potentially the most massive known single-evolution white dwarf and is by far the most massive with a strong association with an open cluster., Comment: 13 pages, 4 figures, 1 table. Accepted for publication in the Astrophysical Journal Letters

This paper presents a numerical method for the simulation of elastic solid materials coupled to fluid inclusions. The application is motivated by the modeling of vascularized tissues and by problems in medical imaging which target the estimation of effective (i.e., macroscale) material properties, taking into account the influence of microscale dynamics, such as fluid flow in the microvasculature. The method is based on the recently proposed Reduced Lagrange Multipliers framework. In particular, the interface between solid and fluid domains is not resolved within the computational mesh for the elastic material but discretized independently, imposing the coupling condition via non-matching Lagrange multipliers. Exploiting the multiscale properties of the problem, the resulting Lagrange multipliers space is reduced to a lower-dimensional characteristic set. We present the details of the stability analysis of the resulting method considering a non-standard boundary condition that enforces a local deformation on the solid-fluid boundary. The method is validated with several numerical examples.

Background Weight loss trajectories after bariatric surgery vary widely between individuals, and predicting weight loss before the operation remains challenging. We aimed to develop a model using machine learning to provide individual preoperative prediction of 5-year weight loss trajectories after surgery. Methods In this multinational retrospective observational study we enrolled adult participants (aged $\ge$18 years) from ten prospective cohorts (including ABOS [NCT01129297], BAREVAL [NCT02310178], the Swedish Obese Subjects study, and a large cohort from the Dutch Obesity Clinic [Nederlandse Obesitas Kliniek]) and two randomised trials (SleevePass [NCT00793143] and SM-BOSS [NCT00356213]) in Europe, the Americas, and Asia, with a 5 year followup after Roux-en-Y gastric bypass, sleeve gastrectomy, or gastric band. Patients with a previous history of bariatric surgery or large delays between scheduled and actual visits were excluded. The training cohort comprised patients from two centres in France (ABOS and BAREVAL). The primary outcome was BMI at 5 years. A model was developed using least absolute shrinkage and selection operator to select variables and the classification and regression trees algorithm to build interpretable regression trees. The performances of the model were assessed through the median absolute deviation (MAD) and root mean squared error (RMSE) of BMI. Findings10 231 patients from 12 centres in ten countries were included in the analysis, corresponding to 30 602 patient-years. Among participants in all 12 cohorts, 7701 (75$\bullet$3%) were female, 2530 (24$\bullet$7%) were male. Among 434 baseline attributes available in the training cohort, seven variables were selected: height, weight, intervention type, age, diabetes status, diabetes duration, and smoking status. At 5 years, across external testing cohorts the overall mean MAD BMI was 2$\bullet$8 kg/m${}^2$ (95% CI 2$\bullet$6-3$\bullet$0) and mean RMSE BMI was 4$\bullet$7 kg/m${}^2$ (4$\bullet$4-5$\bullet$0), and the mean difference between predicted and observed BMI was-0$\bullet$3 kg/m${}^2$ (SD 4$\bullet$7). This model is incorporated in an easy to use and interpretable web-based prediction tool to help inform clinical decision before surgery. InterpretationWe developed a machine learning-based model, which is internationally validated, for predicting individual 5-year weight loss trajectories after three common bariatric interventions., Comment: The Lancet Digital Health, 2023

White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterised by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink toward the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium. Several mechanisms compete with gravitational settling to change a white dwarf's surface composition as it cools, and the fraction of white dwarfs with helium atmospheres is known to increase by a factor ~2.5 below a temperature of about 30,000 K; therefore, some white dwarfs that appear to have hydrogen-dominated atmospheres above 30,000 K are bound to transition to be helium-dominated as they cool below it. Here we report observations of ZTF J203349.8+322901.1, a transitioning white dwarf with two faces: one side of its atmosphere is dominated by hydrogen and the other one by helium. This peculiar nature is likely caused by the presence of a small magnetic field, which creates an inhomogeneity in temperature, pressure or mixing strength over the surface. ZTF J203349.8+322901.1 might be the most extreme member of a class of magnetic, transitioning white dwarfs -- together with GD 323, a white dwarf that shows similar but much more subtle variations. This new class could help shed light on the physical mechanisms behind white dwarf spectral evolution., Comment: 45 pages, 11 figures

AM CVn systems are ultra-compact binaries where a white dwarf accretes from a helium-rich degenerate or semi-degenerate donor. Some AM CVn systems will be among the loudest sources of gravitational waves for the upcoming Laser Interferometer Space Antenna (LISA), yet the formation channel of AM CVns remains uncertain. We report the study and characterisation of a new eclipsing AM CVn, SRGeJ045359.9+622444 (hereafter SRGeJ0453), discovered from a joint SRG/eROSITA and ZTF program to identify cataclysmic variables (CVs). We obtained optical photometry to confirm the eclipse of SRGeJ0453 and determine the orbital period to be $P_\textrm{orb} = 55.0802 \pm 0.0003$ min. We constrain the binary parameters by modeling the high-speed photometry and radial velocity curves and find $M_\textrm{donor} = 0.044 \pm0.024 M_{\odot}$ and $R_\textrm{donor}=0.078 \pm 0.012 R_{\odot}$. The X-ray spectrum is approximated by a power-law model with an unusually flat photon index of $\Gamma\sim 1$ previously seen in magnetic CVs with SRG/eROSITA, but verifying the magnetic nature of SRGeJ0453 requires further investigation. Optical spectroscopy suggests that the donor star of SRGeJ0453 could have initially been a He star or a He white dwarf. SRGeJ0453 is the ninth eclipsing AM CVn system published to date, and its lack of optical outbursts have made it elusive in previous surveys. The discovery of SRGeJ0453 using joint X-ray and optical surveys highlights the potential for discovering similar systems in the near future., Comment: Submitted to ApJ

We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain., Comment: 26 pages, 17 figures. Accepted to OJA

We report on Imaging X-ray polarimetry explorer (IXPE) observations of the Be-transient X-ray pulsar LS V +44 17/RX J0440.9+4431 made at two luminosity levels during the giant outburst in January--February 2023. Considering the observed spectral variability and changes in the pulse profiles, the source was likely caught in supercritical and subcritical states with significantly different emission-region geometry, associated with the presence of accretion columns and hot spots, respectively. We focus here on the pulse-phase-resolved polarimetric analysis and find that the observed dependencies of the polarization degree and polarization angle (PA) on the pulse phase are indeed drastically different for the two observations. The observed differences, if interpreted within the framework of the rotating vector model (RVM), imply dramatic variations in the spin axis inclination, the position angle, and the magnetic colatitude by tens of degrees within the space of just a few days. We suggest that the apparent changes in the observed PA phase dependence are predominantly related to the presence of an unpulsed polarized component in addition to the polarized radiation associated with the pulsar itself. We then show that the observed PA phase dependence in both observations can be explained with a single set of RVM parameters defining the pulsar's geometry. We also suggest that the additional polarized component is likely produced by scattering of the pulsar radiation in the equatorial disk wind., Comment: 10 pages, 7 figures, accepted in A&A

We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post-common envelope carbon-oxygen (CO) white dwarf, and a warm donor ($T_{\rm eff,\,donor}= 16,400\pm1000\,\rm K$). The donor probably formed during a common envelope phase between the CO white dwarf and an evolving giant which left behind a helium star or helium white dwarf in a close orbit with the CO white dwarf. We measure gravitational wave-driven orbital inspiral with $\sim 35\sigma$ significance, which yields a joint constraint on the component masses and mass transfer rate. While the accretion disk in the system is dominated by ionized helium emission, the donor exhibits a mixture of hydrogen and helium absorption lines. Phase-resolved spectroscopy yields a donor radial-velocity semi-amplitude of $771\pm27\,\rm km\, s^{-1}$, and high-speed photometry reveals that the system is eclipsing. We detect a {\it Chandra} X-ray counterpart with $L_{X}\sim 3\times 10^{31}\,\rm erg\,s^{-1}$. Depending on the mass-transfer rate, the system will likely evolve into either a stably mass-transferring helium CV, merge to become an R Crb star, or explode as a Type Ia supernova in the next million years. We predict that the Laser Space Interferometer Antenna (LISA) will detect the source with a signal-to-noise ratio of $24\pm6$ after 4 years of observations. The system is the first \emph{LISA}-loud mass-transferring binary with an intrinsically luminous donor, a class of sources that provide the opportunity to leverage the synergy between optical and infrared time domain surveys, X-ray facilities, and gravitational-wave observatories to probe general relativity, accretion physics, and binary evolution., Comment: 13 pages, 7 figures, 2 tables, submitted to ApJL

Recently, the inner main belt asteroid (152830) Dinkinesh was identified as an additional fly-by target for the Lucy mission. The heliocentric orbit and approximate absolute magnitude of Dinkinesh are known, but little additional information was available prior to its selection as a target. In particular, the lack of color spectrophotometry or spectra made it impossible to assign a spectral type to Dinkinesh from which its albedo could be estimated. We set out to remedy this knowledge gap by obtaining visible wavelength spectra with the Keck telescope on 2022 November 23 and with Gemini-South on 2022 December 27. The spectra measured with the Keck I/Low Resolution Imaging Spectrometer (LRIS) and the Gemini South/Gemini Multi-Object Spectrograph South (GMOS-S) are most similar to the average spectrum of S- and Sq-type asteroids. The most diagnostic feature is the $\approx$15$\pm$1$\%$ silicate absorption feature at $\approx$0.9-1.0~micron. Small S- and Sq-type asteroids have moderately high albedos ranging from 0.17-0.35. Using this albedo range for Dinkinesh in combination with measured absolute magnitude, it is possible to derive an effective diameter and surface brightness for this body. The albedo, size and surface brightness are important inputs required for planning a successful encounter by the Lucy spacecraft., Comment: 7 pages, 1 figure. Accepted for publication in Icarus

Magnetars are the most strongly magnetized neutron stars, and one of the most promising targets for X-ray polarimetric measurements. We present here the first Imaging X-ray Polarimetry Explorer (IXPE) observation of the magnetar 1RXS J170849.0-400910, jointly analysed with a new Swift observation and archival NICER data. The total (energy and phase integrated) emission in the 2-8 keV energy range is linerarly polarized, at a ~35% level. The phase-averaged polarization signal shows a marked increase with energy, ranging from ~20% at 2-3 keV up to ~80% at 6-8 keV, while the polarization angle remain constant. This indicates that radiation is mostly polarized in a single direction. The spectrum is well reproduced by a combination of either two thermal (blackbody) components or a blackbody and a power law. Both the polarization degree and angle also show a variation with the spin phase, and the former is almost anti-correlated with the source counts in the 2-8 keV and 2-4 keV bands. We discuss the possible implications and interpretations, based on a joint analysis of the spectral, polarization and pulsation properties of the source. A scenario in which the surface temperature is not homogeneous, with a hotter cap covered by a gaseous atmosphere and a warmer region in a condensed state, provides a satisfactory description of both the phase- and energy-dependent spectro-polarimetric data. The (comparatively) small size of the two emitting regions, required to explain the observed pulsations, does not allow to reach a robust conclusion about the presence of vacuum birefringence effects., Comment: 14 pages, 5 Figures, accepted for publication on ApJL

In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute orbital period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater due to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modeling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs., Comment: 48 Pages, 12 figures, 2 tables, Published online by Nature on Oct 5, 2022

Cen X-3 is the first X-ray pulsar discovered 50 years ago. Radiation from such objects is expected to be highly polarized due to birefringence of plasma and vacuum associated with propagation of photons in presence of the strong magnetic field. Here we present results of the observations of Cen X-3 performed with the Imaging X-ray Polarimetry Explorer. The source exhibited significant flux variability and was observed in two states different by a factor of ~20 in flux. In the low-luminosity state no significant polarization was found either in pulse phase-averaged (with the 3$\sigma$ upper limit of 12%) or phase-resolved data (the 3$\sigma$ upper limits are 20-30%). In the bright state the polarization degree of 5.8$\pm$0.3% and polarization angle of $49.6\deg\pm1.5\deg$ with significance of about 20$\sigma$ was measured from the spectro-polarimetric analysis of the phase-averaged data. The phase-resolved analysis showed a significant anti-correlation between the flux and the polarization degree as well as strong variations of the polarization angle. The fit with the rotating vector model indicates a position angle of the pulsar spin axis of about 49$\deg$ and a magnetic obliquity of 17$\deg$. The detected relatively low polarization can be explained if the upper layers of the neutron star surface are overheated by the accreted matter and the conversion of the polarization modes occurs within the transition region between the upper hot layer and a cooler underlying atmosphere. A fraction of polarization signal can also be produced by reflection of radiation from the neutron star surface and the accretion curtain., Comment: Submitted to ApJL. 18 pages, 12 figures

Numerical simulations of pulsatile blood flow in an aortic coarctation require the use of turbulence modeling. This paper considers three models from the class of large eddy simulation (LES) models (Smagorinsky, Vreman, $\boldsymbol{\sigma}$-model) and one model from the class of variational multiscale models (residual-based) within a finite element framework. The influence of these models on the estimation of clinically relevant biomarkers used to assess the degree of severity of the pathological condition (pressure difference, secondary flow degree, normalized flow displacement, wall shear stress) is investigated in detail. The simulations show that most methods are consistent in terms of severity indicators such as pressure difference and stenotic velocity. The numerical results indicate that second order velocity elements outperform first order elements in terms of accuracy. Moreover, using second order velocity finite elements, different turbulence models might lead to considerably different results concerning other clinically relevant quantities such as wall shear stresses. These differences may be attributed to differences in numerical dissipation introduced by the turbulence models., Comment: 30 pages, 22 figures. Submitted to International Journal for Numerical Methods in Biomedical Engineering