Papers for Friday, Apr 14 2023

Dark matter subhalos with extended profiles and density cores, and globular stars clusters of mass $10^6-10^8 M_\odot$, that live near the critical curves in galaxy cluster lenses can potentially be detected through their lensing magnification of stars in background galaxies. In this work we study the effect such subhalos have on lensed images, and compare to the case of more well studied microlensing by stars and black holes near critical curves. We find that the cluster density gradient and the extended mass distribution of subhalos are important in determining image properties. Both lead to an asymmetry between the image properties on the positive and negative parity sides of the cluster that is more pronounced than in the case of microlensing. For example, on the negative parity side, subhalos with cores larger than about $50\,$pc do not generate any images with magnification above $\sim 100$ outside of the immediate vicinity of the cluster critical curve. We discuss these factors using analytical and numerical analysis, and exploit them to identify observable signatures of subhalos: subhalos create pixel-to-pixel flux variations of $\gtrsim 0.1$ magnitudes, on the positive parity side of clusters. These pixels tend to cluster around (otherwise invisible) subhalos. Unlike in the case of microlensing, signatures of subhalo lensing can be found up to $1''$ away from the critical curves of massive clusters.

We perform X-ray spectral analyses to derive characteristics (e.g., column density, X-ray luminosity) of $\approx$10,200 active galactic nuclei (AGNs) in the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS), which was designed to investigate the growth of supermassive black holes across a wide dynamic range of cosmic environments. Using physical torus models (e.g., Borus02) and a Bayesian approach, we uncover 22 representative Compton-thick (CT; $N_{\rm H} \;>\; 1.5\times10^{24}\; \rm cm^{-2}$) AGN candidates with good signal-to-noise ratios as well as a large sample of 136 heavily obscured AGNs. We also find an increasing CT fraction (\fct ) from low ($z<0.75$) to high ($z>0.75$) redshift. Our CT candidates tend to show hard X-ray spectral shapes and dust extinction in their SED fits, which may shed light on the connection between AGN obscuration and host-galaxy evolution.

The flavor composition of high-energy neutrinos carries important information about their birth. However, the two most common production scenarios, $pp$ (hadronuclear) and $p\gamma$ (photohadronic) processes, lead to the same flavor ratios when neutrinos and antineutrinos cannot be distinguished. The Glashow resonant interaction $\bar{\nu}_e+e^- \rightarrow W^-$ becomes a window to differentiate the antineutrino contribution from the total diffuse neutrino flux, thus lifting this degeneracy. We examine the power of Glashow resonant events in measuring the fraction of the $\bar{\nu}_e$ flux with current IceCube data, and produce projected sensitivities based on the combined exposure of planned Cherenkov neutrino telescopes around the globe. We find that $pp$ and $p\gamma$ can be distinguished at a 2$\sigma$ significance level in the next decades, in both an event-wise analysis and a more conservative statistical analysis, even with pessimistic assumptions on the spectral index of the astrophysical flux. Finally, we consider the sensitivity of future experiments to mixed production mechanisms.

Observations of debris disks have significantly improved over the past decades, both in terms of sensitivity and spatial resolution. At near-infrared wavelengths, new observing strategies and post-processing algorithms allow us to drastically improve the final images, revealing faint structures in the disks. These structures inform us about the properties and spatial distribution of the small dust particles. We present new $H$-band observations of the disk around HD 129590, which display an intriguing arc-like structure in total intensity but not in polarimetry, and propose an explanation for the origin of this arc. Assuming geometric parameters for the birth ring of planetesimals, our model provides the positions of millions of particles of different sizes to compute scattered light images. We demonstrate that if the grain size distribution is truncated or strongly peaks at a size larger than the radiation pressure blow-out size we are able to produce an arc quite similar to the observed one. If the birth ring is radially narrow, given that particles of a given size have similar eccentricities, they will have their apocenters at the same distance from the star. Since this is where the particles will spend most of their time, this results in a "apocenter pile-up" that can look like a ring. Due to more efficient forward scattering this arc only appears in total intensity observations and remains undetected in polarimetric data. This scenario requires sharp variations either in the grain size distribution or for the scattering efficiencies $Q_\mathrm{sca}$. Alternative possibilities such as a wavy size distribution and a size-dependent phase function are interesting candidates to strengthen the apocenter pile-up. We also discuss why such arcs are not commonly detected in other systems, which can mainly be explained by the fact that most parent belts are usually broad.

We investigated the snapshots of five NICER observations of the black hole transient GX 339-4 when the source transited from the hard state into the soft state during its outburst in 2021. In this paper, we focused our study on the evolution of quasi-periodic oscillations (QPOs) and noise components using power-density spectra. In addition, we derived hardness ratios comparing count rates above and below 2 keV. The evolution from the hard to the soft state was a somewhat erratic process showing several transitions between states that are dominated by top-flat noise and can show type-C QPOs; those that are dominated by red noise and can show type-B QPOs. From the parameters that we studied, we only found a strong correlation between the hardness ratio and the type of QPO observed. This implies that the appearance of type-B QPOs is related to a change in the accretion geometry of the system that also reflects in altered spectral properties. We also observed that the type-B QPO forms from or disintegrates into a broad peaked feature when the source comes out of or goes to the hard-intermediate state, respectively. This implies some strong decoherence in the process that creates this feature.

We present a new reconstruction of the Event Horizon Telescope (EHT) image of the M87 black hole from the 2017 data set. We use PRIMO, a novel dictionary-learning based algorithm that uses high-fidelity simulations of accreting black holes as a training set. By learning the correlations between the different regions of the space of interferometric data, this approach allows us to recover high-fidelity images even in the presence of sparse coverage and reach the nominal resolution of the EHT array. The black hole image comprises a thin bright ring with a diameter of $41.5\pm0.6\,\mu$as and a fractional width that is at least a factor of two smaller than previously reported. This improvement has important implications for measuring the mass of the central black hole in M87 based on the EHT images.

We used the reparameterized Near-Earth Asteroid Thermal Model to model observations of a curated set of over 4000 asteroids from the Wide-field Infrared Survey Explorer in two wavelength bands (W2-3 or W3-4) and compared the results to previous results from all four wavelength bands (W1-4). This comparison was done with the goal of elucidating unique aspects of modeling two-band observations so that any potential biases or shortcomings for planned two-band surveys (e.g., the NASA Near-Earth Object Surveyor Mission) can be anticipated and quantified. The W2-3 two-band fits usually yielded slightly smaller diameters than the four-band fits, with a median diameter difference of -10%, with the 5% and 95% quantiles of the distribution at -32% and -1.5%, respectively. We conducted similar comparisons for W3-4, in part because the longest wavelength bands are expected to provide the best two-band results. We found that the W3-4 two-band diameters are slightly larger than the four-band results, with a median diameter difference of 11% and the 5% and 95% quantiles of the distribution at -2.1% and 26%, respectively. The diameter uncertainty, obtained with bootstrap analysis, is larger by 30% and 35% (median values) for the W2-3 and W3-4 fits, respectively, than for the corresponding four-band fits. Using 23 high-quality stellar occultation diameters as a benchmark, we found that the median errors of W2-3 and W3-4 diameter estimates are -15% and +12%, respectively, whereas the median error of the four-band fits is 9.3%. Although the W2-3 and W3-4 diameters appear to have greater systematic errors and uncertainties than their four-band counterparts, two-band estimates remain useful because they improve upon diameter estimates obtained from visible photometry alone.

Dark matter is required in galaxies at galactocentric radii that are larger than the $a_0$-radius, which is where the gravitational acceleration generated by baryons of the galaxy equals the constant $a_0=1.2\times 10^{-10}$ms$^{-2}$ known as the galactic acceleration scale. It was found previously for massive early-type galaxies that the radial number-density profiles of their globular cluster (GC) systems follow broken power laws and the breaks occur at the $a_0$-radii. We have newly analyzed the distribution of GCs around galaxies in the Fornax cluster in existing photometric catalogs. We found that 1) the coincidence between $a_0$-radii and the break radii of globular cluster systems is valid for early-type galaxies of all masses and, 2) this also applies to the red and blue sub-populations of GCs separately.

