Papers for Tuesday, Feb 04 2025

The process by which a system of non-luminous bodies form around a star is fundamental to understanding the origins of our own solar system and how it fits into the context of other systems we have begun to study around other stars. Some basics of solar system formation have emerged to describe the process by which dust and gas around a newly formed star evolve into what we see today. The combination of occultation observations and the flyby observations by New Horizons of the Cold-Classical Kuiper Belt Object (CCKBO), (498958) Arrokoth, has provided essential new constraints on formation models through its three-dimensional shape. We present a case that an occultation-driven survey of CCKBOs would provide fundamental new insight into solar system formation processes by measuring population-wide distributions of shape, binarity, and spin-pole orientation as a function of size in this primordial and undisturbed reservoir.

Pulsar timing, i.e. the analysis of the arrival times of pulses from a pulsar, is a powerful tool in modern astrophysics. It allows us to measure the time delays of an electromagnetic signal caused by a number of physical processes as the signal propagates from the source to the observer. Joint analysis of an ensemble of pulsars (Pulsar Timing Arrays, PTAs) can be used to address a variety of astrophysical challenges, including the problem of direct detection of space-time metric perturbations, in particular those induced by gravitational waves. Here we present a comprehensive review of the current state of research in the field of pulsar timing, with particular emphasis on recent advancements in the detection of stochastic background of nHz gravitational waves, reported by a number of international collaborations such as NANOGrav (North American Nanohertz Observatory for Gravitational Waves), European Pulsar Timing Array (EPTA), Chinese Pulsar Timing Array (CPTA) and Indian Pulsar Timing Array (InPTA), which are joining their efforts within the International PTA (IPTA). Additionally, this paper reviews contemporary constraints on scalar ultralight matter (pseudoscalar bosons), obtained from timing and polarimetry data of pulsars. The prospects for applying these tools to other problems in fundamental physics and cosmology are explored.

Hyperluminous supersoft X-ray sources, such as bright extragalactic sources characterized by particularly soft X-ray spectra, offer a unique opportunity to study accretion onto supermassive black holes in extreme conditions. Examples of hyperluminous supersoft sources are tidal disruption events, systems exhibiting quasi-periodic eruptions, changing-look AGN, and anomalous nuclear transients. Although these objects are rare phenomena amongst the population of X-ray sources, we developed an efficient algorithm to identify promising candidates exploiting archival observations. In this work, we present the results of a search for hyperluminous supersoft X-ray sources in the recently released Chandra catalog of serendipitous X-ray sources. This archival search has been performed via both a manual implementation of the algorithm we developed and a novel machine-learning-based approach. This search identified a new tidal disruption event, which might have occurred in an intermediate-mass black hole. This event occurred between 2001 and 2002, making it one of the first tidal disruption events ever observed by Chandra.

Astronomical radio interferometers achieve exquisite angular resolution by cross-correlating signal from a cosmic source simultaneously observed by distant pairs of radio telescopes to produce a Fourier-type measurement called a visibility. Million Points of Light (MPoL) is a Python library supporting feed-forward modeling of interferometric visibility datasets for synthesis imaging and parametric Bayesian inference, built using the autodifferentiable machine learning framework PyTorch. Neural network components provide a rich set of modular and composable building blocks that can be used to express the physical relationships between latent model parameters and observed data following the radio interferometric measurement equation. Industry-grade optimizers make it straightforward to simultaneously solve for the synthesized image and calibration parameters using stochastic gradient descent.

The Teledyne COSMOS-66 is a next-generation CMOS camera designed for astronomical imaging, featuring a large-format sensor ($8120 \times 8120$ pixels, each $10 \mu m$), high quantum efficiency, high frame rates, and a correlated multi-sampling mode that achieves low read noise. We performed a suite of bench-top and on-sky tests to characterize this sensor and analyze its suitability for use in astronomical instruments. This paper presents measurements of linearity, conversion gain, read noise, dark current, quantum efficiency, image lag, and crosstalk. We found that the sensor exhibits nonlinear response below 5% of saturation. This nonlinearity is plausibly attributable to the trapping of electrons in each pixel. We developed and implemented a pixel-by-pixel nonlinearity correction, enabling accurate photometric measurements across the dynamic range. After implementing this correction, operating in the correlated multi-sampling mode, the sensor achieved an effective read noise of $2.9 e^-$ and dark current of $0.12 e^-/pix/s$ at $-25^\circ C$. The quantum efficiency exceeded 50% from 250 nm to 800 nm, peaking at 89% at 600 nm. We observed significant optical crosstalk between the pixels, likely caused by photoelectron diffusion. To demonstrate the sensor's astronomical performance, we mounted it on the WINTER 1m telescope at Palomar Observatory. These tests confirmed that the linearity calibration enables accurate stellar photometry and validated our measured noise levels. Overall, the COSMOS-66 delivers similar noise performance to large-format CCDs, with higher frame rates and relaxed cooling requirements. If pixel design improvements are made to mitigate the nonlinearity and crosstalk, then the camera may combine the advantages of low-noise CMOS image sensors with the integration simplicity of large-format CCDs, broadening its utility to a host of astronomical science cases.

Hypervelocity stars (HVSs) are produced by the Hills mechanism when a stellar binary is disrupted by a supermassive black hole (SMBH). The HVS Survey detected 21 unbound B-type main-sequence stars in the Milky Way's outer halo that are consistent with ejection via the Hills mechanism. We revisit the trajectories of these stars in light of proper motions from {\it Gaia} DR3 and modern constraints on the Milky Way -- Large Magellanic Cloud (LMC) orbit. We find that half of the unbound HVSs discovered by the HVS Survey trace back not the Galactic Center, but to the LMC. Motivated by this finding, we construct a forward-model for HVSs ejected from an SMBH in the LMC and observed through the selection function of the HVS Survey. The predicted spatial and kinematic distributions of simulated HVSs are remarkably similar to the observed distributions. In particular, we reproduce the conspicuous angular clustering of HVSs around the constellation Leo. This clustering occurs because HVSs from the LMC are boosted by $\sim300\,{\rm km\,s^{-1}}$ by the orbital motion of the LMC, and stars launched parallel to this motion are preferentially selected as HVS candidates. We find that the birth rate and clustering of LMC HVSs cannot be explained by supernova runaways or dynamical ejection scenarios not involving a SMBH. From the ejection velocities and relative number of Magellanic vs. Galactic HVSs, we constrain the mass of the LMC SMBH to be $10^{5.8^{+0.2}_{-0.4}} M_{\odot}$ ($\simeq 6\times10^5 M_{\odot}$).

Brightest cluster galaxies (BCGs) lie deep within the largest gravitationally bound structures in existence. Though some cluster finding techniques identify the position of the BCG and use it as the cluster center, other techniques may not automatically include these coordinates. This can make studying BCGs in such surveys difficult, forcing researchers to either adopt oversimplified algorithms or perform cumbersome visual identification. For large surveys, there is a need for a fast and reliable way of obtaining BCG coordinates. We propose machine learning to accomplish this task and train a neural network to identify positions of candidate BCGs given no more information than multiband photometric images. We use both mock observations from The Three Hundred project and real ones from the Sloan Digital Sky Survey (SDSS), and we quantify the performance. Training on simulations yields a squared correlation coefficient, R$^2$, between predictions and ground truth of R$^2 \approx 0.94$ when testing on simulations, which decreases to R$^2 \approx 0.60$ when testing on real data due to discrepancies between datasets. Limiting the application of this method to real clusters more representative of the training data, such those with a BCG r-band magnitude $r_{\text{BCG}} \leq 16.5$, yields R$^2 \approx 0.99$. The method performs well up to a redshift of at least $z\approx 0.6$. We find this technique to be a promising method to automate and accelerate the identification of BCGs in large datasets.

In the first TABASCAL paper we showed how to calibrate in the presence of Radio Frequency Interference (RFI) sources by simultaneously isolating the trajectories and signals of the RFI sources. Here we show that we can accurately remove RFI from simulated MeerKAT radio interferometry target data, for a single frequency channel, corrupted by up to 9 simultaneous satellites with average RFI amplitudes varying from weak to very strong (1 - 1000 Jy). Additionally, TABASCAL also manages to leverage the RFI signal-to-noise to phase calibrate the recovered astronomical signal. TABASCAL effectively performs a suitably phased up fringe filter for each RFI source which allows essentially perfect removal of RFI across all strengths. As a result, TABASCAL reaches image noises equivalent to the uncorrupted, no-RFI, case. Consequently, point-source science with TABASCAL almost matches the no-RFI case with near perfect completeness for all RFI amplitudes. In contrast the completeness of AOFlagger and idealised 3$\sigma$ flagging drops below 40% for strong RFI amplitudes where recovered flux errors are $\sim$10x-100x worse than those from TABASCAL. Finally we highlight that TABASCAL works for both static and varying astronomical sources.

The improved sensitivity of interferometric facilities to the 21-cm line of atomic hydrogen (HI) enables studies of its properties in galaxies beyond the local Universe. In this work, we perform a 21 cm line spectral stacking analysis combining the MIGHTEE and CHILES surveys in the COSMOS field to derive a robust HI-stellar mass relation at z=0.36. In particular, by stacking thousands of star-forming galaxies subdivided into stellar mass bins, we optimize the signal-to-noise ratio of targets and derive mean HI masses in the different stellar mass intervals for the investigated galaxy population. We combine spectra from the two surveys, estimate HI masses, and derive the scaling relation log10(MHI) = (0.32 +- 0.04)log10(M*) + (6.65 +- 0.36). Our findings indicate that galaxies at z=0.36 are HI richer than those at z=0, but HI poorer than those at z=1, with a slope consistent across redshift, suggesting that stellar mass does not significantly affect HI exchange mechanisms. We also observe a slower growth rate HI relative to the molecular gas, supporting the idea that the accretion of cold gas is slower than the rate of consumption of molecular gas to form stars. This study contributes to understanding the role of atomic gas in galaxy evolution and sets the stage for future development of the field in the upcoming SKA era.

Stellar mergers and common-envelope evolution are fast (dynamical-timescale) interactions in binary stars that drastically alter their evolution. They are key to understanding a plethora of astrophysical phenomena. Stellar mergers are thought to produce blue straggler stars, blue supergiants, and stars with peculiar rotation and surface chemical abundances. Common-envelope evolution is proposed as a key stage in the formation of gravitational wave sources, X-ray binaries, type Ia supernovae, cataclysmic variables, and other systems. A significant fraction (tens of percent) of binary stars undergo such a phase during their evolution. In this chapter, we first discuss processes leading to a stellar merger or common-envelope phase. We then explain these complex interactions, starting from underlying physical principles like entropy sorting in stellar mergers and the energy formalism in common envelopes. This is followed by a more complete picture revealed by three-dimensional (magneto)hydrodynamical simulations. The outcomes of these interactions are discussed comprehensively and special emphasis is given to the role of magnetic fields. Both stellar mergers and common-envelope evolution remain far from fully understood, and we conclude by highlighting open questions in their study.

We report the discovery of post-starburst ultra-diffuse galaxies (UDGs), identified through spectroscopic analysis with KCWI at the Keck II Telescope. Our analysis is based on a sample of 44 candidate UDGs selected from the Systematically Measuring Ultra-Diffuse Galaxies (SMUDGes) program. Our measured spectroscopic redshifts reveal $\sim 80\%$ of the entire KCWI sample exhibit large physical sizes ($R_{e} \gtrsim 1~{\rm kpc}$) and low surface brightnesses ($24 \lesssim \mu_{0,g} \lesssim 25$ mag arcsec$^{-2}$) which categorize them as UDGs. We find $20\%$ of the confirmed UDG population contain post-starburst (or K+A) features, characterized by minimal to no emission in H$\beta$ indicative of quenched star formation and a predominant presence of spectral A-type stars. Studying the local environments of the post-starburst UDGs, we find that half are isolated systems, including two systems that reside $2-3~R_{\rm vir}$ away from potential nearby massive hosts ($M_{\star} >10^{10}~\mathrm{M}_{\odot}$). Without the influence of external environmental mechanisms, these post-starburst UDGs may represent systems experiencing star formation feedback such that a recent burst may lead to (at least temporary) quenching. Overall, our results highlight the potentially diverse quenching pathways of UDGs in the local Universe.

The peculiar velocity field of the local Universe provides direct insights into its matter distribution and the underlying theory of gravity, and is essential in cosmological analyses for modelling deviations from the Hubble flow. Numerous methods have been developed to reconstruct the density and velocity fields at $z \lesssim 0.05$, typically constrained by redshift-space galaxy positions or by direct distance tracers such as the Tully-Fisher relation, the fundamental plane, or Type Ia supernovae. We introduce a validation framework to evaluate the accuracy of these reconstructions against catalogues of direct distance tracers. Our framework assesses the goodness-of-fit of each reconstruction using Bayesian evidence, residual redshift discrepancies, velocity scaling, and the need for external bulk flows. Applying this framework to a suite of reconstructions -- including those derived from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm and from linear theory -- we find that the non-linear BORG reconstruction consistently outperforms others. We highlight the utility of such a comparative approach for supernova or gravitational wave cosmological studies, where selecting an optimal peculiar velocity model is essential. Additionally, we present calibrated bulk flow curves predicted by the reconstructions and perform a density-velocity cross-correlation using a linear theory reconstruction to constrain the growth factor, yielding $S_8 = 0.69 \pm 0.034$. This result is in significant tension with Planck but agrees with other peculiar velocity studies.

