Papers for Thursday, Sep 28 2023

The number of gravitational arcs systems detected is increasing quickly and should even increase at a faster rate in the near future. This wealth of new gravitational arcs requires the development of a purely automated method to reconstruct the lens and source. A general reconstruction method based on the singular perturbative approach is proposed in this paper. This method generates a lens and source reconstruction directly from the gravitational arc image. The method is fully automated and works in two steps. The first step is to generate a guess solution based on the circular solution in the singular perturbative approach. The second step is to break the sign degeneracy and to refine the solution by using a general source model. The refinement of the solution is conducted step by step to avoid the source-lens degeneracy issue. One important asset of this automated method is that the lens solution is written in universal terms which allows the computation of statistics. Considering the large number of lenses which should be available in the near future this ability to compute un-biased statistics is an important asset.

A long-lived magnetar, potentially originating from a binary neutron star system, has been proposed to explain the extended emission observed in certain short-duration gamma-ray bursts (sGRBs), and is posited as a potential central engine to power the engine-fed kilonovae. Previously, the process by which energy is injected into the surrounding ejecta/jet was widely believed to be nearly isotropic. In this study, we employ special relativity magnetohydrodynamic (SRMHD) simulations to investigate the wind injection process from a magnetar central engine. We explore the dynamics and energy distribution within the system and found that the parameter $\alpha=u_{\rm A}/u_{\rm MWN}$ can be used to indicate the collimation of the magnetar wind energy injection, where $u_{\rm A}$ is the local Alfven four-speed and $u_{\rm MWN}$ is the four-speed of the magnetar wind nebular (MWN) formed from wind-ejecta collision. A significant portion of the injected energy from the magnetar spin-down wind will be channeled to the jet axis due to collimation within the MWN. Achieving isotropic energy injection requires a significantly small $\alpha$ that necessitates either an ultra-relativistic expanding MWN or an extremely low magnetization MWN, both of which are challenging to attain in sGRBs. Consequently, a considerably reduced energy budget (by a factor of up to 10) is anticipated to be injected into the ejecta for engine-fed kilonovae. Engine-fed kilonovae would appear fainter than originally anticipated.

The third data release of Gaia was the first to include orbital solutions assuming non-single stars. Here, we apply the astrometric triage technique of Shahaf et al. 2019 to identify binary star systems with companions that are not single main-sequence stars. Gaia's synthetic photometry of these binaries is used to distinguish between systems likely to have white-dwarf companions and those that may be hierarchical triples. The study uncovered a population of nearly 3200 binaries, characterised by orbital separations on the order of an astronomical unit, in which the faint astrometric companion is probably a white dwarf. Remarkably, over 110 of these systems exhibit significant ultraviolet excess flux, confirming this classification and, in some cases, indicating their relatively young cooling ages. We show that the sample is not easily reproduced by binary population synthesis codes. Therefore, it challenges current binary evolution models, offering a unique opportunity to gain insights into the processes governing white-dwarf formation, binary evolution, and mass transfer.

Ly$\alpha$ emission nebulae are ubiquitous around high-z galaxies and are tracers of the gaseous environment on scales out to >100 kpc. High-z radio galaxies (HzRGs, type-2 radio-loud quasars) host large scale nebulae observed in the ionised gas differ from those seen in other types of high-z quasars. In this work, we exploit MUSE observations of Lya nebulae around eight HzRGs ($2.9<z<4.5$). All the HzRGs have large scale Lya emission nebulae with seven of them extended over 100 kpc at the observed surface brightness (SB) limit. Because the emission line profiles are significantly affected by neutral hydrogen absorbers across the entire nebulae, we perform an absorption correction to infer maps of the intrinsic Lya SB, central velocity and velocity width, all at the last scattering surface of the observed Lya photons. We find: (i) The intrinsic SB radial profiles can be described by an inner exponential and an outter power law; (ii) our HzRGs have higher SB and more asymmetric nebulae than both RL and RQ type-1s; (iii) intrinsic nebula kinematics of four HzRGs show evidence of jet-driven outflows but no general trends for the whole sample; (iv) a relation between the nebula maximum extent and the offset between the AGN and the nebula centroids; (v) an alignment between radio jet position angles and the nebula morphology. All support a scenario where the orientation of the AGN has an impact on the observed nebular morphologies and resonant scattering may affect the shape of the SB profiles, nebular kinematics and relations between the observed Lya morphologies. Furthermore, we find evidence showing that the outskirts of the ionised gas nebulae may be 'contaminated' by Lya photons from nearby emission halos. Overall, this work provides results which allow us to compare Lya nebulae around various classes of quasars at and beyond Cosmic Noon. [Abridged]

We present the OGLE collection of delta Scuti stars in the Large Magellanic Cloud and in its foreground. Our dataset encompasses a total of 15,256 objects, constituting the largest sample of extragalactic delta Sct stars published so far. In the case of 12 delta Sct pulsators, we detected additional eclipsing or ellipsoidal variations in their light curves. These are the first known candidates for binary systems containing delta Sct components beyond the Milky Way. We provide observational parameters for all variables, including pulsation periods, mean magnitudes, amplitudes, and Fourier coefficients, as well as long-term light curves in the I- and V-bands collected during the fourth phase of the OGLE project. We construct the period-luminosity (PL) diagram, in which fundamental-mode and first-overtone delta Sct stars form two nearly parallel ridges. The latter ridge is an extension of the PL relation obeyed by first-overtone classical Cepheids. The slopes of the PL relations for delta Sct variables are steeper than those for classical Cepheids, indicating that the continuous PL relation for first-overtone delta Sct variables and Cepheids is non-linear, exhibiting a break at a period of approximately 0.5 d. We also report the enhancement of the OGLE collection of Cepheids and RR Lyrae stars with newly identified and reclassified objects, including pulsators contained in the recently published Gaia DR3 catalog of variable stars. As a by-product, we estimate the contamination rate in the Gaia DR3 catalogs of Cepheids and RR Lyrae variables.

We derive the stellar population parameters of 11 quiescent ultra-diffuse galaxies (UDGs) from Keck/KCWI data. We supplement these with 14 literature UDGs, creating the largest spectroscopic sample of UDGs to date (25). We find a strong relationship between their $\alpha$-enhancement and their star formation histories: UDGs that formed on very short timescales have elevated [Mg/Fe] abundance ratios, whereas those forming over extended periods present lower values. Those forming earlier and faster are overall found in high-density environments, being mostly early infalls into the cluster. No other strong trends are found with infall times. We analyze the stellar mass-metallicity, age-metallicity and [Mg/Fe]-metallicity relations of the UDGs, comparing them to other types of low mass galaxies. Overall, UDGs scatter around the established stellar mass--metallicity relations of classical dwarfs. We find that GC-rich UDGs have intermediate-to-old ages, but previously reported trends of galaxy metallicity and GC richness are not reproduced with this spectroscopic sample due to the existence of GC-rich UDGs with elevated metallicities. In addition, we also find that a small fraction of UDGs could be 'failed-galaxies', supported by their GC richness, high $\alpha$-abundances, fast formation timescales and that they follow the mass-metallicity relation of z~2 galaxies. Finally, we also compare our observations to simulated UDGs. We caution that there is not a single simulation that can produce the diverse UDG properties simultaneously, in particular the low metallicity failed-galaxy like UDGs.

Aims: We leverage the largest available Atacama Large Millimetre/submillimetre Array (ALMA) survey from the archive (A$^3$COSMOS) to study to study infrared luminosity function and dust-obscured star formation rate density of sub-millimeter/millimeter (sub-mm/mm) galaxies from $z=0.5\,-\,6$. Methods: The A$^3$COSMOS survey utilizes all publicly available ALMA data in the COSMOS field, therefore having inhomogeneous coverage in terms of observing wavelength and depth. In order to derive the luminosity functions and star formation rate densities, we apply a newly developed method that corrects the statistics of an inhomogeously sampled survey of individual pointings to those representing an unbiased blind survey. Results: We find our sample to mostly consist of massive ($M_{\star} \sim 10^{10} - 10^{12}$ $\rm M_{\odot}$), IR-bright ($L_* \sim 10^{11}-10^{13.5} \rm L_{\odot}$), highly star-forming (SFR $\sim 100-1000$ $\rm M_{\odot}$ $\rm yr^{-1}$) galaxies. We find an evolutionary trend in the typical density ($\Phi^*$) and luminosity ($L^*$) of the galaxy population, which decrease and increase with redshift, respectively. Our IR LF is in agreement with previous literature results and we are able to extend to high redshift ($z > 3$) the constraints on the knee and bright-end of the LF, derived by using the Herschel data. Finally, we obtain the SFRD up to $z\sim 6$ by integrating the IR LF, finding a broad peak from $z \sim 1$ to $z \sim 3$ and a decline towards higher redshifts, in agreement with recent IR/mm-based studies, within the uncertainties, thus implying the presence of larger quantities of dust than what is expected by optical/UV studies.

Searches for primordial non-Gaussianity in cosmological perturbations are a key means of revealing novel primordial physics. However, robustly extracting signatures of primordial non-Gaussianity from non-linear scales of the late-time Universe is an open problem. In this paper, we apply k-Nearest Neighbor cumulative distribution functions, kNN-CDFs, to the \textsc{quijote-png} simulations to explore the sensitivity of kNN-CDFs to primordial non-Gaussianity. An interesting result is that for halo samples with $M_h<10^{14}$ M$_\odot$/h, the kNN-CDFs respond to \textit{equilateral} PNG in a manner distinct from the other parameters. This persists in the galaxy catalogs in redshift space and can be differentiated from the impact of galaxy modelling, at least within the halo occupation distribution (HOD) framework considered here. kNN-CDFs are related to counts-in-cells and, through mapping a subset of the kNN-CDF measurements into the count-in-cells picture, we show that our results can be modeled analytically. A caveat of the analysis is that we only consider the HOD framework, including assembly bias. It will be interesting to validate these results with other techniques for modeling the galaxy--halo connection, e.g., (hybrid) effective field theory or semi-analytical methods.

Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfur enrichment in this nearby system and to establish benchmarks for such ratios in metal-poor galaxies. We conducted pointed observations of CS and its isotopologues toward N113, one of the most prominent star-formation regions in the LMC, utilizing the Atacama Pathfinder EXperiment 12~m telescope. We present the first robust detection of C$^{33}$S in the LMC by successfully identifying two C$^{33}$S transitions. Our measurements result in the first direct determination of the $^{34}$S/$^{33}$S isotope ratio, which is 2.0$\pm$0.2. Our comparative analysis indicates that the $^{32}$S/$^{33}$S and $^{34}$S/$^{33}$S isotope ratios are about a factor of 2 lower in the LMC than in the Milky Way. Our findings suggest that the low $^{34}$S/$^{33}$S isotope ratio in the LMC can be attributed to a combination of the age effect, low metallicity, and star formation history.