[ABRIGED] Stellar feedback in high-redshift galaxies plays an important role in the re-ionization epoch of the Universe. Green Pea galaxies (GPs) postulate as favorite local laboratories. However, at their typical redshift of $z\sim0.2$, the most intimate interaction between stars and surrounding ISM cannot be disentangled. Detailed studies of Blue Compact Dwarf galaxies (BCDs) are necessary to anchor our investigations on them. We present here a study in detail UM 462, a BCD with similar properties to GPs uisng high quality optical IFS data with MUSE. Total oxygen abundance by means of the direct method is 12+$\log$(O/H)$\sim$8.02 and homogenous all over the galaxy, in stark contrast with the metallicities derived from several strong line methods. The velocity field for the ionised gas presents a velocity stratification in the area towards the north with redder velocities in the high ionisation lines and bluer velocities in the low ionisation lines. This is the only area with velocity dispersions clearly above the MUSE instrumental width, and it is surrounded by two $\sim$1 kpc-long structures nicknamed \emph{the horns}. We interpret the observational evidence in that area as a fragmented super-bubble fruit of the stellar feedback and it may constitute a preferred channel for LyC photons from the youngest generation of stars to escape. The most recent SF seems to propagate from the outer to the inner parts of the galaxy, and then from east to west. We identified a supernova remnant and Wolf-Rayet stars - as traced by the red bump - that support this picture. The direction of the propagation implies the presence of younger Wolf-Rayet stars at the maximum in H$\alpha$. The ensemble of results exemplifies the potential of 2D detailed spectroscopic studies of dwarf star-forming galaxies at high spatial resolution as key reference for similar studies on primeval galaxies.

N-body simulations are the most powerful method to study the non-linear evolution of large-scale structure. However, they require large amounts of computational resources, making unfeasible their direct adoption in scenarios that require broad explorations of parameter spaces. In this work, we show that it is possible to perform fast dark matter density field emulations with competitive accuracy using simple machine-learning approaches. We build an emulator based on dimensionality reduction and machine learning regression combining simple Principal Component Analysis and supervised learning methods. For the estimations with a single free parameter, we train on the dark matter density parameter, $\Omega_m$, while for emulations with two free parameters, we train on a range of $\Omega_m$ and redshift. The method first adopts a projection of a grid of simulations on a given basis; then, a machine learning regression is trained on this projected grid. Finally, new density cubes for different cosmological parameters can be estimated without relying directly on new N-body simulations by predicting and de-projecting the basis coefficients. We show that the proposed emulator can generate density cubes at non-linear cosmological scales with density distributions within a few percent compared to the corresponding N-body simulations. The method enables gains of three orders of magnitude in CPU run times compared to performing a full N-body simulation while reproducing the power spectrum and bispectrum within $\sim 1\%$ and $\sim 3\%$, respectively, for the single free parameter emulation and $\sim 5\%$ and $\sim 15\%$ for two free parameters. This can significantly accelerate the generation of density cubes for a wide variety of cosmological models, opening the doors to previously unfeasible applications, such as parameter and model inferences at full survey scales as the ESA/NASA Euclid mission.

We present a CO dynamical estimate of the mass of the super-massive black hole (SMBH) in three nearby early-type galaxies: NGC0612, NGC1574 and NGC4261. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3-6 observations of the $^{12}$CO(2-1) emission line with spatial resolutions of $14-58$ pc ($0.01''-0.26''$). We detect disc-like CO distributions on scales from $\lesssim200$ pc (NGC1574 and NGC4261) to $\approx10$ kpc (NGC0612). In NGC0612 and NGC1574 the bulk of the gas is regularly rotating. The data also provide evidence for the presence of a massive dark object at the centre of NGC1574, allowing us to obtain the first measure of its mass, $M_{\rm BH}=(1.0\pm0.2)\times10^{8}$ M$_{\odot}$ (1$\sigma$ uncertainty). In NGC4261, the CO kinematics is clearly dominated by the SMBH gravitational influence, allowing us to determine an accurate black hole mass of $(1.62{\pm 0.04})\times10^{9}$ M$_{\odot}$ ($1\sigma$ uncertainty). This is fully consistent with a previous CO dynamical estimate obtained using a different modelling technique. Signs of non-circular gas motions (likely outflow) are also identified in the inner regions of NGC4261. In NGC0612, we are only able to obtain a (conservative) upper limit of $M_{\rm BH}\lesssim3.2\times10^{9}$ M$_{\odot}$. This has likely to be ascribed to the presence of a central CO hole (with a radius much larger than that of the SMBH sphere of influence), combined with the inability of obtaining a robust prediction for the CO velocity curve. The three SMBH mass estimates are overall in agreement with predictions from the $M_{\rm BH}-\sigma_{\star}$ relation.

Observations of gravitational waves emitted by merging compact binaries have provided tantalising hints about stellar astrophysics, cosmology, and fundamental physics. However, the physical parameters describing the systems, (mass, spin, distance) used to extract these inferences about the Universe are subject to large uncertainties. The current method of performing these analyses requires performing many Monte Carlo integrals to marginalise over the uncertainty in the properties of the individual binaries and the survey selection bias. These Monte Carlo integrals are subject to fundamental statistical uncertainties. Previous treatments of this statistical uncertainty has focused on ensuring the precision of the inferred inference is unaffected, however, these works have neglected the question of whether sufficient accuracy can also be achieved. In this work, we provide a practical exploration of the impact of uncertainty in our analyses and provide a suggested framework for verifying that astrophysical inferences made with the gravitational-wave transient catalogue are accurate. Applying our framework to models used by the LIGO-Virgo-Kagra collaboration, we find that Monte Carlo uncertainty in estimating the survey selection bias is the limiting factor in our ability to probe narrow populations model and this will rapidly grow more problematic as the size of the observed population increases.

We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low energy X-ray properties and the larger AGN population. However some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs.

We report early-time ultraviolet (UV) and optical spectroscopy of the young, nearby Type II supernova (SN) 2022wsp obtained by the Hubble Space Telescope (HST)/STIS at about 10 and 20 days after the explosion. The SN 2022wsp UV spectra are compared to those of other well-observed Type II/IIP SNe, including the recently studied Type IIP SN 2021yja. Both SNe exhibit rapid cooling and similar evolution during early phases, indicating a common behavior among SNe II. Radiative-transfer modeling of the spectra of SN 2022wsp with the TARDIS code indicates a steep radial density profile in the outer layer of the ejecta, a supersolar metallicity, and a relatively high total extinction of E(B-V) = 0.35 mag. The early-time evolution of the photospheric velocity and temperature derived from the modeling agree with the behavior observed from other previously studied cases. The strong suppression of hydrogen Balmer lines in the spectra suggests interaction with a pre-existing circumstellar environment could be occurring at early times. In the SN 2022wsp spectra, the absorption component of the Mg II P Cygni profile displays a double-trough feature on day +10 that disappears by day +20. The shape is well reproduced by the model without fine-tuning the parameters, suggesting that the secondary blueward dip is a metal transition that originates in the SN ejecta.

It has been suggested that the observed flat rotation curves of disk galaxies can be a peculiar effect of General Relativity (GR) rather than evidence for the presence of dark matter (DM) halos in Newtonian gravity. In Ciotti (2022) the problem has been quantitatively addressed by using the well known weak-field, low-velocity gravitomagnetic limit of GR, for realistic exponential baryonic (stellar) disks. As expected, the resulting GR and Newtonian rotation curves are indistinguishable, with GR corrections at all radii of the order of $v^2/c^2\approx 10^{-6}$. Here we list some astrophysical problems that must be faced if the existence of DM halos is attributed to a misinterpretation of weak field effects of GR.

Glacial ice is used as a target material for the detection of ultra-high energy neutrinos, by measuring the radio signals that are emitted when those neutrinos interact in the ice. Thanks to the large attenuation length at radio frequencies, these signals can be detected over distances of several kilometers. One experiment taking advantage of this is the Radio Neutrino Observatory Greenland (RNO-G), currently under construction at Summit Station, near the apex of the Greenland ice sheet. These experiments require a thorough understanding of the dielectric properties of ice at radio frequencies. Towards this goal, calibration campaigns have been undertaken at Summit, during which we recorded radio reflections off internal layers in the ice sheet. Using data from the nearby GISP2 and GRIP ice cores, we show that these reflectors can be associated with features in the ice conductivity profiles; we use this connection to determine the index of refraction of the bulk ice as n=1.778 +/- 0.006.

