Papers for Wednesday, Dec 14 2022

The recently discovered system Gaia 0007-1605 consisting of a white dwarf with a close brown dwarf companion and a distant white dwarf tertiary very much resembles the triple system containing the first transiting planet candidate around a white dwarf ever discovered: WD 1856+534. We have previously argued that the inner binary in WD 1856+534 most likely formed through common envelope evolution but triple star dynamics represent an alternative scenario. Here we analyze different formation scenarios for Gaia 0007-1605. We reconstructed the potential common envelope evolution of the system and find that assuming standard parameters for the energy budget provides a reasonable solution. In agreement with other close white dwarf + brown dwarf binaries, and in contrast to WD 1856+534, no energy sources other than orbital energy during common envelope evolution are required to understand the current configuration of the system. In addition, using analytical prescriptions for triple dynamics, we show that Von Zeipel-Lidov-Kozai oscillations might have trigger tidal migration due to high eccentricity incursions (e \gtrsim 0.997). We conclude that the inner binary in Gaia 0007-1605, as its sibling WD 1856+534, formed either through common envelope evolution, triple dynamics or a combination of both mechanisms.

In this Letter we revisit the relationship between the fraction of the intra-cluster light (ICL) and both the virial mass and the fraction of Early Type Galaxies in the host halo. This is based on a statistically significant and homogeneous sample of 22 groups and clusters of galaxies in the local Universe ($z \leq 0.05$), obtained with the VST Early-type GAlaxy Survey (VEGAS). Taking advantage of the long integration time and large area of the VEGAS images, we are able to map the galaxy outskirts and ICL down to $\mu_g$ $\geq$ 29-30 mag/arcsec$^2$ and out to hundreds of kpc. With this data-set, we have expanded the sample of ICL measurements, doubling the previous measures available from the literature for z $\leq$ 0.05. The main result of this work is the lack of any significant trend between the fraction of ICL and the virial mass of the host environment, covering a wide range of virial masses ( $\sim$ $10^{12.5} \leq M_{vir} \leq 10^{15.5} M_{\odot}$), in agreement with some theoretical studies. Since the new data points are all derived with the same methodology and from the same observational setup, and all have comparable depth, the large observed scatter indicates an intrinsic variation in the ICL fraction.On the other hand, there is a weak relation between the fraction of ICL and the fraction of Early Type Galaxies in the host halo, where a larger fraction of ICL is found in groups and clusters of galaxies dominated by earlier morphological types, indicating a connection between the ICL and the dynamical state of the host system.

Gravitational wave observations of the inspiral of stellar-mass compact objects into massive black holes (MBHs), extreme mass ratio inspirals (EMRIs), enable precision measurements of parameters such as the MBH mass and spin. The Laser Interferometer Space Antenna is expected to detect sufficient EMRIs to probe the underlying source population, testing theories of the formation and evolution of MBHs and their environments. Population studies are subject to selection effects that vary across the EMRI parameter space, which bias inference results if unaccounted for. This bias can be corrected, but evaluating the detectability of many EMRI signals is computationally expensive. We mitigate this cost by (i) constructing a rapid and accurate neural network interpolator capable of predicting the signal-to-noise ratio of an EMRI from its parameters, and (ii) further accelerating detectability estimation with a neural network that learns the selection function, leveraging our first neural network for data generation. The resulting framework rapidly estimates the selection function, enabling a full treatment of EMRI detectability in population inference analyses. We apply our method to an astrophysically motivated EMRI population model, demonstrating the potential selection biases and subsequently correcting for them. Accounting for selection effects, we predict that LISA will measure the MBH mass function slope to a precision of 8.8%, the CO mass function slope to a precision of 4.6%, the width of the MBH spin magnitude distribution to a precision of 10% and the event rate to a precision of 12% with EMRIs at redshifts below z=6.

We present the selection of high-redshift ($z\gtrsim5.7$) radio-loud (RL) quasi-stellar object (QSO) candidates from the combination of the radio Rapid ASKAP Continuum Survey (RACS; at 888 MHz) and the optical/near-infrared Dark Energy Survey (DES). In particular, we selected six candidates brighter than $S_{\rm 888MHz}>1$ mJy beam$^{-1}$ and ${\rm mag}(z_\mathrm{{DES}})<21.3$ using the dropout technique (in the $i$-band). From this sample, we were able to confirm the high-$z$ nature ($z\sim6.1$) of two sources, which are now among the highest-redshift RL QSOs currently known. Based on our Gemini-South/GMOS observations, neither object shows a prominent Ly$\alpha$ emission line. This suggests that both sources are likely to be weak emission-line QSOs hosting radio jets and would therefore further strengthen the potential increase of the fraction of weak emission-line QSOs recently found in the literature. However, further multiwavelength observations are needed to constrain the properties of these QSOs and of their relativistic jets. From the discovery of these two sources, we estimated the space density of RL QSOs in the redshift range $5.9<z<6.4$ to be 0.13$^{+0.18}_{-0.09}$ and found it to be consistent with the expectations based on our current knowledge of the blazar population up to $z\sim5$.

We investigate the accuracy and precision of triaxial dynamical orbit models by fitting two dimensional mock observations of a realistic N-body merger simulation resembling a massive early-type galaxy with a supermassive black hole (SMBH). We show that we can reproduce the triaxial N-body merger remnant's correct black hole mass, stellar mass-to-light ratio and total enclosed mass (inside the half-light radius) for several different tested orientations with an unprecedented accuracy of 5-10%. Our dynamical models use the entire non-parametric line-of-sight velocity distribution (LOSVD) rather than parametric LOSVDs or velocity moments as constraints. Our results strongly suggest that state-of-the-art integral-field projected kinematic data contain only minor degeneracies with respect to the mass and anisotropy recovery. Moroever, this also demonstrates the strength of the Schwarzschild method in general. We achieve the proven high recovery accuracy and precision with our newly developed modeling machinery by combining several advancements: (i) our new semi-parametric deprojection code probes degeneracies and allows to constrain the viewing angles of a triaxial galaxy; (ii) our new orbit modeling code SMART uses a 5-dim orbital starting space to representatively sample in particular near-Keplerian orbits in galaxy centers; (iii) we use a generalised information criterion AICp to optimise the smoothing and to compare different mass models to avoid biases that occur in $\chi^2$-based models with varying model flexibilities.

Hubble Space Telescope imaging for 26 giant early-type galaxies, all drawn from the MAST archive, is used to carry out photometry of their surrounding globular cluster (GC) systems. Most of these targets are Brightest Cluster Galaxies (BCGs) and their distances range from 24 to 210 Mpc. The catalogs of photometry, completed with DOLPHOT, are publicly available. The GC color indices are converted to [Fe/H] through a combination of 12-Gyr SSP (Single Stellar Population) models and direct spectroscopic calibration of the fiducial color index (F475W-F850LP). All the resulting metallicity distribution functions (MDFs) can be accurately matched by bimodal Gaussian functions. The GC systems in all the galaxies also exhibit shallow metallicity gradients with projected galactocentric distance that average $Z \sim R_{gc}^{-0.3}$. Several parameters of the MDFs including the means, dispersions, and blue/red fractions are summarized. Perhaps the most interesting new result is the trend of blue/red GC fraction with galaxy mass, which connects with predictions from recent simulations of GC formation within hierarchical assembly of large galaxies. The observed trend reveals two major transition stages: for low-mass galaxies, the metal-rich (red) GC fraction $f(red)$ increases steadily with galaxy mass, until halo mass $M_h \simeq 3 \times 10^{12} M_{\odot}$. Above this point, more than half the metal-poor (blue) GCs come from accreted satellites and $f(red)$ starts declining. But above a still higher transition point near $M_h \simeq 10^{14} M_{\odot}$, the data hint that $f(red)$ may start to increase again because the metal-rich GCs also become dominated by accreted systems.

We use the Giant Metrewave Radio Telescope (GMRT) Cold-HI AT $z\approx1$ (CAT$z1$) survey, a 510 hr HI 21cm emission survey of galaxies at $z=0.74-1.45$, to report the first measurements of atomic hydrogen (HI) scaling relations at $z\approx1$. We divide our sample of 11,419 blue star-forming galaxies at $z\approx1$ into three stellar mass ($M_*$) subsamples and obtain detections (at $\geq 4\sigma$ significance) of the stacked HI 21cm emission signal from galaxies in all three subsamples. We fit a power-law relation to the measurements of the average HI mass ($M_{HI}$) in the three stellar-mass subsamples to find that the slope of the $M_{HI}-M_{*}$ relation at $z\approx1$ is consistent with that at $z\approx0$. However, we find that the $M_{HI}-M_{*}$ relation has shifted downwards from $z\approx1$ to $z\approx0$, by a factor of $3.54\pm0.48$. Further, we find that the HI depletion timescales ($t_{dep,HI}$) of galaxies in the three stellar-mass subsamples are systematically lower than those at $z\approx0$, by factors of $\approx2-4$. We divide the sample galaxies into three specific star-formation rate (sSFR) subsamples, again obtaining $\geq 4\sigma$ detections of the stacked HI 21cm emission signal in all three subsamples. We find that the relation between the ratio of HI mass to stellar mass and the sSFR evolves between $z\approx1$ and $z\approx0$. Unlike the efficiency of conversion of molecular gas to stars, which does not evolve significantly with redshift, we find that the efficiency with which HI is converted to stars is much higher for star-forming galaxies at $z\approx1$ than those at $z\approx0$.