TRAPPIST-1 d is generally assumed to be at the boundary between a Venus-like world and an Earth-like world, although recently published works on TRAPPIST-1 b and c raise concerns that TRAPPIST-1 d may be similarly devoid of a substantial atmosphere. TRAPPIST-1 d is also relatively understudied in comparison with TRAPPIST-1 e. The latter has generally appeared to be within the habitable zone of most atmospheric modeling studies. Assuming that TRAPPIST-1 d still retains a substantial atmosphere, we demonstrate via a series of 3D general circulation model experiments using a dynamic ocean that the planet could reside within the habitable zone in a narrow parameter space. At the same time, it could also be an exo-Venus or exo-Dead-type world or in transition between between one of these. Studies like this can help distinguish between these types of worlds.

Primordial black holes (PBH) may constitute a considerable fraction of dark matter. In this work we use the recent observations by the LIGO-Virgo-KAGRA (LVK) collaborations to set direct limits on stellar-mass range PBHs. We evaluate the merger rates of PBH binaries. Those type of interactions contribute to what is a minimum of PBH merger rates at low redshifts detectable by LVK. Thus, they allow us to derive what is the most conservative upper limits on the presence of merging PBH binaries in the gravitational-wave observations. We study the cases where PBHs have a monochromatic mass-distribution and the case case with a distribution described by a log-normal function. To derive limits on PBHs, we simulate binaries with masses following the relevant mass-distributions and merger rates as a function of redshift, up to $z\sim 0.8$, relevant for the current LVK observations. In those simulations we also take into account that the LVK observatories measure the detected black hole binaries' masses with a finite resolution. In comparing to the LVK observations, we combine the simulated PBH binaries with the binaries following a power-law mass- and redshift-distribution. The latter binaries dominate the observed population of LVK detections. Our derived limits on the mass fraction of dark matter composed of PBHs is in the range of $10^{-3}$ to O(1), depending on the exact assumptions relating to the PBH binaries properties. For reasonable assumptions on those PBH binaries' properties before their evolution inside dark matter halos, we get that fraction to be in the range of $5\times 10^{-3} - 0.1$, for PBH masses of 5-80 $M_{\odot}$. Our limits even with a conservative choice on evaluating the low-redshift merger rates provide some of the most competitive limits in the mass range of 5-40 $M_{\odot}$. [abridged]

Canonically, a protoplanetary disk is thought to undergo (gravito-)viscous evolution, wherein the angular momentum of the accreting material is transported outwards. However, several lines of reasoning suggest that the turbulent viscosity in a typical protoplanetary disk is insufficient to drive the observed accretion rates. An emerging paradigm suggests that radially extended magnetic disk winds may play a crucial role in the disk evolution. We propose a global model of magnetic wind-driven accretion for evolution of protoplanetary disks, based on the insights gained from local shearing box simulations. Here we develop this model and constrain its parameters with the help of theoretical expectations and comparison with observations. The magnetic wind is characterized with the associated loss of angular momentum and mass, which depend on the local disk conditions and stellar properties. We incorporate the disk winds self-consistently in the code FEOSAD and study formation and long-term evolution of protoplanetary disks. We include disk self-gravity and an adaptive turbulent alpha, while the co-evolution of dust is also considered. Synthetic observations are obtained via radiation thermo-chemical code ProDiMo. The models with inclusion of disk winds satisfy general expectations from both theory and observations. The disk wind parameters can be guided by observational constraints and the synthetic observations resulting from such a model compare favorably with the selected ALMA survey data of Class II disks. The proposed magnetic disk wind model is a significant step forward in the direction of representing a more complete disk evolution, wherein the disk experiences concurrent torques from viscous, gravitational, and magnetic wind processes.

In a previous paper, I described a new way of determining the high-latitude solar rotation rate statistically from space-time maps of polar faculae observed in the 6767 Å continuum by the Michelson Doppler Interferometer (MDI) on the Solar and Heliospheric Observatory (SOHO) Sheeley (2024). Now, I have tested the technique by applying it to simulated images whose faculae have known speeds, and I have been able to recover those speeds with an accuracy better than 0.01 km s$^{-1}$. Repeated measurements of the Sun's polar faculae gave the same high-latitude profile as before, but with a slightly faster synodic rotation rate of 9.$^{\circ}$10 day$^{-1}$ and a rotation period of 39.6 days. Applying this space-time tracking procedure to magnetic flux elements observed with the Helioseismic Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO), I obtained a similar rotation profile with a speed of 9.$^{\circ}$55 day$^{-1}$ and a synodic rotation period of 37.7 days. These rates are comparable to polar rotation rates, obtained by other techniques, but the new latitude profiles are noticeably flatter than the quartic fits to those prior measurements.

The next-generation mm/sub-mm/far-IR astronomy will in part be enabled by advanced digital signal processing (DSP) techniques. The Prime-Cam instrument of the Fred Young Submillimeter Telescope (FYST), featuring the largest array of submillimeter detectors to date, utilizes a novel overlap-channel polyphase synthesis filter bank (OC-PSB) for the AC biasing of detectors, implemented on a cutting-edge Xilinx Radio Frequency System on Chip (RFSoC). This design departs from traditional waveform look-up-table(LUT) in memory, allowing real-time, dynamic signal generation, enhancing usable bandwidth and dynamic range, and enabling microwave kinetic inductance detector (MKID) tracking for future readout systems. Results show that the OC-PSB upholds critical performance metrics such as signal-to-noise ratio (SNR) while offering additional benefits such as scalability. This paper will discuss DSP design, RFSoC implementation, and laboratory performance, demonstrating OC-PSB's potential in submillimeter-wave astronomy MKID readout systems.

This manual is specific to the procedure of reducing galaxy long-slit spectroscopic data from the Boller and Chivens Spectrograph aboard the 90 inch Bok telescope at Kitt Peak National Observatory. The Image Reduction Analysis Facility (IRAF) is utilized to complete the steps of image reduction. This manual will discuss methods of removing instrumental signatures, correcting for radiation events from environmental conditions, and the data frames for wavelength calibration for the particular data presented. As a result, the reduced and calibrated image with a zoomed-in figure will present the h-alpha emission line and continuum of a galaxy.

Since their first detection in 2015, gravitational wave observations have enabled a variety of studies, ranging from stellar evolution to fundamental physics. In this chapter, we focus on their use as "standard sirens", describing the different methodologies that can be adopted to measure cosmological parameters with compact object binaries from ground-based gravitational wave detectors. We cover the three main classes of standard siren measurements, showing how the expansion of the Universe can be constrained through Bayesian statistics both with gravitational wave observations alone and with the aid of electromagnetic emission from the electromagnetic counterpart of gravitational wave events and from galaxies. Finally, we summarize the existing measurements and prospects for future constraints on cosmological parameters.

A reexamination of period finding algorithms is prompted by new large area astronomical sky surveys that can identify billions of individual sources having a thousand or more observations per source. This large increase in data necessitates fast and efficient period detection algorithms. In this paper, we provide an initial description of an algorithm that is being used for detection of periodic behavior in a sample of 1.5 billion objects using light curves generated from Zwicky Transient Facility (ZTF) data (Bellm et al. 2019; Masci et al. 2018). We call this algorithm "Fast Periodicity Weighting" (FPW), derived using a Gaussian Process (GP) formalism. A major advantage of the FPW algorithm for ZTF analysis is that it is agnostic to the details of the phase-folded waveform. Periodic sources in ZTF show a wide variety of waveforms, some quite complex, including eclipsing objects, sinusoidally varying objects also exhibiting eclipses, objects with cyclotron emission at various phases, and accreting objects with complex waveforms. We describe the FPW algorithm and its application to ZTF, and provide efficient code for both CPU and GPU.

This study shows that the sputtering of hydrocarbon water ice leads to the production of mostly CO2, CO, and fragmented hydrocarbons. The onset of sputtered hydrocarbons is immediate, and quickly reaches a steady state, whereas CO2 and CO are formed more gradually. It is found that higher temperatures cause more sputtering, and that there are some notable differences in the distribution of species that are sputtered at different temperatures, indicating local heterogeneity of sputtering yields depending on the surface temperature.

The Swift Burst Alert Telescope (BAT) is a coded aperture gamma-ray instrument with a large field of view that was designed to detect and localize transient events. When a transient is detected, either on-board or externally, the BAT saves time-tagged event (TTE) data which provides the highest quality information of the locations of the photons on the detector plane and their energies. This data can be used to produce spectra, lightcurves, and sky images of a transient event. While these data products are produced by the Swift Data Center and can be produced by current software, they are often preset to certain time and energy intervals which has limited their use in the current time domain and multi-messenger environment. Here, we introduce a new capability for the BatAnalysis python package to download and process TTE data under an open-source pythonic framework that allows for easy interfacing with other python packages. The new capabilities of the BatAnalysis software allows for TTE data to be used by the community in a variety of advanced customized analyses of astrophysical sources which BAT may have TTE data for, such as Fast Radio Bursts (FRBs), Gamma-ray Bursts (GRBs), Low Mass X-ray Binaries (LMXB), Soft Gamma-ray Repeaters, magnetars, and many other sources. We highlight the usefulness of the BatAnalysis package in analyzing TTE data produced by an on-board GRB trigger, a FRB external trigger, a sub-threshold detection of the LMXB EXO 0748-676, and an external trigger of a GRB that BAT detected during a slew.

Multi-messenger observations of gravitational waves and electromagnetic emission from compact object mergers offer unique insights into the structure of neutron stars, the formation of heavy elements, and the expansion rate of the Universe. With the LIGO-Virgo-KAGRA (LVK) gravitational-wave detectors currently in their fourth observing run (O4), it is an exciting time for detecting these mergers. However, assessing whether to follow up a candidate gravitational-wave event given limited telescope time and resources is challenging; the candidate can be a false alert due to detector glitches, or may not have any detectable electromagnetic counterpart even if it is real. GWSkyNet-Multi is a deep learning model developed to facilitate follow-up decisions by providing real-time classification of candidate events, using localization information released in LVK rapid public alerts. Here we introduce GWSkyNet-Multi II, an updated model targeted towards providing more robust and informative predictions during O4 and beyond. Specifically, the model now provides normalized probability scores and associated uncertainties for each of the four corresponding source categories released by the LVK: glitch, binary black hole, neutron star-black hole, and binary neutron star. Informed by explainability studies of the original model, the updated model architecture is also significantly simplified, including replacing input images with intuitive summary values, making it more interpretable. For significant O4 event alerts issued between May 2023 and December 2024, GWSkyNet-Multi II produces a prediction that is consistent with the updated LVK classification for 93% of events. The updated model can be used by the community to help make time-critical follow-up decisions.

We present a comprehensive study of the large-scale structure, jet and outflow morphology, and kinematics of the Class 0/I protostellar binary Ced 110 IRS4, using JWST NIRCam (F150W and F410M) and MIRI MRS observations from the JWST ERC program IceAge, along with ALMA data from the Early Planet Formation in Embedded Disks (eDisk) program. NIRCam images, combined with ALMA continuum and CO data, reveal arc-like structures ($\sim$1100 au), suggesting a dense envelope around the protostars. We detect disk shadows from both protostars in F150W. The MIRI MRS IFU data reveal a jet from both protostars in multiple [Fe II] lines, [Ar II] 6.99 $\mu$m and [Ne II] 12.81 $\mu$m, marking the first detection of a jet from the system. The [Fe II] (5.34 $\mu$m) jet from Ced 110 IRS4A has a width of $\leq$ 51~au at the protostellar location, with a large opening angle of 23\arcdeg{} $\pm$ 4\arcdeg. After inclination correction, the jet velocity is 124 km s$^{-1}$, corresponding to a dynamical timescale of 25 years. The molecular H$_2$ outflow displays a distinct morphology resembling two hemispheres placed back-to-back. The consistent H$_2$ emission extent across transitions, differing from previous observations of protostellar outflows detected with JWST, suggests that MHD disk winds may not drive the observed outflow. We find that the upper limit to the width of the outflow at the protostellar location is 130 $\pm$ 10 au which is smaller than the disk diameter of 183.4 {$\pm$ 0.4 au} but much larger than the width of the [Fe II] jet.

Recently, Syderenko et al. (JCAP, 10: 018, 2016) investigated magnetogenesis and chiral asymmetry in the early hot universe. This study explores the impact of minimally coupling a constant torsion in their cosmological model, suggesting new chiral physics. Physically, this means that if torsion is right chiral, the difference between the number of right and left chiralities does not change. Moreover, the decay of chiral asymmetry depends on torsion chirality. We solve the chiral torsionful dynamo equation for magnetic field seeds. Magnetic helical fields are considered important for chiral fermion asymmetry. Even in $(3+1)$ dimensional spacetime, torsion is highly suppressed beyond inflation (Eur Phys J C 82: 291, 2022). However, torsion of $1\,\mathrm{MeV}$ appears in the early universe. Equations for correlated magnetic field coefficients are solved in terms of torsion. Weak magnetic fields of the order of $10^{-42}$ Gauss are boosted by powerful torsionful dynamo amplification, generating a much stronger magnetic field of the order of $10^{-9}$ Gauss in the present universe. A galactic magnetic field of $10^{-6}$ Gauss in the present universe, with torsion of $10^{-15}$ Gauss, leads us to a galactic dynamo seed of $10^{-9}$ Gauss. We also discuss reheating dynamo regeneration of decaying cosmic magnetic fields during the hadronization era. The relation between the reheating contribution to e-folds and the connection between CMF and temperature squared allows us to obtain dynamo amplification in terms of N-folds of inflation. The main innovation of this work is the exploration of constant torsion in a cosmological model, revealing new chiral physics. This study offers a new perspective on the origin and evolution of magnetic fields in the early universe.

