Papers for Wednesday, Jun 25 2025

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period vs the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher-order numerical solver with adaptive time steps. Specifically, we aim to (i) reproduce a pulsar's high-energy emission maps, namely one with $10\%$ the surface $B$-field strength of Vela, and spectra produced by an independent gyro-centric pulsar emission model; (ii) test convergence of these results to the radiation-reaction limit of Aristotelian Electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel}$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using $10\%$ field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel}$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf{E}\times \mathbf{B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators' model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

As the field of exoplanetary astronomy has matured, there has been growing demand for precise stellar abundances to probe subtle correlations between stellar compositions and planetary demographics. However, drawing population-level conclusions from the disparate measurements in the literature is challenging, with various groups measuring metallicities using bespoke codes with differing line lists, radiative transfer models, and other assumptions. Here we use the neural-net framework uberMS to measure iron abundances and alpha enchriments from high-resolution optical spectra observed by the Tillinghast Reflector Echelle Spectrograph (TRES), a key resource used for the follow-up of candidate exoplanet hosts. To contextualize these measurements and benchmark the performance of our method, we compare to external constraints on metallicity using the Hyades cluster, wide binaries, and asteroids, to external constraints on $T_{\rm eff}$ and $\log g$ using stars with interferometric radii, and to the results of other abundance measurement methods using overlap samples with the APOGEE and SPOCS catalogs, as well as by applying the SPC method directly to TRES spectra. We find that TRES-uberMS provides parameter estimates with errors of roughly 100K in $T_{\rm eff}$, 0.09dex in $\log g$, and 0.04dex in [Fe/H] for nearby dwarf stars. However, [Fe/H] performance is significantly poorer for mid-to-late K dwarfs, with the bias worsening with decreasing $T_{\rm eff}$. Performance is also worse for evolved stars. Our [$\alpha$/Fe] error may be as good as 0.03dex for dwarfs based on external benchmarks, although there are large systematic differences when comparing with specific alpha-element abundances from other catalogs.

Disentangling the rich astrophysical structure of the stochastic gravitational-wave background (SGWB) from its cosmological component is essential for the Laser Interferometer Space Antenna (LISA) to access the physics of the early Universe beyond the reach of any other probe. In this work, we develop an analytical framework to compute the astrophysical contribution to the SGWB arising from an unresolved ensemble of inspiraling and merging black hole binaries. Accounting for various resolvability thresholds, we leverage this framework to predict the amplitude, spectral shape, and detection signal-to-noise ratio of the unresolved background from massive black hole binaries (MBHBs), capturing the possible diversity of its SGWB imprint across a range of astrophysically motivated populations. Through a joint analysis of astrophysical and primordial contributions to the SGWB, we determine the minimum detectable amplitude of the cosmological background across a range of spectral shapes. We demonstrate that, even under optimistic subtraction thresholds, unresolved MBHBs can degrade the detectability of a cosmological signal by multiple orders of magnitude, depending on the spectral shape of the primordial component. Ultimately, the MBHB-induced astrophysical SGWB acts both as a veil and a lens: it imposes a fundamental limit on cosmological sensitivity, yet simultaneously reveals the hidden population of massive black holes beyond the reach of individual detections. Accurate modeling of this late-Universe background is therefore a prerequisite for robust component separation and for realizing LISA's full scientific potential.

Radio broadband spectro-polarimetric observations are sensitive to the spatial fluctuations of the Faraday depth (FD) within the telescope beam. Such FD fluctuations are referred to as "Faraday complexity", and can unveil small-scale magneto-ionic structures in both the synchrotron-emitting and the foreground volumes. We explore the astrophysical origin of the Faraday complexity exhibited by 191 polarised extragalactic radio sources (EGSs) within 5 deg from the Galactic plane in the longitude range of 20-52 deg, using broadband data from the Karl G. Jansky Very Large Array presented by a previous work. A new parameter called the FD spread is devised to quantify the spatial FD fluctuations. We find that the FD spread of the EGSs (i) demonstrates an enhancement near the Galactic mid-plane, most notable within Galactic latitude of +-3 deg, (ii) exhibits hints of modulations across Galactic longitude, (iii) does not vary with the source size across the entire range of 2.5"-300", and (iv) has an amplitude higher than expected from magneto-ionic structures of extragalactic origin. All these suggest that the primary cause of the Faraday complexity exhibited by our target EGSs is <2.5"-scale magneto-ionic structures in the Milky Way. We argue that the anisotropic turbulent magnetic field generated by galactic-scale shocks and shears, or the stellar feedback-driven isotropic turbulent magnetic field, are the most likely candidates. Our work highlights the use of broadband radio polarimetric observations of EGSs as a powerful probe of multi-scale magnetic structures in the Milky Way.

We present new analysis of the CWW 89 system as part of the Orbital Architectures of Transiting Massive Exoplanets And Low-mass stars (OATMEAL) survey. The CWW 89 system is a member of the 2.8 Gyr old Ruprecht 147 (NGC 6774) cluster and features two stars, CWW 89A (EPIC 219388192) and CWW 89B, with the primary hosting a transiting brown dwarf. We use in-transit, highly precise radial velocity measurements with the Keck Planet Finder (KPF) to characterize the Rossiter-McLaughlin (RM) effect and measure the projected spin-orbit obliquity $|\lambda|=1.4\pm2.5^\circ$ and the full 3D spin-orbit obliquity of the brown dwarf to be $\psi=15.1^{+15.0^\circ}_{-10.9}$. This value of $\lambda$ implies that the brown dwarf's orbit is prograde and well-aligned with the equator of the host star, continuing the trend of transiting brown dwarfs showing a preference for alignment ($\lambda \approx 0^\circ$) regardless of the stellar effective temperature. We find that this contrast with the transiting giant planet population, whose spin-orbit alignments depend on host $T_{\rm eff}$, shows an increasingly clear distinction in the formation and orbital migration mechanisms between transiting giant planets and transiting brown dwarfs like CWW 89Ab. For this system in particular, we find it plausible that the brown dwarf may have undergone coplanar high-eccentricity migration influence by CWW 89B.

Active, low-mass stars are widely observed to have radii that are larger than predicted by standard stellar models. Proposed mechanisms for this radius inflation generally involve stellar magnetism, either in the form of added pressure support in the outer layers and/or suppression of convection via starspots. We have assembled a large sample of 261 low-mass stars in the young clusters Upper Scorpius, $\alpha$ Persei, Pleiades, and Praesepe (spanning ages 10--670 Myr) for which the data exist to empirically measure the stellar radii, rotation periods, and starspot covering fractions. We find a clear, strong relationship between the degree of radius inflation and stellar rotation as represented by the Rossby number; this inflation-rotation relationship bears striking resemblance to canonical activity-rotation relationships, including both the so-called linear and saturated regimes. We also demonstrate here for the first time that the radius inflation depends directly on the starspot covering fraction. We furthermore find that the stars' effective temperatures decrease with decreasing Rossby number as well, and that this temperature suppression balances the radius inflation so as to preserve the stellar bolometric luminosity. These relationships are consistent across the age range sampled here, which spans from the pre--main-sequence to the zero-age main sequence. The favorable comparison of our findings to the predictions of modern starspot-based stellar evolution models suggests that, while rotation is clearly the underlying driver, magnetism may be the most likely direct cause of the radius inflation phenomenon.

Nested sampling (NS) is the preferred stochastic sampling algorithm for gravitational-wave inference for compact binary coalenscences (CBCs). It can handle the complex nature of the gravitational-wave likelihood surface and provides an estimate of the Bayesian model evidence. However, there is another class of algorithms that meets the same requirements but has not been used for gravitational-wave analyses: Sequential Monte Carlo (SMC), an extension of importance sampling that maps samples from an initial density to a target density via a series of intermediate densities. In this work, we validate a type of SMC algorithm, called persistent sampling (PS), for gravitational-wave inference. We consider a range of different scenarios including binary black holes (BBHs) and binary neutron stars (BNSs) and real and simulated data and show that PS produces results that are consistent with NS whilst being, on average, 2 times more efficient and 2.74 times faster. This demonstrates that PS is a viable alternative to NS that should be considered for future gravitational-wave analyses.

Understanding the internal structure of core helium burning (CHeB) stars is crucial for evaluating transport processes in nuclear-burning regions, constructing accurate stellar population models, and assessing nucleosynthesis processes that impact the chemical evolution of galaxies. While asteroseismic observations have recently enabled detailed probing of CHeB star interiors, seismic signatures related to structural variations at the boundary between the convective and radiative core, and chemical composition gradients within the radiative core remain underexplored. This paper investigates how such gradients affect the oscillation modes of low-mass CHeB stars, focusing on mixed dipole modes and uncoupled g-modes as diagnostic tools. Using semi-analytical models calibrated with the evolutionary codes $\texttt{BaSTI-IAC}$,$ \texttt{CLES}$, and $\texttt{MESA}$, we examine the impact of density discontinuities and associated structural glitches on mode period spacings. These codes span diverse physical prescriptions, allowing us to isolate robust features relevant for calibration. Our approach enables controlled glitch insertion while preserving a realistic representation of the star. Consistent with prior works, we find that structural glitches introduce periodic components in the period spacings, providing constraints on the location and amplitude of interior variations. We compare models with smooth and sharp transitions, demonstrating how glitch sharpness affects period spacing and mode trapping. Simulations based on four-year $\textit{Kepler}$ data show that our models yield oscillation frequencies closely matching observations. Ultimately, our results offer realistic predictions of how specific structural features affect the power spectral density, validating our theoretical framework and guiding future efforts to interpret glitch signatures in high-precision asteroseismic data.

The presence of magnetic fields in galaxies and their haloes could have important consequences for satellite galaxies moving through the magnetised circumgalactic medium (CGM) of their host. We therefore study the effect of magnetic fields on ram pressure stripping of satellites in the CGM of massive galaxies. We use cosmological `zoom-in' simulations of three massive galaxy haloes ($M_{\rm{200c}} = 10^{12.5-13}$M$_\odot$), each simulated with and without magnetic fields. The fraction of gas retained after infall through the CGM across our full satellite sample shows no population-wide impact of magnetic fields on ram pressure stripping. However, for the most massive satellites, we find that twice as much gas is stripped without magnetic fields. The evolution of a galaxy's stripped tail is also strongly affected. Magnetic fields reduce turbulent mixing, significantly inhibiting the dispersion of metals into the host CGM. This suppressed mixing greatly reduces condensation from the host CGM onto the stripped tail. By studying the magnetic field structure, we find evidence of magnetic draping and attribute differences in the stripping rate to the draping layer. Differences in condensation from the host CGM are attributed to magnetic field lines aligned with the tail suppressing turbulent mixing. We simulate one halo with enhanced resolution in the CGM and show these results are converged with resolution, though the structure of the cool gas in the tail is not. Our results show that magnetic fields can play an important role in ram pressure stripping in galaxy haloes and should be included in simulations of galaxy formation.

The Solar and Heliospheric Observatory (SOHO) Extreme-ultraviolet Imaging Telescope (EIT) has been taking images of the Solar disk and corona in four narrow EUV bandpasses (171Å, 195Å, 284Å, and 304Å) at a minimum cadence of once per day since early 1996. The time series of fully-calibrated EIT images now spans approximately 28 years, from early 1996 to early 2024, covering solar cycles 23, 24, and the beginning of cycle 25. We convert this extensive EIT image archive into a time series of `Sun-as-a-star' light curves in EIT's four bandpasses, providing a long-term record of solar EUV variability. These Sun-as-a-star light curves, available for download from this https URL, trace the Sun as if it were a distant point source, viewed from a fixed perspective. We find that our EUV light curves trace the $\sim$ 11-year solar activity cycle and the $\sim$ 27-day rotation period much better than comparable optical observations. In particular, we can accurately recover the solar rotation period from our 284Ålight curve for 26 out of 28 calendar years of EIT observations (93% of the time), compared to only 3 out of 29 calendar years (10% of the time) of the VIRGO total solar irradiance time series, which is dominated by optical light. Our EIT light curves, in conjunction with Sun-as-a-star light curves at optical wavelengths, will be valuable to those interested in inferring the EUV/UV character of stars with long optical light curves but no intensive UV observations, as well as to those interested in long-term records of solar and space weather.