The fraction of ionizing photons escaping from galaxies, $f_{esc}$, is at the same time a crucial parameter in modelling reionization and a very poorly known quantity, especially at high redshift. Recent observations are starting to constrain the values of $f_{esc}$ in low-z star-forming galaxies, but the validity of this comparison remains to be verified. Applying at high-z the empirical relation between $f_{esc}$ and the UV slope trends derived from the Low-z Lyman Continuum Survey, we use the DELPHI semi-analytical galaxy formation model to estimate the global ionizing emissivity of high-z galaxies, which we use to compute the resulting reionization history. We find that both the global ionizing emissivity and reionization history match the observational constraints. Assuming that the low-z correlations hold during the epoch of reionization, we find that galaxies with $-16 \lesssim M_{UV} \lesssim -13.5$ are the main drivers of reionization. We derive a population-averaged $\langle f_{esc} \rangle \simeq 8\%, 10\%, 20\%$ at z=4.5, 6, 8.

To unravel the dominant cause of the weak emission line in a subset of optically selected radio-quiet 'weak emission line quasars' (WLQs), we have investigated the possibility of an underdeveloped broad line region (BLR). For this, we have modeled spectral energy distributions (SED) of 61 WLQs by using their optical and infrared (IR) photometric observations from SDSS and WISE respectively. SED fit consists of various emission components, including the luminosity from the dusty torus ($L_{tor}$). For comparison with the normal quasar, we have used a control sample of 55 QSOs for each WLQs matching in emission redshift and SDSS r-band. Based on our measurement of $L_{tor}$, we found a decrement of $42\pm2$\% in IR-luminosity in WLQs w.r.t the control sample of normal QSOs. Using $L_{tor}$/$L_{bol}$ as the measure of torus covering factor ($CF_{tor}$) we found a similar decrement in WLQs covering factor, with their $CF_{tor}$ distribution being significantly different w.r.t. the normal QSOs with a KS-test $P_{null}$ of $4.27 \times 10^{-14}$. As dusty torus and BLR covering factors are expected to be of a similar order in AGN, our results suggest that the BLR in the WLQs is underdeveloped and could be a dominant cause of the weakness of their emission line. As a result, our analysis gives support to the models of WLQs based on the evolution scenario being in an early stage of AGNs.

Neutral hydrogen (HI) emission exhibits complex morphology that encodes rich information about the physics of the interstellar medium (ISM). We apply the scattering transform (ST) to characterize HI emission structure via a set of compact and interpretable coefficients, and find a connection between HI emission morphology and HI cold neutral medium (CNM) phase content. Where HI absorption measurements are unavailable, the HI phase structure is typically estimated from the emission via spectral line decomposition. Here we present the first probe of CNM content using measures solely derived from HI emission spatial information. We apply the scattering transform to GALFA-HI data at high Galactic latitudes (|b|>30 deg), and compare the resulting coefficients to CNM fraction measurements derived from archival HI emission and absorption spectra. We quantify the correlation between the ST coefficients and measured CNM fraction (fCNM), and find that HI emission morphology encodes substantial fCNM-correlating information, and that ST-based metrics for small-scale linearity are particularly predictive of fCNM. This is further corroborated by the enhancement of $I_{857}/N_{HI}$ ratio with larger ST measures of small-scale linearity. These results are consistent with the picture that regions with higher CNM content are more populated with small-scale filamentary HI structures. Our work illustrates a physical connection between HI morphology and phase content, and suggests that future phase decomposition methods can be improved by making use of both HI spectral and spatial information.

We present the first systematic broad-band X-ray study of super-Eddington accretion onto SMBHs with simultaneous {\it NuSTAR} and {\it XMM-Newton} or {\it Swift}/XRT observations of a sample of 8 super-Eddington accreting AGN with Eddington ratio $1<\lambda_{\rm Edd}<426$. We find that the SEAMBHs show a steep primary continuum slope as expected for sources accreting in the super Eddington regime, mostly dominated by relativistic reflection. The Iron K$\alpha$ emission lines of the sources of our sample show relativistic broadening. In addition the equivalent widths of the narrow components of the Iron K$\alpha$ lines follow the 'X-ray Baldwin' effect, also known as the 'Iwasawa-Taniguchi' effect. We found a statistically significant correlation between the photon-index of the primary power-law and the Eddington ratio, consistent with past studies. Moreover, as expected for super-Eddington sources, the median value of the reflection fraction of the sources we analysed is a factor $\sim 2$ higher than the median reflection fraction value of the type\,1 AGN from the BASS sample. We are able to estimate the coronal temperature for three sources of our sample: Mrk\,382 ($kT_e=7.8$\,keV), PG\,0026+129 ($kT_e=39$\,keV) and IRAS\,04416+1215 ($kT_e=3$\,keV). Looking at the position of the SEAMBHs sources of our sample in the compactness-temperature diagram it appears that in super-Eddington AGN, as for lower Eddington ratio AGN, the X-ray corona is controlled by pair production and annihilation.

Action-based dynamical modelling, using stars as dynamical tracers, is an excellent diagnostic to estimate the underlying axisymmetric matter distribution of the Milky Way. However, the Milky Way's bar causes non-axisymmetric resonance features in the stellar disc. Using Roadmapping (an action-based dynamical modelling framework to estimate the gravitational potential and the stellar distribution function), we systematically quantify the robustness of action-based modelling in the presence of a bar. We construct a set of test-particle simulations of barred galaxies (with varying bar properties), and apply Roadmapping to different survey volumes (with varying azimuthal position, size) drawn from these barred models. For realistic bar parameters, the global potential parameters are still recovered to within ${\sim \! 1 \! - \! 20}$ percent. However, with increasing bar strength, the best-fit values of the parameters progressively deviate from their true values. This happens due to a combination of radial heating, radial migration, and resonance overlap phenomena in our bar models. Furthermore, the azimuthal location and the size of the survey volumes play important roles in the successful recovery of the parameters. Survey volumes along the bar major axis produce larger (relative) errors in the best-fit parameter values. In addition, the potential parameters are better recovered for survey volumes with larger spatial coverage. As the Sun is located just ${\sim \! 28 \! - \! 33}$ degrees behind the bar's major axis, an estimate for the bar-induced systematic bias -- as provided by this study -- is therefore crucial for future modelling attempts of the Milky Way.

Rocky planets orbiting M-dwarf stars are among the most promising and abundant astronomical targets for detecting habitable climates. Planets in the M-dwarf habitable zone are likely synchronously rotating, such that we expect significant day-night temperature differences, and potentially limited fractional habitability. Previous studies have focused on scenarios where fractional habitability is confined to the substellar or "eye" region, but in this paper we explore the possibility of planets with terminator habitability, defined by the existence of a habitable band at the transition between a scorching dayside and a glacial nightside. Using a global climate model, we show that for water-limited planets it is possible to have scorching temperatures in the "eye" and freezing temperatures on the nightside, while maintaining a temperate climate in the terminator region, due to a reduced atmospheric energy transport. Whereas on water-rich planets, increasing stellar flux leads to increased atmospheric energy transport and a reduction in day-night temperature differences, such that the terminator does not remain habitable once the dayside temperatures approach runaway or moist greenhouse limits. We also show that, while water-abundant simulations may result in larger fractional habitability, they are vulnerable to water loss through cold-trapping on the nightside surface or atmospheric water vapor escape, suggesting that even if planets were formed with abundant water, their climates could become water-limited and subject to terminator habitability.

Ultraviolet (UV) transmission spectroscopy probes atmospheric escape, which has a significant impact on planetary atmospheric evolution. If unaccounted for, interstellar medium absorption (ISM) at the position of specific UV lines might bias transit depth measurements, and thus potentially affect the (non-)detection of features in transmission spectra. Ultimately, this is connected to the so called ``resolution-linked bias'' (RLB) effect. We present a parametric study quantifying the impact of unresolved or unconsidered ISM absorption in transit depth measurements at the position of the MgII h&k resonance lines (i.e. 2802.705 {\AA} and 2795.528 {\AA} respectively) in the near-ultraviolet spectral range. We consider main-sequence stars of different spectral types and vary the shape and amount of chromospheric emission, ISM absorption, and planetary absorption, as well as their relative velocities. We also evaluate the role played by integration bin and spectral resolution. We present an open-source tool enabling one to quantify the impact of unresolved or unconsidered MgII ISM absorption in transit depth measurements. We further apply this tool to a few already or soon to be observed systems. On average, we find that ignoring ISM absorption leads to biases in the MgII transit depth measurements comparable to the uncertainties obtained from the observations published to date. However, considering the bias induced by ISM absorption might become necessary when analysing observations obtained with the next generation space telescopes with UV coverage (e.g. LUVOIR, HABEX), which will provide transmission spectra with significantly smaller uncertainties compared to what obtained with current facilities (e.g. HST).

We present ultraviolet (UV) to near-infrared (NIR) observations and analysis of the nearby Type Ia supernova SN 2021fxy. Our observations include UV photometry from Swift/UVOT, UV spectroscopy from HST/STIS, and high-cadence optical photometry with the Swope 1-m telescope capturing intra-night rises during the early light curve. Early $B-V$ colours show SN 2021fxy is the first "shallow-silicon" (SS) SN Ia to follow a red-to-blue evolution, compared to other SS objects which show blue colours from the earliest observations. Comparisons to other spectroscopically normal SNe Ia with HST UV spectra reveal SN 2021fxy is one of several SNe Ia with flux suppression in the mid-UV. These SNe also show blue-shifted mid-UV spectral features and strong high-velocity Ca II features. One possible origin of this mid-UV suppression is the increased effective opacity in the UV due to increased line blanketing from high velocity material, but differences in the explosion mechanism cannot be ruled out. Among SNe Ia with mid-UV suppression, SNe 2021fxy and 2017erp show substantial similarities in their optical properties despite belonging to different Branch subgroups, and UV flux differences of the same order as those found between SNe 2011fe and 2011by. Differential comparisons to multiple sets of synthetic SN Ia UV spectra reveal this UV flux difference likely originates from a luminosity difference between SNe 2021fxy and 2017erp, and not differing progenitor metallicities as suggested for SNe 2011by and 2011fe. These comparisons illustrate the complicated nature of UV spectral formation, and the need for more UV spectra to determine the physical source of SNe Ia UV diversity.

