Abstract visibility
Abstract text size

Papers for Thursday, May 01 2025

Papers with local authors

Carlos Contreras Peña, Jeong-Eun Lee, Ho-Gyu Lee, Gregory Herczeg, Doug Johnstone, Hanpu Liu, Phillip W. Lucas, Zhen Guo, Michael A. Kuhn, Leigh C. Smith, Mizna Ashraf, Jessy Jose, Sung-Yong Yoon, Sung-Chul Yoon
0 votes
Paper 21 — arXiv:2504.21237
0 votes
Paper 21 — arXiv:2504.21237

Long-lasting episodes of high accretion can strongly impact stellar and planetary formation. However, the universality of these events during the formation of young stellar objects (YSOs) is still under debate. Accurate statistics of strong outbursts (FUors), are necessary to understand the role of episodic accretion bursts. In this work, we search for a population of FUors that may have gone undetected in the past because they either a) went into outburst before the start of modern monitoring surveys and are now slowly fading back into quiescence or b) are slow-rising outbursts that would not commonly be classified as candidate FUors. We hypothesise that the light curves of these outbursts should be well fitted by linear models with negative (declining) or positive (rising) slopes. The analysis of the infrared light curves and photometry of $\sim$99000 YSO candidates from SPICY yields 717 candidate FUors. Infrared spectroscopy of 20 candidates, from both the literature and obtained by our group, confirms that 18 YSOs are going through long-term outbursts and identifies two evolved sources as contaminants. The number of candidate FUors combined with previously measured values of the frequency of FUor outbursts, yield average outburst decay times that are 2.5 times longer than the rise times. In addition, a population of outbursts with rise timescales between 2000 and 5000 days must exist to obtain our observed number of YSOs with positive slopes. Finally, we estimate a mean burst lifetime of between 45 and 100 years.

All other papers

Louise Rebecca (1), Dominic Sebastian (1), C Sivaram (2), Kenath Arun (1) ((1) Department of Physics and Electronics, Christ University, Bangalore, India (2) Indian Institute of Astrophysics, Bangalore, India)

In view of the negative results from various dark matter detection experiments, we had earlier proposed an alternate theoretical framework through Modification of Newtonian Gravity (MONG). Here, the Poison's equation is modified by introducing an additional gravitational self-energy density term along with the usual dark energy density term. In this work we extend this model to account for the presence of low-density gas at high temperatures (10^8 K) in the intra cluster medium (ICM) by estimating the velocities to which particles will be subjected by the modified gravitational force. Considering that the ICM is under the influence of the cluster's gravity, particle velocities of the ions in the ICM must be balanced by the cluster's gravitational force. The particle velocities obtained for various clusters from their temperature profiles match the velocity produced by the MONG gravitational force. Thus, the increase in the gravitational potential at the outskirts of galaxies balances the thermal pressure of the ICM, maintaining hydrostatic equilibrium without invoking DM. The effect of MONG on the angular momentum of galaxies is also studied by obtaining a scaling relation between the angular momentum and the mass of a galaxy. MONG predicts a higher dependence on mass in comparison to the \Lambda-CDM model. This increased dependence on mass compensates for the halo contribution to the angular momentum. The angular momentum from MONG for galaxies from the SPARC database is compared to the halo angular momentum by a Chi-square fit technique. The correlation coefficient is found to be unity, showing a replicable result.

A limit on the expansion parameter a_hNR at which dark matter becomes non-relativistic has been obtained from the observed minimum halo mass hosting Milky Way satellites. This limit is in disagreement with measurements. In the present study, we attempt to understand this disagreement. We find that the limit does not include the following phenomena: non-linear regeneration of the power spectrum of density perturbations, the stripping of galaxy halos by neighboring galaxies, and baryons that act as cold dark matter. Considering these phenomena, we find that there is no longer a significant discrepancy between the limit and the measurements.

In this study, colour-magnitude relations (CMRs) for DA-type white dwarfs in the Sloan Digital Sky Survey (SDSS) photometric system were investigated. For this purpose, the SDSS data for 20,247 white dwarf stars, as provided in the study by Anguiano et al. (2017), were matched with the Gaia third data release (Gaia DR3) catalogue to obtain trigonometric parallax ($\varpi$) data. The SDSS $ugriz$ magnitudes of the stars were corrected for interstellar extinction using dust maps provided for the Milky Way, and distances from the Sun to the stars were calculated. The SDSS magnitudes were thus corrected for the effects of interstellar extinction. For the calibration of the stars, 5,516 white dwarf stars were selected, with apparent magnitudes brighter than $g_0=21$ mag and relative parallax errors measured to better than $\sigma_\varpi/\varpi=0.1$. Subsequently, three separate CMRs were derived for the absolute magnitudes $M_{\rm g}$, $M_{\rm r}$, and $M_{\rm i}$, each calibrated to two-colour indices. The coefficient of determination ($R^2$) of the obtained CMRs are highly reliable in the bf range of 0.86 to 0.95. Moreover, the standard deviations of the differences between the absolute magnitudes obtained from the relations and the original ones of the calibration stars range from 0.26 to 0.37 mag.

Kelsey S. Glazer, Tucker Jones, Yuguang Chen, Ryan L. Sanders, Marusa Bradac, Anthony J. Pahl, Alice E. Shapley, Richard S. Ellis, Michael W. Topping, Naveen A. Reddy

We present a spectral stacking analysis of galaxies at $z\geq6$ observed with the James Webb Space Telescope (JWST). We curate a sample of $64$ galaxies spanning redshifts $z_{\rm spec} = 6.0 - 9.4$ which have NIRSpec medium resolution data. The stacks achieve sufficient signal-to-noise to measure equivalent widths (EW) and velocity centroids ($v_{\rm{cen}}$) of low-ionization species (LIS) absorption features, transmitted Lyman-alpha ($\rm{Ly\alpha}$) emission, and nebular emission lines. Overall, we find our sample has weaker LIS absorption lines ($\rm{EW}(\rm{LIS}) \approx 1 Å$), smaller $v_{\rm{cen,LIS}} \approx -20 \pm 50~ \rm{km} \; \rm{s}^{-1}$, and significantly suppressed $\rm{Ly\alpha}$ emission ($\rm{EW}(\rm{Ly\alpha}) \approx 5~Å$), compared to similar studies undertaken at lower redshift. The weaker LIS absorption may suggest a lower covering fraction of HI and larger escape fraction of ionizing photons from our sample. Additionally, the smaller blueshifted $v_{\rm{cen,LIS}}$ indicates less prevalent or weaker outflows in $z>6$ galaxies. Stacking our sub-sample of $\rm{Ly\alpha}$ emitters (LAEs), we find high EW$(\rm{H}\beta) \approx 170 \pm 4~Å$ and a detection of nebular $\rm{C}\; \rm{IV}$ emission suggesting higher $\xi_{ion}$ in LAEs at $z>6$. This work showcases the enormous potential for stacked JWST spectra revealing properties of galaxies and their diffuse interstellar medium in the epoch of reionization.

Christopher Cain, Matthew McQuinn, Evan Scannapieco, Anson D'Aloisio, Hy Trac

The intergalactic medium (IGM) underwent intense heating that resulted in pressure disequilibrium in the wake of ionization fronts during cosmic reionization. The dynamical relaxation to restore pressure balance may have driven small-scale turbulence and, hence, the amplification of intergalactic magnetic fields. We investigate this possibility using a suite of $\approx 100$ pc resolution radiation-hydrodynamics simulations of IGM gas dynamics. We show that as the spatial resolution improves beyond that achieved with most prior studies, much of the IGM becomes turbulent unless it was pre-heated to $\gg 100~$K before reionization. In our most turbulent simulations, we find that the gas energy spectrum follows the expected $k^{-5/3}$ Kolmogorov scaling to the simulation's resolution, and the eddy turnover time of the turbulence is $< 1$ Gyr at $k \approx 1 ~$kpc$^{-1}$. Turbulence will grow magnetic fields, and we show that the fields grown by reionization-driven turbulence could explain lower limits on IGM B-field strengths from observations of TeV blazars.

R. A. P. Oliveira, R. Poleski, P. Mróz, A. Udalski, J. Skowron, M. Mróz, M. K. Szymański, I. Soszyński, K. Ulaczyk, P. Pietrukowicz, K. Rybicki, P. Iwanek, M. Wrona, M. Gromadzki

Gravitational microlensing depends primarily on the lens mass and presents a larger occurrence rate in crowded regions, which makes it the best tool to uncover the initial mass function (IMF) of low-mass stars in the Galactic bulge. The bulge IMF can be obtained from the luminosity function measured with the Hubble Space Telescope if one knows the statistics of binary stellar systems in the bulge. We aim to analyse a statistically significant number of binary-lens/single-source and single-lens/binary-source events, in order to explore the lower-mass end of the bulge IMF even in unresolved binary systems. This paper deals with events with clearly separated bumps and no caustic crossing or approach, whereas other types will be analysed in following works. A fully-automated approach in the search and modeling of binary events was implemented. Event detection was carried out with a modified version of the algorithm used in previous studies. Model fitting was carried out with Markov chain Monte Carlo and nested sampling methods, in order to find the most probable solution among binary lens or binary source models. We retrieved 107 binary events in Optical Gravitational Lensing Experiment (OGLE) light curves spanning ten years in 9 high-cadence and 112 low-cadence fields towards the bulge. Several criteria were applied to reduce false positives, resulting in 55 most likely binary lenses and 52 binary sources. The tools were effective detecting a bona-fide sample of binary events, with a distribution of Einstein timescales around 35-40 days and flat distributions for mass ratio and source flux ratio. After proper consideration of detection efficiency, the statistics for binary fraction and mass ratio will provide valuable constraints for the bulge IMF.

