Papers for Wednesday, Mar 29 2023

We investigate how feedback and environment shapes the X-ray scaling relations of early-type galaxies (ETGs), especially at the low-mass end. We select central-ETGs from the IllustrisTNG-100 box that have stellar masses $\log_{10}(M_{\ast}/\mathrm{M_{\odot}})\in[10.7, 11.9]$. We derive mock X-ray luminosity ($L_{\mathrm{X, 500}}$) and spectroscopic-like temperature ($T_{\mathrm{sl, 500}}$) of hot gas within $R_{500}$ of the ETG haloes using the MOCK-X pipeline. The scaling between $L_{\mathrm{X, 500}}$ and the total mass within 5 effective radii ($M_{5R_{\rm e}}$) agrees well with observed ETGs from Chandra. IllustrisTNG reproduces the observed increase in scatter of $L_{\mathrm{X, 500}}$ towards lower masses, and we find that ETGs with $\log_{10} (M_{5R_{\rm e}}/\mathrm{M_{\odot}}) \leqslant 11.5$ with above-average $L_{\mathrm{X, 500}}$ experienced systematically lower cumulative kinetic AGN feedback energy historically (vice versa for below-average ETGs). This leads to larger gas mass fractions and younger stellar populations with stronger stellar feedback heating, concertedly resulting in the above-average $L_{\mathrm{X, 500}}$. The $L_{\mathrm{X, 500}}$--$T_{\mathrm{sl, 500}}$ relation shows a similar slope to the observed ETGs but the simulation systematically underestimates the gas temperature. Three outliers that lie far below the $L_{\rm X}$--$T_{\rm sl}$ relation all interacted with larger galaxy clusters recently and demonstrate clear features of environmental heating. We propose that the distinct location of these backsplash ETGs in the $L_{\rm X}$--$T_{\rm sl}$ plane could provide a new way of identifying backsplash galaxies in future X-ray surveys.

Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here, we analyse a massive core in the well-studied infrared dark cloud (IRDC) called the "dragon cloud" (G028.37+00.07 or "Cloud C"). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the center of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, yet the more massive core remains devoid of any star formation indicators. Conclusions: We have presented evidence for the existence of a massive, prestellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular resolution observations to further assess its monolithic and prestellar nature.

When galaxies move through the intracluster medium (ICM) inside galaxy clusters, the ram pressure of the ICM can strip the gas from galaxies. The stripped gas forms tails on the trailing side. These galaxies are hence dubbed ``jellyfish galaxies''. ESO 137-001 is a quintessential jellyfish galaxy located in the nearest rich cluster, the Norma cluster. Its spectacular multiphase tail has complex morphology and kinematics both from the imprinted galaxy's interstellar medium (ISM) and as a result of the interactions between the stripped gas and the surrounding hot plasma, mediated by radiative cooling and magnetic fields. We study the kinematics of the multiphase tail using high-resolution observations of the ionized and the molecular gas in the entire structure. We calculate the velocity structure functions (VSFs) in moving frames along the tail and find that turbulence driven by Kelvin-Helmholtz (KH) instability quickly overwhelms the original ISM turbulence and saturates at $\sim 30$ kpc. There is also a hint that the far end of the tail has possibly started to inherit pre-existing large-scale ICM turbulence likely caused by structure formation. Turbulence measured by the molecular gas is generally consistent with that measured by the ionized gas in the tail but has a slightly lower amplitude. Most of the measured turbulence is below the mean free path of the hot ICM ($\sim 11$ kpc). Using warm/cool gas as a tracer of the hot ICM, we find that the isotropic viscosity of the hot plasma must be suppressed below 0.01 percent Spitzer level.

We report on the identification of a massive ($\sim10^5$ M$_\odot$) sub-structured stellar system in the Galactic Perseus complex likely undergoing hierarchical cluster assembly. Such a system comprises nine star clusters (including the well-known clusters NGC 654 and NGC 663) and an extended and low-density stellar halo. Gaia-DR3 and available spectroscopic data show that all its components are physically consistent in the 6D phase-space (position, parallax, and 3D motion), homogeneous in age (14 $-$ 44 Myr), and chemical content (half-solar metallicity). In addition, the system's global stellar density distribution is that of typical star clusters and shows clear evidence of mass segregation. We find that the hierarchical structure is mostly contracting towards the center with a speed of up to $\simeq4-5$ km s$^{-1}$, while the innermost regions expand at a lower rate (about $\simeq1$ km s$^{-1}$) and are dominated by random motions. Interestingly, this pattern is dominated by the kinematics of massive stars, while low-mass stars ($M<2$ M$_\odot$) are characterized by contraction across the entire cluster. Finally, the nine star clusters in the system are all characterized by a relatively flat velocity dispersion profile possibly resulting from ongoing interactions and tidal heating. We show that the observational results are generally consistent with those found in $N$-body simulations following the cluster violent relaxation phase strongly suggesting that the system is a massive cluster in the early assembly stages. This is the second structure with these properties identified in our Galaxy and, following the nomenclature of our previous work, we named it LISCA II.

We analyse all the available Chandra observations of the Cartwheel Galaxy and its compact group, taken between 2001 and 2008, with the main aim of addressing the variability in the X-ray band for this spectacular collisional ring galaxy. We focus on the study of point-like sources, in particular we are interested in Ultraluminous X-ray sources (ULXs, Lx >= 10^39 erg/s), that we treat as a class. We exploit archival XMM-Newton data to enrich the study of the long-term variability, on timescales of months to years. We find a total of 44 sources in the group area, of which 37 in total are ULXs positionally linked with the galaxies and of which we can study variability. They are 29 in the Cartwheel itself, 7 in G1 and 1 in G3. About one third of these 37 sources show long-term variability, while no variability is detected within the single observations. Of those, 5 ULXs have a transient behaviour with a maximum range of variability (Lmax/Lmin) of about one order of magnitude and are the best candidate neutron stars. The X-ray Luminosity Function (XLF) of the point-like sources remains consistent in shape between the Chandra observations both for the Cartwheel galaxy itself and for G1, suggesting that flux variability does not strongly influence the average properties of the population on the observation timescales.

Wind-fed models offer a unique way to form predictive models of the accretion flow surrounding Sagittarius A*. We present 3D, wind-fed MHD and GRMHD simulations spanning the entire dynamic range of accretion from parsec scales to the event horizon. We expand on previous work by including nonzero black hole spin and dynamically evolved electron thermodynamics. Initial conditions for these simulations are generated from simulations of the observed Wolf-Rayet stellar winds in the Galactic Centre. The resulting flow tends to be highly magnetized ($\beta \approx 2$) with an $\sim$ $r^{-1}$ density profile independent of the strength of magnetic fields in the winds. Our simulations reach the MAD state for some, but not all cases. In tilted flows, SANE jets tend to align with the angular momentum of the gas at large scales, even if that direction is perpendicular to the black hole spin axis. Conversely, MAD jets tend to align with the black hole spin axis. The gas angular momentum shows similar behavior: SANE flows tend to only partially align while MAD flows tend to fully align. With a limited number of dynamical free parameters, our models can produce accretion rates, 230 GHz flux, and unresolved linear polarization fractions roughly consistent with observations for several choices of electron heating fraction. Absent another source of large-scale magnetic field, winds with a higher degree of magnetization (e.g., where the magnetic pressure is 1/100 of the ram pressure in the winds) may be required to get a sufficiently large RM with consistent sign.

The triggering mechanism for the most luminous, quasar-like active galactic nuclei (AGN) remains a source of debate, with some studies favouring triggering via galaxy mergers, but others finding little evidence to support this mechanism. Here, we present deep Isaac Newton Telescope/Wide Field Camera imaging observations of a complete sample of 48 optically-selected type 2 quasars $-$ the QSOFEED sample (L$_{\rm [OIII]}>$10$^{8.5}$$L_{\odot}$; $z < 0.14$). Based on visual inspection by eight classifiers, we find clear evidence that galaxy interactions are the dominant triggering mechanism for quasar activity in the local universe, with 65$^{+6}_{-7}$ per cent of the type 2 quasar hosts showing morphological features consistent with galaxy mergers or encounters, compared with only 22$^{+5}_{-4}$ per cent of a stellar-mass- and redshift-matched comparison sample of non-AGN galaxies $-$ a 5$\sigma$ difference. The type 2 quasar hosts are a factor 3.0$^{+0.5}_{-0.8}$ more likely to be morphologically disturbed than their matched non-AGN counterparts, similar to our previous results for powerful 3CR radio AGN of comparable [OIII] emission-line luminosity and redshift. In contrast to the idea that quasars are triggered at the peaks of galaxy mergers as the two nuclei coalesce, and only become visible post-coalescence, the majority of morphologically-disturbed type 2 quasar sources in our sample are observed in the pre-coalescence phase (61$^{+8}_{-9}$ per cent). We argue that much of the apparent ambiguity that surrounds observational results in this field is a result of differences in the surface brightness depths of the observations, combined with the effects of cosmological surface brightness dimming.

We present GEO-FPT (Geometric Fitted Perturbation Theory), a new model for the galaxy bispectrum anisotropic signal in redshift space, with functional form rooted in perturbation theory. It also models the dependence of the bispectrum with the geometric properties of the triangles in Fourier space, and has a broader regime of validity than state-of-the-art theoretical models based on perturbation theory. We calibrate the free parameters of this model using high-resolution dark matter simulations and perform stringent tests to show that GEO-FPT describes the galaxy bispectrum accurately up to scales of $k\simeq0.12 h{\rm Mpc}^{-1}$ for different cosmological models, as well as for biased tracers of the dark matter field, considering a survey volume of $100$ (Gpc $h^{-1})^3$. In particular, a joint analysis of the power spectrum and bispectrum anisotropic signals, taking into account their full covariance matrix, reveals that the relevant physical quantities -- the BAO peak position (along and across the line-of-sight), and the growth of structure parameters times the amplitude of dark matter fluctuations, $f\sigma_8$-- are recovered in an unbiased way, with an accuracy better than $0.4\%$ and $2\%$ respectively (which is our $2\sigma$ statistical limit of the systematic error estimate). In addition, the bispectrum signal breaks the $f\sigma_8$ degeneracy without detectable bias: $f$ and $\sigma_8$ are recovered with better than $2.7\%$ and $3.8\%$ accuracy respectively (which is our $2\sigma$ statistical limit of the systematic error estimate). GEO-FPT boosts the applicability of the bispectrum signal of galaxy surveys beyond the current limitation of $k\lesssim 0.08\,h$ Mpc$^{-1}$ % and makes the bispectrum a key statistic to unlock the information content from the mildly non-linear regime in the on-going and forthcoming galaxy redshift surveys.

Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity $Z=10^{-11}$) with our population-synthesis code SEVN, and compare them with Population II binary stars ($Z=10^{-4}$). In our models, the lower edge of the pair-instability mass gap corresponds to a BH mass of $\approx{86}$ ($\approx{91}$) M$_\odot$ for single Population III (II) stars. Overall, we find only mild differences between the properties of binary BHs (BBHs) born from Population III and II stars, especially if we adopt the same initial mass function and initial orbital properties. Most BBH mergers born from Population III and II stars have primary BH mass below the pair-instability gap, and the maximum secondary BH mass is $ < 50$ M$_\odot$. Only up to $\approx{3.3}$% ($\approx{0.09}$%) BBH mergers from Population III (II) progenitors have primary mass above the gap. Unlike metal-rich binary stars, the main formation channel of BBH mergers from Population III and II stars involves only stable mass transfer episodes in our fiducial model.

