Locally authored papers of the past 5 days

Supermassive black hole (SMBH) binaries in gaseous and stellar environments are prime targets for next-generation space-based gravitational wave detectors. Yet, realistic accretion conditions under which these binary systems evolve are not fully understood. In this work, we demonstrate the hypermagnetized multi-phase nature of the surrounding accretion flow formed by large-scale feeding from a galaxy background. Our simulations indicate that the hypermagnetized circumbinary disk is eccentric and warped, hosting a hot gas core for a parsec-scale separated binary. We also observe collimated bipolar magnetic tower-like outflows launched from each SMBH.

Large sinusoidal variations in the radio light curves of the blazars PKS J0805$-$0111 and PKS 2131$-$021 have recently been discovered with an 18-year monitoring programme by the Owens Valley Radio Observatory, making these systems strong supermassive black hole binary (SMBHB) candidates. The sinusoidal variations in PKS 2131$-$021 dominate its light curves from 2.7 GHz to optical frequencies. We report sinusoidal variations observed in both objects with the Atacama Cosmology Telescope (ACT) at 95, 147 and 225 GHz consistent with the radio light curves. The ACT 95 GHz light curve of PKS 2131$-$021 agrees well with the contemporaneous 91.5 GHz ALMA light curve and is comparable in quality. Broadband, intermittent, sinusoidal variations are also observed in PKS J0805$-$0111, for which there are no ALMA or other millimetre light curves, showing that PKS 2131$-$021 is not an isolated case and that these three properties could be common in blazar SMBHB phenomenology. In both blazars the sinusoid phase as a function of frequency as well as the achromaticity of the sinusoid amplitudes are consistent with the expected signature of jets in SMBHB systems. Monitoring of ~8000 blazars by the Simons Observatory over the next decade should provide a large number of SMBHB candidates that will shed light on the nature of the nanohertz gravitational-wave background.

Recent multi-wavelength observations of M87* \citep{2024A&A...692A.140A} revealed a high-energy $\gamma$-ray flare without a corresponding millimeter counterpart. We present a theoretical polarimetric study to evaluate the presence and nature of a potential millimeter flare in M87*, using a suite of general relativistic magnetohydrodynamical simulations with varying black hole (BH) spins and magnetic field configurations. We find that the emergence of a millimeter flare is strongly influenced by both spin and magnetic structure, with limited sensitivity to the electron distribution (thermal vs. non-thermal). We model the intensity light curve with a damped random walk (DRW) and compare the characteristic timescale ($\tau$) with recent SMA observations, finding that the simulated $\tau$ exceeds observed values by over an order of magnitude. In a flaring case with BH spin a=+0.5, we identify a distinct millimeter flare followed by an order-of-magnitude flux drop. All Stokes parameters show variability near the flare, including a sign reversal in the electric vector position angle. While most $\beta_m$ modes remain stable, the $EB$-correlation phase is highly sensitive to both the flare peak and decay. We examine polarimetric signatures in photon sub-rings, focusing on modes ns=0 and ns=1. The ns=0 signal closely matches the full image, while ns=1 reveals distinct behaviors, highlighting the potential of space VLBI to isolate sub-ring features. Finally, we analyze the magnetic and velocity field evolution during the flare, finding that magnetic reconnection weakens during the flux decay, and the clockwise velocity flow transitions into an outflow-dominated regime. These results suggest that transient radio variability near flares encodes key information about black hole spin and magnetic field structure, offering a novel probe into the physics of active galactic nuclei.

Thermochemistry, ray-tracing radiation, and radiation-matter interactions are important processes which are computationally difficult to model in astrophysical simulations, addressed by introducing novel algorithms optimized for heterogeneous architectures in the Kratos framework. Key innovations include a stoichiometry-compatible reconstruction scheme for consistent chemical species advection, which ensures element conservation while avoiding matrix inversions, and a LU decomposition method specifically designed for multi-thread parallelization in order to solve stiff thermochemical ordinary differential equations with high efficiency. The framework also implements efficient ray-tracing techniques for radiation transport for radiation-matter interactions. Various verification tests, spanning from chemical advection, combustion, Strömgren spheres, and detonation dynamics, are conducted to demonstrate the accuracy and robustness of Kratos, with results closely matching semi-analytic solutions and benchmarks such as Cantera and the Shock and Detonation Toolbox. The modular design and performance optimizations position it as a versatile tool for studying coupled microphysical processes in the diverse environments of contemporary astrophysical studies.