Filippo Santoliquido, Jacopo Tissino, Ulyana Dupletsa, Marica Branchesi, Jan Harms, Manuel Arca Sedda, Maximilian Dax, Annalena Kofler, Stephen R. Green, Nihar Gupte, Isobel M. Romero-Shaw, Emanuele Berti

The Einstein Telescope (ET), along with other third-generation gravitational wave (GW) detectors, will be a key instrument for detecting GWs in the coming decades. However, analyzing the data and estimating source parameters will be challenging, especially given the large number of expected detections - of order $10^5$ per year - which makes current methods based on stochastic sampling impractical. In this work, we use Dingo-IS to perform Neural Posterior Estimation (NPE) of high-redshift events detectable with ET in its triangular configuration. NPE is a likelihood-free inference technique that leverages normalizing flows to approximate posterior distributions. After training, inference is fast, requiring only a few minutes per source, and accurate, as corrected through importance sampling and validated against standard Bayesian inference methods. To confirm previous findings on the ability to estimate parameters for high-redshift sources with ET, we compare NPE results with predictions from the Fisher information matrix (FIM) approximation. We find that FIM underestimates sky localization errors substantially for most sources, as it does not capture the multimodalities in sky localization introduced by the geometry of the triangular detector. FIM also overestimates the uncertainty in luminosity distance by a factor of $\sim 3$ on average when the injected luminosity distance is $d^{\mathrm{inj}}_{\mathrm{L}} > 10^5~$Mpc, further confirming that ET will be particularly well suited for studying the early Universe.

We report the discovery of eleven high-velocity HI clouds at Galactic latitudes of 25-30 degrees, likely embedded in the Milky Way's nuclear wind. The clouds are detected with deep Green Bank Telescope 21 cm observations of a $3.2^\circ \times 6.2^\circ$ field around QSO 1H1613-097, located behind the northern Fermi Bubble. Our measurements reach $3\sigma$ limits on $ N_{\mathrm{HI}}$ as low as $3.1 \times 10^{17}$ cm$^{-2}$, more than twice as sensitive as previous HI studies of the Bubbles. The clouds span $-180 \leq v_{\mathrm{LSR}} \leq -90$ km/s and are the highest-latitude 21 cm HVCs detected inside the Bubbles. Eight clouds are spatially resolved, showing coherent structures with sizes of 4-28 pc, peak column densities of $\log(N_{\mathrm{HI}}/\mathrm{cm}^2) = 17.9\text{-}18.7$, and HI masses up to 1470 $M_\odot$. Several exhibit internal velocity gradients. Their presence at such high latitudes is surprising, given the short expected survival times for clouds expelled from the Galactic Center. These objects may be fragments of a larger cloud disrupted by interaction with the surrounding hot gas.

Amar Aryan, Ting-Wan Chen, Sheng Yang, James H. Gillanders, Albert K.H. Kong, S. J. Smartt, Heloise F. Stevance, Yi-Jung Yang, Aysha Aamer, Rahul Gupta, Lele Fan, Wei-Jie Hou, Hsiang-Yao Hsiao, Amit Kumar, Cheng-Han Lai, Meng-Han Lee, Yu-Hsing Lee, Hung-Chin Lin, Chi-Sheng Lin, Chow-Choong Ngeow, Matt Nicholl, Yen-Chen Pan, Shashi Bhushan Pandey, Aiswarya Sankar.K, Shubham Srivastav, Guanghui Sun, Ze-Ning Wang

The launch of the Einstein Probe (EP) mission has revolutionized the detection and follow-up observations of fast X-ray transients (FXTs) by providing prompt and timely access to their precise localizations. In the first year of its operation, the EP mission reports the discovery of 72 high signal-to-noise FXTs. Subjected to the visibility in the sky and weather conditions, we search for the optical counterparts of 42 EP-discovered FXTs from the Lulin observatory. We successfully detect the optical counterparts of 12 FXTs, and five of those are first discovered by us from the Lulin observatory. We find that the optical counterparts are generally faint (r>20 mag) and decline rapidly (>0.5 mag per day). We also find that 11 out of 42 FXTs had shown direct evidence of their association with Gamma-Ray Bursts (GRBs) through significant temporal and spatial overlapping. Furthermore, the luminosities and redshifts of FXTs with confirm optical counterparts in our observations are fully consistent with the faintest end of the GRB population. However, the non-detection of any associated optical counterpart with a significant fraction of FXTs suggests that EP FXTs are likely a subset of so-called `Dark' FXTs, similar to `Dark' GRBs. Additionally, the luminosities of two FXTs were also consistent with jetted tidal disruption events (TDEs). However, their luminosities differ significantly from those of typical supernova shock breakout or kilonova emissions. Thus, we conclude that a significant fraction of EP-discovered FXTs are associated with events having relativistic jets; either a GRB or a jetted TDE.

Thais Lemos, Rodrigo Gonçalves, Joel Carvalho, Jailson Alcaniz

Electromagnetic emissions from astrophysical sources at cosmological distances can be used to estimate the photon mass, $m_{\gamma}$. In this paper, we combine measurements of the dispersion measure ($\mathrm{DM}$) of fast radio bursts (FRB) with the luminosity distance from type Ia supernovae (SNe) to investigate update constraints on the photon rest mass. We derive the expression of $\mathrm{DM}$ dependence concerning a non-vanishing photon mass from a cosmological-model independent approach and constrain the parameter $m_{\gamma}$ from measurements of 68 well-localized FRBs and 1048 SNe data from the Pantheon compilation. We consider two scenarios for the baryon fraction in the intergalactic medium ($f_{\mathrm{IGM}}$): one where the value is fixed according to recent reports and another where it is treated as a free parameter, $f_{\mathrm{IGM}} = f_{\mathrm{IGM,0}}$. In the latter case, we find $m_{\gamma} = (29.4_{-15.5}^{+5.80}) \times 10^{-51}$ kg, at $1\sigma$ level. Our results also demonstrate an anticorrelation between $f_{\mathrm{IGM}}$ and $m_{\gamma}$, which highlights the importance of analyzing a larger sample of FRBs for a more comprehensive understanding of their properties.

Analyses of the galaxy N-Point Correlation Functions (NPCFs) have a large number of degrees of freedom, meaning one cannot directly estimate an invertible covariance matrix purely from mock catalogs, as has been the standard approach for the 2PCF and power spectrum. Instead, analyses use templates based on assuming a Gaussian Random Field density with the true, Boltzmann-solver-computed power spectrum. The resulting covariance matrices are sparse but have notable internal structure. To understand this structure better, we seek a fully analytic, closed-form covariance matrix template, using a power law power spectrum $P(k) \propto 1/k$. We obtain a simple closed-form solution for the covariance of the 2PCF, as well as closed-form solutions for the fundamental building blocks of the covariance matrices for the 3PCF, 4PCF, and beyond. We use our results to present a clearer picture of the covariance matrices' structure and sparsity, corresponding to triangular and non-triangular regions. This will be useful in guiding future NPCF analyses with spectroscopic surveys such as DESI, Euclid, Roman, and SPHEREx.

Sam Kumagai, Michael S. Vogeley, Miguel A. Aragon-Calvo, Kelly A. Douglass, Segev BenZvi, Mark Neyrinck

We present DeepVoid, an application of deep learning trained on a physical definition of cosmic voids to detect voids in density fields and galaxy distributions. By semantically segmenting the IllustrisTNG simulation volume using the tidal tensor, we train a deep convolutional neural network to classify local structure using a U-Net architecture for training and prediction. The model achieves a void F1 score of 0.96 and a Matthews correlation coefficient over all structural classes of 0.81 for dark matter particles in IllustrisTNG with interparticle spacing of $\lambda=0.33 h^{-1} \text{Mpc}$. We then apply the machine learning technique of curricular learning to enable the model to classify structure in data with significantly larger intertracer separation. At the highest tracer separation tested, $\lambda=10 h^{-1} \text{Mpc}$, the model achieves a void F1 score of 0.89 and a Matthews correlation coefficient of 0.6 on IllustrisTNG subhalos.

Patrick Horlaville, John J. Ruan, Michael Eracleous, Jaeden Bardati, Jessie C. Runnoe, Daryl Haggard

Supermassive black hole binaries (SMBHBs) at the centers of galaxies emit continuous gravitational waves (GWs) at nanohertz frequencies, and ongoing pulsar timing array (PTA) experiments aim to detect the first individual system. Identifying the exact host galaxy of a SMBHB detected in GWs is paramount for a variety of multi-messenger science cases, but it will be challenging due to the large number of candidate galaxies in the sky localization region. Here, we apply recent insights on the distinct characteristics of SMBHB host galaxies to archival galaxy datasets, to predict which nearby massive galaxies are most likely to host SMBHBs detectable by PTAs. Specifically, we use archival galaxy IFU surveys to search for nearby galaxies with distinct stellar kinematic signatures of SMBHB host galaxies, as informed by cosmological simulations. These distinct stellar kinematic signatures, including slow rotation and strong kinematic/photometric misalignments, are a hallmark of recent major galaxy mergers that led to the formation of SMBHBs in these galaxies. We produce a list of nearby massive galaxies that may currently host SMBHBs, ranked by a combination of their host galaxy stellar kinematic properties and their hypothetical GW strain. We discuss how our ranked list can be used (1) for targeted searches for individual sources of continuous GWs by PTAs, (2) to corroborate candidate SMBHBs identified through other approaches, and (3) to select candidate recoiling AGN and closely-separated (<100 pc) dual AGN for telescope follow-up confirmation.

