2026-02-10 | | Total: 104
We explore, using a state-of-the-art simulation code in 3D and to late enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state. Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear equations of state (EOS) for a 9-$M_{\odot}$/solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino ``light" curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multi-progenitor set of 3D supernova simulations for various EOSes to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest towards a fully quantitative theory of supernova observables.
We have used data collected by the Birmingham Solar-Oscillations Network (BiSON) to perform a helioseismic diagnosis of changes to the Sun's internal structure between four successive solar cycle minima, beginning with the minimum at the end of cycle 21 and ending with the recent minimum at the beginning of cycle 25. The unique duration of the BiSON database makes such a study possible. We used the low-degree BiSON p-mode frequencies to constrain structural changes between minima in the layers above $\approx 0.9 R_{\odot}$. We accomplished this by examining variations in the HeII ionisation zone signature; and by inverting the frequency differences to infer changes in the sound speed. Additionally, we employed frequency differences between various solar models that had subtle modifications to their internal structures to facilitate analysis of the observations. We find evidence for small, but marginally significant, changes in structure between different minima. The HeII signature was larger, and the sound speed in the range $\approx 0.93$ to $0.97 R_{\odot}$ was slightly higher, during the cycle 23/24 minimum, than during the other minima. The cycle 23/24 minimum was the deepest, as measured by proxies of global solar activity. These findings are consistent with magnetic flux levels having been lower in this minimum than the others, resulting in a higher gas pressure, higher temperatures, and higher sound speed. Our results demonstrate the potential of using asteroseismic data to perform similar analyses on other solar-type stars.
The Solar Hemispherical Helicity Rule(s) (SHHR) is a term used to represent a bias observed in proxies for the magnetic helicity in active regions at the solar surface. The SHHR states that predominantly negative magnetic helicity is observed in active regions in the northern hemisphere, whereas predominantly positive is found in the southern. The $Σ$-effect model of \cite{longcope1998flux} is one of the most cited models for the explanation of the SHHR. In this model, the magnetic structures derive the bias in their magnetic helicity from the kinetic helicity of the turbulent convection through which they travel, where the latter is handed owing to the rotational influence of the star. The original paper built an elegant mathematical model for the dynamics of thin flux tubes influenced by parameterized helical turbulence. Here, we attempt to explore the conceptual ideas of this original simplified model using fully-nonlinear, three-dimensional, Cartesian-domain simulations of isolated, finite cross-sectional, twisted magnetic flux structures rising though rotating, overshooting, turbulent compressible convection. We look for evidence of a correlation between the kinetic helicity content of the turbulence and the evolving magnetic helicity of the structures. We find little evidence of such a relationship, and do not even find any clear hemispheric dependence. Although these simulations are far from a perfect representation of the ideas, this work raises many questions about the potential efficacy of the $Σ$-effect in reality.
Granulation is the observable signature of convection in envelopes of low-mass stars, forming the background in stellar power spectra. While well-studied in evolved giants, granulation on the MS has received less attention. We here study and characterise granulation signatures of MS and SGB stars, extending previous studies of giants to provide a continuous physical picture across evolutionary stages. We analyse 753 Kepler short-cadence stars using a Bayesian nested-sampling framework to evaluate three background descriptions and compare model preferences. This yields full posterior distributions for all parameters, enabling robust comparisons across a diverse stellar sample. No universal preference between background models is found. Assuming a Gaussian oscillation envelope, $ν_\mathrm{max}$ estimates are sensitive to model misspecification, with the resulting systematics exceeding the formal uncertainties. The envelope width scales with $ν_\mathrm{max}$ across models and shows a dependence on effective temperature. Total granulation amplitudes in dwarfs broadly follow giant-based scalings, however a decoupling appears between the timescale of the primary granulation and the oscillations for MS stars cooler than the Sun. The prolonged granulation timescale is reproduced by 3D simulations of a K-dwarf, driven by reduced convective velocities due to more efficient convective energy transport in denser envelopes. The prolonged granulation timescale increases the frequency separation to the oscillation excess, potentially aiding seismic detectability, while the reduced convective velocities may influence the excitation of stellar oscillations and relate to the low amplitudes observed in cool dwarfs. Finally, we contribute a dataset linking granulation, oscillations, and stellar parameters, providing a foundation for future investigations into their interdependence across the HR diagram.
We report a significant increase in mid-infrared emission $\leq10$ $μ$m in a transitional disk. The 2024 JWST/MIRI observation of the disk around CVSO 1942 reveals flux increase by a factor of two at $λ\leq10$ $μ$m, compared to the near photospheric flux level observed with Spitzer/IRS in 2005. No significant change in flux at $\gtrsim15$ $μ$m is detected in the spectra. Comparing the MIRI/MRS spectrum and NEOWISE photometry, we found that this $\leq10$ $μ$m flux increase occurs on a timescale of 2 weeks and is consistent with the presence of warm (1,400 K), optically thick, large ($\gtrsim1$ $μ$m) dust grains near the dust sublimation radius. We propose that this rapid dust appearance may indicate in situ dust formation, possibly from planetesimal collisions in the inner disk.
PN K4-47 is a young planetary nebula that exhibits shock-excited bipolar lobes and a complex molecular environment, with the highest number of molecules detected within a planetary nebula. It has been questioned whether K4-47 is indeed a `typical' planetary nebula, or may be more exotic in nature. We examine this question using optical imaging and spectroscopy, and sub-millimetre and radio interferometry. Our observations spatially resolve the sub-millimetre environment of K4-47 for the first time. We find elongated CO (2-1) emission along a similar PA to the optical bipolar outflow. We derive a distance upper limit of 6 kpc to the source. The source hosts a fast ($\sim$350 km s$^{-1}$) bipolar optical outflow, and a slow (~50-60 km s$^{-1}$) molecular outflow along a similar PA. The outflow velocity indicates an age of 336 $\pm$ 119 yr. We also find that the excitation temperature and density of the atomic gas is~20 kK and 2800 cm$^{-3}$, respectively. The elemental and isotopic enrichment of K4-47 infers an AGB progenitor mass of 4-6 M$_{\odot}$, which corresponds to a white dwarf mass of ~1 M$_{\odot}$, following the initial-final mass relation for white dwarfs. We find that the core of K4-47 must contain 10$^{-2}$ M$_{\odot}$ of dust to explain the extinction, and that photoionisation alone cannot explain the excitation of the atomic gas. We instead require an additional heating mechanism, with shocks a likely scenario. It is likely that the progenitor star of K4-47 was a J-type carbon AGB star, which formed the molecular and dusty circumstellar environment. The bipolar outflow is later triggered, punching through the circumstellar environment, producing shocks, and shaping the environment into the elongated structure seen in the sub-millimetre. We therefore classify K4-47 as a genuine, if unusual, planetary nebula.
We report the identification of two previously unrecognized interstellar meteor candidates in the NASA CNEOS fireball database. Using the empirically calibrated low-discrepancy uncertainty model of Peña-Asensio et al.\ (2025) for post-2018 CNEOS velocity accuracy (1$σ$: speed 0.55~km~s$^{-1}$, right ascension 1.35$^\circ$, declination 0.84$^\circ$), we transform CNEOS velocity vectors to heliocentric orbits and assess interstellar candidacy via $10^{6}$-draw Monte Carlo simulations. Two post-2018 events have heliocentric speeds robustly exceeding the Solar System escape speed. CNEOS-22 (2022-07-28; 6.0$^\circ$S, 86.9$^\circ$W; eastern tropical Pacific) has $v_{\rm hel}=46.98$~km~s$^{-1}$, exceeding escape by $Δ= 5.18 \pm 0.60$~km~s$^{-1}$ ($z_Δ=8.7σ$), with interstellar speed $v_{\infty,\odot}=21.5$~km~s$^{-1}$. CNEOS-25 (2025-02-12; 73.4$^\circ$N, 49.3$^\circ$E; Barents Sea) has $v_{\rm hel}=45.63$~km~s$^{-1}$, exceeding escape by $Δ= 3.22 \pm 0.58$~km~s$^{-1}$ ($z_Δ=5.5σ$), with $v_{\infty,\odot}=16.9$~km~s$^{-1}$. For both events, none of $10^{6}$ realizations yield a gravitationally bound orbit ($p_{\rm bound} < 3\times 10^{-6}$). The adopted error model would need to underestimate the true uncertainties by factors of 5--9 for either candidate to be marginally bound.
Microquasars are primary candidates for Galactic PeVatrons, yet their collective contribution to the cosmic ray (CR) ``knee" remains poorly understood. We investigate this contribution by simulating anisotropic diffusive propagation through the Galactic magnetic field (GMF). Our results demonstrate that the GMF establishes a transport regime where magnetic connectivity between sources and the solar neighborhood determines the local flux. Active sources aligned with local GMF lines, such as Cygnus X-1, exhibit significant flux enhancements, while magnetically disconnected sources, such as V616 Mon, are strongly suppressed. By integrating source evolution with anisotropic transport, we show that the observed proton bump at the CR ``knee" is best reproduced by the cumulative contribution of microquasar remnants, which is often dominated by a few nearby or recent events, rather than the active ones alone. We find that a harder injection spectrum allows CRs from remnants to reproduce the PeV bump after propagation, as low-energy CRs have sufficient time to accumulate while high-energy CRs escape the Galactic plane. Our findings suggest that the integrated history of microquasar remnants, governed by the interplay of source age and magnetic connectivity, is the primary driver populating the observed CR ``knee''.
GJ 877 is a bright M dwarf in the solar neighborhood with two currently reported nontransiting exoplanets with periods of $9~\mathrm{d}$ and $21~\mathrm{d,}$ along with an additional unconfirmed signal at $50~\mathrm{d}$. We reanalyzed the system with 101 new HARPS and 12 new ESPRESSO radial velocities (RVs) secured with a cadence to confirm or refute the origin of the $50~\mathrm{d}$ signal. To do so, we searched for signals related to stellar activity in photometric data and spectroscopic indicators. We modeled the stellar activity in the RVs with Gaussian processes (GPs). With the Bayesian analysis, we confirmed a four-planet model, including the two previously known planets at periods of $9.2619\pm0.0005~\mathrm{d}$ and $21.784\pm0.004~\mathrm{d,}$ as well as two newly confirmed exoplanets: an Earth-mass planet, with a $4.42490\pm0.00014~\mathrm{d}$ period and a sub-meter-per-second amplitude, and a super-Earth with a $50.77\pm0.05~\mathrm{d}$ period located in the habitable zone (HZ). This super-Earth is the second closest planet in the HZ, after Proxima Cen b. We found an additional signal in a 2:1 resonance with the $4.4~\mathrm{d}$ planet at $2.21661\pm0.00010~\mathrm{d}$ with an amplitude of $0.37\pm0.09~\mathrm{m/s}$, which could be related to an additional planet. However, other explanations of its origin are also plausible. This signal remains a candidate, as further investigation is required to confirm its true nature. If the signal is caused by a planet, its minimum mass would be half that of Earth. We measured the stellar rotation period with the characteristic periodic timescale of the GP. We found a period of $38.7\pm0.5~\mathrm{d}$, which is consistent with the rotation period determined from photometry and other activity indices.
We present a data-driven reconstruction of the late-time expansion history and its implications for dark-energy dynamics. Modeling the reduced Hubble rate with a node-based Gaussian-process-kernel interpolant, we constrain the reconstruction using CC, Pantheon+ SNIa, BAO data from SDSS and DESI, transversal BAO data, and external $H_0$ priors (SH0ES and H0DN). Assuming GR at the background level, we map the reconstructed kinematics onto a dark-energy fluid and a scalar-field description, yielding the total potential and kinetic contributions that reproduce the inferred $H(z)$. To interpret the reconstruction, we consider both a minimal single-field model (canonical or phantom) and a two-field (quintom) system consisting of one canonical and one phantom scalar field (or families). Within the GR-based effective-fluid mapping, the inferred dark-energy density changes sign for all dataset combinations explored, transitioning from $ρ_{\rm DE}<0$ at higher redshift to $ρ_{\rm DE}>0$ toward the present, and defining a transition redshift $z_\dagger$ by $ρ_{\rm DE}(z_\dagger)=0$. A single canonical scalar cannot realize such a smooth evolution during expansion, whereas a phantom field or a two-field quintom framework can accommodate the required behavior; in particular, the two-field system permits smooth phantom-divide crossings at finite $ρ_{\rm DE}>0$ and distinguishes them from the separate notion of a density zero crossing. The reconstructed kinematics admit intermediate-redshift structure in some combinations, including hints of an additional accelerated-expansion interval around $z\sim 1.7$--$2.3$. The present-day equation of state remains close to a cosmological constant: combinations including supernovae give $w_0\simeq -1$, while combinations without supernovae but with an external $H_0$ prior show only a mild preference for $w_0<-1$ at the $\sim1.5$--$1.7σ$ level.
The black hole X-ray binary 4U 1755$-$338 underwent an outburst in 2020 after 25 years of quiescence. The comprehensive spectral analysis revealed that the system has a low interstellar neutral hydrogen column density of $0.34\pm0.01 \times$10$^{22}$ cm$^{-2}$. The outburst began with a low mass-accretion rate and was characterized as a low-luminosity outburst. The radius of the inner accretion disc remained constant throughout the outburst. Additionally, a growing neutral medium with constant density was detected in the local environment of 4U 1755$-$338.The hardness-intensity diagram (HID) did not follow the standard q-shaped pattern, indicating a non-canonical outburst. Instead, the HID showed a correlated evolution of hardness and source flux, suggesting a thermal disc origin of the flux. A wideband spectral analysis was performed using simultaneous NICER-NuSTAR data in two frameworks, based on kerrbb and bhspec. The results of bhspec (kerrbb) based modeling indicate that 4U 1755$-$338 is a high-inclination system, $67.44_{-3.03}^{+9.75}$ ($75.25_{-4.68}^{+5.59}$) degrees, and harbors a moderately spinning black hole with a spin parameter of $0.78_{-0.14}^{+0.02}$ ($0.50_{-0.43}^{+0.19}$) and a mass of $3.37_{-1.04}^{+0.45} (3.28_{-1.1}^{+1.7})M_{\odot}$ respectively. The inferred key parameters: black hole mass, spin, and system inclination are consistent across both modeling approaches. No reflection features were detected in the spectra of 4U 1755$-$338. The high spectral index, the blackbody nature ($L\propto T^4$) of the hardness ratio, the absence of reflection signatures, and the weak variability in the power density spectra indicate that the source remained in the high/soft state throughout the outburst.
We report results of one of the most sensitive pulsar surveys to date targeting the innermost region of the Galactic Center (GC) using the Robert C. Byrd Green Bank Telescope (GBT) at X-band (8--12GHz) using data from the Breakthrough Listen initiative. In total, we collected 9.5 hr of data covering the wider $\sim 8'$ diameter of the GC bulge, and 11 hr on the inner $1.4'$ region between 2021 May and 2023 December. We conducted a comprehensive Fourier-domain periodicity search targeting both canonical pulsars (CPs) and millisecond pulsars (MSPs), using constant and linearly changing acceleration searches to improve sensitivity to compact binaries. Assuming weak scattering, our searches reached luminosity limits of $L_{\rm min} \approx 0.14~{\rm mJy~kpc^{2}}$ for CPs and $L_{\rm min} \approx 0.26~{\rm mJy~kpc^{2}}$ for MSPs -- sensitive enough to detect the most luminous pulsars expected in the GC. Among 5,282 signal candidates, we identify an interesting 8.19 ms MSP candidate (DM of 2775 pc cm$^{-3}$), persistent in time and frequency across a 1-hr scan at a flux density of $S_{\rm min} \approx 0.007~{\rm mJy}$. We introduce a novel randomization test for evaluating candidate significance against noise fluctuations, including signal persistence via Kolmogorov-Smirnov tests and flux-vs-DM behavior. We are unable to make a definitive claim about the candidate due to a mixed degree of confidence from these tests and, more broadly, its non-detection in subsequent observations. This deepens the ongoing missing pulsar problem in the GC, reinforcing the idea that strong scattering and/or extreme orbital dynamics may obscure pulsar signals in this region.
The upper $100 \, \mathrm{m}$ to $150 \, \mathrm{m}$ of the polar ice sheet, called the firn, has a time-dependent density due to seasonal variations in the surface temperature and snow accumulation. We present RF simulations of an in-ice neutrino-induced radio source that show that these density anomalies create variations in the amplitude and propagation times of radio signals propagating through polar firn at an altitude of ${\sim}3000 \, \mathrm{m}$ above sea level. The received power from signals generated in the ice that refract within the upper ${\sim} 15 \, \mathrm{m}$ firn are subject to a seasonal variation on the order of 10\%. These variations result in an irreducible background uncertainty on the reconstructed neutrino energy and arrival direction for detectors using ice as a detection medium.
Binaries in the cores of globular clusters are known to prevent the gravitational collapse of the cluster, and simulations predict that the core of NGC 6397 contains a large number of white dwarfs (WDs), of which many are expected to be part of a binary system. In this work, we report the discovery of a compact binary system consisting of two WDs in the centre of the Galactic globular cluster NGC 6397. The system, known in the literature as NF1, was observed as part of a MUSE radial-velocity survey aiming at characterizing the binary population in the centre of NGC 6397. The spectral analysis of NF1 provides an effective temperature of 16000 K and a surface gravity (log g) of 5.72 (cgs), which is consistent with an extremely low-mass He-core WD nature. This is further supported by the mass of 0.23 +/- 0.03 Msun obtained from fitting the star's spectral energy distribution using its HST magnitude in various filters. The system has a circular orbit with a period of 0.54 days. The radial velocities show a large semi-amplitude of 200 km/s, implying a minimum mass of 0.78 Msun for the invisible companion, which is likely another WD, or a neutron star if the inclination of the system is smaller than about 50 deg. Some significant residuals in radial velocity remain with our best orbital solution and we tested whether a model with a third body can explain these deviations. While this possibility seems promising, additional measurements are needed to confirm whether the star is actually part of a triple system.
Given the wealth of data provided by Gaia and the upcoming PLATO mission, it is essential to improve stellar models to obtain accurate stellar ages. Our objective is to apply a mapping technique to estimate the age of a system and the initial chemical composition. We also evaluate the influence of observational uncertainties in mass and heavy-element mixtures on results. We applied an inverse calibration method to the evolution of a multiple stellar system, assuming that the stars share the same age and initial chemical composition. This approach determines age, the initial mass fractions of helium ($Y_{ini}$) and heavy elements ($Z_{ini}$), as well as the convective mixing-length parameters ($α_A $ and $α_B$). It uses the observed luminosities ($L_A$ and $L_B$), radii ($R_A$ and $R_B$), and surface chemical compositions ($Z/X_A$ and $Z/X_B$). We used the most recent observational data for $M$, $R$, $L$, and $[Fe/H]$ of $α$ Centauri A and B as input data for our method. We compared two assumptions for the $Z/X$ ratio, following the results for the solar composition. For an assumed high solar $Z/X_\odot =0.0245$, we obtain an age of $7.8 \pm 0.6$ Ga, $Y_{ini} = 0.284 \pm 0.004$, and $Z_{ini} = 0.0335 \pm 0.0015$. For a low solar $Z/X_\odot = 0.0181$, the derived age is $8.7 \pm 0.6$ Ga, $Y_{ini} = 0.267 \pm 0.008$, and $Z_{ini} = 0.025 \pm 0.002$. Observational errors in the stellar masses of $\pm$0.002 lead to an age error of 0.6 Ga. Overshooting of $0.05-0.20H_p$ at the boundary of the convective core increases the age by $0.6-2.1$ Ga. Models with higher $Z/X$ and radiative cores, with ages of $7.2-7.8$ Ga, appear preferable and show better agreement with the observed asteroseismic frequencies.
LOFAR observations of diffuse interstellar polarization at meter wavelengths reveal intricate polarized intensity structures with an unexpected correlation with neutral HI filaments that could not be reproduced in simulations with low cold neutral medium (CNM) abundance. We investigate whether MHD simulations of thermally bi-stable neutral interstellar medium, with a range of CNM fraction, can reproduce the properties of the 3C196 field, the high Galactic latitude test field. Using 50 pc simulations with varying levels of turbulence and compressibility, we generated synthetic 21 cm and synchrotron observations, including instrumental noise and beam effects, for different line-of-sight orientations relative to the magnetic field. We developed MOOSE, a code to generate synthetic synchrotron polarization and Faraday tomography. We also developed a metric based on the HOG algorithm, to quantify the relative contribution of cold and warm neutral medium structures to the Faraday tomographic data. The synthetic observations show levels of polarization intensity and RM values comparable to the 3C196 field, indicating that thermal electrons associated with the neutral HI phase can account for a significant fraction of the synchrotron polarized emission at 100-200 MHz. The simulations consistently reveal a correlation between CNM and Faraday tomographic structures that depends on turbulence level, magnetic field orientation, and observational noise, but only weakly on CNM fraction. We found slightly weaker CNM-Synchrotron polarized emission correlation level than observed in the 3C196 field. These results suggest that low-frequency polarimetric observations provide a valuable probe of magnetic-field morphology in the multi-phase Solar-neighborhood ISM, while simultaneously underscoring the need for improved modeling of the turbulent, multi-phase, and partially ionized interstellar medium.
We present Local Thermodynamic Equilibrium (LTE) estimates of the physical properties and star formation efficiency (SFE) of molecular gas in the Central Molecular Zone (CMZ), using new $^{12}$CO $J=2\to1$ observations from the James Clerk Maxwell Telescope. Combined with CHIMPS2 $^{12}$CO and $^{13}$CO $J=3\to2$, and SEDIGISM $^{13}$CO $J=2\to1$ data, we estimate a median excitation temperature of $T_{\rm ex} = 11$K for $^{13}$CO throughout the CMZ, with peaks exceeding $120$K in the Sgr B1/B2 complex. Cooler gas dominates around Sgr A and nearby clouds. We derive a median H$_{2}$ column-density of $N(\mathrm{H}2) = 2 \times 10^{22}$ cm$^{-2}$ and a total $^{13}$CO-traced gas mass of $M_{\rm gas} = 7 \times 10^6$ M$_\odot$, consistent with previous estimates when accounting for spatial coverage. The instantaneous SFE is assessed using Hi-GAL compact sources detected at 70-$μm$ and 160--500-$μm$. The 70-$μm$-bright SFE, tracing current star formation, is modest overall but elevated in Sgr B1/B2, the Arches cluster, and Sgr C. In contrast, the 160--500-$μm$ SFE, tracing cold pre-stellar gas, is more broadly enhanced, particularly in the dust ridge clouds and towards negative longitudes surrounding Sgr C. The contrasting distributions suggest an evolutionary gradient in SFE, consistent with a transition from dense, cold gas to embedded protostars. Our results imply that the CMZ may be enter a more active phase of star formation, with large reservoirs of gas primed for future activity.
We present a new method for estimating the total energy radiated by stellar flares in broad-band continua, which assumes a constant emitting area but incorporates a time-dependent temperature evolution. This physically motivated approach offers an alternative to the commonly used method that assumes a fixed flare temperature about of 10 000\,K and variable area. By allowing the temperature to vary over time while keeping the emitting area constant, our method captures more realistic flare behaviour. This time-dependent treatment of the flare temperature is supported by numerous solar observations, numerical simulations, and multiwavelength studies of active stars. We demonstrate that using peak flare temperatures estimated from a semi-empirical model grid, rather than assuming an ad-hoc flare temperature value, improves the accuracy of total energy estimates. Although the most precise results still require a multi-band photometry/spectroscopy or independently constrained flare temperatures, our method offers a practical and scalable solution for single-band observations. It is particularly well suited for main-sequence stars of spectral types K4 and later with known effective temperatures. Finally we discuss how the flare continuum behaves under varying chromospheric conditions. Our method improves flare energy estimates by incorporating a physically relevant time-dependent temperature evolution and empirically derived peak temperatures, rather than assuming a constant 10 000\,K value. This modification reduces systematic errors that can reach factors up to ten as compared to previous estimates. We proved this on a sample of 50,320 TESS flares.
The Kelvin-Helmholtz (KH) instability occurs in multiple heliospheric (solar-wind stream interfaces, planetary magnetospheres, cometary tails, heliopause flanks) and interstellar (protoplanetary disks, relativistic jets, neutron star accretion disks) environments. While the KH instability has been well-studied in the magnetohydrodynamic (MHD) limit, only limited studies were performed in the collisionless regime, which is conducive to development of anisotropic pressures. Collisionless plasmas are often described using the Chew Goldberger and Low (CGL) equations which feature an anisotropic pressure tensor. This paper presents a comprehensive analysis of the CGL version of the KH instability using linearised and numerical techniques. We find that the largest growth rates and the greatest incidence of magnetic effects occur in the MHD limit. In the large relaxation time CGL limit, part of the energy goes into the formation of pressure anisotropies, resulting in smaller amounts of energy being available for bending the field lines. Consequently, when we cross-compare CGL and MHD simulations that are otherwise identical, the current densities are largest in the MHD limit, and the largest magnetic islands also form in that limit. Early and late time formation of pressure anisotropies have also been studied. We also find that the strongest trend for forming intermittencies in the flow also occurs in the MHD limit. The paper also discusses possible consequences of our results for turbulence and reconnection in the heliosheath (the layer between the solar wind termination shock and the heliopause).
The solar chromosphere is a region where radiation plays a critical role in energy transfer and interacts strongly with the plasma. In this layer, strong spectral lines, such as the Lyman lines, contribute significantly to radiative energy exchange. Due to the long ionization/relaxation timescale, departures from LTE become significant in the chromosphere. Accurately modeling this layer therefore requires one to solve the non-LTE radiative transfer for the Lyman transitions. We present an updated version of the MURaM code to enable more accurate simulations of chromospheric hydrogen level populations and temperature evolution. In the previous extension, a non-LTE equation of state, collisional transitions of hydrogen, and radiative transitions of non-Lyman lines were already implemented in the code. Building on this, we have now incorporated radiative transfer for the Lyman lines to compute radiative rate coefficients and the associated radiative losses. These were used to solve the population and temperature evolution equations, rendering the system self-consistent. To reduce computational cost, a single-frequency approximation was applied to each line in the numerical solution of the radiative transfer problem. The extended model shows good agreement with reference solutions from the Lightweaver framework, accurately capturing the radiative processes associated with Lyman lines in the chromosphere. The extension brings the simulated hydrogen level populations in the deep chromosphere closer to detailed radiative balance, while those in the upper chromosphere remain significantly out of balance, consistent with the expected conditions in the real solar atmosphere. The extension enables the MURaM code to accurately capture chromospheric dynamics.
Two mysterious unresolved HI structures at velocities of +12 and -6 km/s were discovered in high resolution 21-cm survey data in the direction of a faint white dwarf star. Examination of the HI morphology in this area of sky shows that the star and HI features exists in a large cavity in interstellar HI surrounding the Upper Sco-Cen OB2 Association. The cavity may have been created by an ancient supernova. It is hypothesized that the pair of HI features and filamentary HI structure found in its immediate vicinity may be the remnants of a planetary nebula some 3 x 10^5 years old that have cooled to the point that the gas is neutral and emitting the 21-cm spectral line. This remnant has maintained the morphological characteristics of the original planetary nebula because it expanded into a volume of space relatively devoid of interstellar gas that would otherwise have absorbed any traces of the original nebula.
Context. The origin of the petrologic diversity observed in ordinary chondrites (OCs), the most common meteorites on Earth, remains debated. Competing models invoke either depth-dependent sampling of a single thermally stratified ("onion-shell") parent body or contributions from multiple distinct parent bodies. Aims. We aim to determine which of the two models is preferred for LL chondrites. These are unique among OCs in exhibiting a bimodal petrologic distribution, with most meteorites being LL3 or LL6. Methods. We compare the spectral and mineralogical properties of LL chondrites and corresponding LL-chondrite-like near-Earth objects (NEOs) with their possible sources in the main asteroid belt. We also model the thermal histories of the proposed parent bodies, based on revised estimates of parent-body sizes. Results. The spectral and mineralogical diversity of LL chondrites is consistent with contributions from the bright, S-type component of the Nysa family (NysaS) and the Flora family, with NysaS supplying mainly low-petrologic-type material and Flora higher-grade material. Unequilibrated, LL3 chondrites appear to originate exclusively from NysaS. Similarly, LL-chondrite-like NEOs form two distinct subpopulations consistent with origins in these same families. Conclusions. Our results favour multiple parent bodies for LL chondrites. The petrologic differences between the NysaS and Flora parent bodies indicate that planetesimal accretion within the OC reservoir extended over 0.5-0.7 Myr.
The search for the particle identity of dark matter (DM) continues to be a primary objective in modern physics. In this field, the sub-GeV mass range of DM detection remains a crucial yet challenging window. We investigate synchrotron self-Compton (SSC) emission from electrons and positrons produced by MeV-scale DM annihilation as a novel indirect detection channel. Focusing on the globular cluster Omega Centauri and the sensitivity of the Square Kilometre Array, we derive constraints on the annihilation cross section reaching $\langleσv\rangle \sim 10^{-30}\,\rm{cm}^{3}\,\rm{s}^{-1}$ in the tens-of-MeV range. Furthermore, constraints could even reach below $\langleσv\rangle \sim 10^{-32}\,\rm{cm}^{3}\,\rm{s}^{-1}$ for extreme parameter choices. Remarkably, even under deliberately conservative astrophysical assumptions, this channel outperforms existing indirect limits, establishing SSC emission as a robust probe of sub-GeV DM.
Based on very high-resolution resistive 2D and 3D magnetohydrodynamical (MHD) simulations of current sheets, our findings suggest that the answer to this question is likely no. In contrast, turbulence-mediated reconnection yields significantly faster reconnection rates - about an order of magnitude higher than the so-called universal rate for plasmoid-mediated reconnection in MHD flows ($V_\text{rec}/V_A \sim 0.01$). We conclude that turbulence-driven reconnection is the dominant mechanism responsible for fast reconnection and flares in systems such as accretion flows and relativistic jets in Active Galactic Nuclei (AGNs). In these environments, turbulence is driven by instabilities such as the magneto-rotational instability (MRI), Parker-Rayleigh-Taylor instability (PRTI), and current-driven kink instability (CDKI). Finally, we present 3D General Relativistic MHD simulations of accretion flows that confirm the crucial role of turbulence-mediated reconnection in AGN systems. These findings have important implications for understanding the origin of flares, particle acceleration, and the production of polarized radiation in these extreme environments.
We present a novel framework to estimate the sensitivity and discovery potential of IceCube to axion-like particles (ALPs) produced in core-collapse supernovae (CCSNe), covering ALP masses from 1 MeV to several hundred MeV. A key feature of this work is the explicit handling of the final-state leptons produced in ALP interactions with $^{16}$O nuclei and protons, which can generate Cherenkov light detectable in IceCube. These processes are being fully integrated into a detector-level simulation chain, enabling realistic detector signal modeling beyond existing estimates. The framework enables sensitivity forecasts for both direct detection and constraints based on time delays relative to the neutrino burst, across a range of ALP emission models. This approach may also extend to other MeV-scale dark sector particles. Preliminary sensitivity estimates are in progress and will be presented.