2026-05-15 | | Total: 17
Solar pores possess irregular and evolving boundaries that are often far from the ideal circular flux tubes assumed in many magnetohydrodynamic (MHD) oscillation models. To tackle this irregularity in a consistent way for wave analysis, we introduce a framework that employs the convex hull of the pore boundary - derived from intensity minimum error thresholding - as the domain to perform further analysis. Using the modal assurance criterion, we find the range of pore shapes for which this approximation is valid. We demonstrate the usefulness of this framework by applying it to multi-height, high-cadence observations (4170 angstrom continuum, G-band, Na~\textsc{i}, and Ca~\textsc{ii}~K) of a solar pore, and apply Proper Orthogonal Decomposition of the convex hull to determine wave modes. The fundamental sausage ($m=1$) and kink ($m=2$) mode is found to remain reliable, while higher-order fluting modes ($m\ge3$) are strongly degraded by small-scale boundary irregularity. As expected, sausage-like modes dominate the variance at all heights and exhibit a systematic upward shift in frequency, consistent with freely propagating compressive waves. In contrast, the kink-like motions appear weak, confined to a persistent low-frequency peak, and most plausibly interpreted as a forced response to granular buffeting rather than a propagating mode. Together, these results establish a practical methodology for boundary-mode analysis in real, highly structured pores and provide new constraints on the nature and height evolution of MHD waves in the lower solar atmosphere.
During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions indicating possible energization by magnetic reconnection up to 10s -- 100s keV nucleon$^{-1}$. Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff $E_{\text{max}}$ consistent with in-situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as $E_{\text{max}} \propto (Q/M)^α$ with $α\sim 0.8-1.1$, approximately consistent with PSP measurements in the broader range $α\sim 0.6-1.7$. In the limit where ions are injected at the same energy per nucleon, $α$ can be as low as $\sim 0.3$. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.
Brown dwarfs, being transitional objects between giant planets and low-mass stars, possess dense, cool interiors that provide optimal conditions to explore non-standard physics. Capture and accumulation of dark-matter particles can alter the thermal, structural and dynamic of these substellar objects. We aim to apply a self-consistent two-fluid framework to model the internal structure of self-gravitating brown dwarfs and to quantify how the presence of a dark-matter component modifies their mass--radius relations and dynamical properties. The brown dwarf is modeled as a composite system of a baryonic fluid, described by a polytropic equation of state, and an independent dark-matter fluid. Both components are coupled through their shared gravitational potential in hydrostatic equilibrium. We solve numerically the coupled Lane-Emden equations for a range of dark-matter mass fractions. We find that dark matter accumulating in the core reshapes the baryonic density profile, modifying both the radius and the second-order Love number. Radius and dynamical anomalies in brown dwarfs can serve as diagnostic tools to constrain dark-matter properties. Future high-precision astrometric missions could identify these structural signatures, establishing brown dwarfs as possible detectors of dark matter in the Galaxy.
The NASA TESS mission is conducting a nearly full-sky survey, enabling the photometric characterisation of millions of stars. The forthcoming ESA PLATO mission will provide long-duration, high-precision photometry of tens of thousands of bright stars to be characterised through asteroseismology. The TESS Luminaries Sample is a catalogue of 196 bright naked-eye (V < 6) main-sequence (MS) and sub-giant (SG) stars exhibiting solar-like oscillations. Among them, the subset located within the PLATO Long-duration Observation Phase (LOP) fields constitutes an exceptional set of targets that will be observable by PLATO from the earliest phases of the mission, making them ideal calibrators during commissioning and the first months of science operations. This paper aims to provide an in-depth asteroseismic characterisation of 32 Luminaries stars that fall within the PLATO LOP fields of view. Individual mode parameters were extracted for the first time for 26 of them. We used three independent seismic pipelines, one of which is similar to the algorithms used in the official PLATO pipeline. Statistical criterion were applied to identify the optimal combination of data calibration, observing cadence, and fitting pipeline for each star. For all stars, we derived large and small separations, the asymptotic phase term, radial mode amplitudes, and mean linewidths per order. Comparisons reveal consistent trends in the seismic parameters, confirming the robustness of our analysis. In sub-giants, mixed-mode identification differs in the three pipelines, revealing extraction inconsistencies requiring longer datasets to improve our mode identifications. The Luminaries stars located in the PLATO LOP fields constitute a unique sample that will play a crucial role in validating, calibrating, and optimising PLATO's seismic performance.
The hitherto neglected close binary V486 Car is studied with the aid of newly applied satellite photometry (HIPPARCOS and TESS), high dispersion spectrometry (HERCULES) and ground-based B and V photometry. While the sinusoidal light variations are suggestive of a near-contact system, the stars have only shallow eclipse, so highly confident parametrization becomes challenging. We find: $M_1 = 2.1 \pm 0.1$, $M_2 = 0.4 \pm 0.1$; $R_1 = 3.20 \pm 0.02$, $R_2 = 1.48 \pm 0.01$; (${\odot}$); $T_{e1} = 10000 \pm 500$, $T_{e2} = 6200 \pm 200$ (K); distance = 162 $\pm$ 12 (pc). New times of minima for V486 Car have been examined, including recent observations from TESS. The role of the relatively significant O'Connell effect is examined. As well as the conspicuous asymmetry from the main effect of about 0.036 mag (V), a jitter, with amplitude of about 0.005 V mag and quasi-period of order $\sim$ 10 d is noticed. There is a tendency for such photometric excursions at one maximum to precede those at the other. As well, the O -- C data indicate the presence of a low mass star $\sim$0.3 M$_{\odot}$ in an orbit separated by a few AU from the close binary. More accurate and plentiful spectroscopic data would be requisite for further investigations. A brief discussion reviews possible approaches to understanding the system in the context of near-contact binary scenarios.
We examine the blue straggler star (BSS) populations of six old ($\geq$4 Gyr) open clusters: M67, NGC 188, NGC 6791, Berkeley 32, Berkeley 39, and Trumpler 19. We find that 50% of BSSs have color-magnitude diagram (CMD) locations corresponding to single stars in the final third of their main-sequence lifetimes. This build-up of BSSs near the terminal-age main sequence (TAMS) is primarily, but not solely, driven by more massive BSSs. Eleven of the BSSs have white dwarf companions with measured cooling ages; their evolution age distributions indicate that more massive BSSs typically form far from the zero-age main sequence, whereas lower mass BSSs can form at every evolutionary age. We show that inferred core helium amounts (above primordial) of late-evolution-age BSSs correspond to the core helium fused by cluster main-sequence stars near the turnoffs. We also find that the masses of asymptotic giant branch (AGB) mass-transfer BSSs require evolved main-sequence accretors and conservative mass transfer. These findings indicate that helium enrichment of progenitor accretors leads to the prevalence of BSSs near the TAMS. We further classify the evolutionary stages of the progenitor donors in M67 and NGC 188 and find mass transfer during the AGB accounts for at least half of the BSSs. We trace how the main-sequence binary population of NGC 188 evolves, and find that only 30-40% of interacting binaries create BSSs and that progenitor orbits must change to match current BSS periods.
Evolutionary tracks for $2-8M_\odot$ stars, with metallicities of $Z=0.014$, $0.006$, and $0.002$, including rotation, are computed with Modules for Experiments in Stellar Astrophysics (MESA). We study how rotation impacts the evolutionary properties of classical Cepheids. We examine whether rotation can offer a plausible explanation for the mass discrepancy problem when it is included in the evolutionary code using the fully diffusive approximation for rotationally induced mixing processes. We find that rotation barely influences the appearance and luminosity levels of the blue loops. While luminosity increases with increasing initial rotation rate, the increase does not exceed 0.04 dex, a fraction of the increase resulting from including the main sequence (MS) core overshooting of $0.2H_p$. As a consequence, rotation alone cannot resolve the mass discrepancy problem without simultaneously requiring significant MS core overshooting. Similar to the mass-luminosity relation, the period-radius and period-luminosity relations are barely affected by rotation, while the period-age relation predicts Cepheid ages to be only a few per cent longer compared with models without rotation. The predicted surface rotational velocities are too large compared with observations. These results are in contrast with those obtained with the Geneva code, which implements rotational mixing using the advective-diffusive scheme. In that approach, the luminosity levels of the loops are significantly higher, their luminosity extent increases, and the predicted rotation velocities are lower, compared with MESA models. The differences between the two approaches arise from significantly more efficient rotation-induced mixing during the MS evolution in models computed with the advective-diffusive scheme.
Magnetic switchbacks, large-amplitude, localized Alfvenic like rotations of the solar wind magnetic field, have been the subject of intensive investigation, with approximately 200 refereed papers published in the last decade. Yet, fundamental controversies persist regarding whether switchbacks can be described with Ideal MHD (magnetohydrodynamic) physics or Hall-MHD physics and whether their origin is at the solar surface or in the solar wind. To settle these controversies, we present Parker Solar Probe electric field measurements between 13 and 40 solar radii, which show that switchbacks have non-zero electric fields in the plasma frame, a finding that definitively settles the physics controversy by proving that switchbacks are Hall-MHD, not Ideal MHD, structures. Along with these electric fields, there are enhanced Poynting vectors having three components with similar magnitudes that exist only inside the switchbacks. These facts contradict the view of switchbacks as simple outward-propagating pulses. Together, they resolve one controversy by showing that switchbacks in the young solar wind are a non-MHD process. They contribute to the second (source) controversy by identifying switchbacks as sites of active, in-situ, evolution. These findings provide a new framework for understanding energy transport and dissipation in astrophysical plasmas.
We present non-local thermodynamic equilibrium radiative-transfer modeling of SN2023ixf during the nebular phase out to 1000d, using the same ejecta that matched its photospheric evolution, namely a partially stripped red-supergiant star of initially 15Msun whose terminal explosion yielded ejecta with 7-8Msun, kinetic energy of 1.2e51erg, and 56Ni mass of 0.05Msun, augmented with a cold dense shell (CDS) of 0.2Msun at 8000km/s. Interaction with circumstellar material persists at all epochs, powering the ultraviolet (UV) flux at all times, but dominating the optical only after ~600d. Matching the V-band light curve requires invoking both enhanced gamma-ray escape and dust formation after ~200d, first in the CDS and eventually in the inner ejecta as well. Depending on where they form relative to the dust, emission lines are uniformly attenuated or skewed with a blue-red asymmetry. Our models suggest a rising dust mass (chosen as an C-rich and Si-rich mixture) in the CDS and inner ejecta, possibly reaching 1e-4Msun at 700d, while an external cold dust component is required to match the mid-infrared emission. The UV radiation, largely unaffected by dust, is influenced by the emission and absorption from Fe lines, together with strong, blueshifted emission from Lyalpha and MgII2800, both present at >~200d and with a strengthening fractional flux thereafter. Optical-depth effects play a critical role for the UV flux, and most notably on Lyalpha whose strength depends strongly on the CDS structure (mass and extent) and the treatment of power injection. The CDS is continuously slowing down from 8000km/s at 112d to ~6500km/s at 998d, suggesting a growth in mass of several 0.1Msun. SN2023ixf shares many similarities with SN1993J at 1-3yr, but it is eventually fainter due to dust extinction and cooler (i.e., weak [NII] and no [OIII] lines) likely as a result of greater CDS and ejecta masses.
Binary stars undergoing mass transfer provide unique laboratories for testing stellar evolution. Here, we present a comprehensive photometric and spectroscopic analysis of the semi-detached system ET Cru. Using spectral disentangling, we independently determined the effective temperatures and chemical abundances of both components with high precision, including nine elements (eleven species). We find masses of $13.41\,M_\odot$ and $6.00\,M_\odot$ for the primary and secondary, respectively, with uncertainties of only $\sim$1.3%. The radii are $5.58\,R_\odot$ and $5.68\,R_\odot$, measured to within 0.4% and 0.5%. Surface gravities are constrained to better than 1%, while effective temperatures are determined to within 3-5%. The secondary exhibits extreme chemical anomalies, with severe carbon depletion and nitrogen enrichment far exceeding those reported in classical Algol systems. Multi-wavelength spectral energy distribution modelling yields a distance of $\sim$2.5 kpc, inconsistent with the $Gaia$ DR3 parallax, suggesting systematic astrometric uncertainties in the parallax distance. Together, these results establish ET Cru as a benchmark Algol-type binary, revealing direct spectroscopic evidence of deep CNO-cycle exposure in the donor and confirming the primary star as a rejuvenated gainer. ET Cru thus provides a chemically and dynamically illustrative case for understanding advanced binary interactions and the late evolutionary stages of massive-star evolution.
Procyon is a great system to probe stellar evolution of non-interacting binaries. We present an extensive grid of MESA (Modules for Experiments in Stellar Astrophysics) evolutionary tracks to constrain the evolution of Procyon A and B. We systematically vary the initial parameters of our grid anchored by precise dynamical masses and spectroscopic determinations of effective temperature (T$_{eff}$) and luminosity ($L$) to match the stars' positions on the H-R Diagram. Our goal is two-fold: (i) to quantify how the inferred system age and the progenitor mass of Procyon B depend on metallicity ($Z$), mixing length ($α$), and core overshoot ($β$), and (ii) to determine the best fitting model of Procyon B within a model-based initial-to-final mass relationship (IMFR) for hydrogen-deficient white dwarfs. Our best-fit models reproduce the observed properties for both components, yielding $M_\mathrm{A}=1.487 \pm 0.095$ M$_{\odot}$, $M_\mathrm{B}=0.592\pm 0.082$ M$_{\odot}$, a system age of $2.23 \pm 0.90$ Gyr, and a white dwarf cooling age of $1.20\pm0.49$ Gyr for Procyon B, consistent with independent determinations. Our results point to higher core overshoot than the standard adopted range, with the best fits ranging from $β=0.5-1.0$. From our model grid, we map Procyon B to the initial-to-final mass relationship for H-deficient white dwarfs in the $1.9\!-\!2.6$ M$_{\odot}$ progenitor range. Additionally, we implement the accretion of heavy metals onto the surface of the WD and fit our isotopic abundances to spectroscopic observations. We outline the physics used in our analysis.
The dependence of the stellar initial mass function (IMF) on star-formation environment, particularly at low metallicities and high redshifts, remains poorly constrained. Ultra-faint dwarf galaxies (UFDs) are local fossils of high-redshift galaxies hosting old, metal-poor populations, and their resolved stellar populations provide unique pathways to constrain the sub-solar IMF. We investigate the low-mass IMF in the Boötes I (Boo I) UFD with JWST/NIRCam, leveraging its capability to resolve over 10,000 stars reaching $\lesssim0.15 M_{\odot}$, obtaining one of the largest, deepest resolved stellar samples for UFDs. We explore three different functional forms of the IMF with machine learning and statistical techniques, combining forward modeling of synthetic color-magnitude diagrams with simulation-based inference. We find that a single power-law IMF fails to reproduce the observed luminosity function and also deviates from the canonical Salpeter IMF. Our best-fit broken power-law and lognormal IMF parameters are consistent with the Milky Way within 68% confidence level, providing evidence that star formation at metallicities as low as [Fe/H]$\approx-2.4$ follows a similar IMF as in the Milky Way. By treating Boo I as a local relic analogous to a high-redshift galaxy with a stellar mass of $\lesssim10^5 M_{\odot}$ at $z\gtrsim6$, our results provide evidence for the universality of the IMF across both local and high-redshift environments.
When considering how self-interaction affects an object's motion, it can be convenient to decompose the self-force into conservative and dissipative pieces. As a toy model for understanding such decompositions of the gravitational self-force, we consider objects that do not affect the spacetime, but are instead coupled to a nonlinear scalar field. There is then a standard splitting of the first-order scalar self-force into conservative and dissipative components. Multiple criteria can be used to obtain this splitting, all of which imply the same result. However, the implications of these criteria generically differ at higher orders. Demanding that any reasonable conservative sector be Hamiltonian, we identify multiple possible definitions of the conservative second-order self-force. Motivations for these possibilities and their properties are discussed and relevant Hamiltonians are obtained. We assume the existence of a three-point function with certain properties that is a generalization of the Detweiler-Whiting two-point function. These results apply to the two-body problem but are restricted to unbound scattering trajectories, due to infrared divergences that arise for bound orbits.
Early-time flux excesses in Type Ia supernovae (SNe~Ia) offer a unique insight into their progenitor systems and explosion mechanisms. Although individual early-excess events and larger searches have been reported, demographic studies remain limited by sample size. We present a systematic search for early-time excess emission in a volume-limited sample ($z<0.06$) of SNe~Ia based on the Zwicky Transient Facility Data Release 2 (ZTF DR2). Using ZTF $g$- and $r$-band light curves, we identify candidates showing early-excesses shortly after the explosion time, and we apply conservative coverage and quality requirements to build reliable ``excess'' and ``no-excess'' bump and no-bump catalogs. From an initial sample of 1547 SNe~Ia, our final catalogs contain 42 early-excess and 110 no-excess events. We compare the two populations using SN and host environment parameters from ZTF DR2 and quantify the differences using two-sample statistical tests. We find the strongest differences are in SN light-curve properties: early-excess events have larger SALT2 stretch $x_1$ ($7.91σ$) and larger $r$-band secondary-maximum flux $\mathcal{F}_{r_2}$ ($6.25σ$), while differences in SALT2 color $c$ are weak ($0.57σ$). Early-excess events also favor bluer $(g-z)_{\rm local}$ ($3.41σ$) and lower $\log_{10} (M_*/M_\odot)_{\rm local}$ ($2.73σ$). Our results connect early excesses with SNe~Ia diversity, and motivate further analyses of upcoming larger samples.
A recent analysis of photospheric radius expansion X-ray bursts from the low-mass X-ray binary 4U 1746-37 reported unusually small mass and radius estimates for the neutron star, suggesting it could be a quark star or quark-cluster star. Here, we propose an alternative interpretation: the star's mass and radius could be underestimated from significant blocking of the X-ray flux. Significant blocking factors ($\mathcal{B} \gtrsim 0.8$, reducing the observed flux to $\sim17\%$ of the intrinsic emission) permit neutron star parameters consistent with the canonical values: $M = 1.59 \pm 0.69 M_{\odot}$, $R = 13.0 \pm 5.45\,\mathrm{km}$, or $M = 2.12 \pm 1.08 M_{\odot}$, $R = 9.80 \pm 4.13\,\mathrm{km}$. The blocking factor, which varies with the photospheric radius, provides a natural explanation for the anomalously large peak-to-touchdown flux ratio ($\sim2.0$) and highlights the importance of accounting for geometric system configuration in neutron star mass--radius estimates.
Planets around pulsars were the first discovered exoplanets, found thanks to the extremely precise pulsar timing. Here we suggest that they could also be found through the radio emission generated by the pulsar-planet magnetospheric interaction. We present the results of special relativistic numerical simulations of planets in a pulsar wind of velocity $v=0.985~c$, corresponding to a Lorentz factor $γ=5.795$. Planets, modeled as a perfectly conducting solid surface in an external magnetic field originating from the pulsar, produce radio emission in the extended magnetic structure on the planet's nightside. We find that the planet around a known pulsar, PSR J0636+5129 b, could be detected via its radio emission. We outline the observational prospects for such objects.
The interaction of electrically charged particles with magnetic fields is a fundamental problem in several areas of physics. An example is the motion of energetic particles through a magnetized plasma. The most accurate and reliable way to explore theoretically the interactions between particles and fields is via test-particle simulations. In such simulations one creates the turbulent magnetic field and solves the Newton-Lorentz equation numerically by employing an integration scheme. In the current article we discuss exponential integrators and derive systematically from this the Rodrigues scheme as well as the famous Boris integrator. For an approach where one creates the magnetic field anew at each time step, both integrators are overall comparable. In theory the Rodrigues approach should be more accurate due to the fact that the occurring matrix exponential is evaluated without further approximations. Practically, both methods provide very similar results. It is argued in the current article that a Rodrigues based integrator is a very strong alternative because for the specific problem discussed here, it does not require longer computing times.