A Synchronous Photometry Data Extraction (SPDE) program, performing indiscriminate monitors of all stars appearing at the same field of view of astronomical image, is developed by integrating several Astropy affiliated packages to make full use of time series observed by the traditional small/medium aperture ground-based telescope. The complete full-frame stellar photometry data reductions implemented for the two time series of cataclysmic variables: RX J2102.0+3359 and Paloma J0524+4244 produce 363 and 641 optimal light curves, respectively. A cross-identification with the SIMBAD finds 23 known stars, of which 16 red giant-/horizontal-branch stars, 2 W UMa-type eclipsing variables, 2 program stars, a X-ray source and 2 Asteroid Terrestrial-impact Last Alert System variables. Based on the data productions of the SPDE program, a followup Light Curve Analysis (LCA) program identifies 32 potential variable light curves, of which 18 are from the time series of RX J2102.0+3359, and 14 are from that of Paloma J0524+4244. They are preliminarily separated into periodical, transient, and peculiar types. By querying for the 58 VizieR online data catalogs, their physical parameters and multi-band brightness spanning from X-ray to radio are compiled for future analysis.

We introduce a new scheme based on the marked correlation function to probe gravity using the large-scale structure of the Universe. We illustrate our approach by applying it to simulations of the metric-variation $f(R)$ modified gravity theory and general relativity (GR). The modifications to the equations in $f(R)$ gravity lead to changes in the environment of large-scale structures that could, in principle, be used to distinguish this model from GR. Applying the Monte Carlo Markov Chain algorithm, we use the observed number density and two-point clustering to fix the halo occupation distribution (HOD) model parameters and build mock galaxy catalogues from both simulations. To generate a mark for galaxies when computing the marked correlation function we estimate the local density using a Voronoi tessellation. Our approach allows us to isolate the contribution to the uncertainty in the predicted marked correlation function that arises from the range of viable HOD model parameters, in addition to the sample variance error for a single set of HOD parameters. This is critical for assessing the discriminatory power of the method. In a companion paper we apply our new scheme to a current large-scale structure survey.

A method of near real-time detection and tracking of resident space objects (RSOs) using a convolutional neural network (CNN) and linear quadratic estimator (LQE) is proposed. Advances in machine learning architecture allow the use of low-power/cost embedded devices to perform complex classification tasks. In order to reduce the costs of tracking systems, a low-cost embedded device will be used to run a CNN detection model for RSOs in unresolved images captured by a gray-scale camera and small telescope. Detection results computed in near real-time are then passed to an LQE to compute tracking updates for the telescope mount, resulting in a fully autonomous method of optical RSO detection and tracking. Keywords: Space Domain Awareness, Neural Networks, Real-Time, Object Detection, Embedded Systems.

We have performed a series of two-dimensional two-temperature general relativistic magnetohydrodynamic simulations of magnetized accretion flows initiated from tori with different sizes and poloidal magnetic loop polarities. In these two temperature simulations, we trace the process of heating electrons through turbulence and reconnection, most of the time these electrons are trapped in plasmoids. We found that the accretion process strongly depends on the size of the magnetic loops. The accretion flows never reach the magnetically arrested (MAD) regime in small loop cases. Interaction between magnetic field with different polarities dissipates and decreases the efficiency of magneto-rotational instability. The dependency on the wavelength of the loops places a lower limit on the loop size. In the large loop cases, after reaching a quasi-steady phase, a transition from Standard And Normal Evolution (SANE) flow to MAD flow is observed. The transition of the accretion state and the transition time depends on the initial loop wavelength. The formation of plasmoids strongly depends on the size of the magnetic loops. The frequent magnetic reconnection between the magnetic loops is responsible for the formation of most of the plasmoids. For some plasmoids, Kelvin-Helmholtz and tearing instabilities are coexisting, showing another channel of plasmoid formation. The simulations present that electrons in the plasmoids are well-heated up by turbulent and magnetic reconnection. Different properties of plasmoid formation in different magnetic field configurations provide new insights for the understanding of flaring activity and electron thermodynamics in Sgr A*.

We present a short overview of the TeV-Halos objects as a discovery and a relevant contribution of the High Altitude Water \v{C}erenkov (HAWC) observatory to TeV astrophysics. We discuss history, discovery, knowledge, and the next step through a new and more detailed analysis than the original study in 2017. TeV-Halos will contribute to resolving the problem of the local positron excess observed on the Earth. To clarify the latter, understanding the diffusion process is mandatory.

We report VLA B-configuration observations of the HI 21 cm line on the close disk galaxy pair NGC 5595 and NGC 5597. At the angular resolution of the observations, $\sim7.1'' \times 4.2''$, while most of the HI 21 cm in NGC 5595 and in NGC 5597 has the same extent as the optical disk, we have detected for the first time extended structures (streamers) to the north-east (NE), and south-west (SW) of NGC 5595 with no counterparts in blue, red optical (continuum), 20 cm radio continuum, or H$\alpha$ spectral-line emission. One structure is extended by $\sim 45''$ to the NE with blue-shifted velocities, and the other by $\sim 20''$ to the SW with red-shifted velocities with respect to the systemic velocity. No HI 21 cm emission is detected from the innermost central (nuclear) regions of either galaxy. Lower angular resolution HI 21 cm imaging indicates the non-existence of any intergalactic HI 21 cm gas as tails or bridges between the two galaxies. Our new 20 cm radio continuum emission image of NGC 5597 shows a strong unresolved elongated structure at the central region, in the north-east south-west direction, very similar to the spatial location of the innermost H$\alpha$ spectral line emission. There is no 20 cm continuum emission from its north spiral arm. In NGC 5595, the 20 cm radio continuum image shows no continuum emission from the NE nor the SW extended structures with HI 21 cm emission.

We present polarization observations of the young supernova remnant (SNR) Cas A using the High-resolution Airborne Wideband Camera-Plus (HAWC+) instrument onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization map at 154 microns reveals dust grains with strong polarization fractions (5 - 30 percent), supporting previous measurements made over a smaller region of the remnant at 850 microns. The 154 microns emission and the polarization signal is coincident with a region of cold dust observed in the southeastern shell and in the unshocked central ejecta. The highly polarized far-IR emission implies the grains are large (greater than 0.14 microns) and silicate-dominated. The polarization level varies across the SNR, with an inverse correlation between the polarization degree and the intensity and smaller polarization angle dispersion for brighter SNR emission. Stronger polarization is detected between the bright structures. This may result from a higher collision rate between the gas and dust producing a lower grain alignment efficiency where the gas density is higher. We use the dust emission to provide an estimate of the magnetic field strength in Cas A using the Davis-Chandrasekhar-Fermi method. The high polarization level is direct evidence that grains are highly elongated and strongly aligned with the magnetic field of the SNR. The dust mass from the polarized region is 0.14+-0.04 Msun, a lower limit of the amount of dust present within the ejecta of Cas A. This result strengthens the hypothesis that core-collapse SNe are an important contributor to the dust mass in high redshift galaxies.

Relying on precise observations for the CPD-54 810 binary system, we investigate the robustness of the estimated age and convective core overshooting for a system with both stars on the main sequence (MS). [...] We adopt the SCEPtER pipeline, based on grids of stellar models computed for a different initial chemical composition and convective core overshooting efficiency. The base fit suggests a common age of $3.02 \pm 0.15$ Gyr, in agreement with recent literature. This estimated convective core overshooting parameter is $\beta = 0.09 \pm 0.01$, with a corresponding convective core mass $M_c = 0.059^{+0.017}_{-0.021}$ $M_{\odot}$. The robustness of these estimates were tested assuming a narrow constraint on the helium-to-metal enrichment ratio. The chemical solution of the system changes, but the age and the overshooting parameter are almost unchanged ($3.08^{+0.17}_{-0.14}$ Gyr and $0.09 \pm 0.01$). In a further test, we halved the uncertainty as to the effective temperature of both stars and again the estimated parameter shows only small variations ($3.02 \pm 0.12$ Gyr and $0.09 \pm 0.01$). This low variability suggests that the age of the system with both stars in the MS can be reliably estimated at a 5\% level, but it also indicates that the power of the investigation is probably low. [...] Despite the great increase in the observational constraints' precision, the results support the conclusions of previous theoretical works on the stellar parameter calibration with double MS star binary systems.