The resonantly produced sterile neutrino ($N_1$) of the neutrino minimal standard model ($\nu$MSM) is a compelling dark matter candidate, especially in the reported possible detection of $N_1$ with mass $m_\mathrm{s}=7.1$~keV in X-ray decay. This particle will be accessible to the XRISM X-ray mission over the next 12 months. We revisit the physics behind $N_1$ and the uncertainty in its parameters. We compare predictions for the $m_\mathrm{s}=7.1$keV $N_1$ mixing angle, $\sin^{2}(2\theta)$, and half-mode mass, $M_\mathrm{hm}$, to existing X-ray observations and structure formation constraints. The strongest available constraints rule out $N_1$ as a dark matter candidate, and a more optimistic reading of the data prefers $\sin^{2}(2\theta)=5\times10^{-11}$ and $M_\mathrm{hm}=3.5\times10^{8}$$\mathrm{M}_{\odot}$. We highlight that the most promising upcoming opportunity for a detection is to find a line of velocity dispersion $\sim500$$\mathrm{kms^{-1}}$ in the Virgo cluster with XRISM, and then draw up a list of future objects of study to determine: (i) whether the line is from dark matter generally, and (ii) if from dark matter, whether that candidate is indeed $N_1$.

Supermassive black holes (SMBHs) are a key catalyst of galaxy formation and evolution, leading to an observed correlation between SMBH mass $M_{\rm BH}$ and host galaxy velocity dispersion $\sigma_{\rm e}$. Outside the local Universe, measurements of $M_{\rm BH}$ are usually only possible for SMBHs in an active state: limiting sample size and introducing selection biases. Gravitational lensing makes it possible to measure the mass of non-active SMBHs. We present models of the $z=0.169$ galaxy-scale strong lens Abell~1201. A cD galaxy in a galaxy cluster, it has sufficient `external shear' that a magnified image of a $z = 0.451$ background galaxy is projected just $\sim 1$ kpc from the galaxy centre. Using multi-band Hubble Space Telescope imaging and the lens modeling software $\texttt{PyAutoLens}$ we reconstruct the distribution of mass along this line of sight. Bayesian model comparison favours a point mass with $M_{\rm BH} = 3.27 \pm 2.12\times10^{10}\,$M$_{\rm \odot}$ (3$\sigma$ confidence limit); an ultramassive black hole. One model gives a comparable Bayesian evidence without a SMBH, however we argue this model is nonphysical given its base assumptions. This model still provides an upper limit of $M_{\rm BH} \leq 5.3 \times 10^{10}\,$M$_{\rm \odot}$, because a SMBH above this mass deforms the lensed image $\sim 1$ kpc from Abell 1201's centre. This builds on previous work using central images to place upper limits on $M_{\rm BH}$, but is the first to also place a lower limit and without a central image being observed. The success of this method suggests that surveys during the next decade could measure thousands more SMBH masses, and any redshift evolution of the $M_{\rm BH}$--$\sigma_{\rm e}$ relation. Results are available at https://github.com/Jammy2211/autolens_abell_1201.

Population III (Pop. III) binary stars likely produced the first stellar-born binary black hole (BBH) mergers in the Universe. Here, we quantify the main sources of uncertainty for the merger rate density evolution and mass spectrum of Pop. III BBHs by considering four different formation histories of Pop. III stars and 11 models of the initial orbital properties of their binary systems. The uncertainty on the orbital properties affects the BBH merger rate density by up to two orders of magnitude; models with shorter initial orbital periods lead to higher BBH merger rates, because they favour the merger via stable mass transfer episodes. The uncertainty on the star formation history also has a substantial impact on both the shape and the normalisation of the BBH merger rate density: the peak of the merger rate density shifts from $z\sim{8}$ up to $z\sim{16}$ depending on the assumed star formation rate, while the maximum BBH merger rate density for our fiducial binary population model spans from $\sim{2}$ to $\sim{30}$ Gpc$^{-3}$ yr$^{-1}$. The typical BBH masses are not affected by the star formation rate model and only mildly influenced by the binary population parameters. The primary black holes born from Pop. III stars tend to be rather massive ($30-40$ M$_\odot$) with respect to those born from metal-rich stars ($8-10$ M$_\odot$). However, we expect that Pop. III BBH mergers with primary mass $m_1>60$ M$_\odot$ are rare ($<10^{-2}$ Gpc$^{-3}$ yr$^{-1}$). Finally, we estimate that the Einstein Telescope will detect $10-10^4$ Pop. III BBH mergers per year, depending on the star formation history and binary star properties.

High-energy processes are ubiquitous even in the earliest stages of protostellar evolution. Motivated by the results of our systematic search for intense centimeter radio flares in Young Stellar Objects (YSOs) and by rare findings of strong millimeter-wavelength variability, we have conducted a systematic search for such variability in the Orion Nebula Cluster (ONC) using Atacama Large Millimeter/submillimeter Array (ALMA). Rapid variability on timescales of minutes to hours in the (centimeter)millimeter-wavelength range indicates (gyro)synchrotron radiation. Additionally, mass accretion will also affect the millimeter-wavelength luminosity but typically on longer timescales. Beyond studies of individual YSOs, our characterization of strong millimeter-wavelength variability with ALMA in the ONC sets first systematic constraints on the occurrence of such variability in a large number of YSOs ($\sim$130). We report the discovery of an order of magnitude millimeter-flare within just a few minutes from a known YSO previously reported as a radio flaring source at cm-wavelengths (the "ORBS'' source). We also present an assessment of the systematic variability effects caused by the use of time-sliced imaging of a complex region. These are mostly due to the impact of a changing synthesized beam throughout the observations. We use simulated ALMA observations to reproduce and quantify these effects and set a lower limit for the variability that can be studied using our method in a complex region such as the ONC. Our results demonstrate that the utility of time domain analysis of YSOs extends into the millimeter-wavelength range, potentially interfering with the conversion of observed fluxes into dust masses.

Gas inflow into galaxies should affect the star formation and hence the evolution of galaxies across cosmic time. In this work, we use TNG100 of the IllustrisTNG simulations to understand the role of environment on gas inflow rates in massive galaxies at z >= 2. We divide our galaxies (log(M*/Msolar )>= 10.5) into cluster (log Mhalo/Msolar >= 13) and field (log Mhalo/Msolar < 13) galaxies at z = 2 and further divide into centrals and satellites. We track their gas inflow rates from z = 6 to 2 and find that the total gas inflow rates of satellite galaxies rapidly decline after their infall into cluster halos and as they reach the cluster center. At z = 2, the gas inflow rate of cluster satellite galaxies is correlated with the cluster-centric radii and not the host halo mass. In contrast, the gas inflow rate in centrals is strongly correlated with the host halo mass at z >= 2. Our study indicates that between redshifts 6 to 2, the gas that normally is accreted by the satellite galaxies is redirected to the center of the cluster halo as inflows to the cluster centrals and forming the intra-cluster medium. Our analysis suggest that the inequality of gas accretion between massive satellite and central galaxies is responsible for the starvation of cluster satellite galaxies that evolve into the massive quenched cluster galaxies observed at z<0.5.

Low-luminosity active galactic nuclei are strong sources of X-ray emission produced by Compton scattering originating from the accretion flows surrounding their supermassive black holes. The shape and energy of the resulting spectrum depend on the shape of the underlying electron distribution function (DF). In this work, we present an extended version of the grmonty code, called $\kappa$monty. The grmonty code previously only included a thermal Maxwell J\"utner electron distribution function. We extend the gromty code with non-thermal electron DFs, namely the $\kappa$ and power-law DFs, implement Cartesian Kerr-Schild coordinates, accelerate the code with MPI, and couple the code to the non-uniform AMR grid data from the GRMHD code BHAC. For the Compton scattering process, we derive two sampling kernels for both distribution functions. Finally, we present a series of code tests to verify the accuracy of our schemes. The implementation of non-thermal DFs opens the possibility of studying the effect of non-thermal emission on previously developed black hole accretion models.

The Euclid mission will conduct an extragalactic survey over 15000 deg$^2$ of the extragalactic sky. The spectroscopic channel of the Near-Infrared Spectrometer and Photometer (NISP) has a resolution of $R\sim450$ for its blue and red grisms that collectively cover the $0.93$--$1.89 $\micron;range. NISP will obtain spectroscopic redshifts for $3\times10^7$ galaxies for the experiments on galaxy clustering, baryonic acoustic oscillations, and redshift space distortion. The wavelength calibration must be accurate within $5$\AA to avoid systematics in the redshifts and downstream cosmological parameters. The NISP pre-flight dispersion laws for the grisms were obtained on the ground using a Fabry-Perot etalon. Launch vibrations, zero gravity conditions, and thermal stabilisation may alter these dispersion laws, requiring an in-flight recalibration. To this end, we use the emission lines in the spectra of compact planetary nebulae (PNe), which were selected from a PN data base. To ensure completeness of the PN sample, we developed a novel technique to identify compact and strong line emitters in Gaia spectroscopic data using the Gaia spectra shape coefficients. We obtained VLT/X-SHOOTER spectra from $0.3$ to $2.5$ \micron;for 19 PNe in excellent seeing conditions and a wide slit, mimicking Euclid's slitless spectroscopy mode but with 10 times higher spectral resolution. Additional observations of one northern PN were obtained in the $0.80$--$1.90$ \micron range with the GMOS and GNIRS instruments at the Gemini North observatory. The collected spectra were combined into an atlas of heliocentric vacuum wavelengths with a joint statistical and systematic accuracy of 0.1 \AA in the optical and 0.3 \AA in the near-infrared. The wavelength atlas and the related 1D and 2D spectra are made publicly available.

A major puzzle concerning the wide stellar binaries (semimajor axes $a\gtrsim 10^3$ AU) in the Solar neighborhood is the origin of their observed superthermal eccentricity distribution function (DF), which is well-approximated by $P(e) \propto e^\alpha$ with $\alpha \approx 1.3$. This DF evolves under the combined influence of (i) tidal torques from the Galactic disk and (ii) scattering by passing stars, molecular clouds, and substructure. Recently, Hamilton (2022) (H22) demonstrated that Galactic tides alone cannot produce a superthermal eccentricity DF from an initially isotropic, non-superthermal one, under the restrictive assumptions that the eccentricity DF was initially of power law form and then was rapidly phase-mixed toward a steady state by the tidal perturbation. In this paper we first prove analytically that H22's conclusions are in fact valid at all times, regardless of these assumptions. We then adopt H22's Galactic disk model and numerically integrate the equations of motion for several ensembles of isotropically oriented wide binaries to study the time evolution in detail. We find that even non-power law DFs can be described by an effective power law index $\alpha_\mathrm{eff}$ which accurately characterizes both their initial and final states. Any DF with initial (effective or exact) power law index $\alpha_\mathrm{i}$ is transformed by Galactic tides into another power law with index $\alpha_\mathrm{f} \approx (1+\alpha_\mathrm{i})/2$ on a timescale $\sim 2$ Gyr $(a/10^4\mathrm{AU})^{-3/2}$. In a companion paper, we investigate separately the effect of stellar scattering. As the GAIA data continues to improve, these results will place strong constraints on wide binary formation channels.

We report CFHT/SITELLE imaging Fourier Transform Spectrograph observations of the Brightest Cluster Galaxy (BCG) of galaxy cluster Abell 2390 at z=0.228. The BCG displays a prominent cone of emission in H$\alpha$, H$\beta$, [NII], and [OII] to the North-West with PA = 42$^o$, 4.4 arcsec in length (15.9 kpc), which is associated with elongated and asymmetric Chandra soft X-ray emission. The H$\alpha$ flux map also contains a "hook" of H$\alpha$ and [NII] emission resulting in a broadened northern edge to the cone. Using SITELLE/LUCI software we extract emission line flux, velocity, velocity dispersion, and continuum maps, and utilize them to derive flux ratio maps to determine ionization mechanisms and dynamical information in the BCG's emission line region. The Baldwin-Phillips-Terlevich diagnostics on the BCG cone indicate a composite ionization origin of photoionization due to star formation and shock. Strong LINER-like emission is seen in the nuclear region which hosts an AGN. As Abell 2390 is a cool-core cluster, we suggest that the cooling flow is falling onto the central BCG and interacting with the central AGN. The AGN produces jets that inflate "bubbles" of plasma in the ICM, as is often observed in local galaxy clusters. Furthermore, combining signs of AGN activities from radio, optical emission line and X-ray data over a large range of physical scale, we find evidence for three possible episodes of AGN activity in different epochs associated with the Abell 2390 BCG.

We test the regime of validity of the effective field theory (EFT) of intrinsic alignments (IA) at the one-loop level by comparing with 3D halo shape statistics in N-body simulations. This model is based on the effective field theory of large-scale structure (EFT of LSS) and thus a theoretically well-motivated extension of the familiar non-linear alignment (NLA) model and the tidal-alignment-tidal-torquing (TATT) model. It contains a total of $8$ free bias parameters. Specifically, we measure the dark matter halo shape-shape multipoles $P_{EE}^{(0)}(k), P_{EE}^{(2)}(k), P_{BB}^{(0)}(k), P_{BB}^{(2)}(k)$ as well as the matter-shape multipoles $P_{\delta E}^{(0)}(k), P_{\delta E}^{(2)}(k)$ from the simulations and perform a joint fit to determine the largest wavenumber $k_{\text{max}}$ up to which the theory predictions from the EFT of IA are consistent with the measurements. We find that the EFT of IA is able to describe intrinsic alignments of dark matter halos up to $k_\text{max}=0.30\,h/\text{Mpc}$ at $z=0$. This demonstrates a clear improvement over other existing alignment models like NLA and TATT, which are only accurate up to $k_\text{max}=0.05\,h/\text{Mpc}$ . We examine the posterior distributions of the higher-order bias parameters, and show that their inclusion is necessary to describe intrinsic alignments in the quasi-linear regime. Further, the EFT of IA is able to accurately describe the auto-spectrum of intrinsic alignment B-modes, in contrast to the other alignment models considered.

The Positron Puzzle is a half-century old conundrum about the origin of the Galactic $\gamma$-ray emission line at photon energies of 511 keV, and the shape of its morphology, showing a bulge-to-disk luminosity ratio of $\sim 1$ - unlike any astrophysical source distribution. Positrons that have been cooled to the eV scale capture electrons and form the intermediate bound state of Positronium (Ps) which decays on a nano-second timescale into two or three photons. Assuming the emission to originate from the Galactic bulge, centre, and disk, a visible annihilation rate in the Milky Way of $\sim 5 \times 10^{43}\,\mathrm{e^+\,s^{-1}}$ has to be explained, either by a quasi-steady state of production and annihilation, or by possibly multiple burst-like events that flood the Galaxy with positrons, then fading away on a Myr timescale. In this paper, I will review what the real Positron Puzzle is, where data and simulations have been used inadequately which resulted in false claims and an apparent quandary, what we really know and absolutely not know about the topic, and how this epistemic problem might be advancing.

With long data sets available for asteroseismology from space missions, it is sometimes necessary to deal with time series that have large gaps. This is becoming particularly relevant for TESS, which is revisiting many fields on the sky every two years. Because solar-like oscillators have finite mode lifetimes, it has become tempting to close large gaps by shifting time stamps. Using actual data from the Kepler Mission, we show that this results in artificial structures in the power spectrum that compromise the measurements of mode frequencies and linewidths.

PKS 2005-489 is a well-known, bright southern BL Lac object that has been detected up to TeV energies. In a low-flux state it exhibits the expected multiwavelength double-peaked spectrum in the radio -- $\gamma$-ray band. The high-flux state shows extreme flux variations in the X-ray band with a hardening as well as a peculiar curved feature in the spectrum. Thus far, PKS 2005-489 is the only source to exhibit such a feature. To study the X-ray variability further, we obtained the first hard X-ray spectrum of the source with NuSTAR. We compare quasi-simultaneous radio, optical, UV, soft and hard X-ray, and $\gamma$-ray data of PKS 2005-489 to archival data in order to study its broadband behavior. We find a very consistent quiet state in the SED, with little variation in spectral shape or flux between the 2012 and 2020 data. A possible explanation for the peculiar X-ray spectrum in the flaring state is an additional component in the jet, possibly accelerated via magnetic reconnection, that is not co-spatial to the low-flux state emission region.

Secular gravitational instability (GI) is one promising mechanism for explaining planetesimal formation. The previous studies on secular GI utilized a razor-thin disk model and derived the growth condition in terms of the vertically integrated physical values such as dust-to-gas surface density ratio. However, in weakly turbulent disks where secular GI can operate, a dust disk can be orders of magnitude thinner than a gas disk, and analyses treating the vertical structures are necessary to clarify the interplay of the midplane dust motion and the upper gas motion. In this work, we perform vertically global linear analyses of secular GI with the vertical domain size of a few gas scale heights. We find that dust grains accumulate radially around the midplane while gas circulates over the whole vertical region. We obtain well-converged growth rates when the outer gas boundary is above two gas scale heights. The growth rates are underestimated if we assume the upper gas to be steady and regard it just as the source of external pressure to the dusty lower layer. Therefore, treating the upper gas motion is important even when the dust disk is much thinner than the gas disk. Conducting a parameter survey, we represent the growth condition in terms of the Toomre's $Q$ value for dust and dust-to-gas surface density ratio. The critical dust disk mass for secular GI is $\sim10^{-4}$ stellar mass for the dust-to-gas surface density ratio of 0.01, the Stokes number of 0.1, and dimensionless radial dust diffusivity of $10^{-4}$.

The Ha line, one of the most studied chromospheric diagnostics, is a tracer for magnetic field structures, while its line core intensity provides an estimate of the mass density. The brightness temperatures from Atacama Large Millimetre-submm Array (ALMA) observations provide a complementary view of the activity and the thermal structure of stellar atmospheres. These two diagnostics together can provide insights into the physical properties of stellar atmospheres. In this paper, we present a comparative study between synthetic continuum brightness temperature maps for mm wavelengths (0.3 mm to 8.5 mm) and the width of the Ha 6565{\AA} line. The 3D radiative transfer codes Multi3D and Advanced Radiative Transfer (ART) are used to calculate synthetic spectra for the Ha line and the mm continua respectively, from an enhanced network atmosphere model with non-equilibrium hydrogen ionisation generated with the state-of-the-art 3D rMHD code Bifrost. We use Gaussian Point Spread Function (PSF) for simulating the effect of ALMA's limited spatial resolution and calculate the Ha vs. mm continuum correlations and slopes of scatter plots for the original and degraded resolution of the whole box, quiet sun and enhanced network patches separately. The Ha linewidth and mm brightness temperatures are highly correlated and the correlation is highest at a wavelength 0.8 mm i.e. ALMA Band 7. The correlation increases with decreased resolution. On the other hand, the slopes decrease with increasing wavelength. The degradation of resolution does not have a significant impact on the calculated slopes. With decreasing spatial resolution the standard deviations of the observables, Ha linewidth and brightness temperatures decrease and the correlations between them increase, but the slopes do not change significantly. Hence, these relations may prove useful to calibrate the mm continuum maps observed with ALMA.

Measuring broad emission-line widths in active galactic nuclei (AGN) is not straightforward owing to the complex nature of flux variability in these systems. Line-width measurements become especially challenging when signal-to-noise is low, profiles are narrower, or spectral resolution is low. We conducted an extensive correlation analysis between emission-line measurements from the optical spectra of Markarian 142 (Mrk 142; a narrow-line Seyfert galaxy) taken with the Gemini North Telescope (Gemini) at a spectral resolution of 185.6+\-10.2 km/s and the Lijiang Telescope (LJT) at 695.2+\-3.9 km/s to investigate the disparities in the measured broad-line widths from both telescope data. Mrk~142 posed a challenge due to its narrow broad-line profiles, which were severely affected by instrumental broadening in the lower-resolution LJT spectra. We discovered that allowing the narrow-line flux of permitted lines having broad and narrow components to vary during spectral fitting caused a leak in the narrow-line flux to the broad component, resulting in broader broad-line widths in the LJT spectra. Fixing the narrow-line flux ratios constrained the flux leak and yielded the Hydrogen-beta broad-line widths from LJT spectra $\sim$54\% closer to the Gemini Hydrogen-beta widths than with flexible narrow-line ratios. The availability of spectra at different resolutions presented this unique opportunity to inspect how spectral resolution affected emission-line profiles in our data and adopt a unique method to accurately measure broad-line widths. Reconsidering line-measurement methods while studying diverse AGN populations is critical for the success of future reverberation-mapping studies. Based on the technique used in this work, we offer recommendations for measuring line widths in narrow-line AGN.

Tracking the motions of transient jets launched by low-mass X-ray binaries (LMXBs) is critical for determining the moment of jet ejection, and identifying any corresponding signatures in the accretion flow. However, these jets are often highly variable and can travel across the resolution element of an image within a single observation, violating a fundamental assumption of aperture synthesis. We present a novel approach in which we directly fit a single time-dependent model to the full set of interferometer visibilities, where we explicitly parameterise the motion and flux density variability of the emission components, to minimise the number of free parameters in the fit, while leveraging information from the full observation. This technique allows us to detect and characterize faint, fast-moving sources, for which the standard time binning technique is inadequate. We validate our technique with synthetic observations, before applying it to three Very Long Baseline Array (VLBA) observations of the black hole candidate LMXB MAXI J1803-298 during its 2021 outburst. We measured the proper motion of a discrete jet component to be $1.37\pm0.14$ mas/hr, and thus we infer an ejection date of MJD $59348.08_{-0.06}^{+0.05}$, which occurs just after the peak of a radio flare observed by the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/Sub-Millimeter Array (ALMA), while MAXI J1803-298 was in the intermediate state. Further development of these new VLBI analysis techniques will lead to more precise measurements of jet ejection dates, which, combined with dense, simultaneous multi-wavelength monitoring, will allow for clearer identification of jet ejection signatures in the accretion flow.

The reconstruction of cosmic-ray-induced extensive air showers with a non-imaging Cherenkov detector array requires knowledge of the Cherenkov yield of any given air shower for a given set of shower parameters. Although air showers develop in a stochastic cascade, certain characteristics of the particles in the shower have been shown to come from universal probability distributions, a property known as shower universality. Both the energy and the angular distributions of charged particles within a shower have been parameterized. One can use these distributions to calculate the Cherenkov photon yield as an angular distribution from the Cherenkov cones of charged particles at various stages of shower development. This Cherenkov photon yield can then be tabulated for use in the reconstruction of air showers. In this work, we develop the calculation of both the Cherenkov angular distribution and Cherenkov yield per shower particle, and show how a look-up table was constructed to capture the relevant features of these distributions for general use. We compare the results of our calculations with the results of full, particle-stack, Monte Carlo simulation of the Cherenkov light produced in extensive air showers using CORSIKA-IACT. We make comparisons of both the lateral distribution of the Cherenkov photon flux amongst several detectors and of the arrival-time distribution of the Cherenkov photons in a single detector.

Because of their period-luminosity relation (PLR), classical Cepheids play a key role in the calibration of the extragalactic distance scale and the determination of the Hubble-Lema\^{i}tre constant $H_0$. Recent findings show that the majority of classical Cepheids should be in binary or multiple systems, which might undermine their accuracy, as the extra -- and unaccounted for -- light from the companions of Cepheids causes a shift in the PLR. We quantify this shift using synthetic populations of binary Cepheids that we developed for this purpose, as described in Paper I of this series. We find that while all PLRs are shifted toward brighter values due to the excess light from the companions, the bias in the relative distance modulus between two galaxies hosting binary Cepheids can be either positive or negative, depending on the percentage of binary Cepheids in them. If the binarity percentage in the two galaxies is similar, the effect of binarity is canceled. Otherwise, it introduces a shift in the distance modulus of the order of millimags in the near-infrared passbands and Wesenheit indices, and tens of millimags in the visual domain; its exact value depends on the variant of the synthetic population (a unique combination of metallicity, star formation history, shape and location of the instability strip, and initial parameter distributions). Such shifts in distance moduli to type Ia supernova host galaxies introduce an additional statistical error on $H_0$, which however does not prevent measuring $H_0$ with a precision of 1%.

[Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the steady-state MRI-driven accretion, by employing and improving our previous 1+1D MRI-driven turbulence model. Second, we relax the steady-state accretion assumption in this disk accretion model, and partially couple it to a dust evolution model in order to investigate how the evolution of dust (dynamics and grain growth processes combined) and MRI-driven accretion are intertwined on million-year timescales. Results. Dust coagulation and settling lead to a higher gas ionization degree in the protoplanetary disk, resulting in stronger MRI-driven turbulence as well as a more compact dead zone. On the other hand, fragmentation has an opposite effect because it replenishes the disk in small dust particles. Since the dust content of the disk decreases over million years of evolution due to radial drift, the MRI-driven turbulence overall becomes stronger and the dead zone more compact until the disk dust-gas mixture eventually behaves as a grain-free plasma. Furthermore, our results show that dust evolution alone does not lead to a complete reactivation of the dead zone. Conclusions. The MRI activity evolution (hence the temporal evolution of the MRI-induced $\alpha$-parameter) is controlled by dust evolution and occurs on a timescale of local dust growth, as long as there is enough dust particles in the disk to dominate the recombination process for the ionization chemistry. Once it is no longer the case, it is expected to be controlled by gas evolution and occurs on a viscous evolution timescale.

The observation of short gamma ray bursts (SGRBs) in the TeV energy range plays an important role in understanding the radiation mechanism and probing new areas of physics such as Lorentz invariance violation. However, no SGRB has been observed in this energy range due to the short duration of SGRBs and the weakness of current experiments. New experiments with new technology are required to detect sub-TeV SGRBs. In this work, we observe the very high energy (VHE) $\gamma$-ray emissions from SGRBs and calculate the annual detection rate with the High Altitude Detection of Astronomical Radiation HADAR (HADAR) experiment. First, a set of pseudo-SGRB samples is generated and checked using the observations of Fermi-GBM, Fermi-LAT, and SWIFT measurements. The annual detection rate is calculated from these SGRB samples based on the performance of the HADAR instrument. As a result, the HADAR experiment can detect 0.5 SGRB per year if the spectral break-off of $\gamma$-rays caused by the internal absorption is larger than 100 GeV. For a GRB09010-like GRB in HADAR's view, it should be possible to detect approximately 2000 photons considering the internal absorption. With a time delay assumption due to the Lorentz invariance violation effects, a simulated light curve of GRB090510 has evident energy dependence. We hope that the HADAR experiment can perform the SGRB observations and test our calculations in the future.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) has acquired tens of millions of low-resolution spectra of stars. This paper investigated the parameter estimation problem for these spectra. To this end, we proposed a deep learning model StarGRU network (StarGRUNet). This network was further applied to estimate the stellar atmospheric physical parameters and 13 elemental abundances from LAMOST low-resolution spectra. On the spectra with signal-to-noise ratios greater than or equal to $5$, the estimation precisions are $94$ K and $0.16$ dex on $T_\texttt{eff}$ and $\log \ g$ respectively, $0.07$ dex to $0.10$ dex on [C/H], [Mg/H], [Al/H], [Si/H], [Ca/H], [Ni/H] and [Fe/H], and $0.10$ dex to $0.16$ dex on [O/H], [S/H], [K/H], [Ti/H] and [Mn/H], and $0.18$ dex and $0.22$ dex on [N/H] and [Cr/H] respectively. The model shows advantages over available models and high consistency with high-resolution surveys. We released the estimated catalog computed from about 8.21 million low-resolution spectra in LAMOST DR8, code, trained model, and experimental data for astronomical science exploration and data processing algorithm research respectively.

Cylindrical molecular filaments are observed to be the main sites of Sun-like star formation, while massive stars form in dense hubs, at the junction of multiple filaments. The role of hub-filament configurations has not been discussed yet in relation to the birth environment of the solar system and to infer the origin of isotopic ratios of Short-Lived Radionuclides (SLR, such as $^{26}$Al) of Calcium-Aluminum-rich Inclusions (CAIs) observed in meteorites. In this work, we present simple analytical estimates of the impact of stellar feedback on the young solar system forming along a filament of a hub-filament system. We find that the host filament can shield the young solar system from the stellar feedback, both during the formation and evolution of stars (stellar outflow, wind, and radiation) and at the end of their life (supernovae). We show that the young solar system formed along a dense filament can be enriched with supernova ejecta (e.g., $^{26}$Al) during the formation timescale of CAIs. We also propose that the streamers recently observed around protostars may be channeling the SLR-rich material onto the young solar system. We conclude that considering hub-filament configurations as the birth environment of the Sun is important when deriving theoretical models explaining the observed properties of the solar system.

Novae were discovered to emit transient gamma rays during the period of several days to a few weeks after initial explosion, indicating presence of acceleration processes of particles in their expanding shells. In the case of recurrent novae, electrons can be in principle accelerated in the nova shells for the whole recurrence period of nova producing delayed $\gamma$ ray emission as considered in Bednarek (2022). Here we extend the ideas presented in this article by considering the fate of electrons which diffuse out of the shells of novae supplying fresh relativistic electrons to the recurrent nova super-remnants during the whole active period of nova ($\ge 10^4$ yrs). We develop a model for the acceleration of electrons and their escape from the nova shells. The electrons within the recurrent nova super-remnants produce $\gamma$ rays in the comptonization process of the radiation from the red giant companion and the Cosmic Microwave Background Radiation. As an example, the case of a symbiotic nova RS Oph (with the recurrence period estimated on $\sim$10-50 yrs) is considered in more detail. Predicted $\gamma$-ray emission from the nova super-remnant around RS Oph is discussed in the context of its observability by satellite experiments (i.e. Fermi-LAT) as well as current and future Cherenkov telescopes.

How structures, e.g., magnetic loops, in the upper atmosphere, i.e., the transition region and corona, are heated and sustained is one of the major unresolved issues in solar and stellar physics. Various theoretical and observational studies on the heating of coronal loops have been undertaken. The heating of quiescent loops caused by eruptions is, however, rarely observed. In this study, employing data from the Solar Dynamics Observatory (SDO) and Solar Upper Transition Region Imager (SUTRI), we report the heating of quiescent loops associated with nearby eruptions. In active regions (ARs) 13092 and 13093, a long filament and a short filament, and their overlying loops are observed on 2022 September 4. In AR 13093, a warm channel erupted toward the northeast, whose material moved along its axis toward the northwest under the long filament, turned to the west above the long filament, and divided into two branches falling to the solar surface. Subsequently, the short filament erupted toward the southeast. Associated with these two eruptions, the quiescent loops overlying the long filament appeared in SDO/Atmospheric Imaging Assembly (AIA) high-temperature images, indicating the heating of loops. During the heating, signature of magnetic reconnection between loops is identified, including the inflowing motions of loops, and the formation of X-type structures and newly reconnected loops. The heated loops then cooled down. They appeared sequentially in AIA and SUTRI lower-temperature images. All the results suggest that the quiescent loops are heated by reconnection between loops caused by the nearby warm channel and filament eruptions.

[abridged] We explore the possibility that the dark matter (DM) component in galaxies may originate fractional gravity. In such a framework, the standard law of inertia continues to hold, but the gravitational potential associated to a given DM density distribution is determined by a modified Poisson equation including fractional derivatives (i.e., derivatives of non-integer type), that are meant to describe non-local effects. We derive analytically the expression of the potential that in fractional gravity corresponds to various spherically symmetric density profiles, including the Navarro-Frenk-White (NFW) distribution that is usually exploited to describe virialized halos of collisionless DM as extracted from $N-$body cosmological simulations. We show that in fractional gravity the dynamics of a test particle moving in a cuspy NFW density distribution is substantially altered with respect to the Newtonian case (i.e., basing on the standard Poisson equation), mirroring what in Newtonian gravity would instead be sourced by a density profile with an inner core. We test the fractional gravity framework on galactic scales, showing that: (i) it can provide accurate fits to the stacked rotation curves of galaxies with different properties; (ii) it can reproduce to reasonable accuracy the observed shape and scatter of the radial acceleration relation (RAR); (iii) it can properly account for the universal surface density and the core radius vs. disk scale-length scaling relations. Finally, we discuss the possible origin of the fractional gravity behavior as a fundamental or emerging property of the elusive DM component.

The presence of carbon-chain molecules in the interstellar medium (ISM) has been known since the early 1970s and $>100$ such species have been identified to date, making up $>40\%$ of the total of detected ISM molecules. They are prevalent not only in star-forming regions in our Galaxy, but also in other galaxies. These molecules provide important information on physical conditions, gas dynamics, and evolutionary stages of star-forming regions. More complex species of polycyclic aromatic hydrocarbons (PAHs) and fullerenes (C$_{60}$ and C$_{70}$) have been detected in circumstellar envelopes around carbon-rich Asymptotic Giant Branch (AGB) stars and planetary nebulae, while PAHs are also known to be a widespread component of interstellar dust in most galaxies. Recently, two line survey projects toward the starless core Taurus Molecular Cloud-1 with large single-dish telescopes have detected many new carbon-chain species, including molecules containing benzene rings. These new findings raise fresh questions about carbon-bearing species in the Universe. This article reviews various aspects of carbon-chain molecules, including observational studies, chemical simulations, quantum calculations, and laboratory experiments, and discusses open questions and how they may be answered by future facilities.

An array of large observational programs using ground-based and space-borne telescopes is planned in the next decade. The forthcoming wide-field sky surveys are expected to deliver a sheer volume of data exceeding an exabyte. Processing the large amount of multiplex astronomical data is technically challenging, and fully automated technologies based on machine learning and artificial intelligence are urgently needed. Maximizing scientific returns from the big data requires community-wide efforts. We summarize recent progress in machine learning applications in observational cosmology. We also address crucial issues in high-performance computing that are needed for the data processing and statistical analysis.

The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic winds in driving radial metallicity profiles and setting the mass-metallicity gradient relation (MZGR). We include a comprehensive treatment of these processes by including them as sources that supply mass, metals, and energy to marginally unstable galactic discs in pressure and energy balance. We show that both feedback and accretion that can drive turbulence and enhance metal-mixing via diffusion are crucial to reproduce the observed MZGR in local galaxies. Metal transport also contributes to setting metallicity profiles, but it is sensitive to the strength of radial gas flows in galaxies. While the mass loading of galactic winds is important to reproduce the mass metallicity relation (MZR), we find that metal mass loading is more important to reproducing the MZGR. Specifically, our model predicts preferential metal enrichment of galactic winds in low-mass galaxies. This conclusion is robust against our adopted scaling of the wind mass-loading factor, uncertainties in measured wind metallicities, and systematics due to metallicity calibrations. Overall, we find that at $z \sim 0$, galactic winds and metal transport are more important in setting metallicity gradients in low-mass galaxies whereas star formation feedback and gas accretion dominate setting metallicity gradients in massive galaxies.

Gamma-ray flares from Active Galactic Nuclei (AGN) show substantial variability on ultrafast timescales (i.e. shorter than the light crossing time of the AGN's supermassive black hole). We propose that ultrafast variability is a byproduct of the turbulent dissipation of the jet Poynting flux. Due to the intermittency of the turbulent cascade, the dissipation is concentrated in a set of reconnecting current sheets. Electrons energised by reconnection have a strong pitch angle anisotropy, i.e. their velocity is nearly aligned with the guide magnetic field. Then each current sheet produces a narrow radiation beam, which dominates the emission from the whole jet when it is directed towards the observer. The ultrafast variability is set by the light crossing time of a single current sheet, which is much shorter than the light crossing time of the whole emission region. The predictions of our model are: (i) The bolometric luminosity of ultrafast AGN flares is dominated by the inverse Compton (IC) emission, as the lower energy synchrotron emission is suppressed due to the pitch angle anisotropy. (ii) If the observed luminosity includes a non-flaring component, the variations of the synchrotron luminosity have a small amplitude. (iii) The synchrotron and IC emission are less variable at lower frequencies, as the cooling time of the radiating particles exceeds the light crossing time of the current sheet. Simultaneous multiwavelength observations of ultrafast AGN flares can test these predictions.

Analytical models of magnetised, geometrically thick discs are relevant to understand the physical conditions of plasma around compact objects and to explore its emitting properties. This has become increasingly important in recent years in the light of the Event Horizon Telescope observations of Sgr A* and M87. Models of thick discs around black holes usually consider constant angular momentum distributions and do not take into account the magnetic response of the fluid to applied magnetic fields. We present a generalisation of our previous work on stationary models of magnetised accretion discs with magnetic polarisation (Pimentel et al. 2018). This extension is achieved by accounting for non-constant specific angular momentum profiles, done through a two-parameter ansatz for those distributions. We build a large number of new equilibrium solutions of thick discs with magnetic polarisation around Kerr black holes, selecting suitable parameter values within the intrinsically substantial parameter space of the models. We study the morphology and the physical properties of those solutions, finding qualitative changes with respect to the constant angular momentum tori of (Pimentel et al. 2018). However, the dependences found on the angular momentum distribution or on the black hole spin do not seem to be strong. Some of the new solutions, however, exhibit a local maximum of the magnetisation function, absent in standard magnetised tori. Due to the enhanced development of the magneto-rotational instability as a result of magnetic susceptibility, those models might be particularly well-suited to investigate jet formation through general-relativistic MHD simulations. The new equilibrium solutions reported here can be used as initial data in numerical codes to assess the impact of magnetic susceptibility in the dynamics and observational properties of thick disc-black hole systems.

The bright blazar OJ~287 routinely parades high brightness bremsstrahlung flares, which are explained as being a result of a secondary supermassive black hole (SMBH) impacting the accretion disc of a more massive primary SMBH in a binary system. The accretion disc is not rigid but rather bends in a calculable way due to the tidal influence of the secondary. Below we refer to this phenomenon as a variable disc level. We begin by showing that these flares occur at times predicted by a simple analytical formula, based on general relativity inspired modified Kepler equation, which explains impact flares since 1888. The 2022 impact flare, namely flare number 26, is rather peculiar as it breaks the typical pattern of two impact flares per 12-year cycle. This is the third bremsstrahlung flare of the current cycle that follows the already observed 2015 and 2019 impact flares from OJ~287. It turns out that the arrival epoch of flare number 26 is sensitive to the level of primary SMBH's accretion disc relative to its mean level in our model. We incorporate these tidally induced changes in the level of the accretion disc to infer that the thermal flare should have occurred during July-August 2022, when it was not possible to observe it from the Earth. Thereafter, we explore possible observational evidence for certain pre-flare activity by employing spectral and polarimetric data from our campaigns in 2004/05 and 2021/22. We point out theoretical and observational implications of two observed mini-flares during January-February 2022.

We take advantage of the unprecedented dynamical range provided by the "Cosmic-Zoom" project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the $\Lambda$CDM paradigm when the dark matter is made of weakly interacting massive particles of mass $100\ \mathrm{GeV}$. In particular, we complement previous studies by exploring the MAHs of haloes with mass from $10^8\ h^{-1}\mathrm{M_{\odot}}$ down to Earth mass, $10^{-6}\ h^{-1}\mathrm{M_{\odot}}$. The formation redshift of low-mass haloes anti-correlates weakly with mass, peaking at $z=3$ for haloes of mass $10^{-4}\ h^{-1}\mathrm{M_{\odot}}$. Even lower masses are affected by the free-streaming cutoff in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press-Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological $N$-body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, $M>10^{10}\ h^{-1}\mathrm{M_{\odot}}$, the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than $20$ percent at high redshift. We conclude that EPS theory predicts the hierarchical build-up of dark matter haloes quite well over the entire mass range.

We have implemented ad-hoc modifications to the CORSIKA Monte-Carlo generator which allow us to simultaneously adjust the multiplicity, elasticity and cross-section of hadronic interactions with respect to the predictions of the Sibyll 2.3d interaction model, in order to assess whether a reasonable combination of changes (that is not excluded by current experimental data) could alleviate the observed tension between the model predictions and observed features of extensive air showers induced by ultra-high energy cosmic rays (UHECR). Previously, we have studied the effects of such changes on proton-initiated showers. Because a multitude of experimental data suggest that the primary composition of the UHECR is mixed, we have expanded the modification procedure to include nuclear projectiles in a consistent way based on the superposition model, in a similar manner as was used in the previous studies carried out using one-dimensional simulation methods. As we are using a fully three-dimensional approach, we can quantify the effects of the changes on both longitudinal and lateral features of the showers. With the inclusion of nuclear projectiles, we can study the impact of the changes on observable quantities for realistic primary beams as well as on the determination of the primary composition from data under the assumption of the modified hadronic interactions.

HERMES (High Energy Rapid Modular Ensemble of Satellites) pathfinder is an in-orbit demonstration consisting of a constellation of six 3U nano-satellites hosting simple but innovative detectors for the monitoring of cosmic high-energy transients. The main objective of HERMES Pathfinder is to prove that accurate position of high-energy cosmic transients can be obtained using miniaturized hardware. The transient position is obtained by studying the delay time of arrival of the signal to different detectors hosted by nano-satellites on low Earth orbits. To this purpose, the goal is to achive an overall accuracy of a fraction of a micro-second. In this context, we need to develop novel tools to fully exploit the future scientific data output of HERMES Pathfinder. In this paper, we introduce a new framework to assess the background count rate of a space-born, high energy detector; a key step towards the identification of faint astrophysical transients. We employ a Neural Network (NN) to estimate the background lightcurves on different timescales. Subsequently, we employ a fast change-point and anomaly detection technique to isolate observation segments where statistically significant excesses in the observed count rate relative to the background estimate exist. We test the new software on archival data from the NASA Fermi Gamma-ray Burst Monitor (GBM), which has a collecting area and background level of the same order of magnitude to those of HERMES Pathfinder. The NN performances are discussed and analyzed over period of both high and low solar activity. We were able to confirm events in the Fermi/GBM catalog and found events, not present in Fermi/GBM database, that could be attributed to Solar Flares, Terrestrial Gamma-ray Flashes, Gamma-Ray Bursts, Galactic X-ray flash. Seven of these are selected and analyzed further, providing an estimate of localisation and a tentative classification.

(Abridged) We present an exploration of the internal rotation of GD 278, the first pulsating extremely low-mass white dwarf that shows rotational splittings within its periodogram. We assess the theoretical frequency splittings expected for different rotation profiles and compare them to the observed frequency splittings of GD 278. To this aim, we employ an asteroseismological model representative of the pulsations of this star, obtained by using the LPCODE stellar evolution code. We also derive a rotation profile that results from detailed evolutionary calculations carried out with the MESA stellar evolution code and use it to infer the expected theoretical frequency splittings. We found that the best-fitting solution when assuming linear profiles for the rotation of GD 278 leads to values of the angular velocity at the surface and the center that are only slightly differential, and still compatible with rigid rotation. The values of the angular velocity at the surface and the center for the simple linear rotation profiles and for the rotation profile derived from evolutionary calculations are in very good agreement. Also, the resulting theoretical frequency splittings are compatible with the observed frequency splittings, in general, for both cases. The results obtained from the different approaches followed in this work to derive the internal rotation of GD 278 agree. The fact that they were obtained employing two independent stellar evolution codes gives robustness to our results. Our results suggest only a marginally differential behavior for the internal rotation in GD 278, and considering the uncertainties involved, very compatible with the rigid case, as has been observed previously for white dwarfs and pre-white dwarfs. The rotation periods derived for this star are also in line with the values determined asteroseismologically for white dwarfs and pre-white dwarfs in general.

The LIGO/Virgo collaboration published the catalogs GWTC-1, GWTC-2.1 and GWTC-3 containing candidate gravitational-wave (GW) events detected during its runs O1, O2 and O3. These GW events can be possible sites of neutrino emission. In this paper, we present a search for neutrino counterparts of 90 GW candidates using IceCube DeepCore, the low-energy infill array of the IceCube Neutrino Observatory. The search is conducted using an unbinned maximum likelihood method, within a time window of 1000 s and uses the spatial and timing information from the GW events. The neutrinos used for the search have energies ranging from a few GeV to several tens of TeV. We do not find any significant emission of neutrinos, and place upper limits on the flux and the isotropic-equivalent energy emitted in low-energy neutrinos. We also conduct a binomial test to search for source populations potentially contributing to neutrino emission. We report a non-detection of a significant neutrino-source population with this test.

The smallest extreme-ultraviolet (EUV) brightening events that were detected so far, called campfires, have recently been uncovered by the High Resolution EUV telescope (HRIEUV), which is part of the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter. HRIEUV has a broad bandpass centered at 17.4 nm that is dominated by Fe ix and Fe x emission at about 1 MK. We study the thermal properties of EUI brightening events by simultaneously observing their responses at different wavelengths using spectral data from the Spectral Imaging of the Coronal Environment (SPICE) also on board Solar Orbiter and imaging data from EUI. We studied three EUI brightenings that were identified in HRIEUV data that lie within the small areas covered by the slit of the SPICE EUV spectrometer. We obtained the line intensities of the spectral profiles by Gaussian fitting. These diagnostics were used to study the evolution of the EUI brightenings over time at the different line-formation temperatures. We find that (i) the detection of these EUI brightenings is at the limit of the SPICE capabilities. They could not have been independently identified in the data without the aid of HRIEUV observations. (ii) Two of these EUI brightenings with longer lifetimes are observed up to Ne viii temperatures (0.6 MK). (iii) All of the events are detectable in O vi (0.3 MK), and the two longer-lived events are also detected in other transition region (TR) lines. (iv) In one case, we observe two peaks in the intensity light curve of the TR lines that are separated by 2.7 min for C iii and 1.2 min for O vi. The Ne viii intensity shows a single peak between the two peak times of the TR line intensity. Spectral data from SPICE allow us to follow the thermal properties of EUI brightenings. Our results indicate that at least some EUI brightenings barely reach coronal temperatures.

We study the low corona evolution of the `Cartwheel' coronal mass ejection (CME; 2008-04-09) by reconstructing its 3D path and modeling it with magneto-hydrodynamic simulations. This event exhibits a double-deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The `Cartwheel CME' travels toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the pre-eruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight on the motion exhibited by the 2008-04-09 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR towards the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Using very deep, high spectral resolution data from the SAMI Integral Field Spectrograph we study the stellar population properties of a sample of dwarf galaxies in the Fornax Cluster, down to a stellar mass of $10^{7}$ M$_{\odot}$, which has never been done outside the Local Group. We use full spectral fitting to obtain stellar population parameters. Adding massive galaxies from the ATLAS$^{3D}$ project, which we re-analysed, and the satellite galaxies of the Milky Way, we obtained a galaxy sample that covers the stellar mass range $10^{4}$ to $10^{12} M_{\odot}$. Using this large range we find that the mass - metallicity relation is not linear. We also find that the [$\alpha$/Fe]-stellar mass relation of the full sample shows a U-shape, with a minimum in [$\alpha$/Fe] for masses between $10^{9}-10^{10} M_{\odot}$. The relation between [$\alpha$/Fe] and stellar mass can be understood in the following way: When the faintest galaxies enter the cluster environment, a rapid burst of star formation is induced, after which the gas content is blown away by various quenching mechanisms. This fast star formation causes high [$\alpha$/Fe] values, like in the Galactic halo. More massive galaxies will manage to keep their gas longer and form several bursts of star formation, with lower [$\alpha$/Fe] as a result. For massive galaxies, stellar populations are regulated by internal processes, leading to [$\alpha$/Fe] increasing with mass. We confirm this model by showing that [$\alpha$/Fe] correlates with clustercentric distance in three nearby clusters, and also in the halo of the Milky Way.