Context: Stellar flares have an impact on habitable planets. To relate the observations of the Sun with those of stars, one needs to use a Sun-as-a-star analysis, that is, to degrade the resolution of the Sun to a single point. With the data of the Sun-as-a-star observations, a simulation of solar flares is required to provide a systemic clue for the Sun-as-a-star study. Aims: We aim to explore how the Sun-as-a-star spectrum varies with the flare magnitude and location based on a grid of solar flare models. Methods: Using 1D radiative hydrodynamics modeling and multi-thread flare assumption, we obtained the spectrum of a typical flare with an enhancement of chromospheric lines. Result: The Sun-as-a-star spectrum of the H$\alpha$ line shows enhanced and shifted components, which are highly dependent on the flare magnitude and location. The equivalent width $\Delta\mathrm{EW}$ is a good indicator of energy release. The bisector method can be used to diagnose the sign of the line-of-sight velocity in the flaring atmosphere. For both H$\alpha$ and H$\beta$ lines, the Sun-as-a-star spectrum of a limb flare tends to be wider and shows a dip in the line center. In particular, we propose two quantities to diagnose the magnitude and location of the stellar flares. Besides this, caution must be taken when calculating the radiation energy, since the astrophysical flux-to-energy conversion ratio is dependent on the flare location.

We present a correlation between the stellar metallicities and the mutual inclinations of multi-planet systems hosting short-period small planets (a/Rs<12, Rp<4Re). We analyzed 89 multi-planet systems discovered by Kepler, K2, and TESS, where the innermost planets have periods shorter than 10 days. We found that the mutual inclinations of the innermost two planets are higher and more diverse around metal-rich stars. The mutual inclinations are calculated as the absolute differences between the best-fit inclinations of the innermost two planets from transit modeling, which represent the lower limits of the true mutual inclinations. The mean and variance of the mutual inclination distribution of the metal-rich systems are 3.1+-0.5 and 3.1+-0.4 degrees, while for the metal-poor systems they are 1.3+-0.2 and 1.0+-0.2 degrees. This finding suggests that inner planetary systems around metal-rich stars are dynamically hotter. We summarized the theories that could plausibly explain this correlation, including the influence of giant planets, higher solid densities in protoplanetary disks around metal-rich stars, or secular chaos coupled with an excess of angular momentum deficits. Planet formation and population synthesis models tracking the mutual inclination evolution would be essential to fully understand this correlation.

This work attempts to provide a new interpretation for the hot corona in active galactic nuclei (AGNs). A thin parabolic magnetic reconnection layer, anchored at the innermost disk and extending along the boundary of the magnetic tower for a few tens of gravitational radii, serves as a hard-X-ray source above the disk. Within this reconnection layer, the tearing instability leads to the formation of a chain of plasmoids, which contain relativistic electrons that generate X-ray radiation through inverse-Compton (IC) scattering of soft photons emitted by the accretion disk. Based on previous theoretical works and numerical simulations, we develop a heuristic framework to parameterize the geometry and magnetization of the reconnection layer, as well as to compute both the power of the IC-scattering radiation and the height of the reconnection layer. Our model allows for a quantitative investigation of the relation between the height of the corona and the X-ray radiation luminosity, which can be directly compared against the observed relation from X-ray reverberation mapping of individual AGNs. The theoretical results are in good agreement with the observations of IRAS 13224-3809, indicating the validation of our model.

We derive stellar population parameters and their radial gradients within 0.65 R$_e$ for spatially resolved spectra of 2417 early-type galaxies from the MaNGA survey with stellar velocity dispersions ($\sigma$) between 50 kms$^{-1}$ and 340 kms$^{-1}$. We invert a grid of metallicity-composite stellar population models to find mean age and the abundances of C, N, Na, Mg, and Fe. These models have significant improvements compared to past models, including isochrones that respond to individual element abundances. Globally, age rises with $\sigma$ while [Fe/H] gently falls. Individual light element abundances strengthen with $\sigma$ but strengthen faster for log($\sigma$)$>$2.0. [Fe/H] shows a maximum at log($\sigma$)$\approx$2.0, falling to either side. Light element [X/Fe] anticorrelate with [Fe/H]. Heterogeneity as measured by astrophysical scatter is highest in low-$\sigma$ galaxies, most dramatically for age, Fe, and N. For galaxy-internal parameters, age shows nearly flat radial gradients in low-$\sigma$ galaxies, slightly negative at high $\sigma$. The mean radial gradient in [Fe/H] is negative and light element [X/Fe]s fall. Intrinsic scatter in gradients is highest in high-$\sigma$ galaxies, most dramatically for age and Fe. Evidently, nearly as many galaxies form inside-out as form outside-in. A near-zero radial gradient in age and light elements coupled with a mild [Fe/H] gradient supports the hierarchical merging scenario for ETG evolution. IllustrisTNG hierarchical simulations reproduce the age structure we find, show the abundance slope changes at log($\sigma$)$\approx$2.0 that we observe, and exhibit flat gradients similar to those we derive, although the abundances predicted by IllustrisTNG are significantly higher than our observations overall.

Formaldehyde is a key precursor in the formation routes of many complex organic molecules (COMs) in space. It is also an intermediate step in CO hydrogenation sequence that leads to methanol formation on the surface of interstellar grains in cold dense prestellar cores where pristine ices are formed. Various chemical models successfully reproduce the COMs abundances in cold cores, however, they consistently overpredict the abundance of formaldehyde by an order of magnitude. This results in an inverse H2CO:CH3OH abundance ratios obtained in the astrochemical simulations as compared to the observed values. In this work, we present a homogeneous data set of formaldehyde observational maps obtained towards seven dense cores in the L1495 filament with the IRAM 30 m telescope. Resolving the spatial distribution of the molecules is essential to test the chemical models. We carefully estimate the formaldehyde column densities and abundances to put reliable observational constraints on the chemical models of cold cores. Through numerous tests, we aim to constrain the updated chemical model MONACO to better align with the observed formaldehyde abundance and its ratio to methanol. In particular, we elaborate on the branching ratio of the CH3 + O reaction at low temperatures. The revised MONACO model reproduces abundances of both methanol and formaldehyde within an order of magnitude. However the model tends to overproduce formaldehyde and underpredict methanol. Consequently, the model systematically overestimates the H2CO:CH3OH ratio, although it remains within an order of magnitude of the values derived from observations.

A main source of bias in transmission spectroscopy of exoplanet atmospheres is magnetic activity of the host star in the form of stellar spots, faculae or flares. However, the fact that main-sequence stars have a chromosphere and a corona, and that these optically thin layers are dominated by line emission may alter the global interpretation of the planetary spectrum, has largely been neglected. Using a JWST NIRISS/SOSS data set of hot Jupiter HAT-P-18 b, we show that even at near-IR and IR wavelengths, the presence of these layers leads to significant changes in the transmission spectrum of the planetary atmosphere. Accounting for these stellar outer layers thus improves the atmospheric fit of HAT-P-18 b, and increases its best-fit atmospheric temperature from 536 K to 736 K, a value much closer to the predicted equilibrium temperature of 852 K. Our analysis also decreases the best-fit abundance of CO2 by almost an order of magnitude. The approach provides a new window to the properties of chromospheres/corona in stars other than our Sun.

The International Liquid Mirror Telescope (ILMT) is a 4-m aperture, zenith-pointing telescope with a field-of-view of 22', situated in the foothills of the Himalayas. The telescope operates in continuous survey mode, making it a useful instrument for time-domain astronomy, particularly for detecting transients, variable stars, active galactic nuclei variability, and asteroids. This paper presents the PyLMT transient detection pipeline to detect such transient/varying sources in the ILMT images. The pipeline utilises the image subtraction technique to compare a pair of images from the same field, identifying such sources in subtracted images with the help of convolutional neural networks (CNN) based real/bogus classifiers. The test accuracies determined for the real/bogus classifiers ranged from 94% to 98%. The resulting precision of the pipeline calculated over candidate alerts in the ILMT frames is 0.91. It also houses a CNN-aided transient candidate classifier that classifies the transient/variable candidates based on host morphology. The test accuracy of the candidate classifier is 98.6%. It has the provision to identify catalogued asteroids and other solar system objects using public databases. The median execution time of the pipeline is approximately 29 minutes per image of 17 minutes exposure. Relevant CNNs have been trained on data acquired with the ILMT during the cycle of October-November 2022. Subsequent tests on those images have confirmed the detection of numerous catalogued asteroids, variable stars, and other uncatalogued sources. The pipeline has been operational and has detected 12 extragalactic transients, including 2 new discoveries in the November 2023-May 2024 observation cycle.

The International Liquid Mirror Telescope (ILMT) project was motivated by the need for an inexpensive 4 metre diameter optical telescope that could be devoted entirely to astronomical surveys. Its scientific programmes include the detection and study of transients, variable objects, asteroids, comets, space debris and low surface brightness galaxies. To this end, a collaboration was formed between the Institute of Astrophysics and Geophysics (Liège University, Belgium), several Canadian universities (University of British Columbia, Laval University, University of Montreal, University of Toronto, York University, University of Victoria) and the Aryabhatta Research Institute of Observational Sciences (ARIES, India). After several years of design work in Belgium and construction in India on the ARIES Devasthal site, the telescope saw its first light on 29 April 2022. Its commissioning phase lasted from May 2022 until June 2023 (beginning of the monsoon). The ILMT was inaugurated on 21 March 2023 and has been in regular operation since October 2023. The telescope continuously observes the sky passing at the zenith using the SDSS g', r', and i' filters. This paper describes the ILMT, its operation, performance and shows some initial results.

The recent detections of a large number of candidate active galactic nuclei at high redshift (i.e. $z \gtrsim 4$) has increased speculation that heavy seed massive black hole formation may be a required pathway. Here we re-implement the so-called Lyman-Werner (LW) channel model of Dijkstra et al. (2014) to calculate the expected number density of massive black holes formed through this channel. We further enhance this model by extracting information relevant to the model from the $\texttt{Renaissance}$ simulation suite. $\texttt{Renaissance}$ is a high-resolution suite of simulations ideally positioned to probe the high-$z$ universe. Finally, we compare the LW-only channel against other models in the literature. We find that the LW-only channel results in a peak number density of massive black holes of approximately $\rm{10^{-4} \ cMpc^{-3}}$ at $z \sim 10$. Given the growth requirements and the duty cycle of active galactic nuclei, this means that the LW-only is likely incompatible with recent JWST measurements and can, at most, be responsible for only a small subset of high-$z$ active galactic nuclei. Other models from the literature seem therefore better positioned, at present, to explain the high frequency of massive black holes at high $z$.

(ABRIDGED) We use galaxy-galaxy lensing to investigate how the dark matter (DM) haloes and stellar content of galaxies with $0.012 \leq z \leq 0.32$ and $10 \leq \log_{10}(M_\star/\mathrm{M}_\odot) \leq 12$ change as a result of the merger process. To this end, we construct two samples of galaxies obtained from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), comprising 1 623 post-mergers and $\sim$30 000 non-merging controls, that live in low-density environments to use as our lenses. These samples are weighted to share the same distributions of stellar mass, redshift, and geometric mean distance to a galaxy's three nearest neighbours to ensure differences in the lensing signal are due to the merger process itself. We do not detect a significant difference in the excess surface density profile of post-mergers and non-merging controls with current data. Fitting haloes composed of a point-like stellar mass component and an extended DM structure described by a Navarro-Frenk-White profile to the lensing measurements yields, for both samples, halo masses of $M_\text{halo}\sim 4\times10^{12}\,\mathrm{M}_\odot$ and a moderately negative correlation between $M_\text{halo}$ and concentration $c$. . . . Extreme merger-induced starbursts, in which more than 60 percent of the stars are formed in the burst, are ruled out at the 95 per cent confidence level. Although we do not find merger-induced differences in the DM haloes and stellar content of galaxies with our current data, application of our methods to upcoming surveys that are able to provide samples $\sim$10$\times$ larger than our current catalogue are expected to detect the weak-lensing signatures of mergers and further constrain their properties.

As part of the SIGNALS survey, which comprises a sample of approximately 40 nearby galaxies observed with the Fourier transform spectrometer SITELLE, we present a study of metal mixing in the spiral galaxy NGC 6946. Taking advantage of the blue sensitivity of our setup, we measure the oxygen and nitrogen abundances of 638 H II regions, and focus our analysis on the abundance fluctuations about the radial gradients. We detect an azimuthal variation of about 0.1 dex in these abundances across the NE spiral arm, with the leading edge being more metal-poor than the trailing edge. This result aligns with galaxy simulations, where radial gas flows along the spiral arms lead to dilution on the leading edge and enrichment on the trailing edge, due to the presence of radial metallicity gradients. Our 2D analysis reveals that oxygen and nitrogen exhibit comparable spatial correlation scales, despite the different injection energies and distinct nucleosynthetic origins -- core-collapse supernovae in the case of oxygen and primarily AGB stars for nitrogen. The observed similarity suggests that stellar processes drive these two elements into the ISM over equivalent spatial scales.