It is shown that the optical properties of an irregular porous grain with effective radius $a_{\rm eff}\lesssim 3\lambda$ (where $\lambda$ is the wavelength) can be well approximated by a ``spheroidal analogue'': a spheroid with appropriate axial ratio and size, with a dielectric function obtained from an effective medium theory. Prescriptions for specifying the axial ratio and porosity of the spheroidal analogue, based on simple geometric properties, are given. The accuracy of the spheroidal analogue method is studied for irregular grains with a range of structures and porosities. Different effective medium theories are compared; Bruggeman's theory is found to give the best results. The accuracy of the spheroidal analogue method justifies the use of spheroids for modeling absorption, scattering and polarization by interstellar, circumstellar, or interplanetary dust.

Time-varying dark energy is often modeled in observational analyses through generic parameterizations of its equation of state $w(z)$, which typically use two free parameters $\{w_0, w_a\}$ to span a broad range of behaviors as a function of redshift. However, this broad range of behaviors can only approximately capture the dynamics of any given microphysical theory of dark energy. A complementary approach is to use targeted parameterizations designed to model specific classes of dynamical dark energy with greater precision. Focusing on the class of thawing dark energy, we quantify and compare the precision with which nineteen generic and targeted parameterizations can capture the dynamics of physically motivated thawing quintessence theories. We find that a targeted parameterization derived from a Padé expansion of $w$ is the most reliable of these, producing accurate reconstructions of $w(z)$, the expansion history $H(z)$, and cosmological parameters such as $H_0$ and $\Omega_m$ for a broad range of microphysical theories.

There has been growing evidence that the rich star clusters in the Magellanic Clouds contain significant fractions of rapidly rotating stars. In this work, we aim to constrain these fractions by studying the colour-magnitude diagrams of four star clusters, selected among those with the most striking signatures of fast rotators. Using isochrones derived from PARSEC v2.0 stellar tracks, we generate distinct stellar populations, each covering a limited interval of initial rotation rates $\omega_\mathrm{i}$, referred to as 'Partial Models' (PMs). Using optimization algorithms and Monte Carlo Markov Chains, PMs are combined to create the final best-fitting model. In our analysis, we adopt two key assumptions: a uniform age and an isotropic distribution of stellar spin axes within each cluster. The solutions are allowed to explore the entire range of $\omega_\mathrm{i}$, and different values of age, metallicity, distance and foreground extinction. We find that the rotational velocity distributions in all four clusters reveal a high fraction of stars with $\omega_\mathrm{i}$ close to the break-up value, in all cases. Specifically, the fraction of stars with $\omega_\mathrm{i}>0.7$ exceeds $80\%$ in the clusters NGC 419 of the Small Magellanic Cloud (SMC) and NGC 1831 and NGC 1866 of the Large Magellanic Cloud (LMC). For NGC 2203 of the LMC, this fraction is smaller, although it still exceeds $50\%$, confirming that also this cluster is mainly populated by fast-rotating stars.

Pulsar halos are regions around middle-aged pulsars extending out to tens of parsecs. The large extent of the halos and well-defined central cosmic-ray accelerators make this new class of Galactic sources an ideal laboratory for studying cosmic-ray transport. LHAASO J0621+3755 is a candidate pulsar halo associated with the middle-aged gamma-ray pulsar PSR J0622+3749. We observed LHAASO J0621+3755 with VERITAS and XMM-Newton in the TeV and X-ray bands, respectively. For this work, we developed a novel background estimation technique for imaging atmospheric Cherenkov telescope observations of such extended sources. No halo emission was detected with VERITAS (0.3--10 TeV) or XMM-Newton (2--7 keV) within 1 degree and 10 arcmin around PSR J0622+3749, respectively. Combined with the LHAASO-KM2A and Fermi-LAT data, VERITAS flux upper limits establish a spectral break at ~1--10 TeV, a unique feature compared with Geminga, the most studied pulsar halo. We model the gamma-ray spectrum and LHAASO-KM2A surface brightness as inverse Compton emission and find suppressed diffusion around the pulsar, similar to Geminga. A smaller diffusion suppression zone and harder electron injection spectrum than Geminga are necessary to reproduce the spectral cutoff. A magnetic field <= 1 uG is required by our XMM-Newton observation and synchrotron spectral modeling, consistent with Geminga. Our findings support slower diffusion and lower magnetic field around pulsar halos than the Galactic averages, hinting at magnetohydrodynamic turbulence around pulsars. Additionally, we report the detection of an X-ray point source spatially coincident with PSR J0622+3749, whose periodicity is consistent with the gamma-ray spin period of 333.2 ms. The soft spectrum of this source suggests a thermal origin.