In canonical single-field inflation, the production of primordial black holes (PBH) requires a transient violation of the slow-roll condition. The transient ultra slow-roll inflation is an example of such scenario, and more generally, one can consider the transient constant-roll inflation. We investigate the squeezed bispectrum in the transient constant-roll inflation, and find that the Maldacena's consistency relation holds for a sufficiently long-wavelength mode, whereas it is violated for modes around the peak scale for the non-attractor case. We also demonstrate how the one-loop corrections are modified compared to the case of the transient ultra slow-roll inflation, focusing on representative one-loop terms orgiginating from a time derivative of the second slow-roll parameter in the cubic action. We find that the perturbativity requirement on those terms does not rule out the production of PBH from the transient constant-roll inflation. Therefore, it is a simple counterexample of the recently claimed no-go theorem of PBH production from single-field inflation.

Mini-filament eruptions are one of the most common small-scale transients in the solar atmosphere. However, their eruption mechanisms are still not understood thoroughly. Here, with a combination of 174 A images of high spatio-temporal resolution taken by the Extreme Ultraviolet Imager on board Solar Orbiter and images of the Atmospheric Imaging Assembly on board Solar Dynamics Observatory, we investigate in detail an erupting mini-filament over a weak magnetic field region on 2022 March 4. Two bright ribbons clearly appeared underneath the erupting mini-filament as it quickly ascended, and subsequently, some dark materials blew out when the erupting mini-filament interacted with the outer ambient loops, thus forming a blowout jet characterized by a widening spire. At the same time, multiple small bright blobs of 1-2 Mm appeared at the interaction region and propagated along the post-eruption loops toward the footpoints of the erupting fluxes at a speed of ~ 100 km/s. They also caused a semi-circular brightening structure. Based on these features, we suggest that the mini-filament eruption first experiences internal and then external reconnection, the latter of which mainly transfers mass and magnetic flux of the erupting mini-filament to the ambient corona.