We explore the potential of precision cosmological data to study non-minimal dark sectors by updating the cosmological constraint on the mirror twin Higgs model (MTH). The MTH model addresses the Higgs little hierarchy problem by introducing dark sector particles. In this work, we perform a Bayesian global analysis that includes the latest cosmic shear measurement from the DES three-year survey and the Planck CMB and BAO data. In the early Universe, the mirror baryon and mirror radiation behave as dark matter and dark radiation, and their presence modifies the Universe's expansion history. Additionally, the scattering between mirror baryon and photon generates the dark acoustic oscillation process, suppressing the matter power spectrum from the cosmic shear measurement. We demonstrate how current data constrain these corrections to the $\Lambda$CDM cosmology and find that for a viable solution to the little hierarchy problem, the proportion of MTH dark matter cannot exceed about $30\%$ of the total dark matter density, unless the temperature of twin photon is less than $30\%$ of that of the standard model photon. While the MTH model is presently not a superior solution to the observed $H_0$ tension compared to the $\Lambda$CDM+$\Delta N_{\rm eff}$ model, we demonstrate that it has the potential to alleviate both the $H_0$ and $S_8$ tensions, especially if the $S_8$ tension persists in the future and approaches the result reported by the Planck SZ (2013) analysis. In this case, the MTH model can relax the tensions while satisfying the DES power spectrum constraint up to $k \lesssim 10~h\rm {Mpc}^{-1}$. If the MTH model is indeed accountable for the $S_8$ and $H_0$ tensions, we show that the future China Space Station Telescope (CSST) can determine the twin baryon abundance with a $10\%$ level precision.

We present a parameterization for the eddy diffusion profile of gas giant exoplanets based on physical phenomena and we explore how the parameterized eddy profile impacts the chemical composition, the thermal structure, the haze microphysics, and the transit spectra of 8 hot-Jupiters. Our eddy parameterization depends on the planetary intrinsic temperature (T$_{int}$ ), we thus evaluate how the increase of this parameter to values higher than those typically used ($\sim$100K) impacts the atmospheric structure and composition. Our investigation demonstrates that despite the strong impact of T$_{int}$ on the chemical composition of the deep atmosphere, the upper atmosphere is not affected for T$_{eq}$ $>$ 1300 K owing to high altitude quench levels at these conditions. Below this threshold, however, the larger atmospheric temperatures produced by increasing T$_{int}$ affect the quenched chemical composition. Our parameterization depends on two parameters, the eddy magnitude at the radiative-convective boundary (K$_0$) and the corresponding magnitude at the homopause (K$_{top}$). We demonstrate that, when using common K$_0$ and K$_{top}$ values among most of the different planet cases studied, we derive transit spectra consistent with Hubble Space Telescope observations. Moreover, our simulations show that increasing the eddy profile enhances the photochemical production of haze particles and reduces their average radius, thus providing a steeper UV-Visible slope. Finally, we demonstrate for WASP-39b that the James Webb Space Telescope observations provide improved constraints for the hazes and clouds and we show that both components seem necessary to interpret the combined transit spectrum from HST and JWST observations.

The nature of dark matter (DM) is one of the most relevant questions in modern astrophysics. We present a brief overview of recent results that inquire into a possible fermionic quantum nature of the DM particles, focusing mainly on the interconnection between the microphysics of the neutral fermions {and the macrophysical structure of galactic halos, including their formation both in the linear and non-linear cosmological regimes. We discuss the general relativistic Ruffini-Arg\"uelles-Rueda (RAR) model of fermionic DM in galaxies, its applications to the Milky Way, the possibility that the Galactic center harbors a DM core instead of a supermassive black hole (SMBH), the S-cluster stellar orbits with an in-depth analysis of the S2's orbit including precession, the application of the RAR model to other galaxy types (dwarf, elliptic, big elliptic and galaxy clusters), and universal galaxy relations. All the above focusing on the model parameters constraints, most relevant to the fermion mass. We also connect the RAR model fermions with particle physics DM candidates, self-interactions, and galactic observables constraints. The formation and stability of core-halo galactic structures predicted by the RAR model and their relation to warm DM cosmologies are also treated. Finally, we briefly discuss how gravitational lensing, dynamical friction, and the formation of SMBHs can also probe the DM nature.

Astrometric surveys can be used to constrain the stochastic gravitational wave background (SGWB) at very low frequencies. We use proper motion data provided by Gaia DR3 to fit a generic dipole+quadrupole field. We analyse several quasar-based datasets and discuss their purity and idoneity to set constraints on gravitational waves. For the cleanest dataset, we derive an upper bound on the (frequency-integrated) energy density of the SGWB $h_{70}^2\Omega_{\rm GW}\lesssim 0.087$ for $4.2\times 10^{-18}~\mathrm{Hz}\lesssim f\lesssim 1.1\times 10^{-8}~\mathrm{Hz}$. We also reanalyse previous VLBI-based data to set the constraint $h_{70}^2\Omega_{\rm GW}\lesssim 0.024$ for $5.8\times 10^{-18}~\mathrm{Hz}\lesssim f\lesssim 1.4\times 10^{-9}~\mathrm{Hz}$ under the same formalism, standing as the best astrometric constraint on GWs. Based on our results, we discuss the potential of future Gaia data releases to impose tighter constraints.

Common envelope events have been associated with the formation of a planetary nebulae since its proposition more than forty five years ago. However, until recently there have been doubts as to whether a common envelope while the donor is ascending the red giant branch, rather than the subsequent asymptotic red giant branch, would result in a planetary nebula. There is now strong theoretical and observational evidence to suggest that some planetary nebulae are, indeed, the products of common envelope phases which occurred while the nebular progenitor was on the red giant branch. The characterisation of these systems is challenging but has the potential to reveal much about the common envelope -- a critical evolutionary phase in the formation of a plethora of interesting astrophysical phenomena.

We use the latest constraints on the population of stellar origin binary black holes (SOBBH) from LIGO/Virgo/KAGRA (LVK) observations, to estimate the stochastic gravitational wave background (SGWB) they generate in the frequency band of LISA. We account for the faint and distant binaries, which contribute the most to the SGWB, by extending the merger rate at high redshift assuming it tracks the star formation rate. We adopt different methods to compute the SGWB signal: an analytical evaluation, Monte Carlo sums over SOBBH population realisations, and a method that accounts for the role of the detector by simulating LISA data and iteratively removing resolvable signals until only the confusion noise is left, allowing for the extraction of both the expected SGWB and the number of resolvable SOBBHs. Since the latter are few for SNR thresholds larger than five, we confirm that the spectral shape of the SGWB in the LISA band follows the analytical prediction of a power law. We infer the probability distribution of the SGWB amplitude from the LVK GWTC-3 posterior of the binary population model; its interquartile range of $h^2\Omega_\mathrm{GW}(f=3\times10^{-3}\,\mathrm{Hz}) \in [5.65,\,11.5]\times10^{-13}$ is in agreement with most previous estimates. We perform a MC analysis to assess LISA's capability to detect and characterise this signal. Accounting for both the instrumental noise and the galactic binaries foreground, with four years of data, LISA will be able to detect the SOBBH SGWB with percent accuracy, narrowing down the uncertainty on the amplitude by one order of magnitude with respect to the range of possible amplitudes inferred from the population model. A measurement of this signal by LISA will help to break the degeneracy among some of the population parameters, and provide interesting constraints, in particular on the redshift evolution of the SOBBH merger rate.

The atmospheric composition of planets is determined by the chemistry of the disks in which they form. Studying the gas-phase molecular composition of disks thus allows us to infer what the atmospheric composition of forming planets might be. Recent observations of the IRS 48 disk have shown that (asymmetric) dust traps can directly impact the observable chemistry, through radial and vertical transport, and the sublimation of ices. The asymmetric HD 142527 disk provides another good opportunity to investigate the role of dust traps in setting the disk's chemical composition. In this work, we use archival ALMA observations of the HD 142527 disk to obtain an as large as possible molecular inventory, which allows us to investigate the possible influence of the asymmetric dust trap on the disk's chemistry. We present the first ALMA detections of [C I], 13C18O, DCO+, H2CO and additional transition of HCO+ and CS in this disk. In addition, we have acquired upper limits for non-detected species such as SO and CH3OH. For the majority of the observed molecules, a decrement in the emission at the location of the dust trap is found. For the main CO isotopologues continuum over-subtraction likely causes the observed asymmetry, while for CS and HCN we propose that the observed asymmetries are likely due to shadows cast by the misaligned inner disk. As the emission of the observed molecules is not co-spatial with the dust trap and no SO or CH3OH are found, thermal sublimation of icy mantles does not appear to play a major role in changing the gas-phase composition of the outer disk in HD 142527 disk. Using our observations of 13C18O and DCO+ and a RADMC-3D model, we determine the CO snowline to be located beyond the dust traps, favouring cold gas-phase formation of H2CO, rather than the hydrogenation of CO-ice and subsequent sublimation.

Detecting planetary signatures in radial velocity time-series of young stars is challenging due to their inherently strong stellar activity. However, it is possible to learn information about the properties of the stellar signal by using activity indicators measured from the same stellar spectra used to extract radial velocities. In this manuscript, we present a reanalysis of spectroscopic HARPS data of the young star K2-233, which hosts three transiting planets. We perform a multidimensional Gaussian Process regression on the radial velocity and the activity indicators to characterise the planetary Doppler signals. We demonstrate, for the first time on a real dataset, that the use of a multidimensional Gaussian Process can boost the precision with which we measure the planetary signals compared to a one-dimensional Gaussian Process applied to the radial velocities alone. We measure the semi-amplitudes of K2-233 b, c, and d as 1.31(-0.74)(+0.81), 1.81(-0.67)(+0.71), and 2.72(-0.70)(+0.66) m/s, which translates into planetary masses of 2.4(-1.3)(+1.5), 4.6(-1.7)(+1.8), and 10.3(-2.6)(+2.4), respectively. These new mass measurements make K2-233 d a valuable target for transmission spectroscopy observations with JWST. K2-233 is the only young system with two detected inner planets below the radius valley and a third outer planet above it. This makes it an excellent target to perform comparative studies, to inform our theories of planet evolution, formation, migration, and atmospheric evolution.

As one of the possible extensions of Einstein's General Theory of Relativity, it has been recently suggested that the presence of spacetime torsion could solve problems of the very early and the late-time universe undergoing accelerating phases. In this paper, we use the latest observations of high-redshift data, coming from multiple measurements of quasars and baryon acoustic oscillations, to phenomenologically constrain such cosmological model in the framework of Einstein-Cartan (EC) endowed with spacetime torsion. Such newly compiled quasar datasets in the cosmological analysis is crucial to this aim, since it will extend the Hubble diagram to high-redshift range in which predictions from different cosmologies can be distinguished. Our results show that out of all the candidate models, the torsion plus cosmological constant model is strongly favoured by the current high-redshift data, where torsion itself would be expected to yield the current cosmic acceleration. Specially, in the framework of Friedmann-like cosmology with torsion, the determined Hubble constant is in very good agreement with that derived from the Planck 2018 CMB results. On the other hand, our results are compatible with zero spatial curvature and there is no significant deviation from flat spatial hypersurfaces. Finally, we check the robustness of high-redshift observations by placing constraints on the torsion parameter $\alpha$, which is strongly consistent with other recent works focusing on torsion effect on the primordial helium-4 abundance.

Cosmic strings, if they exist, source nonlinear and non-Gaussian perturbations all the way back to the time of equal matter and radiation (and earlier). Here, we compute the mass function of halos seeded by a scaling distribution of cosmic string loops, and we compare the results with the predictions of the standard Gaussian $\Lambda$CDM model. Assuming a simple linear relation between stellar mass and halo mass, we also compute the stellar mass function. The contribution of cosmic strings dominates at sufficiently high redshifts $z > z_c$ where $z_c$ depends on the mass of the halo and on the mass per unit length $\mu$ of the strings and is of the order $z_c \sim 12$ for $G\mu = 10^{-8}$. We find that strings with this value of $G\mu$ can explain the preliminary JWST data on the high redshift stellar mass density. Based on an extreme value statistic, we find that the mass of the heaviest expected string-seeded galaxy for the current JWST sky coverage is compatible with the heaviest detected galaxy. Given the uncertainties in the interpretation of the JWST data, we discuss predictions for higher redshift observations.

We study the abundance patterns and the radial gradients of s-process elements (Y, Zr, Ba, La and Ce), r-process elements (Eu) and mixed-process elements (Mo, Nd and Pr) in the Galactic thin disc by means of a detailed two-infall chemical evolution model for the Milky Way with state-of-the-art nucleosynthesis prescriptions. We consider r-process nucleosynthesis from merging neutron stars (MNS), magneto-rotational supernovae (MR-SNe) and s-process synthesis from low- and intermediate- mass stars (LIMS) and rotating massive stars. The predictions of our model are compared with data from the sixth data release of the Gaia-ESO survey, from which we consider 62 open clusters with age > 0.1 Gyr and 1300 Milky Way disc field stars. We conclude that: i) the [Eu/Fe] vs. [Fe/H] is reproduced by both a prompt and a delayed source, but the quick source completely dominates the Eu production; ii) rotation in massive stars contribute substantially to the s-process elements of the first peak, but MNS and MR-SNe are necessary in order to reproduce the observations; iii) due to the adopted yields, our model overpredicts Pr and underpredicts Nd, while the [Mo/Fe] vs. [Fe/H] is nicely reproduced. For the radial gradients, we conclude that: i) our predicted slope of the [Fe/H] gradient is in agreement with the one observed in open clusters by Gaia-ESO and other high-resolution spectroscopic surveys. ii) The predicted slope of the [Eu/H] radial gradient is steeper than the observed one, independently on how quick the production of Eu is. We discuss the possible causes of this discrepancy in terms of both different Galaxy formation scenarios and stellar radial migration effects. iii) For all the elements belonging to the second s-process peak (Ba, La, Ce) as well as for Pr, we predict a plateau at low Galactocentric distances, which is probably due to the enhanced enrichment from LIMS in the inner regions.