Cool outflows are now commonly observed in galaxies, but their physical origin and driving mechanism remain unclear. Active galactic nucleus (AGN) feedback can potentially accelerate cool galactic outflows via cosmic rays (CR) and radiation pressure on dust. Here we investigate the relative importance of CR and radiation feedback in AGNs, and we analyse the physical conditions for outflow launching as a function of the black hole accretion flow mode. We assume CRs from AGN jet origin and consider the analogy with Galactic X-ray binaries, whereby the jet is prominent at low accretion rates (hard state) and quenched at high accretion rates (soft state). We show that CR-driven outflows can be powered at low accretion rates and at large radii, whereas radiation pressure-driven outflows dominate at high accretion rates and small radii. Thus the two AGN feedback mechanisms -- CRs and radiation pressure on dust -- may play complementary roles in driving cool outflows on galactic scales. The transition from radiation pressure-driven outflows at higher accretion rates to CR-driven outflows at lower accretion rates likely corresponds to a transition in the underlying accretion flow modes (from a radiatively efficient accretion disc to a radiatively inefficient jet-dominated flow) over cosmic time.

We present a study of nearby dwarf galaxies selected from the ALFALFA blind HI survey. A primary goal of the project was to utilize a non-standard selection method with the hope of detecting previously unrecognized extremely metal-poor (XMP) galaxies. The study was motivated by the recent discovery of two XMP galaxies $-$ Leo P and Leoncino $-$ which were both originally found via the ALFALFA survey. We have obtained narrowband H$\alpha$ images for 42 dwarf systems, many of which are located in the local void in front of the Pisces-Perseus Supercluster. Spectra for eleven of the best candidates resulted in the determination of metal abundances for ten of the systems. None were found to be extremely metal poor, although one system (AGC 123350) was found to have an oxygen abundance of log(O/H)+12 = 7.46, or $\sim$6\% solar. One of the galaxies in our sample exhibits a high oxygen abundance for its luminosity, suggesting the possibility that it may have a tidal origin.

Machine learning has been successfully applied in varied field but whether it is a viable tool for determining the distance to molecular clouds in the Galaxy is an open question. In the Galaxy, the kinematic distance is commonly employed as the distance to a molecular cloud. However, there is a problem in that for the inner Galaxy, two different solutions, the ``Near'' solution, and the ``Far'' solution, can be derived simultaneously. We attempted to construct a two-class (``Near'' or ``Far'') inference model using a Convolutional Neural Network (CNN), a form of deep learning that can capture spatial features generally. In this study, we used the CO dataset toward the 1st quadrant of the Galactic plane obtained with the Nobeyama 45-m radio telescope (l = 62-10 degree, |b| < 1 degree). In the model, we applied the three-dimensional distribution (position-position-velocity) of the 12CO (J=1-0) emissions as the main input. The dataset with ``Near'' or ``Far'' annotation was made from the HII region catalog of the infrared astronomy satellite WISE to train the model. As a result, we could construct a CNN model with a 76% accuracy rate on the training dataset. By using the model, we determined the distance to molecular clouds identified by the CLUMPFIND algorithm. We found that the mass of the molecular clouds with a distance of < 8.15 kpc identified in the 12CO data follows a power-law distribution with an index of about -2.3 in the mass range of M >10^3 Msun. Also, the detailed molecular gas distribution of the Galaxy as seen from the Galactic North pole was determined.

During solar flares a tremendous amount of magnetic energy is released and transported through the Sun's atmosphere and out into the heliosphere. Despite over a century of study, many unresolved questions surrounding solar flares are still present. Among those are how does the solar plasma respond to flare energy deposition, and what are the important physical processes that transport that energy from the release site in the corona through the transition region and chromosphere? Attacking these questions requires the concert of advanced numerical simulations and high spatial-, temporal-, and spectral- resolution observations. While flares are 3D phenomenon, simulating the NLTE flaring chromosphere in 3D and performing parameter studies of 3D models is largely outwith our current computational capabilities. We instead rely on state-of-the-art 1D field-aligned simulations to study the physical processes that govern flares. Over the last decade, data from the Interface Region Imaging Spectrograph (IRIS) have provided the crucial observations with which we can critically interrogate the predictions of those flare loop models. Here in Paper 2 of a two-part review of IRIS and flare loop models, I discuss how forward modelling flares can help us understand the observations from IRIS, and how IRIS can reveal where our models do well and where we are likely missing important processes, focussing in particular on the plasma properties, energy transport mechanisms, and future directions of flare modelling.

TOI-2076 b is a sub-Neptune-sized planet ($R=2.39 \pm 0.10 \mathrm{R_\oplus}$) that transits a young ($204 \pm 50 \mathrm{MYr}$) bright ($V = 9.2$) K-dwarf hosting a system of three transiting planets. Using spectroscopic observations with the NEID spectrograph on the WIYN 3.5 m Telescope, we model the Rossiter-McLaughlin effect of TOI-2076 b, and derive a sky-projected obliquity of $\lambda=-3_{-15}^{+16\:\circ}$. Using the size of the star ($R=0.775 \pm0.015 \mathrm{R_\odot}$), and the stellar rotation period ($P_{\mathrm{rot}}=7.27\pm0.23$ days), we estimate a true obliquity of $\psi=18_{-9}^{+10\:\circ}$ ($\psi < 34^\circ$ at 95% confidence), demonstrating that TOI-2076 b is on a well-aligned orbit. Simultaneous diffuser-assisted photometry from the 3.5 m Telescope at Apache Point Observatory rules out flares during the transit. TOI-2076 b joins a small but growing sample of young planets in compact multi-planet systems with well-aligned orbits, and is the fourth planet with an age $\lesssim 300$ Myr in a multi-transiting system with an obliquity measurement. The low obliquity of TOI-2076 b and the presence of transit timing variations in the system suggest the TOI-2076 system likely formed via convergent disk migration in an initially well-aligned disk.

We develop a "dual cone" model for millimeter wavelength emission near a spinning black hole. The model consists of optically thin, luminous cones of emission, centered on the spin axis, which are meant to represent jet walls. The resulting image is dominated by a thin ring. We first consider the effect of black hole's spin on the image, and show that the dominant effect is to displace the ring perpendicular to the projection of the spin axis on the sky by $2 a_* \sin i + \mathcal{O}(a_*^3)$. This effect is lower order in $a_*$ than changes in the shape and size of the photon ring itself, but is undetectable without a positional reference. We then show that the centerline of the jet can provide a suitable reference: its location is exactly independent of spin if the observer is outside the cone, and nearly independent of spin if the observer is inside the cone. If astrophysical uncertainties can be controlled for, then spin displacement is large enough to be detectable by future space VLBI missions. Finally, we consider ring substructure in the dual cone model and show that features in total intensity are not universal and depend on the cone opening angle.

The aim of this work is to provide new insights on the dynamics associated to the resonances which arise as a consequence of the coupling of the effect due to the oblateness of the Earth and the Solar Radiation Pressure (SRP) effect for an uncontrolled object with moderate to high area-to-mass ratio. Analytical estimates for the location of the resulting resonant equilibrium points are provided, together with formulas to compute the maximum amplitude of the corresponding variation in the eccentricity, as a function of the initial conditions of the object and of its area-to-mass ratio. The period of the variations of the eccentricity and inclinations due to such resonances is estimated using classical formulas. A classification based on the strength of the SRP resonances is provided. The estimates presented in the paper are validated using numerical tools, including the use of Fast Lyapunov Indicators to draw phase portraits and bifurcation diagrams. Many FLI maps depicting the location and overlapping of SRP resonances are presented. The results from this paper suggest that SRP resonances could be modeled in the context of either the Extended Fundamental Model by [1] or the Second Fundamental Model by [2]. [1] S. Breiter. Extended fundamental model of resonance. Celestial Mechanics and Dynamical Astronomy, 85:209-218, 03 (2003). doi: 10.1023/A:1022569419866 [2] J. Henrard and A. Lemaitre. A second fundamental model for resonance. Celestial Mechanics, 30:197-218, (1983). doi: 10.1007/BF01234306

We use the MaNGA integral-field spectroscopic survey of low-redshift galaxies to compare the stellar populations of the bulge and disc components, identified from their Sersic profiles, for various samples of galaxies. Bulge dominated regions tend to be more metal-rich and have slightly older stellar ages than their associated disc dominated regions. The metallicity difference is consistent with the deeper gravitational potential in bulges relative to discs, which allows bulges to retain more of the metals produced by stars. The age difference is due to star formation persisting longer in discs relative to bulges. Relative to galaxies with lower stellar masses, galaxies with higher stellar masses tend to have bulge dominated regions that are more metal-rich and older (in light-weighted measurements) than their disc dominated regions. This suggests high-mass galaxies quench from the inside out, while lower-mass galaxies quench across the whole galaxy simultaneously. Early-type galaxies tend to have bulge dominated regions the same age as their disc dominated regions, while late-type galaxies tend to have disc dominated regions significantly younger than their bulge dominated regions. Central galaxies tend to have a greater metallicity difference between their bulge dominated regions and disc dominated regions than satellite galaxies at similar stellar mass. This difference may be explained by central galaxies being subject to mergers or extended gas accretion bringing new, lower-metallicity gas to the disc, thereby reducing the average metallicity and age of the stars; quenching of satellite discs may also play a role.

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region being one of the youngest (<10 Myr), closest (280-320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, PanSTARRS). Methods. We use Bayesian methodologies accounting for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbound stars. Their mass distributions are broadly compatible with that of Chabrier for masses >0.1 $M_\odot$ and do not show hints of over-abundance of low-mass stars in NGC1333 with respect to IC348. These groups' ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than in previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 $M_\odot$), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star-formation scenario.