The Color Magnitude Diagram (CMD) morphology of the "extended" main sequence turnoff (eMSTO) and upper main sequence (MS) of the intermediate age ($\lesssim 2$ Gyr) Large Magellanic Cloud Cluster NGC 1783 shows the presence of a small group of UV-dim stars, that, in the ultraviolet Hubble Space Telescope filters, are located at colors on the red side of the typical "fan" shape displayed by the eMSTO. We model the UV-dim stars by assuming that some of the stars which would intrinsically be located on the left side of the eMSTO are obscured by a ring of dust due to grain condensation at the periphery of the excretion disc expelled when they spin at the high rotation rates typical of stars in the Be stage. A reasonably low optical depth at 10$\mu$ is necessary to model the UV-dim group. Introduction of dust in the interpretation of the eMSTO may require a substantial re-evaluation of previous conclusions concerning the role of age and/or rotation spreads in the MC clusters: the entire eMSTO can be populated by dusty stars, and the reddest UV-dim stars simply represents the tail of the distribution with both maximum obscuration and the dust ring seen along the line of sight. The model stars having higher rotational projected velocity ($v \sin$ i) are predicted to be preferentially redder than the slowly-rotating stars. The mass loss responsible for the dust may also cause the non-monotonic distribution of stars in the upper main sequence, with two peaks and gaps showing up in the UV CMD.