Mars is the next milestone in human exploration. However, there are still several challenges that must be assessed to ensure appropriate conditions in a future settlement. Communications services will be essential for this task, providing not only a link between Earth and Mars but also supporting Martian weather forecasting and any potential rescue missions. These applications require a robust, high data rate communications network that allows for rapid response, remote sensing and public engagement. This research aims to study the feasibility of ground-to-satellite (and vice versa) optical communications during extreme Martian weather conditions, focusing on the link between a ground station on the surface of Mars and a satellite orbiting the planet. Long-lasting and expansive Martian dust storms, particularly common in the southern hemisphere, pose a considerable challenge when considering the feasibility of optical communications with Mars due to their significant impact in terms of signal attenuation and scattering. The methodology of this study is based on a computer simulation of the system featuring the characterisation of the Martian atmosphere and optical link to measure the attenuation and undesired effects suffered by the data signal when applying different environmental configuration parameters. The flexibility of the approach allows for the prediction of communications link quality in extreme cases such as global dust storms. The simulation is based on atmospheric data from the Mars Reconnaissance Orbiter's Mars Climate Sounder instrument and considers the recently launched Laser Communication Relay Demonstration (LCRD). The extreme conditions during dust storms in the southern polar-hood region lead to the proposal of a new communications network architecture to ensure connectivity during these events.

Extraterrestrial intelligences are speculated to surround stars with structures to collect their energy or to signal distant observers. If they exist, these most likely are megaswarms, vast constellations of satellites (elements) in orbit around the hosts. Although long-lived megaswarms are extremely powerful technosignatures, they are liable to be subject to collisional cascades once guidance systems start failing. The collisional time is roughly an orbital period divided by the covering fraction of the swarm. Structuring the swarm orbits does not prolong the initial collisional time as long as there is enough randomness to ensure collisions, although it can reduce collision velocities. I further show that once the collisional cascade begins, it can develop extremely rapidly for hypervelocity collisions. Companion stars or planets in the stellar system induce perturbations through the Lidov-Kozai effect among others, which can result in orbits crossing within some millions of years. Radiative perturbations, including the Yarkovsky effect, also can destabilize swarms. Most megaswarms are thus likely to be short-lived on cosmic timescales without active upkeep. I discuss possible mitigation strategies and implications for megastructure searches.

Rapid optical transient events can be hard to detect because of the limited number of photons they produce. I discuss a method of inferring the presence of fast, chaotic variability in photometry using the normalized autocorrelation function, what is called $g^{(2)}$ in quantum optics. The variability's signature is a bump in the function at short lags. No periodicity is needed for the method to work. Versions of this method are attested in stellar variability studies, but its uses in some other subfields apparently have not been realized. I calculate expected signal-to-noise ratios with shot noise and scintillation. This method could be used to find unknown phenomena, particularly sub-millisecond optical variability. I present simple models of three example use cases: a flickering artificial "lantern" near a host sun, optical microbursts from the Crab pulsar, and frequent irregular transits of a star by cometary bodies.

We present a novel way to fit solar flare X-ray spectra that offers more sensitivity to physical flare parameters than traditional approaches to spectroscopy. We decouple physically distinct emission types in solar flare X-ray spectra using timing behaviors, a technique we call time-decomposed spectroscopy. By fitting the shapes of particular time series to others across a time interval, we extract X-ray emission of distinct physical origins before any forward modeling. We perform spectroscopy on the original and time-decomposed spectra and find good agreement, with physical reasons for the few disagreements. In general, the time decomposition technique provides more precise results than those of traditional spectroscopy. The thermal and nonthermal energies are better constrained by more than an order of magnitude using the time decomposition approach, relative to traditional spectroscopy. We explain mathematically how and why the technique works using the multifractal formalism, and show that fractality is one way to choose appropriate light curves for this analysis. We speculate about applications to different wavelengths for solar data analysis, as well as applications to other physics subfields and domains.

Using three-dimensional general relativistic magnetohydrodynamic simulations with electron and proton thermodynamics, we probe the inner radial and vertical structure of weakly magnetized geometrically thin accretion discs around rapidly spinning black holes. We find that the thin, cold disc transitions to a thick, hot accretion flow at a radius dependent on the mass accretion rate. At high accretion rates, the disc truncates close to the innermost stable circular orbit $r\approx2r_g$, demonstrating that even in the canonical thin disc model, the plunging region should be treated with two-temperature physics. At intermediate accretion rates, the transition radius moves outward by a factor of two to $r\approx 5r_g$, forming a radiatively inefficient inner flow. The simulations also reveal extended cooling along the surface of the disc out to $\sim10r_g$, with 40% of the total cooling at intermediate accretion rates occurring above the disc body. These results have implications for X-ray binary state transitions and the physical origin of the X-ray corona.

Johanna Vazquez, Kathleen Barger, Frances Cashman, Andrew Fox, Bart Wakker, Felix Lockman, Alex Hill, Suraj Poudel, April Horton, Jaq Hernandez, Matthew Nuss, Alice Blake, Lauren Corlies, Molly Peeples

The Smith Cloud is a high-velocity cloud (HVC) on its final approach to the Milky Way that shows evidence of interaction with the Galaxy's disk. We investigate the metallicity and gas-phase chemical depletion patterns in this HVC using UV absorption-line observations toward two background QSOs taken with the Hubble Space Telescope (HST)/Cosmic Origin Spectrograph (COS) and H I 21-cm emission-line observations taken with Green Bank Telescope (GBT). We find evidence of silicon gas-phase depletion with [Si/S]$=\,-0.72^{+0.26}_{-0.24}$ and [Si/O]$_{3\sigma}\,\lesssim\,-0.05$, implying the presence of dust within the Smith Cloud. Because dust is galactic in origin, this HVC could trace the return leg of a Galactic fountain or a dwarf galaxy that passed through the Galactic plane.

H.E.S.S. Collaboration, F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, M. Bouyahiaoui, F. Bradascio, R. Brose, A. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, T. Bylund, S. Casanova, J. Celic, M. Cerruti, T. Chand, S. Chandra, A. Chen, J. Chibueze, O. Chibueze, G. Cotter, P. Cristofari, J. Damascene Mbarubucyeye, I.D. Davids, J. Devin, J. Djuvsland, A. Dmytriiev, K. Egberts, S. Einecke, S. Fegan, G. Fontaine, S. Funk, S. Gabici, J.F. Glicenstein, J. Glombitza, P. Goswami, G. Grolleron, L. Haerer, B. Heß, J.A. Hinton, W. Hofmann, T. L. Holch, M. Holler, D. Horns, Zhiqiu Huang, M. Jamrozy, F. Jankowsky, A. Jardin-Blicq, E. Kasai, K. Katarzyński, R. Khatoon, B. Khélifi, Nu. Komin, K. Kosack, D. Kostunin, A. Kundu, R.G. Lang, S. Le Stum, A. Lemière, M. Lemoine-Goumard, J.-P. Lenain, F. Leuschner, A. Luashvili, J. Mackey, D. Malyshev, V. Marandon, G. Martí-Devesa, R. Marx, A. Mehta, A. Mitchell, R. Moderski, L. Mohrmann, A. Montanari, M. de Naurois, J. Niemiec, P. O'Brien, L. Olivera-Nieto, E. de Ona Wilhelmi, M. Ostrowski, S. Panny, M. Panter, U. Pensec, G. Pühlhofer, M. Punch, A. Quirrenbach, S. Ravikularaman, M. Regeard

The results of the first extragalactic gamma-ray survey by the High Energy Stereoscopic System (H.E.S.S.) are presented. The survey comprises 2720 hours of very high-energy gamma-ray observations of the extragalactic sky, recorded with H.E.S.S. from 2004 up to the end of 2012. These data have been re-analysed using a common consistent set of up-to-date data calibration and analysis tools. From this analysis, a list of 23 detected objects, predominantly blazars, was obtained. This catalogue was assessed in terms of the source class populations that it contains. The level of source parameter bias for the blazar sources, probed by this observational dataset, was evaluated using Monte-Carlo simulations. Spectral results obtained with the H.E.S.S. data were compared with the \textit{Fermi}-LAT catalogues to present the full gamma-ray picture of the detected objects. Lastly, this unique dataset was used to assess the contribution of BL Lacertae objects and flat-spectrum radio quasars to the extragalactic gamma-ray background light at several hundreds of giga-electronvolts. These results are accompanied by the release of the high-level data to the astrophysical community.

AX J145732-5901 is an unidentified X-ray source discovered in the ASCA Galactic plane survey. Its extended nature and heavily absorbed X-ray spectrum suggest that AX J145732-5901 is a cluster of galaxies behind the Galactic plane. However, due to limited photon statistics, the spectral shape was not well examined. Using the results of the Galactic ridge X-ray emission and Cosmic X-ray background studies based on the Suzaku observations, we reanalyzed the ASCA data of AX J145732-5901. We confirmed that the source is more extended than the point spread function and the angular size is 14'x10'. The spectrum was heavily absorbed by interstellar matter equivalent to an N_{H} of ~10^{23} cm^{-2} and the emission line feature was confirmed. The spectrum was represented by a thin thermal plasma model with a temperature of 2.6 keV and a redshift of 0.12. Assuming the redshift value, the X-ray luminosity is calculated to be 2.6x10^{44} erg s^{-1} in the 1-10 keV energy band. The observational results indicate that AX J145732-5901 is a cluster of galaxies behind the Galactic plane.