We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin for Saturn's remarkably young rings. Such an event could have been triggered a few hundred million years ago by resonant instabilities in a previous satellite system. Using high-resolution smoothed particle hydrodynamics simulations, we find that this kind of impact can produce a wide distribution of massive objects and scatter material throughout the system. This includes the direct placement of pure-ice ejecta onto orbits that enter Saturn's Roche limit, which could form or rejuvenate rings. In addition, fragments and debris of rock and ice totalling more than the mass of Enceladus can be placed onto highly eccentric orbits that would intersect with any precursor moons orbiting in the vicinity of Mimas, Enceladus, or Tethys. This could prompt further disruption and facilitate a collisional cascade to distribute more debris for potential ring formation, the re-formation of the present-day moons, and evolution into an eventual cratering population of planeto-centric impactors.

Context: Circumstellar debris disks provide insight into the formation and early evolution of planetary systems. Resolved belts in particular help to locate planetesimals in exosystems, and can hint at the presence of disk-sculpting exoplanets. Aims: We study the circumstellar environment of HD 112810 (HIP 63439), a mid-F type star in the Sco-Cen association with a significant infrared excess indicating the presence of a circumstellar debris disk. Methods: We collected five high-contrast observations of HD 112810 with VLT/SPHERE. We identified a debris disk in scattered light, and found that the debris signature is robust over a number of epochs and a variety of reduction techniques. We modelled the disk, accounting for self-subtraction and assuming that it is optically thin. Results: We find a single-belt debris disk, with a radius of 118$\pm$9au and an inclination angle of ${75.7}^{+1.1}_{-1.3}$$\deg$. This is in good agreement with the constraints from SED modelling and from a partially-resolved ALMA image of the system. No planets are detected, though planets below the detection limit ($\sim$2.6M$_\textrm{J}$ at a projected separation of 118au) could be present and could have contributed to sculpting the ring of debris. Conclusions: HD 112810 adds to the growing inventory of debris disks imaged in scattered light. The disk is faint, but the radius and the inclination of the disk are promising for follow-up studies of the dust properties.

A new era of lunar exploration has begun with participation of all major space agencies. This activity brings opportunities for revolutionary science experiments and observatories on the Moon. The idea of a lunar gravitational-wave detector was already proposed during the Apollo program. The key characteristic of the Moon is that it is seismically extremely quiet. It was also pointed out that the permanently shadowed regions at the lunar poles provide ideal conditions for gravitational-wave detection. In recent years, three different detector concepts were proposed with varying levels of technological complexity and science potential. In this paper, we confront the three concepts in terms of their observational capabilities based on a first more detailed modeling of instrumental noise. We identify important technological challenges and potential show-stoppers.

The degree of mass loss, i.e. the fraction of stars lost by globular clusters, and specifically by their different populations, is still poorly understood. Many scenarios of the formation of multiple stellar populations, especially the ones involving self-enrichment, assume that the first generation (FG) was more massive at birth than now to reproduce the current mass of the second generation (SG). This assumption implies that, during their long-term evolution, clusters lose around $90\%$ of the FG. We have tested whether such strong mass loss could take place in a massive globular cluster orbiting the Milky Way at $4\ {\rm kpc}$ from the centre and composed of two generations. We perform a series of $N$-body simulations for ${12\ \rm Gyr}$ to probe the parameter space of internal cluster properties. We have derived that, for an extended FG and a low-mass second one, the cluster loses almost $98\%$ of its initial FG mass and the cluster mass can be as much as 20 times lower after a Hubble time. Furthermore, under these conditions, the derived fraction of SG stars, $f_{\rm enriched}$, falls in the range occupied by observed clusters of similar mass ($\sim 0.6-0.8$). In general, the parameters that affect the most the degree of mass loss are the presence or not of primordial segregation, the depth of the central potential, $W_{0,FG}$, the initial mass of the SG, $M^{ini}_{SG}$, and the initial half-mass radius of the SG, $r_{h,SG}$. Higher $M^{ini}_{SG}$ have not been found to imply higher final $f_{\rm enriched}$ due to the deeper cluster potential well which slows down mass loss.

Astronomy observation is difficult in urban environments due to the background noise generated by human activities. Consequently, promoting astronomy in metropolitan areas is challenging. In this work, we propose a low-cost, educational experiment called Wok the Hydrogen (WTH) that offers opportunities for scientific observation in urban environments, specifically the observation of the $21$ cm ($f_{21} = 1420.4$ MHz) emission from neutral hydrogen in the Milky Way. We demonstrate how to construct a radio telescope using kitchenware, along with additional electronic equipment that can be easily purchased online. The total system cost is controlled within 150 dollars. We also outline the subsequent data analysis procedures for deriving the recession velocity of galactic hydrogen from the raw data. The system was tested on the campus of the Hong Kong University of Science and Technology, which is located approximately 2 km northeast of the nearest residential area with a population of 0.4 million and about 10 km east of the downtown area with a population of 2 million. We show that a precision of $\Delta v \approx \pm 20$ km s$^{-1}$ can be achieved for determining the recession velocity of neutral hydrogen with this relatively simple setup, and the precision can be further improved with longer exposure time.

In this article, we employ a machine learning (ML) approach for the estimations of four fundamental parameters, namely, the Hubble constant ($H_0$), matter ($\Omega_{0m}$), curvature ($\Omega_{0k}$) and vacuum ($\Omega_{0\Lambda}$) densities of non-flat $\Lambda$CDM model. We use $53$ Hubble parameter values measured by differential ages (DA) and baryon acoustic oscillations (BAO) techniques in the redshift interval $0.07 \leq z \leq 2.36$. We create an artificial neural network (called ParamANN) and train it with simulated values of $H(z)$ using various sets of $H_0$, $\Omega_{0m}$, $\Omega_{0k}$, $\Omega_{0\Lambda}$ parameters chosen from different and sufficiently wide prior intervals. We use a correlated noise model in the analysis. We demonstrate accurate validation and prediction by ParamANN. ParamANN provides an excellent cross-check for the validity of the $\Lambda$CDM model and alleviates the Hubble tension problem which has been reported earlier in the literature. We obtain $H_0 = 66.11 \pm 2.59$ $\rm{kmMpc^{-1}sec^{-1}}$, $\Omega_{0m} = 0.3359 \pm 0.0814$, $\Omega_{0k} = 0.0237 \pm 0.1248$ and $\Omega_{0\Lambda} = 0.6405 \pm 0.0861$ by using the trained network. These parameter values agree very well with the results of global CMB observations of Planck collaboration.

It is widely believed that the binary neutron star merger GW190425 produced a black hole promptly upon merger. Motivated by the potential association with the fast radio burst FRB 20190425A, which took place 2.5~hours after the merger, we revisit the question of the outcome of GW190425 by means of numerical relativity simulations. We show that current laboratory and astrophysical constraints on the equation of state of dense matter do not rule out the formation of a long-lived remnant. However, the formation of a stable remnant would have produced a bright kilonova, in tension with upper limits by ZTF at the location and time of FRB 20190425A. Moreover, the ejecta would have been optically thick to radio emission for days to months, preventing a putative FRB from propagating out. The predicted dispersion measure is also several orders of magnitude larger than that observed for FRB 20190425A. Our results indicate that FRB 20190425A and GW190425 are not associated. However, we cannot completely rule out the formation of a long-lived remnant, due to the incomplete coverage of the relevant sky regions. More observations of GW190425-like events, including potential upper limit, have the potential to constrain nuclear physics. To this aim, it is important that follow-up observational campaigns of gravitational wave events are informed by the properties of the source, such as their chirp mass, and we urge the LIGO-Virgo-KAGRA collaboration to promptly release them publicly.

Merging black holes (BH) are expected to produce remnants with large dimensionless spin parameters ($a_{\rm spin} \sim 0.7$). However, gravitational wave (GW) observations with LIGO/Virgo suggest that merging BH are consistent with modestly positive but not high spin ($a_{\rm spin} \sim 0.2$), causing tension with models suggesting that high mass mergers are produced by hierarchical merger channels. Some BH also show evidence for strong in-plane spin components. Here we point out that \emph{spin down} of BH due to eccentric prograde post-merger orbits within the gas of an active galactic nucleus (AGN) disk can yield BH with masses in the upper mass gap, but only modestly positive $a_{\rm spin}$, and thus observations of BH with low spin \emph{do not} rule out hierarchical models. We also point out that the fraction of BBH mergers with significant in-plane spin components is a strong test of interactions between disk binary black holes (BBH) and nuclear spheroid orbiters. Spin magnitude and spin tilt constraints from LIGO/Virgo observations of BBH are an excellent test of dynamics of black holes in AGN disks, disk properties and the nuclear clusters interacting with AGN.

Models with an interaction between dark energy and dark matter have already been studied for about twenty years. However, in this paper, we provide for the first time a general analytical solution for models with an energy transfer given by $\mathcal{E} = 3H(\xi_1 \rho_c + \xi_2 \rho_d)$. We also use a new set of age-redshift data for 114 old astrophysical objects (OAO) and constrain some special cases of this general energy transfer. We use a method inspired on artificial intelligence, known as Chaos Quantum-behaved Particle Swarm Optimization (CQPSO), to explore the parameter space and search the best fit values. We test this method under a simulated scenario and also compare with previous MCMC results and find good agreement with the expected results.

We observed Sedna, Gonggong, and Quaoar with the NIRSpec instrument on the James Webb Space Telescope (JWST). All three bodies were observed in the low-resolution prism mode at wavelengths spanning 0.7 to 5.2 {\mu}m. Quaoar was also observed at 10x higher spectral resolution from 0.97 to 3.16 {\mu}m using medium-resolution gratings. Sedna's spectrum shows a large number of absorption features due to ethane (C2H6), as well as acetylene (C2H2), ethylene (C2H4), H2O, and possibly minor CO2. Gonggong's spectrum also shows several, but fewer and weaker, ethane features, along with stronger and cleaner H2O features and CO2 complexed with other molecules. Quaoar's prism spectrum shows even fewer and weaker ethane features, the deepest and cleanest H2O features, a feature at 3.2 {\mu}m possibly due to HCN, and CO2 ice. The higher-resolution medium grating spectrum of Quaoar reveals several overtone and combination bands of ethane and methane (CH4). Spectra of all three objects show steep red spectral slopes and strong, broad absorptions between 2.7 and 3.6 {\mu}m indicative of complex organic molecules. The suite of light hydrocarbons and complex organic molecules are interpreted as the products of irradiation of methane. The differences in apparent abundances of irradiation products are likely due to their distinctive orbits, which lead to different timescales of methane retention and to different charged particle irradiation environments. In all cases, however, the continued presence of light hydrocarbons implies a resupply of methane to the surface. We suggest that these three bodies have undergone internal melting and geochemical evolution similar to the larger dwarf planets and distinct from all smaller KBOs.