Very long baseline radio interferometry (VLBI) with ground-based observatories is limited by the size of Earth, the geographic distribution of antennas, and the transparency of the atmosphere. In this whitepaper, we present Capella, a tentative design of a space-only VLBI system. Using four small (<500 kg) satellites on two orthogonal polar low-Earth orbits, and single-band heterodyne receivers operating at frequencies around 690 GHz, the interferometer is able to achieve angular resolutions of approximately 7 microarcsec. Within a total observing time of three days, a near-complete uv plane coverage can be reached, with a 1-sigma point source sensitivity as good as about 6~mJy for an instantaneous bandwidth of 1 GHz. The required downlink data rates of >10 Gbps can be reached through near-infrared laser communication; depending on the actual downlink speed, one or multiple ground communication stations are necessary. We note that all key technologies required for the Capella system are already available, some of them off-the-shelf. Data can be correlated using dedicated versions of existing Fourier transform (FX) software correlators; dedicated routines will be needed to handle the effects of orbital motion, including relativistic corrections. With the specifications assumed in this whitepaper, Capella will be able to address a range of science cases, including: photon rings around supermassive black holes; the acceleration and collimation zones of plasma jets emitted from the vicinity of supermassive black holes; the chemical composition of accretion flows into active galactic nuclei through observations of molecular absorption lines; mapping supermassive binary black holes; the magnetic activity of stars; and nova eruptions of symbiotic binary stars - and, like any substantially new observing technique, has the potential for unexpected discoveries.

(Abridged) The star-formation rate (SFR) quantitatively describes the star-formation process in galaxies. Current ways to calibrate this rate do not usually employ observational methods accounting for the low-mass end of stellar populations as their signatures are too weak. Accessing the bulk of protostellar activity within galactic star-forming regions can be achieved by tracing signposts of ongoing star formation. One such signpost is molecular outflows, which are bright in molecular emission. We propose to utilize the protostellar outflow emission as a tracer of the SFR. In this work, we introduce a novel version of the galaxy-in-a-box model, which can be used to relate molecular emission from star formation in galaxies with the SFR. We measured the predicted para-H2O emission at 988 GHz and corresponding SFRs for galaxies with LFIR = $10^8$ - $10^{11}$ L$_\odot$ in a distance-independent manner, and compared them with expectations from observations. We evaluated the derived results by varying the star formation efficiency, the free-fall time scaling factor, and the initial mass function. For the chosen H2O transition, relying on the current Galactic observations and star formation properties, we are underestimating the total galactic emission, while overestimating the SFRs, particularly for more starburst-like configurations. The current version of the galaxy-in-a-box model accounts for a limited number of processes and configurations, that is, it focuses on ongoing star formation in massive young clusters in a spiral galaxy. Therefore, the inferred results, which underestimate the emission and overestimate the SFR, are not surprising: known sources of emission are not included in the model. To improve the results, the next version of the model needs to include a more detailed treatment of the entire galactic ecosystem and other processes that would contribute to the emission.