Based on the first 13.4 years of Fermi science data in the energy range from 300 MeV to 500 GeV, we discovered a bright GeV gamma-ray source with a $\sim$5.64 $\sigma$ detection, named Fermi~J1242.5+3236, which has an offset of about $0.0996^{\circ}$ from a nearby star-forming galaxy NGC 4631. When using the 12 yrs' data, the detection significance of Fermi~J1242.5+3236 is about 4.72 $\sigma$. Fermi~J1242.5+3236 is a steady point source without significant temporal variability and has a hard gamma-ray photon index of about $-$1.60$\pm 0.24$. The spatial offset and the hard gamma-ray spectrum disfavors this source as the diffuse gamma-ray emission from this galaxy. This new source might has a possible origin of an unidentified background blazar, which is more likely a high-synchrotron-peaked blazar (HSP) for its hard gamma-ray photon index. A follow-up optical observation would help distinguish origin of Fermi~J1242.5+3236.

With unparalleled rotational stability, millisecond pulsars (MSPs) serve as ideal laboratories for numerous astrophysical studies, many of which require precise knowledge of the distance and/or velocity of the MSP. Here, we present the astrometric results for 18 MSPs of the "MSPSR$\pi$" project focusing exclusively on astrometry of MSPs, which includes the re-analysis of 3 previously published sources. On top of a standardized data reduction protocol, more complex strategies (i.e., normal and inverse-referenced 1D interpolation) were employed where possible to further improve astrometric precision. We derived astrometric parameters using sterne, a new Bayesian astrometry inference package that allows the incorporation of prior information based on pulsar timing where applicable. We measured significant ($>3\,\sigma$) parallax-based distances for 15 MSPs, including $0.81\pm0.02\,$kpc for PSR J1518+4904 -- the most significant model-independent distance ever measured for a double neutron star system. For each MSP with a well-constrained distance, we estimated its transverse space velocity and radial acceleration. Among the estimated radial accelerations, the updated ones of PSR J1012+5307 and PSR J1738+0333 impose new constraints on dipole gravitational radiation and the time derivative of Newton's gravitational constant. Additionally, significant angular broadening was detected for PSR J1643-1224, which offers an independent check of the postulated association between the HII region Sh 2-27 and the main scattering screen of PSR J1643-1224. Finally, the upper limit of the death line of $\gamma$-ray-emitting pulsars is refined with the new radial acceleration of the hitherto least energetic $\gamma$-ray pulsar PSR J1730-2304.

Recent observations show that the metallicity of the broad line region ($Z_{\rm BLR}$) in active galactic nuclei (AGNs) is solar-to-supersolar, which is positively correlated with the mass of supermassive black holes ($M_{\rm BH}$) and does not evolve with redshift up to $z \sim 7$. We revisit the $M_{\rm BH}-Z_{\rm BLR}$ correlation with more AGNs with $M_{\rm BH}\sim 10^{6-8} M_{\odot}$ and find that the positive correlation become flat in low-mass range. It is known that outer part of accretion disks is gravitationally unstable and can fragment into stars. Considering the star formation and supernovae (SNe) in the outer AGN disk, we calculate the metal enrichment and find that positive $M_{\rm BH}-Z_{\rm BLR}$ correlation can be roughly reproduced if the stellar mass distribution is ``top-heavy". We find that the observed BLR size is more or less similar to the self-gravity radius of the AGN disk, which suggests that the BLR may be closely correlated with the underlying accretion process.

Dual active galactic nuclei (AGN), which are the manifestation of two actively accreting supermassive black holes (SMBHs) hosted by a pair of merging galaxies, are a unique laboratory for studying the physics of SMBH feeding and feedback during an indispensable stage of galaxy evolution. In this work, we present NOEMA CO(2-1) observations of seven kpc-scale dual-AGN candidates drawn from a recent Chandra survey of low-redshift, optically classified AGN pairs. These systems are selected because they show unexpectedly low 2-10 keV X-ray luminosities for their small physical separations signifying an intermediate-to-late stage of merger. Circumnuclear molecular gas traced by the CO(2-1) emission is significantly detected in 6 of the 7 pairs and 10 of the 14 nuclei, with an estimated mass ranging between $(0.2 - 21) \times10^9\rm~M_{\odot}$. The primary nuclei, i.e., the ones with the higher stellar velocity dispersion, tend to have a higher molecular gas mass than the secondary. Most CO-detected nuclei show a compact morphology, with a velocity field consistent with a kpc-scale rotating structure. The inferred hydrogen column densities range between $5\times10^{21} - 2\times10^{23}\rm~cm^{-2}$, but mostly at a few times $10^{22}\rm~cm^{-2}$, in broad agreement with those derived from X-ray spectral analysis. Together with the relatively weak mid-infrared emission, the moderate column density argues against the prevalence of heavily obscured, intrinsically luminous AGNs in these seven systems, but favors a feedback scenario in which AGN activity triggered by a recent pericentric passage of the galaxy pair can expel circumnuclear gas and suppress further SMBH accretion.

JWST/MIRI imaging of the nearby galaxies IC 5332, NGC 628, NGC 1365 and NGC 7496 from PHANGS reveals a richness of gas structures that in each case form a quasi-regular network of interconnected filaments, shells and voids. We examine whether this multi-scale network of structure is consistent with the fragmentation of the gas disk through gravitational instability. We use FilFinder to detect the web of filamentary features in each galaxy and determine their characteristic radial and azimuthal spacings. These spacings are then compared to estimates of the most Toomre-unstable length (a few kpc), the turbulent Jeans length (a few hundred pc) and the disk scale height (tens of pc) reconstructed using PHANGS-ALMA observations of the molecular gas as a dynamical tracer. Our analysis of the four galaxies targeted in this work indicates that Jeans-scale structure is pervasive. Future work will be essential for determining how the structure observed in gas disks impacts not only the rate and location of star formation but also how stellar feedback interacts positively or negatively with the surrounding multi-phase gas reservoir.

We assess a neural network (NN) method for reconstructing 3D cosmological density and velocity fields (target) from discrete and incomplete galaxy distributions (input). We employ second-order Lagrangian Perturbation Theory to generate a large ensemble of mock data to train an autoencoder (AE) architecture with a Mean Squared Error (MSE) loss function. The AE successfully captures nonlinear features arising from gravitational dynamics {and} the discreteness of the galaxy distribution. In comparison, the traditional Wiener filter (WF) reconstruction exhibits a stronger suppression of the power on smaller scales and contains regions of unphysical negative densities. In the density reconstruction, the reduction of the AE MSE relative to the WF is $\sim 15 \%$, whereas for the velocity reconstruction a relative reduction of up to a factor of two can be achieved. The AE is particularly superior to the WF towards the tails of the density and velocity distributions. In fact, trained with an MSE loss, any NN estimate approaches the unbiased mean of the underlying target given the input. This implies a slope of unity in the linear regression of the true on the NN-reconstructed field. Only for the special case of Gaussian fields, the NN and WF estimates are equivalent. Nonetheless, we also recover a linear regression slope of unity for the WF with non-Gaussian fields.

High-redshift primordial galaxies have recently been found with evolved stellar populations and complex star-formation histories reaching back to 250 Myr after the Big Bang. Their intense bursts of star-formation appear to be interspersed with sustained periods of strong quenching, however the processes underlying this evolutionary behaviour remain unclear. Unlike later epochs, galaxies in the early Universe are not located in large associations like clusters. Instead, they co-evolve with their developing circumgalactic halo as relatively isolated ecosystems. Thus, the mechanisms that could bring about the downfall of their star-formation are presumably intrinsic, and feedback processes associated with their intense starburst episodes likely play an important role. Cosmic rays are a viable agent to deliver this feedback, and could account for the star-formation histories inferred for these systems. The cosmic ray impact on galaxies may be investigated using the wealth of multi-wavelength data soon to be obtained with the armada of new and upcoming facilities. Complementary approaches to probe their action across the electromagnetic spectrum can be arranged into a distance ladder of cosmic ray feedback signatures. With a clear understanding of how cosmic ray activity in primordial systems can be traced, it will be possible to extend this ladder to high redshifts and map-out the role played by cosmic rays in shaping galaxy evolution over cosmic time.

Context. The Sun acts as a cornerstone of stellar physics. Thanks to spectroscopic, helioseismic and neutrino flux observations, we can use the Sun as a laboratory of fundamental physics in extreme conditions. The conclusions we draw are then used to inform and calibrate evolutionary models of all other stars in the Universe. However, solar models are in tension with helioseismic constraints. The debate on the ``solar problem'' has hitherto led to numerous publications discussing potential issues with solar models and abundances. Aims. Using the recently suggested high-metallicity abundances for the Sun, we investigate whether standard solar models, as well as models with macroscopic transport reproducing the solar surface lithium abundances and analyze their properties in terms of helioseismic and neutrino flux observations. Methods. We compute solar evolutionary models and combine spectroscopic and helioseismic constraints as well as neutrino fluxes to investigate the impact of macroscopic transport on these measurements. Results. When high-metallicity solar models are calibrated to reproduce the measured solar lithium depletion, tensions arise with respect to helioseismology and neutrino fluxes. This is yet another demonstration that the solar problem is also linked to the physical prescriptions of solar evolutionary models and not to chemical composition alone. Conclusions. A revision of the physical ingredients of solar models is needed in order to improve our understanding of stellar structure and evolution. The solar problem is not limited to the photospheric abundances if the depletion of light elements is considered. In addition, tighter constraints on the solar beryllium abundance will play a key role in the improvement of solar models.

After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multi-dimensional simulations. Obtaining nucleosynthesis yields based on such multi-dimensional core-collapse supernova (CCSN) simulations, however, is not straightforward and frequently tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamics simulation. Reconstruction can be done in a number of ways and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to a better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction, to the point when nuclear statistical equilibrium last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing.