We perform an MMT/Hectospec redshift survey of the North Ecliptic Pole Wide (NEPW) field covering 5.4 square degrees, and use it to estimate the photometric redshifts for the sources without spectroscopic redshifts. By combining 2572 newly measured redshifts from our survey with existing data from the literature, we create a large sample of 4421 galaxies with spectroscopic redshifts in the NEPW field. Using this sample, we estimate photometric redshifts of 77755 sources in the band-merged catalog of the NEPW field with a random forest model. The estimated photometric redshifts are generally consistent with the spectroscopic redshifts, with a dispersion of 0.028, an outlier fraction of 7.3%, and a bias of -0.01. We find that the standard deviation of the prediction from each decision tree in the random forest model can be used to infer the fraction of catastrophic outliers and the measurement uncertainties. We test various combinations of input observables, including colors and magnitude uncertainties, and find that the details of these various combinations do not change the prediction accuracy much. As a result, we provide a catalog of 77755 sources in the NEPW field, which includes both spectroscopic and photometric redshifts up to z~2. This dataset has significant legacy value for studies in the NEPW region, especially with upcoming space missions such as JWST, Euclid, and SPHEREx.

We report the discovery of ZTF J172132.75+445851.0 as a new possible VY Sculptoris variable, with an orbital period of 0.109765426(44) days. Survey observations from 2005 to the present reveal significant activity in the system, with brightness variations ranging between 18.1 and 21.3 (zr) magnitudes, including deep fades exceeding an amplitude of 3 magnitudes.

We present a systematic search for 1696 Green Pea (GP) galaxy candidates in the southern hemisphere selected from the Dark Energy Survey Data Release 2 (DES DR2) and provide preliminary results from spectroscopic follow-up observations of 26 targets chosen among them. Our selection criteria include the colors in $gri$-bands and compact morphology in the color composite images. The multi-wavelength spectral energy distribution fitting shows that the selected GP candidates exhibit star formation rates up to several tens $\mathrm{M}_\odot\,\mathrm{yr}^{-1}$. With the mean stellar mass of $\mathrm{log}\,M_*/\mathrm{M}_\odot=8.6$, GP candidates are located at roughly 1 dex above the main sequence of star-forming galaxies at $z\sim0.3$. Spectroscopic follow-up observations of the GP candidates with Gemini/GMOS are underway. All 26 targets are spectroscopically confirmed to be at $z=0.3-0.41$ and have [OIII] equivalent width larger than 85$Å$, classified to be starbursts with low to moderate dust attenuation. These confirmed GPs show a lower metallicity offset from the mass-metallicity relation of local star-forming galaxies, indicating that GPs are less chemically evolved systems at their early stage of evolution.

We reanalyze the spectral lag data for 56 GRBs in the cosmological rest frame to search for Lorentz Invariance Violation using frequentist inference. For this purpose, we use the technique of profile likelihood to deal with the nuisance parameters corresponding to a constant astrophysical spectral lag in the GRB rest frame and the unknown intrinsic scatter, while the parameter of interest is the energy scale ($E_{QG}$) for LIV. With this method, we do not obtain a global minimum for $\chi^2$ as a function of $E_{QG}$ up to the Planck scale. Therefore, we can set one-sided lower limits on $E_{QG}$ in a seamless manner. The 95% c.l. lower limits which we obtain on $E_{QG}$ are then given by $E_{QG}\geq 1.22 \times 10^{15}$ GeV and $E_{QG}\geq 6.64\times 10^{5}$ GeV, for linear and quadratic LIV, respectively. Therefore, this work represents yet another proof of principles application of profile likelihood in the search for LIV using GRB spectral lags.

Models of highly sub-Eddington accretion onto black holes commonly use a single fluid model for the collisionless, near-horizon plasma. These models must specify an equation of state. It is common to use an ideal gas with $p = (\gamma - 1) u$ and $\gamma = 4/3, 13/9,$ or $5/3$, but these produce significantly different outcomes. We discuss the origins of this discrepancy and the assumptions underlying the single fluid model. The main result of this investigation is that under conditions relevant to low luminosity black hole accretion the best choice of single fluid adiabatic index is close to but slightly less than $5/3$. Along the way we provide a simple equilibrium model for the relation between the ion-to-electron dissipation ratio and the ion to electron temperature ratio and explore the implications for electron temperature fluctuations in Event Horizon Telescope sources.

In this study, we train score-based diffusion models to super-resolve gigaparsec-scale cosmological simulations of the 21-cm signal. We examine the impact of network and training dataset size on model performance, demonstrating that a single simulation is sufficient for a model to learn the super-resolution task regardless of the initial conditions. Our best-performing model achieves pixelwise $\mathrm{RMSE}\sim0.57\ \mathrm{mK}$ and dimensionless power spectrum residuals ranging from $10^{-2}-10^{-1}\ \mathrm{mK^2}$ for $128^3$, $256^3$ and $512^3$ voxel simulation volumes at redshift $10$. The super-resolution network ultimately allows us to utilize all spatial scales covered by the SKA1-Low instrument, and could in future be employed to help constrain the astrophysics of the early Universe.

Extreme mass-ratio inspirals (EMRIs), consisting of a secondary (stellar mass) black hole (BH) orbiting around a supermassive BH, are one of the primary targets for future spaceborne gravitational wave (GW) detectors. The spin of the secondary BH encodes the formation history of the stellar mass BH and the formation process of the EMRI. In this work, we construct a kludge EMRI waveform model taking the secondary spin into account and forecast the measurement precision of the secondary spin by future spaceborne GW detectors with the Fisher information matrix. We find the secondary spin should be measured with a reasonably good precision for generic eccentric and inclined EMRIs, though the measurement precision largely degrades for low-eccentricity and nearly equatorial EMRIs. As an example of its astrophysical applications, we propose that the secondary spin can be used for distinguishing dry (loss cone) EMRIs (where the secondary BHs were born in the collapse of individual massive stars and are of low spin) and Hills EMRIs (where the secondary BHs are remnants of massive star binaries and the secondary spins follow a bimodal distribution).

The nearest spiral galaxy, M31, exhibits a kinematically hot stellar disc, a global star formation episode ~2-4 Gyr ago, and conspicuous substructures in its stellar halo, suggestive of a recent accretion event. Recent chemodynamical measurements in the M31 disc and inner halo can be used as additional constraints for N-body hydrodynamical simulations that successfully reproduce the disc age-velocity dispersion relation and star formation history, together with the morphology of the inner halo substructures. We combine an available N-body hydrodynamical simulation of a major merger (mass ratio 1:4) with a well-motivated chemical model to predict abundance distributions and gradients in the merger remnant at z=0. We computed the projected phase space and the [M/H] distributions for the substructures in the M31 inner halo, i.e. the GS, the NE-, W- Shelves. We compare these chemodynamical properties of the simulated M31 remnant with recent measurements for the M31 stars in the inner halo. This major merger model predicts (i) distinct multiple components within each of the substructure; (ii) a high mean metallicity and large spread in the GS, NE- and W- Shelves, explaining various photometric and spectroscopic metallicity measurements; (iii) simulated phase space diagrams that qualitatively reproduce various features identified in the projected phase space of the substructures in published data from the DESI; (iv) a large distance spread in the GS, as suggested by previous tip of the RGB measurements, and (v) phase space ridges caused by several wraps of the secondary, as well as up-scattered main M31 disc stars, that also have plausible counterparts in the observed phase spaces. These results provide further independent arguments for a major satellite merger in M31 ~3 Gyr ago and a coherent explanation for many of the observational results that make M31 look so different from the MW.

In this study, we investigate the effect of Poisson noise, originating from the discrete distribution of stellar-mass primordial black holes (PBHs), on the abundance of ultradense dark matter halos (UDMHs). We incorporate the contribution of PBH shot noise into the power spectrum and modify the primordial power spectrum to calculate the differential mass function of UDMHs across various models, while accounting for crucial physical and geometrical factors. We further compare our results with observational constraints on the abundance of compact dark matter and the mass distribution of PBHs derived from the thermal history of the early Universe. Our findings demonstrate that the mass of PBHs contributing to Poisson noise plays a pivotal role in determining the abundance of UDMHs. Moreover, we show that Poisson noise from lighter PBHs strengthens the single-component dark matter scenario, whereas the corresponding effect from heavier PBHs supports the multi-component dark matter scenario.

The discovery that we live in an accelerating universe changed drastically the paradigm of physics and introduced the concept of \textit{dark energy}. In this work, we present a brief historical description of the main events related to the discovery of cosmic acceleration and the basic elements of theoretical and observational aspects of dark energy. Regarding the historical perspective, we outline some of the key milestones for tracing the journey from Einstein's proposal of the cosmological constant to the type Ia supernovae results. Conversely, on the theoretical/observational side, we begin by analyzing cosmic acceleration within the context of the standard cosmological model, i.e., in terms of the cosmological constant. In this case, we show how a positive cosmological constant drives accelerated expansion and discuss the main observational aspects, such as updated results and current cosmological tensions. We also explore alternative descriptions of dark energy, encompassing dynamic and interacting dark energy models.

Active Galactic Nuclei (AGN) feedback is a key physical mechanism proposed to regulate star formation, primarily in massive galaxies. In particular, cosmic rays associated with AGN jets have the potential to efficiently suppress cooling flows and quench star formation. The locus of cosmic ray production and their coupling to gas play a crucial role in the overall self-regulation process. To investigate this in detail, we conduct high-resolution, non-cosmological MHD simulations of a massive $10^{14} {\rm M_\odot}$ halo using the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We explore a variety of AGN jet feedback scenarios with cosmic rays, examining different values for the cosmic ray energy fraction in jets, cosmic ray coupling sites (in the black hole vicinity versus at the large-scale jet-driven shock front), and jet precession parameters. Our findings indicate that when cosmic rays are injected near the black hole, they efficiently inhibit black hole accretion by suppressing the density before the jet propagates out to large radii. As a result, this leads to episodic black hole accretion, with the jet not having sufficient energy flux to reach large radii and impact cooling flows. Conversely, if the cosmic rays are injected at the jet-driven shock front, not only does the jet sustain a higher overall energy flux for an extended period, but it also disperses cosmic rays out to larger radii, more effectively suppressing the cooling flow. Furthermore, the period and angle of jet precession can influence the position of shock fronts. We identify an optimal range of jet precession periods ($\sim$ tens of Myr) that generates shocks at the inner circumgalactic medium, where cooling flows are most severe. We report that this specific configuration offers the most effective scenario for cosmic rays at the shock front to suppress the cooling flow and star formation.

We present results from a suite of 3D high-resolution hydrodynamic simulations of supernova-driven outflows from galactic disc regions with a range of gas surface density, metallicity, and supernova scale height. We use this suite to quantify how outflow properties -- particularly the loading factors for mass, metallicity, and energy -- vary with these parameters. We find that the winds fall into three broad categories: steady and hot, multiphase and moderately bursty, and cool and highly bursty. The first of these is characterised by efficient metal and energy loading but weak mass loading, the second by moderate loading of mass, metals, and energy, and the third by negligible metal and energy loading but substantial mass loading. The most important factor in determining the kind of wind a galaxy will produce is the ratio of supernova to gas gas scale heights, with the latter set by a combination of supernova rate, metallicity-dependent cooling rate, and the gravitational potential. These often combine in counterintuitive ways -- for example increased cooling causes cold clouds to sink into the galactic midplane more rapidly, lowering the volume-filling factor of dense gas and making the environment more favourable for strong winds. Our findings suggest that the nature of galactic winds is likely highly sensitive to phenomena such as runaway stars occuring at a large height and dense gas and are poorly captured in most simulations, and that metal loading factors for type Ia supernovae may be substantially larger than those for type II, with important implications for galactic chemical evolution.

We discuss the turbulent structure and dynamics of the upper solar convection zone using a 3D radiative hydrodynamic simulation model at 45 degrees latitude. The model reveals the self-formation of meridional flows, the leptocline, and the radial differential rotation. Unlike previous studies, the model shows a complex variation of the characteristic scales of turbulent flows with depth. In particular, an increase in the characteristic convective scale is trackable within an individual snapshot up to a depth of 7 Mm, near the bottom of the hydrogen ionization zone, where turbulent flows become weaker and more homogeneous. However, the turbulent spectra show an increase in scale with depth and a qualitative change in convective patterns below 7 Mm (near the bottom of the leptocline), suggesting changes in the diffusivity properties and energy exchange among different scales.

The solar magnetic field, thought to be generated by the motion of plasma within the Sun, alternates on the order of 11-year cycles and is incompletely understood. Industries rely on accurate forecasts of solar activity, but can solar cycles be predicted? And how well did predictions perform for the current cycle? The answer is less-than-stellar... Of more than 100 predictions, most underestimated the amplitude (peak sunspot number) of cycle 25, some massively so. Fewer predictions were made for the timing of solar maximum, but timing predictions seem to have performed better than amplitude predictions. Reasons for inaccurate prediction are suggested, and perspectives are given on how future studies might improve upon the extant literature.