Quasi-periodic pulsations (QPPs) in the Balmer continuum of solar white-light flares (WLFs) are rarely reported, and accurately pinpointing the spatial source of flaring QPPs remains a significant challenge. We present spatiotemporal characteristics of QPPs of an X2.8 two-ribbon solar WLF (SOL2023-12-14T17:02), which was well observed by the White-light Solar Telescope (WST) aboard the Advanced Space-based Solar Observatory, with high-cadence imaging (1--2 s) in the Balmer continuum at 3600 Å. Combined with additional multi-instrument data, we find that the enhancement of the WLF in both Balmer and Paschen continua shows strong spatiotemporal correlation with hard X-ray (HXR) emissions. Notably, the pulses in the WST Balmer continuum exhibited a near-zero time lag with most HXR pulses, whereas soft X-ray and extreme ultraviolet emissions showed a lag of 2--3 s. Interestingly, quasi-harmonic QPPs with periods of $\sim$11 s and $\sim$20 s were observed in multiple wavelengths in the rising phase of the white-light continuum. Furthermore, we employed Fourier transform to spatially locate the QPPs around 11 and 20 s, revealing that they primarily originated from the east flare ribbon, which exhibited the most substantial continuum enhancement. More interestingly, we find that the west ribbon contributed significantly to the 11-second QPP but had a weaker contribution to the 20-second QPP. Moreover, the occurrence of quasi-harmonic QPPs is temporally coincident with the rapid elongation and separation motions of flare ribbons. Possible mechanisms for the quasi-harmonic QPPs have been discussed. These observations provide valuable insights into QPP modeling for solar and stellar flares.

Deuteration of hydrocarbon material, including polycyclic aromatic hydrocarbons (PAHs), has been proposed to account for the low gas-phase abundances of D in the interstellar medium. JWST spectra of four star-forming regions in M51 show an emission feature, with central wavelength $\sim$4.647$\mu$m and FWHM 0.0265$\mu$m, corresponding to the C-D stretching mode in aliphatic hydrocarbons. The emitting aliphatic material is estimated to have (D/H)$_{\rm aliph}\approx 0.17\pm0.02$ -- a factor $\sim$$10^4$ enrichment relative to the overall interstellar medium (ISM). On $\sim$$50\,$pc scales, deuteration levels toward four H$\,$II regions in M51 are 2-3 times higher than in the Orion Bar photodissociation region (PDR), with implications for the processes responsible for the formation and evolution of hydrocarbon nanoparticles, including PAHs. The deuteration of the aliphatic material is found to anticorrelate with helium ionization in the associated H$\,$II, suggesting that harsh FUV radiation may act to lower the deuteration of aliphatics in PDRs near massive stars. No evidence is found for deuteration of aromatic material, with (D/H)$_{\rm arom} \lesssim 0.016$: deuteration of the aliphatic material exceeds that of the aromatic material by at least a factor 10. The observed levels of deuteration may account for the depletion of D observed in the Galactic interstellar medium. If so, the $4.65\mu$m feature may be detectable in absorption.

We present a new dynamical measurement of the supermassive black hole mass and intrinsic shape of the stellar halo of the massive radio galaxy NGC 315 as part of the MASSIVE survey. High signal-to-noise ratio spectra from integral-field spectrographs at the Gemini and McDonald Observatories provide stellar kinematic measurements in $304$ spatial bins from the central ${\sim}0.3''$ out to $30''$. Using ${\sim} 2300$ kinematic constraints, we perform triaxial stellar orbit modeling with the TriOS code and search over ${\sim}$15,000 galaxy models with a Bayesian scheme to simultaneously measure six mass and intrinsic shape parameters. NGC 315 is triaxial and highly prolate, with middle-to-long and short-to-long axis ratios of $p=0.854$ and $q=0.833$ and a triaxiality parameter of $T=0.89$. The black hole mass inferred from our stellar kinematics is $M_\mathrm{BH} = \left(3.0 {\pm} 0.3\right) {\times} 10^{9}\ M_\odot$, which is higher than $M_\mathrm{BH}=(1.96^{+0.30}_{-0.13}) {\times} 10^{9} M_\odot$ inferred from CO kinematics (scaled to our distance). When the seven galaxies with $M_\mathrm{BH}$ measurements from both stellar and CO kinematics are compared, we find an intrinsic scatter of 0.28 dex in $M_\mathrm{BH}$ from the two tracers and do not detect statistically significant biases between the two methods in the current data. The implied black hole shadow size (${\approx} 4.7\, \mu{\rm as}$) and the relatively high millimeter flux of NGC 315 makes this galaxy a prime candidate for future horizon-size imaging studies.