PSR J1528-3146 is a 60.8 ms pulsar orbiting a heavy white dwarf (WD) companion, with an orbital period of 3.18 d. This work aimed at characterizing the pulsar's astrometric, spin and orbital parameters by analyzing timing measurements conducted at the Parkes, MeerKAT and Nan\c{c}ay radio telescopes over almost two decades. The measurement of post-Keplerian perturbations to the pulsar's orbit can be used to constrain the masses of the two component stars of the binary, and in turn inform us on the history of the system. We analyzed timing data from the Parkes, MeerKAT and Nan\c{c}ay radio telescopes collected over $\sim$16 yrs, obtaining a precise rotation ephemeris for PSR J1528-3146. A Bayesian analysis of the timing data was carried out to constrain the masses of the two components and the orientation of the orbit. We further analyzed the polarization properties of the pulsar, in order to constrain the orientations of the magnetic axis and of the line-of-sight with respect to the spin axis. We measured a significant rate of advance of periastron for the first time, and put constraints on the Shapiro delay in the system and on the rate of change of the projected semi-major axis of the pulsar's orbit. The Bayesian analysis yielded measurements for the pulsar and companion masses of respectively $M_p = 1.61_{-0.13}^{+0.14}$ M$_\odot$ and $M_c = 1.33_{-0.07}^{+0.08}$ M$_\odot$ (68\% C.L.), confirming that the companion is indeed massive. This companion mass as well as other characteristics of PSR J1528$-$3146 make this pulsar very similar to PSR J2222-0137, a 32.8 ms pulsar orbiting a WD whose heavy mass ($\sim 1.32$ M$_\odot$) was unique among pulsar-WD systems until now. Our measurements therefore suggest common evolutionary scenarios for PSRs J1528-3146 and J2222-0137.

GRB 211211A is a rare burst with a genuinely long duration, yet its prominent kilonova association provides compelling evidence that this peculiar burst was the result of a compact binary merger. However, the exact nature of the merging objects, whether they were neutron star pairs, neutron star-black hole systems, or neutron star-white dwarf systems, remains unsettled. This Letter delves into the rarity of this event and the possibility of using current and next-generation gravitational wave detectors to distinguish between the various types of binary systems. Our research reveals an event rate density of $\gtrsim 5.67^{+13.04}_{-4.69} \times 10^{-3}\ \rm Gpc^{-3} yr^{-1}$ for GRB 211211A-like GRBs, which is significantly smaller than that of typical long and short GRB populations. We further calculated that if the origin of GRB 211211A is a result of a neutron star-black hole merger, it would be detectable with a significant signal-to-noise ratio, given the LIGO-Virgo-KAGRA designed sensitivity. On the other hand, a neutron star-white dwarf binary would also produce a considerable signal-to-noise ratio during the inspiral phase at decihertz and is detectable by next-generation space-borne detectors DECIGO and BBO. However, to detect this type of system with millihertz space-borne detectors like LISA, Taiji, and TianQin, the event must be very close, approximately 3 Mpc in distance or smaller.

An anti-correlation between the central density of the dark matter halo ($\rho_{150,\ {\rm DM}}$) and the pericentric distances ($r_{p}$) of the Milky Way's (MW's) dwarf spheroidal galaxies (dSphs) has been reported in the literature. The existence and origin of such anti-correlation is however controversial, one possibility being that only the densest dSphs can survive the tidal field towards the centre of our Galaxy. In this work, we place particular emphasis on quantifying the statistical significance of such anti-correlation, by using available literature data in order to explore its robustness under different assumptions on the MW gravitational potential, and for various derivations of $\rho_{150}$ and $r_{p}$. We consider models in which the MW is isolated and has a low ($8.8\times10^{11}\,M_{\odot}$) and high ($1.6\times10^{12}\, M_{\odot}$) halo mass, respectively, as well as configurations in which the MW's potential is perturbed by a Large Magellanic Cloud (LMC) infall. We find that, while data generally support models in which the dSphs' central DM density decreases as a function of their pericentric radius, this anti-correlation is statistically significant at $3\sigma$ level only in $\sim$12$\%$ of the combinations of $\rho_{150}$ and $r_{p}$ explored. Moreover, including the impact of the LMC's infall onto the MW weakens or even washes away this anti-correlation, with respect to models in which the MW is isolated. Our results suggest that the strength and existence of such anti-correlation is still debatable: exploring it with high-resolution simulations including baryonic physics and different DM flavours will help us to understand its emergence.

Bipolar Magnetic Regions (BMRs) provide crucial information about solar magnetism. They exhibit varying morphology and magnetic properties throughout their lifetime, and studying these properties can provide valuable insights into the workings of the solar dynamo. The majority of previous studies have counted every detected BMR as a new one and have not been able to study the full life history of each BMRs. To address this issue, we have developed an Automatic Tracking Algorithm (AutoTAB) for BMRs, that tracks the BMRs for their entire lifetime or throughout their disk passage. AutoTAB uses the binary maps of detected BMRs to automatically track the regions. This is done by differentially rotating the binary maps of the detected regions and checking for overlaps between them. In this first article of this project, we provide a detailed description of the working of the algorithm and evaluate its strengths and weaknesses. We also compare its performance with other existing tracking techniques. AutoTAB excels in tracking even for the small features and it successfully tracks 9152 BMRs over the last two solar cycles (1996-2020), providing a comprehensive dataset that depicts the evolution of various properties for each tracked region. The tracked BMRs follow familiar properties of solar cycles except for these small BMRs that appear at all phases of the solar cycle and show weak latitudinal dependency, which is represented through the butterfly diagram. Finally, we discuss the possibility of adapting our algorithm to other datasets and expanding the technique to track other solar features in the future.

We combine the Planck-SZ2 galaxy cluster catalogue with near-infrared photometry of galaxies from the VISTA Hemisphere Survey to identify candidate brightest cluster galaxies (BCGs) in 306 massive clusters in the Southern skies at redshifts of $z>0.1$. We find that 91% of these clusters have at least one candidate BCG within the 95% confidence interval on the cluster centers quoted by the Planck collaboration, providing reassurance that our analyses are statistically compatible, and find 92% to be reasonable candidates following a manual inspection. We make our catalog publicly available to assist colleagues interested in multi-wavelength studies of cluster cores, and the search for gravitationally lensed explosive transients in upcoming surveys including the Legacy Survey of Space and Time by the Vera C. Rubin Observatory.

Multiwavelength studies indicate that nuclear activity and bulge properties are closely related, but the details remain unclear. To study this further, we combine $Hubble~Space~Telescope$ bulge structural and photometric properties with 1.5 GHz, $e$-MERLIN nuclear radio continuum data from the LeMMINGs survey for a large sample of 173 `active' galaxies (LINERs and Seyferts) and `inactive' galaxies (H IIs and absorption line galaxies, ALGs). Dividing our sample into active and inactive, they define distinct (radio core luminosity)$-$(bulge mass), L_R,core-M_*,bulge, relations, with a mass turnover at M_*, bulge ~ 10^(9.8 +- 0.3) M_sun (supermassive black hole mass M_BH ~ 10^(6.8 +- 0.3) M_sun), which marks the transition from AGN-dominated nuclear radio emission in more massive bulges to that mainly driven by stellar processes in low-mass bulges. None of our 10/173 bulgeless galaxies host an AGN. The AGN fraction increases with increasing M_*, bulge such that f_optical_AGN $\propto$ M_*,bulge^(0.24 +- 0.06) and f_radio_AGN $\propto$ M_*,bulge^(0.24 +- 0.05). Between M_*,bulge ~ 10^8.5 and 10^11.3 M_sun, f_optical_AGN steadily rises from 15 +- 4 to 80 +- 5 per cent. We find that at fixed bulge mass, the radio loudness, nuclear radio activity and the (optical and radio) AGN fraction exhibit no dependence on environment. Radio-loud hosts preferentially possess an early-type morphology than radio-quiet hosts, the two types are however indistinguishable in terms of bulge S\'ersic index and ellipticity, while results on the bulge inner logarithmic profile slope are inconclusive. We finally discuss the importance of bulge mass in determining the AGN triggering processes, including potential implications for the nuclear radio emission in nearby galaxies.

We present pure spectroscopic constraints on the UV luminosity functions and cosmic star formation rate (SFR) densities from 25 galaxies at $z_\mathrm{spec}=8.61-13.20$. By reducing the JWST/NIRSpec spectra taken in multiple programs of ERO, ERS, GO, and DDT with our analysis technique, we independently confirm 16 galaxies at $z_\mathrm{spec}=8.61-11.40$ including new redshift determinations, and a bright interloper at $z_\mathrm{spec}=4.91$ that was claimed as a photometric candidate at z~16. In conjunction with nine galaxies at redshifts up to $z_\mathrm{spec}=13.20$ in the literature, we make a sample of 25 spectroscopically-confirmed galaxies in total and carefully derive the best estimates and lower limits of the UV luminosity functions. These UV luminosity function constraints are consistent with the previous photometric estimates within the uncertainties and indicate mild redshift evolution towards z~12 showing tensions with some theoretical models of rapid evolution. With these spectroscopic constraints, we obtain firm lower limits of the cosmic SFR densities and spectroscopically confirm a high SFR density at z~12 beyond the constant star-formation efficiency models, which supports earlier claims from the photometric studies. While there are no spectroscopically-confirmed galaxies with very large stellar masses violating the $\Lambda$CDM model due to the removal of the bright interloper, we confirm star-forming galaxies at $z_\mathrm{spec}=11-13$ with stellar masses much higher than model predictions. Our results indicate possibilities of high star-formation efficiency (>5%), hidden AGN, top-heavy initial mass function (possibly with Pop-III), and large scatter/variance. Having these successful and unsuccessful spectroscopy results, we suggest observational strategies for efficiently removing low redshift interlopers for future JWST programs.