How do the characteristics of starspots influence the triggering of stellar flares? Here we investigate the activity of two K-type stars, similar in every way from mass to rotation periods and planetary systems. Both stars exhibit about a hundred spots, however, Kepler-411 produced 65 superflares, while Kepler-210 presented none. The spots of both stars were characterized using the planetary transit mapping technique, which yields the intensity, temperature, and radius of starspots. The average radius was $(17\pm7) \times 10^3$ km and $(58 \pm 23) \times 10^3$ km, while the intensity ratio with respect to the photosphere was $(0.35\pm0.24)$ $I_{c}$ and $(0.64\pm0.15)$ $I_{c}$, and the temperature was $(3800 \pm 700)$ K and $(4180 \pm 240)$ K for spots of Kepler-411 and Kepler-210, respectively. Therefore, spots on the star with no superflares, Kepler-210, are mostly larger, less dark, and warmer than those on the flaring star, Kepler-411. This may indicate magnetic fields with smaller magnitude and complexity of the spots on Kepler-210 when compared to those on Kepler-411. Thus, starspot area appears not to be the main culprit of superflares triggering. Perhaps the magnetic complexity of active regions is more important.

The young planetary nebula NGC 6302 is known to exhibit Raman-scattered He II features at 6545 and 4851 Angstrom. These features are formed through inelastic scattering of He II$\lambda\lambda$ 1025 and 972 with hydrogen atoms in the ground state, for which the cross sections are $1.2 \times 10^{-21}$ and $1.4\times 10^{-22} {\rm\ cm^2}$, respectively. We investigate the spectrum of NGC 6302 archived in the ESO Science Portal. Our Gaussian line fitting analysis shows that the Raman-scattered He II features are broader and more redshifted than the hypothetical model Raman features that would be formed in a cold static H I medium. We adopt a simple scattering geometry consisting of a compact He II emission region surrounded by a H I medium to perform Monte Carlo simulations using the radiative transfer code ${\it STaRS}$. Our simulations show that the H I region is characterized by the H I column density $N_{\rm HI}=3\times 10^{21}{\rm\ cm^{-2}}$ with the random speed component $v_{\rm ran}=10{\rm\ km\ s^{-1}}$ expanding with a speed $v_{\rm exp}= 13{\rm\ km\ s^{-1}}$ from the He II emission region. Based on our best fit parameters, we estimate the H I mass of the neutral medium $M_{\rm HI} \simeq 1.0\times 10^{-2}\ {\rm M_\odot}$, pointing out the usefulness of Raman He II spectroscopy as a tool to trace H I components.