We present a study of the third star orbiting around known contact eclipsing binary J04+25 using spectra from the LAMOST medium-resolution survey (MRS) and publicly available photometry. This is a rare case of a hierarchical triple, where the third star is significantly brighter than the inner contact subsystem. We successfully extracted radial velocities for all three components, using the binary spectral model in two steps. Third star radial velocities have high precision and allow direct fitting of the orbit. The low precision of radial velocity measurements in the contact system is compensated by large number statistics. We employed a template matching technique for light curves to find periodic variation due to the light time travel effect (LTTE) using several photometric datasets. Joint fit of third star radial velocities and LTTE allowed us to get a consistent orbital solution with $P_3=941.40\pm0.03$ day and $e_3=0.059\pm0.007$. We made estimations of the masses $M_{\rm 12,~3}\sin^3{i_3}=1.05\pm0.02,~0.90\pm0.02~M_\odot$ in a wide system and discussed possible determination of an astrometric orbit in the future data release of Gaia. Additionally, we propose an empirical method for measuring a period and minimal mass of contact systems, based on variation of the projected rotational velocity ($V\sin{i}$) from the spectra.

Ursa Major III/UNIONS 1 (UMa3/U1) is the faintest Milky Way satellite discovered to date, exhibiting a half-light radius of 3 $\pm$ 1 pc and an absolute V-band magnitude of +2.2 $\pm$ 0.4. Previous studies suggest UMa3/U1 is a dwarf galaxy, based on its large internal velocity dispersion and the improbability (indicated by dynamical cluster simulations) of its long-term survival if it were a dark-matter-free star cluster. In this paper, we model the evolution of UMa3/U1 as a star cluster using collisional N-body simulations that include a description of stellar evolution and the external tidal field of the Milky Way, with some simulations including primordial binaries. We find that UMa3/U1 has a substantial remaining lifetime of 2.7 $\pm$ 0.4 Gyr, primarily due to the retention of compact stellar remnants within the cluster. This retention is facilitated by mass segregation and the preferential loss of low-mass stars. Furthermore, we demonstrate that the observed large velocity dispersion of UMa3/U1 can be successfully reproduced. These results support the possibility that UMa3/U1 is a self-gravitating star cluster. Our simulations reveal that modelling UMa3/U1 as a dark matter free star cluster produces a markedly altered present-day mass function, driven by a strong depletion of low-mass stars. However, the degree of mass segregation among the visible stars is not statistically significant. We therefore recommend that future observations of UMa3/U1 and other very small Milky Way satellites focus on measuring their present-day mass functions to determine their nature.

Mariko Kimura, Takayuki Hayashi, Yuuki Wada, Wataru Iwakiri, Shigeyuki Sako, Martina Veresvarska, Simone Scaringi, Noel Castro-Segura, Christian Knigge, Keith C. Gendreau, Zaven Arzoumanian

We present our X-ray and optical observations performed by NICER, NuSTAR, and Tomo-e Gozen during the 2023 outburst in the intermediate polar GK Persei. The X-ray spectrum consisted of three components: blackbody emission of several tens of eVs from the irradiated white-dwarf surface, a source possibly including several emission lines around 1 keV, and multi-temperature bremsstrahlung emission from the accretion column. The 351.3-s white-dwarf spin pulse was detected in X-rays, and the observable X-ray flux from the column drastically decreased at the off-pulse phase, which suggests that the absorption of the column by the accreting gas called the curtain was the major cause of the pulse. As the system became brighter in optical, the column became fainter, the pulse amplitude became higher, and the energy dependence of pulses became weaker at $<$8~keV. These phenomena could be explained by the column's more pronounced absorption by the denser curtain as mass accretion rates increased. The blackbody and line fluxes rapidly decreased at the optical decline, which suggests the expansion of the innermost disk edge with decreasing accretion rates. The electron scattering or the column geometry may be associated with almost no energy dependence of high-energy pulses. The irradiated vertically-thick structure at the disk may generate optical QPOs with a period of $\sim$5700 s.

The cosmological tensions between early- and late-Universe probes, particularly the Hubble tension ($H_0$) and $S_8$ discrepancy, challenge the validity of the standard ${\rm{\Lambda}CDM}$ model. Motivated by these tensions, we perform a comprehensive joint analysis of three representative dark energy models - ${\rm{\Lambda}CDM}$, the Chevallier-Polarski-Linder (CPL) parametrization, and the Phenomenologically Emergent Dark Energy (PEDE) model - using the latest observational datasets: baryon acoustic oscillation (BAO) measurements from DESI Data Releases 1 and 2 (DR1/DR2), the Pantheon Plus sample of Type Ia supernovae (SNe Ia), and time-delay cosmography from TDCOSMO lensing. Our multi-probe approach breaks key degeneracies among cosmological parameters ($H_0$, $r_d$, and $M_B$) and provides robust constraints on dark energy dynamics. The CPL model yields a statistically significant improvement over ${\rm{\Lambda}CDM}$, with $\Delta \chi^2 \approx -3.6$ for DESI DR2+Pantheon Plus+TDCOSMO, favoring a quintessence-like behavior ($w_0=-0.87^{+0.045 }_{-0.045} $, $w_a=-0.41^{+ 0.28}_{-0.28}$ at 1$\sigma$ confidence level). In contrast, the PEDE model exhibits severe tension with observations, yielding $\Delta \chi^2 \approx +53.1$ (DR1) and $\Delta \chi^2 \approx +132.3$ (DR2), despite its potential to marginally alleviate the $H_0$ tension. DESI DR2 tightens constraints on dynamical dark energy by $\sim$40\%, reinforcing evidence for redshift-evolving $w(z)$. Remarkably, our results demonstrate that such data combination can achieve precision comparable to Planck CMB measurements for dynamical dark energy studies, while offering complementary advantages in probing the late-time universe. This synergy between different observational data significantly enhances our ability to constrain dark energy properties.

Primordial power spectra with low power at long wavelengths can alleviate lensing anomaly. However the extent to which data favours such a primordial spectra is not clear. In this work, we investigate power suppression and related mitigation of lensing anomaly with the help of phenomenological models which are valid over scales of interest. We consider simple extensions to nearly scale invariant power spectra such as those which includes running and running of running of spectral index. We perform Bayesian analysis of these models, which are agnostic about power suppression, with various data sets and show that data tend to choose parameters which leads to power suppression at low multipoles. We then analyse the significance of these findings using information criteria. Further, we investigate the ability of near-ultimate future CMB missions such as ECHO to put tighter constraints on these models. We conclude that we can make stronger conclusions about the presence of power suppression in the future by studying such simple phenomenological models.

Fenja Schweder, Sebastian Trujillo-Gomez, Kai Polsterer

Astronomical data is rich in volume, information and facets. Although this offers multiple research perspectives, processing the data remains a challenge. Infrastructures for analyzing, inspecting, exploring and communicating with data are mandatory. To address this issue, we introduce Jasmine, the JAvaScript Multimodal INformation Explorer. Jasmine allows users to open different data viewer modals that show a specific data point from a set. The viewer currently supports image data, as well as point cloud objects. Users can decide on which information about the data point they like to have displayed. Point clouds are interactive and allow for zooming, tossing, and turning. Picking a data point is enabled by providing a structured view of the set, arranged by a key property. This arrangement is achieved by autoencoding.

This paper describes advances in solar magnetography and developments in instrumental techniques of polarimetry and spectroscopy made at Paris-Meudon observatory in the second half of the twentieth century. The adventure started from Lyot expertise and extended progressively to the measurement of vector magnetic fields using various and improving polarimetric techniques (such as beam exchange or grid) or new spectroscopic methods (such as the MSDP imaging slicer), at Meudon and Pic du Midi, ending by the achievement of the state-of-the-art optimized and polarization free telescope THEMIS in 1999.

Rapid uniformly-rotating neutron stars are expected to be formed for instance in the collapse of some massive stars, the accretion of some compact object binaries, and some double neutron star mergers. The huge amount of the rotational energy has been widely believed to be the source of some cosmic gamma-ray bursts and superluminous supernovae. Benefited from the tight constraints on the equation of state of the neutron star matter set by the latest multi-messenger data, the chiral effective field theory ($\chi$EFT) and perturbative quantum chromodynamics (pQCD), here we present the maximum gravitational mass as well as the kinetic rotational energy for a neutron star at a given spin period. Our nonparametric EOS analysis reveals that the critical Keplerian configurations ($\Omega_{\rm kep}^{\rm crit}=1.00\pm0.07\times 10^{4}~ {\rm rad/s}$) can sustain maximum gravitational masses of $M_{\rm kep}^{\rm crit}=2.76^{+0.11}_{-0.09} M_\odot$ with corresponding rotational energy reaching $E_{\rm rot,kep}^{\rm crit}=2.38^{+0.25}_{-0.24}\times 10^{53}$. However, the maximum rotational energy that can be feasibly extracted from a neutron star is limited to $1.40^{+0.15}_{-0.13}\times 10^{53}$ erg, which holds for a baryon mass of $2.68^{+0.10}_{-0.09}M_\odot$. All these parameters, obtained via the nonparametric reconstruction of the equation of state, are at the $68.3\%$ confidence level and the adoption of a quarkonic model yields rather similar results. These findings are found to have already set some intriguing constraints on the millisecond magnetar interpretation of some exciting data.