Direct Collapse scenario to form the supermassive black hole (SMBH) seeds offers the most promising way to explain the origin of quasars at $z>7$. Assuming atomic primordial gas, can Ly$\alpha$ photons escape from the central regions of the collapse and serve as a diagnostic for the detection of these pre-SMBH objects? Models assuming spherical collapse have found these photons to be trapped and destroyed. We use Ly$\alpha$ radiation transfer within the inflow-outflow geometry, based on earlier zoom-in cosmological modeling involving radiation transfer and magnetic forces. Adopting geometry that includes ongoing disky and spherical accretion, and formation of bi-conical outflow funnel, we obtain formation of a dense radiatively-driven expanding shell. The Ly$\alpha$ transfer is performed through a Monte Carlo algorithm, accounting for destruction of Ly$\alpha$ photons and emergence of two-photon emission. We find that substantial fraction of Ly$\alpha$ photons can escape through the funnel, and calculate the line profiles, the line peak velocity shift, asymmetry and cuspiness, by varying basic model parameters. The escaping Ly$\alpha$ emission is anisotropic, and sensitive to overall inflow-outflow geometry. The escaping fraction of Ly$\alpha$ radiation exceeds 95% from a $z=10$ pre-SMBH object -- in principle detectable by the JWST NIRSpec in the MOS mode, during $\sim 10^4$seconds for a $10\sigma$ signal-to-noise ratio. Moreover, comparisons with line shapes from high-$z$ galaxies and quasars allow to separate them from pre-SMBH objects based on the line asymmetry -- the pre-SMBH line show a profound asymmetry and extended red tail.

The new generation of optical spectrographs (i.e., WEAVE, 4MOST, DESI, and WST) offer unprecedented opportunities for statistically studying the star formation histories of galaxies. However, these observations are not easily comparable to predictions from cosmological simulations. Our goal is to build a reference framework for comparing spectroscopic observations with simulations and test tools for deriving stellar population properties of galaxies. We focus on the observational strategy of the Stellar Population at Intermediate Redshift Survey (StePS) with the WEAVE instrument. We generate mock datasets of ~750 galaxies at redshifts z = 0.3, 0.5, and 0.7 using the TNG50 simulation, perform radiative transfer with SKIRT, and analyze the spectra with pPXF as if they were real observations. We present the methodology to generate these datasets and provide an initial exploration of stellar population parameters (i.e., mass-weighted ages and metallicities) and star formation histories for three galaxies at z = 0.7 and their descendants at z = 0.5 and 0.3. We find good agreement between the mock spectra and intrinsic ages in TNG50 (average difference $0.2\pm0.3$ Gyr) and successfully recover their star formation histories, especially for galaxies form the bulk of their stars on short timescales and at early epochs. We release these datasets, including multi-wavelength imaging and spectra, to support forthcoming WEAVE observations.

Some galaxies such as the Milky Way and Andromeda display coherently rotating satellite planes, posing tensions with cosmological simulations. NGC 2750 has emerged as an additional candidate system hosting a co-rotating group of galaxies. We aim to assess the presence of a coherent satellite plane around NGC 2750 by identifying new candidate dwarf galaxies and low surface brightness features. We conducted deep, wide-field photometric observations of NGC 2750 using the Large Binocular Telescope in the g- and r-bands. Standard data reduction techniques were applied to enhance the detection of low-surface-brightness features down to about 31 mag/arcsec^2 in r. Our observations led to the discovery of six new candidate dwarf galaxies, including one with properties consistent with an ultra-diffuse galaxy. We also identified tidal features around NGC 2750, indicating past interactions with nearby satellites. The spatial distribution of satellites suggests a moderate flattening, further supported by the newly identified candidates. Follow-up spectroscopic measurements will be critical in confirming or challenging the strong kinematic coherence observed previously. The luminosity function of NGC 2750 reveals an excess of bright satellites compared to similar systems, adding to the growing tension between observed satellite populations and cosmological simulations.

Precise localizations of a small number of repeating fast radio bursts (FRBs) using very long baseline interferometry (VLBI) have enabled multiwavelength follow-up observations revealing diverse local environments. However, the 2--3\% of FRB sources that are observed to repeat may not be representative of the full population. Here we use the VLBI capabilities of the full CHIME Outriggers array for the first time to localize a nearby (40 Mpc), bright (kJy), and apparently one-off FRB source, FRB 20250316A, to its environment on 13-pc scales. We use optical and radio observations to place deep constraints on associated transient emission and the properties of its local environment. We place a $5\sigma$ upper limit of $L_{\mathrm{9.9~\mathrm{GHz}}} < 2.1\times10^{25}~\mathrm{erg~s^{-1}~Hz^{-1}}$ on spatially coincident radio emission, a factor of 100 lower than any known compact persistent radio source associated with an FRB. Our KCWI observations allow us to characterize the gas density, metallicity, nature of gas ionization, dust extinction and star-formation rate through emission line fluxes. We leverage the exceptional brightness and proximity of this source to place deep constraints on the repetition of FRB 20250316A, and find it is inconsistent with all well-studied repeaters given the non-detection of bursts at lower spectral energies. We explore the implications of a measured offset of 190$\pm20$ pc from the center of the nearest star-formation region, in the context of progenitor channels. FRB 20250316A marks the beginning of an era of routine localizations for one-off FRBs on tens of mas-scales, enabling large-scale studies of their local environments.

We present deep James Webb Space Telescope near-infrared imaging to search for a quiescent or transient counterpart to FRB 20250316A, which was precisely localized with the CHIME/FRB Outriggers array to an area of $11\times13$ pc in the outer regions of NGC 4141 at $d\approx40$ Mpc. Our F150W2 image reveals a faint source near the center of the FRB localization region ("NIR-1"; $M_{\rm F150W2}\approx-2.5$ mag; probability of chance coincidence $\approx0.36$), the only source within $\approx2.7\sigma$. We find that it is too faint to be a globular cluster, young star cluster, red supergiant star, or a giant star near the tip of the red giant branch (RGB). It is instead consistent with a red giant near the RGB "clump" or a massive ($\gtrsim20$ M$_{\odot}$) main sequence star, although the latter explanation is less likely. The source is too bright to be a supernova remnant, Crab-like pulsar wind nebula, or isolated magnetar. Alternatively, NIR-1 may represent transient emission, namely a dust echo from an energetic outburst associated with the FRB, in which case we would expect it to fade in future observations. We explore the stellar population near the FRB and find that it is composed of a mix of young massive stars ($\sim10-100$ Myr) in a nearby HII region that extends to the location of FRB 20250316A, and old evolved stars ($\gtrsim$ Gyr). The overlap with a young stellar population, containing stars of up to $\approx20$ M$_\odot$, may implicate a neutron star / magnetar produced in the core collapse of a massive star as the source of FRB 20250316A.

Water ice is a fundamental building material of comets and other bodies in the outer solar system. Yet, the properties of cometary water ice are challenging to study, due to its volatility and the typical distances at which comets are observed. Cometary outbursts, impulsive mass-loss events that can liberate large amounts of material, offer opportunities to directly observe and characterize cometary water ice. We present a study of comet 243P/NEAT, instigated by a $-3$ mag outburst that occurred in December 2018. Optical images and a 251-day lightcurve were examined to characterize the outburst and the comet's quiescent activity. Variations in the quiescent lightcurve appear to be dominated by coma asymmetries, rather than changing activity levels as the comet approached and receded from the Sun. Furthermore, the lightcurve shows evidence for 1 to 2 additional small outbursts ($-0.3$ mag) occurring in September 2018. The large December 2018 outburst likely ejected water ice grains, yet no signatures of ice were found in color photometry, a color map, nor a near-infrared spectrum. We discuss possible dynamical and thermal reasons for this non-detection. In this context, we examined the comae of comets 103P/Hartley 2 and C/2013 US$_{10}$ (Catalina), and show that a one-to-one mapping between continuum color and the presence of water ice cannot be supported. We also discuss possible causes for the large outburst, and find that there is an apparent grouping in the kinetic energy per mass estimates for the outbursts of 5 comets.

Combinations of the most recent CMB, BAO, and SNeIa datasets, when analyzed using the CPL parametrization, $w(a) = w_0 + (1 - a) w_a$, exclude $\Lambda$CDM at $\gtrsim\!3\sigma$ in favor of a dark energy equation of state (EoS) parameter that crosses the phantom divide. We confirm this behavior and show that it persists when DESI BAO data are replaced by SH0ES $H_0$ measurements, despite the known tension between these probes in the presence of CMB data. In both cases, the constraints favor a transition from an early-time phantom-like phase to a late-time quintessence-like phase, with the crossing occurring at different redshifts depending on the dataset combination. The probability that a phantom divide line (PDL) crossing does not occur within the expansion history is excluded at significance levels ranging between $3.1\sigma$-$5.2\sigma$. To investigate whether the apparent PDL crossing is a genuine feature preferred by the data or an artifact of the linear form of the CPL parametrization, we isolate the PDL crossing feature by introducing two modified versions of CPL that explicitly forbid it: CPL${}_{>a_\mathrm{c}}$ and CPL${}_{<a_\mathrm{c}}$. These models are physically motivated in that they phenomenologically capture the behavior of thawing and freezing scalar field scenarios. While previous studies have primarily considered thawing quintessence as a non-crossing alternative, we find that a freezing phantom-like model is the only one capable of performing comparably to CPL -- and only in a few cases. Nevertheless, across all dataset combinations considered, the standard CPL model consistently provides the best fit, strongly indicating that the data genuinely favor a PDL crossing.

21 cm intensity mapping (HI IM) can efficiently map large cosmic volumes with good redshift resolution, but systematics and foreground contamination pose major challenges for extracting accurate cosmological information. Cross-correlation with galaxy surveys offers an efficient mitigation strategy, as both datasets have largely uncorrelated systematics. We evaluate the detectability of the 21 cm signal from the BINGO radio telescope by cross-correlating with the LSST photometric survey, given their strong overlap in area and redshift. Using lognormal simulations, we model the cosmological signal in the BINGO frequency range (980 - 1260 MHz), incorporating thermal noise, foregrounds, and cleanup. The LSST simulations include uncertainties in photometric redshift (photo-z) and galaxy number density in the first three redshift intervals (mean redshift approximately equal to 0.25, 0.35, 0.45), corresponding to the expected performance after 10 years of the survey. We show that photo-z errors significantly increase the noise in the cross-correlation, reducing its statistical significance to levels comparable to those of the autocorrelation. Still, the HI signal remains detectable through the cross-correlation, even with photo-z uncertainties similar to those of the LSST. Our results corroborate the feasibility of this approach under realistic conditions and motivate further refinement of the current analysis methods.