We model helium-rich stars with solar metallicity ($X=0.7,\:Z=0.02$) progenitors that evolve to form AM Canum Venaticorum systems through a helium-star formation channel, with the aim to explain the observed properties of Gaia14aae and ZTFJ1637+49. We show that semi-degenerate, H-exhausted ($X\leq 10^{-5}$), He-rich ($Y\approx0.98$) donors can be formed after a common envelope evolution (CEE) phase if either additional sources of energy are used to eject the common envelope, or a different formalism of CEE is implemented. We follow the evolution of such binary systems after the CEE phase using the Cambridge stellar evolution code, when they consist of a He-star and a white dwarf accretor, and report that the mass, radius, and mass-transfer rate of the donor, the orbital period of the system, and the lack of hydrogen in the spectrum of Gaia14aae and ZTFJ1637+49 match well with our modelled trajectories wherein, after the CEE phase Roche lobe overflow is governed not only by the angular momentum loss (AML) owing to gravitational wave radiation ($\mathrm{AML_{GR}}$) but also an additional AML owing to $\alpha-\Omega$ dynamos in the donor. This additional AML is modelled with our double-dynamo (DD) model of magnetic braking in the donor star. We explain that this additional AML is just a consequence of extending the DD model from canonical cataclysmic variable donors to evolved donors. We show that none of our modelled trajectories match with Gaia14aae or ZTFJ1637+49 if the systems are modelled only with $\mathrm{AML_{GR}}$.

An analysis of the CH molecule in non-local thermodynamic equilibrium (NLTE) is performed for the physical conditions of cool stellar atmospheres typical of red giants (log g = 2.0, Teff = 4500 K) and the Sun. The aim of the present work is to explore whether the G-band of the CH molecule, which is commonly used in abundance diagnostics of Carbon-Enhanced Metal-Poor (CEMP) stars, is sensitive to NLTE effects. LTE and NLTE theoretical spectra are computed with the MULTI code. We use one-dimensional (1D) LTE hydrostatic MARCS model atmospheres with parameters representing eleven red giant stars with metallicities ranging from [Fe/H] = -4.0 to [Fe/H] = 0.0 and carbon-to-iron ratios [C/Fe] = 0.0, +0.7, +1.5, and +3.0. The CH molecule model is represented by 1981 energy levels, 18377 radiative bound-bound transitions, and 932 photo-dissociation reactions. The rates due to transitions caused by collisions with free electrons and hydrogen atoms are computed using classical recipes. Our calculations suggest that NLTE effects in the statistical equilibrium of the CH molecule are significant and cannot be neglected for precision spectroscopic analysis of C abundances. The NLTE effects in the G-band increase with decreasing metallicity. We show that the C abundances are always under-estimated if LTE is assumed. The NLTE corrections to C abundance inferred from the CH feature range from +0.04 dex for the Sun to +0.21 dex for a red giant with metallicity [Fe/H] = -4.0. Departures from the LTE assumption in the CH molecule are non-negligible and NLTE effects have to be taken into account in the diagnostic spectroscopy based on the CH lines. We show here that the NLTE effects in the optical CH lines are non-negligible for the Sun and red giant stars, but further calculations are warranted to investigate the effects in other regimes of stellar parameters.

The James Webb Space Telescope (JWST), which has recently become operational, is capable of detecting objects at record-breaking redshifts, $z \gtrsim 15$. This is a crucial advance for observational cosmology, as at these redshifts the differences between alternative cosmological models manifest themselves in the most obvious way. In recent years, some observational hints have emerged indicating that the Standard Cosmological Model could require correcting. One of these hints is related to the discovery of remote galaxies whose redshifts correspond to the very young Universe (less than one billion years after the Big Bang) but which are similar to nearby galaxies. The issue is that such galaxies in the early Universe do not have enough time to evolve into something similar to the late-Universe galaxies. JWST observations of high-redshift objects are expected to shed light on the origin of this issue. Here we provide results on performing the ``angular diameter -- redshift'' cosmological test for the first JWST observation data. We compare this result with predictions of the standard $\Lambda$CDM cosmological model and some static cosmological models, including Zwicky's ``tired-light'' model. The latter is currently assumed to be ruled out by observations. We challenge this assumption and show that a static model can provide a natural and straightforward way of solving the puzzle of the well-evolved galaxies and better agreements with the results of the JWST ``angular diameter -- redshift'' test at higher redshifts than the correcting evolution model within the $\Lambda$CDM framework. We discuss several cosmological tests that will be important for further research on the possibility of revising the expanding Universe paradigm.

Context. The details of cosmic-ray transport have a strong impact on galaxy evolution. The peak of the cosmic-ray energy distribution is observable in the radio continuum using the electrons as proxy. Aims. We measure the length that the cosmic-ray electrons (CRE) are transported during their lifetime in the nearby galaxy M 51 across one order of magnitude in cosmic-ray energy (approximately 1-10 GeV). To this end we use new ultra-low frequency observations from the LOw Frequency ARay (LOFAR) at 54 MHz and ancillary data between 144 and 8350 MHz. Methods. As the the CRE originate from supernova remnants, the radio maps are smoothed in comparison to the distribution of the star formation. By convolving the map of the star-formation rate (SFR) surface density with a Gaussian kernel, we can linearise the radio-SFR relation. The best-fitting convolution kernel is then our estimate of the CRE transport length. Results. We find that the CRE transport length increases at low frequencies, as expected since the CRE have longer lifetimes. The CRE transport length is $l_{\rm CRE} = \sqrt{4Dt_{\rm syn}}$, where $D$ is the isotropic diffusion coefficient and $t_{\rm syn}$ is the CRE lifetime as given by synchrotron and inverse Compton losses. We find that the data can be well fitted by diffusion, where $D=(2.14\pm 0.13) \times 10^{28}~\rm cm^2\,s^{-1}$. With $D\propto E^{0.001\pm 0.185}$, the diffusion coefficient is independent of the CRE energy $E$ in the range considered. Conclusions. Our results suggest that the transport of GeV-cosmic ray electrons in the star-forming discs of galaxies is governed by energy-independent diffusion.

We present a study of synthetic observations of polarized dust emission at 353 GHz as seen by an observer within a cavity. The cavity's physical properties are comparable to the Local Bubble in which the Solar System is embedded. We find that the combination of the local density enhancement and the coherent magnetic field in the walls of the cavity makes the Local Bubble a translucent polarization filter to the emission coming from beyond its domains. This underlines the importance of studying the Local Bubble in detail, as it is a foreground screen to Galactic astronomy and cosmology. We find that the magnetic field lines inferred from synthetic dust polarization data are qualitatively in agreement, at low and middle latitudes (|b|<60 deg), with the all-sky maps of polarized emission at 353 GHz by the Planck satellite. This is in agreement with current models of the magnetic field structure in the Local Bubble walls. However, we also find that the polarized signal coming from the selected cavity is insufficient to explain the regular polarization orientation patterns observed by Planck toward the polar Galactic caps. This may be due to the advanced stage of evolution of the candidate bubble and to the limited simulation domain. Finally, we show that from our synthetic polarization maps, it is difficult to distinguish between an open and a closed Galactic cap using the inferred magnetic field morphology alone.

The era of Gravitational-Wave (GW) astronomy will grant the detection of the astrophysical GW background from unresolved mergers of binary black holes, and the prospect of probing the presence of primordial GW backgrounds. In particular, the low-frequency tail of the GW spectrum for causally-generated primordial signals (like a phase transition) offers an excellent opportunity to measure unambiguously cosmological parameters as the equation of state of the universe, or free-streaming particles at epochs well before recombination. We discuss whether this programme is jeopardised by the uncertainties on the astrophysical GW foregrounds that coexist with a primordial background. We detail the motivated assumptions under which the astrophysical foregrounds can be assumed to be known in shape, and only uncertain in their normalisation. In this case, the sensitivity to a primordial signal can be computed by a simple and numerically agile procedure, where the optimal filter function subtracts the components of the astrophysical foreground that are close in spectral shape to the signal. We show that the degradation of the sensitivity to the signal in presence of astrophysical foregrounds is limited to a factor of a few, and only around the frequencies where the signal is closer to the foregrounds. Our results highlight the importance of modelling the contributions of eccentric or intermediate-mass black hole binaries to the GW background, to consolidate the prospects to perform precision cosmology with primordial GW backgrounds.

The deepest all-sky survey available in the $\gamma$-ray band - the last release of the Fermi-LAT catalogue (4FGL-DR3) based on the data accumulated in 12 years, contains more than 6600 sources. The largest population among the sources is blazar subclass - 3743, $60.1\%$ of which are classified as BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs), while the rest are listed as blazar candidates of uncertain type (BCU) as their firm optical classification is lacking. The goal of this study is to classify BCUs using different machine learning algorithms which are trained on the spectral and temporal properties of already classified BL Lacs and FSRQs. Artificial Neural Networks, \textit{XGBoost} and LightGBM algorithms are employed to construct predictive models for BCU classification. Using 18 input parameters of 2219 BL Lacs and FSRQs, we train (80\% of the sample) and test (20\%) these algorithms and find that LightGBM model, state-of-the-art classification algorithm based on gradient boosting decision trees, provides the highest performance. Based on our best model, we classify 825 BCUs as BL Lac candidates and 405 as FSRQ candidates, however, 190 remain without a clear prediction but the percentage of BCUs in 4FGL is reduced to 5.1\%. The $\gamma$-ray photon index, synchrotron peak frequency, and high energy peak frequency of a large sample are used to investigate the relationship between FSRQs and BL Lacs (LBLs, IBLs, and HBLs).

In this paper, we forecast the expected detection rates and redshift distributions of gravitationally lensed gravitational waves (GWs) from three different mass distributions of primordial black holes (PBHs) and two stellar formation models of astrophysical black holes (ABHs) in the context of DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) and it's smaller scale version B-DECIGO. It suggests that DECIGO will be able to detect $10^4-10^5$ GW signals from such binary black holes (BBHs) each year and the event rate distributions for PBHs will differ from those for ABHs due to their different merger rate with respect to redshift. The large number of event rates make $5-100$ detections of lensed GW signals being possible. After considering the gravitational lensing effect, the difference between the detection rates and distributions for PBHs and ABHs will be more significant. Therefore, this can be served as a complementary method to distinguish PBHs from ABHs.