Coronal heating has puzzled solar physicists for decades. The question of why the Sun's upper atmosphere is significantly hotter than its lower atmosphere remains a key mystery. It is commonly believed that the source of coronal heating comes from the Sun's magnetic field, and more complex magnetic dynamics is more efficient in heating. In an earlier work we studied the secondary (or finer-scale) flux emergence identified in five ephemeral regions (ERs), selected during the last solar minimum (Yang et al. 2024). Here we further explore the atmospheric response to the secondary flux emergences (SFEs) that were identified in the first paper. We further reveal that approximately 80 percent of the 172 identified SFEs are associated with atmospheric heating events. The heating is most likely associated with magnetic reconnection involved in the SFE. Overall, a solar quiet region is heated by several hundred thousand degrees, during flux emergence of an ER.

A popular class of models for interpreting quasi-periodic X-ray eruptions from galactic nuclei (QPEs) invoke collisions between an object on an extreme mass ratio inspiral (EMRI) and an accretion disk around a supermassive black hole. There are strong links between QPE systems and those disks which formed following a tidal disruption event (TDE), and at least two events (AT2019qiz and AT2022upj) are known to have occurred following an otherwise typical TDE. We show that the fact that these disks were formed following a TDE strongly constrains their properties, more so than previous models have assumed. Models based on steady-state AGN-like disks have mass contents which grow strongly with size $M_{\rm disk}\propto R_{\rm out}^{7/2}$ and do not conserve the mass or angular momentum of the disrupted star. A very different scaling must be satisfied by a TDE disk in order to conserve the disrupted stars angular momentum, $M_{\rm disk} \propto R_{\rm out}^{-1/2}$. These constraints substantially change the predicted scaling relationships between QPE observables (luminosity, duration, energy, temperature) and the QPE period. They also allow QPE observables to be written in terms of the properties of the two stars assumed to be involved (the one tidally disrupted and the one on an EMRI), making plausibility tests of these models possible. We show that these modifications to the disk structure imply that (i) QPEs cannot be powered by collisions between an orbiting black hole and a TDE disk, (ii) QPEs also cannot be powered by collisions between the surface of a stellar EMRI and a TDE disk. A framework in which the collisions are between a TDE disk and a star which has puffed up to fill its Hills sphere with a trailing debris stream (as seen in recent simulations) cannot be ruled out from the data, and should be the focus of further study.

The nature of dark matter is one of the most fundamental questions in cosmology. Using the cosmic microwave background (CMB), type Ia supernova (SN) and DESI's new measurements of baryon acoustic oscillations (BAO), we find the robust $\sim2\,\sigma$ evidences of the evolution of dark matter in the dynamical dark matter (DDM) model, $\omega_{dm}(a)=\omega_{dm0}+\omega_{dma}(1-a)$. Based on CMB data, we find a very strong linear relation $\omega_{dma}=-\omega_{dm0}$, inducing the single-parameter DDM model, $\omega_{dm}(a)=\omega_{dm}a$, where the $\sim2\,\sigma$ DDM evidences is well captured and even strengthened. We demonstrate that there are beyond $2\,\sigma$ evidences of the coexistence of DDM and dynamical dark energy using the combinations of CMB, DESI BAO and Pantheon+ SN data. In such models, at a beyond $5\,\sigma$ confidence level, we verify that the universe remains in a matter-dominated state for a substantial period in the past, accelerate in the distant future and finally becomes completely dominated by dark matter. We propose that the ultimate fate of the universe is the ``Super Rip'' induced by dark matter with an extremely negative pressure. Our findings fundamentally challenge the prevailing understanding of cosmic acceleration and deepen our insight into the universe's evolution.

Fenfen Yin, Liang Xiao, Wenying Du, Zhujun Jiang, Zhiwei Min, Jaime Forero-Romero, Jiacheng Ding, Le Zhang, Xiao-Dong Li

We investigate the potential of $\beta$-cosmic-web weighted angular correlation functions to improve the cosmological constraints. Using SDSS DR12 CMASS-NGC galaxies and simulated mock catalogs with $\Omega_m$ varying in the range of 0.25-0.40, we quantify the discriminative power of different statistics via $\Delta \chi^2$ measurements. Our results demonstrate significant improvements when incorporating weighted statistics. Especially, adding the $\bar{D}_{\rm nei}$-weighting statistics enhances $\Delta \chi^2$ by 39%-130%, while adding the $1/\bar{D}_{\rm nei}$-weighted statistics yields 229%-336% gains over solely using the traditional angular statistics. These findings align with 3D correlation function studies \cite{Yin+etal+2024}, confirming the superior performance of $\beta$-cosmic-web weighted statistics. The thin-redshift-slice approach makes our method particularly relevant for slitless surveys (such as Euclid, CSST) where redshift errors challenge traditional 3D analyses. This work also establishes the first theoretical framework for marked statistics in 2D angular clustering.

Spencer J. Magnall, Christian Ecker, Luciano Rezzolla, Paul D. Lasky, Simon R. Goode

Gravitational-wave astronomy shows great promise in determining nuclear physics in a regime not accessible to terrestrial experiments. We introduce physics-informed priors constrained by nuclear theory and perturbative Quantum Chromodynamics calculations, as well as astrophysical measurements of neutron-star masses and radii. When these priors are used in gravitational-wave astrophysical inference, we show a significant improvement on nuclear equation of state constraints. Applying these to the first observed gravitational-wave binary neutron-star merger GW170817, the constraints on the radius of a $1.4\,M_\odot$ neutron star improve from $R_{1.4} ={12.54^{+1.05}_{-1.54}} \, {\rm km}$ to $R_{1.4} = 12.11^{+0.91}_{-1.11} \,{\rm km}$ and those on the tidal deformability from $\tilde{\Lambda}_{1.186} < 720$ to $\tilde{\Lambda}_{1.186} = 384^{+306}_{-158}$ ($90\%$ confidence intervals) at the events measured chirp mass $\mathcal{M}=1.186\,M_\odot$. We also show these priors can be used to perform model selection between binary neutron star and neutron star-black hole mergers; in the case of GW190425, the results provide only marginal evidence with a Bayes factor $\mathcal{BF}=1.33$ in favour of the binary neutron star merger hypothesis. Given their ability to improve the astrophysical inference of binary mergers involving neutron stars, we advocate for these physics-informed priors to be used as standard in the literature and provide open-source code for reproducibility and adaptation of the method.

Chin-Ping Hu, Zorawar Wadiasingh, Wynn C. G. Ho, Matthew G. Baring, George A. Younes, Teruaki Enoto, Sebastien Guillot, Tolga Guver, Marlon L. Bause, Rachael Stewart, Alex Van Kooten, Chryssa Kouveliotou

During the 2022 outburst of SGR 1935+2154, a Fast-Radio-Burst-like event (FRB 20221014A) and X-ray activities occurred between two spin-up glitches, suggesting these glitches may connect to multiwavelength phenomenology. However, the mechanisms altering the magnetar's magnetosphere to enable radio emission remain unclear. This study presents high-cadence NICER and NuSTAR observations revealing spectral changes in burst and persistent emission. Hardness ratio and spectral analysis reveal significant changes during an "intermediate flare" 2.5 hours before FRB 20221014A. This 80-second flare, releasing $>(6.3\pm0.2)\times10^{40}$ erg, coincides with a rapid spectral softening in both burst and persistent emission and a notable decrease in burst occurrence rate. The intermediate flare is bright enough to be detected if placed at a few Mpc, and would appear as a fast X-ray transient. This implies that the connection between magnetar X-ray activity and FRBs can be observed in the local Universe. Post-flare burst spectra peak near 5 keV, resembling the characteristics of the FRB-associated X-ray burst of 2020. Such change persisted for a few hours, implying magnetospheric evolution on similar timescales. However, no radio emission was detected from post-flare bursts, suggesting that FRB emission requires conditions beyond peculiar short bursts. The burst waiting times exhibit a broken power-law distribution, likely resulting from contamination by enhanced persistent emission. Although the bursts appear randomly distributed in the spin phase, the hardness ratio profile as a function of spin phase follows that of the persistent emission, indicating that X-ray bursts originate at low altitudes.

Matúš Labaj, Sean M. Ressler, Michal Zajaček, Tomáš Plšek, Bart Ripperda, Florian Peißker

Intermediate-mass black holes (IMBHs) are a missing link in black hole demographics, with only tentative observational evidence to date. Dense stellar clusters such as IRS 13E near the Galactic Center are promising IMBH hosts, where accretion is likely driven by winds from nearby Wolf-Rayet (WR) stars. Yet, the dynamics of such wind-fed systems remain largely unexplored. We investigate how high-velocity stellar winds, magnetic fields, and metallicity-dependent radiative cooling influence gas dynamics and black hole accretion in compact WR clusters. Using three-dimensional (magneto)hydrodynamic simulations, we model each WR star as a source of mass, momentum, energy, and magnetic flux, and include a cooling function that depends on chemical abundance. We compare isotropic versus disk-like stellar distributions to explore the impact of cluster geometry. Across all models, we find that the accretion rate onto the IMBH is suppressed by up to five orders of magnitude relative to the total stellar mass-loss rate. Turbulent, shock-heated outflows driven by wind-wind collisions dominate the flow, expelling most injected gas. While enhanced cooling in high-metallicity runs promotes the formation of dense clumps, these structures are typically unable to reach the black hole. The system's integrated X-ray luminosity is dominated by colliding WR winds, masking the IMBH's radiative signature. Accretion occurs in short-lived, quasi-periodic episodes triggered by close stellar passages, but even these flares remain difficult to detect against the luminous wind background. Our results naturally explain the low detectability of IMBHs in compact WR clusters and provide theoretical predictions to guide future X-ray and infrared observational strategies.