The total lunar eclipse on March 14, 2025 UT occurs nearly exactly 521 years (one Hypersaros) after a similar eclipse on March 1, 1504 UT that is renowned for its importance to the voyage of Columbus to Jamaica. Eclipses separated by a Hypersaros have similar depths, appear very close to the same location in the sky, and occur at nearly the same time of year. This paper summarizes the results from a search for analogous cycles within the Five Millennium Catalogs of Lunar and Solar Eclipses. Under the two simple constraints of similar eclipse dates relative to the vernal equinox and similar paths of the Moon through the Earth's shadow, the most common time intervals between lunar eclipses separated by less than 1000 years are the 521-year Hypersaros and a 633-yr period of the Icosa-Inex-Triple-Saros (IITS). Notable cycles at longer periods occur at 1154, 1284, 1787, 1917, and 2308 years.

We investigate the time delay incurred during ultra-high energy cosmic ray (UHECR) propagation over cosmological distances and its potential impact on the correlation between UHECR directions of arrival and sources such as Active Galactic Nuclei (AGNs), the UHECR chemical composition, and extragalactic magnetic field constraints. We propagate particles in different magnetic field configurations, spanning over an extended range of particle Larmor radii and magnetic field coherence lengths, also including attenuation losses. We find that UHECR delays could easily be comparable to (and longer than) AGN duty cycles, effectively erasing the correlation between known AGNs and UHECR anisotropies. We finally consider how strong constraints on the chemical composition of the heaviest UHECRs could enable a better characterization of extragalactic magnetic fields.

Recent pulsar timing array (PTA) observations have detected nanohertz gravitational waves, likely originating from massive black hole binaries (MBHBs). The detected amplitude is unexpectedly higher than inferred from the electromagnetic measurements. We present new gravitational wave background (GWB) results from the ASTRID simulation. Its large volume and on-the-fly dynamical friction for MBHs provide new insights into the MBHB population, offering a more accurate assessment of its contribution to the observed GWB. ASTRID predicts a GWB from MBHBs of $h_c=2.8\times10^{-15}$, or $\sim45\%$ of the observed amplitude at $\sim 4\,{\rm nHz}$ and $h_c=2.5\times10^{-16}$ ($5\%$) with $h_c\propto f^{-1.6}$ at $\sim 30\,{\rm nHz}$. These predictions remain below current PTA constraints but align with previous empirical models based on the observed MBH mass functions. By comparison, TNG300 with post-processed MBH dynamics yields a range between $70-90\%$ ($20\% - 30\%$) of the observed levels at low (high) frequencies. At low frequencies, ASTRID predicts that the bulk of the GWB originates from MBHB with masses $M_{\rm tot}=1-3\times 10^9\,M_\odot$ peaking at $z\approx 0.3$, consistent with TNG300. Notably, both simulations predict significant GWB contribution from minor mergers ($q<0.2$) by up to $\sim 40\%$. By tracing the full merger trees of local MBHs in ASTRID, we show that they generate GWs at $\sim 10\%-80\%$ of the maximum signal assuming no accretion and recent equal-mass mergers. Finally, we demonstrate the importance of on-the-fly MBH dynamics, the lack of which leads to $3- 5$ times excessive mass growth by merger, and a similar boost to the GWB prediction.

The central regions of the globular cluster Omega Centauri ($\omega$ Cen) have been extensively studied, but its outer regions and tidal structure have been less so. Gaia's astrometry uncovered substantial tidal substructure associated with $\omega$ Cen, yet the lack of chemical tagging makes these associations tenuous. In this paper, we utilise the Gaia-Synthetic CaHK-band photometry, metallicities from the Pristine survey and Gaia's astrometry to explore up to a clustercentric radius of 5 degrees from $\omega$ Cen. We identify $\omega$ Cen-like stars based on proper motion, colour-magnitude and colour-colour space, exploring the morphology, and stellar populations of the outer regions. Our probabilistic approach recovers the tidal tails of $\omega$ Cen, and we investigate the metallicity distribution of $\omega$ Cen ranging from a radius of 15 arcmin to the tidal radius, and beyond into the tidal tails. We present (1) two components between 15 arcmin and tidal radius at -1.83 and -1.45 dex which are also the dominant populations within 15 arcmin, and (2) the first evidence that the same two populations in the outer regions of the cluster are present outside the tidal radius and into the tidal tails. These populations are mixed about the stream, and are typically amongst the faintest stars in our sample; all indicating that the tidal tails are made of tidally stripped $\omega$ Cen stars.

This study analyzes the multi-wavelength flaring activity of the distant flat spectrum radio quasar (FSRQ) OP 313 (z=0.997) during November 2023 to March 2024, using data from Fermi-Large Area Telescope, Swift X-ray Telescope, and Ultraviolet and Optical Telescope. The analysis highlights two significant very high energy(VHE) detection epochs and GeV gamma-ray ?aring episodes, providing insight into jet emission processes and radiative mechanisms. Key findings include broadband spectral energy distribution (SED) evolution, including enigmatic X-ray spectral changes. Modeling of the multi-wavelength SED with a one-zone leptonic radiative processes attributes the emissions to synchrotron radiation, Synchrotron Self-Compton (SSC), and External Compton (EC) mechanisms, with torus photons as the primary source for EC processes. The results suggest that the gamma-ray emitting region lies outside the broad-line region but within the dusty torus. Furthermore, we find that the radiated power is significantly smaller than the total jet power, suggesting that most of the bulk energy remains within the jet even after passing through the blazar emission zone. These findings advance our understanding of particle acceleration, jet dynamics, and photon field interactions in FSRQs.

Context. The HI distribution at high Galactic latitudes is found to be filamentary and closely related to the far infrared (FIR) in caustics with coherent velocity structures. These structures trace the orientation of the magnetic field lines. Aims. Recent absorption observations with the Australian SKA Pathfinder Telescope have led to major improvements in our understanding of the physical properties of the cold neutral medium (CNM) at high Galactic latitudes. We use these results to explore how far the physical state of the CNM may be related with caustics in HI and FIR. Methods. We trace filamentary FIR and HI structures and probe the absorption data for coincidences in position and velocity. Results. 57 percent of the absorption positions are associated with known FIR/HI caustics, filamentary dusty structures with a coherent velocity field. The remaining part of the absorption sample is coincident in position and velocity with genuine HI filaments that are closely related to the FIR counterparts. Thus, within the current sensitivity limitations, all of the positions with HI absorption lines are associated with filamentary structures in FIR and/or HI. We summarize physical parameters for then CNM along filaments in the framework of filament velocities that have been determined from a Hessian analysis of FIR and HI emission data. Velocity deviations between absorption components and filament velocities are due to local turbulence and we determine for the observed CNM an average turbulent velocity dispersion of 2.48 < delta_vturb < 3.9 km/s. The CNM has a mean turbulent Mach number of Mt = 3.4 +/- 1.6 km/s. Conclusions. Most, if not all, of the CNM in the diffuse ISM at high Galactic latitudes is located in filaments, identified as caustics with the Hessian operator.

K-dwarf stars are promising targets in the exploration of potentially habitable planets. Their properties, falling between G and M dwarfs, provide an optimal trade-off between the prospect of habitability and ease of detection. The KOBE experiment is a blind-search survey exploiting this niche, monitoring the radial velocity of 50 late-type K-dwarf stars. It employs the CARMENES spectrograph, with an observational strategy designed to detect planets in the habitable zone of their system. In this work, we exploit the KOBE data set to characterize planetary signals in the K7V star HIP 5957 (KOBE-1) and to constrain the planetary population within its habitable zone. We used 82 CARMENES spectra over a time span of three years. We employed a GLS periodogram to search for significant periodic signals that would be compatible with Keplerian motion on KOBE-1. We carried out a model comparison within a Bayesian framework to ensure the significance of the planetary model over alternative configurations of lower complexity. We also inspected two available TESS sectors in search of planetary signals. We identified two signals: at 8.5d and 29.7d. We confirmed their planetary nature through ruling out other non-planetary configurations. Their minimum masses are 8.80+/-0.76ME and 12.4+/-1.1ME, corresponding to absolute masses within the planetary regime at a high certainty (>99.7%). By analyzing the sensitivity of the CARMENES time series to additional signals, we discarded planets above 8.5ME within the habitable zone. We identified a single transit-like feature in TESS, whose origin is still uncertain, but still compatible within 1sigma with a transit from planet c. We have explored future prospects for characterizing this system, concluding that nulling interferometry with the LIFE mission could be capable of directly imaging both planets and characterizing their atmospheres in future studies.

Magnetic reconnection is one of the fundamental dynamical processes in the solar corona. The method of studying reconnection in active region-scale magnetic fields generally depends on non-local methods (i.e. requiring information across the magnetic field under study) of magnetic topology, such as separatrix skeletons and quasi-separatrix layers. The theory of General Magnetic Reconnection is also non-local, in that its measure of the reconnection rate depends on determining the maxima of integrals along field lines. In this work, we complement the above approaches by introducing a local theory of magnetic reconnection, that is one in which information about reconnection at a particular location depends only on quantities at that location. The theory connects the concept of the field line slippage rate, relative to ideal motion, to the underlying local geometry of the magnetic field characterized in terms of the Lorentz force and field-aligned current density. It is argued that the dominant non-ideal term for the solar corona, discussed in relation to this new theory, is mathematically equivalent to the anomalous resistivity employed by many magnetohrdrodynamic simulations. However, the general application of the theory is adaptable to the inclusion of other non-ideal terms, which may arise from turbulence modelling or the inclusion of a generalized Ohm's law. The theory is illustrated with two examples of coronal magnetic fields related to flux ropes: an analytical model and a nonlinear force-free extrapolation. In terms of the latter, the slippage rate corresponds to the reconnection which would happen if the given (static) force-free equilibrium were the instantaneous form of the magnetic field governed by an Ohm's law with non-ideal terms.

Context. Protostellar jets driven by massive protostars are collimated outflows producing high-speed shocks through dense interstellar medium. Fast shocks can accelerate particles up to relativistic energies via diffusive shock acceleration, producing non-thermal emission that can originate ${\gamma}$-ray photons. HH 80-81 is one of the most powerful collimated protostellar jets in our galaxy, with non-thermal emission detected in radio, X-ray, and ${\gamma}$-ray bands. Characterize the ${\gamma}$-ray emission originated by the accelerated particles of the region is crucial for demonstrating the capability of protostars to accelerate cosmic rays. Aims. Our goal is to determine the particle distribution that is producing the ${\gamma}$-ray spectrum of HH 80-81 in order to ascertain the leptonic/hadronic origin of the ${\gamma}$-ray emission. We aim at associating the high-energy emission in the region with the HH 80-81 system, characterize its spectrum, and elaborate emission models based on what we expect from the diffusive shock acceleration. Methods. We use the 15 yr database provided by the Fermi-LAT satellite to study the high-energy emission of the jet, spanning from 300 MeV to 100 GeV. In addition, we perform a source association based on positional arguments. Then, we employ the naima and Gamera softwares to analyze the possible mechanisms that are producing ${\gamma}$-rays considering the ambient conditions. We perform a radiative fitting and study the nature of the particles behind the ${\gamma}$-ray emission. Results. By analyzing all the candidates to produce the ${\gamma}$-ray emission that we detect, we conclude that HH 80-81 is the most probable candidate to explain the ${\gamma}$-ray emission in the region. The detected spectrum can be explained by both hadronic and leptonic particle components.

The main goal of this work was to measure the masses and radii of white dwarfs that belong to widely separated, common proper motion binaries with non-degenerate companions. These can be assessed, independently from theoretical mass-radius relations, through measurements of gravitational redshifts and photometric radii. We studied 50 white dwarfs with hydrogen-dominated atmospheres, performing a detailed analysis of high-resolution (R ~ 18,500) spectra by means of state-of-the-art grids of synthetic models and specialized software. Hence, we measured accurate radial velocities from the H-alpha and H-beta line-cores, thus obtaining the white dwarf gravitational redshifts. Jointly with a photometric analysis that is formalized by a Bayesian inference method, we measured precise white dwarf radii that allowed us to directly measure the white dwarf masses from their gravitational redshifts. The distributions of measured masses and radii agree within 6% (at the 1-sigma level) from the theoretical mass-radius relation, thus delivering a much smaller scatter in comparison with previous analyses that used gravitational redshift measurements from low-resolution spectra. A comparison against model-dependent spectroscopic estimates produces a larger scatter of 15% on the mass determinations. We find an agreement within ~10% from previous model-based, photometric mass estimates from the literature. Combining gravitational redshift measurements and photometric analysis of white dwarfs delivers precise and accurate, empirical estimates of their masses and radii. This work confirms the reliability of the theoretical mass-radius relation from the lightest to the heaviest white dwarfs in our sample (0.38-1.3 Msun). [abridged]

Dynamical tide consists of various waves that can resonate with orbital motion. We test this coupling of dynamical tide and orbital motion using a simple two-dimensional shallow water model, which can be applied to a rocky planet covered with thin ocean or atmosphere. Then we take the earth-moon system as a fiducial model to calculate the tidal resonances and orbital evolution. We find that tidal dissipation can even increase with increasing orbital separation because of the coupling of dynamical tide and orbital motion. We draw the conclusion that the coupling is not negligible to study the orbital evolution on secular timescale.