Type Ia-CSM supernovae (SNe) are a rare and peculiar subclass of thermonuclear SNe characterized by emission lines of hydrogen or helium, indicative of a high-density circumstellar medium (CSM). Their implied mass-loss rates of $\sim 10^{-4}-10^{-1}$ M$_{\odot}$ yr$^{-1}$ (assuming $\mathrm{ \sim 100 \ km\ s^{-1}}$ winds) from optical observations are generally in excess of values observed in realistic SN Ia progenitors. In this paper, we present an independent study of CSM densities around a sample of 29 archival Ia-CSM SNe using radio observations with the Very Large Array at 6 GHz. Motivated by the late ($\sim$2 yr) radio detection of the Ia-CSM SN 2020eyj, we observed old ($>$1 yr) SNe where we are more likely to see the emergent synchrotron emission that may have been suppressed earlier by free-free absorption by the CSM. We do not detect radio emission down to 3$\sigma$ limits of $\sim$35 $\mu$Jy in our sample. The only radio-detected candidate in our sample, SN 2022esa, was likely mis-classified as a Ia-CSM with early spectra, and appears more consistent with a peculiar Ic based on later-epochs. Assuming a wind-like CSM with temperatures between $2 \times 10^4$ K and $10^5$ K, and magnetic field-to-shock energy fraction ($\epsilon_B$) = $0.01-0.1$, the radio upper limits rule out mass-loss rates between $\sim 10^{-4}-10^{-2}$ M$_{\odot}$ yr$^{-1}$ (100 km s$^{-1}$)$^{-1}$. This is somewhat in tension with the estimates from optical observations, and may indicate that more complex CSM geometries and/or lower values of $\epsilon_B$ may be present.

In this work, we derive systemic velocities and subsequently orbits for 8456 RR~Lyrae stars. We identify interlopers from other Milky Way (MW) structures, which amount to 22 percent of the total sample. Most interlopers are associated with the halo, with the remainder linked to the Galactic disk. We confirm the previously reported lag in the rotation curve of bulge RR~Lyrae stars regardless of the removal of interlopers. Metal-rich RR~Lyrae stars' rotation patterns are consistent with that of non-variable metal-rich giants, following the MW bar, while metal-poor stars exhibit slower rotation. The analysis of orbital parameter space is used to distinguish bulge stars that, in the bar reference frame, have prograde orbits from those in retrograde orbits. We classify the prograde stars into orbital families and estimate the chaoticity (in the form of frequency drift) of their orbits. RR~Lyrae stars with banana-like orbits have a bimodal distance distribution, similar to the distance distribution seen in the metal-rich red clump stars. The fraction of stars with banana-like orbits decreases linearly with metallicity, as does the fraction of stars on prograde orbits (in the bar reference frame). The retrograde moving stars (in the bar reference frame) form a centrally concentrated nearly spherical distribution. Analyzing an $N$-body+SPH simulation, we find that some stellar particles in the central parts oscillate between retrograde and prograde orbits and only a minority stays prograde over a long period of time. Based on the simulation, the ratio between prograde and retrograde stellar particles seems to stabilize within a couple of gigayears after bar formation. The non-chaoticity of retrograde orbits and their high numbers can explain some of the spatial and kinematical features of the MW bulge that have been often associated with a classical bulge.

The new Starlink V2 satellites incorporate improvements to the chassis brightness through dielectric mirrors, off-pointing solar arrays, and black paint on exposed components. For the general case in which the reflectivities are initially unknown, we simulate LSST operations and repeated photometry of every satellite in simulated model constellations. We derive a brightness model of the Starlink V2 satellite and study the simulated apparent brightness as a function of the satellite position relative to the observer and the sun. We find that the V2 Starlink satellites appear brightest at two distinct positions in the sky: when oriented toward the sun at low elevations where light is specularly reflected, and nearly overhead where the satellite is closest to the observer. A simulation of Starlink V2 satellites at 550 km height distributed across a series of Walker constellations with varying inclinations was analyzed to study the impact on the LSST observations. Some bright satellites will be visible in LSST observations. For every thousand V1.5 Starlink satellites imaged by LSST in the first hour of a summer night, we find 1.2 of them will appear brighter than 7 AB magnitude. By comparison, for every thousand V2 Starlink satellites observed, we find only 0.93 of them will appear this bright. The off-pointed solar array and reduced diffuse reflection of the chassis mitigate the brightness. Finally, we simulate lowering this Walker constellation to 350km. Only 0.56 V2 Starlink satellites per thousand brighter than 7 AB magnitude will be observed in the first hour at this height. This is a 40% reduction in number of bright satellites entering the focal plane compared to the constellation at 550km height. We find that a combination of factors yield an apparent surface brightness of these satellites for LSST operations only 5% brighter than at 550km orbit.

The large-angle polarization anisotropies observed in the cosmic microwave background (CMB) are generated by Thomson scattering of CMB photons off free electrons in the post-recombination Universe. In the standard $\Lambda$ cold dark matter cosmological model, the free electron density becomes large at redshifts $z \lesssim 10$ due to the formation of the first stars, leading to reionization of the intergalactic medium. In this work, we use \emph{Gaussian processes} to give a model-independent reconstruction of the reionization history of the Universe constrained by the \textit{Planck} CMB measurements. This method gives a consistent reconstruction of the standard reionization at $z \lesssim 10$, and strongly constrains any additional high-$z$ reionization. We use our results to construct a new derived parameter expressing the high-$z$ contribution to the CMB optical depth, $\tau_{\mathrm{highz}}$. The posterior distribution on $\tau_{\mathrm{highz}}$ can be used to set accurate constraints on models for exotic energy injection, which we demonstrate for the case of decaying dark matter with particle mass in the range $\mathcal{O}(1\,\text{MeV})$. A companion paper demonstrates the use of this methodology to test ensembles of multi-axion models. Our results can be recast to different models using the provided data and code: \href{this https URL}{this http URL\_reio\_gpr}.

The cosmic time dilation observed in Type Ia supernova light curves suggests that the passage of cosmic time varies throughout the evolution of the Universe. This observation implies that the rate of proper time is not constant, as assumed in the standard FLRW metric, but instead is time-dependent. Consequently, the commonly used FLRW metric should be replaced by a more general framework, known as the Conformal Cosmology (CC) metric, to properly account for cosmic time dilation. The CC metric incorporates both spatial expansion and time dilation during cosmic evolution. As a result, it is necessary to distinguish between comoving and proper (physical) time, similar to the distinction made between comoving and proper distances. In addition to successfully explaining cosmic time dilation, the CC metric offers several further advantages: (1) it preserves Lorentz invariance, (2) it maintains the form of Maxwell's equations as in Minkowski space-time, (3) it eliminates the need for dark matter and dark energy in the Friedmann equations, and (4) it successfully predicts the expansion and morphology of spiral galaxies in agreement with observations.

How much of Io's SO$_2$ atmosphere is driven by volcanic outgasing or sublimation of SO$_2$ surface frost is a question with a considerable history. We develop a time dependent surface temperature model including thermal inertia and the exact celestial geometry to model the radiation driven global structure and temporal evolution of Io's atmosphere. We show that many observations can be explained by assuming a purely sublimation driven atmosphere. We find that a thermal diffusivity $\alpha=2.41\times10^{-7}$ m$^2$s$^{-1}$ yields an averaged atmospheric SO$_2$ column density decreasing by more than one order of magnitude from the equator to the poles in accordance with the observed spatial variations of Io's column densities. Our model produces a strong day-night-asymmetry with modeled column density variations of almost two orders of magnitude at the equator as well as a sub-anti-Jovian hemisphere asymmetry, with maximum dayside column densities of $3.7\times10^{16}$ cm$^{-2}$ for the sub-Jovian and $8.5\times10^{16}$ cm$^{-2}$ for the anti-Jovian hemisphere. Both are consistent with the observed temporal and large-scale longitudinal variation of Io's atmosphere. We find that the diurnal variations of the surface temperature affect the subsurface structure up to a depth of 0.6m. Furthermore, we quantify seasonal effects with Io having a northern summer close to perihelion and a northern winter close to aphelion. Finally, we found that at Io's anomalous warm polar regions a conductive heat flux of at least 1.2 Wm$^{-2}$ is necessary to reach surface temperatures consistent with observations.

Changing-look active galactic nuclei, or CL-AGN, are AGN which appear to transition between Seyfert Type 1 and 2 over periods of months to years. Several mechanisms to trigger these transitions have been proposed, but we have yet to conclusively determine their cause. Recent studies suggest CL-AGN are hosted primarily in galaxies which are shutting down star formation (Dodd et al. 2021; Liu et al. 2021, Wang et al. 2023), which may indicate a link between galaxy quenching and changing look events. We use Prospector stellar population synthesis software (Leja et al. 2017; Johnson & Leja 2017; Johnson et al. 2021) to model non-parametric star formation histories for 39 CL-AGN host galaxies. We find that $43^{+13}_{-12}\%$ of our gold sample CL-AGN at z < 0.15 are star forming, while $29^{+13}_{-10}\%$ fall in the Green Valley of the stellar mass-sSFR diagram. At z > 0.15, $57^{+13}_{-18}\%$ of CL-AGN in the gold sample are star-forming and $29^{+19}_{-14}\%$ are in the Green Valley. CL-AGN hosts have similar star formation properties to the host galaxies of Seyfert 1 and 2 AGN at z < 0.15 and to Seyfert 2 AGN at z > 0.15. We find no statistically significant differences in the star formation properties of turn-on and turn-off CL-AGN. We also find no evidence for rapid quenching in the Green Valley CL-AGN. We conclude that CL-AGN state transitions are not associated with the formation history of CL-AGN host galaxies on large spatial scales, implying CL-AGN state transitions may instead result from nuclear-scale or accretion disk effects.

In this work, we present high-resolution spectroscopic observations for six metal-poor stars with [Fe/H]<-3 (including one with [Fe/H]<-4), selected using narrow-band Ca II HK photometry from the DECam MAGIC Survey. The spectroscopic data confirms the accuracy of the photometric metallicities and allows for the determination of chemical abundances for 16 elements, from carbon to barium. The program stars have chemical abundances consistent with this metallicity range. A kinematic/dynamical analysis suggests that all program stars belong to the distant Milky Way halo population (heliocentric distances 35 < dhelio/kpc < 55), including three with high-energy orbits that might have been associated with the Magellanic system and one, J0026-5445, having parameters consistent with being a member of the Sagittarius stream. The remaining two stars show kinematics consistent with the Gaia-Sausage/Enceladus dwarf galaxy merger. J0433-5548, with [Fe/H]=-4.12, is a carbon-enhanced ultra metal-poor star, with [C/Fe]=+1.73. This star is believed to be a bona fide second-generation star, and its chemical abundance pattern was compared with yields from metal-free supernova models. Results suggest that J0433-5548 could have been formed from a gas cloud enriched by a single supernova explosion from a ~11Mo star in the early universe. The successful identification of such objects demonstrates the reliability of photometric metallicity estimates, which can be used for target selection and statistical studies of faint targets in the Milky Way and its satellite population. These discoveries illustrate the power of measuring chemical abundances of metal-poor Milky Way halo stars to learn more about early galaxy formation and evolution.

The last decade has witnessed significant progress in our understanding of the growth of super-massive black holes (SMBHs). It is now clear that an Active Galactic Nucleus (AGN: the observed manifestation of a growing SMBH) is an "event" within the broader lifecycle of a galaxy, which can significantly influence the shape and evolution of the galaxy itself. Our view of the obscuring medium that affects the observed properties of an AGN has also undergone a revolution, and we now have a more physical understanding of the connection between the fuelling of (and feedback from) the SMBH and the broader host-galaxy and larger-scale environment. We have a greater understanding of the physics of SMBH accretion, can identify AGNs out to z = 8-10 witnessing the very earliest phases of SMBH growth, and have a more complete census of AGN activity than ever before. This great progress has been enabled by new innovative facilities, an ever-increasing quantity of multi-wavelength data, the exploitation and development of new techniques, and greater community-wide engagement. In this article we review our understanding of AGNs and the growth of SMBHs, providing an update of the earlier Alexander and Hickox (2012) review. Using citation-network analyses we also show where this review fits within the broader black-hole research literature and, adopting the previous article as a snapshot of the field over a decade ago, identify the drivers that have enabled the greatest scientific progress.