Birendra Chhotaray, Gaurava K. Jaisawal, Sachindra Naik, Arghajit Jana

We present the results of the broadband timing and spectral analysis of the poorly understood SMC pulsar RX J0032.9-7348 (= SXP 7.02) using NuSTAR and NICER observations during its X-ray brightening in 2024. Our timing analysis revealed a pulsation period of approximately 7.02 s in the X-ray light curve. The pulse profile obtained in the broad energy range is double-peaked and asymmetric in nature and shows moderate variation with the energy. An absorbed power-law model describes the 0.5-8 keV NICER spectra well. The 3-50 keV NuSTAR spectrum is best described with an absorbed power-law modified with a high-energy cutoff model. We find no evidence of iron or cyclotron line features in the energy spectrum. During our observation period, the 0.5-50 keV luminosity varies in the range of $\sim 8\times10^{36} - 4\times10^{37}$ erg s$^{-1}$. We also discuss the dependence of spectral parameters on the rotational phase of the pulsar through phase-resolved spectroscopy.

F. Onori, M. Nicholl, P. Ramsden, S. McGee, R. Roy, W. Li, I. Arcavi, J. P. Anderson, E. Brocato, M. Bronikowski, S. B. Cenko, K. Chambers, T. W. Chen, P. Clark, E. Concepcion, J. Farah, D. Flammini, S. González-Gaitán, M. Gromadzki, C. P. Gutiérrez, E. Hammerstein, K. R. Hinds, C. Inserra, E. Kankare, A. Kumar, L. Makrygianni, S. Mattila, K. K. Matilainen, T. E. Müller-Bravo, T. Petrushevska, G. Pignata, S. Piranomonte, T. M. Reynolds, R. Stein, Y. Wang, T. Wevers, Y. Yao, D. R. Young

We present the results from our multi-wavelength monitoring campaign of the transient AT2022wtn, discovered by the Zwicky Transient Facility in the nucleus of SDSSJ232323.79+104107.7, the less massive galaxy in an active merging pair with a mass ratio of ~10:1. AT2022wtn shows spectroscopic and photometric properties consistent with a X-ray faint N-strong TDE-H+He with a number of peculiarities. Specifically, a 30-days long plateau at maximum luminosity, a corresponding dip in temperature and the development of a double-horned N III+ He II line profile. Strong and time-evolving velocity offsets in the TDE broad emission lines and the detection of a transient radio emission, indicate the presence of outflows. Overall, the observed properties are consistent with the full disruption of a low-mass star by a ~10$^{6}$ M$_{\odot}$ SMBH followed by an efficient disk formation and the launch of a quasi-spherical reprocessing envelope of fast expanding outflowing material. The observed differences between the He II and the Hydrogen and N III lines can be explained either with a spatial separation of the lines emitting region or with a late-time reveal of shocks from the returning debris streams, as the photosphere recedes. Finally, we present an extensive analysis of the hosting environment and discuss the implications for the discovery of two TDEs in interacting galaxy pairs, finding indication for an over-representation of TDEs in these systems. The AT2022wtn host galaxy properties suggest that it is in the early stages of the merger, therefore we may be witnessing the initial enhanced rate of TDEs in interacting galaxies before the post-starburst phase.

Jingbo Sun, Hengxiao Guo, Wenwen Zuo, Paulina Lira, Minfeng Gu, Philip G. Edwards, Shu Wang, Jamie Stevens, Tao An, Samuzal Barua, Zhen-yi Cai, Haicheng Feng, Alok C. Gupta, Luis C. Ho, Dragana Ilić, Andjelka B. Kovačević, ShaSha Li, Mar Mezcua, Luka Č. Popović, Paula Sánchez-Sáez, Mouyuan Sun, Rongfeng Shen, Vivian U, Oliver Vince, Junxian Wang, Xuebing Wu, Zhefu Yu, Zhenya Zheng

The search for robust evidence of intermediate-mass black holes (IMBHs) is crucial for understanding black hole seeding process and the formation of supermassive black holes in the early Universe. NGC 4395 and POX 52 are two prototypical IMBH hosts, both exhibiting multi-line evidence of low-mass black hole activity. Here, we report the first detection of mid-infrared (MIR) lags in response to optical variability, with measurements of $3.0^{+2.4}_{-1.9}$ days for NGC 4395 and $35.2^{+14.2}_{-11.7}$ days for POX~52 at $3.4$ $\mu$m, respectively, using archival optical data and observations from the Wide-field Infrared Survey Explorer (WISE). This detection provides the first reverberation evidence of low-mass black hole activity in POX 52. The time lags of these two low-mass, low-luminosity active galactic nuclei (AGNs) generally follow the extent of the $R_{\rm dust}-L_{\rm 5100}$ relation found in higher-mass AGNs. Based on an empirical relation between the broad-line region and dusty torus size, we constrain the black hole mass of POX 52 to log($M_{\rm BH}$/$M_\odot$) = 5.5 $\pm$ 0.37 (systemic and statistical errors), confirming its IMBH nature. Furthermore, long-term optical continuum monitoring of POX 52 reveals a mild inter-band lag of $\lesssim$ 1 day. However, no significant intranight variability was detected during its one-night, high-cadence monitoring, which we attribute to the longer duty cycle of fast variability in POX 52 compared to that in NGC 4395.

In Paper I (Rowan-Robinson 2024), models derived in 2009 to fit mid-infrared (8-24 micron) source counts from the IRAS, ISO and Spitzer missions, were found to provide an excellent fit to deep counts at 7.7-21 mu with JWST, demonstrating that the evolution of dusty star-forming galaxies is well understood. Here the treatment of optical spectral energy distributions (SEDs) is improved and the counts are extended to 5.6 mu and optical wavelengths. The models proved a good fit to the latest, deeper, JWST counts. The models are also extended to radio and X-ray wavelengths. Predicted redshift distributions are given for a range of wavelengths and flux-densities.

In an era when we are charting multiple planets per system, one might wonder the extent to which "missing" (or failing to detect) a planet can skew our interpretation of the system architecture. We address this question with a simple experiment: starting from a large, homogeneous catalog, we remove planets and monitor how several well-defined metrics of the system architecture change. We first perform this test on a catalog of observed exoplanets. We then repeat our test on a catalog of synthetic planetary systems with underlying hyperparameters that have been fit to reproduce the observed systems as faithfully as possible (though imperfectly). For both samples, we find that the failure to detect one or more planets tends to create more irregularly spaced planets, whereas the planet mass similarity and coplanarity are essentially unaffected. One key difference between the synthetic and observed data sets is that the observed systems have more evenly spaced planets than the observation-bias-applied synthetic systems. Since our tests show that detection bias tends to increase irregularity in spacing, the even spacing in the observed planetary systems is likely astrophysical rather than the result of the Kepler missions' inherent detection biases. Our findings support the interpretation that planets in the same system have similar sizes and regular spacing and reinforce the need to develop an underlying model of planetary architectures that reproduces these observed patterns.

Luis Welbanks, Matthew C. Nixon, Peter McGill, Lana J. Tilke, Lindsey S. Wiser, Yoav Rotman, Sagnick Mukherjee, Adina Feinstein, Michael R. Line, Sara Seager, Thomas G. Beatty, Darryl Z. Seligman, Vivien Parmentier, David Sing

Claims of detections of gases in exoplanet atmospheres often rely on comparisons between models including and excluding specific chemical species. However, the space of molecular combinations available for model construction is vast and highly degenerate. Only a limited subset of these combinations is typically explored for any given detection. As a result, apparent detections of trace gases risk being artifacts of incomplete modeling rather than robust identification of atmospheric constituents, especially in the low signal-to-noise regime. We illustrate these challenges using the sub-Neptune K2-18~b, where recent claims of a potential biosignature detection vanish when the considered model space is expanded. We show that numerous alternative models without potential biosignature gases provide equivalent or better fits to the observations. We demonstrate that the significance of a claimed detection relies on the choice of models being compared, and that model preference does not necessarily imply the presence of any specific gas.

Under $\Lambda$CDM, recent baryon acoustic oscillation (BAO) distance measures from DESI, which favor a low matter density $\Omega_m$, are in moderate $2-3\sigma$ tension with cosmic microwave background (CMB) observations. This tension appears alternately as a preference for the sum of neutrino masses dropping below the $\sum m_\nu = 0.06$eV value required by neutrino oscillation measurements to formally negative values; a discrepant value of $\Omega_m$ at 0.06eV; or preference for dynamical dark energy beyond $\Lambda$CDM. We show that this tension largely arises from the CMB lensing constraints on the calibration of the sound horizon for geometric measurements and relies on the measurement of the reionization optical depth $\tau$ from large-angle CMB polarization to set the lensing amplitude. Dropping these constraints removes the neutrino tension at $\sum m_\nu=0.06$eV entirely, favoring $\tau = 0.091\pm 0.011$ in $\Lambda$CDM. Beyond $\Lambda$CDM, it brings the preference for $w_0-w_a$ dynamical dark energy to below $95\%$ CL. We explore the freedom in interpreting the low-$\ell$ EE polarization constraint due to analysis choices and reionization modeling beyond the standard step-function assumption and find that this drops the neutrino tension in $\Lambda$CDM to below $95\%$ CL. Alternately, this raising of $\tau$ can also be achieved by the same reduction in large-scale curvature fluctuations that also ameliorates the low-$\ell$ temperature anomaly.

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value $f_\nu$ for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor $(1-f_\nu)$. Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of $f_\nu$ over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is however inadequate to estimate $f_\nu$: in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of $f_\nu$ to jump from $f_\nu=0$ to $f_\nu\sim 10^5$. While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential as the interactions used in simulations avoid direct calculations of $f_\nu$. We are however reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of $f_\nu$ can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