This article explores a groundbreaking update to the cosmography of the local Universe within z = 0.1, incorporating galaxy peculiar velocity datasets from the first data releases of WALLABY, FAST and DESI surveys. The galaxies with peculiar velocities currently selected in each survey is 655, 4796 and 4191 respectively. The new CF4++ compendium enables a more comprehensive study of the nearby Universe bulk flow dynamics. This analysis reveals that the dynamical scale of homogeneity is not yet reached in the interval [200-300] Mpc/h from the observer. This new data also refines the structure of local superclusters, revealing more spherical shapes and more clearly defined boundaries for key regions such as Great Attractor (Laniakea) and Coma. Very few measurements make a big difference in revealing the hidden Vela supercluster.

Pulsar Wind Nebulae (PWNe) dominate the galactic gamma-ray sky at very high energies, and are major contributors to the leptonic cosmic ray flux. However, whether or not pulsars also accelerate ions to comparable energies is not yet experimentally confirmed. We aim to constrain the birth period and pair-production multiplicity for a set of pulsars. In doing so, we aim to constrain the proportion of ions in the pulsar magnetosphere and hence the proportion of ions that could enter the pulsar wind. We estimate possible ranges of the value of the average pair production multiplicity for a sample of 26 pulsars in the Australia Telescope National Facility (ATNF) catalogue, which have also been observed by the High Energy Stereoscopic System (H.E.S.S.) telescopes. We then derive lower limits for the pulsar birth periods and average pair production multiplicities for a subset of these sources where the extent of the pulsar wind nebula and surrounding supernova shell have been measured in the radio. We also derive curves for the average pair production multiplicities as a function of birth period for sources recently observed by the Large High Altitude Air Shower Observatory (LHAASO). We show that there is a potential for hadrons entering the pulsar wind for most of the H.E.S.S. and LHAASO sources we consider, dependent upon the efficiency of luminosity conversion into particles. We also present estimates of the pulsar birth period for six of these sources, which all fall into the range of $\simeq$10-50 ms.

Recent studies have focused on how spinning black holes (BHs) within a binary system containing a strongly magnetized neutron star, then immersed in external magnetic fields, can acquire charge through mechanisms like the Wald process and how this charge could power pulsar-like electromagnetic radiation. Those objects called ``Black hole pulsar'' mimic the behaviour of a traditional pulsar, and they can generate electromagnetic fields, such as magnetic dipoles. Charged particles within an accretion disk around the black hole would then be influenced not only by the gravitational forces but also by electromagnetic forces, leading to different geometries and dynamics. In this context, we focus here on the interplay of the magnetic dipole and the accretion disk. We construct the equilibrium structures of non-conducting charged perfect fluids orbiting Kerr black holes under the influence of a dipole magnetic field aligned with the rotation axis of the BH. The dynamics of the accretion disk in such a system are shaped by a complex interplay between the non-uniform, non-Keplerian angular momentum distribution, the black hole's induced magnetic dipole, and the fluid's charge. We show how these factors jointly influence key properties of the disk, such as its geometry, aspect ratio, size, and rest mass density.

We append an additional fifteen years (2009-2024) to the Chandra X-ray light curve of M31*, the supermassive black hole at the center of M31, the Andromeda galaxy. Extending and expanding on the work in Li et al. 2011, we show that M31* has remained in an elevated X-ray state from 2006 through at least 2016 (when regular Chandra monitoring ceased) and likely through 2024, with the most recent observations still showing an elevated X-ray flux. We identify one moderate flare in 2013 where the other nuclear X-ray sources are in low-flux states, making that flare a valuable target for followup with multiwavelength and multimessenger archival data. We extract a mostly uncontaminated spectrum for M31* from this observation, showing that its X-ray properties are similar to those observed at Sgr A* in its quiescent state by Baganoff et al. 2003. Furthermore, we find no substantial change in the source's hardness ratio in the 2006 and 2013 flares compared to the post-2006 elevated state, suggesting the these flares are increases in the regular X-ray emission mechanisms instead of entirely new emission modes. Our extended light curve for M31* provides valuable context for multimessenger or multiwavelength observations of nearby supermassive black holes.

The Sun's polar magnetic field is pivotal in understanding solar dynamo processes and forecasting future solar cycles. However, direct measurements of the polar field is only available since the 1970s. The chromospheric Ca II K polar network index (PNI; the fractional area of the chromospheric network regions above a certain latitude) has recently emerged as a reliable proxy for polar magnetic fields. In this study, we derive PNI estimates from newly calibrated, rotation-corrected Ca II K observations from the Kodaikanal Solar Observatory (1904-2007) and modern data from the Rome Precision Solar Photometric Telescope (2000-2022). We use both of those Ca II K archives to identify polar network regions with an automatic adaptive threshold segmentation technique and calculate the PNI. The PNI obtained from both the archives shows a significant correlation with the measured polar field from WSO (Pearson correlation coefficient r > 0.93) and the derived polar field based on an Advective Flux Transport Model (r > 0.91). The PNI series also shows a significant correlation with faculae counts derived from Mount Wilson Observatory observations (r > 0.87) for both KoSO and Rome-PSPT data. Finally, we use the PNI series from both archives to reconstruct the polar magnetic field over a 119-year-long period, which includes last 11 solar cycles (Cycle 14-24). We also obtain a relationship between the amplitude of solar cycles (in 13-month smoothed sunspot number) and the strength of the reconstructed polar field at the preceding solar cycle minimum to validate the prediction of the ongoing solar cycle, Cycle 25.

A transiting planet was recently discovered around a star in the Taurus star-forming region, IRAS 04125+2902, making it one of the youngest known planets. The discovery paper cited two age estimates for IRAS 04125+2902, one based on a comparison to two sets of model isochrones in the Hertzsprung-Russell (H-R) diagram and a second age reported by an earlier study for a putative population in Taurus that includes IRAS 04125+2902 (D4-North). However, the model isochrones in question differ significantly for most temperatures and luminosities of young low-mass stars, and do not reproduce the observed sequences for the TW Hya and 32 Ori associations (10 and 21 Myr). Meanwhile, as found in my previous work, D4-North is a collection of field stars and fragments of several distinct Taurus groups and older associations, so its quoted age is not meaningful. The true parent population for IRAS 04125+2902 is a small group that is ~35 pc behind the L1495 and B209 clouds (B209N). I have analyzed the age of B209N through a comparison to TW Hya and 32 Ori. The M star sequences in the latter two associations have the same shapes, but the sequence for B209N is flatter, indicating that >M4 stars at ages of <10 Myr fade more quickly than stars at earlier types and older ages. For the one member of B209N that is earlier than M4 (IRAS 04125+2902), I estimate an age of 3.0+/-0.4 Myr based on its offsets from TW Hya and 32 Ori, which by happenstance is similar to the value derived through the comparison to model isochrones.

We present MauveSim, the instrument simulator software for Mauve, the latest mission from Blue Skies Space dedicated to time-domain stellar astronomy. The tool is designed to generate simulated stellar spectra, enabling the assessment of various scientific objectives, as well as determining limiting magnitudes and conducting signal-to-noise (S/N) analyses. MauveSim functions as an end-to-end simulator that takes an input stellar spectrum-either observed or synthetic-and produces a simulated observation based on the instrument's performance and characteristics. The results of MauveSim have been validated against instrument performance data from extensive ground testing campaigns, ensuring that the software reflects the most up-to-date understanding of the payload performance. Accessible to all scientists involved in the mission, MauveSim serves as a crucial tool for target selection and observation planning.

We present a new model for the TeV afterglow of GRB 221009A. The rapid increase of the TeV flux in the very early phase is reproduced by the magnetic acceleration. We consider the change in the radial structure of the circumstellar medium from homogeneous to wind-like to describe the breaks in the TeV light curve. Our results imply a highly magnetized ejecta with a significantly thick width, making the deceleration time around 400 s for observers. In our model, no early jet break is required.

Following the discovery of the brightest high-energy neutrino sources in the sky, the further detection of fainter sources is more challenging. A natural solution is to combine fainter source candidates, and instead of individual detections, aim to identify and learn about the properties of a larger population. Due to the discreteness of high-energy neutrinos, they can be detected from distant very faint sources as well, making a statistical search benefit from the combination of a large number of distant sources, a called deep-stacking. Here we show that a Bayesian framework is well-suited to carry out such statistical probes, both in terms of detection and property reconstruction. After presenting an introductory explanation to the relevant Bayesian methodology, we demonstrate its utility in parameter reconstruction in a simplified case, and in delivering superior sensitivity compared to a maximum likelihood search in a realistic simulation.

The possibility of generating a magnetic field by dynamo effect with anisotropic electrical conductivity rather than turbulent flow has been demonstrated theoretically (Plunian & Alboussière 2020) and experimentally (Alboussière et al. 2022). If the electrical conductivity is anisotropic, the electrical currents will flow preferentially in certain directions rather than others, and a simple differential rotation will suffice to generate a magnetic field. In a galaxy with spiral arms, it is reasonable to assume that the electrical conductivity will be twice larger along the arms than in the perpendicular direction, suggesting the possibility of an anisotropic dynamo. However, a further geometrical criterion must be satisfied to obtain a dynamo (Plunian & Alboussière 2022). It is given by $\Omega' \cdot\sin p > 0$, where $p$ is the pitch angle of the spiral arms, with $p \in[-\frac{\pi}{2}, \frac{\pi}{2}]$, and $\Omega'$ is the radial shear of the angular velocity. We find that the usual spiral arms galaxies, which satisfy $|\Omega'|<0$, do not satisfy this dynamo condition because they have trailing arms instead of leading arms. Even the galaxy NGC 4622, which has both trailing and leading arms, does not satisfy this dynamo condition either. This is confirmed by numerical simulations of the induction equation. Thus, for all the spiral arms galaxies known to date, the anisotropy of the spiral arms cannot explain the existence of galactic magnetic fields until further notice.

This study investigates the dispersion of magnetohydrodynamic waves influenced by thermal misbalance in a cylindrical configuration with a finite axial magnetic field within solar coronal plasmas. Specifically, it examines how thermal misbalance, characterized by two distinct timescales directly linked to the cooling and heating functions, influences the dispersion relation. This investigation is a key approach for understanding non-adiabatic effects on the behaviour of these waves. Our findings reveal that the effect of thermal misbalance on fast sausage and kink modes, consistent with previous studies on slabs, is small but slightly more pronounced than previously thought. The impact is smaller at long-wavelength limits but increases at shorter wavelengths, leading to higher damping rates. This minor effect on fast modes occurs despite the complex interaction of thermal misbalance terms within the dispersion relation, even at low-frequency limits defined by the characteristic timescales. Additionally, a very small amplification is observed, indicating a suppressed damping state for the long-wavelength fundamental fast kink mode. In contrast, slow magnetoacoustic modes are significantly affected by thermal misbalance, with the cusp frequency shifting slightly to lower values, which is significant for smaller longitudinal wavenumbers. This thermal misbalance likely accounts for the substantial attenuation observed in the propagation of slow magnetoacoustic waves within the solar atmosphere. The long-wavelength limit leads to an analytical expression that accurately describes the frequency shifts in slow modes due to misbalance, closely aligning with both numerical and observational results.

The paper considers the implementation of the fast multipole method (FMM) in the PHANTOM code for the calculation of forces in a self-gravitating system. The gravitational interaction forces are divided into short-range and long-range interactions depending on the value of the tree opening parameter of the hierarchical kd-tree. It is demonstrated that Newton's third law holds for any pair of cells of the kd-tree engaged in mutual interaction. However, for the entire system a linear momentum is not conserved. As a result, there is an unphysical force that causes the center of mass to migrate. For example, for a pair of neutron stars, the displacement of the system's center of mass is found to be comparable to the radii of the objects at times of a few tens of Keplerian revolutions. This displacement cannot be reduced by increasing the number of particles for values of the tree opening parameter greater than 0.2. For smaller values, the time required for the calculation is significantly longer.

In this work, the ability of rare VHE gamma ray selection with neural network methods is investigated in the case when cosmic radiation flux strongly prevails (ratio up to {10^4} over the gamma radiation flux from a point source). This ratio is valid for the Crab Nebula in the TeV energy range, since the Crab is a well-studied source for calibration and test of various methods and installations in gamma astronomy. The part of TAIGA experiment which includes three Imaging Atmospheric Cherenkov Telescopes observes this gamma-source too. Cherenkov telescopes obtain images of Extensive Air Showers. Hillas parameters can be used to analyse images in standard processing method, or images can be processed with convolutional neural networks. In this work we would like to describe the main steps and results obtained in the gamma/hadron separation task from the Crab Nebula with neural network methods. The results obtained are compared with standard processing method applied in the TAIGA collaboration and using Hillas parameter cuts. It is demonstrated that a signal was received at the level of higher than 5.5{\sigma} in 21 hours of Crab Nebula observations after processing the experimental data with the neural network method.