GX 301-2 is a high-mass X-ray binary (HMXB) with strong stellar outflows. The evolution of these binaries can be closely related with the interstellar environment due to strong wind interactions. We try to constrain the energy injected in the interstellar medium by GX 301-2 through stellar wind using HAWK-I and Herschel data. We analysed HAWK-I images in four different filters (Br$\gamma$, H$_2$, J, and Ks) and tried to retrieve signatures of the impact of GX 301-2 on its environment. We used Herschel data to outline the interstellar medium and the Gaia DR3 catalogue to infer the proper motion of GX 301-2. Finally, we estimated the energy injected in the interstellar medium since the first supernova event of the HMXB. Using both HAWK-I and Herschel images, we deduce an approximation of the total mass injected from GX~301-2 in the interstellar medium of $M_{\rm inj} = 3.05 ^{+0.05}_{-0.03} 10^{-2} M_{\odot}$.

It was recently proposed that there exists a "gateway" in the orbital parameter space through which Centaurs transition to Jupiter-family Comets (JFCs). Further studies have implied that the majority of objects that eventually evolve into JFCs should leave the Centaur population through this gateway. This may be naively interpreted as gateway Centaurs being pristine progenitors of JFCs. This is the point we want to address in this work. We show that the opposite is true: gateway Centaurs are, on average, more thermally processed than the rest of the population of Centaurs crossing Jupiter's orbit. Using a dynamically-validated JFC population, we find that only $\sim 20\%$ of Centaurs pass through the gateway prior to becoming JFCs, in accordance with previous studies. We show that more than half of JFC dynamical clones entering the gateway for the first time have already been JFCs -they simply avoided the gateway on their first pass into the inner solar system. By coupling a thermal evolution model to the orbital evolution of JFC dynamical clones, we find a higher than 50\% chance that the layer currently contributing to the observed activity of gateway objects has been physically and chemically altered, due to previously sustained thermal processing. We further illustrate this effect by examining dynamical clones that match the present-day orbits of 29P/Schwassmann-Wachmann 1, P/2019 LD2 (ATLAS), and P/2008 CL94 (Lemmon).

One of the common approximations in long-term evolution studies of small bodies is the use of circular orbits averaging the actual eccentric ones, facilitating the coupling of processes with very different timescales, such as the orbital changes and the thermal processing. Here we test a number of averaging schemes for elliptic orbits in the context of the long-term evolution of comets, aiming to identify the one that best reproduces the elliptic orbits' heating patterns and the surface and subsurface temperature distributions. We use a simplified thermal evolution model applied on simulated comets both on elliptic and on their equivalent averaged circular orbits, in a range of orbital parameter space relevant to the inner solar system. We find that time averaging schemes are more adequate than spatial averaging ones. Circular orbits created by means of a time average of the equilibrium temperature approximate efficiently the subsurface temperature distributions of elliptic orbits in a large area of the orbital parameter space, rendering them a powerful tool for averaging elliptic orbits.

Solar flares are transient yet dramatic events in the atmosphere of the Sun, during which a vast amount of magnetic energy is liberated. This energy is subsequently transported through the solar atmosphere or into the heliosphere, and together with coronal mass ejections flares comprise a fundamental component of space weather. Thus, understanding the physical processes at play in flares is vital. That understanding often requires the use of forward modelling in order to predict the hydrodynamic and radiative response of the solar atmosphere. Those predictions must then be critiqued by observations to show us where our models are missing ingredients. While flares are of course 3D phenomenon, simulating the flaring atmosphere including an accurate chromosphere with the required spatial scales in 3D is largely beyond current computational capabilities, and certainly performing parameter studies of energy transport mechanisms is not yet tractable in 3D. Therefore, field-aligned 1D loop models that can resolve the relevant scales have a crucial role to play in advancing our knowledge of flares. In recent years, driven in part by the spectacular observations from the Interface Region Imaging Spectrograph (IRIS), flare loop models have revealed many interesting features of flares. For this review I highlight some important results that illustrate the utility of attacking the problem of solar flares with a combination of high quality observations, and state-of-the-art flare loop models, demonstrating: (1) how models help to interpret flare observations from IRIS, (2) how those observations show us where we are missing physics from our models, and (3) how the ever increasing quality of solar observations drives model improvements. Here in Paper 1 of this two part review I provide an overview of modern flare loop models, and of electron-beam driven mass flows during solar flares.

We report the detection of a high density of redshift $z\approx 10$ galaxies behind the foreground cluster Abell 2744, selected from imaging data obtained recently with NIRCam onboard JWST by three programs -- GLASS-JWST, UNCOVER, and DDT-2756. To ensure robust estimates of the lensing magnification $\mu$, we use an improved version of our model that exploits the first epoch of NIRCam images and newly obtained MUSE spectra, and avoids regions with $\mu>5$ where the uncertainty may be higher. We detect seven bright $z\approx 10$ galaxies with demagnified rest-frame $-22 \lesssim M_{\rm UV}\lesssim -19$ mag, over an area of $\sim37$ sq. arcmin. Taking into account photometric incompleteness and the effects of lensing on luminosity and cosmological volume, we find that the density of $z\approx 10$ galaxies in the field is about $10\times$ ($3\times$) larger than the average at $M_{UV}\approx -21 (-20)$ mag reported so far. The density is even higher when considering only the GLASS-JWST data, which are the depeest and the least affected by magnification and incompleteness. The GLASS-JWST field contains 5 out of 7 galaxies, distributed along an apparent filamentary structure of 2 Mpc in projected length, and includes a close pair of candidates with $M_{\rm UV}< -20$ mag having a projected separation of only 16 kpc. These findings suggest the presence of a $z\approx 10$ overdensity in the field. In addition to providing excellent targets for efficient spectroscopic follow-up observations, our study confirms the high density of bright galaxies observed in early JWST observations, but calls for multiple surveys along independent lines of sight to achieve an unbiased estimate of their average density and a first estimate of their clustering.

ESO 137-G006 is the brightest cluster galaxy (BCG) of the cool-core and dynamically young Norma cluster. We discover an atomic hydrogen (HI) absorption line associated with this BCG using the Australia Telescope Compact Array. We estimate a gas column density of $ \approx (1.3 \pm 0.2) \times 10^{20}\,T_{\rm{spin}}$ atoms cm$^{-2}$ with spin temperature, $T_{\rm{spin}} \leq 194$ K, consistent with the HI properties of other early-type galaxies and cool-core cluster BCGs. The relationship between the presence of cold gas and a cluster cooling flow is unclear. Our results support the scenario that ESO 137-G006 may be a recent arrival to the cluster centre and not the original BCG. This scenario is consistent with the observed spatial alignment of the BCG's wide-angle tail radio lobes with Norma's X-ray sub-cluster and the significant line-of-sight velocity offset between the mean velocity of Norma and that of the BCG.

We have derived luminosity functions, and set constraints on the UV luminosity and SFR density from z~17 to z~8, using the three most-studied JWST/NIRCam data sets, the SMACS0723, GLASS Parallel, and CEERS fields. We first used our own selections on two independent reductions of these datasets using the latest calibrations. 18 z~8, 12 z~10, 5 z~13, and 1 z~17 candidate galaxies are identified over these fields in our primary reductions, with a similar number of candidates in our secondary reductions. We then use these two reductions, applying a quantitative discriminator, to segregate the full set of z>~8 candidates reported over these fields from the literature, into three different samples, ``robust,'' ``solid,'' and ``possible''. Using all of these samples we then derive UV LF and luminosity density results at $z\geq8$, finding substantial differences. For example, including the full set of ``solid'' and ``possible'' z>~12 candidates from the literature, we find UV luminosity densities which are ~7x and ~20x higher than relying on the ``robust'' candidates alone. These results indicate the evolution of the UV LF and luminosity densities at z>~8 is still extremely uncertain, emphasizing the need for spectroscopy and deeper NIRCam+optical imaging to obtain reliable results. Nonetheless, even with the very conservative ``robust'' approach to selections, both from our own and those of other studies, we find the luminosity density from luminous (M(UV)<-19) galaxies to be ~2x larger than is easily achievable using constant star-formation efficiency models, similar to what other early JWST results have suggested.

This paper presents the results of a search for neutrinos that are spatially and temporally coincident with 22 unique, non-repeating Fast Radio Bursts (FRBs) and one repeating FRB (FRB121102). FRBs are a rapidly growing class of Galactic and extragalactic astrophysical objects that are considered a potential source of high-energy neutrinos. The IceCube Neutrino Observatory's previous FRB analyses have solely used track events. This search utilizes seven years of IceCube's cascade events which are statistically independent of the track events. This event selection allows probing of a longer range of extended timescales due to the low background rate. No statistically significant clustering of neutrinos was observed. Upper limits are set on the time-integrated neutrino flux emitted by FRBs for a range of extended time-windows.

Detection of massive binary black hole (BBH) mergers at high redshifts is a target for LISA space mission. While the individual masses of a BBH merger are redshifted, the mass ratio of BBH mergers is independent of their redshift. Therefore, if there is an independent correlation between the mass ratio and redshift, such a relationship can be used to i) infer the redshift of the merging binaries, and together with the luminosity distance measurement ($D_L$), constrain the expansion rate of the universe at high redshifts $H(z)$, and ii) constrain models of supermassive black hole seed formation in the universe assuming a fixed cosmology. We discuss why there is an expected relation between the mass ratio of the massive BBHs with their redshift and show the clues for this relation by analyzing cosmological hydrodynamical simulations of BBH mergers. This approach opens up the possibility of directly measuring the expansion rate at redshift $z \approx 2$ with LISA for the first time. Moreover, we discover a trend between seed mass and mass ratio of massive BBHs which by itself is a major result that could be exploited to constrain the formation scenarios of supermassive BH seeds.