We present hydrodynamic and magnetohydrodynamic (MHD) simulations of sub galactic regions including photoionising and supernova feedack. We aim to improve the initial conditions of our region extraction models by including an initial population of stars. We also investigate the reliability of extracting regions in simulations, and show that with a good choice of region, results are comparable with using a larger region for the duration of our simulations. Simulations of star formation on molecular cloud scales typically start with a turbulent cloud of gas, from which stars form and then undergo feedback. In reality, a typical cloud or region within a galaxy may already include, or reside near some population of stars containing massive stars undergoing feedback. We find the main role of a prior population is triggering star formation, and contributing to gas dynamics. Early time supernova from the initial population are important in triggering new star formation and driving gas motions on larger scales above 100 pc, whilst the ionising feedback contribution from the initial population has less impact, since many members of the initial population have cleared out gas around them in the prior model. In terms of overall star formation rates though, the initial population has a relatively small effect, and the feedback does not for example suppress subsequent star formation. We find that MHD has a relatively larger impact than initial conditions, reducing the star formation rate by a factor of 3 at later times.

The self-interacting dark matter can affect various cosmological processes. Such interactions can be number conserving (\emph{e.g.} $2 \rightarrow 2$) or number violating (\emph{e.g.} $3 \rightarrow 2,\,4 \rightarrow 2$ etc.). The latter processes where three (or more) dark matter particles undergo self-annihilation/scattering to produce less number of dark matter is termed as ``Cannibalism'' process. In this work, the self-interaction of dark matter and the strength of such interactions are investigated in the light of experimental results of the global 21-cm spectrum of neural hydrogen from the era of cosmic dawn. From the present work, it appears that $2\rightarrow 2$ process is much more dominant over the $3\rightarrow 2$ process. It is also found that such interactions affect the dark matter-baryon elastic scattering cross-section. The study also indicates the presence of multi component dark matter of different mass range in the Universe.

The Tip of the Red Giant Branch (TRGB) provides a luminous standard candle for constructing distance ladders to measure the Hubble constant. In practice its measurements via edge-detection response (EDR) are complicated by the apparent fuzziness of the tip and the multi-peak landscape of the EDR. As a result, it can be difficult to replicate due to a case-by-case measurement process. Previously we optimized an unsupervised algorithm, Comparative Analysis of TRGBs (CATs), to minimize the variance among multiple halo fields per host without reliance on individualized choices, achieving state-of-the-art $\sim$ $<$ 0.05 mag distance measures for optimal data. Further, we found an empirical correlation at 5$\sigma$ confidence in the GHOSTS halo survey between our measurements of the tip and their contrast ratios (ratio of stars 0.5 mag just below and above the tip), useful for standardizing the apparent tips at different host locations. Here, we apply this algorithm to an expanded sample of SN Ia hosts to standardize these to multiple fields in the geometric anchor, NGC 4258. In concert with the Pantheon$+$ SN Ia sample, this analysis produces a (baseline) result of $H_0= 73.22 \pm 2.06$ km/s/Mpc. The largest difference in $H_0$ between this and similar studies employing the TRGB derives from corrections for SN survey differences and local flows used in most recent SN Ia compilations but which were absent in earlier studies. SN-related differences total $\sim$ 2.0 km/s/Mpc. A smaller share, $\sim$ 1.4 km/s/Mpc, results from the inhomogeneity of the TRGB calibration across the distance ladder. We employ a grid of 108 variants around the optimal TRGB algorithm and find the median of variants is $72.94\pm1.98$ km/s/Mpc with an additional uncertainty due to algorithm choices of 0.83 km/s/Mpc. None of these TRGB variants result in $H_0$ less than 71.6 km/s/Mpc.

The tip of the red giant branch (TRGB) is a standard candle that can be used to help refine the determination of the Hubble constant. $Gaia$ Data Release 3 (DR3) provides synthetic photometry constructed from low-resolution BP/RP spectra for Milky Way field stars that can be used to directly calibrate the luminosity of the TRGB in the Johnson-Cousins I band, where the TRGB is least sensitive to metallicity. We calibrate the TRGB luminosity using a two-dimensional maximum likelihood algorithm with field stars and $Gaia$ synthetic photometry and parallaxes. For a high-contrast and low-contrast break (characterized by the values of the contrast parameter $ R$ or the magnitude of the break $ \beta $), we find $M^{TRGB}_I$ =$-4.02$ and $-3.92$ mag respectively, or a midpoint of $-3.970$ $^{+0.042} _{-0.024}$ (sys) $\pm$ $0.062$ (stat) mag. This measurement improves upon the TRGB measurement from Li et al. (2022), as the higher precision photometry based on $ Gaia $ DR3 allows us to constrain two additional free parameters of the luminosity function. We also investigate the possibility of using $Gaia$ DR3 synthetic photometry to calibrate the TRGB luminosity with $\omega$ Centauri, but find evidence of blending within the inner region for cluster member photometry that precludes accurate calibration with $Gaia$ DR3 photometry. We instead provide an updated TRGB measurement of $m^{TRGB}_I$ = $ 9.82 \pm 0.04$ mag in $\omega$ Centauri using ground-based photometry from the most recent version of the database described in Stetson et al. (2019), which gives $M^{TRGB}_I$ = $-3.97$ $\pm$ $0.04$ (stat) $\pm$ 0.10 (sys) mag when tied to the $Gaia$ EDR3 parallax distance from the consensus of Vasiliev & Baumgardt (2021), Soltis et al. (2021), and Ma\'{i}z Apell\'{a}niz et al. (2022a).

In this work, we present a simple interpolation methodology for spectroscopic time series, based on conventional interpolation techniques (random forests) implemented in widely-available libraries. We demonstrate that our existing library of simulations is sufficient for training, producing interpolated spectra that respond sensitively to varied ejecta parameter, post-merger time, and viewing angle inputs. We compare our interpolated spectra to the AT2017gfo spectral data, and find parameters similar to our previous inferences using broadband light curves. However, the spectral observations have significant systematic short-wavelength residuals relative to our models, which we cannot explain within our existing framework. Similar to previous studies, we argue that an additional blue component is required. We consider a radioactive heating source as a third component characterized by light, slow-moving, lanthanide-free ejecta with $M_{\rm th} = 0.003~M_\odot$, $v_{\rm th} = 0.05$c, and $\kappa_{\rm th} = 1$ cm$^2$/g. When included as part of our radiative transfer simulations, our choice of third component reprocesses blue photons into lower energies, having the opposite effect and further accentuating the blue-underluminosity disparity in our simulations. As such, we are unable to overcome short-wavelength deficits at later times using an additional radioactive heating component, indicating the need for a more sophisticated modeling treatment.

Cosmological solitonic objects such as monopoles, cosmic strings, domain walls, oscillons and Q-balls often appear in theories of the early Universe. We demonstrate that such scenarios are generically accompanied by a novel production source of gravitational waves stemming from soliton isocurvature perturbations. The resulting induced universal gravitational waves (UGWs) reside at lower frequencies compared to gravitational waves typically associated with soliton formation. We show that UGWs from axion-like particle (ALP) oscillons, originating from ALP misalignment, extend the frequency range of produced gravitational waves by more than two orders of magnitude regardless of the ALP mass and decay constant and can be observable in upcoming gravitational wave experiments. UGWs open a new route for gravitational wave signatures in broad classes of cosmological theories.