The triatomic hydrogen ion H3+ is one of the most important species for the gas phase chemistry of the interstellar medium. Observations of H3+ are used to constrain important physical and chemical parameters of interstellar environments. However, the temperatures inferred from the two lowest rotational states of H3+ in diffuse lines of sight - typically the only ones observable - appear consistently lower than the temperatures derived from H2 observations in the same sightlines. All previous attempts at modelling the temperatures of H3+ in the diffuse interstellar medium failed to reproduce the observational results. Here we present new studies, comparing an independent master equation for H3+ level populations to results from the Meudon PDR code for photon dominated regions. We show that the populations of the lowest rotational states of H3+ are strongly affected by the formation reaction and that H3+ ions experience incomplete thermalisation before their destruction by free electrons. Furthermore, we find that for quantitative analysis more than two levels of H3+ have to be considered and that it is crucial to include radiative transitions as well as collisions with H2. Our models of typical diffuse interstellar sightlines show very good agreement with observational data, and thus they may finally resolve the perceived temperature difference attributed to these two fundamental species.

We use subhalo abundance and age distribution matching to create magnitude-limited mock galaxy catalogs at $z\sim0.43$, $0.52$, and $0.63$ with $z$-band and $3.4$ micron $W1$-band absolute magnitudes and ${r-z}$ and ${r-W1}$ colors. From these magnitude-limited mocks we select mock luminous red galaxy (LRG) samples according to the $(r-z)$-based (optical) and $(r-W1)$-based (infrared) selection criteria for the LRG sample of the Dark Energy Spectroscopic Instrument (DESI) Survey. Our models reproduce the number densities, luminosity functions, color distributions, and projected clustering of the DESI Legacy Surveys that are the basis for DESI LRG target selection. We predict the halo occupation statistics of both optical and IR DESI LRGs at fixed cosmology, and assess the differences between the two LRG samples. We find that IR-based SHAM modeling represents the differences between the optical and IR LRG populations better than using the $z$-band, and that age distribution matching overpredicts the clustering of LRGs, implying that galaxy color is uncorrelated with halo age in the LRG regime. Both the optical and IR DESI LRG target selections exclude some of the most luminous galaxies that would appear to be LRGs based on their position on the red sequence in optical color-magnitude space. Both selections also yield populations with a non-trivial LRG-halo connection that does not reach unity for the most massive halos. We find the IR selection achieves greater completeness ($\gtrsim 90\%$) than the optical selection across all redshift bins studied.