We present a hybrid method for reconstructing the primordial density from late-time halos and galaxies. Our approach involves two steps: (1) apply standard Baryon Acoustic Oscillation (BAO) reconstruction to recover the large-scale features in the primordial density field and (2) train a deep learning model to learn small-scale corrections on partitioned subgrids of the full volume. At inference, this correction is then convolved across the full survey volume, enabling scaling to large survey volumes. We train our method on both mock halo catalogs and mock galaxy catalogs in both configuration and redshift space from the Quijote $1(h^{-1}\,\mathrm{Gpc})^3$ simulation suite. When evaluated on held-out simulations, our combined approach significantly improves the reconstruction cross-correlation coefficient with the true initial density field and remains robust to moderate model misspecification. Additionally, we show that models trained on $1(h^{-1}\,\mathrm{Gpc})^3$ can be applied to larger boxes--e.g., $(3h^{-1}\,\mathrm{Gpc})^3$--without retraining. Finally, we perform a Fisher analysis on our method's recovery of the BAO peak, and find that it significantly improves the error on the acoustic scale relative to standard BAO reconstruction. Ultimately, this method robustly captures nonlinearities and bias without sacrificing large-scale accuracy, and its flexibility to handle arbitrarily large volumes without escalating computational requirements makes it especially promising for large-volume surveys like DESI.

The SiTian Project represents a groundbreaking initiative in astronomy, aiming to deploy a global network of telescopes, each with a 1-meter aperture, for comprehensive time-domain sky surveys. The network's innovative architecture features multiple observational nodes, each comprising three strategically aligned telescopes equipped with filters. This design enables three-color (g, r, i) channel imaging within each node, facilitating precise and coordinated observations. As a pathfinder to the full-scale project, the Mini-SiTian Project serves as the scientific and technological validation platform, utilizing three 30-centimeter aperture telescopes to validate the methodologies and technologies planned for the broader SiTian network. This paper focuses on the development and implementation of the Master Control System (MCS),and the central command hub for the Mini-SiTian array. The MCS is designed to facilitate seamless communication with the SiTian Brain, the project's central processing and decision-making unit, while ensuring accurate task allocation, real-time status monitoring, and optimized observational workflows. The system adopts a robust architecture that separates front-end and back-end functionalities.A key innovation of the MCS is its ability to dynamically adjust observation plans in response to transient source alerts, enabling rapid and coordinated scans of target sky regions...(abridged)

The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-time inflationary physics, and underpinned by gravitation through general relativity. There have always been open questions about the soundness of the foundations of the standard model. However, recent years have shown that there may also be questions from the observational sector with the emergence of differences between certain cosmological probes. In this White Paper, we identify the key objectives that need to be addressed over the coming decade together with the core science projects that aim to meet these challenges. These discordances primarily rest on the divergence in the measurement of core cosmological parameters with varying levels of statistical confidence. These possible statistical tensions may be partially accounted for by systematics in various measurements or cosmological probes but there is also a growing indication of potential new physics beyond the standard model. After reviewing the principal probes used in the measurement of cosmological parameters, as well as potential systematics, we discuss the most promising array of potential new physics that may be observable in upcoming surveys. We also discuss the growing set of novel data analysis approaches that go beyond traditional methods to test physical models. [Abridged]

White-light observations from the WISPR instrument on NASA's Parker Solar Probe recently revealed the presence of a narrow, dense dust trail close to the orbit of asteroid 3200 Phaethon. Although Geminid-related, it aligns imperfectly with Phaethon's orbit and known Geminid meteoroid orbits. To address the nature of this dust trail, we performed a detailed comparison between the WISPR trail observations and several well-developed Geminid models. Simulating these models in the WISPR field of view visually demonstrates that the WISPR trail almost certainly represents the true ``density core'' of the Geminid stream. Trends in model trail width and offset from Phaethon's orbit, both as functions of true anomaly, agree with observations to varying extents. All the models, however, place their apparent core interior to the parent orbit due to Poynting-Robertson forces, contradictory to the WISPR trail which is exterior to Phaethon's orbit. Therefore, Phaethon's current orbit likely does not represent the orbit of the system parent, which most probably had a larger semi-major axis. These findings provide new initial conditions for future Geminid models, with WISPR identifying the Geminid core's position.

Magnetic reconnection is a key process that drives the energy release in solar flares. This process can occur at multiple locations along the coronal loop. The reconnection generates energetic electrons capable of exciting wave modes and emissions as they propagate through the loop. In this follow-up study, we investigate the influence of the injection site location of these energetic electrons - either at the looptop (LT) or at the leg of the loop around a footpoint (FP) - on the excitation of wave modes especially the second harmonic emissions (X2) in coronal loops. Our simulations reveal that the injection location significantly impacts the spatial distribution and intensity of excited wave modes. When electrons are injected at the LT, electromagnetic X2, and Z modes dominate along the loop, with minimal excitation of Langmuir waves (Yousefzadeh et al. 2021; 2022). Conversely, the present study reveals that injection close to FP leads to a strong Langmuir wave excitation throughout the loop, particularly as electrons ascend toward the LT. We find that X2 and Z modes are consistently excited at the injection site with different intensities, regardless of the injection location. However, electron injection near the FP scenario creates favorable conditions for significant Langmuir wave generation, potentially leading to plasma emission under specific circumstances. These findings emphasize the importance of electron injection location in determining the properties of the excited and emitted waves in solar coronal loops.