This study analyzes Io's thermally detected volcanic outbursts and mini-outbursts, generally called bright transient eruptions. We examine their evolving characteristics over the history of outburst observations between the Voyager flybys in 1978 and 2022. We catalog, compare, and interpret the data of these bright transient eruptions from several spacecraft flybys and numerous ground-based observation campaigns. To test the spatiotemporal behavior of these events, we compare them to a population of randomly spaced, stochastic events with an equal likelihood of occurrence anywhere on Io's surface. We find that the aggregate of all outbursts is consistent with a random distribution across Io, whereas mini-outbursts strongly prefer the trailing hemisphere (180 to 360 W). On shorter timescales, however, outbursts show a significant change in spatiotemporal behavior before and after the year 2012. Outbursts from 1995 to 2007 favor the northern leading hemisphere, while outbursts from 2013 to 2021 favor the southern trailing hemisphere. These temporally separated clusters of outbursts are remarkably similar to Io's two primary mountainous regions, indicating that outbursts may be related to mountain-forming activity. These trends show how bright transient eruptions are distinct from Io's other forms of volcanism. These could be essential constraints to assess models of Io's interior heat transport between tidal generation and volcanic distribution.

TRAPPIST-1e is a potentially habitable terrestrial exoplanet orbiting an ultra-cool M Dwarf star and is a key target for observations with the James Webb Space Telescope (JWST). One-dimensional photochemical modelling of terrestrial planetary atmospheres has shown the importance of the incoming stellar UV flux in modulating the concentration of chemical species, such as O$_3$ and H$_2$O. In addition, three-dimensional (3D) modelling has demonstrated anisotropy in chemical abundances due to transport in tidally locked exoplanet simulations. We use the Whole Atmosphere Community Climate Model Version 6 (WACCM6), a 3D Earth System Model, to investigate how uncertainties in the incident UV flux, combined with transport, affect observational predictions for TRAPPIST-1e (assuming an initial Earth-like atmospheric composition). We use two semi-empirical stellar spectra for TRAPPIST-1 from the literature. The UV flux ratio between them can be as large as a factor of 5000 in some wavelength bins. Consequently, the photochemically-produced total O$_3$ columns differ by a factor of 26. Spectral features of O$_3$ in both transmission and emission spectra vary between these simulations (e.g. differences of 19 km in transmission spectra effective altitude for O$_3$ at 0.6 $\mu$m). This leads to potential ambiguities when interpreting observations, including overlap with scenarios that assume alternative O$_2$ concentrations. Hence, to achieve robust interpretations of terrestrial exoplanetary spectra, characterisation of the UV spectra of their host stars is critical. In the absence of such stellar measurements, atmospheric context can still be gained from other spectral features (e.g. H$_2$O), or by comparing direct imaging and transmission spectra in conjunction.