OVz stars are identified by their optical spectra, in which the He II 4686 absorption line is stronger than any other He line. This spectral characteristic suggests that these stars may be less luminous, less evolved, and located closer to the zero-age main sequence (ZAMS). Although their evolutionary status is still debated, OVz stars are key candidates for understanding the earliest stages of massive star evolution. In this study, we analyze OVz stars in the star-forming regions NGC 346 in the Small Magellanic Cloud (SMC) and 30 Doradus in the Large Magellanic Cloud (LMC), aiming to determine their physical properties and investigate the role of low metallicity in the emergence of the Vz peculiarity. We identified a sample of OVz and OV stars in NGC 346 using spectra obtained with the Magellan Echellette (MagE) spectrograph at Las Campanas Observatory. Spectral classification was based on the equivalent width ratios of key He lines. Quantitative spectroscopic analyses were carried out using the IACOB-BROAD and IACOB-GBAT/FASTWIND tools. For consistency, previously identified OV and OVz stars in 30 Doradus were reclassified following the same criteria.

The ultraviolet and optical wavelength ranges have proven to be a key addition to infrared observations of exoplanet atmospheres, as they offer unique insights into the properties of clouds and hazes and are sensitive to signatures of disequilibrium chemistry. Here we present the 0.2-0.8 $\mu$m transmission spectrum of the Teq = 2000 K Jupiter KELT-7b, acquired with HST WFC3/UVIS G280 as part of the HUSTLE Treasury program. We combined this new spectrum with the previously published HST WFC3/IR G141 (1.1-1.7 $\mu$m) spectrum and Spitzer photometric points at 3.6$\mu$m and 4.5$\mu$m, to reveal a generally featureless transmission spectrum between 0.2 and 1.7 $\mu$m, with a slight downward slope towards bluer wavelengths, and a asymmetric water feature in the 1.1-1.7 $\mu$m band. Retrieval models conclude that the 0.2 - 1.7$\mu$m spectrum is primarily explained by a high H- abundance ($\sim 10^{-5}$), significantly above the equilibrium chemistry prediction ($\sim 10^{-12}$), suggesting disequilibrium in KELT-7b's upper atmosphere. Our retrievals also suggest the presence of bright inhomogeneities in the stellar surface, and tentative evidence of CO2 at the Spitzer wavelengths. We demonstrate that with the UV-optical coverage provided by WFC3 UVIS/G280, we are able to confirm the presence and constrain the abundance of H-, and obtain evidence for bright stellar inhomogeneities that would have been overlooked using infrared data alone. Observations redward of 1$\mu$m with JWST should be able to further constrain the abundance of H-, as well as confirm the presence of CO2 inferred by the two Spitzer datapoints.

The $\nu$HDM is the only cosmological model based on Milgromian Dynamics (MOND) with available structure formation simulations. While MOND accounts for galaxies, with a priori predictions for spirals and ellipticals, a light sterile neutrino of 11 eV can assist in recovering scaling relations on the galaxy-cluster scales. In order to perform MONDian cosmological simulations in this theoretical approach, initial conditions derived from a fit to the angular power spectrum of Cosmic Microwave Background (CMB) fluctuations are required. In this work, we employ CosmoSIS to perform a Bayesian study of the $\nu$HDM model. Using the best-fit values of the posterior, the CMB power spectrum is reevaluated. The excess of power in the transfer function implies a distinct evolution scenario, which can be used further as an input for a set of hydro-dynamical calculations. The resulting values H0 $\approx$ 56 km/s/Mpc and ${\Omega}_{m_{0}} \approx 0.5$ are far from agreement with respect to the best fit ones in the canonical Cold Dark Matter model, but may be significant in MONDian cosmology. The assumed Planck CMB initial conditions are only valid for the $\Lambda$CDM cosmology. This work constitutes a first step in an iterative procedure needed to disentangle the model dependence of the derived initial density and velocity fields.

Using relative stellar astrometry for the detection of coherent gravitational wave sources is a promising method for the microhertz range, where no dedicated detectors currently exist. Compared to other gravitational wave detection techniques, astrometry operates in an extreme high-baseline-number and low-SNR-per-baseline limit, which leads to computational difficulties when using conventional Bayesian search techniques. We extend a technique for efficiently searching pulsar timing array datasets through the precomputation of inner products in the Bayesian likelihood, showing that it is applicable to astrometric datasets. Using this technique, we are able to reduce the total dataset size by up to a factor of $\mathcal{O}(100)$, while remaining accurate to within 1% over two orders of magnitude in gravitational wave frequency. Applying this technique to simulated astrometric datasets for the Kepler Space Telescope and Nancy Grace Roman Space Telescope missions, we obtain forecasts for the sensitivity of these missions to coherent gravitational waves. Due to the low angular sky coverage of astrometric baselines, we find that coherent gravitational wave sources are poorly localized on the sky. Despite this, from $10^{-8}$ Hz to $10^{-6}$ Hz, we find that Roman is sensitive to coherent gravitational waves with an instantaneous strain above $h_0 \simeq 10^{-11.4}$, and Kepler is sensitive to strains above $h_0 \simeq $ $10^{-12.4}$. At this strain, we can detect a source with a frequency of $10^{-7}$ Hz and a chirp mass of $10^9$ $M_\odot$ at a luminosity distance of 3.6 Mpc for Kepler, and 0.3 Mpc for Roman.

We present the numerical lunar time ephemeris LTE440 based on the definition of Lunar Coordinate Time (TCL) given by the International Astronomical Union (IAU) in IAU 2024 Resolution II. LTE440 can be used to obtain the numerical transformation between TCL and Solar System Barycentric Dynamical Time (TCL-TDB) or Solar System Barycentric Coordinate Time (TCL-TCB). The theoretical model, numerical method, and performance of LTE440 are discussed. The secular drifts between TCL and TDB and between TCL and TCB are respectively estimated as $\left< d\mathrm{TCL}/d\mathrm{TDB}\right>=1-6.798\,355\,238\times10^{-10}$, and $\left< d\mathrm{TCL}/d\mathrm{TCB}\right>=1-1.482\,536\,216\,67\times10^{-8}$. The most significant periodic terms have amplitudes of about 1651 microseconds for the annual term and 126 microseconds for the monthly term. The precision of LTE440 is estimated to be at the level of several picoseconds. LTE440 is freely available at this https URL.

The distribution of galaxies in the Zone of Avoidance (ZoA) is incomplete due to the presence of our own Galaxy. Our research focused on the identification and characterisation of galaxies in the ZoA, using the new near-infrared data from the VVVX survey in the regions that cover the southern Galactic disc. We used our previously-established procedure based on photometric and morphological criteria to identify galaxies. The large data volume collected by the VVVX required alternatives to visual inspection, including artificial intelligence techniques, such as classifiers based on neural networks. The VVV NIR galaxy catalogue is presented, covering the southern Galactic disc, significantly extending the vision down to $K^0_s=16$ mag throughout the ZoA. This catalogue provides positions, photometric and morphological parameters for a total of 167,559 galaxies with their probabilities determined by the CNN and XGBoost algorithms based on image and photometric data, respectively. The construction of the catalogue involves the employment of optimal probability criteria. 14% of these galaxies were confirmed by visual inspection or by matching with previous catalogues. The peculiarities exhibited by distinct regions across the Galactic disc, along with the characteristics of the galaxies, are thoroughly examined. The catalogue serves as a valuable resource for extragalactic studies within the ZoA, providing a crucial complement to the forthcoming radio catalogues and future surveys utilizing the Vera C.~Rubin Observatory and the Nancy Grace Roman Space Telescope. We present a deep galaxy map, covering a 1080 sq. deg. region, which reveals that the apparent galaxy density is predominantly influenced by foreground extinction from the Milky Way. However, the presence of intrinsic inhomogeneities, potentially associated with candidate galaxy groups or clusters and filaments, is also discernible.

Motivated by the discovery of a pulsar in the direction of the old open cluster NGC 6791, we conducted a search for radio pulsars in archival Parkes observations targeting similar old open clusters. We reprocessed 224 observations totalling 75.02 hours from four clusters: Theia 1661, NGC 6259, Pismis 3, and Trumpler 20. Our analysis identified five known pulsars and three new rotating radio transient (RRAT) candidates. By comparing the measured dispersion measures (DMs) with the expected DM values for each cluster derived from YMW16 and NE2001 models, we conclude that most detected sources are likely background pulsars. However, RRAT J1749-25 in Theia 1661 and RRAT J1237-60 in Trumpler 20 have DMs reasonably close to their respective clusters, suggesting possible membership. The association between PSR J1750-2536 and Theia 1661 remains ambiguous due to its intermediate DM. These candidate cluster-associated neutron stars warrant follow-up with more sensitive telescopes such as MeerKAT or the SKA, potentially offering valuable insights into neutron star retention mechanisms and evolution in open cluster environments.

The Galactic black hole (BH) X-ray binary 4U 1630-47 went into a new outburst in 2021 after $\sim$ 600 days from its 2020 outburst. We perform a detailed analysis of quasi-periodic oscillations and spectral evolutions during its 2021 outburst based on \textit{Insight}-HXMT observations. The main science aims to study the reflection features evolution of this accreting black hole using the observations of detecting quasi-periodic oscillations (QPOs) and quasi-regular modulations (QRMs). The QPOs frequencies evolve from $\sim 1.6 - 3.6$ Hz, and QRMs have low frequencies around 0.05 - 0.07 Hz. The reflection fraction varies during the outburst and has a positive correlation with the hardness ratio when QPOs are detected. The centroid frequency of QPOs is anti-correlated to the reflection fraction. This is consistent with the prediction of precessing inner flow model and provides evidence for a geometrical origin of QPOs. The centroid frequency of QRMs also shows an anti-correlation to the reflection fraction, but the hardness ratio shows no relation to the reflection fraction during the period. We suggest that QRMs may have a different origin from QPOs and be caused by instabilities in the corona.

The fundamental nature of dark matter (DM) remains unknown, with significant uncertainties in its density profile. DM environments surrounding massive binary black holes (BBHs) modify their orbital dynamics, thereby altering gravitational wave (GW) emissions. For BBH systems at the Galactic Center, dynamical friction induced by DM spikes could produce detectable deviations in GW spectra, potentially observable by future space-based detectors. To address the uncertainties in the Galactic Center's DM profile, we systematically examine two scenarios: the generalized Navarro-Frenk-White (gNFW) profile and its post-spike modification. We investigate the evolutionary effects of DM dynamical friction and accretion on the eccentricity and semi-latus rectum of secondary black holes (BHs) in elliptical orbits. By constructing orbital models with varying initial eccentricities across the mass-semi-latus rectum parameter space and utilizing 30 years of simulated pulsar timing array data from the Square Kilometer Array (SKA), we identify detectable parameter regimes of DM effects and employ these GW observational signatures to constrain different DM density profiles. Our analysis reveals that among gNFW profiles ($\gamma=2,1.5,1,0.5$), only $\gamma=2$ produces significant detectable signatures. The formation of DM spikes further enhances these observable waveform deviations for all gNFW slopes.