Computation, if treated as a set of physical processes that act on information represented by states of matter, encompasses biological systems, digital systems, and other constructs, and may be a fundamental measure of living systems. The opportunity for biological computation, represented in the propagation and selection-driven evolution of information-carrying organic molecular structures, has been partially characterized in terms of planetary habitable zones based on primary conditions such as temperature and the presence of liquid water. A generalization of this concept to computational zones is proposed, with constraints set by three principal characteristics: capacity, energy, and instantiation (or substrate). Computational zones naturally combine traditional habitability factors, including those associated with biological function that incorporate the chemical milieu, constraints on nutrients and free energy, as well as element availability. Two example applications are presented by examining the fundamental thermodynamic work efficiency and Landauer limit of photon-driven biological computation on planetary surfaces and of generalized computation in stellar energy capture structures (a.k.a. Dyson structures). It is shown that computational zones involving nested structures or substellar objects could manifest unique observational signatures as cool far-infrared emitters. While this is an entirely hypothetical example, its simplicity offers a useful, complementary introduction to computational zones.

We present the detection of ethanol (C2H5OH), acetone (CH3COCH3), and propanal (C2H5CHO) toward the cyanopolyyne peak of TMC-1. These three O-bearing complex organic molecules are known to be present in warm interstellar clouds, but had never been observed in a starless core. The addition of these three new pieces to the puzzle of complex organic molecules in cold interstellar clouds stresses the rich chemical diversity of cold dense cores in stages prior to the onset of star formation. The detections of ethanol, acetone, and propanal were made in the framework of QUIJOTE, a deep line survey of TMC-1 in the Q band that is being carried out with the Yebes 40m telescope. We derive column densities of (1.1 +/- 0.3)e12 cm-2 for C2H5OH, (1.4 +/- 0.6)e11 cm-2 for CH3COCH3, and (1.9 +/- 0.7)e11 cm-2 for C2H5CHO. The formation of these three O-bearing complex organic molecules is investigated with the aid of a detailed chemical model which includes gas and ice chemistry. The calculated abundances at a time around 2e5 yr are in reasonable agreement with the values derived from the observations. The formation mechanisms of these molecules in our chemical model are as follows. Ethanol is formed on grains by addition of atomic carbon on methanol followed by hydrogenation and non-thermal desorption. Acetone and propanal are produced by the gas-phase reaction between atomic oxygen and two different isomers of the C3H7 radical, where the latter follows from the hydrogenation of C3 on grains followed by non-thermal desorption. A gas-phase route involving the formation of (CH3)2COH+ through several ion-neutral reactions followed by its dissociative recombination with electrons do also contribute to the formation of acetone.

Lithium-rich giant stars are rare and their existence challenges our understanding of stellar structure and evolution. We profit from the high-quality sample gathered with HARPS and UVES, in order to search for Li-rich giants and to identify the Li enrichment mechanisms responsible. We derive stellar parameters for 247 stars belonging to 32 open clusters, with 0.07 Ga < ages < 3.6 Ga. We employed the spectral synthesis technique code FASMA for the abundance analysis of 228 stars from our sample. We also determined ages, distances, and extinction using astrometry and photometry from Gaia and PARSEC isochrones to constrain their evolutionary stage. Our sample covers a wide range of stellar masses from 1 to more than 6 solar masses where the majority of the masses are above 2 solar masses. We have found 14 canonical Li-rich giant stars which have experienced the first dredge-up. This corresponds to 6% of our total sample, which is higher than what is typically found for field stars. Apart from the canonical limit, we use the maximum Li abundance of the progenitor stars as a criterion for Li enrichment. We find Li enhancement also among eight stars which have passed the first dredge up and show strong Li lines based on the fact that stars at the same evolutionary stage in the same cluster have significantly different Li abundances. We confirm that giants with higher Li abundance correspond to a higher fraction of fast-rotating giants, suggesting a connection between Li enhancement and stellar rotation as predicted by stellar models. Our Li-rich giants are found in various evolutionary stages implying that no unique Li production mechanism is responsible for Li enrichment but rather different intrinsic or external mechanisms can be simultaneously at play.

The goal of this thesis is twofold; introduce the fundamentals of Bayesian inference and computation focusing on astronomical and cosmological applications, and present recent advances in probabilistic computational methods developed by the author that aim to facilitate Bayesian data analysis for the next generation of astronomical observations and theoretical models. The first part of this thesis familiarises the reader with the notion of probability and its relevance for science through the prism of Bayesian reasoning, by introducing the key constituents of the theory and discussing its best practices. The second part includes a pedagogical introduction to the principles of Bayesian computation motivated by the geometric characteristics of probability distributions and followed by a detailed exposition of various methods including Markov chain Monte Carlo (MCMC), Sequential Monte Carlo (SMC), and Nested Sampling (NS). Finally, the third part presents two novel computational methods (Ensemble Slice Sampling and Preconditioned Monte Carlo) and their respective software implementations (zeus and pocoMC). [abridged]

We present a new catalogue of the high-mass X-ray binaries (HMXBs) in the Galaxy improving upon the most recent such catalogue. We include new HMXBs discovered since aforementioned publication and revise the classification for several objects previously considered HMXBs or candidates. The catalogue includes both basic information such as source names, coordinates, types, and more detailed data such as distance and X-ray luminosity estimates, binary system parameters and other characteristic properties of 169 HMXBs, together with appropriate references to the literature. Finding charts in several bands from infra-red to hard X-rays are also included for each object. The aim of this catalogue is to provide the reader a list of all currently known Galactic HMXBs with some basic information on both compact objects and non-degenerate counterpart properties (where available). We also include objects tentatively classified as HXMBs in the literature and give a brief motivation for the classifcation in each relevant case. The catalogue is compiled based on a search of known HMXBs and candidates in all commonly available databases and literature published before 31 October 2022. Relevant properties in the optical and other bands were collected for all objects either from the literature or using the data provided by large-scale surveys. In the later case, the counterparts in each individual survey were found by cross-correlating positions of identified HMXBs with relevant databases. An up-to date catalogue of Galactic HMXBs is presented to facilitate research in this area. An attempt was made to collect a larger set of relevant HMXB properties in a more uniform way compared to previously published works.

We present a new catalogue of low-mass X-ray binaries (LMXBs) in the Galaxy. The catalogue contains source names, coordinates, source types, fluxes, distances, system parameters, and other characteristic properties of 348 LMXBs, including LMXBs that were newly discovered or re-classified since the latest releases of the catalogues by Liu et al. (2007) and Ritter and Kolb (2003). The aim of this catalogue is to provide a list of all currently known Galactic objects identified as LMXBs with some basic information on each system (including X-ray and optical/IR properties where possible). Literature published before March 2023 has, as far as possible, been taken into account when compiling this information. References for all reported properties as well as object finding charts in several energy bands are provided as part of the catalogue. We plan to update the catalogue regularly, in particular to reflect new objects discovered in the ongoing large scale surveys such as Gaia and eROSITA.

Dusty magnetized structures observable in the far-infrared (FIR) at high Galactic latitudes are ubiquitous and found to be closely related to HI filaments with coherent velocity structures. Considering dimensionless morphological characteristics based on Minkowski functionals, we determine the distribution of filamentarities $F$ and aspect ratios $A$ for these structures. Our data are based on Planck FIR and HI4PI HI observations. Filaments have previously been extracted by applying the Hessian operator. We trace individual filamentary structures along the plane of the sky and determine $A$ and $F$. Filaments in the diffuse interstellar medium (ISM) are seldom isolated structures, but are rather part of a network of filaments with a well-defined, continuous distribution in $A$ and $F$. This distribution is self-replicating, and the merger or disruption of individual filamentary structures leads only to a repositioning of the filament in $A$ and $F$ without changing the course of the distribution. FIR and HI filaments identified at high Galactic latitudes are a close match to model expectations for narrow filaments with approximately constant widths. This distribution is continuous without clear upper limits on the observed aspect ratios. Filaments are associated with enhanced column densities of CO-dark $H_2$. Radial velocities along the filaments are coherent and mostly linear with typical dispersions of $\Delta v_{\mathrm{LSR}} = 5.24 $ km/s. The magnetic field strength in the diffuse turbulent ISM scales with hydrogen volume density as $B \propto n_{\mathrm{H}}^{0.58} $. At high Galactic latitudes, we determine an average turbulent magnetic field strength of $\langle \delta B \rangle = 5.3 ~\mu$G and an average mean strength of the magnetic field in the plane of the sky of $\langle B_{\mathrm{POS}} \rangle = 4.4 ~\mu$G.