The relationship between the integrated H$\beta$ line luminosity and the velocity dispersion of the ionized gas of HII galaxies and giant HII regions represents an exciting standard candle that presently can be used up to redshifts $z\sim~4$. Locally it is used to obtain precise measurements of the Hubble constant by combining the slope of the relation obtained from nearby ($z<0.2$) HII galaxies with the zero point determined from giant HII regions belonging to an `anchor sample' of galaxies for which accurate redshift-independent distance moduli are available (e.g Cepheids, TRGBs). Through this chapter, we present a general description of the method and results obtained so far in the determination of the local value of the Hubble constant and cosmological constraints. We account for the main systematic effects associated with the method and the possibility of improvement in the future. The discussion presented here is the product of an independent approach to the standard methods used in the literature to measure distances and, although at present less precise than the latest SNIa results, it is amenable to substantial improvement.

In this contribution we study the possibility of the formation of cosmic ray ensembles (CRE) created by the interaction of ultra-high energy (UHE) photons with the magnetic field of the Sun. The lack of observation of those UHE and the difficulties for their identification given the current methodologies motivates this study. We performed simulations using the PRESHOWER program in order to simulate the expected extensive air showers which might be spatially correlated generated upon entering the Earth's atmosphere. We found characteristic features like very thing and extremely elongates cascades of secondary photons with their corresponding energies spanning the entire cosmic range spectrum. Shower footprints are as large as hundreds of kilometres. An application of this study is the scenario of gamma-ray emission from the vicinity of the Sun as a result of ultra-high energy photon cascading in the solar magnetic field in order to understand recent observations made by the HAWC and Fermi-LAT observatories.

The planetary infrared excess (PIE) technique has the potential to efficiently detect and characterize the thermal spectra of both transiting and non-transiting exoplanets. However, the technique has not been evaluated on multiplanet systems. We use the TRAPPIST-1 system as our test bed to evaluate PIE's ability to resolve multiple planets. We follow the unfolding discoveries in the TRAPPIST-1 system and examine the results from the PIE technique at every stage. We test the information gained from observations with JWST and next-generation infrared observatories like the proposed MIRECLE mission concept. We find that even in the case where only the star is known, the PIE technique would infer the presence of multiple planets in the system. The precise number inferred is dependent on the wavelength range of the observation and the noise level of the data. We also find that in such a tightly packed, multiplanet system such as TRAPPIST-1, the PIE technique struggles to constrain the semi-major axis beyond prior knowledge. Despite these drawbacks and the fact that JWST is less sensitive to the fluxes from planets g and h, with strong priors in their orbital parameters we are able to constrain their equilibrium temperatures. We conclude that the PIE technique may enable the discovery of unknown exoplanets around solar-neighborhood M dwarfs and could characterize known planets around them.

Galaxy clusters, which underwent a recent ($\leq3$ Gyr) major merger, offer a harsher environment due to the global hydrodynamical disturbance and the merger-shock heated ICM. However, the aftermath of such extreme cluster interactions on the member galaxy properties is not very well constrained. We explore the integrated star formation properties of galaxies through galaxy colours, as well as morphology buildup in three nearby ($0.04<z<0.07$) young ($\sim$0.6-1 Gyr) post-merger clusters -- A3667, A3376 and A168 -- and 7 relaxed clusters, to disentangle merger-induced post-processing signatures from the expected effects due to high-density cluster environments. Exploiting the optical spectroscopy and photometry from the OmegaWINGS survey, we find that post-merger clusters are evolved systems demonstrating uniform spiral fractions, uniform fraction of blue galaxies and constant scatter in the colour-magnitude relations, a regularity that is absent in dynamically relaxed clusters. While no clear merger-induced signatures were revealed in the global colours of galaxies, we conclude that different global star formation histories of dynamically relaxed clusters lead to considerable scatter in galaxy properties, resulting in the pre-merger cluster environment to potentially contaminate any merger-induced signal in galaxy properties. We discover red spirals to be common to both post-merger and relaxed clusters while post-merger clusters appear to host a non-negligible population of blue early-type galaxies. We propose that while such merging cluster systems absorb extra cosmic web populations hitherto not part of the original merging subclusters, a $\sim$ 1 Gyr timescale is possibly insufficient to result in changes in global colours and morphologies of galaxies.

As comets journey into the inner solar system, they deliver particulates and volatile gases into their comae that reveal the most primitive materials in the solar system. Cometary dust particles provide crucial information for assessing the physico-chemical conditions in the outer disk from which they formed. Compared to the volatiles and soluble organics, the refractory dust particles are more robust and may be traceable to other small bodies. Using data from the Spitzer Heritage Archive, we present thermal dust models of 57 observations of 33 comets observed spectroscopically with the NASA Spitzer Space Telescope. This comet spectral survey offers the opportunity to study comets with data from the same instrument, reduced by the same methods, and fitted by the same thermal model using the same optical constants. The submicron dust tends to be dominated by amorphous carbon, and the submicron silicate mass tends to be dominated by amorphous silicate materials. We discuss the implications of these findings as they relate to Mg-rich crystalline silicates, which are high-temperature condensates, as well as to potential ion irradiation of amorphous Mg:Fe silicates prior to their incorporation into comets. These results impact our understandings of the protoplanetary disk conditions of planetesimal formation. Lastly, we cannot definitively conclude that a distinct difference exists in the dust composition between Oort cloud and Jupiter-family comet dynamical population as a whole.

Isotopic anomalies provide a means of probing the materials responsible for the formation of terrestrial planets. By analyzing new iron isotopic anomaly data from Martian meteorites and drawing insights from published data for O, Ca, Ti, Cr, Fe, Ni, Sr, Zr, Mo, Ru, and Si, we scrutinize potential changes in the isotopic composition of the material accreted by Mars and Earth during their formation. A Principal Component Analysis of isotopic anomalies in meteorites identifies three main clusters (forming the three parts of the isotopic trichotomy): CI, CC=CM+CO+CV+CR, and NC=EH+EL+H+L+LL. Our results suggest that Earth is primarily an isotopic mixture of ~92% E, 6 % CI, and <2% COCV and O. Mars, on the other hand, appears to be a mixture of ~65% E, 33% O, and <2% CI and COCV. We establish that Earth's CI contribution substantially increased during the latter half of its accretion. Mars began accreting a mix of O and E but predominantly accreted E later. Mars' changing isotopic makeup during accretion can be explained if it underwent gas-driven type I migration from its origin near the O - E boundary to a location well within the E region during the first few million years of solar system history. Earth's later increased CI contribution may be attributed to the stochastic impact of an interloper carbonaceous embryo that moved inside the inner solar system region while nebular gas was still present, and subsequently participated in the stage of chaotic growth. The recent findings of Si isotopic anomalies in enstatite chondrites when compared to terrestrial rocks likely stems from insufficient correction for high-temperature equilibrium isotopic fractionation. With appropriate adjustments for this influence, both the silicate Earth and enstatite chondrites exhibit comparable Si isotopic anomalies, reaffirming a genetic link between them.

A wide class of scalar field models including Quintessence and K-essence have the attractive property of tracker regimes, where the energy density stored in the field evolves so as to mimic that of the dominant background component for a period of time. During this evolution, for a brief period of time there is an increase in the energy density of the field as it spirals in towards it's attractor solution. We show that when the peak of this energy density occurs around the epoch of equality, we can address a key requirement of early dark energy (EDE), postulated as a solution to the Hubble tension. In particular we demonstrate how this can occur in a wide class of Quintessence, axion and K-essence models, before showing that the Quintessence models suffer in that they generally lead to sound speeds incompatible with the requirements of EDE, whereas the K-essence and axion models can do a better job of fitting the data.

Galaxy properties primarily depend on their host halo mass. Halo mass, in turn, depends on the cosmic web environment. We explore if the effect of the cosmic web on galaxy properties is entirely transitive via host halo mass, or if the cosmic web has an effect independent of mass. The secondary galaxy bias, sometimes referred to as "galaxy assembly bias", is the beyond-mass component of the galaxy-halo connection. We investigate the link between the cosmic web environment and the secondary galaxy bias in simulations. We measure the secondary galaxy bias through the following summary statistics: projected two-point correlation function, $\wprp$, and counts-in-cylinders statistics, $\Pncic$. First, we examine the extent to which the secondary galaxy bias can be accounted for with a measure of the environment as a secondary halo property. We find that the total secondary galaxy bias preferentially places galaxies in more strongly clustered haloes. In particular, haloes at fixed mass tend to host more galaxies when they are more strongly associated with nodes or filaments. This tendency accounts for a significant portion, but not the entirety, of the total secondary galaxy bias effect. Second, we quantify how the secondary galaxy bias behaves differently depending on the host halo proximity to nodes and filaments. We find that the total secondary galaxy bias is relatively stronger in haloes more associated with nodes or filaments. We emphasise the importance of removing halo mass effects when considering the cosmic web environment as a factor in the galaxy-halo connection.

The Vera C. Rubin Legacy Survey of Space and Time will discover thousands of microlensing events across the Milky Way Galaxy, allowing for the study of populations of exoplanets, stars, and compact objects. It will reach deeper limiting magnitudes over a wider area than any previous survey. We evaluate numerous survey strategies simulated in the Rubin Operation Simulations (OpSims) to assess the discovery and characterization efficiencies of microlensing events. We have implemented three metrics in the Rubin Metric Analysis Framework: a discovery metric and two characterization metrics, where one estimates how well the lightcurve is covered and the other quantifies how precisely event parameters can be determined. We also assess the characterizability of microlensing parallax, critical for detection of free-floating black hole lenses, in a representative bulge and disk field. We find that, given Rubin's baseline cadence, the discovery and characterization efficiency will be higher for longer duration and larger parallax events. Microlensing discovery efficiency is dominated by observing footprint, where more time spent looking at regions of high stellar density including the Galactic bulge, Galactic plane, and Magellanic clouds, leads to higher discovery and characterization rates. However, if the observations are stretched over too wide an area, including low-priority areas of the Galactic plane with fewer stars and higher extinction, event characterization suffers by > 10%, which could impact exoplanet, binary star, and compact object events alike. We find that some rolling strategies (where Rubin focuses on a fraction of the sky in alternating years) in the Galactic bulge can lead to a 15-20% decrease in microlensing parallax characterization, so rolling strategies should be chosen carefully to minimize losses.

We explore the kinematic and chemical properties of Monoceros stellar overdensity by combining data from 2MASS, WISE, APOGEE, and \text{Gaia}. Monoceros is a structure located towards the Galactic anticenter and close to the disk. We identified that its stars have azimuthal velocity in the range of $200 < v_{\phi}\,{\rm(km\,s^{-1})}< 250$. Combining their kinematics and spatial distribution, we designed a new method to select stars from this overdensity. This method allows us to easily identify the structure in both hemispheres and estimate their distances. Our analysis was supported by comparison with simulated data from the entire sky generated by $\texttt{Galaxia}$ code. Furthermore, we characterized, for the first time, the Monoceros overdensity in several chemical-abundance spaces. Our results confirm its similarity to stars found in the thin disk of the Galaxy and suggest an \textit{in situ} formation. Furthermore, we demonstrate that the southern (Mon-S) and northern (Mon-N) regions of Monoceros exhibit indistinguishable chemical compositions.

Previous work by our group has shown the potential for an aerogravity assist at Triton using a trailing ballute to capture a probe into orbit about Neptune. The current work extends that study by using the LOFTID aeroshell configuration, a flight-proven inflatable, to perform the aeromaneuver. Numerical simulations were carried out beginning at the atmospheric interface for a range of entry velocities (3.0-14.0 km/s) and angles (45-56 degrees). The spin-stabilized vehicle was assumed to fly at zero angle of attack, producing no lift, and is capable of reaching outbound conditions sufficient to capture into orbit about Neptune without violating the vehicle's aerothermal limits or penetrating too deeply into Triton's atmosphere.

The nature of the first seeds of supermassive black holes (SMBHs) is currently unknown, with postulated initial masses ranging from $\sim10^5~M_{\odot}$ to as low as $\sim10^2~M_{\odot}$. However, most existing cosmological simulations resolve BHs only down to $\sim10^5-10^6~M_{\odot}$. In this work, we introduce a novel sub-grid BH seed model that is directly calibrated from high resolution zoom simulations that can trace the formation and growth of $\sim 10^3~M_{\odot}$ seeds forming in halos with pristine, star-forming gas. We trace the BH growth along merger trees until their descendants reach masses of $\sim10^4$ or $10^5~M_{\odot}$. The descendants assemble in galaxies with a broad range of properties (e.g., halo masses $\sim10^7-10^9~M_{\odot}$) that evolve with redshift and are sensitive to seed parameters. The results are used to build a new stochastic seeding model that directly seeds these descendants in lower resolution versions of our zoom region. Remarkably, we find that by seeding the descendants simply based on total galaxy mass, redshift and an environmental richness parameter, we can reproduce the results of the detailed gas based seeding model. The baryonic properties of the host galaxies are well reproduced by the mass-based seeding criterion. The redshift-dependence of the mass-based criterion captures the influence of halo growth, star formation and metal enrichment on seed formation. The environment based seeding criterion seeds the descendants in rich environments with higher numbers of neighboring galaxies. This accounts for the impact of unresolved merger dominated growth of BHs, which produces faster growth of descendants in richer environments with more extensive BH merger history. Our new seed model will be useful for representing a variety of low mass seeding channels within next generation larger volume uniform cosmological simulations.

The Asteroid Terrestrial-impact Last Alert System (ATLAS) observes the visible sky every night in search of dangerous asteroids. With four soon five) sites ATLAS is facing new challenges for scheduling observations and linking detections to identify moving asteroids. Flexibility in coping with diverse observation sites and times of detections that can be linked is critical, as is optimization of observing time for coverage versus depth. We present new algorithms to fit orbits rapidly to sky-plane observations, and to test and link sets of detections to find the ones which belong to moving objects. The PUMA algorithm for fitting orbits to angular positions on the sky executes in about a millisecond, orders of magnitude faster than the methods currently in use by the community, without sacrifice in accuracy. The PUMALINK algorithm to find linkages among sets of detections has similarities to other approaches, notably HelioLinC, but it functions well at asteroid ranges of a small fraction of an AU. Candidate linkages are checked by the PUMA library to test that the detections correspond to a real orbit, even at close range, and the false alarm rate is manageable. We present the results of tests of PUMALINK on three datasets which illustrate PUMALINK's effectiveness and economy: 2 weeks of all ATLAS detections over the sky, 2 weeks of special ATLAS opposition observations with long exposure time, and 2 weeks of simulated LSST asteroid observations. Testing only pairs of nights PUMALINK achieves approximately 90% detection probability for real objects while keeping the false alarm rate below 10%. Both numbers improve greatly when PUMALINK is given a third night.

We consider the dynamics of an equatorial explosion powered by a millisecond magnetar formed from the core collapse of a massive star. We study whether these outflows -- generated by a priori magneto-centrifugally-driven, relativistic magnetar winds -- might be powerful enough to produce an ultra-relativistic blade ("lamina") that successfully carves its way through the dense stellar interior. We present high-resolution numerical special-relativistic hydrodynamic simulations of axisymmetric centrifugally-driven explosions inside a star and follow the blast wave propagation just after breakout. We estimate the engine requirements to produce ultra-relativistic lamina jets and comment on the physicality of the parameters considered. We find that sufficiently collimated -- half-opening angle $\theta_r \leq 0.2^\circ$ -- laminas successfully break out of a compact progenitor at ultra-relativistic velocities ($\Gamma_{\rm core} \gtrsim 30$) and extreme isotropic energies ($E_{k,\rm iso} \sim 5 \times 10^{52}\text{erg}$) within a few percent of the typical spin-down period for a millisecond magnetar. The various phases of these ultra-thin outflows such as collimation shocks, Kelvin-Helmholtz instabilities, and lifetime are discussed and we speculate on the observational signatures echoed by this outflow geometry.

Search for extraterrestrial intelligence (SETI) has been mainly focused on nearby stars and their planets in recent years. Barnard's star is the second closest star system to the sun and the closest star in the FAST observable sky which makes the minimum Equivalent Isotropic Radiated Power (EIRP) required for a hypothetical radio transmitter from Barnard's star to be detected by FAST telescope a mere 4.36x10^8 W. In this paper, we present the Five-hundred-meter Aperture Spherical radio Telescope (FAST) telescope as the most sensitive instrument for radio SETI observations toward nearby star systems and conduct a series of observations to Barnard's star (GJ 699). By applying the multi-beam coincidence matching (MBCM) strategy on the FAST telescope, we search for narrow-band signals (~Hz) in the frequency range of 1.05-1.45 GHz, and two orthogonal linear polarization directions are recorded. Despite finding no evidence of radio technosignatures in our series of observations, we have developed predictions regarding the hypothetical extraterrestrial intelligence (ETI) signal originating from Barnard's star. These predictions are based on the star's physical properties and our observation strategy.

Although it is generally believed that the solar photosphere is not magnetically force-free owning to its high plasma $\beta$, the estimations of force-freeness using observed magnetograms have produced disputable results. Some studies confirmed that the photosphere is largely not force-free whereas some authors argued that the photosphere is not far away from being force-free. In a previous paper of ours we demonstrated that, due to the fact that the noise levels of the transverse field in the magnetograms are much larger than those of the vertical field, wrong judgements on the force-freeness could be made: a truly force-free field could be judged as being not-force-free and a truly not-force-free field could be judged as being force-free. Here in this letter we propose an approach to overcome this serious problem. By reducing the spatial resolution to lower the noise level, the heavy influence of the measurement noise on the force-freeness judgement can be significantly suppressed. We first use two analytical solutions to show the success and the effectiveness of this approach. Then we apply this new approach to two large datasets of active region magnetograms, obtained with the HMI/SDO and SP/Hinode, respectively. Our analysis shows that the photospheric magnetic fields are actually far away from being force-free. Particularly and most notably, the mean value of Fz/Fp (where Fz the net Lorentz force in the vertical direction and Fp the total Lorentz force) is as low as -0.47, with more than 98% of the active regions having |Fz/Fp|>0.1, when using the SP/Hinode magnetograms of true field strength.

The mass ($M_{\rm TOV}$) and radius ($R_{\rm TOV}$) of the maximum-mass nonrotating neutron star (NS) play a crucial role in constraining the elusive equation of state (EOS) of cold dense matter and in predicting the fate of remnants from binary neutron star (BNS) mergers. In this study, we introduce a novel method to deduce these parameters by examining the mergers of second-generation (2G) BHs with NSs. These 2G BHs are assumed to originate from supramassive neutron stars (SMNSs) formed in BNS mergers. Since the properties of the remnant BHs arising from the collapse of SMNSs follow a universal relation governed by $M_{\rm TOV}$ and $R_{\rm TOV}$, we anticipate that by analyzing a series ($\sim 100$ detections) of mass and spin measurements of the 2G BHs using the third-generation ground-based gravitational wave detectors, $M_{\rm TOV}$ and $R_{\rm TOV}$ can be determined with a precision of $\sim 0.01M_\odot$ and $\sim 0.6$ km, respectively.

We report the measurement of rest-frame UV size and morphology of H$\alpha$-emission-selected star-forming galaxies (HAEs) in four protoclusters at z $\sim$ 2 (PKS 1138-262, USS 1558-003, PHz G237.0+42.5, and CC 2.2) using archival Hubble Space Telescope Advanced Camera Surveys (HST/ACS) F814W data. We compare the measurement of 122 HAEs in protoclusters detected by HST/ACS to a coeval comparison field sample of 436 HAEs. We find the size distributions of protocluster and field HAEs are similar with typical half-light radius of $\sim$ 2.5 kpc. At fixed stellar mass, there is no significant difference between HAE in protocluster and in field, which is also supported by stacking analyses. This result suggests that the environment does not significantly affect the size of galaxies during the star-forming phase at this epoch. Based on S\'ersic index and non-parametric morphologies, HAE morphologies in both environments at $z\sim2$ in rest-frame UV are consistent with disk-like star-forming galaxies, although we also find $29\% \pm 4\%$ HAEs have peculiar or disturbed morphologies. The fraction of disturbed galaxies is higher in protocluster environment, with $39 \pm 8 \%$ protocluster HAEs showing disturbed morphologies, compared to $26\pm4\%$ in the comparison field. The apparent disturbed morphologies are correlated with higher star-formation activity and may be caused by either in situ giant clumps or mergers.

Context. The Doppler shift predicted by general relativity for light escaping a gravitational potential has been observed on Earth as well as in the direction of various stars and galaxy clusters at optical wavelengths. Aims. Observing the gravitational redshift in the X-ray band within galaxy clusters could provide information on their properties and, in particular, their gravitational potential. We present a feasibility study of such a measurement, using the capabilities of the next-generation European X-ray observatory Athena. Methods. We used a simple generalized Navarro-Frenk-White potential model along with a beta-model for the density of baryonic matter, which sets the emission to provide an estimation of the observed redshift in the simplest of cases. We generated mock observations with the Athena X-ray Integral Field Unit (X-IFU) for a nearby massive cluster, while seeking to recover the gravitational redshift along with other properties of the toy model cluster. Results. We investigated the observability of the gravitational redshift in an idealized test case of a nearby massive cluster with the Athena X-IFU instrument, as well as its use in probing the properties of the potential well. We were also able to constrain the mass to a 20 % level of precision and the cosmological redshift to less than 1%, within a simplified and idealized observational framework. More refined simulations accounting for further effects such as the internal gas motions and the actual shape of the potential well are required to fully investigate the feasibility of measuring the gravitational redshift for a single target or statistically over a sample of galaxy clusters.

The sources of the astrophysical flux of high-energy neutrinos detected by IceCube are still largely unknown, but searches for temporal and spatial correlation between neutrinos and electromagnetic radiation are a promising approach in this endeavor. All major imaging atmospheric Cherenkov telescopes (IACTs) - FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program of target-of-opportunity observations of neutrino alerts issued by IceCube. These programs use several complementary neutrino alert streams. A publicly distributed alert stream is formed by individual high-energy neutrino candidate events of potentially astrophysical origin, such as IceCube-170922A (which could be linked to the flaring blazar TXS\,0506+056). A privately distributed alert stream is formed by clusters of neutrino events in time and space around either pre-selected gamma-ray sources or anywhere in the sky. Here, we present joint searches for multi-wavelength emission associated with a set of IceCube alerts, both private and public, received through mid-January 2021. We will give an overview of the programs of the participating IACTs. We will showcase the various follow-up and data analysis strategies employed in response to the different alert types and various possible counterpart scenarios. Finally, we will present results from a combined analysis of the VHE gamma-ray observations obtained across all involved instruments, as well as relevant multi-wavelength data.

Context. In the scope of space weather forecasting, it is crucial to be able to more reliably predict the arrival time, speed, and magnetic field configuration of coronal mass ejections (CMEs). From the time a CME is launched, the dominant factor influencing all of the above is the interaction of the interplanetary CME (ICME) with the ambient plasma and interplanetary magnetic field. Aims. Due to a generally anisotropic heliosphere, differently oriented ICMEs may interact differently with the ambient plasma and interplanetary magnetic field, even when the initial eruption conditions are similar. For this, we examined the possible link between the orientation of an ICME and its propagation in the heliosphere (up to 1 AU). Methods. We investigated 31 CME-ICME associations in the period from 1997 to 2018. The CME orientation in the near-Sun environment was determined using an ellipse-fitting technique applied to single-spacecraft data from SOHO/LASCO C2 and C3 coronagraphs. In the near-Earth environment, we obtained the orientation of the corresponding ICME using in situ plasma and magnetic field data. The shock orientation and nonradial flows in the sheath region for differently oriented ICMEs were investigated. In addition, we calculated the ICME transit time to Earth and drag parameter to probe the overall drag force for differently oriented ICMEs. The drag parameter was calculated using the reverse modeling procedure with the drag-based model. Results. We found a significant difference in nonradial flows for differently oriented ICMEs, whereas a significant difference in drag for differently oriented ICMEs was not found.

Every population of small bodies in the Solar system contains a sizable fraction of multiple systems. Among these, the Jupiter Trojans have the lowest number of known binary systems and the least characterized. We aim at characterizing the reported binary system (17365) Thymbraeus, one of the only seven multiple systems known among Jupiter Trojans. We conducted light curves observing campaigns in 2013, 2015, and 2021 with ground-based telescopes. We model these lightcurves using dumbbell equilibrium figures. We show that Thymbraeus is unlikely a binary system. Its light curves are fully consistent with a bilobated shape: a dumbbell equilibrium figure. We determine a low density of 830 +/- 50 kg.m-3 , consistent with the reported density of other Jupiter Trojan asteroids and small Kuiper-belt objects. The angular velocity of Thymbraeus is close to fission. If separated, its components would become a similarly-sized double asteroid such as the other Jupiter Trojan (617) Patroclus.

The dynamics of accreting and outgoing flows around compact objects depends crucially on the strengths and configurations of the magnetic fields therein, especially of the large-scale fields that remain coherent beyond turbulence scales. Possible origins of these large-scale magnetic fields include flux advection and disc dynamo actions. However, most numerical simulations have to adopt an initially strong large-scale field rather than allow them to be self-consistently advected or amplified, due to limited computational resources. The situation can be partially cured by using sub-grid models where dynamo actions only reachable at high resolutions are mimicked by artificial terms in low-resolution simulations. In this work, we couple thin-disc models with local shearing-box simulation results to facilitate more realistic sub-grid dynamo implementations. For helical dynamos, detailed spatial profiles of dynamo drivers inferred from local simulations are used, and the nonlinear quenching and saturation is constrained by magnetic helicity evolution. In the inner disc region, saturated fields have dipole configurations and can reach $\beta\simeq 0.1$ to $100$, with correlation lengths $\simeq h$ in the vertical direction and $\simeq 10h$ in the radial direction, where $h$ is the disc scale height. The dynamo cycle period is $\simeq 40$ orbital time scale, compatible with previous global simulations. Additionally, we explore two dynamo mechanisms which do not require a net kinetic helicity and have only been studied in shearing-box setups. We show that such dynamos are possible in thin accretion discs, but produce field configurations that are incompatible with previous results. We discuss implications for future general-relativistic magnetohydrodynamics simulations.

The atmosphere of a terrestrial planet that is replenished with secondary gases should have accumulated hydrogen-rich gas from its protoplanetary disk. Although a giant impact blows off a large fraction of the primordial atmosphere of a terrestrial planet in the late formation stage, the remaining atmosphere can become water-rich via chemical reactions between hydrogen and vaporized core material. We find that a water-rich post-impact atmosphere forms when a basaltic or CI chondrite core is assumed. In contrast, little post-impact water is generated for an enstatite chondrite core. We investigate the X-ray- and UV-driven mass loss from an Earth-mass planet with an impact-induced multi-component H$_2$$-$He$-$H$_2$O atmosphere for Gyrs. We show that water is left in the atmosphere of an Earth-mass planet when the low flux of escaping hydrogen cannot drag water upward via collisions. For a water-dominated atmosphere to form, the atmospheric mass fraction of an Earth-mass planet with an oxidizing core after a giant impact must be less than a few times 0.1%. We also find that Earth-mass planets with water-dominated atmospheres can exist at semimajor axes ranging from a few times 0.1 au to a few au around a Sun-like star depending on the mass loss efficiency. Such planets are important targets for atmospheric characterization in the era of JWST. Our results indicate that efficient mixing between hydrogen and rocky components during giant impacts can play a role in the production of water in an Earth-mass planet.

Coherent radiation of pulsars, magnetars, and fast radio bursts could, in theory, be interpreted as radiation from solitons and soliton-like waves. The solitons are meant to contain a large number of electric charges confined on long time-scales and may radiate strongly by coherent curvature emission. However, solitons are also known to undergo a wave collapse, which may cast doubts on the correctness of the soliton radio emission models of neutron stars. We investigate the evolution of the caviton type of solitons self-consistently formed by the relativistic streaming instability and compare their apparent stability in 1D calculations with more generic 2D cases, in which the solitons are seen to collapse. Three representative cases of beam Lorentz factors and plasma temperatures are studied to obtain soliton dispersion properties. We utilized 1D electrostatic and 2D electromagnetic relativistic particle-in-cell simulations at kinetic microscales. We found that no solitons are generated by the streaming instability in the 2D simulations. Only superluminal L-mode (relativistic Langmuir) waves are produced during the saturation of the instability, but these waves have smaller amplitudes than the waves in the 1D simulations. The amplitudes tend to decrease after the instability has saturated, and only waves close to the light line, $\omega=c k$, remain. Solitons in the 1D approach are stable for $\gamma_\mathrm{b}\gtrsim60$ but disappear for low beam Lorentz factor $\gamma_\mathrm{b}<6$. Our examples show that the superluminal soliton branch that is formed in 1D simulations will not be generated by the relativistic streaming instability when more dimensional degrees of freedom are present. The soliton model can, therefore, not be used to explain the coherent radiation of pulsars, magnetars, and fast radio bursts - unless there are alternative generation mechanisms.

The tidal disruption event AT2018hyz was a regular optically detected one with no special prompt features. However, it suddenly displayed a fast-rising radio flare almost three years after the disruption. The flare is most naturally interpreted as arising from an off-axis relativistic jet. We didn't see the jet at early times as its emission was relativistically beamed away from us. However, we could see the radiation once the jet has slowed down due to interaction with the surrounding matter. Analysis of the radio data enabled estimates of the jet's kinetic energy and opening angle, as well as the conditions (size and magnetic field) within the radio-emitting region. We show here that such a jet satisfies the Hillas condition for the acceleration of UHECRs to the highest energies. We also show that the rate and total power of this event are consistent with the observed luminosity density of UHECRs. These results strongly support earlier suggestions that TDEs are the sources of UHECRs.

V570 Per is a binary star system containing two F-type stars in a 1.90 d period circular orbit. It shows shallow partial eclipses that were discovered from its Hipparcos light curve. We present an analysis of this system based on two sectors of high-quality photometry from the NASA Transiting Exoplanet Survey Satellite (TESS) mission, and published spectroscopic light ratio and radial velocity measurements. We find masses of 1.449 +/- 0.006 and 1.350 +/- 0.006 Msun, and radii of 1.538 +/- 0.035 and 1.349 +/- 0.032 Rsun. The radius measurements are set by the spectroscopic light ratio and could be improved by obtaining a more precise light ratio. The eclipses in the TESS data arrived 660 +/- 30 s later than expected, suggesting the presence of a faint third body on a wider orbit around the eclipsing system. Small trends in the residuals of the fit to the TESS light curve are attributed to weak starspots. The distance to the system is close to the Gaia DR3 value, but the Gaia spectroscopic orbit is in moderate disagreement with the results from the published ground-based data.

The Roche potential is the sum of the gravitational and rotational potentials experienced by a massless body rotating alongside two massive bodies in a circular orbit. The Lagrangian points are five stationary points in the Roche potential. The positions of two of the Lagrangian points (L4 and L5) are fixed. The other three (L1, L2 and L3) are along the line joining the two masses: their positions depend on the mass ratio, $q$, and can be calculated numerically by finding the roots of a quintic polynomial. Analytical approximations to their positions are useful in several situations, but existing ones are designed for small mass ratios. We present new approximations valid for all mass ratios from zero to unity: \begin{eqnarray*} x_{\rm L1} & = & 1 - \frac{q^{0.33071}}{0.51233\,q^{0.49128} + 1.487864} \\ x_{\rm L2} & = & 1 + \frac{q^{0.8383} + 2.891\,q^{0.3358}}{1.525\,q^{0.848} + 4.046596} \\ x_{\rm L3} & = & -1 + \frac{q^{1.007}}{1.653\,q^{0.9375} + 1.66308} \end{eqnarray*} in a rotating frame of reference where the more massive body is at $x=0$ and the less massive body at $x=1$. The three approximations are precise to $6 \times 10^{-5}$ for all mass ratios.

We make use of the deepest JWST/MIRI image at 5.6 um, obtained in the Hubble eXtreme Deep Field (XDF), to constrain the role of strong Ha emitters (HAEs) in Cosmic Reionization at z~7-8. Our sample of bright (M(UV) < -20 mag) HAEs is comprised of young (<30 Myr) galaxies with low stellar masses (<= 10^9 Msun). They span a wide range of UV-beta slopes, with a median beta = -2.22+-0.35, which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency (xi_ion,0) of these sources (assuming f_esc,LyC = 0), which yields a median value log10(xi_ion,0/(Hz erg^(-1))) = 25.54(+0.09, -0.10). We show that xi_ion,0 positively correlates with EW0(Ha) and specific star formation rate (sSFR). Instead xi_ion,0 weakly anti-correlates with stellar mass and beta. Based on the beta values, we estimate f_esc,LyC=0.07(+0.03, -0.02), which results in log10(xi_ion/(Hz erg^(-1))) = 25.59 (+0.06, -0.04). By considering this result along with others from the literature, we find a mild evolution of xi_ion with redshift. Finally, we assess the impact of strong HAEs during Cosmic Reionization at z~7-8. We find that our HAEs do not need high values of f_esc, rel (only 6-10%) to be able to reionize their surrounding intergalactic medium. They have N_dot_ion = 10^(50.43+-0.3) s^(-1)Mpc^(-3) and contribute more than a factor of two in terms of emitted ionizing photons per comoving volume compared to non-Ha emitters in the same redshift bin, suggesting that strong, young, and low stellar-mass emitters could have played a central role during the Epoch of Reionization.

Intended and unintended radio emissions from satellites can interfere with sensitive radio telescopes in the frequency ranges of key experiments in astrophysics and cosmology. We detect strong intended and unintended electromagnetic radiation from Starlink satellites at the site of the future SKA-Low facility in Western Australia, using an SKA-low prototype station known as the Engineering Development Array version 2 (EDA2). We aim to show that Starlink satellites are easily detectable utilising a configuration of low frequency radio antennas representative of an SKA-Low 'station' and that our results complement similar findings with the LOFAR telescope. Utilising the EDA2 at frequencies of 137.5 MHz and 159.4 MHz, we detect trains of Starlink satellites on 2023-03-17/18 and 2021-11-16/17, respectively, via the formation of all-sky images with a frequency resolution of 0.926 MHz and a time resolution of 2 s. Time differencing techniques are utilised to isolate and characterise the transmissions from Starlink and other satellites. We observe Starlink satellites reaching intensities of $10^6$ Jy/beam, with the detected transmissions exhibiting a range of behaviours, from periodic bursts to steady transmission. The results are notable because they demonstrate that Starlink satellites are detected in the SKA-Low frequency range, transmitting both intentionally and unintentionally. Follow-up work and discussion are needed to identify the cause of this unintentional radiation as it has the potential to interfere with SKA-Low science. Our results indicate that both intended and unintended radiation from Starlink satellites will be detrimental to key SKA science goals without mitigation. Continued conversation with SpaceX could potentially result in future mitigations which the EDA2 instrument could efficiently monitor and characterise at the SKA-Low site.

We present the first deep X-ray observations of a luminous FBOT AT2018cow at $\sim 3.7\,\rm{yr}$ since discovery. These observations revealed the presence of a luminous X-ray source with $L_{\rm x}\approx 4\times 10^{38}\,\rm{erg\,s^{-1}}$ at the location of AT2018cow. The very soft X-ray spectrum and sustained luminosity are clearly distinct from the spectral and temporal behavior of AT2018cow in the first $\sim100$ days of evolution, and signal the emergence of a new emission component at late times. We interpret these findings in the context of the late-time panchromatic emission from AT2018cow, which includes the detection of persistent, slowly-fading UV emission with $\nu L_{\nu}\approx 10^{39}\,\rm{erg\,s^{-1}}$. Similar to previous works, (and in strict analogy with arguments used for Ultra-Luminous X-ray sources --ULXs), we find that these late-time observations are consistent with thin-disks around Intermediate Mass Black Holes (IMBHs, with $M_{BH}\approx 10^3-10^4\,\rm{M_{\odot}}$) accreting at sub-Eddington rates. However, differently from previous studies, we find that smaller-mass BHs with $M_{BH}\approx 10-100\,\rm{M_{\odot}}$ accreting at $\gtrsim$ the Eddington rate cannot be ruled out, and in fact provide a natural explanation for the inferred compact size ($R_{\rm out}\approx 40\,R_{\odot}$) of the accretion disk years after the optical flare. Most importantly, irrespective of the accretor mass, our study lends support to the hypothesis that LFBOTs are accretion-powered phenomena and that, specifically, LFBOTs constitute electromagnetic manifestations of super-Eddington accreting systems that evolve to $\lesssim$ Eddington over a $\approx 100$\,days time scale.

Gravity modes (g modes), mixed gravito-acoustic modes (mixed modes), and gravito-inertial modes (gi modes) possess unmatched properties as probes for stars with radiative interiors. The structural and dynamical constraints that they are able to provide cannot be accessed by other means. While they provide precious insights into the internal dynamics of evolved stars as well as massive and intermediate-mass stars, their non-detection in main sequence (MS) solar-type stars make them a crucial missing piece in our understanding of angular momentum transport in radiative zones and stellar rotational evolution. In this work, we aim to apply certain analysis tools originally developed for helioseismology in order to look for g-mode signatures in MS solar-type stars. We select a sample of the 34 most promising MS solar-type stars with Kepler four-year long photometric time series. All these stars are well-characterised late F-type stars with thin convective envelopes, fast convective flows, and stochastically excited acoustic modes (p modes). For each star, we compute the background noise level of the Fourier power spectrum to identify significant peaks at low frequency. After successfully detecting individual peaks in 12 targets, we further analyse four of them and observe distinct patterns of surrounding peaks with a low probability of being noise artifacts. Comparisons with the predictions from reference models suggest that these patterns are compatible with the presence of non-asymptotic low-order pure g modes, pure p modes, and mixed modes. Given their sensitivity to both the convective core interface stratification and the coupling between p- and g-mode resonant cavities, such modes are able to provide strong constraints on the structure and evolutionary states of the related targets. [abridged]

In this work, we study the dynamics of two less massive objects moving around a central massive object, which are all embedded within a thin accretion disc. In addition to the gravitational interaction between these objects, the disc-object interaction is also crucial for describing the long-term dynamics of the multi-body system, especially in the regime of mean-motion resonances. We point out that near the resonance the density waves generated by the two moving objects generally coherently interfere with each other, giving rise to extra angular momentum fluxes. The resulting backreaction on the objects is derived within the thin-disc scenario, which explicitly depends on the resonant angle. With this density-wave mediated interaction included, we find that a system initially locked into the mean-motion resonance either asymptotes to a quai-stationary fixed point or automatically exits the resonance with large amplitude circulations. We have performed hydrodynamical simulations with planets embedded within a thin accretion disc and have found signatures of interfering density waves from the evolution of planet eccentricities. By including the type-I migration torques in the evolution of a pair of planets, we show that the eccentricity-damping effect contributed by the interfering densities may increase the period ratio of the planets when they are trapped in mean-motion resonances. This may explain the $1\%-2\%$ offset (for the period ratios) from the exact resonance values as observed in Kepler multi-planet systems.

Recently, cosmography emerged as a valuable tool to effectively describe the vast amount of astrophysical observations without relying on a specific cosmological model. Its model-independent nature ensures a faithful representation of data, free from theoretical biases. Indeed, the commonly assumed fiducial model, the $\Lambda$CDM, shows some shortcomings and tensions between data at late and early times that need to be further investigated. In this paper, we explore an extension of the standard cosmological model by adopting the $f(z)$CDM approach, where $f(z)$ represents the cosmographic series characterizing the evolution of recent universe driven by dark energy. To construct $f(z)$, we take into account the Pad\'e series, since this rational polynomial approximation offers a better convergence at high redshifts than the standard Taylor series expansion. Several orders of such an approximant have been proposed in previous works, here we want to answer the questions: What is the impact of the cosmographic series choice on the parameter constraints? Which series is the best for the analysis? So, we analyse the most promising ones by identifying which order is preferred in terms of stability and goodness of fit. Theoretical predictions of the $f(z)$CDM model are obtained by the Boltzmann solver code and the posterior distributions of the cosmological and cosmographic parameters are constrained by a Monte Carlo Markov Chains analysis. We consider a joint data set of cosmic microwave background temperature measurements from the Planck collaboration, type Ia supernovae data from the latest Pantheon+ sample, baryonic acoustic oscillations and cosmic chronometers data. In conclusions, we state which series can be used when only late time data are used, while which orders has to be considered in order to achieve the necessary stability when large redshifts are considered.

Star formation histories (SFH) of early (6$<z<$12) galaxies have been found to be highly stochastic in both simulations and observations, while at $z\lesssim$6 the presence of a main sequence (MS) of star-forming galaxies imply secular processes at play. In this work, we aim at characterising the SFH variability of early galaxies as a function of their stellar mass and redshift. We use the JADES public catalogue and derive the physical properties of the galaxies as well as their SFH using the spectral energy distribution modelling code CIGALE. To this aim, we implement a non-parametric SFH with a flat prior allowing for as much stochasticity as possible. We use the SFR gradient, an indicator of the movement of galaxies on the SFR-$M_\ast$ plane, linked to the recent SFH of galaxies. This dynamical approach of the relation between the SFR and stellar mass allows us to show that, at $z>9$, 87% of massive galaxies, ($\log(M_\ast/M_\odot)\gtrsim$9), have SFR gradients consistent with a stochastic star-formation activity during the last 100 Myr, while this fraction drops to 15% at $z<7$. On the other hand, we see an increasing fraction of galaxies with a star-formation activity following a common stream on the SFR-$M_\ast$ plane with cosmic time, indicating that a secular mode of star-formation is emerging. We place our results in the context of the observed excess of UV emission as probed by the UV luminosity function at $z\gtrsim10$, by estimating $\sigma_{UV}$, the dispersion of the UV absolute magnitude distribution, to be of the order of 1.2mag and compare it with predictions from the literature. In conclusion, we find a transition of star-formation mode happening around $z\sim9$: Galaxies with stochastic SFHs dominates at $z\gtrsim9$, although this level of stochasticity is too low to reach those invoked by recent models to reproduce the observed UV luminosity function.

Line-intensity mapping (LIM) experiments coming online now will survey fluctuations in aggregate emission in the [C II] ionized carbon line from galaxies at the end of reionization. Experimental progress must be matched by theoretical reassessments of approaches to modelling and the information content of the signal. We present a new model for the halo-[C II] connection, building upon results from the FIRE simulations suggesting that gas mass and metallicity most directly determine [C II] luminosity. Applying our new model to an ensemble of peak-patch halo lightcones, we generate new predictions for the [C II] LIM signal at $z\gtrsim6$. We expect a baseline 4000-hour LIM survey from the CCAT facility to have the fundamental sensitivity to detect the [C II] power spectrum at a significance of $4\sigma$ at $z\sim6$, with an extended or successor Stage 2 experiment improving significance to $36\sigma$ at $z\sim6$ and achieving $8\sigma$ at $z\sim7.5$. Cross-correlation through stacking, simulated against a mock narrow-band Lyman-break galaxy survey, would yield a strong detection of the radial profile of cosmological [C II] emission surrounding star-forming galaxies. We also analyse the role of a few of our model's parameters through the pointwise relative entropy (PRE) of the distribution of [C II] intensities. While the PRE signature of different model parameters can become degenerate or diminished after factoring in observational distortions, various parameters do imprint themselves differently on the one-point statistics of the intrinsic signal. Further work can pave the way to access this information and distinguish different sources of non-Gaussianity in the [C II] LIM observation.

The temperature of the solar atmosphere increases from thousands to millions of degrees moving from the lower layer, the chromosphere, to the outermost one, the corona, while density drops by several orders of magnitude. Such a phenomenon is called temperature inversion and how it happens is still largely unknown. We argue that temperature fluctuations in the chromosphere play a key role, as suggested by the study of a kinetic model of a plasma confined in a semicircular tube subjected to the gravity of the Sun and in contact with a thermostat at its feet, mimicking a coronal loop anchored in the chromosphere. Collisions are neglected in the corona, with a sharp transition to a fully collisional chromosphere. Numerical simulations and analytical calculations show that suitable fluctuations of the thermostat temperature drive the plasma towards a non-thermal stationary state with temperature and density profiles strikingly similar to those observed in the atmosphere of the Sun, suggesting this mechanism may significantly contribute to coronal heating.

The cosmic acceleration problem remains one of the most significant challenges in cosmology. One of the proposed solutions to this problem is the modification of gravity on large scales. In this paper, we explore the well-known $\mu$ and $\Sigma$ parametrization scenarios, and confront them with observational data, including the cosmic microwave background (CMB) radiation from the Wilkinson Microwave Anisotropy Probe (WMAP), Atacama Cosmology Telescope (ACT), and South Pole Telescope (SPT), as well as large-scale structure data from the Sloan Digital Sky Survey (SDSS: BAO+RSD) and Pantheon Supernovae (SN) catalog. We employ a Bayesian framework to constrain the model parameters and discuss the implications of our results on the viability of modified gravity theories. Our analysis reveals the strengths and limitations of the $\mu$ and $\Sigma$ parametrization and provides valuable insights into the nature of gravity on cosmological scales. From the joint analysis ACT + WMAP + SDDS + SN, we find $\mu -1 = 0.02 \pm 0.19$ and $\Sigma -1 = 0.021 \pm 0.068$ at 68\% CL. In light of the SPT + WMAP + SDDS + SN, we find $\mu -1 = 0.07 \pm 0.18$ and $\Sigma -1 = -0.009^{+0.078}_{-0.11}$ at 68\% CL. We also discuss and present several other results involving these datasets. In all the analysis carried out, we did not find any deviations from the theory of general relativity. Our results represent an observational update on the well-known $\mu$-$\Sigma$ parameterization in view of current CMB data, independent and competitive with the constraints obtained with the Planck data.

Next-generation gravitational-wave detectors will provide unprecedented sensitivity to inspiraling binary neutron stars and black holes, enabling detections at the peak of star formation and beyond. However, the signals from these systems will last much longer than those in current detectors, and overlap in both time and frequency, leading to increased computational cost to search for them with standard matched filtering analyses, and a higher probability that they are observed in the presence of non-Gaussian noise. We therefore present a method to search for gravitational waves from compact binary inspirals in next-generation detectors that is computationally efficient and robust against gaps in data collection and noise non-stationarities. Our method, based on the Hough Transform, finds tracks in the time/frequency plane of the detector that uniquely describe specific inspiraling systems. We find that we could detect $\sim 5$ overlapping, intermediate-strength signals without a sensitivity loss. Additionally, we demonstrate that our method can enable multi-messenger astronomy: using only low frequencies ($2-20$ Hz), we could warn astronomers $\sim 2.5$ hours before a GW170817-like merger at 40 Mpc and provide a sky localization of $\sim 20$ deg$^2$ using only one interferometer. Additionally, assuming that primordial black holes (PBHs) exist, we derive projected constraints on the fraction of dark matter they could compose, $f_{\rm PBH}\sim 10^{-6}-10^{-4}$, for $\sim 1-0.1M_\odot$ equal-mass systems, respectively, using a rate suppression factor $f_{\rm sup}=2.5\times 10^{-3}$. Our method only incurs a strain sensitivity loss of a factor of a few at binary neutron star masses compared to the matched filter, which may be reduced depending on available computational power.

We present constraints on the amplitude of local Primordial Non-Gaussianities (PNG), $f_{\rm NL}$, using the quasar sample in the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 16. We analyze the power spectrum monopole, testing for the presence of scale dependent galaxy bias induced by local PNG. Our analysis makes use of optimal redshift weights that maximize the response of the quasar sample to the possible presence of non zero PNG. We find $-4<f_{\rm NL}<27$ at $68\%$ confidence level, which is among the strongest bounds with Large Scale Structure data. The optimal analysis reduces the error bar by about 10% compared to the standard one, an amount which is smaller than what was expected from a Fisher matrix analysis. This, and the reduced improvement over previous releases of the same catalog, suggest the presence of still unknown systematic effects in the data. If the quasars have a lower response to local PNG, our optimal constraint becomes $-23<f_{\rm NL}<21$ at $68\%$, with an improvement of $30\%$ over standard analyses. We also show how to use the optimal weights to put data-driven priors on the sample's response to local PNG.

The no-scale flipped SU(5) superstring framework constitutes a very promising paradigm for physics below the Planck scale providing us with a very rich cosmological phenomenology in accordance with observations. In particular, it can accommodate Starobinsky-like inflation, followed by a reheating phase, which is driven by a light "flaton" field, and during which the GUT phase transition occurs. In this Letter, we extract for the first time a gravitational-wave (GW) signal which naturally arises in the context of the flipped SU(5) cosmological phenomenology and is related to the existence of an early matter era (eMD) driven by the flaton field. Specifically, we study GWs non-linearly induced by inflationary perturbations and which are abundantly produced during a sudden transition from the flaton-driven eMD era to the late-time radiation-dominated era. Remarkably, we find a GW signal with a characteristic peak frequency $f_\mathrm{GW,peak}$ depending only on the string slope $\alpha'$ and reading as $f_\mathrm{GW,peak} \propto 10^{-9} \left(\frac{\alpha'}{\alpha'_*}\right)^4 \mathrm{Hz}$, where $\alpha'_{*}$ is the fiducial string slope being related directly to the reduced Planck scale $M_\mathrm{Pl}$ as $\alpha'_{*} = 8/M^2_\mathrm{Pl}$. Interestingly enough, $f_\mathrm{GW,peak}$ lies within the $\mathrm{nHz}$ frequency range; hence rendering this primordial GW signal potentially detectable by SKA, NANOGrav and PTA probes at their very low frequency region of their detection bands.

Gamma-ray burst (GRB) afterglow light curves and spectra provide information about the density of the environment, the energy of the explosion, the properties of the particle acceleration process, and the structure of the decelerating jet. Due to the large number of parameters involved, the model can present a certain degree of parameter degeneracy. In this paper, we present synthetic photometric observations of GRB afterglows and model them using the Markov Chain Monte Carlo (MCMC) method. This method has emerged as the preferred approach for analysing and interpreting data in astronomy. We show that, depending on the choice of priors, the parameter degeneracy can go unnoticed by the MCMC method. Furthermore, we apply the MCMC method to analyse the GRB~170817A afterglow. We find that there is a complete degeneracy between the energy of the explosion, the density of the environment, and the microphysical parameters describing the particle acceleration process, which cannot be determined by the afterglow light curve alone. Our results emphasise the importance of gaining a deep understanding of the degeneracy properties which can be present in GRB afterglows models, as well as the limitations of the MCMC method.

We compute the density of a spin-$\frac32$ particle, the raritron, produced at the end of inflation due to gravitational interactions. We consider a background inflaton condensate as the source of this production, mediated by the exchange of a graviton. This production greatly exceeds the gravitational production from the emergent thermal bath during reheating. The relic abundance limit sets an absolute minimum mass for a stable raritron, though there are also model dependent constraints imposed by unitarity. We also examine the case of gravitational production of a gravitino, taking into account the goldstino evolution during reheating. We compare these results with conventional gravitino production mechanisms.

In this work we explore a numerical technique, based on the spherical harmonic decomposition and the discretization of the radial coordinate through \v{C}eby\v{s}\"ev polynomial interpolation, for the computation of quasi-bound states of linear massive scalar and vector perturbations in spinning black hole spacetimes in General Relativity. The aim is studying black hole superradiant instabilities, an energy-extraction mechanism triggered by the presence of massive bosonic fields near black holes, which finds wide applications in constraining scenarios beyond Standard Model and General Relativity. This method does not rely on any separation ans\"atze, thus it can have wide applications. Consequently we extend the technique so that it can be applied also to the computation of massive tensor quasi-bound states in spinning black holes in General Relativity, whose separability ansatz is currently unknown. We also apply it to spinning black holes in scalar-tensor theory non-linearly interacting with plasma, wherein the massless scalar perturbations acquires an effective mass, finding a novel way for constraining scalar-tensor theories.

Extended galaxy surveys revealed that at a large-scale the Universe consists of matter concentrations in the form of galactic clusters, filaments and vast regions devoid of galaxies. In this paper we shall investigate the problem of structure formation in terms of nonlinear dynamics of a self-gravitating N-body system. We reformulate the dynamics of a self-gravitating N-body system in terms of a geodesic flow on a curved Riemannian manifold of dimension 3N equipped by the Maupertuis's metric. The regions of negative sectional curvatures are responsible for the exponential instability of geodesic trajectories and for the chaotic behaviour of the system and the relaxation phenomena, while the regions of the phase space of positive sectional curvatures are responsible for the generation of geodesic focusing and gravitational caustics, the regions of space where the density of matter - particle and galaxies - is larger than in the ambient space of the Universe. We investigate the structure formation dynamics by analysing the stability of geodesic trajectories of N-body systems by means of the Jacobi deviation equation and the geodesic focusing, generation of conjugate points and caustics by means of the Raychaundhuri equation. By solving these equations for an N-body system we estimated the characteristic relaxation time scales of stars in galaxies in globular clusters and demonstrated the geodesic focusing phenomenon in self-gravitating N-body systems that leads to the generation of caustics. The gravitational caustics are space regions where the density of matter is higher than the average density in the surrounding Universe. These regions can represent galaxies, galactic clusters and filaments, and the regions of lower density between caustics represent a relatively empty space, the voids.

Decades of analytic and computational work have demonstrated that a charge immersed in a hot plasma is screened. For both Abelian and non-Abelian interactions, the characteristic screening length $1/m_D$ is set by the so-called Debye mass $m_D \sim g_s T$, proportional to the plasma temperature $T$ and the dimensionless gauge coupling $g_s$. One of the most interesting naturally occurring examples is the quark-gluon plasma (QGP) that filled the early universe prior to the QCD confinement phase transition at $t_{\rm QCD} \sim 10^{-5}\,{\rm s}$. During this early epoch, regimes of strong spacetime curvature are of significant cosmological interest, such as near primordial black holes (PBHs). However, the typical description of Debye screening only applies within Minkowski spacetime, and is therefore insufficient to describe the dynamics of charged plasmas near PBHs or other primordial features. We construct an effective field theory for soft modes of the gauge field $A_\mu^a$ to give a full description of Debye screening in non-Abelian plasmas within arbitrary curved spacetimes, recovering a temperature-dependent Debye mass that exhibits gravitational redshift. We then apply our results to some scenarios of cosmological interest: an expanding FLRW universe and the vicinity of a PBH immersed in a hot QGP.

Motivated by the increasing interest in finding physically viable rotating sources, we present a new class of anisotropic rotating solutions. The energy-momentum tensor compatible with the metric is composed of anisotropic matter with a non-vanishing energy flow around the symmetry axis and vanishing viscosity. The new class of solutions can be used to find new possible sources for the Kerr metric, to obtain new regular black hole solutions and to study galaxies with a central rotating black hole and an halo of dark matter. As an example, we obtain a 5-parameter class of solutions representing a two-way traversable wormhole smoothly matched to the Kerr one and satisfying all energy conditions outside the wormhole for a wide range of parameters, in particular for compact objects. Finally, with a simple modification of the aforementioned solution, we obtain a source for Kerr metric with a throat geometry, non-representing a two-way traversable wormhole and satisfying all energy conditions.

Neutrinos are amongst the most abundant particles in the universe. The fact that they are massive particles proves that the Standard Model is an incomplete theory and needs to be extended. This thesis focuses on the study of neutrino properties from the available data from laboratory experiments and cosmological observations. It also presents sensitivity studies for future neutrino experiments. The first part covers the status of the determination of the oscillation parameters in the three-neutrino framework, the mass ordering and the absolute mass scale. A second part is devoted to other neutrino properties predicted in neutrino mass models. The topics covered include the study of spin-flavour precession in solar neutrinos from non-zero neutrino magnetic moments, the prospects for neutrino non-standard interactions with quarks from the solar sector and the expected limits on CPT violation from a separate analysis of neutrino and antineutrino data. Regarding cosmology, the impact of non-standard interactions with electrons and a non-unitary three-neutrino mixing matrix in the process of neutrino decoupling are addressed. Finally, a third part explores possible connections between dark matter and neutrinos. The existence of dark matter has been proven indirectly via its gravitational effects. This thesis presents two studies that explore how neutrinos could probe the fraction of dark matter that could exist in the form of primordial black holes and the experimental signatures expected from a hypothetical coupling between neutrinos and an ultralight scalar dark matter candidate.

Recently, the scalar-tensor representation of $f (R,T)$ gravity was used to explore gravitationally induced particle production/annihilation. Using the framework of irreversible thermodynamics of open systems in the presence of matter creation/annihilation, the physical and cosmological consequences of this setup were investigated in detail. In this paper, we test observationally the scalar-tensor representation of $f(R,T)$ gravity in the context of the aforementioned framework, using the Hubble and Pantheon+ measurements. The best fit parameters are obtained by solving numerically the modified Friedmann equations of two distinct cosmological models in scalar tensor $f(R, T)$ gravity, corresponding to two different choices of the potential, and by performing a Markov Chain Monte Carlo analysis. The best parameters are used to compute the cosmographic parameters, i.e., the deceleration, the jerk and the snap parameters. Using the output resulting from the Markov Chain Monte Carlo analysis, the cosmological evolution of the creation pressure and of the matter creation rates are presented for both models. To figure out the statistical significance of the studied scalar-tensor $f(R,T)$ gravity, the Bayesian and the corrected Akaike information criteria are used. The latter indicates that the first considered model in scalar tensor $f(R,T)$ gravity is statistically better than $\Lambda$CDM, i.e., it is more favored by observations. Besides, a continuous particle creation process is present in Model 1. On the other hand, for large redshifts, in Model 2 the particle creation rate may become negative, thus indicating the presence of particle annihilation processes. However, both models lead to an accelerating expansion of the Universe at late times, with a deceleration parameter equivalent to that of the $\Lambda$CDM model.

The prospect of unprecedented high-quality data of gravitational waves in the upcoming decades demands a theoretical effort to optimally study and analyze the signals that next generation detectors will provide. Here we study the gravitational wave emission and related dynamics during the inspiralling phase of the Circular Restricted Three Body Problem, a modification of the conventional binary scenario in which a small third object co-rotates with the parent binary system. Specifically, we obtain analytic expressions for the emitted power, frequency variation and other dynamical variables that describe the evolution of the system. As a key highlight, we find that the presence of the third body actually slows down the coalescence of the binary, which can be partially interpreted as an effective rescaling of the binary's chirp-mass. Our analysis assumes semi-Keplerian orbits for the particles and a highly mass asymmetric parent binary needed for the stability of orbits.

Orbital eccentricity is a crucial physical effect to unveil the origin of compact-object binaries detected by ground- and spaced-based gravitational-wave (GW) observatories. Here, we perform for the first time a Bayesian inference study of inspiral-merger-ringdown eccentric waveforms for binary black holes with non-precessing spins using two (instead of one) eccentric parameters: eccentricity and relativistic anomaly. We employ for our study the multipolar effective-one-body (EOB) waveform model SEOBNRv4EHM, and use initial conditions such that the eccentric parameters are specified at an orbit-averaged frequency. We show that this new parametrization of the initial conditions leads to a more efficient sampling of the parameter space. We also assess the impact of the relativistic-anomaly parameter by performing mock-signal injections, and we show that neglecting such a parameter can lead to significant biases in several binary parameters. We validate our model with mock-signal injections based on numerical-relativity waveforms, and we demonstrate the ability of the model to accurately recover the injected parameters. Finally, using standard stochastic samplers employed by the LIGO-Virgo-KAGRA Collaboration, we analyze a set of real GW signals observed by the LIGO-Virgo detectors during the first and third runs. We do not find clear evidence of eccentricity in the signals analyzed, more specifically we measure $e^{\text{GW150914}}_{\text{gw, 10Hz}}= 0.08^{+0.09}_{-0.06}$, $e^{\text{GW151226}}_{\text{gw, 20Hz}}= {0.04}^{+0.05}_{-0.04} $, and $e^{\text{GW190521}}_{\text{gw, 5.5Hz}}= 0.15^{+0.12}_{-0.12}$.

An analysis of contemporary Cosmology is presented, with the aim of identifying the elements present in it according to the scientific program structure created by I. Lakatos. We look at some modern controversies from this point of view and clarify the meaning of issues related to them within this context.

We study the vacuum dynamics of pseudo-Nambu-Goldstone bosons (pNGBs) for $SO(N+1) \rightarrow SO(N)$ spontaneous and explicit symmetry breaking. We determine the magnitude of explicit symmetry breaking consistent with an EFT description of the effective potential at zero and finite temperatures. We expose and clarify novel additional vacuum transitions that can arise for generic pNGBs below the initial scale of $SO(N+1) \rightarrow SO(N)$ spontaneous symmetry breaking, which may have phenomenological relevance. In this respect, two phenomenological scenarios are analyzed: thermal and supercooled dark sector pNGBs. In the thermal scenario the vacuum transition is first-order but very weak. For a supercooled dark sector we find that, depending on the sign of the explicit symmetry breaking, one can have a symmetry-restoring vacuum transition $SO(N-1) \rightarrow SO(N)$ which can be strongly first-order, with a detectable stochastic gravitational wave background signal.