L. Monaco, E. Caffau, P. Molaro, P. Bonifacio, G. Cescutti

Aims. The metal-poor star BPM 3066 belongs to the retrograde halo and presents unexpectedly strong spectral features of lithium. To gain insight into the origin of this peculiar abundance, we investigate the chemistry and kinematic properties of this star. Methods. We performed a local thermodynamic equilibrium chemical abundance analysis of UVES/VLT high-resolution spectra of BPM 3066 using ATLAS9 and ATLAS12 model atmospheres and the MyGIsFOS code. We further characterised the orbital properties of the star by integrating its orbit and analysing its integrals of motion using the galpy code. Results. The star BPM 3066 shows an exceptional overabundance of both lithium and beryllium. The abundances are A(Li) = 3.0 and A(Be) = 2.1, which are respectively about 0.8 and 2.2 dex higher than the Li and Be abundances expected at [Fe/H] = -1.5, the metallicity of the star. The observed ratio 7Li/9Be is 7.9, which is close to that expected from a synthesis by spallation processes. Overabundances of Si, Al, and of the neutron capture elements Sr,Y, Zr, and Ba are also measured. Kinematically, BPM 3066 has an eccentric, strongly retrograde orbit, confined to a height lower than 1 kpc from the galactic plane, and it is a candidate member of the Sequoia/Thamnos accreted galaxy. Conclusions. The processes leading to the 7Li and 9Be synthesis could have occurred in the environment of a hypernova. This is supported by some abundance anomalies like the high value of Si, [Si/Fe]=1.2 and [Si/O]=1.1. However, the simultaneous high values of N, Na, Al, Sc, Ti, and Cu are at odds with the expectations from a hypernova. Alternatively, the abundances of BPM 3066 could result from the engulfing of rocky planets that were rich in spallated Li and Be. In both cases, it is remarkable that such an extreme abundance pattern has been found in a star belonging to the Sequoia/Thamnos accreted galaxy.

Franc O, Pavlos Protopapas, Dominic W. Pesce, Angelo Ricarte, Sheperd S. Doeleman, Cecilia Garraffo, Lindy Blackburn, Mauricio Santillana

Using very long baseline interferometry, the Event Horizon Telescope (EHT) collaboration has resolved the shadows of two supermassive black holes. Model comparison is traditionally performed in image space, where imaging algorithms introduce uncertainties in the recovered structure. Here, we develop a deep learning framework to perform parameter inference in visibility space, directly using the data measured by the interferometer without introducing potential errors and biases from image reconstruction. First, we train and validate our framework on synthetic data derived from general relativistic magnetohydrodynamics (GRMHD) simulations that vary in magnetic field state, spin, and $R_\mathrm{high}$. Applying these models to the real data obtained during the 2017 EHT campaign, and only considering total intensity, we do not derive meaningful constraints on either of these parameters. At present, our method is limited both by theoretical uncertainties in the GRMHD simulations and variation between snapshots of the same underlying physical model. However, we demonstrate that spin and $R_\mathrm{high}$ could be recovered using this framework through continuous monitoring of our sources, which mitigates variations due to turbulence. In future work, we anticipate that including spectral or polarimetric information will greatly improve the performance of this framework.

H. A. Adarsha (MCNS, India), Chandrachur Chakraborty (MCNS, India), Sudip Bhattacharyya (TIFR, India)

Sub-solar mass black holes could show up in gravitational observations in future and near-solar mass black holes might have been involved in the events GW190425 and GW190814. Since they cannot form from the stellar evolution, their creation requires exotic mechanisms. One such mechanism involves the capture of dark matter particles by stellar objects and their thermalization. When the criterion for the collapse of these dark matter particles is satisfied, a tiny endoparasitic black hole (EBH) forms and then it accretes matter from the host. The EBH may transmute the host into a black hole of nearly the same mass as the host or lesser, depending on the type of accretion. We examine this complex and poorly explored accretion mechanism, considering the effects of rotation and viscosity but ignoring some other effects, such as those of pressure and magnetic field, as the first step. Using a general framework to assess the effects of rotation and viscosity on accretion, we show that the accretion could be stalled in some white dwarfs, but not in neutron stars. The stalled accretion should cause an opening in the host's polar regions, the extent of which depends on the mass and spin of the host.

We present StarryStarryProcess, a novel hierarchical Bayesian framework for mapping stellar surfaces using exoplanet transit light curves. While previous methods relied solely on stellar rotational light curves--which contain limited information about spot properties -- our approach leverages planetary transits as probes of stellar surfaces. When a planet crosses a spot during transit, it creates a distinctive change in the light curve that directly reveals spot properties. Our model integrates planetary transit modeling with stellar variability analysis by combining the spherical harmonic surface map representation from starry, the probabilistic approach to spot properties of StarryProcess, and a comprehensive transit model that accounts for spot-crossing events during transits. We demonstrate through synthetic data experiments that our model successfully recovers spot distributions, stellar orientation, and spot physical properties. We extend the framework to handle evolving stellar surfaces through time-dependent modeling. Applying our method to TESS observations of TOI-3884, we find evidence for high-latitude spot concentrations and significant spin-orbit misalignment. The transit-based approach overcomes fundamental limitations of previous models by providing constraints on spot properties that would remain hidden in the null space of rotational light curves alone. This methodology enables more accurate exoplanet characterization by disentangling stellar activity due to starspots from planetary signals while simultaneously providing insights into stellar magnetic activity patterns. The whole paper is reproducible, and can be found by clicking the GitHub icon.

E. Petrova, T. Van Doorsselaere, M. van Noort, D. Berghmans, J. S. Castellanos Durán

Context. Fine-scale structures of the solar chromosphere, particularly fibrils, are known to host various types of magnetohydrodynamic (MHD) waves that can transport energy to the corona. In particular, absorption features observed in the H{\alpha} channel have been widely detected that exhibit transverse oscillations. Aims. We aimed to detect a high-frequency transverse oscillation in fibrils. Methods. We conducted a case study on a high-frequency transverse oscillation in a chromospheric fibril. A chromospheric fibril was observed on 24 August 2018, in the H{\alpha} spectral line, with the prototype Microlensed Hyperspectral Imager (MiHI) at the Swedish 1- meter Solar Telescope. The MiHI instrument is an integral field spectrograph capable of achieving ultra-high resolution simultaneously in the spatial, temporal, and spectral domains. Results. The detected oscillation characteristics include a period of 15 s and a displacement amplitude of 42 km. Using the bisector method, we derived Doppler velocities and determined that the polarisation of the oscillation was elliptical. Conclusions. The energy contained in the oscillation ranges from 390 to 2300 W/m2, which is not sufficient to balance radiative losses of the chromosphere.

Motivated by the recent discovery of situations where cosmological fluctuations recohere during inflation, we investigate the relationship between quantum recoherence (late-time purification after a transient phase of decoherence), adiabaticity, and Markovianity. To that end, we study a simple setup of two linearly-coupled harmonic oscillators, and compute the purity of one oscillator when the interaction is switched off. We find that there exists a critical value for the coupling strength below which the purity oscillates and above which it decays exponentially. This decay cannot be captured with perturbation theory, hence, decoherence is always a non-perturbative phenomenon. When the interaction is turned off, the purity either freezes to its value prior to the turn-off, or it smoothly goes back to a value very close to one (recoherence). This depends on the rate at which the turn-off occurs. We thus develop a new adiabatic-expansion scheme and find complete recoherence at any finite order in the inverse turn-off time. Therefore, decoherence is always a non-adiabatic effect. The critical value of the turn-off time above which recoherence takes place is then expressed in terms of the other time scales of the problem. Finally, we show that even though the dynamics of the system is never Markovian, when decoherence takes place, it can always be interpreted as resulting from a locally Markovian process. We introduce a new measure of Markovianity dubbed the Bures velocity and use it to optimise Markovian approximations.

We suggest that a gauge theory CP-violating phase could be degenerate with a magnetic dual of a $4$-form flux. Discrete discharge of this flux by membrane nucleations could reduce the total CP-violating phase to below $10^{-10}$ well before BBN if the charge and the tension of the membranes are in the $\sim keV$ range.

Using the Kadanoff--Baym formalism, we perform a detailed study of neutrino-antineutrino synchrotron emission from strongly magnetized, dense quark matter under conditions relevant to compact stars. Starting from an exact expression for the emission rate that fully accounts for Landau-level quantization of quarks, we derive an approximate formula applicable in the regime where quark chemical potentials are much larger than all other relevant energy scales. We demonstrate that the emission rate is largely controlled by a single dimensionless ratio between two low-energy scales: the Landau-level spacing at the Fermi surface, $|e_f B|/\mu_{f}$, and the temperature of the quark matter, $T$. When the ratio $|e_f B|/(\mu_{f} T)$ is small, many closely spaced Landau levels contribute to the emission. In the opposite limit, when the ratio is large, the rate is dominated by transitions between adjacent levels and is exponentially suppressed. Our results show that, even in the presence of the strongest magnetic fields expected in compact stars, the synchrotron emission remains suppressed by more than three orders of magnitude compared to the direct Urca process. This implies that such emission is unlikely to play any substantial role in the cooling of magnetized quark stars, at least those made of unpaired quark matter phases.

Francisco R. Candón, Pablo Casaseca, Maurizio Giannotti, Mathieu Kaltschmidt, Jaime Ruz, Julia K. Vogel