The fast and accurate estimation of planetary mass-loss rates is critical for planet population and evolution modelling. We use machine learning (ML) for fast interpolation across an existing large grid of hydrodynamic upper atmosphere models, providing mass-loss rates for any planet inside the grid boundaries with superior accuracy compared to previously published interpolation schemes. We consider an already available grid comprising about 11000 hydrodynamic upper atmosphere models for training and generate an additional grid of about 250 models for testing purposes. We develop the ML interpolation scheme (dubbed "atmospheric Mass Loss INquiry frameworK"; MLink) using a Dense Neural Network, further comparing the results with what was obtained employing classical approaches (e.g. linear interpolation and radial basis function-based regression). Finally, we study the impact of the different interpolation schemes on the evolution of a small sample of carefully selected synthetic planets. MLink provides high-quality interpolation across the entire parameter space by significantly reducing both the number of points with large interpolation errors and the maximum interpolation error compared to previously available schemes. For most cases, evolutionary tracks computed employing MLink and classical schemes lead to comparable planetary parameters at Gyr-timescales. However, particularly for planets close to the top edge of the radius gap, the difference between the predicted planetary radii at a given age of tracks obtained employing MLink and classical interpolation schemes can exceed the typical observational uncertainties. Machine learning can be successfully used to estimate atmospheric mass-loss rates from model grids paving the way to explore future larger and more complex grids of models computed accounting for more physical processes.

Atmospheric characterization of Earth-like exoplanets through reflected light spectroscopy is a key goal for upcoming direct imaging missions. A critical challenge in this endeavor is the accurate determination of planetary mass, which may influence the measurement of atmospheric compositions and the identification of potential biosignatures. In this study, we used the Bayesian retrieval framework ExoReL$^\Re$ to investigate the impact of planetary mass uncertainties on the atmospheric characterization of terrestrial exoplanets observed in reflected light. Our results indicate that precise prior knowledge of the planetary mass can be crucial for accurate atmospheric retrievals if clouds are present in the atmosphere. When the planetary mass is known within 10\% uncertainty, our retrievals successfully identified the background atmospheric gas and accurately constrained atmospheric parameters together with clouds. However, with less constrained or unknown planetary mass, we observed significant biases, particularly in the misidentification of the dominant atmospheric gas. For instance, the dominant gas was incorrectly identified as oxygen for a modern-Earth-like planet or carbon dioxide for an Archean-Earth-like planet, potentially leading to erroneous assessments of planetary habitability and biosignatures. These biases arise because, the uncertainties in planetary mass affect the determination of surface gravity and atmospheric scale height, leading the retrieval algorithm to compensate by adjusting the atmospheric composition. Our findings emphasize the importance of achieving precise mass measurements-ideally within 10\% uncertainty-through methods such as extreme precision radial velocity or astrometry, especially for future missions like the Habitable Worlds Observatory.

The nature of dark matter, which constitutes 27% of the Universe's matter-energy content, remains one of the most challenging open questions in physics. Over the past two decades, the DAMA/LIBRA experiment has reported an annual modulation in the detection rate of $\approx$250 kg of NaI(Tl) detectors operated at the Gran Sasso Laboratory, which the collaboration interprets as evidence of the galactic dark matter detection. However, this claim has not been independently confirmed and is refuted under certain dark matter particle and halo model scenarios. Therefore, it is crucial to perform an experiment with the same target material. The ANAIS experiment uses 112.5 kg of NaI(Tl) detectors at the Canfranc Underground Laboratory and it has been collecting data since August 2017 to model-independently test the DAMA/LIBRA result. This article presents the results of the annual modulation analysis corresponding to six years of ANAIS-112 data. Our results, the most sensitive to date with the same target material, NaI(Tl), are incompatible with the DAMA/LIBRA modulation signal at a 4$\sigma$ confidence level. Such a discrepancy strongly challenges the DAMA/LIBRA dark matter interpretation and highlights the need to address systematic uncertainties affecting the comparison, particularly those related to the response of detectors to nuclear recoils, which may require further characterization of the DAMA crystals.

V498 Hya (SDSS J084555.07+033929.2) was identified as a short-period cataclysmic variable (CV) by the Catalina Real-Time Transient Survey (CRTS) in 2008. The superhump period was measured during the detected single superoutburst of V498 Hya. The quiescent spectrum subsequently taken by the \SDSSV\ Milky Way Mapper survey suggested that the CV donor may be a brown dwarf. We present time-resolved follow-up spectroscopy of V498 Hya in quiescence, obtained with the GTC OSIRIS spectrograph, from which we derived the 86.053 min spectroscopic period, systemic radial velocity, and the gravitational redshift of the Mg II line. We also modeled the spectral energy distribution to constrain the system parameters, including the > 0.82 Ms mass of the white dwarf and the best-fit value 0.043 +/- 0.004 Ms of the donor star mass. This combination of parameters implies that V498 Hya has evolved past the period minimum and is a relatively rare ``period bouncer''.

Type Ia supernovae are the established `standard candle' in the construction of the Hubble diagram out to high luminosity distances. Since the Hubble constant that best fits observations of these supernovae often turns out to be high compared to fits to other data, they are currently being investigated for possible systematic effects, with many studies focusing on the calibration of the distance ladder in the local Universe. Here we present a simulation-based assessment of another type of systematic effect, related to the chance that the line of sight to a distant supernova passes close to a foreground galaxy. We consider two cases separately: First, the foreground galaxy may block the line of sight so that the supernova is not observed. Since foreground galaxies are correlated with overdensities that typically magnify the flux of background sources, this effect leads to a systematic removal of lensed supernovae from the sample, biasing the high-redshift Hubble diagram towards demagnified (fainter) supernovae. Second, if the supernova can be observed, its proximity to the foreground galaxy can lead to an incorrect host assignment, especially if the true host has a low surface brightness. Since foreground galaxies are typically found at lower redshifts, this effect introduces another systematic bias. The probability of line-of-sight alignments with foreground galaxies increases with redshift and therefore affects distant supernovae more strongly. We find that both effects are small, but the effect of host misidentification should be included in the systematic error budget at current levels of measurement precision.

High-resolution spectroscopy has provided a wealth of information about the climate and composition of ultra-hot Jupiters. However, the 3D structure of their atmospheres makes observations more challenging to interpret, necessitating 3D forward-modeling studies. In this work, we model phase-dependent thermal emission spectra of the archetype ultra-hot Jupiter WASP-76b to understand how the line strengths and Doppler shifts of Fe, CO, H$_2$O, and OH evolve throughout the orbit. We post-process outputs of the SPARC/MITgcm global circulation model with the 3D Monte-Carlo radiative transfer code gCMCRT to simulate emission spectra at 36 orbital phases. We then cross-correlate the spectra with different templates to obtain CCF and $K_{\text{p}}$$-$$V_{\text{sys}}$ maps. For each species, our models produce consistently negative $K_{\text{p}}$ offsets in pre- and post-eclipse, which are driven by planet rotation. The size of these offsets is similar to the equatorial rotation velocity of the planet. Furthermore, we demonstrate how the weak vertical temperature gradient on the nightside of ultra-hot Jupiters mutes the absorption features of CO and H$_2$O, which significantly hampers their detectability in pre- and post-transit. We also show that the $K_{\text{p}}$ and $V_{\text{sys}}$ offsets in pre- and post-transit are not always a measure for the line-of-sight velocities in the atmosphere. This is because the cross-correlation signal is a blend of dayside emission and nightside absorption features. Finally, we highlight that the observational uncertainty in the known orbital velocity of ultra-hot Jupiters can be multiple km/s, which makes it hard for certain targets to meaningfully report absolute $K_{\text{p}}$ offsets.

We present deep ALMA observations targeting the [CII]$_{158\mu m}\,$ line in JADES-GS-z14-0, the most distant known galaxy at z=14.1793. We do not detect the [CII]$_{158\mu m}\,$ line in our deep observations, implying a luminosity of $<$6$\times10^7$ L$_{\odot}$ (3$\sigma$) for the target. Comparing this with the detected [OIII]$_{88\mu m}\,$ line, we constrain the [OIII]/[CII] ratio to be $>$3.5, significantly improving our probe of the ionization parameter $U$. The observed ratio is higher than analogues in the local universe, but consistent with galaxies at $z\approx6$-9. Through ISM modeling, we infer extreme ionizing conditions with log(U)$>-$2.0, likely requiring a young stellar population. Our modeling also indicates a relatively low gas density ($51_{-32}^{+116}$ cm$^{-3}$), significantly lower than expected from lower redshift trends. We infer a relatively high gas-phase metallicity (16$\pm$6\% solar) consistent with previous results and implying a rapid build-up of metals. Finally, using [CII]$_{158\mu m}\,$ as a molecular/cold gas mass tracer, we infer a low gas fraction ($f_\mathrm{gas} < 0.77$), consistent with previous estimates of the dynamical mass from [OIII]$_{88\mu m}\,$. Combined with the low observed gas density, lack of dust and high ionization parameter, this suggests strong feedback processes are playing an important role in the evolution of this galaxy. Our observations show that JADES-GS-z14-0 is a rapidly evolving galaxy with extreme ISM conditions, shedding light on the earliest phases of galaxy formation.

We show that in realistic models where primordial black holes are formed due to the collapse of sizeable inflationary perturbations, their initial spatial clustering beyond Poisson distribution does not play any role in the binary mergers, including sub-solar primordial black holes, responsible for the gravitational waves detectable by LIGO-Virgo-KAGRA. This is a consequence of the existing FIRAS CMB distortion constraints on the relevant scales. This conclusion might not hold for lighter masses potentially accessible by future gravitational wave observations.

We use the IllustrisTNG cosmological hydrodynamical simulation to study the rotation curves of galaxies in the local universe. To do that, we first select the galaxies with 9.4 $<$ $\log{(M_\mathrm{star}/M_\odot)}$ $<$ 11.5 to make a sample comparable to that of SDSS/MaNGA observations. We then construct the two-dimensional line-of-sight velocity map and conduct the fit to determine the rotational velocity and the slope of the rotation curve in the outer region ($R_\mathrm{t}<r<3\times r_\mathrm{half,*}$). The outer slopes of the simulated galaxies show diverse patterns that are dependent on morphology and stellar mass. The outer slope increases as galaxies are more disky, and decreases as galaxies are more massive, except for the very massive early-type galaxies. The outer slope of the rotation curves shows a correlation with the dark matter fraction, slightly better than for the gas mass fraction. Our study demonstrates that the observed dependence of galaxy rotation curves on morphology and stellar mass can be successfully reproduced in cosmological simulations, and provides a hint that dark matter plays an important role in shaping the rotation curve. The sample of simulated galaxies in this study could serve as an important testbed for the subsequent study tracing galaxies back in time, enabling a deeper understanding of the physical origin behind the diverse rotation curves.

X-ray observing facilities, such as the Chandra X-ray Observatory and the eROSITA, have detected millions of astronomical sources associated with high-energy phenomena. The arrival of photons as a function of time follows a Poisson process and can vary by orders-of-magnitude, presenting obstacles for common tasks such as source classification, physical property derivation, and anomaly detection. Previous work has either failed to directly capture the Poisson nature of the data or only focuses on Poisson rate function reconstruction. In this work, we present Poisson Process AutoDecoder (PPAD). PPAD is a neural field decoder that maps fixed-length latent features to continuous Poisson rate functions across energy band and time via unsupervised learning. PPAD reconstructs the rate function and yields a representation at the same time. We demonstrate the efficacy of PPAD via reconstruction, regression, classification and anomaly detection experiments using the Chandra Source Catalog.

We demonstrate that in the presence of a light scalar spectator field, vacuum transitions taking place during inflation can produce large, potentially detectable non-Gaussian signatures in the primordial curvature perturbation. Such transitions are common in theories with multiple scalar fields when the potential has several minima. Our computation proceeds by numerically finding the instanton solution that describes quantum tunnelling between vacuum states in a de Sitter background, calculating its dependence on the spectator field and, thereby, its effect on the expansion of space. For a scenario with Higgs inflation, we obtain the non-Gaussianity parameter $f_\mathrm{NL} \sim O(10)$ and study its parameter dependence.

Nearby extremely metal-poor galaxies (XMPs) allow us to study primitive galaxy formation and evolution in greater detail than is possible at high redshift. This work, for the first time, promotes the use of convolutional neural networks (CNNs) to efficiently search for XMPs in multi-band imaging data based on their predicted N2 index (N2\,$\equiv\log$\{\rNii/\Ha\}). We developed a sequential characterisation pipeline, composed of three CNN procedures: (i) a classifier for metal-poor galaxies, (ii) a classifier for XMPs, and (iii) an N2 predictor. The pipeline is applied to over 7.7 million SDSS DR17 imaging data without SDSS spectroscopy. The predicted N2 values are used to select promising candidates for observations. This approach was validated by new observations of 45 candidates with redshifts less than 0.065 using the 2.54~m Isaac Newton Telescope (INT) and the 4.1~m Southern Astrophysical Research (SOAR) Telescope between 2023 and 2024. All 45 candidates are confirmed to be metal-poor, including 28 new discoveries. There are 18/45 galaxies lacking detectable \rNii\ lines ($S/N<2$); for these, we report $2\sigma$ upper limits on their oxygen abundance. Our XMPs have estimated oxygen abundances of $7.1\leq$\OH$\leq8.7$ ($2\sigma$ upper limit), based on the N2 index, and 21 of them with estimated metallicity $<0.1~Z_{\odot}$. Additionally, we identified 4 potential candidates of low-metallicity AGNs at $\lesssim0.1Z_{\odot}$. Finally, we found that our observed samples are mostly brighter in the $g-$band compared to other filters, similar to blueberry (BB) galaxies, resembling green pea galaxies and high-redshift Ly$\alpha$ emitters.

CubeSat technology is an emerging alternative to large-scale space telescopes due to its short development time and cost-effectiveness. MeVCube is a proposed CubeSat mission to study the least explored MeV gamma-ray sky, also known as the `MeV gap'. Besides being sensitive to a plethora of astrophysical phenomena, MeVCube can also be important in the hunt for dark matter. If dark matter is made up of evaporating primordial black holes, then it can produce photons in the sensitivity range of MeVCube. Besides, particle dark matter can also decay or annihilate to produce final state gamma-ray photons. We perform the first comprehensive study of dark matter discovery potential of a near-future MeVCube CubeSat mission. In all cases, we find that MeVCube will have much better discovery reach compared to existing limits in the parameter space. This may be an important step towards discovering dark matter through its non-gravitational interactions.

