Papers for Tuesday, May 23 2023

High-redshift quasars ($z\gtrsim6$), powered by black holes (BHs) with large inferred masses, imply rapid BH growth in the early Universe. The most extreme examples have inferred masses of $\sim \! 10^9\,$M$_\odot$ at $z = 7.5$ and $\sim \! 10^{10}\,$M$_\odot$ at $z = 6.3$. Such dramatic growth via gas accretion likely leads to significant energy input into the quasar host galaxy and its surroundings, however few theoretical predictions of the impact of such objects currently exist. We present zoom-in simulations of a massive high-redshift protocluster, with our fiducial FABLE model incapable of reproducing the brightest quasars. With modifications to this model to promote early BH growth, such as earlier seeding and mildly super-Eddington accretion, such `gargantuan' BHs can be formed. With this new model, simulated host dust masses and star formation rates are in good agreement with existing JWST and ALMA data from ultraluminous quasars. We find the quasar is often obscured as it grows, and that strong, ejective feedback is required to have a high probability of detecting the quasar in the rest-frame UV. Fast and energetic quasar-driven winds expel metal-enriched gas, leading to significant metal pollution of the circumgalactic medium (CGM) out to twice the virial radius. As central gas densities and pressures are reduced, we find weaker signals from the CGM in mock X-ray and Sunyaev-Zeldovich maps, whose detection - with proposed instruments such as Lynx, and even potentially presently with ALMA - can constrain quasar feedback.

The detections of four apparently young radio pulsars in the Milky Way globular clusters are difficult to reconcile with standard neutron star formation scenarios associated with massive star evolution. Here we discuss formation of these young pulsars through white dwarf mergers in dynamically-old clusters that have undergone core collapse. Based on observed properties of magnetic white dwarfs, we argue neutron stars formed via white dwarf merger are born with spin periods of roughly $10-100\,$ms and magnetic fields of roughly $10^{11}-10^{13}\,$G. As these neutron stars spin down via magnetic dipole radiation, they naturally reproduce the four observed young pulsars in the Milky Way clusters. Rates inferred from $N$-body cluster simulations as well as the binarity, host cluster properties, and cluster offsets observed for these young pulsars hint further at a white dwarf merger origin. These young pulsars may be descendants of neutron stars capable of powering fast radio bursts analogous to the bursts observed recently in a globular cluster in M81.

The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy and dark matter. However, the development of high-fidelity theoretical models that cover a wide range of scales and redshifts is a necessary stepping-stone. Here, we present public high-resolution weak lensing maps on the light cone, generated using the $N$-body simulation suite AbacusSummit in the Born approximation, and accompanying weak lensing mock catalogues, tuned via fits to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument (DESI). Available in this release are maps of the cosmic shear, deflection angle and convergence fields at source redshifts ranging from $z = 0.15$ to 2.45 with $\Delta z = 0.05$ as well as CMB convergence maps ($z \approx 1090$) for each of the 25 ${\tt base}$-resolution simulations ($L_{\rm box} = 2000\,h^{-1}{\rm Mpc}$, $N_{\rm part} = 6912^3$) as well as for the two ${\tt huge}$ simulations ($L_{\rm box} = 7500\,h^{-1}{\rm Mpc}$, $N_{\rm part} = 8640^3$) at the fiducial AbacusSummit cosmology ($Planck$ 2018). The pixel resolution of each map is 0.21 arcmin, corresponding to a HEALPiX $N_{\rm side}$ of 16384. The sky coverage of the ${\tt base}$ simulations is an octant until $z \approx 0.8$ (decreasing to about 1800 deg$^2$ at $z \approx 2.4$), whereas the ${\tt huge}$ simulations offer full-sky coverage until $z \approx 2.2$. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper-Suprime Cam weak lensing datasets. The produced mock catalogues are validated against theoretical predictions for various clustering and lensing statistics such as galaxy clustering multipoles, galaxy-shear and shear-shear, showing excellent agreement.

Dust absorption is invoked in a number of contexts for hiding a star that has survived some sort of transient event from view. Dust formed in a transient is expanding away from the star and, in spherical models, the mass and energy budgets implied by a high optical depth at late times make such models untenable. Concentrating the dust in a disk or torus can in principle hide a source from an equatorial observer using less mass and so delay this problem. However, using axisymmetric dust radiation transfer models with a range of equatorial dust concentrations, we find that this is quite difficult to achieve in practice. The polar optical depth must be either low or high to avoid scattering optical photons to equatorial observers. Most of the emission remains at wavelengths easily observed by JWST, and the equatorial brightness is reduced by at most a factor of ~2 compared to isotropic emission even for equatorial (visual) optical depths of 1000. It is particularly difficult to hide a source with silicate dusts because the absorption feature near 10\ microns frequently leads to the emission being concentrated just bluewards of the feature, near 8 microns.

For 100 years since galaxies were found to be flying apart from each other, astronomers have been trying to determine how fast. The expansion, characterized by the Hubble constant, H0, is confused locally by peculiar velocities caused by gravitational interactions, so observers must obtain accurate distances at significant redshifts. Very nearby in our Galaxy, accurate distances can be determined through stellar parallaxes. There is no good method for obtaining galaxy distances that is applicable from the near domain of stellar parallaxes to the far domain free from velocity anomalies. The recourse is the distance ladder involving multiple methods with overlapping domains. Good progress is being made on this project, with satisfactory procedures and linkages identified and tested across the necessary distance range. Best values of H0 from the distance ladder lie in the range 73 - 75 km/s/Mpc. On the other hand, from detailed information available from the power spectrum of fluctuations in the cosmic microwave background, coupled with constraints favoring the existence of dark energy from distant supernova measurements, there is the precise prediction that H0 = 67.4 to 1%. If it is conclusively determined that the Hubble constant is well above 70 km/s/Mpc as indicated by distance ladder results then the current preferred LambdaCDM cosmological model based on the Standard Model of particle physics may be incomplete. There is reason for optimism that the value of the Hubble constant from distance ladder observations will be rigorously defined with 1% accuracy in the near future.

One important source of systematics in galaxy redshift surveys comes from the estimation of the galaxy window function. Up until now, the impact of the uncertainty in estimating the galaxy window function on parameter inference has not been properly studied. In this paper, we show that the uncertainty and the bias in estimating the galaxy window function will be salient for ongoing and next-generation galaxy surveys using a simulation-based approach. With a specific case study of cross-correlating Emission-line galaxies from the DESI Legacy Imaging Surveys and the Planck CMB lensing map, we show that neural network-based regression approaches to modelling the window function are superior in comparison to linear regression-based models. We additionally show that the definition of the galaxy overdensity estimator can impact the overall signal-to-noise of observed power spectra. Finally, we show that the additive biases coming from the window functions can significantly bias the modes of the inferred parameters and also degrade their precision. Thus, a careful understanding of the window functions will be essential to conduct cosmological experiments.

Ultra-fast outflows (UFOs) have been revealed in a large number of active galactic nuclei (AGN) and are regarded as promising candidates for AGN feedback on the host galaxy. The nature and launching mechanism of UFOs are not yet fully understood. Here we perform a time- and flux-resolved X-ray spectroscopy on four XMM-Newton observations of a highly accreting narrow-line Seyfert 1 (NLS1) galaxy, Mrk 1044, to study the dependence of the outflow properties on the source luminosity. We find that the UFO in Mrk 1044 responds to the source variability quickly and its velocity increases with the X-ray flux, suggesting a high-density ($10^{9}-4.5\times10^{12}\,\mathrm{cm}^{-3}$) and radiatively driven outflow, launched from the region within a distance of $98-6600\, R_\mathrm{g}$ from the black hole. The kinetic energy of the UFO is conservatively estimated ($L_\mathrm{UFO}\sim4.4\%L_\mathrm{Edd}$), reaching the theoretical criterion to affect the evolution of the host galaxy. We also find emission lines, from a large-scale region, have a blueshift of $2700-4500$ km/s in the spectra of Mrk 1044, which is rarely observed in AGN. By comparing with other sources, we propose a correlation between the blueshift of emission lines and the source accretion rate, which can be verified by a future sample study.

We study tidal dissipation in hot Jupiter host stars due to the nonlinear damping of tidally driven $g$-modes, extending the calculations of Essick & Weinberg (2016) to a wide variety of non-solar type hosts. This process causes the planet's orbit to decay and has potentially important consequences for the evolution and fate of hot Jupiters. Previous studies either only accounted for linear dissipation processes or assumed that the resonantly excited primary mode becomes strongly nonlinear and breaks as it approaches the stellar center. However, the great majority of hot Jupiter systems are in the weakly nonlinear regime in which the primary mode does not break but instead excites a sea of secondary modes via three-mode interactions. We simulate these nonlinear interactions and calculate the net mode dissipation for stars that range in mass from $0.5 M_\odot \le M_\star \le 2.0 M_\odot$ and in age from the early main sequence to the subgiant phase. For stars with $M_\star \lesssim 1.0 M_\odot$ of nearly any age, we find that the orbital decay time is $\lesssim 100 \textrm{ Myr}$ for orbital periods $P_{\rm orb} \lesssim 1 \textrm{ day}$. For $M_\star \gtrsim 1.2 M_\odot$, the orbital decay time only becomes short on the subgiant branch, where it can be $\lesssim 10 \textrm{ Myr}$ for $P_{\rm orb} \lesssim 2 \textrm{ days}$ and result in significant transit time shifts. We discuss these results in the context of known hot Jupiter systems and examine the prospects for detecting their orbital decay with transit timing measurements.

Galactic science encompasses a wide range of subjects in the study of the Milky Way and Magellanic Clouds, from Young Stellar Objects to X-ray Binaries. Mapping these populations, and exploring transient phenomena within them, are among the primary science goals of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). While early versions of the survey strategy dedicated relatively few visits to the Galactic Plane region, more recent strategies under consideration envision higher cadence within selected regions of high scientific interest. The range of galactic science presents a challenge in evaluating which strategies deliver the highest scientific returns. Here we present metrics designed to evaluate Rubin survey strategy simulations based on the cadence of observations they deliver within regions of interest to different topics in galactic science, using variability categories defined by timescale. We also compare the fractions of exposures obtained in each filter with those recommended for the different science goals. We find that the baseline_v2.x simulations deliver observations of the high-priority regions at sufficiently high cadence to reliably detect variability on timescales >10 d or more. Follow-up observations may be necessary to properly characterize variability, especially transients, on shorter timescales. Combining the regions of interest for all the science cases considered, we identify those areas of the Galactic Plane and Magellanic Clouds of highest priority. We recommend that these refined survey footprints be used in future simulations to explore rolling cadence scenarios, and to optimize the sequence of observations in different bandpasses.

Stripped-envelope supernovae (SESNe) are a subclass of core-collapse supernovae that are deficient in hydrogen (SN~IIb, SN~Ib) and possibly helium (SN~Ic) in their spectra. Their progenitors are likely stripped of this material through a combination of stellar winds and interactions with a close binary companion, but the exact ejecta mass ranges covered by each subtype and how it relates to the zero-age main-sequence progenitor mass is still unclear. Using a combination of semi-analytic modeling and numerical simulations, we discuss how the properties of SESN progenitors can be constrained through different phases of the bolometric light curve. We find that the light curve rise time is strongly impacted by the strength of radioactive nickel mixing and treatment of helium recombination. These can vary between events and are often not accounted for in simpler modeling approaches, leading to large uncertainties in ejecta masses inferred from the rise. Motivated by this, we focus on the late time slope, which is determined by gamma-ray leakage. We calibrate the relationship between ejecta mass, explosion energy, and gamma-ray escape time $T_0$ using a suite of numerical models. Application of the fitting function we provide to bolometric light curves of SESNe should result in ejecta masses with approximately 20\% uncertainty. With large samples of SESNe coming from current and upcoming surveys, our methods can be utilized to better understand the diversity and origin of the progenitor stars.

The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet as well as thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL's VPIC (vector particle-in-cell) code, we have run simulations of such mildly-relativistic, collisionless, weakly-magnetized plasmas and compute the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton to electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers.

Evidence of the negative impact of light pollution on ecosystems is increasing every year. Its monitoring and study requires the identification, characterisation and control of the emitting sources. This is the case of urban centres with outdoor lighting that spills light outside the place it is intended to illuminate. The quantity and nature of the pollutant (artificial light at night) depends on the lamps used and how they are positioned. This is important because a greater proportion of blue light means a greater scattering effect. In this study, we analysed the emissions of 100 urban centres in the north of Granada province (Spain), using International Space Station (ISS) images from 2012 and 2021, in order to compare the results with public lighting inventories and verify the validity of these data for characterising night-time lighting emissions. Using inference and cluster analysis techniques, we confirmed an overall increase in emissions and a shift in their colour towards blue, consistent with the results of the lighting inventory analysis. We concluded that it is possible to use ISS imagery to characterise artificial light emissions and the lighting that causes them, none the less there are a number of inherent problems with the data and the way it was collected that require the results to be interpreted with caution.

Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are of particular interest since they are thin (dynamically cold) and therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter models. However, current samples are limited to the Milky Way. With its large field of view, deep imaging sensitivity, and high angular resolution, the upcoming {\it Nancy Grace Roman Space Telescope} ({\it Roman}); presents a unique opportunity to significantly increase the number of observed streams and gaps. This paper presents a first exploration of the prospects for detecting gaps in streams in M31 and other nearby galaxies with resolved stars. We simulate the formation of gaps in a Palomar-5-like stream and generate mock observations of these gaps together with background stars in M31 and foreground Milky Way stellar fields. We assess {\it Roman}'s ability to detect gaps out to 10~Mpc through visual inspection and with the gap-finding tool {\it FindTheGap}. We conclude that gaps of $\approx 1.5$~kpc in streams that are created from subhalos of masses $\geq5 \times 10^6$ {\Msun} are detectable within a 2--3~Mpc volume in exposures of 1000s--1~hour. This volume contains $\approx$ 200 galaxies. Large samples of stream gaps in external galaxies will open up a new era of statistical analyses of gap characteristics in stellar streams and help constrain dark matter models.

The origin and evolution of gas in debris disks is still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only few C observations available. We create an overview of the C and CO content of debris disk gas and use it test state-of-the-art models. We use new and archival ALMA observations of CO and CI emission, complemented by CII data from Herschel, for a sample of 14 debris disks. This expands the number of disks with ALMA measurements of both CO and CI by ten disks. We present new detections of CI emission towards three disks: HD 21997, HD 121191 and HD 121617. We use a simple disk model to derive gas masses and column densities. We find that current state-of-the-art models of secondary gas production overpredict the neutral carbon content of debris disk gas. This does not rule out a secondary origin, but might indicate that the models require an additional C removal process. Alternatively, the gas might be produced in transient events rather than a steady-state collisional cascade. We also test a primordial gas origin by comparing our results to a simplified thermo-chemical model. This yields promising results, but more detailed work is required before a conclusion can be reached. Our work demonstrates that the combination of C and CO data is a powerful tool to advance our understanding of debris disk gas.

Using Asteroid Terrestrial-impact Last Alert System (ATLAS) and Zwicky Transient Facility (ZTF) data, I found that the SW Sex star V1315 Aql entered a low state early in 2023. As far as I know, this is the first such an event since the discovery of this object with observations dating back to 1948. This object is renowned for its nova shell and the nova explosion was estimated to occur 500-1200 yr ago, although no direct detection of a nova eruption was present in the historical record. The present low state probably occurred as the mass-transfer rate from the nova-irradiated secondary decreased secularly. The present low state would provide an opportunity to study the secondary in detail, which has been hampered by the luminous accretion disk in the past. I also provide the recent light curve of the classical nova LV Vul, which now shows high and low states 50 yr after the nova eruption.

By performing ideal magnetohydrodynamical (MHD) simulations with weak vertical magnetic fields in unstratified cylindrical shearing boxes with modified boundary treatment, we investigate MHD turbulence excited by magnetorotational instability. The cylindrical simulation exhibits extremely large temporal variation in the magnetic activity compared to the simulation in a normal Cartesian shearing box, although the time-averaged field strengths are comparable in the cylindrical and Cartesian setups. Detailed analysis of the terms describing magnetic-energy evolution with ``triangle diagrams'' surprisingly reveals that in the cylindrical simulation the compression of toroidal magnetic field is unexpectedly as important as the winding due to differential rotation in amplifying magnetic fields and triggering intermittent magnetic bursts, which are not seen in the Cartesian simulation. The importance of the compressible amplification is also true for a cylindrical simulation with tiny curvature; the evolution of magnetic fields in the nearly Cartesian shearing box simulation is fundamentally different from that in the exact Cartesian counterpart. The radial gradient of epicycle frequency, $\kappa$, which cannot be considered in the normal Cartesian shearing box model, is the cause of this fundamental difference. An additional consequence of the spatial variation of $\kappa$ is continuous and ubiquitous formation of narrow high(low)-density and weak(strong)-field localized structures; seeds of these ring-gap structures are created by the compressible effect and subsequently amplified and maintained under the marginally unstable condition regarding ``viscous-type'' instability.

Using solar wind observation near PSP perihelions as constraints, we have investigated the parameters in various PFSS model methods. It's found that the interplanetary magnetic field extrapolation with source surface height $R_\mathrm{SS} = 2\,Rs$ is better than that with $R_\mathrm{SS} = 2.5\,Rs$. HMI and GONG magnetograms show similar performance in the simulation of magnetic field variation, but the former appears to have a slight advantage in reconstruction of intensity while the latter is more adaptable to sparser grids. The finite-difference method of constructing eigenvalue problem for potential field can achieve similar accuracy as analytic method and greatly improve the computational efficiency. MHD modeling performs relatively less well in magnetic field prediction, but it is able to provide rich information about solar-terrestrial space.

Gamma-Ray Bursts (GRBs), being observed at high redshift (z = 9.4), vital to cosmological studies and investigating Population III stars. To tackle these studies, we need correlations among relevant GRB variables with the requirement of small uncertainties on their variables. Thus, we must have good coverage of GRB light curves (LCs). However, gaps in the LC hinder the precise determination of GRB properties and are often unavoidable. Therefore, extensive categorization of GRB LCs remains a hurdle. We address LC gaps using a 'stochastic reconstruction,' wherein we fit two pre-existing models (Willingale 2007; W07 and Broken Power Law; BPL) to the observed LC, then use the distribution of flux residuals from the original data to generate data to fill in the temporal gaps. We also demonstrate a model-independent LC reconstruction via Gaussian Processes. At 10% noise, the uncertainty of the end time of the plateau, its correspondent flux, and the temporal decay index after the plateau decreases, on average, by 33.3% 35.03%, and 43.32%, respectively for the W07, and by 33.3%, 30.78%, 43.9% for the BPL. The slope of the plateau decreases by 14.76% in the BPL. After using the Gaussian Process technique, we see similar trends of a decrease in uncertainty for all model parameters for both the W07 and BPL models. These improvements are essential for the application of GRBs as standard candles in cosmology, for the investigation of theoretical models and for inferring the redshift of GRBs with future machine learning analysis.

Motivated by the recent discovery of a low surface brightness diffuse emission, a supernova remnant (SNR) candidate, surrounding the young pulsar PSR~J0837--2454, we carry out a likelihood analysis of the $\gamma$-ray data obtained by the \emph{Fermi} Gamma-ray Space Telescope from August 2008 to November 2022. Using a 2D Gaussian spatial template, we detect a significant extended $\gamma$-ray emission with a 68\% containment radius of $\sim1^{\circ}.8$, which is spatially coincident with the new SNR candidate at $\sim12\sigma$ confidence level. The spectrum of the extended $\gamma$-ray emission, obtained in the energy range of 0.1-500.0 GeV, shows a significant spectral curvature at $\sim$1 GeV, with a log-parabola spectral shape. Several scenarios, such as the SNR, pulsar wind nebula, and pulsar halo, are discussed as the potential origins of the extended $\gamma$-ray emission, and our model fitting results are preferred for the SNR scenario.

The investigation of the data for quasi-periodic pulsations observed in the X-ray spectra of the accreting millisecond pulsar XTEJ 1807-294 allows some conclusions to be made about its main parameters - mass and angular momentum. Seven different geodesic models - namely RP, RP1, RP2, TP, TP1, WD and TD are applied in attempt to assess their ability to describe the properties of the central neutron star.

In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass function agree well with results from high-resolution N-body simulations over the full trained redshift range ($10 \le z \le 0$). Furthermore, we assess the temporal consistency of our SR model by constructing halo merger trees. We examine progenitors, descendants and mass growth along the tree branches. All statistical indicators demonstrate the ability of our SR model to generate satisfactory high-resolution simulations based on low-resolution inputs.

This manuscript analyzes lunar lander soil erosion models and trajectory models to calculate how much damage will occur to spacecraft orbiting in the vicinity of the Moon. The soil erosion models have considerable uncertainty due to gaps in our understanding of the basic physics. The results for ~40 t landers show that the Lunar Orbital Gateway will be impacted by 1000s to 10,000s of particles per square meter but the particle sizes are very small and the impact velocity is low so the damage will be slight. However, a spacecraft in Low Lunar Orbit that happens to pass through the ejecta sheet will sustain extensive damage with hundreds of millions of impacts per square meter: although they are small, they are in the hypervelocity regime, and exposed glass on the spacecraft will sustain spallation over 4% of its surface.

The recent discovery of a kilonova from the long duration gamma-ray burst, GRB 211211A, challenges classification schemes based on temporal information alone. Gamma-ray properties of GRB 211211A reveal an extreme event, which stands out among both short and long GRBs. We find very short variations (few ms) in the lightcurve of GRB 211211A and estimate ~1000 for the Lorentz factor of the outflow. We discuss the relevance of the short variations in identifying similar long GRBs resulting from compact mergers. Our findings indicate that in future gravitational wave follow-up campaigns, some long duration GRBs should be treated as possible strong gravitational wave counterparts.

We explore the use of Gibbs sampling in recovering the red noise power spectral density, the red noise Fourier coefficients, and the white noise parameters for the case of single pulsar analyses, and illustrate its effectiveness using the NANOGrav 11-year data set. We find that Gibbs sampling noise modeling (GM) is superior to the current standard Bayesian techniques (SM) for single pulsar analyses, yielding model parameter posteriors with significantly higher computational efficiency and fidelity. Furthermore, the output of GM contains posteriors for the Fourier coefficients that can be used to characterize the underlying red noise process of any pulsar's timing residuals which are absent in current implementations of SM. Through simulations, we demonstrate the potential for such coefficients to recover the overall shape of the spatial cross-correlations between pulsar pairs produced by gravitational waves.

We present time-frequency analysis, based on the Hilbert-Huang transform (HHT), of the evolution on the low-frequency quasi-periodic oscillations (LFQPOs) observed in the black hole X-ray binary MAXI J1820+070. Through the empirical mode decomposition (EMD) method, we decompose the light curve of the QPO component and measure its intrinsic phase lag between photons from different energy bands. We find that the QPO phase lag is negative (low energy photons lag behind high energy photons), meanwhile the absolute value of the lag increases with energy. By applying the Hilbert transform to the light curve of the QPO, we further extract the instantaneous frequency and amplitude of the QPO. Compared these results with those from the Fourier analysis, we find that the broadening of the QPO peak is mainly caused by the frequency modulation. Through further analysis, we find that these modulations could share a common physical origin with the broad-band noise, and can be well explained by the internal shock model of the jet.

We revisit the growth of curvature perturbations in non-minimal curvaton scenario with a non-trivial field metric $\lambda(\phi)$ where $\phi$ is an inflaton field, and incorporate the effect from the non-uniform onset of curvaton's oscillation in terms of an axion-like potential. The field metric $\lambda(\phi)$ plays a central role in the enhancement of curvaton field perturbation $\delta\chi$, serving as an effective friction term which can be either positive or negative, depending on the shape of $\lambda(\phi)$, namely the first derivative $\lambda_{,\phi}$. Our analysis reveals that $\delta\chi$ undergoes the superhorizon growth when the condition $\eta_\text{eff} \equiv - 2 \sqrt{2\epsilon} M_\text{Pl} { \lambda_{,\phi} \over \lambda} < -3$ is satisfied. This is analogous to the mechanism responsible for the amplification of curvature perturbations in the context of ultra-slow-roll inflation, namely the growing modes dominate curvature perturbations. As a case study, we examine the impact of a Gaussian dip in $\lambda(\phi)$ and conduct a thorough investigation of both the analytical and numerical aspects of the inflationary dynamics. Our findings indicate that the behavior of the curvaton perturbation during inflation is not solely determined by the depth of the dip in $\lambda(\phi)$. Rather, the shape of the dip also plays a significant role, a feature that has not been previously highlighted in the literature. Utilizing the $\delta \mathcal{N}$ formalism, we derive analytical expressions for both the final curvature power spectrum and the non-linear parameter $f_\text{NL}$ in terms of an axion-like curvaton's potential leading to the non-uniform curvaton's oscillation. Additionally, the resulting primordial black hole abundance and scalar-induced gravitational waves are calculated, which provide observational windows for PBHs.

Metal-poor stars are key to our understanding of the early stages of chemical evolution in the Universe. New multi-filter surveys, such as the Southern Photometric Local Universe Survey (S-PLUS), are greatly advancing our ability to select low-metallicity stars. In this work, we analyse the chemodynamical properties and ages of 522 metal-poor candidates selected from the S-PLUS data release 3. About 92% of these stars were confirmed to be metal-poor ([Fe/H] $\leq -1$) based on previous medium-resolution spectroscopy. We calculated the dynamical properties of a subsample containing 241 stars, using the astrometry from Gaia Data Release 3. Stellar ages are estimated by a Bayesian isochronal method formalized in this work. We analyse the metallicity distribution of these metal-poor candidates separated into different subgroups of total velocity, dynamical properties, and ages. Our results are used to propose further restrictions to optimize the selection of metal-poor candidates in S-PLUS. The proposed astrometric selection ($\mathrm{parallax}>0.85$ mas) is the one that returns the highest fraction of extremely metal-poor stars (16.3% have [Fe/H] $\leq -3$); the combined selection provides the highest fraction of very metal-poor stars (91.0% have [Fe/H] $\leq -2$), whereas the dynamical selection (eccentricity > 0.35 and diskness < 0.75) is better for targetting metal-poor (99.5% have [Fe/H] $\leq -1$). Using only S-PLUS photometric selections, it is possible to achieve selection fractions of 15.6%, 88.5% and 98.3% for metallicities below $-$3, $-$2 and $-$1, respectively. We also show that it is possible to use S-PLUS to target metal-poor stars in halo substructures such as Gaia-Sausage/Enceladus, Sequoia, Thamnos and the Helmi stream.

We report on the measurement of the thermal Sunyaev-Zel'dovich (tSZ) Effect in the circumgalactic medium (CGM) of 641,923 galaxies with $\rm M_\star$=$\rm 10^{9.8-11.3}M_\odot$ at $z<$0.5, pushing the exploration of tSZ Effect to lower-mass galaxies compared to previous studies. We cross-correlate the galaxy catalog of $WISE$ and $SuperCosmos$ with the Compton-$y$ maps derived from the combined data of $Atacama$ $Cosmology$ $Telescope$ and $Planck$. We improve on the data analysis methods (correcting for cosmic infrared background and Galactic dust, masking galaxy clusters and radio sources, stacking, aperture photometry), as well as modeling (taking into account beam smearing, "two-halo" term, zero-point offset). We have constrained the thermal pressure in the CGM of $\rm M_\star$=$\rm 10^{10.6-11.3}M_\odot$ galaxies for a generalized NFW profile and provided upper limits for $\rm M_\star$=$\rm 10^{9.8-10.6}M_\odot$ galaxies. The relation between $\rm M_{500}$ (obtained from an empirical $\rm M_\star$-$\rm M_{200}$ relation and a concentration factor) and $\rm \tilde Y^{sph}_{R500}$ (a measure of the thermal energy within R$_{500}$) is $>$2$\sigma$ steeper than the self-similarity and the deviation from the same that has been reported previously in higher mass halos. We calculate the baryon fraction of the galaxies, $f_b$, assuming the CGM to be at the virial temperature that is derived from $\rm M_{200}$. $f_b$ exhibits a non-monotonic trend with mass, with $\rm M_\star$=$\rm 10^{10.9-11.2}M_\odot$ galaxies being baryon sufficient.

We perform general relativistic radiation magnetohydrodynamics (MHD) simulations of super-Eddington accretion flows around a neutron star with a dipole magnetic field for modeling the galactic ultra-luminous X-ray source (ULX) exhibiting X-ray pulsations, Swift J0243.6+6124. Our simulations show the accretion columns near the magnetic poles, the accretion disk outside the magnetosphere, and the outflows from the disk. It is revealed that the effectively optically thick outflows, consistent with the observed thermal emission at $\sim10^7$ K, are generated if the mass accretion rate is much higher than the Eddington rate $\dot{M}_{\rm Edd}$ and the magnetospheric radius is smaller than the spherization radius. In order to explain the blackbody radius ($\sim 100-500$ km) without contradicting the reported spin period ($9.8~{\rm s}$) and spin-up rate ($\dot{P}=-2.22\times10^{-8}~{\rm s~s^{-1}}$), the mass accretion rate of $(200-1200)\dot{M}_{\rm Edd}$ is required. Since the thermal emission was detected in two observations with $\dot{P}$ of $-2.22\times10^{-8}~{\rm s~s^{-1}}$ and $-1.75\times10^{-8}~{\rm s~s^{-1}}$ but not in another with $\dot{P}=-6.8 \times10^{-9}~{\rm s~s^{-1}}$, the surface magnetic field strength of the neutron star in Swift J0243.6+6124 is estimated to be between $3\times10^{11}~{\rm G}$ and $4\times10^{12}~{\rm G}$. From this restricted range of magnetic field strength, the accretion rate would be $(200-500)\dot{M}_{\rm Edd}$ when the thermal emission appears and $(60-100)\dot{M}_{\rm Edd}$ when it is not detected. Our results support the hypothesis that the super-Eddington phase in the 2017-2018 giant outburst of Swift J0243.6+6124 is powered by highly super-Eddington accretion flows onto a magnetized neutron star.

All known Small Solar System Bodies have diameters between a few meters and a few thousands of kilometers. Based on the collisional evolution of Solar System Bodies, a larger number of asteroids with diameters down to $\sim 2$ m is thought to exist. As all Solar System Bodies, Small Bodies can be passive sources of high-energy gamma rays, produced by the interaction of energetic cosmic rays impinging on their surfaces. Since the majority of known asteroids are in orbits between Mars and Jupiter (in a region known as the Main Belt), we expect them to produce a diffuse emission close to the ecliptic plane. In this work we have studied the gamma-ray emission coming from the ecliptic using the data collected by the Large Area Telescope onboard the Fermi satellite. We have fit the results with simulations of the gamma-ray intensity at source level (calculated with the software FLUKA) to constrain the Small Solar System Bodies population. Finally, we have proposed a model describing the distribution of asteroid sizes and we have used the LAT data to constrain the gamma-ray emission expected from this model and, in turn, on the model itself.

Black holes with masses in excess of several billion solar masses have been found at redshifts 6-7.5, when the universe was less than 1 Gyr old. The existence of such supermassive black holes already in place at such early epochs has been challenging for theoretical models and distinguishing between different scenarios has prompted the search for their progenitors at earlier epochs. Here we present an extensive analysis of the JWST-NIRSpec spectrum (from the JADES survey) of GN-z11, an exceptionally luminous galaxy at z=10.6, revealing the detection of the high ionization [NeIV]$\lambda$2423 transition and semi-forbidden nebular lines tracing gas densities higher than $\rm 10^{10}~cm^{-3}$, typical of the Broad Line Region of Active Galactic Nuclei (AGN). These spectral features indicate that, in addition to star formation, GN-z11 also hosts an accreting black hole. We do not exclude a contribution from extreme stellar populations, however Wolf Rayet stars alone cannot account for many of the spectral properties. The spectrum also reveals a deep and blueshifted CIV$\lambda$1549 absorption trough, tracing an outflow with a velocity of $\sim 800-1000$ km/s, higher than typically observed in starburst galaxies, hence likely driven by the AGN. Assuming local virial scaling relations, we derive a black hole mass of $\rm \log{(M_{BH}/M_{\odot})}=6.2\pm 0.3$, accreting at about 5 times the Eddington rate. While super-Eddington accretion is probably episodic, if it has been occurring for the previous $\sim 100$ Myr, then the black hole could have potentially originated even from a stellar mass seed at z$\sim$12-15. We finally discuss that our finding naturally explains the high luminosity of GN-z11 and can also provide an explanation for its exceptionally high nitrogen abundance.

The recent discovery of new Active Galactic Nuclei (AGN) at z > 4 with JWST is revolutionising the black hole (BH) landscape at cosmic dawn, unveiling for the first time accreting BHs with masses of 10^6 - 10^7 Msun. To date, the most distant reside in CEERS-1019 at z=8.7 and GNz11 at z=10.6. Given the high rate of newly discovered high-z AGNs, more than 10 at z > 4, we wonder: are we really surprised to find them in the nuclei of z = 5 - 11 galaxies? Can we use the estimated properties to trace their origin? In this work, we predict the properties of 4 < z < 11 BHs and their host galaxies considering an Eddington-limited (EL) and a super-Eddington (SE) BH accretion scenario, using the Cosmic Archaeology Tool (CAT), a semi-analytical model for the formation and evolution of z > 4 AGNs and galaxies. We then calculate the transmitted spectral energy distribution of CAT synthetic candidates, representative of the estimated BH properties in CEERS-1019 and GNz11. We find that the estimated luminosity of high-z JWST detected AGNs are better reproduced by the SE model, where BHs descend from efficiently growing light and heavy seeds. Conversely, the host galaxy stellar masses are better matched in the EL model, in which all the systems detectable with JWST surveys JADES and CEERS appear to be descendants of heavy BH seeds. Our study suggests an evolutionary connection between systems similar to GNz11 at z=10.6 and CEERS-1019 at z=8.7 and supports the interpretation that the central point source of GNz11 could be powered by a super-Eddington (lambda_Edd = 2 - 3) accreting BH with mass 1.5 10^6 Msun, while CEERS-1019 harbours a more massive BH, with M_{BH} = 10^7 Msun, accreting sub-Eddington (lambda_Edd = 0.45 - 1), with a dominant emission from the host galaxy.

In this letter, we provide the first observational evidence of substantial collisionless damping (CD) modulation in the magnetohydrodynamic (MHD) turbulence cascade in Earth's magnetosheath using four Cluster spacecraft. Plasma turbulence is primarily shaped by the forcing on large scales and damping on small scales. Based on an improved compressible MHD decomposition algorithm, our observations demonstrate that CD enhances the anisotropy of compressible MHD modes due to their strong pitch angle dependence. The wavenumber distributions of slow modes are more stretched perpendicular to the background magnetic field ($\mathbf{B_0}$) under CD modulation compared to Alfv\'en modes. In contrast, fast modes are subject to a more significant CD modulation. Fast modes exhibit a scale-independent, slight anisotropy above the CD truncation scales, and their anisotropy increases as the wavenumbers fall below the CD truncation scales. As a result, CD affects the relative energy fractions in total compressible modes. Our findings take a significant step forward in comprehending the functions of CD in truncating the compressible MHD turbulence cascade and the consequential energy anisotropy in the wavevector space.

We present radio very long baseline interferometry (VLBI) and X-ray studies of the starburst galaxy NGC 3628. The VLBI observation at 1.5 GHz reveals seven compact (0.7$-$7 parsec) radio sources in the central $\sim$250 parsec region of NGC 3628. Based on their morphology, high radio brightness temperatures ($10^5-10^7$ K), and steep radio spectra, none of these seven sources can be associated with active galactic nuclei (AGNs); instead, they can be identified as supernova remnants (SNRs), with three of them appearing consistent with partial shells. Notably, one of them (W2) is likely a nascent radio supernova and appears to be consistent with the star formation rate of NGC 3628 when assuming a canonical initial mass function. The VLBI observation provides the first precise measurement of the diameter of the radio sources in NGC 3628, which allow us to fit a well-constrained radio surface brightness - diameter ($\Sigma-D$) correlation by including the detected SNRs. Furthermore, future VLBI observations can be conducted to measure the expansion velocity of the detected SNRs. In addition to our radio VLBI study, we analyze Chandra and XMM-Newton spectra of NGC 3628. The spectral fitting indicates that the SNR activities could well account for the observed X-ray emissions. Along with the Chandra X-ray image, it further reveals that the X-ray emission is likely maintained by the galactic-scale outflow triggered by SN activities. These results provide strong evidence that SN-triggered activities play a critical role in generating both radio and X-ray emissions in NGC 3628 and further suggest that the galaxy NGC 3628 is in an early stage of starbursts.

We study radio and X-ray emissions from IMBHs and explore the unified model for accretion and ejection processes. The radio band survey of IMBH (candidate) hosted galaxies indicates that only a small fraction ($\sim$0.6\%) of them are radio-band active. In addition, very long baseline interferometry observations reveal parsec-scale radio emission of IMBHs, further resulting in a lower fraction of actively ejecting objects (radio emission is produced by IMBHs other than hosts), which is consistent with a long quiescent state in the evolution cycle of IMBHs. Most (75\%, i.e., 3 out of 4 samples according to a recent mini-survey) of the radio-emitting IMBHs are associated with radio relics and there is also evidence of dual radio blobs from episodic ejecting phases. Taking the radio emission and the corresponding core X-ray emission of IMBH, we confirm a universal fundamental plane relation (FMP) of black hole activity. Furthermore, state transitions can be inferred by comparing a few cases in XRBs and IMBHs in FMP, i.e., both radio luminosity and emission regions evolve along these state transitions. These signatures and evidence suggest an analogy among all kinds of accretion systems which span from stellar mass to supermassive black holes, hinting at unified accretion and ejection physics. To validate the unified model, we explore the correlation between the scale of outflows (corresponding to ejection powers) and the masses of central engines; it shows that the largest scale of outflows $\hat{LS}_\mathrm{out}$ follows a power-law correlation with the masses of accretors $M_\mathrm{core}$, i.e., $\log{\hat{LS}_\mathrm{out}} = (0.73\pm0.01)\log{M_\mathrm{core}} - (3.34\pm0.10)$. In conclusion, this work provides evidence to support the claim that the ejection (and accretion) process behaves as scale-invariant and their power is regulated by the masses of accretors.

Main-sequence bolometric corrections (BC) and a standard BC-Teff relation are produced for TESS wavelengths using published physical parameters and light ratios from SED models of 209 detached double-lined eclipsing binaries. This and previous five-band (Johnson B, V, Gaia G, GBP, GRP) standard BC-Teff relations are tested by recovering luminosity (L) of the most accurate 341 single host stars (281 MS, 40 subgiants, 19 giants and one PMS). The recovered L of photometry is compared to L from published R and Teff. A very high correlation ($R^2$ = 0.9983) is achieved for this mixed sample. Error histograms of recovered and calculated L show peaks at 2 and 4 per cent respectively. The recovered L and the published Teff} were then used in $L = 4 \pi R^2 \sigma Teff^4$ to predict the standard R of the host stars. Comparison between the predicted and published R of all luminosity classes are found successful with a negligible offset associated with the giants and subgiants. The peak of the predicted R errors is found at 2 per cent, which is equivalent to the peak of the published R errors. Thus, a main-sequence BC-Teff relation could be used in predicting both L and R of a single star at any luminosity class, but this does not mean BC-Teff relations of all luminosity classes are the same because luminosity information could be more constrained by star's apparent magnitude $\xi$ than its BC since $m_{Bol} = \xi + BC_\xi$.

In cosmological simulations without baryons, the relation between the specific angular momentum $j_{\rm h}$ and mass $M_{\rm h}$ of galactic dark matter halos has the well-established form $j_{\rm h} \propto M_{\rm h}^{2/3}$. This is invariably adopted as the starting point in efforts to understand the analogous relation between the specific angular momentum $j_{\ast}$ and mass $M_{\ast}$ of the stellar parts of galaxies, which are often re-expressed relative to the corresponding halo properties through the retention fractions $f_j = j_{\ast} / j_{\rm h}$ and $f_M = M_{\ast} / M_{\rm h}$. An important caveat here is that the adopted $j_{\rm h} \propto M_{\rm h}^{2/3}$ relation could, in principle, be modified by the gravitational back-reaction of baryons on dark matter (DM). We have tested for this possibility by comparing the $j_{\rm h}$-$M_{\rm h}$ relations in the IllustrisTNG100 and TNG50 simulations that include baryons (full-physics runs) with their counterparts that do not (DM-only runs). In all cases, we find scaling relations of the form $j_{\rm h} \propto M_{\rm h}^{\alpha}$, with $\alpha \approx 2/3$ over the ranges of mass and redshift studied here: $M_{\rm h} \geq 10^{10} \, M_{\odot}$ and $0 \leq z \leq 2$. The values of $\alpha$ are virtually identical in the full-physics and DM-only runs at the same redshift. The only detectable effect of baryons on the $j_{\rm h}$-$M_{\rm h}$ relation is a slightly higher normalization, by 12%-15% at $z=0$ and by 5% at $z=2$. This implies that existing estimates of $f_j$ based on DM-only simulations should be adjusted downward by similar amounts. Finally, we discuss briefly some implications of this work for studies of galaxy formation.

Hub-filament systems are suggested to be the birth cradles of high-mass stars and clusters. We apply the FILFINDER algorithm to the integrated intensity maps of the 13CO (3-2) line to identify filaments in the G333 complex, and extract the velocity and intensity along the filament skeleton from moment maps. Clear velocity and density fluctuations are seen along the filaments, allowing us to fit velocity gradients around the intensity peaks. The velocity gradients fitted to the LAsMA data and ALMA data agree with each other over the scales covered by ALMA observations in the ATOMS survey. Changes of velocity gradient with scale indicate a ''funnel'' structure of the velocity field in PPV space, indicative of a smooth, continuously increasing velocity gradient from large to small scales, and thus consistent with gravitational acceleration. The typical velocity gradient corresponding to a 1 pc scale is ~1.6km/s/pc. Assuming free-fall, we estimate a kinematic mass within 1 pc of ~1190 M$_\odot$, which is consistent with typical masses of clumps in the ATLASGAL survey. We find direct evidence for gravitational acceleration from comparison of the observed accelerations to those predicted by free-fall onto dense hubs. On large scales, we find that the inflow may be driven by the larger scale structure, consistent with hierarchical structure in the molecular cloud and gas inflow from large to small scales. The hub-filament structures at different scales may be organized into a hierarchical system extending up to the largest scales probed, through the coupling of gravitational centers at different scales. We argue that the ''funnel'' structure in PPV space can be an effective probe for the gravitational collapse motions in molecular clouds. The large scale gas inflow is driven by gravity, implying that the molecular clouds in G333 complex may be in the state of global gravitational collapse.

Midplane heating induced by disk accretion plays a key role in determining the disk temperature particularly at the inner disk midplane where planets form. However, the efficiency of accretion heating has been not well constrained by observations. We construct two-dimensional models of the Class II disk around CW Tau, taking into account the midplane heating. The models are compared with the ALMA dust continuum observations at Bands 4, 6, 7 and 8, with an angular resolution of 0.1 arcsec. The observed brightness temperatures are almost wavelength-indenpendent at $\lesssim$10 au. We find that if the maximum dust size $a_{\rm max}$ is $\lesssim100~{\rm \mu m}$, the brightness temperatures predicted by the model exceed the observed values, regardless of the efficiency of accretion heating. The low observed brightness temperatures can be explained if millimeter scattering reduces the intensity. If the disk is passive, $a_{\rm max}$ needs to be either $\sim150~{\rm \mu m}$ or $\gtrsim$ few ${\rm cm}$. The accretion heating significantly increases the brightness temperature particularly when $a_{\rm max}\lesssim300~{\rm \mu m}$, and hence $a_{\rm max}$ needs to be either $\sim300~{\rm \mu m}$ or $\gtrsim$ few ${\rm cm}$. The midplane temperature is expected to be $\sim$1.5-3 times higher than the observed brightness temperatures, depending on the models. The dust settling effectively increases the temperature of the dust responsible for the millimeter emission in the active disk, which makes the model with $300~{\rm \mu m}$-sized dust overpredicts the brightness temperatures when strong turbulence is absent. Porous dust (porosity of 0.9) makes the accretion heating more efficient so that some sort of reduction in accretion heating is required. Future longer wavelength and higher angular resolution observations will help us constrain the heating mechanisms of the inner protoplanetary disks.

The photoionization model of narrow-line regions (NLRs) in active galactic nuclei (AGNs) has been investigated for decades. Many published models are restricted to simple linear scaling abundance relations, dust-free assumption, uniform AGN radiation field, and using one specific photoionization code, which restricts them from providing a satisfactory prediction on a broad range of AGN observations. Through a comprehensive investigation, here we present how the choice of abundance scaling relations, dust inclusion, AGN radiation fields, and different photoionization codes CLOUDY and MAPPINGS affect the predictions on the strength of strong UV, optical, and infrared emission lines. We find the dust-depleted radiation pressure-dominated AGN model built with the latest non-linear abundance sets and photoionization code MAPPINGS V are consistent with AGN observations across a broad range of wavelengths. We also assess new potential HII-AGN separation diagrams in the optical and UV wavelengths.

Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core-collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper-motion data from Gaia EDR3 and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions ($\beta\approx -0.4$$\pm$$0.1$) in the outskirts. We also robustly measure a dark central mass of roughly $800\pm300 \,$M$_{\odot}$, but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), or a dark population of stellar remnants of extent $\approx 0.016\,\rm pc \simeq 3300\,AU$. However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended ($\sim 0.034\, \rm{pc} \approx 7000\,\rm AU$). We use Monte Carlo $N$-body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e., IMBH vs. remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core.

We present the optical photometric and spectroscopic analysis of two type Iax SNe 2018cni and 2020kyg. SN 2018cni is a bright type Iax SN (M$_{V,peak}$ = $-$17.81$\pm$0.21 mag) whereas SN 2020kyg (M$_{V,peak}$ = $-$14.52$\pm$0.21 mag) is a faint one. We derive $^{56}$Ni mass of 0.07 and 0.002 M${_\odot}$, ejecta mass of 0.48 and 0.14 M${_\odot}$ for SNe 2018cni and 2020kyg, respectively. A combined study of the bright and faint type Iax SNe in $R/r$- band reveals that the brighter objects tend to have a longer rise time. However, the correlation between the peak luminosity and decline rate shows that bright and faint type Iax SNe exhibit distinct behaviour. Comparison with standard deflagration models suggests that SN 2018cni is consistent with the deflagration of a CO white dwarf whereas the properties of SN 2020kyg can be better explained by the deflagration of a hybrid CONe white dwarf. The spectral features of both the SNe point to the presence of similar chemical species but with different mass fractions. Our spectral modelling indicates stratification at the outer layers and mixed inner ejecta for both the SNe.

Cosmic rays (CRs) are charged particles that arrive at Earth isotropically from all directions and interact with the atmosphere. The presence of a spectral knee feature seen in the CR spectrum at $\sim$3 PeV energies is an evidence that astrophysical objects within our Galaxy, which are known as 'Galactic PeVatrons', are capable of accelerating particles to PeV energies. Scientists have been trying to identify the origin of Galactic CRs and have been looking for signatures of Galactic PeVatrons through neutral messengers. Recent advancements in ground-based $\gamma$-ray astronomy have led to the discovery of 12 Galactic sources emitting above 100 TeV energies, and even the first time detection of PeV photons from the direction of the Crab Nebula and the Cygnus region. These groundbreaking discoveries have opened up the field of ultra-high energy (UHE, E$>$100 TeV) $\gamma$-ray astronomy, which can help us explore the high energy frontiers of our Galaxy, hunt for PeVatron sources, and shed light on the century-old problem of the origin of CRs. This review article provides an overview of the current state of the art and potential future directions for the search for Galactic PeVatrons using ground-based $\gamma$-ray observations.

We determine the atomic hydrogen (HI) to halo mass relation (HIHM) using Arecibo Legacy Fast ALFA survey HI data at the location of optically selected groups from the Galaxy and Mass Assembly (GAMA) survey. We make direct HI detections for 37 GAMA groups. Using HI group spectral stacking of 345 groups, we study the group HI content as function of halo mass across a halo mass range of $10^{11} - 10^{14.7}\text{ M}_\odot$. We also correct our results for Eddington bias. We find that the group HI mass generally rises as a function of halo mass from $1.3\%$ of the halo mass at $10^{11.6} \text{M}_\odot$ to $0.4\%$ at $10^{13.7} \text{M}_\odot$ with some indication of flattening towards the high-mass end. Despite the differences in optical survey limits, group catalogues, and halo mass estimation methods, our results are consistent with previous group HI-stacking studies. Our results are also consistent with mock observations from SHARK and IllustrisTNG.

The optical spectrum of WR 138 exhibits emission lines typical of a WN6o star and absorption lines from a rapidly-rotating OB star. Using a large set of spectroscopic data, we establish a new orbital solution of the WN6o star based on the radial velocities of highly-ionized nitrogen lines. We show that the WN6o star moves on a 4.3 yr orbit with a comparatively low eccentricity of 0.16. The radial velocities of the OB star display considerable scatter. Our best estimates of the velocities of He I absorption lines result in a mass-ratio of $m_{\rm WN6o}/m_{\rm OB} = 0.53 \pm 0.09$. We disentangle the spectra of the two stars and derive a projected rotational velocity of $v\,\sin{i} = 350 \pm 10$ km s$^{-1}$ for the OB star. Our best orbital parameters, combined with the Gaia parallax of WR 138, are at odds with a previous interferometric detection of the companion, suggesting that there is either a bias in this detection or that WR 138 is actually a triple system.

The growth-rate $f\sigma_8(z)$ of the large-scale structure of the Universe is an important dynamic probe of gravity that can be used to test for deviations from General Relativity. However, in order for galaxy surveys to extract this key quantity from cosmological observations, two important assumptions have to be made: i) a fiducial cosmological model, typically taken to be the cosmological constant and cold dark matter ($\Lambda$CDM) model and ii) the modeling of the observed power spectrum from H$\alpha$ emitters, especially at non-linear scales, which is particularly dangerous as most models used in the literature are phenomenological at best. In this work, we propose a novel approach involving convolutional neural networks (CNNs), trained on the Quijote N-body simulations, to predict $f\sigma_8(z)$ directly and without assuming a model for the non-linear part of the power spectrum, thus avoiding the second of the aforementioned assumptions. We find that the predictions for the value of $f\sigma_8$ from the CNN are in excellent agreement with the fiducial values, while the errors are within a factor of order unity from those of the traditionally optimistic Fisher matrix approach, assuming an ideal fiducial survey matching the specifications of the Quijote simulations. Thus, we find the CNN reconstructions provide a viable alternative in order to avoid the theoretical modeling of the non-linearities at small scales when extracting the growth-rate.

We use the hemisphere comparison method to test the isotropy of the SnIa absolute magnitudes of the Pantheon+ and SH0ES samples in various redshift/distance bins. We compare the identified levels of anisotropy in each bin with Monte-Carlo simulations of corresponding isotropised data to estimate the frequency of such levels of anisotropy in the context of an underlying isotropic cosmological. We find that the identified levels of anisotropy in all bins are consistent with the Monte-Carlo isotropic simulated samples. However, in the real samples for both the Pantheon+ and the SH0ES cases we find sharp changes of the level of anisotropy occuring at distances less than $40Mpc$. For the Pantheon+ sample we find that the redshift bin $[0.005,0.01]$ is significantly more anisotropic than the other 5 redshift bins considered. For the SH0ES sample we find a sharp drop of the anisotropy level at distances larger than about $30Mpc$. These anisotropy transitions are relatively rare in the Monte-Carlo isotropic simulated data and occur in $2\%$ of the SH0ES simulated data and at about $7\%$ of the Pantheon+ isotropic simulated samples. This effect is consistent with the experience of an off center observer in a $30Mpc$ bubble of distinct physics or systematics.

Massive neutrinos are well-known to cause a characteristic suppression in the growth of structures at scales below the neutrino free-streaming length. A detailed understanding of this suppression is essential in the era of precision cosmology we are entering into, enabling us to better constrain the total neutrino mass and possibly probe (beyond)-$\Lambda$CDM cosmological model(s). Instead of the usual N-body simulation or Boltzmann solver, in this article we consider a two-fluid framework at the linear scales, where the neutrino fluid perturbations are coupled to the CDM (+ baryon) fluid via gravity at redshifts of interest. Treating the neutrino mass fraction $f_\nu$ as a perturbative parameter, we find a fully analytic solution to the system with redshift-dependent neutrino free-streaming length in $\Lambda$CDM background. The perturbative scale-dependent solution is shown to be in excellent agreement with numerical solution of the two-fluid equations valid to all orders in $f_{\nu}$, and also agrees with results from {\texttt{CLASS}} to a good accuracy. We further generalize the framework to incorporate different evolving dark energy backgrounds and found sub-percent level differences in the suppression, all of which lie within the observational uncertainty of BOSS-like surveys. We also present a brief discussion on the prospects of the current analysis in the context of upcoming missions.

The Cherenkov Telescope Array Observatory (CTAO) is a next-generation facility for ground-based very high energy gamma ray astronomy. CTAO will be operated as an open observatory. With two sites, in the northern and southern hemispheres, the Cherenkov Telescope Array CTA will provide full-sky coverage, improving sensitivity by an order of magnitude over current instruments, with a wide gamma ray energy coverage from 20 GeV to 300 TeV. CTA will use telescope arrays composed of three types of telescopes, optimized to cover different energy ranges. The large telescopes covering the lowest energies provide rapid slewing capability, for follow-up of transients. Key Science Projects (KSPs) are developed to form a significant part of the CTAO observing program during the first decade of operation, providing legacy data sets such as surveys or deep observations of key targets.

We search the archival Zwicky Transient Facility public survey for rapidly evolving transient (RET) candidates based on well-defined criteria between 2018 May and 2021 December. The search yielded 19 bona-fide RET candidates, corresponding to a discovery rate of $\sim 5.2$ events per year. Even with a Galactic latitude cut of $20^\circ$, 8 of the 19 events ($\sim 42$%) are Galactic, including one with a light-curve shape closely resembling that of the GW170817 kilonova (KN). An additional event is a nova in M31. Four out of the 19 events ($\sim 21$%) are confirmed extragalactic RETs (one confirmed here for the first time) and the origin of 6 additional events cannot be determined. We did not find any extragalactic events resembling the GW170817 KN, from which we obtain an upper limit on the volumetric rate of GW170817-like KNe of $R \le$ 2400 Gpc$^{-3}$ yr$^{-1}$ (95% confidence). These results can be used for quantifying contaminants to RET searches in transient alert streams, specifically when searching for kilonovae independently of gravitational-wave and gamma-ray-burst triggers.

(abridged) We present simulations of halo formation and evolution in scalar field dark matter (SFDM) cosmologies in the Thomas-Fermi regime, aka ``SFDM-TF", where a strong repulsive 2-particle self-interaction (SI) is included, being a valuable alternative to CDM, with the potential to resolve its ``cusp-core" problem. In general, SFDM behaves like a quantum fluid. Previous literature has presented two fluid approximations for SFDM-TF, as well as simulations of halo formation. These results confirmed earlier expectations and are generally in mutual agreement, but discrepancies were also reported. Therefore, we perform dedicated 3D cosmological simulations for the SFDM-TF model, applying both fluid approximations, as well as for CDM. Our results are very well in accordance with previous works and extend upon them, in that we can explain the reported discrepancies as a result of different simulation setups. We find some interesting details: The evolution of both SFDM-TF and CDM halos follows a 2-stage process. In the early stage, the density profile in the center becomes close to a $(n=1.5)$-polytropic core, dominated by an "effective" velocity-dispersion pressure $P_{\sigma}$ which is common to both dark matter models. Consecutively, for CDM halos, the core transitions into a central cusp. In SFDM-TF halos, the additional pressure $P_\text{SI}$ due to SI determines the second stage of the evolution, where the central region follows closely a $(n=1)$-polytropic core, embedded in a nearly isothermal envelope, i.e. the outskirts are similar to CDM. We also encounter a new effect, namely a late-time expansion of both polytropic core plus envelope, because the size of the almost isothermal halo envelope is affected by the expansion of the background universe. So, an initial primordial core of $\sim 100$ pc can evolve into a larger core of $\gtrsim 1$ kpc, even without feedback from baryons.

Model-independent mass limits assess the robustness of current cosmological measurements of the neutrino mass scale. Consistency between high-multipole and low-multiple Cosmic Microwave Background observations measuring such scale further valuate the constraining power of present data. We derive here up-to-date limits on neutrino masses and abundances exploiting either the Data Release 4 of the Atacama Cosmology Telescope (ACT) or the South Pole Telescope polarization measurements from SPT-3G, envisaging different non-minimal background cosmologies and marginalizing over them. Both the most constraining and marginalized bounds are competitive with those found with Planck data: we obtain $\sum m_\nu <0.139$ eV and $N_{\textrm{eff}}= 2.82\pm 0.25$ in a dark energy quintessence scenario, both at $95\%$ CL. These limits translate into $\sum m_\nu <0.20$ eV and $N_{\textrm{eff}}= 2.79^{+0.30}_{-0.28}$ after marginalizing over a plethora of well-motivated fiducial models. Our findings reassess both the strength and the reliability of cosmological neutrino mass constraints.

We study the clustering of primordial black holes (PBHs) and axion miniclusters produced in the model proposed to explain the LIGO/Virgo events or the seeds of the supermassive black holes (SMBHs) in arXiv:2006.13137. It is found that this model predicts large isocurvature perturbations due to the clustering of PBHs and axion miniclusters, from which we obtain stringent constraints on the model parameters. Specifically, for the axion decay constant $f_a=10^{16}~\mathrm{GeV}$, which potentially accounts for the seeds of the SMBHs, the PBH fraction in dark matter should be $f_\mathrm{PBH}\lesssim7\times 10^{-10}$. Assuming that the mass of PBHs increases by more than a factor of $\mathcal{O}(10)$ due to accretion, this is consistent with the observed abundance of SMBHs. On the other hand, for $f_a=10^{17}~\mathrm{GeV}$ required to produce PBHs of masses detected in the LIGO/Virgo, the PBH fraction should be $f_\mathrm{PBH}\lesssim6\times 10^{-8}$, which may be too small to explain the LIGO/Virgo events, although there is a significant uncertainty in calculating the merger rate in the presence of clustering.

We present an X-ray spectro-polarimetric analysis of the bright Seyfert galaxy IC 4329A. The Imaging X-ray Polarimetry Explorer (IXPE) observed the source for ~500 ks, supported by XMM-Newton (~60 ks) and NuSTAR (~80 ks) exposures. We detect polarisation in the 2-8 keV band with 2.97 sigma confidence. We report a polarisation degree of $3.3\pm1.1$ per cent and a polarisation angle of $78\pm10$ degrees (errors are 1 sigma confidence). The X-ray polarisation is consistent with being aligned with the radio jet, albeit partially due to large uncertainties on the radio position angle. We jointly fit the spectra from the three observatories to constrain the presence of a relativistic reflection component. From this, we obtain constraints on the inclination angle to the inner disc (< 39 degrees at 99 per cent confidence) and the disc inner radius (< 11 gravitational radii at 99 per cent confidence), although we note that modelling systematics in practice add to the quoted statistical error. Our spectro-polarimetric modelling indicates that the 2-8 keV polarisation is consistent with being dominated by emission directly observed from the X-ray corona, but the polarisation of the reflection component is completely unconstrained. Our constraints on viewer inclination and polarisation degree tentatively favour more asymmetric, possibly out-flowing, coronal geometries that produce more highly polarised emission, but the coronal geometry is unconstrained at the 3 sigma level.

We present High Sensitivity Array observation of the water megamasers of NGC 1068. We obtain absolute astrometry with 0.3 mas precision that confirms the association of the disk masers with the nuclear radio continuum source S1. The new observations reveal two new blueshifted groups of disk masers. We also detect the 22 GHz continuum on short interferometric baselines. The position-velocity diagram of the disk masers shows a curve consistent with a nonaxisymmetric distribution of maser spots. The curve is probably the result of spiral arms with a constant pitch angle of roughly 5 degrees. The disk kinematics are consistent with Keplerian rotation and low turbulent speeds. The inferred central mass is 17 million solar masses. On the basis of disk stability arguments, the mass of the molecular disk is roughly 110 thousand solar masses. The disk masers further resolve into filamentary structures suggesting an ordered magnetic field threading the maser disk. The magnetic field strengths must be greater than 1.6 mG to withstand turbulent motions in the partially ionized molecular gas. We note apparent asymmetries in the molecular disk that might be explained by anisotropic heating by a misaligned inner accretion disk. The new observations also detect the fainter jet masers north of the disk masers. The distribution and kinematics of the jet masers are consistent with an expanding ring of molecular gas.

We present a first study based on the analysis of the DEep Spectra of Ionized REgions Database (DESIRED). This is a compilation of 190 high signal-to-noise ratio optical spectra of HII regions and other photoionized nebulae, mostly observed with 8-10m telescopes and containing $\sim$29380 emission lines. We find that the electron density --$n_{\rm e}$-- of the objects is underestimated when [SII] $\lambda6731/\lambda6716$ and/or [OII] $\lambda3726/\lambda3729$ are the only density indicators available. This is produced by the non-linear density dependence of the indicators in the presence of density inhomogeneities. The average underestimate is $\sim 300$ cm$^{-3}$ in extragalactic HII regions, introducing systematic overestimates of $T_{\rm e}$([OII]) and $T_{\rm e}$([SII]) compared to $T_{\rm e}$([NII]). The high-sensitivity of [OII] $\lambda\lambda7319+20+30+31/\lambda\lambda3726+29$ and [SII] $\lambda\lambda4069+76/\lambda\lambda6716+31$ to density makes them more suitable for the diagnosis of the presence of high-density clumps. If $T_{\rm e}$([NII]) is adopted, the density underestimate has a small impact in the ionic abundances derived from optical spectra, being limited to up to $\sim$0.1 dex when auroral [SII] and/or [OII] lines are used. However, these density effects are critical for the analysis of infrared fine structure lines, such as those observed by the JWST in local star forming regions, implying strong underestimates of the ionic abundances. We present temperature relations between $T_{\rm e}$([OIII]), $T_{\rm e}$([ArIII]), $T_{\rm e}$([SIII]) and $T_{\rm e}$([NII]) for the extragalactic HII regions. We confirm a non-linear dependence between $T_{\rm e}$([OIII])-$T_{\rm e}$([NII]) due to a more rapid increase of $T_{\rm e}$([OIII]) at lower metallicities.

Aims: The study at hand aims to describe the distribution of atomic carbon, C0, throughout the envelope, in support of an improved understanding of its photo-chemistry. Additionally, we also briefly discuss the observation of [CII] emission towards the star. Methods: We obtain spectra of the [CI] $\mathrm{^3P_1} \rightarrow \mathrm{^3P_0}$ fine structure line at projected distances of up to 78" from the star. The line profiles are characterized by both direct fitting of Gaussian components, and by modeling the observed line of the [CI] triplet. We also report the detection of the $\mathrm{^2P_{3/2}} \rightarrow \mathrm{^2P_{1/2}}$ line from the C+ fine structure singlet at the central position and at 32" from the star. Results: The overall picture of the [CI] emission from IRC +10216 agrees with more limited previous studies. The satisfying agreement between the observed and modeled line profiles, with emission at the systemic velocity appearing beyond one beam from the star, rules out that the C0 is located in a thin shell. Given that the bond energy of CO falls only 0.1 eV below the ionization threshold of C0, the absence of observable [CII] emission from sightlines beyond a projected distance of $\sim 10^{17}$ cm from the star (adopting a distance of 130 pc) does not contradict a scenario where the bulk of C0 is located between that of CO and C+, as expected for an external FUV radiation field. This conjecture is also corroborated by a model in which the C0 shell is located farther outside, failing to reproduce the [CI] line profiles at intermediate sky-plane distances from the star. Comparing a photo-chemical model adopted from literature with the simplifying assumption of a constant C0 abundance with respect to the $\mathrm{H}_2$ density, we constrain the inner boundary of the [CI] emitting shell, located at $\sim 10^{16}$ cm from the star.

One of the major sources of perturbation in the solar cycle amplitude is believed to be the emergence of anomalous active regions which do not obey Hale's polarity law and Joy's law of tilt angles. Anomalous regions containing high magnetic flux that disproportionately impact the polar field are sometimes referred to as "rogue regions". In this study -- utilizing a surface flux transport model -- we analyze the large-scale dipole moment build-up due to the emergence of anomalous active regions on the solar surface. Although these active regions comprise a small fraction of the total sunspot number, they can substantially influence the magnetic dipole moment build-up and subsequent solar cycle amplitude. Our numerical simulations demonstrate that the impact of "Anti-Joy" regions on the solar cycle is similar to those of "Anti-Hale" regions. We also find that the emergence time, emergence latitude, relative number and flux distribution of anomalous regions influence the large-scale magnetic field dynamics in diverse ways. We establish that the results of our numerical study are consistent with the algebraic (analytic) approach to explaining the Sun's dipole moment evolution. Our results are relevant for understanding how anomalous active regions modulate the Sun's large-scale dipole moment build-up and its reversal timing within the framework of the Babcock-Leighton dynamo mechanism -- now believed to be the primary source of solar cycle variations.

We present MeerKAT Fornax Survey atomic hydrogen (HI) observations of the dwarf galaxies located in the central ~2.5 x 4 deg$^2$ of the Fornax galaxy cluster. The HI images presented in this work have a $3\sigma$ column density sensitivity between 2.7 and 50 x 10$^{18}$ cm$^{-2}$ over 25 km s$^{-1}$ for spatial resolution between 4 and 1 kpc. We are able to detect an impressive MHI = 5 x 10$^{5}$ Msun 3$\sigma$ point source with a line width of 50 km s$^{-1}$ at a distance of 20 Mpc. We detect HI in 17 out of the 304 dwarfs in our field -- 14 out of the 36 late type dwarfs (LTDs), and 3 of the 268 early type dwarfs (ETDs). The HI-detected LTDs have likely just joined the cluster and are on their first infall as they are located at large clustocentric radii, with comparable MHI and mean stellar surface brightness at fixed luminosity as blue, star-forming LTDs in the field. The HI-detected ETDs have likely been in the cluster longer than the LTDs and acquired their HI through a recent merger or accretion from nearby HI. Eight of the HI-detected LTDs host irregular or asymmetric HI emission and disturbed or lopsided stellar emission. There are two clear cases of ram-pressure shaping the HI, with the LTDs displaying compressed HI on the side closest to the cluster centre and a one-sided, starless tail pointing away from the cluster centre. The HI-detected dwarfs avoid the most massive potentials, consistent with massive galaxies playing an active role in the removal of HI. We create a simple toy model to quantify the timescale of HI stripping in the cluster. We find that a MHI = 10$^{8}$ Msun dwarf will be stripped in ~ 240 Myr. The model is consistent with our observations, where low mass LTDs are directly stripped of their HI from a single encounter and more massive LTDs can harbour a disturbed HI morphology due to longer times or multiple encounters being required to fully strip their HI.

We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) method to match the morphology in a 3D reconstruction of the dust density (3D dust) and the distribution of the gas tracers' emission. The result of the HOG analysis is a map of the relationship between the distances and radial velocities. The HOG comparison between the 3D dust and the HI emission indicates a morphological match at the distance of Taurus but an anti-correlation between the dust density and the HI emission, which uncovers a significant amount of cold HI within the Taurus MC. The HOG between the 3D dust and $^{12}$CO reveals a pattern in radial velocities and distances that is consistent with converging motions of the gas in the Taurus MC, with the near side of the cloud moving at higher velocities and the far side moving at lower velocities. This convergence of flows is likely triggered by the large-scale gas compression caused by the interaction of the Local Bubble and the Per-Tau shell, with Taurus lying at the intersection of the two bubble surfaces.

Previous STIS long-slit observations of Eta Carinae identified numerous absorption features in both the stellar spectrum, and in the adjacent nebular spectra, along our line-of-sight. The absorption features became temporarily stronger when the ionizing FUV radiation field was reduced by the periastron passage of the secondary star. Subsequently, dissipation of a dusty structure in our LOS has led to a long-term increase in the apparent magnitude of \ec, an increase in the ionizing UV radiation, and the disappearance of absorptions from multiple velocity-separated shells extending across the foreground Homunculus lobe. We use HST/STIS spectro-images, coupled with published infrared and radio observations, to locate this intervening dusty structure. Velocity and spatial information indicate the occulter is ~1000 au in front of Eta Carinae. The Homunculus is a transient structure composed of dusty, partially-ionized ejecta that eventually will disappear due to the relentless rain of ionizing radiation and wind from the current binary system along with dissipation and mixing with the ISM. This evolving complex continues to provide an astrophysical laboratory that changes on human timescales.

In this work we present, for the first time, infrared spectra of Titan from the Spitzer Space Telescope ($2004-2009$). The data are from both the short wavelength-low resolution (SL, $5.13-14.29\mathrm{\mu m}, R\sim60-127$) and short wavelength-high resolution channels (SH, $9.89 - 19.51\mathrm{\mu m}, R\sim600$) showing the emissions of CH$_{4}$, C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, C$_{2}$H$_{6}$, C$_{3}$H$_{4}$, C$_{3}$H$_{6}$, C$_{3}$H$_{8}$, C$_{4}$H$_{2}$, HCN, HC$_{3}$N, and CO$_{2}$. We compare the results obtained for Titan from Spitzer to those of the Cassini Composite Infrared Spectrometer (CIRS) for the same time period, focusing on the $16.35-19.35\mathrm{\mu m}$ wavelength range observed by the SH channel but impacted by higher noise levels in CIRS observations. We use the SH data to provide estimated haze extinction cross-sections for the $16.67-17.54\mathrm{\mu m}$ range that are missing in previous studies. We conclude by identifying spectral features in the $16.35-19.35\mathrm{\mu m}$ wavelength range, including two prominent emission features at 16.39 and $17.35\mathrm{\mu m}$, that could be analyzed further through upcoming James Webb Space Telescope Cycle 1 observations with the Mid-Infrared Instrument ($5.0-28.3\mathrm{\mu m}, R\sim1500-3500$). We also highlight gaps in current spectroscopic knowledge of molecular bands, including candidate trace species such as C$_{60}$ and detected trace species such as C$_{3}$H$_{6}$, that could be addressed by theoretical and laboratory study.

Giant pulses emitted by PSR B1937+21 are bright, intrinsically impulsive bursts. Thus, the observed signal from a giant pulse is a noisy but direct measurement of the impulse response from the ionized interstellar medium. We use this fact to detect 13,025 giant pulses directly in the baseband data of two observations of PSR B1937+21. Using the giant pulse signals, we model the time-varying impulse response with a sparse approximation method, in which the time dependence at each delay is decomposed in Fourier components, thus constructing a wavefield as a function of delay and differential Doppler shift. We find that the resulting wavefield has the expected parabolic shape, with several diffuse structures within it, suggesting the presence of multiple scattering locations along the line of sight. We also detect an echo at a delay of about 2.4 ms, over 1.5 times the rotation period of the pulsar, which between the two observations moves along the trajectory expected from geometry. The structures in the wavefield are insufficiently sparse to produce a complete model of the system, and hence the model is not predictive across gaps larger than about the scintillation time. Nevertheless, within its range, it reproduces about 75% of the power of the impulse response, a fraction limited mostly by the signal-to-noise ratio of the observations. Furthermore, we show that by deconvolution, using the model impulse response, we can successfully recover the intrinsic pulsar emission from the observed signal.

SNe Ia play a key role in the fields of astrophysics and cosmology. It is widely accepted that SNe Ia arise from thermonuclear explosions of WDs in binaries. However, there is no consensus on the fundamental aspects of the nature of SN Ia progenitors and their explosion mechanism. This fundamentally flaws our understanding of these important astrophysical objects. We outline the diversity of SNe Ia and the proposed progenitor models and explosion mechanisms. We discuss the recent theoretical and observational progress in addressing the SN Ia progenitor and explosion mechanism in terms of the observables at various stages of the explosion, including rates and delay times, pre-explosion companion stars, ejecta-companion interaction, early excess emission, early radio/X-ray emission from CSM interaction, surviving companions, late-time spectra and photometry, polarization signals, and SNR properties, etc. Despite the efforts from both the theoretical and observational side, the questions of how the WDs reach an explosive state and what progenitor systems are more likely to produce SNe Ia remain open. No single published model is able to consistently explain all observational features and the full diversity of SNe Ia. This may indicate that either a new progenitor paradigm or the improvement of current models is needed if all SNe Ia arise from the same origin. An alternative scenario is that different progenitor channels and explosion mechanisms contribute to SNe Ia. In the next decade, the ongoing campaigns with the JWST, Gaia and the ZTF, and upcoming extensive projects with the LSST and the SKA will allow us to conduct not only studies of individual SNe Ia in unprecedented detail but also systematic investigations for different subclasses of SNe Ia. This will advance theory and observations of SNe Ia sufficiently far to gain a deeper understanding of their origin and explosion mechanism.

The singlet scalar Higgs portal model provides one of the simplest explanations of dark matter in our Universe. Its Higgs resonant region, $m_\text{DM}\approx m_h/2$, has gained particular attention, being able to reconcile the tension between the relic density measurement and direct detection constraints. Interestingly, this region is also preferred as an explanation of the Fermi-LAT $\gamma$-ray Galactic center excess. We perform a detailed study of this model using $\gamma$-ray data from the Galactic center and from dwarf spheroidal galaxies and combine them with cosmic-ray antiproton data from the AMS-02 experiment that shows a compatible excess. In the calculation of the relic density, we take into account effects of early kinetic decoupling relevant for resonant annihilation. The model provides excellent fits to the astrophysical data either in the case the dark matter candidate constitutes all or a subdominant fraction of the observed relic density. We show projections for future direct detection and collider experiments to probe these scenarios.

We study the impact of renormalization group effects on QCD axion phenomenology. Focusing on the DFSZ model, we argue that the relevance of running effects for the axion couplings crucially depends on the scale where the heavier Higgs doublet, charged under the Peccei-Quinn symmetry, is integrated out. We study the impact of these effects on astrophysical and cosmological bounds as well as on the sensitivity of helioscopes experiments such as IAXO and XENONnT, showing that they can be sizable even in the most conservative case in which the two Higgs doublets remain as light as the TeV scale. We provide simple analytical expressions that accurately fit the numerical solutions of the renormalization group equations as a function of the mass scale of the heavy scalars.

We continue our work on a proposal for what dark matter could be, namely that the dark matter consists of essentially macroscopic objects built from ordinary matter. The only element of new physics is that there should exist several types or phases of vacuum. Then the dark matter particles are bubbles of a new type of vacuum filled with ordinary matter - say diamond - under high pressure and the whole bubble is supposedly contaminated by the materials in the extra galactic or intragalactic space able to form dust. The dust has a lower than dimension 3 structure and may well be long chains of atoms. We begin speculations on identifying the phase transition between the vacua, which we require in our model, with a (second order?) transition revealed in the Columbia plot with quark masses as axes. In the corner of small quark masses there is a true phase transition as temperature is raised, while outside this region there is instead no genuine phase transition under variation of the temperature but rather only a crossover.

An impressive feature of loop quantum gravity (LQG) is that it can elegantly resolve both the big bang and black hole singularities. By using the Newman-Janis algorithm, a regular and effective rotating self-dual black hole(SDBH) metric could be constructed, which alters the Kerr geometry with a polymeric function $P$ from the quantum effects of LQG geometry. In this paper, we investigate its impact on the frequency characteristics of the X-ray quasi-periodic oscillations(QPOs) from 5 X-ray binaries and contrast it with the existing results of the orbital, periastron precession and nodal precession frequencies within the relativistic precession model. We apply a Monte Carlo Markov Chain (MCMC) simulation to examine the possible LQG effects on the X-ray QPOs. We found that the best constraint result for the rotating self-dual geometry from LQG came from the QPOs of X-ray binary GRO J1655-40, which establish an upper bound on the polymeric function $P$ less than $8.6\times 10^{-4}$ at 95\% confidence level. This bound leads to a restriction on the polymeric parameter $\delta$ of LQG to be 0.24.

In this paper, the orbits of a charged particle near the event horizon of a magnetized black hole are investigated. For a static black hole of mass $M$ immersed in a homogeneous magnetic field $B$, the dimensionless parameter $b=eBGM/ (mc^4)$ controls the radius of the circular orbits and determines the position of the innermost stable circular orbit (ISCO), where $m$ and $e$ are the mass and charge of the particle. For large values of the parameter $b$, the ISCO radius can be very close to the gravitational radius. We demonstrate that the properties of such orbits can be effectively and easily found by using a properly constructed ``near-horizon approximation''. In particular, we show that the effective potential (which determines the position of the orbit) can be written in a form which is invariant under rescaling of the magnetic field, and as a result is universal in this sense. We also demonstrate that in the near-horizon approximation, the particle orbits are stationary worldlines in Minkowski spacetime. We use this property to solve the equation describing slow changes in the distance of the particle orbit from the horizon, which arise as a result of the electromagnetic field radiated by the particle itself. This allows us to evaluate the life-time of the particle before it reaches the ISCO and ultimately falls into the black hole.

For the light relativistic dark matter (DM) boosted by high energy cosmic ray, its scattering cross section with the nucleon is sensitively dependent on the momentum-transfer and such an dependence is caused by the mediator in the scattering. For puffy DM particle with a size, the momentum-transfer dependence can also arise from the DM radius effect. All these momentum-transfer dependences should be considered. In this note we study the direct detection limits on the cosmic ray-boosted puffy DM for a simplified model with a light mediator. For comparison, we first re-derive the direct detection limits on the cosmic ray-boosted point-like DM. We display the limits on various planes of parameters and find that the limits for the cosmic ray-boosted puffy DM are stronger than for the point-like DM.

The screening condition in neutron star core has been formulated as equality of velocities of superconducting protons and the electrons $\mathbf{v}_p=\mathbf{u}_e$ at wavenumbers $q\ll\lambda^{-1}$ ($\lambda$ is the London depth) and has been used to derive the force between the electronic flow past the flux tube, which has astrophysical applications. By calculating the current-current response, I find that $\mathbf{v}_p\neq\mathbf{u}_e$ for $l^{-1}<q\ll\lambda^{-1}$ ($l$ is the electron mean free path) at typical realistic parameters. Therefore, the momentum exchange between the electrons and the flux tubes in the core of neutron stars remains an open question.

In this paper, we investigate the Blandford-Znajek (BZ) process within the framework of Einsteinian cubic gravity (ECG). To analytically study the BZ process using the split monopole configuration, we construct a slowly rotating black hole in ECG up to cubic order in small spin, considering the leading order in small coupling constant of higher curvature terms. By deriving the magnetosphere solution around the black hole, we determine the BZ power up to the second relative order in spin. The BZ power is modified by the coupling constant compared to Kerr black hole. Although the general nature of the BZ process in ECG remains unchanged at the leading order in spin, the coupling constant introduces modification at the second relative order in spin. Therefore, we anticipate that it is feasible to discern general relativity from higher derivative gravities by examining the BZ power in rapidly rotating black holes.

Black hole images are theoretically predicted (under mild astrophysical assumptions) to display a stack of lensed "photon rings" that carry information about the underlying spacetime geometry. Despite vigorous efforts, no such ring has been observationally resolved thus far. However, planning is now actively under way for space missions targeting the first (and possibly the second) photon rings of the supermassive black holes M87* and Sgr A*. In this work, we study interferometric photon ring signatures in time-averaged images of Kerr black holes surrounded by different astrophysical profiles. We focus on the first, most easily accessible photon ring, which has a larger width-to-diameter ratio than subsequent rings and whose image consequently lacks a sharply defined diameter. Nonetheless, we show that it does admit a precise angle-dependent diameter in visibility space, for which the Kerr metric predicts a specific functional form that tracks the critical curve. We find that a measurement of this interferometric ring diameter is possible for most astrophysical profiles, paving the way for precision tests of strong-field general relativity via near-future observations of the first photon ring.

Spacecraft and drones aimed at exploring our solar system are designed to operate in conditions where the smart use of onboard resources is vital to the success or failure of the mission. Sensorimotor actions are thus often derived from high-level, quantifiable, optimality principles assigned to each task, utilizing consolidated tools in optimal control theory. The planned actions are derived on the ground and transferred onboard where controllers have the task of tracking the uploaded guidance profile. Here we argue that end-to-end neural guidance and control architectures (here called G&CNets) allow transferring onboard the burden of acting upon these optimality principles. In this way, the sensor information is transformed in real time into optimal plans thus increasing the mission autonomy and robustness. We discuss the main results obtained in training such neural architectures in simulation for interplanetary transfers, landings and close proximity operations, highlighting the successful learning of optimality principles by the neural model. We then suggest drone racing as an ideal gym environment to test these architectures on real robotic platforms, thus increasing confidence in their utilization on future space exploration missions. Drone racing shares with spacecraft missions both limited onboard computational capabilities and similar control structures induced from the optimality principle sought, but it also entails different levels of uncertainties and unmodelled effects. Furthermore, the success of G&CNets on extremely resource-restricted drones illustrates their potential to bring real-time optimal control within reach of a wider variety of robotic systems, both in space and on Earth.

Nobel Prize laureate P.J.E. Peebles [24] has emphasized the importance and difficulties of studying the large scale clustering of matter in cosmology. Nonlinear gravitational instability plays a central role in understanding the clustering of matter and the formation of nonlinear structures in the universe and stellar systems. However, there is no rigorous result on the nonlinear analysis of this instability except for some particular exact solutions without pressure, and numerical and phenomenological approaches. Both Rendall [26] and Mukhanov [21] have highlighted the challenge posed by nonlinear gravitational instability with effective pressure. This has been a longstanding open problem in astrophysics for over a century since the occurrence of linearized Jeans instabilities in Newtonian universes in 1902. This article contributes to a fully nonlinear analysis of the gravitational instability for the Euler-Possion system which models expanding Newtonian universes with inhomogeneous density and pressure. The exponential or finite-time increasing blowups of the density contrast $\varrho$ can be determined, which may account for the considerably faster growth rate of nonlinear structures observed in astrophysics than that suggested by the classical Jeans instability. We believe this is the first rigorous result for the nonlinear Jeans instability with effective pressure and the method is concise and robust.