Recent studies used Illustris TNG50 galaxies to study bar formation, evolution, and properties like length, strength, and pattern speed. In simulations, these are typically derived from particle positions and mass distribution, neglecting stellar light and extinction effects. However, observational studies indicate that bar appearance depends on wavelength. To test whether this dependence exists in TNG50 at $z \sim 0$, we analysed 50 strongly barred galaxies using mock images from SKIRT radiative transfer simulations covering infrared to ultraviolet filters (Spitzer 3.6 um, SDSS i, r, g, S-PLUS J0378, GALEX NUV, and GALEX FUV). Bar ellipticity and length were measured via ellipse fitting. Ellipticity generally increases by 6 percent from 3.6 um to g band, and by 9 percent to J0378 band. On average, TNG50 bars cannot be said to be longer in bluer filters when the entire sample is used. However, the trend is detected when only star-forming galaxies are considered. In this star-forming subsample, bar length increases by 10 percent from 3.6 um to g band, and by 17 percent to J0378 band; moreover, the bar can appear up to 20${-}$30 percent longer in bluer mocks than in the mass map. Over 90 percent of bars vanish in UV due to minimal emission by dominant stellar populations. We reproduced the bar properties morphological dependence phenomenon using age-filtered mass maps, where older stars form shorter, rounder bars, and younger stars generate longer, more elliptical ones. TNG50 bars exhibit a wavelength-dependent trend similar to observations: bars appear more elliptical and longer in bluer filters, with this effect being stronger in star-forming galaxies.

We investigate the formation of X-ray emission lines in the wind of the Wolf-Rayet (WR) companion in Cyg X-3 by analyzing their orbital dynamics using Chandra High Energy Transmission Grating (HEG) observations during a hypersoft state. Our goal is to constrain the X-ray transparency of the recently discovered funnel-like structure surrounding the compact star, as revealed by X-ray polarimetry. All lines exhibit sinusoidal orbital modulation, with the velocity amplitude generally increasing and the orbital phase with the highest blueshift generally decreasing for ions with a higher ionisation potential. The {Fe}{xxvi}-line displays velocity extremes at phase 0.25 (blueshift) and 0.75 (redshift), indicating that the line-emitting region is close to the compact component (disc or corona) and thus reflects the orbital motion. The {Fe}{xxv}-line shows a complex behaviour that cannot be fully resolved with the Chandra/HEG resolution. Other lines display velocity extremes scattered around phase 0.5 (blueshift) and 0.0 (redshift), with velocity amplitudes of 100--300 km/s, suggesting their origin in the WR stellar wind between the two components along 3D \xi-surfaces. Parts of the emission lines of {Ar}{xviii} and {Ca}{xx} originated around the compact star (disc or corona). The recent polarisation funnel-modelling is consistent with the present results during the hypersoft state.

AI-enhanced approaches are becoming common in astronomical data analysis, including in the galaxy morphological classification. In this study we develop an approach that enhances galaxy classification by incorporating an image denoising pre-processing step, utilizing the U-Net Variational Autoencoder (VAE) architecture and effectively mitigating noise in galaxy images and leading to improved classification performance. Our methodology involves training U-Net VAEs on the EFIGI dataset. To simulate realistic observational conditions, we introduce artifacts such as projected stars, satellite trails, and diffraction patterns into clean galaxy images. The denoised images generated are evaluated using Peak Signal-to-Noise Ratio (PSNR) and Structural Similarity Index (SSIM), to quantify the quality improvements. We utilize the denoised images for galaxy classification tasks using models such as DenseNet-201, ResNet50, VGG16 and GCNN. Simulations do reveal that, the models trained on denoised images consistently outperform those trained on noisy images, thus demonstrating the efficiency of the used denoising procedure. The developed approach can be used for other astronomical datasets, via refining the VAE architecture and integrating additional pre-processing strategies, e.g. in revealing of gravitational lenses, cosmic web structures.

We analyze the impact of the decaying magnetic turbulence of primordial magnetic fields (PMFs) and ambipolar diffusion on the ionization and thermal history of the Dark Ages Universe ($30\le z\le300$), and its imprint on the spectral profile of the global signal in the 21 cm hydrogen line. The heating function caused by decaying magnetic turbulence monotonically decreases after cosmological recombination, its amplitude strongly depends on the strength of the PMFs $B_0$ and weakly depends on the spectral index of the initial power spectrum of PMFs $n_{\rm B}$. The heating function caused by ambipolar diffusion, in contrast, noticeably depends on the spectral index in the range $-3\lesssim n_{\rm B}\lesssim-1$, but is subdominant in the Dark Ages epoch for PMF models with $B_0\lesssim0.5$ nG. We computed the ionization and thermal history of intergalactic gas from the cosmological recombination up to the end of the Dark Ages epoch for a range of PMF parameters $0.05\lesssim B_0\lesssim0.5$ nG, $-2.9\lesssim n_{\rm B}\lesssim4$ and show the essentially distinguished thermal evolution from one in the $\Lambda$CDM model. We also show that the profile of the redshifted 21 cm hydrogen line is very sensitive to the PMF parameters from this range and can be used for their constraints in conjunction with other observational data.

Line resolved X-ray spectra of outflows from X-ray binaries are interesting since they provide quantifiable measures of the accreted material on to the compact object (black hole or neutron star), which can not be observed directly in the accretion disk. One such measurement that has been largely overlooked is that of the elemental abundances, which potentially provide insights into the origin of the ejected material. Using the Chandra/HETG grating spectrometer we measure and present elemental abundances in four low-mass X-ray binaries. We compare two measurement methods. One is by fitting line series of individual ions and reconstructing the absorption measure distribution (AMD), and the other is a global fit with one or two individual ionization components. All outflows feature a steep AMD strongly favoring high ionization degrees. The present abundances are consistent with previous works suggesting the abundances in the outflows are non-solar. We find a tentative trend of increasing abundances with atomic number, which fits some core-collapse supernova models, but no exact match to a specific one.

We present International LOw Frequency ARray (LOFAR) Telescope (ILT) observations of the Crab Nebula, the remnant of a core-collapse supernova explosion observed by astronomers in 1054. The field of the Crab Nebula was observed between 120 and 168 MHz as part of the LOFAR Two Meter Sky Survey (LoTSS), and the data were re-processed to include the LOFAR international stations to create a high angular resolution ($0.43'' \times 0.28''$) map at a central frequency of 145 MHz. Combining the ILT map with archival centimeter-range observations of the Nebula with the Very Large Array (VLA) and LOFAR data at 54 MHz, we become sensitive to the effects of free-free absorption against the synchrotron emission of the pulsar wind nebula. This absorption is caused by the ionised filaments visible in optical and infrared data of the Crab Nebula, which are the result of the pulsar wind nebula expanding into the denser stellar ejecta that surrounds it and forming Rayleigh-Taylor fingers. The LOFAR observations are sensitive to two components of these filaments: their dense cores, which show electron densities of $\gtrsim1,000$ cm$^{-3}$, and the diffuse envelopes, with electron densities of $\sim50-250$ cm$^{-3}$. The denser structures have widths of $\sim0.03$ pc, whereas the diffuse component is at one point as large as $0.2$ pc. The morphology of the two components is not always the same. These finding suggests that the layered temperature, density, and ionisation structure of the Crab optical filaments extends to larger scales than previously considered.

Metal-bearing species in diffuse or molecular clouds are often overlooked in astrochemical modeling except for the charge exchange process. However, catalytic cycles involving these metals can affect the abundance of other compounds. We prepared a comprehensive chemical network for Na, Mg, Al, Fe, K, and Si-containing species. Assuming water as the major constituent of interstellar ice in dark clouds, quantum chemical calculations were carried out to estimate the binding energy of important metallic species, considering amorphous solid water as substrate. Significantly lower binding energies (approximately 5 to 6 times) were observed for Na and Mg, while the value for Fe was roughly 4 times higher than what was used previously. Here, we calculated binding energy values for Al and K, for which no prior guesses were available. The total dipole moments and enthalpies of formation for several newly included species are unknown. Furthermore, the assessment of reaction enthalpies is necessary to evaluate the feasibility of the new reactions under interstellar conditions. These parameters were estimated and subsequently integrated into models. Some additional species that were not included in the UMIST/KIDA database have been introduced. The addition of these new species, along with their corresponding reactions, appears to significantly affect the abundances of related species. Some key reactions that significantly influence general metal-related chemistry include: $\rm{M^+ + H_2 \rightarrow M{H_2}^+ + h\nu}$, $\rm{MH + O \rightarrow MO + H}$ ($\rm{M \ = \ Fe, \ Na, \ Mg, \ Al, \ K}$), and $\rm{M_1^+ + M_2H \rightarrow M_1H + M_2^+}$ (where $\rm{M_1 \neq M_2}$, $M_1, \ M_2 \ = \ Na, \ Mg, \ Al, \ K, \ Fe$). Significant changes were observed in magnesium and sodium-bearing cyanides, isocyanides, and aluminum fluoride when additional reaction pathways were considered.

Bright points (BPs) are ubiquitous, small-scale energetic events with multithermal signatures, typically observed in the chromosphere and linked to both photospheric structure and coronal composition. Their evolution is influenced by various physical processes, including plasma dynamics and magnetic interactions. This paper examines BP evolution using a large statistical sample, focusing on when they reach maximum values for key attributes, and explores differences between BPs in the "Active Quiet Sun" (AQS, above the network) and the "True Quiet Sun" (TQS, above the internetwork). Observed attributes include maximum brightness (total and intrinsic), plane-of-sky (POS) speed, travel distance, acceleration, and apparent size (POS area). BPs can reach maximum brightness and size at almost any time during their lifetime, likely due to complex chromospheric interactions. AQS and TQS BPs show similar behaviour overall and tend to reach maximum POS speed around the midpoint of their lifetimes. Positive acceleration usually occurs early, while negative acceleration is more common near the end. Preliminary results suggest two distinct BP regimes with differing relationships between intrinsic brightness and area. We interpret these trends as evidence that BPs, whether due to plasma motion or a heating event, follow arched paths -- likely along short magnetic loops. In this model, the midpoint of a BP's life corresponds to the crest of the arch, producing the greatest POS speeds. Reconnection may occur at a footpoint, with the BP moving along the loop before returning to the photosphere. Deviations from this expected evolution may result from complex chromospheric interactions or BPs with highly non-linear POS motions.

Context. Cosmic voids are vast underdense regions in the cosmic web that encode crucial information about structure formation, the composition of the Universe, and its expansion history. Due to their lower density, these regions are less affected by non-linear gravitational dynamics, making them suitable candidates for analysis using semi-analytic methods. Aims. We assess the accuracy of the PINOCCHIO code, a fast tool for generating dark matter halo catalogs based on Lagrangian Perturbation Theory, in modeling the statistical properties of cosmic voids. We validate this approach by comparing the resulting void statistics measured from PINOCCHIO to those obtained from N-body simulations. Methods. We generate a set of simulations using PINOCCHIO and OpenGADGET3, assuming a fiducial cosmology and varying the resolution. For a given resolution, the simulations share the same initial conditions between the different simulation codes. Snapshots are saved at multiple redshifts for each simulation and post-processed using the watershed void finder VIDE to identify cosmic voids. For each simulation code, we measure the following statistics: void size function, void ellipticity function, core density function, and the void radial density profile. We use these statistics to quantify the accuracy of PINOCCHIO relative to OpenGADGET3 in the context of cosmic voids. Results. We find agreement for all void statistics at better than 2{\sigma} between PINOCCHIO and OpenGADGET3, with no systematic difference in redshift trends. This demonstrates that the PINOCCHIO code can reliably produce void statistics with high computational efficiency compared to full N-body simulations.