We establish the general conditions under which evolution in the laws of physics and matter creation or destruction are closely intertwined. They make use of global time variables canonically dual to the constants of Nature. Such times flow at a rate determined by what can be interpret as the chemical potential of the fundamental constants (in analogy with phenomenological clocks based on isentropic fluids). The general condition for violations of energy conservation is then that a matter parameter evolves as a function of a gravity clock or vice-versa. This framework can be envisaged as the environment within which a natural selection scenario operates, powered by random mutations in the values of the constants of nature (or indeed any other variability in the laws in terms of the times defined above). The prize function is the creation of matter, followed by its preservation. This can be accomplished in an environment where diffeomorphism invariance is among the possible theories, with mutations modelled, for example, on the absorbing Markov chain. In such a set-up the diffeormorphism invariant state with fixed constants (or any nearby state) should be the absorbing state. John Wheeler's ``higgledy-piggledy'' chaotic cosmic start therefore finds a realization in this model, where its own demise and the establishment of order and seemingly immutable laws is also a predection of the model.

We revisit the Fractional Holographic Dark Energy (FHDE) model to reconstruct it by means of dynamic candidates such as ($i$) Quintessence, ($ii$) K-essence, ($iii$) Dilaton, ($iv$) Yang-Mills condensate, ($v$) DBI-essence, and ($vi$) Tachyonic fields in a flat Friedmann-Robertson-Walker (FRW) Universe. In particular, the dark-energy possibilities ($i$)-($vi$) are formulated through suitable field descriptions. Being concrete, we establish a comprehensive correspondence between FHDE and suitable scalar and gauge field frameworks that co-substantiate our investigation and subsequent discussion. In more detail, we methodically compute the corresponding Equation of State (EoS) parameters and field (kinetic and potential) features for the fractional parameter ($\alpha$) range, viz. $1<\alpha\leq2$. Conclusively, our results show that the modifications brought by the fractional features satisfactorily enable late-time cosmic acceleration, together with avoiding quantum instabilities by preventing the EoS from entering the phantom divide i.e., $\omega(z)\rightarrow-\infty$, which is a common issue in standard scalar field models without fractional dynamics (e.g., K-essence field). Our findings further indicate that fractional calculus attributes can be significant in addressing the challenges of dark-energy models by offering a robust framework to prospect late-time acceleration and properly fitting observational constraints. Notably, we find that as the fractional features start to dominate, the EoS parameter of all the effective field configurations asymptotically approaches a $\Lambda$CDM behaviour in the far-future limit $z\rightarrow-1$. In summary, the recent perspective introduced by FHDE \citep{Trivedi:2024inb} can indeed be cast as a promising aspirant through the use of prominent field frameworks.

In this paper, we investigate static spherically symmetric solutions in the context of Conformal Killing Gravity, a recently proposed modified theory of gravity that offers a new approach to the cosmological constant problem. Coupling this new theory with nonlinear electrodynamics, we derive the corresponding field equations and study their behavior under different parameter choices. We analyze three different models, each focusing on different key parameters. Our results reveal a rich causal structure with multiple horizons and transitions between extreme and non-extreme solutions depending on the parameter values. Moreover, we compute the nonlinear Lagrangian density for each model and find that it agrees with Maxwell theory in the limit $F \rightarrow 0$. We also confirm the existence of a central curvature singularity via the Kretschmann scalar. To connect our theoretical results with observational prospects, we study the black hole shadows associated with each model. The analysis shows that the calculated shadow size and shape of the three proposed models are consistent with the data for the supermassive object at the center of our galaxy and are therefore possible candidates for modeling this structure.

The one-loop corrections (1LC) to the power spectrum of scalar perturbations, arising from cubic interactions in the single-field E-type $\alpha$-attractor models of inflation and primordial black hole (PBH) production, are numerically calculated. The results demonstrate the 1LC contributes merely a few percent to the tree-level power spectrum. The model parameters are chosen to predict the PBH masses in the asteroid-mass range, while maintaining consistency with the cosmic microwave background (CMB) observations within 1$\sigma$ confidence levels, and obeying the upper limits on $\mu$-distortions. The PBH formed on scales smaller than the inflation scale can constitute a significant fraction of the present dark matter (DM). The PBH-induced gravitational waves (GW) may be detectable by the future space-based gravitational interferometers. We also consider a reconstruction of the scalar potential from possible GW observations and present a numerical approach tested in the parameter space of the model.

We present novel implementations of Starobisky-like inflation within Supergravity adopting Kahler potentials for the inflaton which parameterize hyperbolic geometries known from the T-model inflation. The associated superpotentials are consistent with an R and a global or gauge U(1)_X symmetries. The inflaton is represented by a gauge singlet or non-singlet superfield and is accompanied by a gauge-singlet superfield successfully stabilized thanks to its compact contribution into the total Kahler potential. Keeping the Kahler manifold intact, a conveniently violated shift symmetry is introduced which allows for a slight variation of the predictions of Starobinsky inflation: The (scalar) spectral index exhibits an upper bound which lies close to its central observational value whereas the constant scalar curvature of the inflaton-sector Kahler manifold increases with the tensor-to-scalar ratio.

The Cold Dark Matter (CDM) hypothesis accurately predicts large-scale structure formation and fits the Cosmic Microwave Background temperature fluctuations (CMB). However, observations of the inner regions of dark matter halos and dwarf galaxy satellites have consistently posed challenges to CDM. On the other hand, the Modified Newtonian Dynamics (MOND) hypothesis can explain galactic phenomena but fails to account for the complex shape of the CMB and matter power spectra. CDM and MOND are effective in nearly mutually exclusive regimes, prompting the question: Is there a physical mechanism where CDM and MOND share a common origin? Q-balls, which are localized, non-topological solitons, can be a bridge between the two hypotheses. Q-balls formed in the early Universe can mimic CDM at cosmological scales. Interestingly, Q-balls can exhibit MOND-like behavior in the late Universe at galactic scales, providing a unified framework. Specifically, we demonstrate that millicharged composite Q-balls formed from complex scalar fields, decoupled from the background radiation, can naturally arise during the radiation-dominated epoch. From the matter-radiation equality, we also obtain the mass of Q-balls to be $1~{eV}$, which are much smaller than the electron mass. Using the constraints from the invisible decay mode of ortho-positronium, we obtain $Q < 3.4 \times 10^{-5}$. We also establish an upper bound on the number density of Q-balls, which depends on the charge of the Q-ball and the small initial charge asymmetry. Furthermore, we demonstrate that the MOND naturally emerges at the galactic scale within the framework of our Q-ball model.

Gravitational waves (GWs) emanated by stellar mass compact binary coalescences (CBCs), and lensed by galaxy- or cluster-scale lenses, will produce two or more copies of the GW signal. These will have identical phase evolution but differing amplitudes. Such lensing signatures are expected to be detected by the end of the LIGO-Virgo-Kagra's (LVK's) fifth observing run (O5). In this work, we propose a novel $\chi_{\mathrm{lens}}^2$ statistic to segregate pairs of detected GW events as either lensed or unlensed, using templates typically used in GW searches. The statistic is an application of the generalized $\chi^2$ discriminator described in \citet{dhurandhar2017}, tailored to probe the similarity (or lack thereof) between the phase evolutions of two CBC signals. We assess the performance of $\chi_{\mathrm{lens}}^2$ on a realistic astrophysical dataset of lensed and unlensed CBCs detectable in O4, assuming a single LIGO-like detector at design sensitivity. We find that we can correctly identify lensed events with efficiencies comparable to existing Bayesian and machine learning methods. Evaluating $\chi_{\mathrm{lens}}^2$ is orders of magnitude faster than Bayesian methods. Moreover, the statistics of $\chi_{\mathrm{lens}}^2$, in stationary Gaussian noise, are fully understood, in contrast to machine learning methods. $\chi_{\mathrm{lens}}^2$ can, therefore, be used to rapidly and accurately weed out the vast majority of unlensed candidate pairs and identify lensed pairs.

We present the Minimally-Implicit Runge-Kutta (MIRK) methods for the numerical evolution of the resistive relativistic magnetohydrodynamic (RRMHD) equations, following the approach proposed by Komissarov (2007) of an augmented system of evolution equations to numerically deal with constraints. Previous approaches rely on Implicit-Explicit (IMEX) Runge-Kutta schemes; in general, compared to explicit schemes, IMEX methods need to apply the recovery (which can be very expensive computationally) of the primitive variables from the conserved ones in numerous additional times. Moreover, the use of an iterative process for the recovery could have potential convergence problems, increased by the additional number of required loops. In addition, the computational cost of the previous IMEX approach in comparison with the standard explicit methods is much higher. The MIRK methods are able to deal with stiff terms producing stable numerical evolutions, minimize the number of recoveries needed in comparison with IMEX methods, their computational cost is similar to the standard explicit methods and can actually be easily implemented in numerical codes which previously used explicit schemes. Two standard numerical tests are shown in the manuscript.

Gravitational-wave (GW) ringdown signals from black holes (BHs) encode crucial information about the gravitational dynamics in the strong-field regime, which offers unique insights into BH properties. In the future, the improving sensitivity of GW detectors is to enable the extraction of multiple quasi-normal modes (QNMs) from ringdown signals. However, incorporating multiple modes drastically enlarges the parameter space, posing computational challenges to data analysis. Inspired by the $F$-statistic method in the continuous GW searches, we develope an algorithm, dubbed as FIREFLY, for accelerating the ringdown signal analysis. FIREFLY analytically marginalizes the amplitude and phase parameters of QNMs to reduce the computational cost and speed up the full-parameter inference from hours to minutes, while achieving consistent posterior and evidence. The acceleration becomes more significant when more QNMs are considered. Rigorously based on the principle of Bayesian inference and importance sampling, our method is statistically interpretable, flexible in prior choice, and compatible with various advanced sampling techniques, providing a new perspective for accelerating future GW data analysis.

We investigated the generation of the $\alpha$ and $\beta$ effects in a rotating spherical plasma system with oppositely polarized kinetic helicity in the northern and southern hemispheres and examined their contributions to the induction of magnetic fields. We found that the $\alpha$ effect is relatively small, and its sign depends on the polarization of kinetic helicity. In contrast, the $\beta$ effect remains negative regardless of the sign of kinetic helicity. Despite its small magnitude, the $\alpha$ effect plays a crucial role in determining the polarity of helical magnetic structures, while a negative $\beta$ indicates energy diffusion from turbulent regions into the large-scale magnetic field. We derived the $\alpha$ and $\beta$ effects with oppositely polarized kinetic helicity using different approaches, incorporating large-scale magnetic data and turbulent kinetic data. These were used to reproduce the large-scale magnetic field and compare it with DNS results. In the kinematic regime, where the magnetic field strength is weak, our results align well; however, in regions with strong nonlinear magnetic effects, the magnetic field reproduced using turbulent kinetic data diverges. This divergence is attributed to insufficient quenching of the $\beta$ effect, suggesting that including the second-moment terms of velocity in the magnetic field effect would improve the accuracy of the $\beta$ coefficient. In this study, we considered the case of a rotating plasma sphere with $Pr_M = 1$ and low Reynolds numbers. However, in reality, Reynolds numbers are much higher, and $Pr_M$ is much less than 1, which necessitates further studies on this topic. We plan to address this in future research.

The match, which is defined as the the similarity between two waveform templates, is a fundamental calculation in computationally expensive gravitational-wave data-analysis pipelines, such as template bank generation. In this paper we introduce LearningMatch, a Siamese neural network that has learned the mapping between the parameters, specifically $\lambda_{0}$ (which is proportional to the chirp mass), $\eta$ (symmetric mass ratio), and equal aligned spin ($\chi_{1}$ = $\chi_{2}$), of two gravitational-wave templates and the match. The trained Siamese neural network, called LearningMatch, can predict the match to within $3.3\%$ of the actual match value. For match values greater than 0.95, a trained LearningMatch model can predict the match to within $1\%$ of the actual match value. LearningMatch can predict the match in 20 $\mu$s (mean maximum value) with Graphical Processing Units (GPUs). LearningMatch is 3 orders of magnitudes faster at determining the match than current standard mathematical calculations that involve the template being generated.