The tidal response of compact objects in an inspiraling binary system is measured by a set of tidal Love and dissipation numbers imprinted in the gravitational waveforms. While a four-dimensional black hole in vacuum within General Relativity has vanishing Love numbers, a black hole in alternative theories of gravity can acquire non-vanishing Love numbers. The dissipation numbers may quantify Planckian corrections at the horizon scale. These properties will allow a test of classical theories of gravity in the strong-field regime with gravitational-wave observation. Since black holes are not in the exact vacuum environment in astrophysical situations, the following question arises: can the environment affect the tidal response? In this paper, we investigate the stability of the tidal response of a Schwarzschild black hole for frequency-dependent tidal-field perturbations against a small modification of the background. Our analysis relies on the scattering theory, which overcomes difficulties in defining the relativistic tidal Love numbers. The tidal Love and dissipation numbers can be extracted from the phase shift of sufficiently low-frequency scattering waves. We show that the tidal Love numbers are sensitive to the property of the modification. Therefore, we need careful consideration of the environment around the black hole in assessing the deviation of the underlying theory of gravity from General Relativity with the Love numbers. The modification has less impact on the dissipation numbers, indicating that quantifying the existence of the event horizon with them is not spoiled. We also demonstrate that in a composite system, i.e., a compact object with environmental effects, the Love and dissipation numbers are approximately determined by the sum of the numbers of each component.

We entertain the novel possibility that long range forces may lead to violations of accidental symmetries, in particular baryon number. Employing an ultralight scalar, with a mass $\ll$ eV, we illustrate that this scenario can lead to vastly disparate nucleon lifetimes, in different astronomical objects. Such a long range interaction can yield a number of potentially observable effects, such as a flux of neutrinos at $\gtrsim 10$ MeV from the Sun and heating of old neutron stars. We examine the prospects for constraining this scenario, with current and future astrophysical data, and find that neutron star heating provides the strongest present and near term bounds. Simple extensions of our setup allow for the ultralight scalar to constitute the dark matter of the Universe. This suggests that matter-enhanced baryon number violation can be a signal of ultralight dark matter, which has apparently been overlooked, so far.

Dark matter (DM) could be a nonthermal relic that freezes in from extremely weak, sub-Hubble annihilation and decay of Standard Model (SM) particles. The case of Dirac DM freezing in via a dark photon mediator is a well-studied benchmark for DM direct detection experiments. Here, we extend prior work to take into account the possibility that DM is pseudo-Dirac with a small mass splitting. If the mass splitting is greater than twice the electron mass but less than the dark photon mass, there will be distinct cosmological signatures. The excited state $\chi_2$ is initially produced in equal abundance to the ground state $\chi_1$. Subsequently, the excited state population decays at relatively late cosmological times, primarily via the three-body process $\chi_2 \rightarrow \chi_1 e^+ e^-$. This process injects energetic electrons into the ambient environment, providing observable signatures involving Big Bang nucleosynthesis, cosmic microwave background spectral distortions and anisotropies, and the Lyman-$\alpha$ forest. Furthermore, the ground state particles that are populated from the three-body decay receive a velocity kick, with implications for DM clustering on small scales. We find that cosmological probes and accelerator experiments are highly complementary, with future coverage of much of the parameter space of the model.

When choosing between competing symbolic models for a data set, a human will naturally prefer the "simpler" expression or the one which more closely resembles equations previously seen in a similar context. This suggests a non-uniform prior on functions, which is, however, rarely considered within a symbolic regression (SR) framework. In this paper we develop methods to incorporate detailed prior information on both functions and their parameters into SR. Our prior on the structure of a function is based on a $n$-gram language model, which is sensitive to the arrangement of operators relative to one another in addition to the frequency of occurrence of each operator. We also develop a formalism based on the Fractional Bayes Factor to treat numerical parameter priors in such a way that models may be fairly compared though the Bayesian evidence, and explicitly compare Bayesian, Minimum Description Length and heuristic methods for model selection. We demonstrate the performance of our priors relative to literature standards on benchmarks and a real-world dataset from the field of cosmology.

We argue that, as long as relativistic quantum particles are in point, the variable $y=E/E_p$ of the rainbow functions pair $g_{_{0}} (y)$ and $g_{_{1}} (y)$ should be fine tuned into $y=|E|/E_p$, where $E_p$ is the Planck's energy scale. Otherwise, the rainbow functions will be only successful to describe the rainbow gravity effect on relativistic quantum particles and the anti-particles will be left unfortunate. Under such fine tuning, we consider Klein-Gordon (KG) particles in cosmic string rainbow gravity spacetime in a non-uniform magnetic field (i.e., $\mathbf{B}=\mathbf{\nabla }\times \mathbf{A}=\frac{3}{2}B_{\circ }r\,\hat{z}$ ). Then we consider KG-particles in cosmic string rainbow gravity spacetime in a uniform magnetic field (i.e., $\mathbf{B}=\mathbf{\nabla }\times \mathbf{A}=\frac{1}{2}B_{\circ }\,\hat{z}$ ). Whilst the former effectively yields KG-oscillators, the later effectively yields KG-Coulombic particles. We report on the effects of rainbow gravity on both KG-oscillators and Coulombic particles using four pairs of rainbow functions: (i) $% g_{_{0}}\left( y\right) =1$, $g_{_{1}}\left( y\right) =\sqrt{1-\epsilon y^{2}% }$, (ii) $g_{_{0}}\left( y\right) =1$, $g_{_{1}}\left( y\right) =\sqrt{% 1-\epsilon y}$, (iii) $g_{_{0}}\left( y\right) =g_{_{1}}\left( y\right) =\left( 1-\epsilon y\right) ^{-1}$, and (iv) $g_{_{0}}\left( y\right) =\left( e^{\epsilon y}-1\right) /\epsilon y$, $g_{_{1}}\left( y\right) =1$, where $y=|E|/E_p$ and $\epsilon$ is the rainbow parameter. It is interesting to report that, all KG particles' and anti-particles' energies are symmetric about $E=0$ value (a natural relativistic quantum mechanical tendency), and a phenomenon of energy states to fly away and disappear from the spectrum is observed for the rainbow functions pair (iii) at $\gamma=\epsilon m/E_p=1$.

We investigate the effects of the low-scale cosmological first-order phase transitions on the neutrino decoupling and constrain the PT parameters with the cosmological observations of big bang nucleosynthesis and cosmic microwave background. We consider the phase transitions that occur at the MeV-scale which can produce stochastic gravitational wave background to be probed by pulsar timing array experiments. We find that the phase transition can modify the effective number of neutrinos and the primordial nucleosynthesis. In turn, the cosmological observations can exclude slow and strong phase transitions around the MeV scale.

We study the evolution of magnetic fields coupled with chiral fermion asymmetry in the framework of chiral magnetohydrodynamics with zero initial total chirality. The initial magnetic field has a turbulent spectrum peaking at a certain characteristic scale and is fully helical with positive helicity. The initial chiral chemical potential is spatially uniform and negative. We consider two opposite cases where the ratio of the length scale of the chiral plasma instability (CPI) to the characteristic scale of the turbulence is smaller and larger than unity. These initial conditions might be realized in cosmological models such as certain types of axion inflation. The magnetic field and chiral chemical potential evolve with inverse cascading in such a way that the magnetic helicity and chirality cancel each other at all times. The CPI time scale is found to determine mainly the time when the magnetic helicity spectrum attains negative values at high wave numbers. The turnover time of the energy-carrying eddies, on the other hand, determines the time when the peak of the spectrum starts to shift to smaller wave numbers via an inverse cascade. The onset of helicity decay is determined by the time when the chiral magnetic effect becomes efficient at the peak of the initial magnetic energy spectrum. When spin flipping is important, the chiral chemical potential vanishes and the magnetic helicity becomes constant, which leads to a faster increase of the correlation length, as expected from magnetic helicity conservation. This also happens when the initial total chirality is imbalanced. Our findings have important implications for baryogenesis after axion inflation.

We build an analytical framework to study the observability of anisotropies and a net chiral polarization of the Stochastic Gravitational Wave Background (SGWB) with a generic network of ground-based detectors. We apply this formalism to perform a Fisher forecast of the performance of a network consisting of the current interferometers (LIGO, Virgo and KAGRA) and planned third-generation ones, such as the Einstein Telescope and Cosmic Explorer. Our results yield limits on the observability of anisotropic modes, spanning across noise- and signal-dominated regimes. We find that if the isotropic component of the SGWB has an amplitude close to the current limit, third-generation interferometers with an observation time of $10$ years can measure multipoles (in a spherical harmonic expansion) up to $\ell = 8$ with ${\cal O }\left( 10^{-3} - 10^{-2} \right)$ accuracy relative to the isotropic component, and an ${\cal O }\left( 10^{-3} \right)$ amount of net polarization. For weaker signals, the accuracy worsens as roughly the inverse of the SGWB amplitude.