We explore the diagnostic potential of the H$\alpha$ line for probing the chromospheric magnetic field using a realistic 3D radiative magnetohydrodynamic (rMHD) model. The Stokes profiles of the H$\alpha$ line are synthesized through full 3D radiative transfer under the field-free approximation, alongside the Ca II 8542 Å and Fe I 6173 Å lines for comparison. The line-of-sight (LOS) magnetic fields are inferred using the weak field approximation (WFA) for theH$\alpha$ and Ca II 8542 Å lines, while the Fe I 6173 Å line is analyzed through Milne-Eddington inversion techniques. The comparison between the inferred LOS magnetic field maps and the magnetic fields in the rMHD model revealed that the H$\alpha$ line core primarily probes the chromospheric magnetic field at log tau_500 = -5.7, which corresponds to higher layers than the Ca II 8542 Å line core, which is most sensitive to conditions at log tau_500 = -5.1. On average, the Stokes V profiles of the H$\alpha$ line core form 500 km higher than those of the Ca II 8542 Å line core. The H$\alpha$ polarization signals persist after adding noise, and with noise at the level of 10^-3 Ic, most simulated magnetic structures remain visible. These findings suggest that spectropolarimetric observations of the H$\alpha$ line can provide complementary insights into the stratification of the magnetic field at higher altitudes, especially when recorded simultaneously with widely used chromospheric diagnostics such as the Ca II 8542 Å line.