Quasiperiodic eruptions (QPEs) from low-mass galaxy centres may result from accretion from a white dwarf in a very eccentric orbit about the central massive black hole. Evolution under gravitational radiation losses reduces the separation and eccentricity. I note that below a critical eccentricity $e_{\rm crit} \simeq 0.97$, the accretion disc's viscous timescale at pericentre passage is probably longer than the orbital period $P$, and periodic eruptive behaviour is no longer possible. These QPE descendant systems (QPEDs) are then likely to produce quasiperiodic oscillations (QPOs) rather than eruptions, varying more smoothly over the orbital cycle, with duty cycles $\sim1$. I identify 2XMM J123103.2+110648 ($P \simeq 3.9$~hr) and (more tentatively) RE J1034+396 ($P \simeq 1$~hr) as candidate systems of this type, and find agreement with their deduced eccentricities $e < e_{\rm crit}$. The absence of eruptions and the lower accretion luminosities resulting from the smaller gravitational radiation losses may make QPED systems harder to discover. Ultimately they must evolve to have viscous times much longer than the orbital period, and either remain steady, or possibly have infrequent but large outbursts. The latter systems would be massive analogues of the soft X-ray transients produced by low stellar-mass X-ray binaries.

The TESS mission detected a companion orbiting TIC 71268730, categorized it as a planet candidate, and designated the system TOI-5375. Our follow-up analysis using radial velocity data from the Habitable-zone Planet Finder (HPF), photometric data from Red Buttes Observatory (RBO), and speckle imaging with NN-EXPLORE Exoplanet Stellar Speckle Imager (NESSI) determined that the companion is a very low mass star (VLMS) near the hydrogen-burning mass limit with a mass of 0.080$\pm{0.002} M_{\Sun}$ ($83.81\pm{2.10} M_{J}$), a radius of 0.1114$^{+0.0048}_{-0.0050} R_{\Sun}$ (1.0841$^{0.0467}_{0.0487} R_{J}$), and brightness temperature of $2600\pm{70}$ K. This object orbits with a period of 1.721553$\pm{0.000001}$ days around an early M dwarf star ($0.62\pm{0.016}M_{\Sun}$). TESS photometry shows regular variations in the host star's TESS light curve, which we interpreted as activity-induced variation of $\sim$2\%, and used this variability to measure the host star's stellar rotation period of 1.9716$^{+0.0080}_{-0.0083}$ days. The TOI-5375 system provides tight constraints on stellar models of low-mass stars at the hydrogen-burning limit and adds to the population in this important region.

We have investigated, both analytically and numerically, accreting supermassive black hole binaries as they inspiral due to gravitational radiation to elucidate the decoupling of binaries from their disks and inform future multi-messenger observations of these systems. Our numerical studies evolve equal-mass binaries from initial separations of $100 GM/c^2$ until merger, resolving scales as small as $\sim0.04 GM/c^2$, where $M$ is the total binary mass. Our simulations accurately capture the point at which the orbital evolution of each binary decouples from that of their circumbinary disk, and precisely resolve the flow of gas throughout the inspiral. We demonstrate analytically and numerically that timescale-based predictions overestimate the binary separations at which decoupling occurs by factors of $\sim3$, and illustrate the utility of a velocity-based decoupling criterion. High-viscosity ($\nu\gtrsim0.03 GM/c$) circumbinary systems decouple late ($a_b\lesssim 15 GM/c^2$) and have qualitatively similar morphologies near merger to circumbinary systems with constant binary separations. Lower-viscosity circumbinary disks decouple earlier and exhibit qualitatively different accretion flows, which lead to precipitously decreasing accretion onto the binary. If detected, such a decrease may unambiguously identify the host galaxy of an ongoing event within a LISA error volume. We illustrate how accretion amplitude and variability evolve as binaries gradually decouple from their circumbinary disks, and where decoupling occurs over the course of binary inspirals in the LISA band. We show that, even when dynamically negligible, gas may leave a detectable imprint on the phase of gravitational waves.

The radiation environment over the African continent, at aviation altitudes, remains mostly uncharacterized and unregulated. In this paper we present initial measurements made by a newly developed active dosimeter on-board long-haul flights between South Africa and Germany. Based on these initial tests, we believe that this low-cost and open-source dosimeter is suitable for continued operation over the Africa continent and can provide valuable long-term measurements to test dosimteric models and inform aviation policy

We develop a formalism to describe electron ejections from graphene-like targets by dark matter (DM) scattering for general forms of scalar and spin 1/2 DM-electron interactions and compare their applicability and accuracy within the density functional theory (DFT) and tight binding (TB) approaches. This formalism allows for accurate prediction of the daily modulation signal expected from DM in upcoming direct detection experiments employing graphene sheets as the target material. A key result is that the physics of the graphene sheet and that of the DM and the ejected electron factorise, allowing for the rate of ejections from all forms of DM to be obtained with a single graphene response function. We perform a comparison between the TB and DFT approaches to modeling the initial state electronic wavefunction within this framework, with DFT emerging as the more self-consistent and reliable choice due to the challenges in the embedding of an appropriate atomic contribution into the TB approach.

We use a formalism that describes electron ejections from graphene-like targets by dark matter (DM) scattering for general forms of scalar and spin 1/2 DM-electron interactions in combination with state-of-the-art density functional calculations to produce predictions and reach estimates for various possible carbon-based detector designs. Our results indicate the importance of a proper description of the target electronic structure. In addition, we find a strong dependence of the predicted observed signal for different DM candidate masses and interaction types on the detailed geometry and design of the detector. Combined with directional background vetoing, these dependencies will enable the identification of DM particle properties once a signal has been established.

The collection of gravitational waves (GWs) that are either too weak or too numerous to be individually resolved is commonly referred to as the gravitational-wave background (GWB). A confident detection and model-driven characterization of such a signal will provide invaluable information about the evolution of the Universe and the population of GW sources within it. We present a new, user-friendly Python--based package for gravitational-wave data analysis to search for an isotropic GWB in ground--based interferometer data. We employ cross-correlation spectra of GW detector pairs to construct an optimal estimator of the Gaussian and isotropic GWB, and Bayesian parameter estimation to constrain GWB models. The modularity and clarity of the code allow for both a shallow learning curve and flexibility in adjusting the analysis to one's own needs. We describe the individual modules which make up {\tt pygwb}, following the traditional steps of stochastic analyses carried out within the LIGO, Virgo, and KAGRA Collaboration. We then describe the built-in pipeline which combines the different modules and validate it with both mock data and real GW data from the O3 Advanced LIGO and Virgo observing run. We successfully recover all mock data injections and reproduce published results.

We explore the impact of choosing different sets of Time-Delay Interferometry (TDI) variables for detecting and reconstructing Stochastic Gravitational Wave Background (SGWB) signals and estimating the instrumental noise in LISA. Most works in the literature build their data analysis pipelines relying on a particular set of TDI channels, the so-called AET variables, which are orthogonal under idealized conditions. By relaxing the assumption of a perfectly equilateral LISA configuration, we investigate to which degree these channels remain orthogonal and compare them to other TDI channels. We show that different sets of TDI variables are more robust under perturbations of the perfect equilateral configuration, better preserving their orthogonality and, thus, leading to a more accurate estimate of the instrumental noise. Moreover, we investigate the impact of considering the noise levels associated with each instrumental noise source to be independent of one another, generalizing the analysis from two to twelve noise parameters. We find that, in this scenario, the assumption of orthogonality is broken for all the TDI variables, leading to a misestimation of measurement error for some of the noise parameters. Remarkably, we find that for a flat power-law signal, the reconstruction of the signal parameters is nearly unaffected in these various configurations.

The near-zone ``Love'' symmetry resolves the naturalness issue of black hole Love number vanishing with $\text{SL}\left(2,\mathbb{R}\right)$ representation theory. Here, we generalize this proposal to $5$-dimensional asymptotically flat and doubly spinning (Myers-Perry) black holes. We consider the scalar response of Myers-Perry black holes and extract its static scalar Love numbers. In agreement with the naturalness arguments, these Love numbers are, in general, non-zero and exhibit logarithmic running unless certain resonant conditions are met; these conditions include new cases with no previously known analogs. We show that there exist two near-zone truncations of the equations of motion that exhibit enhanced $\text{SL}\left(2,\mathbb{R}\right)$ Love symmetries that explain the vanishing of the static scalar Love numbers in the resonant cases. These Love symmetries can be interpreted as local $\text{SL}\left(2,\mathbb{R}\right)\times\text{SL}\left(2,\mathbb{R}\right)$ near-zone symmetries spontaneously broken down to global $\text{SL}\left(2,\mathbb{R}\right)\times U\left(1\right)$ symmetries by the periodic identification of the azimuthal angles. We also discover an infinite-dimensional extension of the Love symmetry into $\text{SL}\left(2,\mathbb{R}\right)\ltimes\hat{U}\left(1\right)_{\mathcal{V}}^2$ that contains both Love symmetries as particular subalgebras, along with a family of $\text{SL}\left(2,\mathbb{R}\right)$ subalgebras that reduce to the exact near-horizon Myers-Perry black hole isometries in the extremal limit. Finally, we show that the Love symmetries acquire a geometric interpretation as isometries of subtracted (effective) black hole geometries that preserve the internal structure of the black hole and interpret these non-extremal $\text{SL}\left(2,\mathbb{R}\right)$ structures as remnants of the enhanced isometry of the near-horizon extremal geometries.

We discuss the electromagnetic collisions of high energy protons, pions and kaons with atmospheric nuclei. In particular, we use the equivalent photon approximation to estimate (i) the diffractive collisions where the projectile scatters inelastically off a nucleus, and (ii) the usual radiative processes (bremsstrahlung, pair production and photonuclear interactions) of these charged hadrons in the air. We then include the processes in the simulator AIRES and study how they affect the longitudinal development of extensive air showers. For $10^{9-11}$ GeV proton primaries we find that they introduce a very small reduction (below 1%) in the average value of both $X_{\rm max}$ and $\Delta X_{\rm max}$. At a given shower age (relative slant depth from $X_{\rm max}$), these electromagnetic processes do not change significantly the number of muons or the total energy carried by electrons and photons, increased by just 1% the muon-to-($\gamma+e$) signal at the ground level.