Sunquakes are seismic emissions visible on the solar surface, associated with some solar flares. Although discovered in 1998, they have only recently become a more commonly detected phenomenon. Despite the availability of several manual detection guidelines, to our knowledge, the astrophysical data produced for sunquakes is new to the field of Machine Learning. Detecting sunquakes is a daunting task for human operators and this work aims to ease and, if possible, to improve their detection. Thus, we introduce a dataset constructed from acoustic egression-power maps of solar active regions obtained for Solar Cycles 23 and 24 using the holography method. We then present a pedagogical approach to the application of machine learning representation methods for sunquake detection using AutoEncoders, Contrastive Learning, Object Detection and recurrent techniques, which we enhance by introducing several custom domain-specific data augmentation transformations. We address the main challenges of the automated sunquake detection task, namely the very high noise patterns in and outside the active region shadow and the extreme class imbalance given by the limited number of frames that present sunquake signatures. With our trained models, we find temporal and spatial locations of peculiar acoustic emission and qualitatively associate them to eruptive and high energy emission. While noting that these models are still in a prototype stage and there is much room for improvement in metrics and bias levels, we hypothesize that their agreement on example use cases has the potential to enable detection of weak solar acoustic manifestations.

We have analysed Optical Gravitational Lensing Experiment photometry for first overtone classical Cepheids in the Large and Small Magellanic Clouds in search for additional periodicities beyond radial modes. We have used standard consecutive prewhitening technique in some cases followed by time-dependent prewhitening. We report new candidates for double-mode radial pulsations. However, majority of signals we have detected cannot be interpreted in terms of radial modes. We report 516 double-periodic stars with period ratios, $P_{\rm x}/P_{\rm 1O}$, in the range 0.60 and 0.65. We study the properties of this class and implications for model explaining these periodicities. We also report 28 stars in which additional variability is of longer period, below radial fundamental mode, with median period ratio, $P_{\rm 1O}/P_{\rm x}$, of 0.684. This class is an analogue of a class known in RR Lyrae stars. Hundreds of other signals were detected that cannot be attributed to radial modes or the above-mentioned classes. Statistical properties of these signals are analysed. We suggest that majority of these signals correspond to non-radial modes. In particular, a significant fraction can be attributed to non-radial modes of moderate degrees, tightly connected to a class with period ratios in between 0.60 and 0.65. In tens of stars, close to radial mode frequency, relatively large-amplitude and coherent signals are observed, that may represent yet another class. In 27 stars periodic modulation of pulsation was detected. Differences in additional frequency content between the two Clouds are discussed.

Large imaging surveys, such as LSST, rely on photometric redshifts and tomographic binning for 3x2pt analyses that combine galaxy clustering and weak lensing. In this paper, we propose a method for optimizing the tomographic binning choice for the lens sample of galaxies. We divide the CosmoDC2 and Buzzard simulated galaxy catalogs into a training set and an application set, where the training set is non-representative in a realistic way, and then estimate photometric redshifts for the application sets. The galaxies are sorted into redshift bins covering equal intervals of redshift or comoving distance, or with an equal number of galaxies in each bin, and we consider a generalized extension of these approaches. We find that placing equal numbers of galaxies in each bin produces the highest signal to noise of the initial binning choices, but that the choice of bin edges can be further optimized. We then train a neural network classifier to identify galaxies that are either highly likely to have accurate photometric redshift estimates, or highly likely to be sorted into the correct redshift bin. The neural network classifier is used to remove poor redshift estimates from the sample, and the results are compared to the case when none of the sample is removed. We find that the choice of bin edges has a larger impact than the sample selection with the NNCs, but the NNCs are able to recover ~70% of the loss in signal to noise that occurs when a non-representative training sample is used.

One of the most mysterious astrophysical states is the common envelope (CE) phase of binary evolution, in which two stars are enshrouded by the envelope shed by one of them. Interactions between the stars and the envelope shrinks the orbit. The CE can lead to mergers or to a subsequent phase of interactions. Mergers may involve any combination of two compact objects and/or stars. Some involving white dwarfs, may produce Type Ia supernovae, while merging neutron stars may yield gamma-ray bursts, and merging compact objects of all kinds produce gravitational radiation. Since CEs can arise from a variety of different initial conditions, and due to the complexity of the processes involved, it is difficult to predict their end states. When many systems are being considered, as in population synthesis calculations, conservation principles are generally employed. Here we use angular momentum in a new way to derive a simple expression for the final orbital separation. This method provides advantages for the study of binaries and is particularly well-suited to higher order multiples, now considered to be important in the genesis of potential mergers. Here we focus on CEs in binaries, and the follow-up paper extends our formalism to multiple star systems within which a CE occurs.

Intermediate mass stars (IMSs) represent the link between low-mass and high-mass stars and cover a key mass range for giant planet formation. In this paper, we present a spectroscopic survey of 241 young IMS candidates with IR-excess, the most complete unbiased sample to date within 300 pc. We combined VLT/X-Shooter spectra with BVR photometric observations and Gaia DR3 distances to estimate fundamental stellar parameters such as Teff, mass, radius, age, and luminosity. We further selected those stars within the intermediate mass range 1.5 <= Mstar/Msun <= 3.5 and discarded old contaminants. We used 2MASS and WISE photometry to study the IR-excesses of the sample, finding 92 previously unidentified stars with IR-excess. We classified this sample into 'protoplanetary', 'hybrid candidates' and 'debris' discs based on their observed fractional excess at 12microns, finding a new population of 17 hybrid disc candidates. We studied inner disc dispersal timescales for {\lambda} < 10{\mu}m and found very different trends for IMSs and low mass stars (LMSs). IMSs show excesses dropping fast during the first 6 Myrs independently of the wavelength, while LMSs show consistently lower fractions of excess at the shortest wavelengths and increasingly higher fractions for longer wavelengths, with slower dispersal rates. In conclusion, this study demonstrates empirically that IMSs dissipate their inner discs very differently than LMSs, providing a possible explanation for the lack of short period planets around IMSs.

Upcoming surveys of cosmic structures will probe scales close to the cosmological horizon, which opens up new opportunities for testing the cosmological concordance model to high accuracy. In particular, constraints on the squeezed bispectrum could rule out the single-field hypothesis during inflation. However, the squeezed bispectrum is also sensitive to dynamical effects of general relativity as well as interactions of matter with residual radiation from the early Universe. In this paper, we present a relativistic simulation pipeline that includes these relativistic effects consistently. We produce light cones and calculate the observed number counts of cold dark matter for five redshift bins between $z=0.55$-$2.25$. We compare the relativistic results against reference Newtonian simulations by means of angular power- and bispectra. We find that the dynamical relativistic effects scale roughly inversely proportional to the multipole in the angular power spectrum, with an amplitude of $0.5\%$ to $5\%$ of the total power. By using a smoothing method applied to the binned bispectrum we detect the Newtonian bispectrum with very high significance. The purely relativistic part of the matter bispectrum, obtained by subtracting the Newtonian bispectrum from the relativistic one, is detected with a significance of $\sim 3\,\sigma$, mostly limited by cosmic variance. Our relativistic pipeline for modelling ultra-large scales yields gauge-independent results as we compute observables consistently on the past light cone, while the Newtonian treatment employs approximations that leave some residual gauge dependence. A gauge-invariant approach is required in order to meet the expected level of precision of forthcoming probes of cosmic structures on ultra-large scales.

The understanding of novae, the thermonuclear eruptions on the surfaces of white dwarf stars in binaries, has recently undergone a major paradigm shift. Though the bolometric luminosity of novae was long thought to be solely attributed to runaway nuclear burning, recent GeV gamma-ray observations have supported the notion that a significant portion of the luminosity could come from radiative shocks. More recently, observations of novae have lent evidence that these shocks are acceleration sites for hadrons for at least some types of novae. In this scenario, a flux of neutrinos may accompany the observed gamma rays. As the gamma rays from most novae have only been observed up to a few GeV, novae have previously not been considered as targets for neutrino telescopes, which are most sensitive at and above TeV energies. Here, we present the first search for neutrinos from novae with energies between a few GeV and 10 TeV using IceCube-DeepCore, a densely instrumented region of the IceCube Neutrino Observatory with a reduced energy threshold. We search both for a correlation between gamma-ray and neutrino emission as well as between optical and neutrino emission from novae. We find no evidence for neutrino emission from the novae considered in this analysis and set upper limits for all gamma-ray detected novae.

The Oort cloud is thought to be a reservoir of icy planetesimals and the source of long-period comets (LPCs) implanted from the outer Solar System during the time of giant planet formation. The abundance of rocky ice-free bodies is a key diagnostic of Solar System formation models as it can distinguish between ``massive" and ``depleted" proto-asteroid belt scenarios and thus disentangle competing planet formation models. Here we report a direct observation of a decimeter-sized ($\sim2$ kg) rocky meteoroid on a retrograde LPC orbit ($e \approx 1.0$, i = $121^{\circ}$). During its flight, it fragmented at dynamic pressures similar to fireballs dropping ordinary chondrite meteorites. A numerical ablation model fit produces bulk density and ablation properties also consistent with asteroidal meteoroids. We estimate the flux of rocky objects impacting Earth from the Oort cloud to be $1.08^{+2.81}_{-0.95} \mathrm{meteoroids/10^6 km^2/yr}$ to a mass limit of 10 g. This corresponds to an abundance of rocky meteoroids of $\sim6^{+13}_{-5}$\% of all objects originating in the Oort cloud and impacting Earth to these masses. Our result gives support to migration-based dynamical models of the formation of the Solar System which predict that significant rocky material is implanted in the Oort cloud, a result not explained by traditional Solar System formation models.