We present a new method to calculate the polarised synchrotron emission of radio AGN sources using magnetic field information from 3-dimensional relativistic magnetohydrodynamical (RMHD) simulations. Like its predecessor, which uses pressure as a proxy for the magnetic field, this method tracks the spatially resolved adiabatic and radiative loss processes using the method adapted from the Radio AGN in Semi-analytic Environments formalism. Lagrangian tracer particles in RMHD simulations carried out using the PLUTO code are used to track the fluid quantities of each `ensemble of electrons' through time to calculate the radio emissivity ex-situ. By using the magnetic field directly from simulations, the full set of linear Stokes parameters I, Q, and U can be calculated to study the synthetic radio polarisation of radio AGN sources. We apply this method to a suite of RMHD simulations to study their polarisation properties. The turbulent magnetic field present in radio lobes influences the emission, causing a complex clumpy structure that is visible at high resolution. Our synthetic polarisation properties are consistent with observations; we find that the fractional polarisation is highest (approximately 50 percent) at the lobe edges. We show that for the same source, the integrated and mean fractional polarisation depends on viewing angle to the source. At oblique viewing angles the behaviour of the integrated and mean fractional polarisation over time depends on the morphology of the jet cocoon. Using Faraday rotation measures, we reproduce known depolarisation effects such as the Laing-Garrington depolarisation asymmetry in jets angled to the line of sight. We show that the hotspots and hence the Fanaroff-Riley classification become less clear with our new, more accurate method.

Common envelope (CE) evolution is largely governed by the drag torque applied on the inspiralling stellar components by the envelope. Previous work has shown that idealized models of the torque based on a single body moving in rectilinear motion through an unperturbed atmosphere can be highly inaccurate. Progress requires new models for the torque that account for binarity. Toward this end we perform a new 3D global hydrodynamic CE simulation with the mass of the companion point particle set equal to the mass of the asymptotic giant branch star core particle to maximize symmetry and facilitate interpretation. First, we find that a region around the particles of a scale comparable to their separation contributes essentially all of the torque. Second, the density pattern of the torque-dominating gas and, to an extent, this gas itself, is roughly in corotation with the binary. Third, approximating the spatial distribution of the torquing gas as a uniform-density prolate spheroid whose major axis resides in the orbital plane and lags the line joining the binary components by a constant phase angle reproduces the torque evolution remarkably well, analogous to studies of binary supermassive black holes. Fourth, we compare the torque measured in the simulation with the predictions of a model that assumes two weak point-mass perturbers undergoing circular motion in a uniform background without gas self-gravity, and find remarkable agreement with our results if the background density is taken to be equal to a fixed fraction ($\approx0.44$) of the density at the spheroid surface. Overall, this work makes progress toward developing simple time-dependent models of the CE phase.

We self-consistently determine the 3D density distribution of NGC 551's stellar disk and study observational signatures of two-component stellar disks. Assuming baryonic disks are in hydrostatic equilibrium, we solved the Poisson-Boltzmann equation to estimate 3D density distribution. We used integral-field spectroscopic observations to estimate stellar velocity dispersion and built a 3D dynamical model using these density solutions and the observed rotation curve. We generated simulated surface brightness maps and compared them with observations to verify modeling consistency. The dynamical model was inclined to 90° to produce an edge-on surface density map, which we investigated by fitting different 2D functions and plotting vertical cuts in logarithmic scale. We estimated vertical stellar velocity dispersion using an iterative method, obtaining results consistent with the Disk Mass Survey formalism. Through dynamical modeling, we produced moment maps that reasonably matched observations. We examined the simulated edge-on model by taking vertical cuts and decomposing them into multiple Gaussian components. We find that artificial double Gaussian components arise due to line-of-sight integration effects, even for single-component disks. This indicates that decomposing vertical intensity cuts into multiple Gaussian components is unreliable for multicomponent disks. Instead, an up-bending break visible in logarithmic-scale vertical cuts serves as a more reliable indicator for two-component disks. We performed 2D fitting on the edge-on surface density map using the product of a scaled modified Bessel function and $sech^2$ function to estimate structural parameters. These traditional methods systematically underestimate the scale length and flattening ratio. Therefore, we suggest using detailed modeling to accurately deduce stellar disk structural parameters.

We present a detailed X-ray/UV high-time resolution monitoring of the final reflaring phase of the accreting millisecond pulsar SAX~J1808.4$-$3658. During its 2022 outburst, we obtained simultaneous XMM-Newton and Hubble Space Telescope (HST) observations. We detected coherent X-ray pulsations down to a 0.5-10~keV luminosity of $L_{X(low)\,0.5-10} \simeq 6.21^{+0.20}_{-0.15} \times 10^{34} \, d^2_{3.5}\,\mathrm{erg \, s^{-1}}$, among the lowest ever observed in this source. The uninterrupted coverage provided by XMM-Newton enabled a detailed characterisation of the spectral and temporal evolution of the source X-ray emission as the flux varied by approximately one order of magnitude. At the lowest flux levels, we observed significant variations in pulse amplitude and phase correlated with the X-ray source flux. We found a sharp phase jump of $\sim 0.4$ cycles, accompanied by a doubling of the pulse amplitude and a softening of the X-ray emission. We interpret the changes in the X-ray pulse profiles as drifts of emission regions on the neutron star surface due to an increase of the inner disk radius occurring when the mass accretion rate decreases. The phase evolution was consistent with a magnetospheric radius scaling as $R_{m} \propto \dot{M}^{\Lambda}$, with $\Lambda = -0.17(9)$, in broad agreement with theoretical predictions. Simultaneous HST observations confirmed the presence of significant UV pulsations. The measured pulsed luminosity $-$ $L_{pulsed}^{UV} = (9 \pm 2) \times 10^{31} \, \text{erg} \, \text{s}^{-1}$ $-$ was approximately half that observed during the 2019 outburst, but the pulsed X-ray to UV luminosity ratio simultaneously measured remained consistent. Yet, such a UV luminosity exceeds the predictions of standard emission models, as further confirmed by the shape of the pulsed spectral energy distribution.

Dynamical interactions between stars and the super massive black hole Sgr A* at the Galactic Centre (GC) may eject stars into the Galactic halo. While recent fast ejections by Sgr A* have been identified in the form of hypervelocity stars (hundreds to thousands km/s), it is also expected that the stellar halo contains slower stars, ejected over the last few billion years. In this study, we use the first data release of DESI to search for these slower GC ejecta, which are expected to stand out from the stellar halo population for their combined high metallicity (${\rm [Fe/H]}\gtrsim0$) and small values of their vertical angular momentum ($L_Z$), whose distribution should peak at zero. Our search does not yield a detection, but allows us to place an upper limit on the ejection rate of stars from the GC of $\sim2.8\times10^{-3}$ yr$^{-1}$ over the past ~5 Gyr, which is ejection model independent. This implies that our result can be used to put constraints on different ejection models, including that invoking mergers of Sgr A* with other massive black holes in the last last few billion years.

The nature of dark energy remains one of the most profound mysteries in modern cosmology. One intriguing proposal is that black holes (BHs) could be the astrophysical source of dark energy through a cosmological coupling mechanism, and strong evidence has been claimed via analyzing the growth of the black hole masses in the red-sequence elliptical galaxies at redshifts $\leq 2.5$. In this work, with a group of very high redshift AGNs detected by the James Webb Space Telescope (JWST) in the red-sequence elliptical galaxies, we show that the possibility of BHs being the astrophysical source of dark energy has been rejected at a confidence level exceeding 10$\sigma$. Moreover, it turns out that the Little Red Dots recently discovered by JWST, characterized by the low accretion rates, can naturally evolve into the red-sequence elliptical galaxies hosting the relatively low mass black holes at the redshifts of $\sim 0.7-2.5$, without the need of black hole cosmological coupling.

We present the first broad-band spectral and timing study of the Be/X-ray pulsar XTE J0111.2$-$7317 (SXP31.0) during its first major outburst since its discovery in 1998. This giant Type II outburst, observed in April-June 2025, marks the source's return to activity after nearly three decades of quiescence. Utilizing our NuSTAR observations together with Swift/XRT and SRG/ART-XC ones, we trace the evolution of the outburst to a peak bolometric luminosity of $L_{\rm bol} = 3.6 \times 10^{38}$ erg s$^{-1}$. The broadband spectra are well described by an absorbed cutoff power law, two blackbody components (hot and soft), and a narrow Fe K$\alpha$ line. No cyclotron absorption features are detected in either the phase-averaged or phase-resolved spectra in the 5-50 keV band. Most notably, we report the discovery of a previously undetected quasi-periodic oscillation (QPO) at $0.8 \pm 0.1$ mHz, characterized by a fractional root-mean-square amplitude of 18% at a super-Eddington bolometric luminosity of $L_{\rm bol} = 2.5 \times 10^{38}$ erg s$^{-1}$. In contrast, the previously reported 1.27 Hz QPO is not detected. While the sub-mHz QPO is present, the pulsed fraction (PF) is low in soft X-rays, consistent with other super-Eddington pulsars exhibiting mHz QPOs, but rises above 20 keV reaching 35%. The QPO vanishes in subsequent observations coinciding with a sharp increase in the PF and a distinct change in pulse profile morphology. It is not observed in any follow-up observations at luminosities above or below its initial detection, suggesting it is a transient phenomenon.

We investigate broad emission-line wings, reaching velocities up to 800 km/s, observed in 26 galaxies with Lyman continuum (LyC) observations, primarily from the Low-redshift Lyman Continuum Survey (LzLCS). Using Magellan/MIKE, VLT/X-shooter, and WHT/ISIS high-resolution spectroscopy, we show that this fast gas appears to probe the dominant feedback mechanisms linked to LyC escape. We find that in 14 galaxies, the wings are best fit with power laws with slopes -3.5 to -1.6, with four others best fit by Gaussians of width ~300 km/s; the remaining eight show ambiguous wing morphologies. Gaussian wings are found only at low O32 = [O III]5007/[O II]3726,3729 and high metallicity, while power-law wings span the full range of these parameters. The general evidence suggests a dual-mode paradigm for LyC escape: radiation-driven superwinds traced by power-law wings and supernova-driven feedback traced by Gaussian wings. For the former, the < 3 Myr-old, pre-supernova stellar population correlates with more luminous, faster winds. The data also show that radiation-driven wind parameters like wind luminosity and power-law slope depend on the UV luminosity more than the optically thick covering fraction, consistent with ``picket-fence" radiative transfer. Observed wing slope values flatten with both escaping LyC luminosity and higher extinction, while still preserving the anticorrelation between these two quantities. Additionally, the differential between red and blue slopes implies that extinction and dense gas are centrally concentrated relative to the wind emission. Overall, our results show that power-law emission-line wings probe LyC-driven winds and LyC escape in metal-poor starbursts.

Low-m inertial modes have been recently discovered in the Sun's high-latitude regions. In this study, we characterize the observational properties of the m = 1 mode by analyzing time-distance subsurface flow maps. Synoptic flow maps, constructed from daily subsurface flow maps using a tracking rate corresponding to the rotation at latitude 65 degrees, are filtered in both the spherical harmonic and Fourier domains to retain only the m = 1 mode and its dominant frequencies. Our analysis reveals a power distribution that is significantly stronger in the northern polar region. The mode's power exhibits an anti-correlation with solar activity, remaining strong and persistent during the solar activity minimum and becoming weaker and more fragmented during the solar maximum. Magnetic flux transported from low to high latitudes influences both the mode's power and lifetime, enhancing its power and shortening its lifetime upon arrival. The phases of the m = 1 mode in the northern and southern polar regions are near-antisymmetric for most of the time with short deviations. We also compute zonal and meridional phase velocities of the mode and find that it exhibits significantly less differential rotation than its surrounding plasma. The meridional phase velocity, comprising both the local plasma's meridional flow and the mode's intrinsic phase motion, is directed poleward below latitude 70 degrees and equatorward above this latitude. These observational findings underscore the need for a deeper understanding of the internal dynamics of the low-m modes, which may offer valuable insights into the structure and dynamics of the solar interior.

The outer regions of the protoplanetary disc surrounding the T Tauri star HD 143006 show rings, dust asymmetries and shadows. Whilst rings and dust asymmetries can arise from companions and other mechanisms, shadows and misaligned discs in particular are typically attributed to the presence of misaligned planets or stellar-mass companions. To understand the mechanisms that drive these traits, the innermost regions of discs need to be studied. Using CHARA/MIRCX and VLTI/PIONIER, we observed the sub-au region of HD 143006. We constrain the orientation of the inner disc of HD 143006 and probe whether a misalignment between the inner and outer disc could be the cause of the shadows. Modelling the visibilities using a geometric model, the inclination and position angle are found to be $i=22^\circ\pm 3^\circ$ and $\mathrm{PA}=158^\circ\pm 8^\circ$ respectively, with an inner dust sublimation radius of $\sim0.04$ au. The inner disc is misaligned by $39^\circ\pm4^\circ$ with respect to the outer disc, with the far side of the inner disc to the east and the far side of the outer disc to the west. We constrain $h/R$ (scattering surface/radius of scattered light) of the outer disc at $18$ au to be about $13\%$ by calculating the offset between the shadow position and the central star. No companion was detected, with a magnitude contrast of $4.4$ in the H-band and placing an upper mass limit of $0.17 M_\odot$ at separations of $0-8$ au. Therefore, we cannot confirm or rule out that a low-mass star or giant planet is responsible for the misalignment and dust sub-structures.

Aims. We investigate the impact of ambipolar diffusion on the development of the Rayleigh-Taylor instability (RTI) with an oblique magnetic field in the incompressible limit. Methods. We developed a general bi-fluid framework comprising charges and neutrals with the specific feature of differing gravity between charges and neutrals. We derived the perturbed magnetohydrodynamic (MHD) equations and obtained an analytic dispersion relation for an oblique magnetic field. The growth rate was then evaluated using a numerical integration of the dispersion relation. In particular, we focussed on the anisotropy in the mode growth induced ambipolar diffusion. Results. In contrast to the case of a magnetic field within the interface, an oblique magnetic field is much less restrictive with respect to the wavenumbers that can develop; rather than a sharp cut-off, it results in a selection of preferred scales for mode growth. The presence of a second neutral fluid that is insensitive to gravity tends to amplify the anisotropy in the possible direction of the instability development. In particular, we show that this effect is maximal when the coupling is in an intermediate range between full and no charge-neutral coupling, specifically in the region where the ambipolar diffusion is highest.

Increasing resolution and coverage of astrophysical and climate data necessitates increasingly sophisticated models, often pushing the limits of computational feasibility. While emulation methods can reduce calculation costs, the neural architectures typically used--optimised via gradient descent--are themselves computationally expensive to train, particularly in terms of data generation requirements. This paper investigates the utility of the Extreme Learning Machine (ELM) as a lightweight, non-gradient-based machine learning algorithm for accelerating complex physical models. We evaluate ELM surrogate models in two test cases with different data structures: (i) sequentially-structured data, and (ii) image-structured data. For test case (i), where the number of samples $N$ >> the dimensionality of input data $d$, ELMs achieve remarkable efficiency, offering a 100,000$\times$ faster training time and a 40$\times$ faster prediction speed compared to a Bi-Directional Recurrent Neural Network (BIRNN), whilst improving upon BIRNN test performance. For test case (ii), characterised by $d >> N$ and image-based inputs, a single ELM was insufficient, but an ensemble of 50 individual ELM predictors achieves comparable accuracy to a benchmark Convolutional Neural Network (CNN), with a 16.4$\times$ reduction in training time, though costing a 6.9$\times$ increase in prediction time. We find different sample efficiency characteristics between the test cases: in test case (i) individual ELMs demonstrate superior sample efficiency, requiring only 0.28% of the training dataset compared to the benchmark BIRNN, while in test case (ii) the ensemble approach requires 78% of the data used by the CNN to achieve comparable results--representing a trade-off between sample efficiency and model complexity.

Currently, achieving a balance between computational efficiency, accuracy, and numerical stability in CME simulations, particularly in the sub-Alfv{'e}nic coronal region, remains a significant challenge. This paper aims to address the challenge by integrating observational data and developing advanced numerical algorithms, focusing on reproducing large-scale CME evolutions that are consistent with observations in the coronal region. Based on the recently developed fully implicit thermodynamic MHD coronal model (Wang et al. 2025a), we further use an observation-based RBSL flux rope to trigger a CME event during CR 2111. Additionally, we improve the temporal accuracy using a 2nd-order accurate ESDIRK2 method, with the intermediate stage solutions computed by the 2nd-order accurate BDF2 pseudo-time marching method. To enhance the numerical stability of ESDIRK2, we apply approximate linearisation in the implicitly solved intermediate stages. Furthermore, we adjust the time-evolving magnetic field B1 to zero at the end of each physical time step to further validate the extended magnetic field decomposition approach proposed by (Wang et al. 2025a). It is noticed that the model successfully reproduces the CME evolution consistent with white-light coronagraph observations, enables faster-than-real-time CME propagation simulations from solar surface to 0.1 AU using only a few dozen CPU cores on approximately 1 million grid cells, and remains numerically stable in CME simulations involving low-\b{eta} regions. The simulation results show that this novel MHD coronal model, combined with an observation-based magnetic flux rope, is sufficiently numerically stable and computationally efficient to reproduce real CME events propagating through the sub-Alfv{é}nic coronal region. Thus, the observation-based CME model is well suited for practical applications in daily space weather forecasting.

Galaxy cluster cosmology relies on complete and pure samples spanning a large range of masses and redshifts. In Xu et al. (2018) and Xu et al. (2022), we discovered an apparently new population of galaxy groups and clusters with, on average, flatter X-ray surface brightness profiles than known clusters; this cluster population was missed in previous cluster surveys. The discovery of such a new class of objects could have a significant impact on cosmological applications of galaxy clusters. In this work we use a subsample of these systems to assess whether they belong to a new population. We follow up three of these galaxy groups and clusters with high-quality XMM-Newton observations. We produce clean images and spectra and use them for model fitting. We also identify known galaxies, groups, or clusters in the field. The observations reveal that all three systems are composed of multiple groups each, either at the same or at different redshifts. In total, we characterize nine groups. We measure flat surface brightness profiles with slope parameter $\beta < 0.6$; i.e, less than the canonical $\beta = 2/3$. For the two main central groups, we even measure $\beta < 0.4$. When the fluxes for the three observations are split up across the nine identified groups, none of them exceeds the typical flux limit adopted in previous RASS cluster catalogs, $\approx 3 \times 10^{-12}\,\mathrm{erg s^{-1}cm^{-2}}$ in the 0.1$-$2.4 keV energy band. The observations reveal that groups with flat surface brightness profiles exist. Whether they form a new, separate population requires additional follow-up observations of further systems from the Xu et al. sample, given the complexity we have discovered. Such extended low surface brightness systems, as well as multiple systems and projection effects, need to be taken into account when determining selection functions of group and cluster samples.

[Abridged] We analyse the dark matter (DM) halo properties of 127 0.3<z<1.5 star-forming galaxies (SFGs) down to low stellar masses (7<log(Mstar/Msun)<11), using data from the MUSE Hubble Ultra Deep Field Survey and photometry from HST and JWST. We employ a 3D forward modelling approach to analyse the morpho-kinematics of our sample, enabling measurement of individual rotation curves out to 2-3 times the effective radius. We perform a disk-halo decomposition with a 3D parametric model that includes stellar, gas, and DM components, with pressure support corrections. We validate our methodology on mock data cubes generated from idealised disk simulations. We select the best-fitting DM model among six density profiles, including the Navarro-Frenk-White and the generalised alpha-beta-gamma profile of Di Cintio et al. (2014, DC14). Our Bayesian analysis shows that DC14 performs as well as or better than the other profiles in >65% of the sample. We find that the kinematically inferred stellar masses agree with values from SED fitting. We find that 89% of galaxies have DM fractions >50%. For 70% of SFGs, we infer a DM inner slope, gamma < 0.5, indicating cored DM profiles, but no correlation is found between gamma and star formation rate of the sample. The stellar- and concentration-mass relations agree with theoretical expectations, but with larger scatter. We confirm the anticorrelation between halo scale radius and DM density. The halo scale radii and DM surface densities increase with Mstar, while DM densities stay constant. We find tentative evidence of an evolution of the DM density with z, which suggests that the DM halos of intermediate-z systems are denser than those of local galaxies. In contrast, the halo scale radii are z-invariant.

The existence of high galactic latitude molecular clouds has been known for several decades, and studies of their dust and gas distributions reveal complicated morphological structures. Their dynamics involve turbulence even in the absence of internal energy sources such as stars. We study in detail two such clouds, MBM 3 and MBM 16, trying to recover the geometric structure and topology of the gas distribution. In particular, we address the evidence of superthermal asymmetric atomic and molecular line profiles as a result of filament superposition combined with turbulent motions. We use a variety of spectroscopic and imaging archival observations of the gas and dust components. The spectroscopic data set comprises HI 21 cm, 12CO, 13CO, and CH line profiles. We also use archival infrared images to study the dust distribution and temperature. To understand the topology of MBM 3 and MBM 16 we compare molecular and atomic spectra, along with profile decomposition of the HI 21 cm line. Standard tools such as Structure Functions of velocity centroids are used to characterise the turbulence in MBM 3, and channel maps and position-velocity diagrams are employed for elucidating the filament topology of both clouds. The unusually large linewidths previously reported for MBM 3 are due to superposition of individual filaments whose superthermal linewidths are about 1 km/s. In MBM 16, the cloud appears to decompose into two adjacent structures with similar properties. The filaments have a high aspect ratio, with lengths of about 1 pc and widths of about 0.1 pc. In general, the molecular gas is embedded within more extended neutral hydrogen structures. Velocity gradients found within these structures are not necessarily dynamical, convergent flows. Projection effects and topology of the driving flows produce signatures that mimic velocity shears even if they are simply distortions of ordered gas.

The O'Connell effect is a phenomenon in eclipsing binary (EB) systems that consists of unequal maxima in a light curve when out of eclipse. Despite being known for decades and with several theories proposed over the years, this effect is still not fully understood. Our goal is to find different O'Connell effect properties using a multiband approach, compare them with each other, and find correlations between these properties and the physical parameters of the systems. We present the analysis of 14 new EBs that show the O'Connell effect using multiband data extracted from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and Zwicky Transient Facility (ZTF) all-sky surveys. We measured the difference in maximum amplitudes (delta m) alongside different light curve features in different passbands via a new modeling process that uses Gaussian fits. We created a brand-new phenomenological model for O'Connell effect systems whose analysis had previously been hampered by overfitting. Although the magnitude of the O'Connell effect seems to be more pronounced at shorter effective wavelengths, supporting the idea that cool starspots cause the effect, a conclusive correlation is not found. On the other hand, we do find strong correlations between the magnitude of the effect and both the temperature and the period, both of which are inconsistent with a previous study. We also find that in systems that show both positive and negative O'Connell effects, there are different correlations with the aforementioned parameters. We conclude that, even though starspots may be one cause of the O'Connell effect, it is likely a multipronged phenomenon; for instance, the physical interaction between the components of close binary stars may be another important factor.

Millimeter (mm) emission from F - M dwarfs (cool stars) primarily traces chromospheric activity, with thermal emission thought to dominate in quiescence. Despite the high chromospheric activity, the quiescent mm spectral fluence (mm-S($\nu$)) of young (< 1 Gyr) M dwarfs (dMs) remain largely unexplored. We present the quiescent mm-S($\nu$) of a young dM, ADLeo, observed around 94 GHz using the Northern Extended Millimetre Array (NOEMA). The observed quiescent mm-S($\nu$) exceeds the thermal flux density from a 1D chromospheric model, constrained by optical-UV spectroscopic data, by up to a factor of 7. This indicates a quasi-steady non-thermal emission powered by supra-thermal electrons unlike in old (> 1 Gyr) cool stars, whose quiescent mm-S($\nu$) generally agree with 1D thermal models. The mm-brightness temperature spectral index ($\alpha_{mm}$; $T_B(\nu)\propto \nu^{- \alpha_{mm}}$) of AD Leo deviates by a factor of 3 from the $\alpha_{mm}$ - $T_{eff}$ scaling law for old sun-like stars (Mohan, A., et al., 2022), while UV Ceti, an older M6V star, follows the trend. Also, we report a double-hump flare with second-scale variability in flux density and spectral index, and a frequency-rising nature with brightness increasing with frequency. The flare resemble certain solar events, but is unlike the second-scale events reported in dMs. The non-thermal flare humps suggest multiple injections of accelerated electrons. The mean flare luminosity (2 - 5 $\times 10^{15} erg s^{-1} Hz^{-1}$) and duration ($18\pm 2$ s) are comparable to flares reported in AU Mic and Proxima Cen, but 100 - 1000 times weaker than the minutes-long dM flares observed by the South Pole Telescope.

External photoevaporation of protoplanetary discs, by massive O stars in stellar clusters, is thought to be a significant process in the evolution of a disc. It has been shown to result in significant mass loss and disc truncation, ultimately reducing the lifetime of the discs, and possibly affecting potential planet populations. It is a well-studied process in the Orion Nebula Cluster (ONC) where the cometary morphology of proplyds is spatially resolvable due to its proximity to Earth. However, we need to study external photoevaporation in additional stellar clusters to better understand its prevalence and significance more globally. Unfortunately, more massive stellar clusters where the majority of stars form are much farther away than the ONC. In these more distant clusters the proplyds are spatially unresolvable with current facilities, hence the cometary morphology is not a useful identification of external photoevaporation. Therefore, in order to identify and interpret external photoevaporation, the only observations we have are of spatially unresolved emission lines. To resolve this issue we have used the CLOUDY code to develop an approximate general model of the emission lines emanating from the hot ionized wind of a proplyd. We have used the model to determine which line ratios are most sensitive to the distance from an OB star, and found that the most sensitive line ratios vary by multiple orders of magnitude over an FUV field of between 10$^3$ G$_0$ to 10$^6$ G$_0$. By identifying spatial gradients of line ratios in stellar clusters, we can identify regions of ongoing external photoevaporation.

Building on previous developments of a harmonic decomposition framework for computing the three-point correlation function (3PCF) of projected scalar fields over the sky, this work investigates how much cosmological information is contained in these higher-order statistics. We perform a forecast to determine the number of harmonic multipoles required to capture the full information content of the 3PCF in the context of galaxy weak lensing, finding that only the first few multipoles are sufficient to capture the additional cosmological information provided by the 3PCF. This study addresses a critical practical question: to what extent can the high-dimensional 3PCF signal be compressed without significant loss of cosmological information? Since the different multipoles contain highly redundant information, we apply a principal component analysis (PCA) which further reduces its dimensionality and preserving information. We also account for non-linear parameter degeneracies using the DALI method, an extension of Fisher forecasting that includes higher-order likelihood information. Under optimistic settings, we find that the 3PCF improves considerably the constraining power of the 2PCF for the matter abundance $\Omega_m$ and fluctuations amplitude $\sigma_8$ by about 18% for galaxy maps at redshift $z=0.5$. However, for the amplitude $S_8$ and the equation of state parameter of dark energy, $w_0$, the improvements are very modest. Thanks to the harmonic basis and PCA, the method can be implemented with relatively low computational cost, providing a practical tool for current and future photometric surveys.

The universe's large-scale structure forms a vast, interconnected network of filaments, sheets, and voids known as the cosmic web. For decades, astronomers have observed that the orientations of neighboring galaxy clusters within these elongated structures are often aligned over separations of tens of Mpc. Using the largest available catalog of galaxy clusters, we show for the first time that clusters orientations are correlated over even larger scales, up to 200-300 comoving Mpc, and such alignments are seen to redshifts of at least z = 1. Comparison with numerical simulations suggests that coherent structures on similar scales may be expected in LCDM models.

Strongly supercooled first-order phase transitions have been proposed as a primordial black hole (PBH) production mechanism. While previous works rely on simplified models with limited thermodynamic precision, we stress that reliable theoretical PBH predictions require precise nucleation dynamics within realistic, classically conformal extensions of the Standard Model. By employing high-temperature dimensional reduction and computing the one-loop fluctuation determinants, we provide a state-of-the-art thermodynamic analysis and estimate the corresponding PBH abundance for classically conformal gauge-Higgs theories. Accounting for constraints from successful percolation and QCD chiral symmetry breaking, the parameter space where PBHs are viable dark matter candidates is severely limited.

This paper describes how intentional and unintentional radio emission from airplanes is recorded with the Radio Neutrino Observatory Greenland (RNO-G). We characterize the received signals and define a procedure to extract a clean set of impulsive signals. These signals are highly suitable for instrument calibration, also for future experiments. A set of signals is used to probe the timing precision of RNO-G in-situ, which is found to match expectations. We also discuss the impact of these signals on the ability to detect neutrinos with RNO-G.

In this work, we carry out a comprehensive perturbative analysis of four cosmological models featuring a sign-switching cosmological constant. Among these, we include the well-known $\Lambda_{\rm s}$CDM model, characterised by an abrupt transition from a negative to a positive cosmological constant. We also consider the L$\Lambda$CDM model, which exhibits a generalised ladder-step evolution, as well as the SSCDM and ECDM models, both of which undergo a smooth sign change at comparable redshifts. We solve the linear cosmological perturbation equations from the radiation-dominated era, imposing initial adiabatic conditions for matter and radiation, for modes well outside the Hubble radius in the early Universe. We analyse the behaviour of the matter density contrast, the gravitational potential, the linear growth rate, the matter power spectrum, and the $f\sigma_8$ evolution . These results are contrasted with predictions from the standard $\Lambda$CDM model and are confronted with observational data.

Inspired by the amplitude bootstrap program, the spirit of this work is to constrain the space of consistent inflationary correlation functions - specifically, the bispectrum of curvature perturbations - using fundamental principles such as unitarity, locality, analyticity, and symmetries. To this end, we assume a setup for inflation in which de Sitter isometries are only mildly broken by the slow roll of the inflaton field, and study the bispectrum imprinted by a generic hidden sector during inflation. Assuming that the hidden sector's contributions to primordial non-Gaussianity are dominated by the exchange of a scalar operator (which does not preclude high-spin UV completions), we derive nontrivial positivity constraints on the resulting bispectrum $B(k_1,k_2,k_3)$. In particular, we show that $B$ must be negative in a certain region around the equilateral configuration. For instance, for isosceles triangles (with $k_2=k_3$) this region is given by $0.027\lesssim k_3/k_1\leq 2$. Furthermore, we demonstrate that unitarity imposes upper and lower bounds on the bispectrum shape, thereby carving out a Bispectrum Island where consistent shapes in our setup can reside. We complement our analysis by contemplating alternative setups where the coupling to the hidden sector is allowed to strongly break de Sitter boosts. We also identify situations that would push the bispectrum off the island and the profound physical features they would reveal.

We present a measurement of the mean number of muons with energies larger than 500 GeV in near-vertical extensive air showers initiated by cosmic rays with primary energies between 2.5 PeV and 100 PeV. The measurement is based on events detected in coincidence between the surface and in-ice detectors of the IceCube Neutrino Observatory. Air showers are recorded on the surface by IceTop, while a bundle of high-energy muons ("TeV muons") from the shower can subsequently produce a track-like event in the IceCube in-ice array. Results are obtained assuming the hadronic interaction models Sibyll 2.1, QGSJet-II.04, and EPOS-LHC. The measured number of TeV muons is found to be in agreement with predictions from air-shower simulations. The results have also been compared to a measurement of low-energy muons by IceTop, indicating an inconsistency between the predictions for low- and high-energy muons in simulations based on the EPOS-LHC model.

We extend the arguments of Maldacena and Núñez to include higher co-dimension brane setups and derive a new no-go theorem. Specifically, we show that under reasonable assumptions on the energy-momentum conservation and the bulk curvature, co-dimension-two branes fail to support stable de Sitter solutions. For co-dimensions higher than two embedded in a compact internal space, we show that negative tension sources would be required. This result places strong constraints on the viability of higher-dimensional braneworld models as a means to obtain de Sitter space within string theory.

It has recently been revealed that, in curved black-hole spacetimes, non-minimally coupled massive Proca fields may be characterized by the existence of poles in their linearized perturbation equations and may therefore develop exponentially growing instabilities. Interestingly, recent numerical computations [H. W. Chiang, S. Garcia-Saenz, and A. Sang, arXiv:2504.04779] have provided compelling evidence that the onset of monopole instabilities in the composed black-hole-field system is controlled by the dimensionless physical parameter $\mu r_-$, where $\mu$ is the proper mass of the non-minimally coupled Proca field and $r_-\equiv (-2\alpha)^{1/3}r_{\text{H}}$ is the radial location of the pole [here $\alpha$ is the non-minimal coupling parameter of the Einstein-Proca theory and $r_{\text{H}}$ is the radius of the black-hole horizon]. In the present paper we use {\it analytical} techniques in order to explore the physical properties of critical (marginally-stable) composed Schwarzschild-black-hole-nonminimally-coupled-monopole-Proca-field configurations. In particular, we derive a remarkably compact analytical formula for the discrete spectrum $\{\mu(r_{\text{H}},r_-;n) \}^{n=\infty}_{n=1}$ of Proca field masses which characterize the critical black-hole-monopole-Proca-field configurations in the dimensionless regime ${r_- -r_{\text{H}}\over{r_{\text{H}}}}\ll1$ of near-horizon poles. The physical significance of the analytically derived resonance spectrum stems from the fact that the critical field mass $\mu_{\text{c}}\equiv\mu(r_{\text{H}},r_-;n=1)$ marks the onset of instabilities in the Schwarzschild-black-hole-nonminimally-coupled-monopole-Proca-field system. In particular, composed black-hole-linearized-Proca-field configurations in the small-mass regime $\mu\leq\mu_{\text{c}}$ of the Proca field are stable.