A finite axion-nucleon coupling enables the production of axions in stellar environments via the thermal excitation and subsequent de-excitation of the $^{57}$Fe isotope. Given its low-lying excited state at 14.4 keV, $^{57}$Fe can be efficiently excited in the hot cores of supergiant stars, possibly leading to axions emission. The conversion of these axions into photons in the Galactic magnetic field results in a characteristic 14.4 keV line, potentially detectable by hard X-ray telescopes such as NASA's Nuclear Spectroscopic Telescope Array (NuSTAR). In this work, we present the first constraints on axion-nucleon couplings derived from \textsc{NuSTAR} observations of Betelgeuse and discuss the potential insights that could be gained from detecting this line in other nearby supergiants. Our results establish significantly more stringent bounds than those obtained from solar observations, setting a limit of $|g_{a\gamma} g_{aN}^{\mathrm{eff}}| < (1.2 - 2.7) \times 10^{-20}$ GeV$^{-1}$ for $m_a \lesssim 10^{-10}$ eV.

In this paper, we introduce a new interacting mechanism within the dark sector, encompassing both dark energy and dark matter, while grounding our analysis in the familiar framework of the $\mathrm{\Lambda CDM}$ model augmented by baryons and radiation components, including photons and neutrinos. The interaction between dark energy and dark matter is confined to the perturbative level. One significant advantage of this proposal is that all geometric probes yield constraints consistent with the so-called vanilla model or the extended vanilla model, where dark energy has a constant equation of state, $w_{x}$. However, the introduction of this new interacting mechanism affects several theoretical signatures, involving contrast dark matter and dark energy densities. We perform an exploratory analysis of those effects in the CMB power spectra, matter-power spectra, and redshift space distortions. For instance, it allows for a decrease/increment in the integrated Sachs-Wolfe (ISW) effect depending on the value taken by the interaction coupling. This effect could be observationally detected by looking for a cross-correlation between the ISW temperature fluctuations and the distribution of galaxies or quasars. At late times, the interaction in the dark sector becomes very effective, affecting the non-linear scale of structure formation. We discuss how the estimators $f\sigma_{8}(z)$ and $S_{8}(z)$ are affected by different interacting couplings; indicating that $f\sigma_{8}(z)$ can show a relative change of up to $15\%$ compared to the concordance model at low redshifts. Finally, we show how the various terms in the dark energy pressure perturbation (both adiabatic and non-adiabatic) are relevant for different scales, demonstrating the absence of large-scale instabilities.

This paper develops a novel hybrid approach for estimating the mixture model of $t$-factor analyzers (MtFA) that employs multivariate $t$-distribution and factor model to cluster and characterize grouped data. The traditional estimation method for MtFA faces computational challenges, particularly in high-dimensional settings, where the eigendecomposition of large covariance matrices and the iterative nature of Expectation-Maximization (EM) algorithms lead to scalability issues. We propose a computational scheme that integrates a profile likelihood method into the EM framework to efficiently obtain the model parameter estimates. The effectiveness of our approach is demonstrated through simulations showcasing its superior computational efficiency compared to the existing method, while preserving clustering accuracy and resilience against outliers. Our method is applied to cluster the Gamma-ray bursts, reinforcing several claims in the literature that Gamma-ray bursts have heterogeneous subpopulations and providing characterizations of the estimated groups.

When testing general relativity (GR) with gravitational wave observations, parametrized tests of deviations from the expected strong-field source dynamics are one of the most widely used techniques. We present an updated version of the parametrized framework with the state-of-art IMRPhenomX waveform family. Our new framework incorporates deviations in the dominant mode as well as in the higher-order modes of the waveform. We demonstrate that the missing physics of either higher-order modes or precession in the parametrized model can lead to a biased conclusion of false deviation from GR. Our new implementation mitigates this issue and enables us to perform the tests for highly asymmetric and precessing binaries without being subject to systematic biases due to missing physics. Finally, we apply the improved test to analyze events observed during the second half of the third observing run of LIGO and Virgo (O3b). We provide constraints on GR deviations by combining O3b results with those from previous observation runs. Our findings show no evidence for violations of GR.

The dynamics of the electroweak phase transition in the early universe has profound implications for cosmology and particle physics. We systematically study the steady-state dynamics of bubble walls in scenarios where the transition is first order within three representative beyond the Standard Model frameworks, characterised by the presence of an additional scalar in different electroweak representations. Focusing on the local thermal equilibrium regime, we numerically solve the coupled scalar and hydrodynamic equations to extract key properties of the phase transition front: the wall velocity, width, plasma and field profiles. Remarkably, we find a near-universal behaviour across models when expressed in terms of thermodynamic quantities, that can be captured by simple fitting functions, useful for phenomenological applications. These results also provide an upper bound on the bubble velocity and represent the first necessary step for the full inclusion of out-of-equilibrium effects.

Since the recent announcements of evidence for a stochastic gravitational wave background from several pulsar timing array collaborations, much effort has been devoted to explore features beyond the fiducial Hellings-Downs background including those arising in modified gravity theories and deterministic gravitational wave signals. Inspired by previous studies, we propose a method to efficiently screen these models using likelihood reweighting based on the fiducial model. In order to alleviate the well-known unstable weight estimates in vanilla importance sampling, we implement reweighting for the second time making use of the kernel density estimation of the previously reweighted samples. We tested this method by analyzing three simulated datasets with an injected sinusoid signal applied to all pulsars. It is found that likelihood reweighting not only gives results compatible with those from full Bayesian analyses when the signal is subdominant, but is also able to recover the signal posterior to a reasonable accuracy in the presence of a rather strong signal. Given samples from the fiducial model, this method could bring an at least $\mathcal{O}(10)$-time speedup in analyzing new models.

Observations of microlensed gravitational waves (GWs) emanated by compact binary coalescences (CBCs) are essential for studying the mass density distribution in the universe, including black holes and dark matter halos. However, no confident detection of microlensed GWs have been reported to date. There are two important challenges in the identification of microlensed GWs. The first is that the source waveform and lens structure models are not known a-priori. The second is that certain classes of unlensed GWs could mimic microlensed GWs, resulting in undesirable false alarms. In this work, we propose to use the Kramers-Kronig relation for gravitational lensing systems. We argue that such systems are essentially linear response systems obeying causality, where KK relation must hold. The power of this method lies in the fact that microlensed GWs, regardless of the lens structure, must obey KK relation, while unlensed GW events are not in general expected to obey it. This, in principle, allows us to identify microlensed GWs while dismissing microlensing mimickers. We provide the first important steps towards a methodology that exploits KK relation, and test its usefulness under idealized conditions.

NGC 1068 is the brightest extragalactic source in high-energy neutrinos as seen by IceCube, yet the accompanying gamma-ray flux is orders of magnitude weaker. It has been argued that this indicates that the bulk of neutrinos and gamma rays are emitted in the innermost vicinity of the central supermassive black hole, which is transparent to neutrinos, but opaque to gamma rays. Even in such extreme scenarios for the acceleration of cosmic rays, astrophysical models typically overestimate the low-energy gamma-ray flux and/or require some fine-tuning in the physical parameters. Here we suggest instead that the dark matter surrounding the supermassive black hole may absorb the gamma rays, inducing the observed deficit. We show that for a dark matter-photon scattering cross section in the range $\sigma_{\rm DM-\gamma}/m_{\rm DM} \simeq 10^{-28}-10^{-30}$ cm$^2$/GeV, Fermi-LAT measurements can be well reconciled with IceCube data. We also present some simple particle physics examples that achieve the correct spectral energy dependence while respecting complementary constraints.

Diffusion in Yukawa crystals is stochastic due to the thermally activated formation of vacancy-interstitial pairs, which have poor statistics in simulations. This makes it difficult to argue if Yukawa crystals exhibit normal diffusion, or if they could be subdiffusive or superdiffusive. To resolve this, we run a long molecular dynamics simulation of an idealized Yukawa crystal for a billion timesteps. We find no evidence of anomalous diffusion in the pure crystal, but also caution readers against overinterpreting this result as real crystals have complicated structures including grains and defects.

In numerical simulations of core-collapse supernova and binary neutron stars mergers information about the energetics and composition of matter is implemented via external tables covering the huge ranges of thermodynamic conditions explored during the astrophysical evolution. More than 120 general purpose equation of state tables have been contributed so far. Not all of them comply with current constraints from theoretical and experimental nuclear physics and astrophysical observations of neutron stars. Systematic investigations of the role that dense matter properties play in the evolution of these astrophysical phenomena require that more equation of state tables are provided. We build a set of general purpose equation of state tables. At zero temperature, they comply with all currently accepted constraints, including ab initio chiral effective field theory calculations of pure neutron and symmetric nuclear matter. This set is designed to explore a wide variety of the behaviors of the effective masses as functions of density, which is reflected into a wide range of thermal behaviors. We employ Brussels extended Skyrme interactions generated by means of Bayesian inference techniques. An extended nuclear statistical equilibrium model is developed for modeling sub-saturated inhomogeneous nuclear matter. We study the properties of sub-saturated inhomogeneous nuclear matter over wide ranges of density, temperature and proton fraction. We analyze the mechanisms of transition to homogeneous matter and estimate the transition density. Our key results include the presence of a think $^{14}$He layer in the inner crusts of (neo-)neutron stars, significant abundance of other exotic isotopes of H and He in warm and neutron rich matter and a detailed study of the thermodynamic stability of cold stellar matter. The equation of state tables will be publicly available in the Compose online database.