In this work, we focus on plasma discharges produced between two electrodes with a high potential difference, resulting in the ionization of the neutral particles supporting a current in the gaseous medium. At low currents and low temperatures, this process can create luminescent emissions: the so-called glow and corona discharges. The parallel plate geometry used in Townsend's (1900) theory lets us develop a theoretical formalism, with explicit solutions for the critical voltage effectively reproducing experimental Paschen curves. However, most discharge processes occur in non-parallel plate geometries, such as discharges between grains or ice particles in multiphase flows. Here, we propose a generalization of the classic parallel plate configurations to concentric spherical and coaxial cylindrical geometries in Earth, Mars, Titan, and Venus atmospheres. In a spherical case, a small radius effectively represents a sharp tip rod, while larger, centimeter-scale radii represents rounded or blunted tips. Similarly, in a cylindrical case, a small radius would correspond to a thin wire. We solve continuity equations in the gap and estimate a critical radius and minimum breakdown voltage that allows ionization of neutral gas and formation of a glow discharge. We show that glow coronae form more easily in Mars's low-pressure, $CO_2$-rich atmosphere than in Earth's high-pressure atmosphere. Additionally, we present breakdown criteria for Titan and Venus. We further demonstrate that critical voltage minima occur at 0.5 cm$\cdot$Torr for all three investigated geometries, suggesting easier initiation around millimeter-size particles in dust and water clouds and could be readily extended to examine other multiphase flows with inertial particles.

Many recent works have shown that the angular resolution of ground-based detectors is too poor to characterize the anisotropies of the stochastic gravitational-wave background (SGWB). For this reason, we asked ourselves if a constellation of space-based instruments could be more suitable. We consider the Laser Interferometer Space Antenna (LISA), a constellation of multiple LISA-like clusters, and the Deci-hertz Interferometer Gravitational-wave Observatory (DECIGO). Specifically, we test whether these detector constellations can probe the anisotropies of the SGWB. For this scope, we considered the SGWB produced by two astrophysical sources: merging compact binaries and a recently proposed scenario for massive black-hole seed formation through multiple mergers of stellar remnants. We find that measuring the angular power spectrum of the SGWB anisotropies is almost unattainable. However, it turns out that it could be possible to probe the SGWB anisotropies through cross-correlation with the CMB fluctuations. In particular, we find that a constellation of two LISA-like detectors and CMB-S4 can marginally constrain the cross-correlation between the CMB lensing convergence and the SGWB produced by the black hole seed formation process. Moreover, we find that DECIGO can probe the cross-correlation between the CMB lensing and the SGWB from merging compact binaries.

In this letter we uncover a new parametric resonance of axionic cosmic strings. This process is triggered by the presence on the string of internal mode excitations that resonantly amplify the amplitude of transverse displacements of the string. We study this process by running numerical simulations that demonstrate the existence of this phenomenon in a $(3+1)$ dimensional lattice field theory and compare the results with the analytic expectations for the effective Lagrangian of the amplitude of these modes and their interactions. Finally, we also analyze the massless and massive radiation produced by these excited strings and comment on its relevance for the interpretation of the results of current numerical simulations of axionic comic string networks.

An isotropic stochastic background of nanohertz gravitational waves creates excess residual power in pulsar-timing-array datasets, with characteristic inter-pulsar correlations described by the Hellings-Downs function. These correlations appear as nondiagonal terms in the noise covariance matrix, which must be inverted to obtain the pulsar-timing-array likelihood. Searches for the stochastic background, which require many likelihood evaluations, are therefore quite computationally expensive. We propose a more efficient method: we first compute approximate posteriors by ignoring cross correlations, and then reweight them to exact posteriors via importance sampling. We show that this technique results in accurate posteriors and marginal likelihood ratios, because the approximate and exact posteriors are similar, which makes reweighting especially accurate. The Bayes ratio between the marginal likelihoods of the exact and approximate models, commonly used as a detection statistic, is also estimated reliably by our method, up to ratios of at least $10^6$.

Classical novae are thermonuclear explosions in stellar binary systems, and important sources of $^{26}$Al and $^{22}$Na. While gamma rays from the decay of the former radioisotope have been observed throughout the Galaxy, $^{22}$Na remains untraceable. The half-life of $^{22}$Na (2.6 yr) would allow the observation of its 1.275 MeV gamma-ray line from a cosmic source. However, the prediction of such an observation requires good knowledge of the nuclear reactions involved in the production and destruction of this nucleus. The $^{22}$Na($p,\gamma$)$^{23}$Mg reaction remains the only source of large uncertainty about the amount of $^{22}$Na ejected. Its rate is dominated by a single resonance on the short-lived state at 7785.0(7) keV in $^{23}$Mg. In the present work, a combined analysis of particle-particle correlations and velocity-difference profiles is proposed to measure femtosecond nuclear lifetimes. The application of this novel method to the study of the $^{23}$Mg states, combining magnetic and highly-segmented tracking gamma-ray spectrometers, places strong limits on the amount of $^{22}$Na produced in novae, explains its non-observation to date in gamma rays (flux < 2.5x$10^{-4}$ ph/(cm$^2$s)), and constrains its detectability with future space-borne observatories.

In this paper, we investigate the structure of hybrid stars consisting of hadrons (neutrons, protons, sigmas, lambdas), leptons (electrons, muons), and quarks (up, down, strange). We use a relativistic mean-field (RMF) model namely the Sigma-omega-rho model for the hadronic phase and the MIT bag model as well as the NJL model for the quark phase. In addition, Maxwell and Gibbs conditions are employed to investigate the hadron-Quark phase transition. Finally, by obtaining the mass-radius relation, $ M (M_{sun}) \leqslant 2.07 $ is predicted for such hybrid stars.

We investigate the entanglement due to geometric corrections in particle creation during inflation. To do so, we propose a single-field inflationary scenario, nonminimally coupled to the scalar curvature of spacetime. We require particle production to be purely geometric, setting to zero the Bogolubov coefficients and computing the $S$ matrix associated to spacetime perturbations, which are traced back to inflaton fluctuations. The corresponding particle density leads to a nonzero entanglement entropy whose effects are investigated at primordial time of Universe evolution. The possibility of modeling our particle candidate in terms of dark matter is discussed. The classical back-reaction of inhomogeneities on the homogeneous dynamical background degrees of freedom is also studied and quantified in the slow-roll regime.

Superconducting cosmic strings can exhibit longitudinal, pinching instabilities in some regions of the parameter space. We make predictions about the onset of this instability using the thin string approximation (TSA) and develop an improved analysis that remains applicable for small wavelength perturbations, where the TSA breaks down. We use simulations of perturbed strings to assess the accuracy of the TSA, test the predictions of our new analysis and demonstrate an improvement over previous methods in the literature. Notably, it appears that the instabilities are typically present for a larger range of magnetic strings than previously expected, and we show examples of pinching instabilities also occurring in electric strings. However, both our simulations and predictions agree that strings near the chiral limit are free from pinching instabilities and in particular our results support our previously published claim that vortons can be stable to all classical perturbations if they are sufficiently large.

The development and adoption of artificial intelligence (AI) technologies in space applications is growing quickly as the consensus increases on the potential benefits introduced. As more and more aerospace engineers are becoming aware of new trends in AI, traditional approaches are revisited to consider the applications of emerging AI technologies. Already at the time of writing, the scope of AI-related activities across academia, the aerospace industry and space agencies is so wide that an in-depth review would not fit in these pages. In this chapter we focus instead on two main emerging trends we believe capture the most relevant and exciting activities in the field: differentiable intelligence and on-board machine learning. Differentiable intelligence, in a nutshell, refers to works making extensive use of automatic differentiation frameworks to learn the parameters of machine learning or related models. Onboard machine learning considers the problem of moving inference, as well as learning, onboard. Within these fields, we discuss a few selected projects originating from the European Space Agency's (ESA) Advanced Concepts Team (ACT), giving priority to advanced topics going beyond the transposition of established AI techniques and practices to the space domain.

We develop a new perturbation method for determining a class of time transfer functions in a stationary spacetime when its metric is a small deformation of a background metric for which the time transfer functions are known in a closed form. The perturbation terms are expressed as line integrals along the null geodesic paths of the background metric. Unlike what happens with the other procedures proposed until now, the time transfer functions obtained in this way are completely free of unbounded terms and do not generate any enhancement in the light travel time. Our procedure proves to be very efficient when the background metric is a linearized Schwarzschild-like metric. Its application to an isolated body slowly rotating about an axis of symmetry leads to integrals which can be calculated with any symbolic computer program. Explicit expressions are obtained for the mass dipole and quadrupole moments and for the leading gravitomagnetic term induced by the spin of the body. A brief numerical discussion is given for the 2002 Cassini experiment.

This thesis presents research exploring aspects of inflation and cosmology in the context of inflation models in which an inflaton is non-minimally coupled to the Ricci scalar, or is considered in conjunction with a term quadratic in the Ricci scalar. We consider a $\phi^{2}$ Palatini inflation model in $R^{2}$ gravity and investigate whether this model can overcome some of the problems of the original $\phi^{2}$ chaotic inflation model. We investigate the compatibility of this model with the observed CMB when treated as an effective theory of inflation in quantum gravity by examining the constraints on the model parameters arising due to Planck-suppressed potential corrections and reheating. Additionally, we consider two possible reheating channels and assess their viability in relation to the constraints on the size of the coupling to the $R^{2}$ term. We present an application of the Affleck-Dine mechanism, in which quadratic $B$-violating potential terms generate the asymmetry, with a complex inflaton as the Affleck-Dine field. We derive the $B$ asymmetry generated in the inflaton condensate analytically and numerically. We use the present-day asymmetry to constrain the size of the $B$-violating mass term and derive an upper bound on the inflaton mass in order for the Affleck-Dine dynamics to be compatible with non-minimally coupled inflation in the metric and Palatini formalisms. We demonstrate the existence of a new class of inflatonic Q-balls in a non-minimally coupled Palatini inflation model, through an analytical derivation of the Q-ball equation and numerical confirmation of the existence of solutions, and derive a range of the inflaton mass squared within which the model can inflate and produce Q-balls. We derive analytical estimates of the properties of these Q-balls, explore the effects of curvature, and discuss observational signatures of the model.