Papers for Friday, Sep 16 2022

The discovery of two galaxy candidates at $z\sim 13$, HD1 and HD2, laid the foundation for a new race to study the early Universe. These galaxies exhibit large UV luminosities, and previous investigations suggested that they are either powered by a supermassive black hole or by an extreme, transient burst of star formation. Given their uncertain nature, the critical question we investigate is whether these sources could be lensed by a hitherto undetected, faint foreground galaxy. We find that at the current limiting magnitude with which HD1 and HD2 were imaged, there is only a $7.39\%$ probability they are lensed and that the hypothetical lensing galaxy was too faint to be detected. Meanwhile, with the limiting magnitudes of the Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) imaging capabilities, the probability drops precipitously to $0.058\%$ and $0.0012\%$, respectively. We further find it unlikely that the luminosities of both sources can be accounted for by lensing that produces a single, resolved image with sufficiently high magnification. Alternatively, in the unlikely low-probability event that their brightness does indeed result from lensing, there is a $30.9 \%$ probability that the lensing galaxy is too faint to be observable at the current limiting magnitude. Future HST (JWST) imaging will drop this probability to $0.245 \%$ ($0.0025 \%$). In summary, while deep imaging with HST and JWST is required to discard the lensing hypothesis entirely, it is unlikely that the exceptional luminosity of the two $z \sim 13$ sources can be accounted for by gravitational lensing.

We present new diagnostic metrics to probe the dynamical state of galaxy clusters. These novel metrics rely on the computation of the power spectra of the matter and gas distributions and their cross-correlation derived from cluster observations. This analysis permits us to cross-correlate the fluctuations in the matter distribution, inferred from high-resolution lensing mass maps derived from Hubble Space Telescope (HST) data, with those derived from the emitted X-ray surface brightness distribution of the hot Intra-Cluster medium (ICM) from the Chandra X-ray Observatory (CXO). These methodological tools allow us to quantify with unprecedented resolution the coherence with which the gas traces the mass and interrogate the assumption that the gas is in hydro-static equilibrium with the underlying gravitational potential. We characterize departures from equilibrium as a function of scale with a new gas-mass coherence parameter. The efficacy of these metrics is demonstrated by applying them to the analysis of two representative clusters known to be in different dynamical states: the massive merging cluster Abell 2744, from the HST Frontier Fields (HSTFF) sample, and the dynamically relaxed cluster Abell 383, from the Cluster Lensing And Supernova Survey with Hubble (CLASH) sample. Using lensing mass maps in combination with archival Chandra data, and simulated cluster analogs available from the OMEGA500 suite, we quantify the fluctuations in the mass and X-ray surface brightness and show that new insights into the dynamical state of the clusters can be obtained from our gas-mass coherence analysis.

Any successful model of dark matter must explain the diversity of observed Milky Way (MW) satellite density profiles, from very dense ultrafaints to large, low density satellites such as Crater~II that appear to be larger their anticipated host dark matter haloes. We consider cold dark matter (CDM), warm dark matter (WDM, 3.3keV thermal relic power spectrum), and a self-interacting dark matter model (SIDM) that induces gravothermal collapse in low mass subhaloes. Predictions for these density profiles are complicated by the limitations of simulation resolution in the stripping of subhaloes by the MW system, therefore we make predictions for satellite properties in these three models using $N$-body simulations combined with a semi-analytic halo stripping algorithm. We find that most CDM and WDM subhaloes of mass $>10^{8}$$M_{\odot}$ are large enough after stripping to fit most satellites; however, the required amount of stripping often requires a stronger tidal field than is available on the subhalo's orbit. The lower concentrations of WDM subhaloes enable this model to explain the required satellite masses with less stripping than is necessary for CDM, and is thus consistent with orbits of larger pericentres. SIDM cores offer the best fits to massive, low density satellites at the expense of predicting many $>10^{9}$$M_{\odot}$ subhaloes to host low density satellites with no observed analogue. We conclude that an SIDM model must have a very high velocity-dependent cross-section in order to match all satellites, and that WDM offers a marginally better fit than CDM to the MW satellite mass function.

The recent JWST discovery of a population of super-early (redshift $z> 10$), relatively massive (stellar mass $M_* = 10^{8-9} M_{\odot}$) and evolved (metallicity $Z \approx 0.1 Z_{\odot}$) galaxies, which nevertheless show blue ($\beta \simeq -2.6$) spectra, and very small dust attenuation ($A_{\rm V} \leq 0.02$), challenges our interpretation of these systems. To solve the puzzle we propose two solutions in which dust is either (a) ejected by radiation pressure, or (b) segregated with respect to UV-emitting regions. We clarify the conditions for which the two scenarios apply, and show that they can be discriminated by ALMA observations, such as the recent non-detection of the $88\mu m$ dust continuum in GHZ2 ($z\simeq 12$) favouring dust ejection.

We introduce the star formation and Supernova (SN) feedback model of the SATIN (Simulating AGNs Through ISM with Non-Equilibrium Effects) project to simulate the evolution of the star forming multi-phase interstellar medium (ISM) of entire disk galaxies. This galaxy-wide implementation of a successful ISM feedback model naturally covers an order of magnitude in gas surface density, shear and radial motions. It is implemented in the adaptive mesh refinement code RAMSES at a peak resolution of 9 pc. New stars are represented by star cluster (sink) particles with individual SN delay times for massive stars. With SN feedback, cooling and gravity, the galactic ISM develops a realistic three-phase structure. The star formation rates naturally follow observed scaling relations for the local Milky Way gas surface density. SNe drive additional turbulence in the warm (300 K < $T$ < 10$^4$ K) gas and increase the kinetic energy of the cold gas, cooling out of the warm phase. The majority of the gas leaving the galactic ISM is warm and hot with mass loading factors of $3 \le \eta \le 10$. While the hot gas is leaving the system, the warm and cold gas falls back onto the disc in a galactic fountain flow.

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $\Omega_{\rm m}$ and $\sigma_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $\Omega_{\rm m}$ and $\sigma_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $\Omega_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

Spherical stellar systems have weakly-damped response modes. The dipole modes are seiche modes. The quadrupole are zero pattern-speed prolate modes, the stable precursors to the radial-orbit instability (ROI). We demonstrate that small deviations in the distribution function (DF) can destabilise the dipole modes and describe the newly identified instabilities in NFW-like dark-matter (DM) halos and power-law spherical systems. The modes were identified in N-body simulations using multivariate singular spectrum analysis (MSSA) and corroborated using linear-response theory. The new mode peaks inside the half-mass radius but has a pattern speed typical of an outer-halo orbit. As it grows, the radial angle of the eccentric orbits that make up the mode correlate and lose angular momentum by a resonant couple to outer-halo orbits. This leads to an unsteady pattern with a density enhancement that swings from one side of the halo to another along a diameter, like the orbits that comprise the instability. In this way, the dipole mode is similar to the ROI. Since the DF found in Nature is unlikely to be smooth and isotropic with $df(E)/dE<0$ necessary for Antonov stability, these modes may be ubiquitous albeit slowly growing. Halos that are less extended than NFW, such as the Hernquist model, tend to be stable to this dipole instability. We present the critical stability exponents for one- and two-power models. These different critical outer power-law exponents illustrate that the gravitational coupling between the inner and outer DM halo depends on the global shape of density profile.

We present laboratory rotational spectroscopy of five isomers of cyanoindene (2-, 4-, 5-, 6-, and 7-cyanoindene) using a cavity Fourier-transform microwave spectrometer operating between 6-40 GHz. Based on these measurements, we report the detection of 2-cyanoindene (1H-indene-2-carbonitrile; 2-C$_9$H$_7$CN) in GOTHAM line survey observations of the dark molecular cloud TMC-1 using the Green Bank Telescope at centimeter wavelengths. Using a combination of Markov Chain Monte Carlo (MCMC), spectral stacking, and matched filtering techniques, we find evidence for the presence of this molecule at the 6.3$\sigma$ level. This provides the first direct observation of the ratio of a cyano-substituted polycyclic aromatic hydrocarbon (PAH) to its pure hydrocarbon counterpart, in this case indene, in the same source. We discuss the possible formation chemistry of this species, including why we have only detected one of the isomers in TMC-1. We then examine the overall hydrocarbon:CN-substituted ratio across this and other simpler species, as well as compare to those ratios predicted by astrochemical models. We conclude that while astrochemical models are not yet sufficiently accurate to reproduce absolute abundances of these species, they do a good job at predicting the ratios of hydrocarbon:CN-substituted species, further solidifying -CN tagged species as excellent proxies for their fully-symmetric counterparts.

Some white dwarfs accreting from non-degenerate companions show anomalous carbon and nitrogen abundances in the photospheres of their stellar components which have been postulated to be descendants of supersoft X-ray binaries. Therefore the carbon-to-nitrogen ratio can provide constraints upon their past evolution. We fit far ultraviolet spectroscopy of the cataclysmic variable HS0218+3229 taken with the Cosmic Origins Spectrograph using Markov Chain Monte Carlo. While some parameters depend upon the amount of reddening, the carbon-to-nitrogen ratio is about one tenth of the Solar value ($log C/N=-0.53^{+0.13}_{-0.14}$ and $-0.58^{+0.16}_{-0.15}$ for almost no reddening and E(B-V)=0.065, respectively, which are consistent within the uncertainties). We also provide estimates of the silicon and aluminum abundances, and upper limits for iron and oxygen. Using the measured parameters of HS0218+3229 we reconstruct its past using evolutionary simulations with MESA. We implemented Gaussian process fits to the MESA grid in order to determiner the most likely initial binary configuration of HS0218+3229. We found that an initial mass of the donor of $M_{\rm donor;i}=0.90-0.98,\mathrm{M}_{\odot}$ and an initial orbital period of $P_{\rm orb;i}=2.88$ days ($P_{\rm orb;i}=3.12-3.16$ days) for an assumed white dwarf mass of $M_{\mathrm{WD}}=0.83\,\mathrm{M}_{\odot}$ ($M_{\mathrm{WD}}=0.60\,\mathrm{M}_{\odot}$) are needed to replicate the measured parameters. These configurations imply that the system did not go through a phase of quasi-steady hydrogen-burning on the white dwarf's surface. However, it could have experienced a phase of thermal timescale mass transfer in the past if the initial mass ratio was $\geq1.5$. We predict that HS0218+3229 will evolve into a CV with a period below the $\simeq80$\,min period minimum for normal CVs, displaying helium and hydrogen in its spectrum.

We perform an analysis of the luminosity functions (LFs) of two types of ringed galaxies -- polar-ring galaxies and collisional ring galaxies -- using data from the Sloan Digital Sky Survey (SDSS). Both classes of galaxies were formed as a result of interaction with their environment and they are very rare objects. We constructed LFs of galaxies by different methods and found their approximations by the Schechter function. The luminosity functions of both types of galaxies show a systematic fall-off at low luminosities. The polar structures around bright ($M_r \leq -20^m$) and red ($g-r > +0.8$) galaxies are about twice as common as around blue ones. The LF of collisional rings is shifted towards brighter luminosities compared to polar-ring galaxies. We analysed the published data on the ringed galaxies in several deep fields and confirmed the increase in their volume density with redshift: up to z$\sim$1 their density grows as $(1+z)^m$, where $m \gtrsim 5$.

Using a sample of 92 Galactic Be stars, we compare inclination angles (the angle between a star's rotation axis and the line-of-sight) determined from H$\alpha$ emission line profile fitting to those determined by the spectroscopic signature of gravitational darkening. We find good agreement: 70% of the sample (64 out of 92 stars) is consistent with zero difference between the two methods using $1\sigma$ errors, and there is a strong linear correlation coefficient between the two methods of $r=+0.63\pm0.05$. There is some evidence that the H$\alpha$ profile fitting method overestimates the inclination angle for $i\lesssim 25^\circ$, perhaps due to the neglect of incoherent electron scattering on the H$\alpha$ line widths, while the gravitational darkening method underestimates the inclination angle for $i\gtrsim70^\circ$, perhaps due to the neglect of disk radiative transfer effects on the optical spectrum. Overall, it is demonstrated that a single H$\alpha$ spectrum of modest resolution and SNR can be used to extract a useful estimate for the inclination angle of an individual Be star. This allows equatorial rotation velocities for individual Be stars to be derived from $v \sin i $ measurements and will allow Be stars to be used to search for correlated spin axes in young, open clusters if unbiased (with respect to inclination) samples of Be stars are used.

Of the more than $3{,}000$ radio pulsars currently known, only ${\sim}300$ are in binary systems, and only five of these consist of young pulsars with massive non-degenerate companions. We present the discovery and initial timing, accomplished using the Canadian Hydrogen Intensity Mapping Experiment telescope (CHIME), of the sixth such binary pulsar, PSR J2108+4516, a $0.577$-s radio pulsar in a 269-day orbit of eccentricity 0.09 with a companion of minimum mass ${\lesssim}11$ M$_{\odot}$. Notably, the pulsar undergoes periods of substantial eclipse, disappearing from the CHIME $400{-}800$ MHz observing band for a large fraction of its orbit, and displays significant dispersion measure and scattering variations throughout its orbit, pointing to the possibility of a circumstellar disk or very dense stellar wind associated with the companion star. Sub-arcsecond resolution imaging with the Karl G. Jansky Very Large Array unambiguously demonstrates that the companion is a bright, $V \simeq 11$ OBe star, EM* UHA 138, located at a distance of $3.26(14)$ kpc. With further multi-wavelength followup, PSR J2108+4516 promises to serve as another rare laboratory for the exploration of companion winds, circumstellar disks, and short-term evolution through extended-body orbital dynamics.

We compare the two largest galaxy morphology catalogues, which separate early and late type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of $\sim$21 mag. The methodologies used for the construction of the catalogues include differences such as the cutout sizes, the labels used for training, and the input to the CNN - monochromatic images versus $gri$-band normalized images. In addition, one catalogue is trained using bright galaxies observed with DES ($i<18$), while the other is trained with bright galaxies ($r<17.5$) and `emulated' galaxies up to $r$-band magnitude $22.5$. Despite the different approaches, the agreement between the two catalogues is excellent up to $i<19$, demonstrating that CNN predictions are reliable for samples at least one magnitude fainter than the training sample limit. It also shows that morphological classifications based on monochromatic images are comparable to those based on $gri$-band images, at least in the bright regime. At fainter magnitudes, $i>19$, the overall agreement is good ($\sim$95\%), but is mostly driven by the large spiral fraction in the two catalogues. In contrast, the agreement within the elliptical population is not as good, especially at faint magnitudes. By studying the mismatched cases we are able to identify lenticular galaxies (at least up to $i<19$), which are difficult to distinguish using standard classification approaches. The synergy of both catalogues provides an unique opportunity to select a population of unusual galaxies.

We found a broad absorption line (BAL) outflow in the VLT/UVES spectrum of the quasar SDSS J235702.54-004824.0, in which we identified four subcomponents. We measured the column densities of the ions in one of the subcomponents ($v$ = -1600 km s$^{-1}$), which include O I and Fe II. We found the kinetic luminosity of this component to be at most ~2.4% of the quasar's Eddington luminosity. This is near the amount required to contribute to AGN feedback. We also examined the time-variability of a C IV mini-BAL found at $v$ = -8700 km s$^{-1}$, which shows a shallower and narrower absorption feature attached to it in previous SDSS observations from 2000 and 2001, but not in the spectra from 2005 and onwards.

We investigate the resolved kinematics of the molecular gas, as traced by ALMA in CO (2-1), of 25 cluster member galaxies across three different clusters at a redshift of $z\sim1.6$. This is the first large-scale analysis of the molecular gas kinematics of cluster galaxies at this redshift. By separately estimating the rotation curve of the approaching and receding side of each galaxy via kinematic modeling, we quantify the difference in total circular velocity to characterize the overall kinematic asymmetry of each galaxy. 3/14 of the galaxies in our sample that we are able to model have similar degrees of asymmetry as that observed in galaxies in the field at similar redshift. However, this leaved 11/14 galaxies in our sample with significantly higher asymmetry, and some of these galaxies have degrees of asymmetry of up to $\sim$50 times higher than field galaxies observed at similar redshift. Some of these extreme cases also have one-sided tail-like morphology seen in the molecular gas, supporting a scenario of tidal and/or ram pressure interaction. Such stark differences in the kinematic asymmetry in clusters versus the field suggest the evolutionary influence of dense environments, established as being a major driver of galaxy evolution at low-redshift, is also active in the high-redshift universe.

Much remains to be understood about the nature of exoplanets smaller than Neptune, most of which have been discovered in compact multi-planet systems. With its inner ultra-short period planet b aligned with the star and two larger outer planets d-c on polar orbits, the multi-planet system HD 3167 features a peculiar architecture and offers the possibility to investigate both dynamical and atmospheric evolution processes. To this purpose we combined multiple datasets of transit photometry and radial velocimetry (RV) to revise the properties of the system and inform models of its planets. This effort was spearheaded by CHEOPS observations of HD 3167b, which appear inconsistent with a purely rocky composition despite its extreme irradiation. Overall the precision on the planetary orbital periods are improved by an order of magnitude, and the uncertainties on the densities of the transiting planets b and c are decreased by a factor of 3. Internal structure and atmospheric simulations draw a contrasting picture between HD 3167d, likely a rocky super-Earth that lost its atmosphere through photo-evaporation, and HD 3167c, a mini-Neptune that kept a substantial primordial gaseous envelope. We detect a fourth, more massive planet on a larger orbit, likely coplanar with HD 3167d-c. Dynamical simulations indeed show that the outer planetary system d-c-e was tilted, as a whole, early in the system history, when HD 3167b was still dominated by the star influence and maintained its aligned orbit. RV data and direct imaging rule out that the companion that could be responsible for the present-day architecture is still bound to the HD\,3167 system. Similar global studies of multi-planet systems will tell how many share the peculiar properties of the HD3167 system, which remains a target of choice for follow-up observations and simulations.

Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computationally challenging. We introduce Spectroscopic r-Process Abundance Retrieval for Kilonovae (SPARK), a modular framework to perform Bayesian inference on kilonova spectra with the goals of inferring elemental abundance patterns and identifying absorption features at early times. SPARK inputs an atomic line list and abundance patterns from reaction network calculations into the TARDIS radiative transfer code. It then performs fast Bayesian inference on observed kilonova spectra by training a Gaussian process surrogate for the approximate posteriors of kilonova ejecta parameters, via active learning. We use the spectrum of GW170817 at 1.4 days to perform the first inference on a kilonova spectrum, and recover a complete abundance pattern. Our inference shows that this ejecta was generated by an r-process with either (1) high electron fraction Y_e ~ 0.35 and high entropy s/k_B ~ 25, or, (2) a more moderate Y_e ~ 0.30 and s/k_B ~ 14. These parameters are consistent with a shocked, polar dynamical component, and a viscously-driven outflow from a remnant accretion disk, respectively. We also recover previous identifications of strontium absorption at ~8000 AA, and tentatively identify yttrium and/or zirconium at < 4500 AA. Our approach will enable computationally-tractable inference on the spectra of future kilonovae discovered through multi-messenger observations.

We present the Distant Giants Survey, a three-year radial velocity (RV) campaign to measure the conditional occurrence of distant giant planets (M_p ~ 0.3 - 13 M_J, P ~100 days - 12 years) in systems hosting close-in small planets, P(DG|CS). For the past two years, we have monitored 47 Sun-like stars hosting small transiting planets (R_p < 10 R_E) detected by TESS. We present the selection criteria used to assemble our sample and report the discovery of two distant giant planets, TOI-1669 c and TOI-1694 c. For TOI-1669 c we find that Msin i = 0.602 +/- 0.078 M_J and P = 501 +/- 16 days, while for TOI-1694 c, Msin i = 0.929 +/- 0.037 M_J and P = 393.2 +/- 4.3 days. We also measured a true mass of M = 31.1 +/- 2.1 M_E for the 3.8-day transiting planet TOI-1694 b. Planets with mass ~30 M_E are rare at periods <~10 days, suggesting that TOI-1694 b may have reached its current location through dynamical interactions with its outer neighbor. At the end of the Distant Giants Survey, we will incorporate TOI-1669 c and TOI-1694 c into our calculation of P(DG|CS), a crucial statistic for understanding the relationship between outer giants and small inner companions.

We use AGN with $\rm L_{X} \sim 10^{42.5-44}\,erg\,s^{-1}$, from the COSMOS-Legacy survey that lie within the UltraVISTA region and cross match them with the LEGA-C catalogue. The latter provides measurements of the calcium break, D$_n$4000, and H$_\delta$ Balmer line that allow us to study the stellar populations of AGN and compare them with a galaxy reference catalogue. Our samples consist of 69 AGN and 2176 non-AGN systems, within $\rm 0.6<z<1.3$, that satisfy the same photometric selection criteria. We construct the SEDs of both population and use the CIGALE code to investigate the effect of the two indices in the SED fitting process. Our analysis shows that the inclusion of D$_n$4000 and H$_\delta$ allows CIGALE to constrain better the ages of the stellar populations. Furthermore, we find an increase of the estimated stellar masses by, on average, $\sim 0.2$ dex, in particular for systems with young stars (D$_n$4000$\,<1.5$), when the two indices are included in the SED fitting. We then compare the D$_n$4000 and H$_\delta$ of AGN with sources in the reference catalogue, accounting for the different stellar mass of the two populations. Low to moderate L$_X$ AGN tend to reside in galaxies with older stellar populations and are less likely to have experienced a recent star formation burst, compared to galaxies in the control sample. We, also, compare the two populations as a function of their morphology (bulge-dominated, BD, vs. non-BD) and compactness. A similar fraction of AGN and non-AGN systems are classified as non-BD ($\sim 70\%$). Our analysis shows that BD AGN tend to have younger stellar populations compared to BD, non-AGN systems. On the other hand, non-BD AGN have, on average, older stellar populations and are less likely to have experienced a burst compared to non-BD sources in the reference sample. Furthermore, AGN tend to prefer more compact systems compared to non-AGN.

The growing set of gravitational-wave sources is being used to measure the properties of the underlying astrophysical populations of compact objects, black holes and neutron stars. Most of the detected systems are black hole binaries. While much has been learned about black holes by analyzing the latest LIGO-Virgo-KAGRA (LVK) catalog, GWTC-3, a measurement of the astrophysical distribution of the black hole spin orientations remains elusive. This is usually probed by measuring the cosine of the tilt angle ($\cos\tau$) between each black hole spin and the orbital angular momentum, $\cos\tau=+1$ being perfect alignment. Abbott et al. (2021e) has modeled the $\cos\tau$ distribution as a mixture of an isotropic component and a Gaussian component with mean fixed at $+1$ and width measured from the data. In this paper, we want to verify if the data require the existence of such a peak at $\cos\tau=+1$. We use various models for the astrophysical tilt distribution and find that a) Augmenting the LVK model such that the mean of the Gaussian is not fixed at $+1$ returns results that strongly depend on priors. If we allow $\mu>+1$ then the resulting astrophysical $\cos\tau$ distribution peaks at $+1$ and looks linear, rather than Gaussian. If we constrain $-1\leq \mu\leq+1$ the Gaussian component peaks at $\mu=0.47^{+0.47}_{-1.04}$ (median and 90\% symmetric credible interval). Two other 2-component mixture models yield $\cos\tau$ distributions that either have a broad peak centered at $0.20^{+0.21}_{-0.18}$ or a plateau that spans the range $[-0.5, +1]$, without a clear peak at $+1$. b) All of the models we considered agree on the fact that there is \textit{no} excess of black hole tilts at around $-1$. c) While yielding quite different posteriors, the models considered in this work have Bayesian evidences that are the same within error bars.

We systematically search for quasi-periodic oscillatory (QPO) signals on the monthly time scale among the 1525 sources given in the Fermi Large Area Telescope Light Curve Repository (LCR). We find a transient QPO of 31.3$\pm$1.8 days in the gamma-ray band light curve of the TeV blazar S5 0716+714, which has 7 cycles (MJD 55918-56137) for the first time by Weighted wavelet Z-transform (WWZ) and Lomb-Scargle periodogram (LSP) methods. Monte Carlo simulations based on power spectral density (PSD) and probability distribution function (PDF) were used to evaluate the confidence level of the QPO, and the result is $\sim 4.1 \sigma$. Seasonal Autoregressive Integrated Moving Average (SARIMA) modeling of the light curve revealed it is a significant physical QPO. the physical models to explain the sporadic QPO of the month time scale in blazar were discussed. Our studies indicate that the helical jet model and blob move helically in a curved jet model to properly explain this kind of transient QPO.

The amplitude of the metagalactic ultraviolet background (UVB) at large-scales is impacted by two factors. First, it naturally attenuates at scales larger than mean-free-path of UVB photons due to the absorption by neutral intergalactic medium. Second, there are discrete and rare ionizing sources distributing in the Universe, emitting the UVB photons, and thus enhancing the local UVB amplitude. Therefore, for cosmological probe that is sensitive to the UVB amplitude and capable of detecting the large scale like Lyman-$\alpha$ forest spectrum, the fluctuation due to the clustering of ionizing sources becomes a significant factor for Lyman-$\alpha$ flux transmission and leave imprints on Lyman-$\alpha$ flux power spectrum at these large scales. In this work, we make use of a radiative transfer model that parametrizes the UVB source distribution by its bias $b_{\rm j}$ and shot noise $\overline{n}_{\rm j}$. We estimate the constraints on this model through the cross-correlation between Lyman-$\alpha$ forest survey and galaxy survey, using the DESI Lyman-$\alpha$ forest survey and the Roman Space Telescope emission line galaxy survey as an example. We show the detection sensitivity improvement for UVB parameters from disjoint to maximal overlap of DESI+Roman survey strategy. We also show that the degeneracy of two ionizing source parameters can be broken by increasing the overlapping survey area. Our results motivate survey strategies more dedicated to probe the UVB large-scale fluctuations.

Top-shaped asteroids have been observed among near-Earth asteroids. About half of them are reported to have moons (on the order of $\sim 1$wt.\% of the top-shaped primary) and many of them have an equatorial ridge. A recent study has shown that the enigmatic top-shaped figure of asteroids (e.g., Ryugu, Bennu, and Didymos) could result from an axisymmetric landslide of the primary during a fast spin-up near the breakup rotation period. Such a landslide would inevitably form a particulate disk around an asteroid with a short timescale ($\sim 3$ hours). However, the long-term full dynamical evolution is not investigated. Here, we perform a continuous simulation ($\sim 700$ hours) that investigates the sequence of events from the surface landslide that forms a top-shaped asteroid and a particulate disk to disk evolution. We show that the disk quickly spreads and produces moons (within $\sim 300$ hours). The mass of the formed moon is consistent with what is observed around the top-shaped asteroids. We also demonstrate that an equatorial ridge is naturally formed because a fraction of the disk particles re-accretes selectively onto the equatorial region of the primary. We envision that Ryugu and Bennu could once have an ancient moon that was later lost due to a successive moon's orbital evolution. Alternatively, at the top-shaped asteroid that has a moon, such as Didymos, no significant orbital evolution of the moon has occurred that would result in its loss. Our study would also be qualitatively applicable to any rubble-pile asteroids near the breakup rotation period.

In observational astronomy, noise obscures signals of interest. Large-scale astronomical surveys are growing in size and complexity, which will produce more data and increase the workload of data processing. Developing automated tools, such as convolutional neural networks (CNN), for denoising has become a promising area of research. We investigate the feasibility of CNN-based self-supervised learning algorithms (e.g., Noise2Noise) for denoising astronomical images. We experimented with Noise2Noise on simulated noisy astronomical data. We evaluate the results based on the accuracy of recovering flux and morphology. This algorithm can well recover the flux for Poisson noise ( $98.13${\raisebox{0.5ex}{\tiny$^{+0.77}_{-0.90} $}$\large\%$}) and for Gaussian noise when image data has a smooth signal profile ($96.45${\raisebox{0.5ex}{\tiny$^{+0.80}_{-0.96} $}$\large\%$}).

Studies and reconstructions of past solar activity require data on sunspots as well as faculae/plage and network. Such data are also important for understanding the magnetic activity and variability of the Sun and Sun-like stars. The longest available direct faculae/plage datasets are white-light facular and Ca II K observations going back to 1874 and 1892, respectively. Prior to that time, the only direct data available are for sunspots. We reassess the relationship between plage areas and sunspot records (areas and numbers) since 1892, to allow reconstructions of facular/plage areas which can be employed for studies going further back in time, i.e. over the period when solely sunspot observations are available. We use the plage areas derived from 38 consistently processed Ca II K archives as well as the plage area composite based on these archives. We find the relationship between plage and sunspot areas and numbers to be well represented by a power law function. We further find that the relationships depend on the bandwidth and the solar cycle strength. The reconstructions with a power law relationship are in good agreement with the original plage area series, whereas employment of a cycle-strength-dependent relationship improves the reconstructions only marginally. Performing the same analysis on other previously published plage area series, usually derived from a single archive with diverse processing techniques, returns different results when using different time series. This highlights the importance of applying a consistent processing to the various archives and demonstrates the uncertainties introduced by using previously published series. Our results have implications for past solar activity and irradiance reconstructions as well as for stellar activity studies, which sometimes assume a linear dependence between plage and sunspot areas.

Emission-line galaxies (ELGs) are crucial in understanding the formation and evolution of galaxies, while little is known about their variability. Here we report the study on the optical variability of a sample of ELGs selected in the COSMOS field, which has narrow-band observations in two epochs separated by $\gtrsim$ 12 years. This sample was observed with Suprime-Cam (SC) and Hyper Suprime-Cam (HSC) on the $Subaru$ telescope in NB816 and $i'/i$ bands, respectively. After carefully removing the wing effect of a narrow-band filter, we check the optical variability in a sample of 181 spectroscopically confirmed ELGs. We find that 0 (0/68) Ha emitters, 11.9% (5/42) [OIII] emitters, and 0 (0/71) [OII] emitters show significant variability ($|\Delta m_{NB}| \geq 3\,\sigma_{\Delta m_{NB,AGN}} = 0.20\, mag$) in the two-epoch narrow-band observations. We investigate the presence of active galactic nucleus (AGN) in this variable ELG (var-ELG) sample with three methods, including X-ray luminosity, mid-infrared activity, and radio-excess. We find zero bright AGN in this var-ELG sample, but cannot rule out the contribution from faint AGN. We find that SNe could also dominate the variability of the var-ELG sample. The merger morphology shown in the HST/F814W images of all the var-ELG sample is in agreement with the enhancement of star formation, i.e., the SNe activity.

SOXS is the new spectrograph for the ESO NTT telescope able to cover the optical and NIR bands thanks to two different arms: the UV-VIS (350-850 nm), and the NIR (800-2000 nm). In this article, we describe the final design of the visible camera cryostats, the test facilities for the CCD characterization, and the first results with the scientific detector. The UV-VIS detector system is based on a e2v CCD 44-82, a custom detector head coupled with the ESO Continuous Flowing Cryostat (CFC) cooling system and the New General Detector Controller (NGC) developed by ESO. The laboratory facility is based on an optical bench equipped with a Xenon lamp, filter wheels to select the wavelength, an integrating sphere, and a calibrated diode to measure the flux. This paper outlines the visible camera cryostat, the test facilities for the CCD characterization and the first results with the scientific detector in the laboratory and after the integration to the instrument.

Black hole X-ray binary (BHXB) GX 339-4 produces the most frequent X-ray outbursts (every 2-4 years) among known X-ray binary sources in our galaxy. Here we present the study of the evolution of GX 339-4 over the years of observations using the Hardness Intensity Diagram (HID), i.e., the q diagram. We present an analysis of two outbursts of the source observed by the Monitor of All-sky X-ray Image (MAXI) telescope aboard the International Space Station (ISS).

The next-generation CMB satellite missions are expected to provide robust constraints on a wide range of cosmological parameters with unprecedented precision. But these constraints on the parameters could weaken if we do not attribute dark energy to a cosmological constant. The cosmological models involving interaction between dark energy and dark matter can give rise to comparable energy densities at the present epoch, thereby alleviating the so-called cosmic coincidence problem. In the present paper, we perform a forecast analysis to test the ability of the future generation high-sensitive Cosmic Microwave Background (CMB), and Baryon Acoustic Oscillation (BAO) experiments to constrain phenomenological interacting dark energy models. We consider cosmic variance limited future CMB polarization experiment PICO along with BAO information from the DESI experiment to constrain the parameters of the interacting dark sector. Based on the stability of the cosmological perturbations, we consider two possibilities for the interaction scenario. We investigate the impact of both the coupling constant and the equation of state parameter of dark energy on the CMB temperature power spectrum, matter power spectrum, and $f\sigma_8$. We have used simulated temperature and polarization data from PICO within the multipole ranges ($\ell = 2 - 4000$), and as expected, we do see PICO alone produces better constraints than Planck on the $\Lambda$CDM parameters. With the integration of the PICO and DESI missions, we observe a significant improvement in the constraints on several cosmological parameters, especially the equation of state parameter of dark energy. However, we note that additional data is required to constrain a small positive coupling constant.

Sausage modes are one leading mechanism for interpreting short period quasi-periodic pulsations (QPPs) of solar flares. Forward modeling their radio emission is crucial for identifying sausage modes observationally and for understanding their connections with QPPs. Using the numerical output from three-dimensional magnetohydrodynamic (MHD) simulations, we forward model the gyrosynchrotron (GS) emission of flare loops modulated by sausage modes and examine the influence of loop fine structures. The temporal evolution of the emission intensity is analyzed for an oblique line of sight crossing the loop center. We find that the low- and high-frequency intensities oscillate in-phase at the period of sausage modes for models with or without fine structures. For low-frequency emissions where the optically thick regime arises, the modulation magnitude of the intensity is dramatically reduced by the fine structures at some viewing angles. On the contrary, for high-frequency emissions where the optically thin regime holds, the effect of fine structures or viewing angle is marginal. Our results show that the periodic intensity variations of sausage modes are not wiped out by the fine structures, and sausage modes remains a promising candidate mechanism for QPPs even when flare loops are fine-structured.

Recent CMB observations have resulted in very precise observational data. A robust and reliable CMB reconstruction technique can lead to efficient estimation of the cosmological parameters. We demonstrate the performance of our methodology using simulated temperature and polarization observations using cosmic variance limited future generation PRISM satellite mission. We generate samples from the joint distribution by implementing the CMB inverse covariance weighted internal-linear-combination (ILC) with the Gibbs sampling technique. We use the Python Sky Model (PySM), d4f1s1 to generate the realistic foreground templates. The Synchrotron is parametrized by a spatially varying spectral index, whereas thermal dust is described as two component dust model. We estimate the marginalized densities of CMB signal ${\bf S}$ and theoretical angular power spectrum $C_{\ell}$ utilizing the samples from the entire posterior distribution. The best-fit cleaned CMB map and the corresponding angular power spectrum are consistent with the CMB realization and the sky $C_{\ell}$ implying an efficient foreground minimized reconstruction. The likelihood function $P(C_{\ell}|{\bf D}$) estimated by making use of the Blackwell-Rao estimator is used for the estimation of the cosmological parameters. Our methodology can estimate tensor to scalar ratio $r\ge 0.0075$. Our current work demonstrates an analysis pipeline starting from the reliable estimation of CMB signal and its angular power spectrum to the case of cosmological parameter estimation using the foreground model independent Gibbs-ILC method.

We present development progress of the scheduler for the Son Of X-Shooter (SOXS) instrument at the ESO-NTT 3.58 meter telescope. SOXS will be a single object spectroscopic facility, consisting of a two-arms high-efficiency spectrograph covering the spectral range 350-2000 nanometer with a mean resolving power R$\approx$4500. SOXS will be uniquely dedicated to the UV-visible and near infrared follow up of astrophysical transients, with a very wide pool of targets available from the streaming services of wide-field telescopes, current and future. This instrument will serve a variety of scientific scopes in the astrophysical community, with each scope eliciting its specific requirements for observation planning, that the observing scheduler has to meet. Due to directions from the European Southern Observatory (ESO), the instrument will be operated only by La Silla staff, with no astronomer present on the mountain. This implies a new challenge for the scheduling process, requiring a fully automated algorithm that should be able to present the operator not only with and ordered list of optimal targets, but also with optimal back-ups, should anything in the observing conditions change. This imposes a fast-response capability to the scheduler, without compromising the optimization process, that ensures good quality of the observations. In this paper we present the current state of the scheduler, that is now almost complete, and of its web interface.

SOXS (SOn of X-Shooter) is a medium resolution (~4500) wide-band (0.35 - 2.0 {\mu}m) spectrograph which passed the Final Design Review in 2018. The instrument is in the final integration phase and it is planned to be installed at the NTT in La Silla by next year. It is mainly composed of five different optomechanical subsystems (Common Path, NIR spectrograph, UV-VIS spectrograph, Camera, and Calibration) and other mechanical subsystems (Interface flange, Platform, cable corotator, and cooling system). A brief overview of the optomechanical subsystems is presented here as more details can be found in the specific proceedings while a more comprehensive discussion is dedicated to the other mechanical subsystems and the tools needed for the integration of the instrument. Moreover, the results obtained during the acceptance of the various mechanical elements are presented together with the experiments performed to validate the functionality of the subsystems. Finally, the mechanical integration procedure is shown here, along with all the modifications applied to correct the typical problems happening in this phase.

We report the implemented architecture for monitoring the health and the quality of the Son Of X-Shooter (SOXS) spectrograph for the New Technology Telescope in La Silla at the European Southern Observatory. Briefly, we report on the innovative no-SQL database approach used for storing time-series data that best suits for automatically triggering alarm, and report high-quality graphs on the dashboard to be used by the operation support team. The system is designed to constantly and actively monitor the Key Performance Indicators (KPI) metrics, as much automatically as possible, reducing the overhead on the support and operation teams. Moreover, we will also detail about the interface designed to inject quality checks metrics from the automated SOXS Pipeline (Young et al. 2022).

Particle acceleration mechanisms in supermassive black hole jets, such as shock acceleration, magnetic reconnection, and turbulence, are expected to have observable signatures in the multi-wavelength polarization properties of blazars. The recent launch of the Imaging X-ray Polarimetry Explorer (IXPE) enables us, for the first time, to use polarization in the X-ray band (2-8 keV) to probe the properties of the jet synchrotron emission in high-frequency-peaked BL Lac objects (HSPs). We report the discovery of X-ray linear polarization (degree $\Pi_{\rm x}=15\pm$2\% and electric-vector position angle $\Psi_{\rm x}=35^\circ\pm4^\circ$) from the jet of the HSP Mrk~421 in an average X-ray flux state. At the same time, the degree of polarization at optical, infrared, and millimeter wavelengths was found to be lower by at least a factor of 3. During the IXPE pointing, the X-ray flux of the source increased by a factor of 2.2, while the polarization behavior was consistent with no variability. The higher level of $\Pi_{\rm x}$ compared to longer wavelengths, and the absence of significant polarization variability, suggest a shock as the most likely X-ray emission site in the jet of Mrk 421 during the observation. The multiwavelength polarization properties are consistent with an energy-stratified electron population, where the particles emitting at longer wavelengths are located farther from the acceleration site, where they experience a more disordered magnetic field.

We present the progresses of the simulation tools, the Exposure Time Calculator (ETC) and End-to-End simulator (E2E), for the Son Of X-Shooter (SOXS) instrument at the ESO-NTT 3.58-meter telescope. The SOXS will be a single object spectroscopic facility, made by a two-arms high-efficiency spectrograph, able to cover the spectral range 350-2000 nanometer with a mean resolving power R$\approx$4500. While the purpose of the ETC is the estimate, to the best possible accuracy, of the Signal-to-Noise ratio (SNR), the E2E model allows us to simulate the propagation of photons, starting from the scientific target of interest, up to the detectors. We detail the ETC and E2E architectures, computational models and functionalities. The interface of the E2E with external simulation modules and with the pipeline are described, too. Synthetic spectral formats, related to different seeing and observing conditions, and calibration frames to be ingested by the pipeline are also presented.

Local primordial non-Gaussianity (PNG) is a promising observable of the underlying physics of inflation, characterized by a parameter denoted by $f_{\rm NL}$. We present the methodology to measure local $f_{\rm NL}$ from the Dark Energy Survey (DES) data using the 2-point angular correlation function (ACF) via the induced scale-dependent bias. One of the main focuses of the work is the treatment of the integral constraint. This condition appears when estimating the mean number density of galaxies from the data and is especially relevant for PNG analyses, where it is found to be key in obtaining unbiased $f_{\rm NL}$ constraints. The methods are analysed for two types of simulations: $\sim 246$ GOLIAT N-body simulations with non-Gaussian initial conditions $f_{\rm NL}$ equal to -100 and 100, and 1952 Gaussian ICE-COLA mocks with $f_{\rm NL}=0$ that follow the DES angular and redshift distribution. We use the GOLIAT mocks to asses the impact of the integral constraint when measuring $f_{\rm NL}$. We obtain biased PNG constraints when ignoring the integral constraint, $f_{\rm NL} = -2.8\pm1.0$ for $f_{\rm NL}=100$ simulations, and $f_{\rm NL}=-10.3\pm1.5$ for $f_{\rm NL}=-100$ simulations, whereas we recover the fiducial values within $1\sigma$ when correcting for the integral constraint with $f_{\rm NL}=97.4\pm3.5$ and $f_{\rm NL}=-95.2\pm5.4$, respectively. We use the ICE-COLA mocks to validate our analysis in a DES-like setup, finding it to be robust to different analysis choices: best-fit estimator, the effect of integral constraint, the effect of BAO damping, the choice of covariance, and scale configuration. We forecast a measurement of $f_{\rm NL}$ within $\sigma(f_{\rm NL})=31$ when using the DES-Y3 BAO sample, with the ACF in the $1\ {\rm deg}<\theta<20\ {\rm deg}$ range.

Jellyfish galaxies are an intriguing snapshot of galaxies undergoing ram-pressure stripping (RPS) in dense environments, showing spectacular star-forming knots in their disks and tails. We study the ionized gas properties of five jellyfish galaxies in massive clusters with Gemini GMOS/IFU observations: MACSJ0916-JFG1 ($z=0.330$), MACSJ1752-JFG2 ($z=0.353$), A2744-F0083 ($z=0.303$), MACSJ1258-JFG1 ($z=0.342$), and MACSJ1720-JFG1 ($z=0.383$). BPT diagrams show that various mechanisms (star formation, AGN, or mixed effects) are ionizing gas in these galaxies. Radial velocity distributions of ionized gas seem to follow disk rotation of galaxies, with the appearance of a few high-velocity components in the tails as a sign of RPS. Mean gas velocity dispersion is lower than 50 \kms~in most star-forming regions except near AGNs or shock-heated regions, indicating that the ionized gas %in these star-forming regions is dynamically cold. Integrated star formation rates (SFRs) of these galaxies range from $7~{\rm M_{\odot}~{\rm yr^{-1}}}$ to $35~{\rm M_{\odot}~{\rm yr^{-1}}}$ and the tail SFRs are from $0.6~{\rm M_{\odot}~{\rm yr^{-1}}}$ to $16~{\rm M_{\odot}~{\rm yr^{-1}}}$, which are much higher than those of other jellyfish galaxies in the local universe. These high SFR values imply that RPS triggers intense star formation activity in these extreme jellyfish galaxies. The phase-space diagrams demonstrate that the jellyfish galaxies with higher stellar masses and higher host cluster velocity dispersion are likely to have more enhanced star formation activity. The jellyfish galaxies in this study have similar gas kinematics and dynamical states to those in the local universe, but they show a much higher SFR.

SOXS (Son Of X-Shooter) is a single object spectrograph offering a simultaneous spectral coverage from U- to H-band, built by an international consortium for the 3.58-m ESO New Technology Telescope at the La Silla Observatory. It is designed to observe all kind of transients and variable sources discovered by different surveys with a highly flexible schedule maintained by the consortium, based on the Target of Opportunity concept. SOXS is going to be a fundamental spectroscopic partner for any kind of imaging survey, becoming one of the premier transient follow-up instruments in the Southern hemisphere. This paper gives an updated status of the project, when the instrument is in the advanced phase of integration and testing in Europe, prior to the activities in Chile.

The Son-Of-XShooter (SOXS) is a single object spectrograph (UV-VIS & NIR) and acquisition camera scheduled to be mounted on the ESO 3.58-m New Technology Telescope at the La Silla Observatory. Although the underlying data reduction processes to convert raw detector data to fully-reduced science ready data are complex and multi-stepped, we have designed the SOXS Data Reduction pipeline with the core aims of providing end-users with a simple-to-use, well-documented command-line interface while also allowing the pipeline to be run in a fully automated state; streaming reduced data into the ESO Science Archive Facility without need for human intervention. To keep up with the stream of data coming from the instrument, there is the requirement to optimise the software to reduce each observation block of data well within the typical observation exposure time. The pipeline is written in Python 3 and has been built with an agile development philosophy that includes CI and adaptive planning.

Observations have revealed that the elemental abundances of carbon and oxygen in the warm molecular layers of some protoplanetary disks are depleted compared to those is the interstellar medium by a factor of ~10-100. Meanwhile, little is known about nitrogen. To investigate the time evolution of nitrogen, carbon, and oxygen elemental abundances in disks, we develop a one-dimensional model that incorporates dust settling, turbulent diffusion of dust and ices, as well as gas-ice chemistry including the chemistry driven by stellar UV/X-rays and the galactic cosmic rays. We find that gaseous CO in the warm molecular layer is converted to CO2 ice and locked up near the midplane via the combination of turbulent mixing (i.e., the vertical cold finger effect) and ice chemistry driven by stellar UV photons. On the other hand, gaseous N2, the main nitrogen reservoir in the warm molecular layer, is less processed by ice chemistry, and exists as it is. Then the nitrogen depletion occurs solely by the vertical cold finger effect of N2. As the binding energy of N2 is lower than that of CO and CO2, the degree of nitrogen depletion is smaller than that of carbon and oxygen depletion, leading to higher elemental abundance of nitrogen than that of carbon and oxygen. This evolution occurs within 1 Myr and proceeds further, when the $\alpha$ parameter for the diffusion coefficient is ~0.001. Consequently, the N2H+/CO column density ratio increases with time. How the vertical transport affects the midplane ice composition is briefly discussed.

The Son Of X-Shooter (SOXS) is a single object spectrograph offering simultaneous spectral coverage in UV-VIS (350-850 nm) and NIR (800-2000 nm) wavelength regimes with an average of R close to 4500 for a 1 slit. SOXS also has imaging capabilities in the visible wavelength regime. It is designed and optimized to observe all kinds of transients and variable sources. The final destination of SOXS is the Nasmyth platform of the ESO NTT at La Silla, Chile. The SOXS consortium has a relatively large geographic spread, and therefore the Assembly Integration and Verification (AIV) of this medium-class instrument follows a modular approach. Each of the five main sub-systems of SOXS, namely the Common Path, the Calibration Unit, the Acquisition Camera, the UV-VIS Spectrograph, and the NIR Spectrograph, are undergoing (or undergone) internal alignment and testing in the respective consortium institutes. INAF-Osservatorio Astronomico di Padova delivers the Common Path sub-system, the backbone of the entire instrument. We report the Common Path internal alignment starting from the assembly of the individual components to the final testing of the optical quality, and the efficiency of the complete sub-system.

The Son Of X-Shooter (SOXS) is a two-channel spectrograph along with imaging capabilities, characterized by a wide spectral coverage (350nm to 2000nm), designed for the NTT telescope at the La Silla Observatory. Its main scientific goal is the spectroscopic follow-up of transients and variable objects. The UV-VIS arm, of the Common Path sub-system, is characterized by the presence of a powered Atmospheric Dispersion Corrector composed (ADC) by two counter-rotating quadruplets, two prisms, and two lenses each. The presence of powered optics in both the optical groups represents an additional challenge in the alignment procedures. We present the characteristics of the ADC, the analysis after receiving the optics from the manufacturer, the emerging issues, the alignment strategies we followed, and the final results of the ADC in dispersion and optical quality.

The Son Of X-Shooter (SOXS) is a single object spectrograph, built by an international consortium for the 3.58-m ESO New Technology Telescope at the La Silla Observatory, ranging from 350 to 2000 nm. In this paper, we present the progress in the AIT phase of the Near InfraRed (NIR) arm. We describe the different AIT phases of the cryo, vacuum, opto-mechanics and detector subsystems, that finally converged at the INAF-OAB premises in Merate (Italy), where the NIR spectrograph is currently being assembled and tested, before the final assembly on SOXS.

SOXS (Son Of X-Shooter) is a single object spectrograph built by an international consortium for the ESO NTT telescope. SOXS is based on the heritage of the X-Shooter at the ESO-VLT with two arms (UV-VIS and NIR) working in parallel, with a Resolution-Slit product of about 4500, capable of simultaneously observing over the entire band the complete spectral range from the U- to the H-band. SOXS will carry out rapid and long-term Target of Opportunity requests on a variety of astronomical objects. The SOXS vacuum and cryogenic control system has been designed to evacuate, cool down and maintain the UV-VIS detector and the entire NIR spectrograph to their operating temperatures. The design chosen allows the two arms to be operated independently. This paper describes the final design of the cryo-vacuum control system, its functionalities and the tests performed in the integration laboratories.

SOXS (SOn of X-Shooter) is a high-efficiency spectrograph with a mean Resolution-Slit product of 3500 over the entire band capable of simultaneously observing the complete spectral range 350-2000 nm. It consists of three scientific arms (the UV-VIS Spectrograph, the NIR Spectrograph and the Acquisition Camera) connected by the Common Path system to the NTT, and the Calibration Unit. We present an overview of the flow from the scientific to the technical requirements, and the realization of the sub-systems. Further, we give an overview of the methodologies used for planning and managing the assembly of the sub-systems, their integration and tests before the acceptance of the instrument in Europe (PAE) along with the plan for the integration of SOXS to the NTT. SOXS could be used as an example for the system engineering of an instrument of moderate complexity, with a large geographic spread of the team.

Recent studies based on numerical models of the Local Group predict the existence of field haloes and galaxies that have visited two distinct galaxies in the past, called Hermeian haloes. This work presents an analysis of the Hermeian haloes population in two high-resolution dark matter-only N-body simulations from the MultiDark suit (the ESMDPL, VSMDPL). Hermeian haloes make up from 0.5 to 2.5 per cent of the total number of field haloes depending on their mass. Furthermore, the results of our study suggest that at a sufficiently high resolution simulation, Hermeian haloes may be found around almost every halo, making them interesting for studies of matter exchange between galaxies. We find that about half (22 out of 49) of the selected Local Group analogues contain Hermeian haloes that passed through the haloes of both the Milky Way and M31 if the distance between the two main haloes is below 1 $h^{-1} \;\mathrm{Mpc}$; this fraction drops to one-fifth (24 out of 108) for distances of up to 1.5 $h^{-1} \;\mathrm{Mpc}$. We confirm earlier findings that unlike other field haloes, Hermeians are grouped along the line connecting the primary hosts of Local Group-like systems, which should facilitate their identification in observations. The vast majority of the Hermeian haloes, whose second target is a Milky Way analogue, are currently moving away from it with increased velocity compared to remaining field halo populations. Interestingly, the obtained data admits that NGC 3109 could have passed through the Andromeda Galaxy and the Milky Way earlier.

We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_\phi$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_\phi$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_\phi$ in a similar manner, roughly preserving the shape of the linear $b_\phi(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_\phi$ in opposite directions. This induces significant changes to the $b_\phi(b_1)$ relation, with higher concentration halos having higher amplitude of $b_\phi(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_\phi \approx 3$, $b_\phi \approx -1$ and $b_\phi \approx 9$, respectively. Varying the halo concentration can make $b_\phi$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_\phi(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $\sigma_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_\phi(b_1)$ relation of halos and galaxies.

Binding energies (BEs) are one of the most important parameters for astrochemical modeling determining, because they govern whether a species stays in the gas-phase or is frozen on the grain surfaces. It is currently known that, in the denser and colder regions of the interstellar medium, sulphur is severely depleted in the gas phase. It has been suggested that it may be locked into the grain icy mantles. However, which are the main sulphur carriers is still a matter of debate. This work aims at establishing accurate BEs of 17 sulphur-containing species on two validated water ice structural models, the proton-ordered crystalline (010) surface and an amorphous water ice surface. We adopted Density Functional Theory (DFT)-based methods (the hybrid B3LYP-D3(BJ) and the hybrid meta-GGA M06-2X functionals) to predict structures and energetics of the adsorption complexes. London's dispersion interactions are shown to be crucial for an accurate estimate of the BEs due to the presence of the high polarizable sulphur element. While on the crystalline model the adsorption is restricted to a very limited number of binding sites with single valued BEs, on the amorphous model several adsorption structures are predicted, giving a BE distribution for each species. With the exception of few cases, both experimental and other computational data are in agreement with our calculated BE values. A final discussion on how useful the computed BEs are with respect to the snow lines of the same species in protoplanetary disks is provided

Several recent works have focused on the search for bright, high-z quasars (QSOs) in the South. Among them, the QUasars as BRIght beacons for Cosmology in the Southern hemisphere (QUBRICS) survey has now delivered hundreds of new spectroscopically confirmed QSOs selected by means of machine learning algorithms. Building upon the results obtained by introducing the probabilistic random forest (PRF) for the QUBRICS selection, we explore in this work the feasibility of training the algorithm on synthetic data to improve the completeness in the higher redshift bins. We also compare the performances of the algorithm if colours are used as primary features instead of magnitudes. We generate synthetic data based on a composite QSO spectral energy distribution. We first train the PRF to identify QSOs among stars and galaxies, then separate high-z quasar from low-z contaminants. We apply the algorithm on an updated dataset, based on SkyMapper DR3, combined with Gaia eDR3, 2MASS and WISE magnitudes. We find that employing colours as features slightly improves the results with respect to the algorithm trained on magnitude data. Adding synthetic data to the training set provides significantly better results with respect to the PRF trained only on spectroscopically confirmed QSOs. We estimate, on a testing dataset, a completeness of ~86% and a contamination of ~36%. Finally, 207 PRF-selected candidates were observed: 149/207 turned out to be genuine QSOs with z > 2.5, 41 with z < 2.5, 3 galaxies and 14 stars. The result confirms the ability of the PRF to select high-z quasars in large datasets.

Motivated by the measured velocity profile of the M87 jet using the KVN and VERA Array (KaVA) by Park et al. indicating that the starting position of the jet acceleration is farther from the central engine of the jet than predicted in general relativistic magnetohydrodynamic simulations, we explore how to mitigate the apparent discrepancy between the simulations and the KaVA observation. We use a semi-analytic jet model proposed by Tomimatsu and Takahashi. consistently solving the trans-magnetic field structure but neglecting any dissipation effects. By comparing the jet model with the observed M87 jet velocity profile, we find that the model can reproduce the logarithmic feature of the velocity profile, and can fit the observed data when choosing $c/(100r_{g}) \lesssim \Omega_{F} \lesssim c/(70r_{g})$ where $r_{g}$ is the gravitational radius. While a total specific energy (${\cal E}$) of the jet changes the terminal bulk Lorentz factor of the jet, a slower angular velocity of the black hole magnetosphere (funnel region) ($\Omega_{F}$) makes a light-cylinder radius ($r_{\rm lc}$) larger and it consequently pushes out a location of a starting point of the jet acceleration. Using the estimated $\Omega_{F}$ we further estimate the magnetic field strength on the event horizon scale in M87 by assuming Blandford-Znajek (BZ) process is in action. The corresponding magnetic flux threading the event horizon of M87 is in good agreement with a magnetically arrested disc (MAD) regime.

In this paper I present a method to enhance the search sensitivity for long transient Gravitational Waves produced by Neutron Star instabilities. This method consists in a selective image filter, called Triangular Filter, to be applied to data spectrograms. It is shown that thanks to this implementation a 20% gain in sensitivity is achievable.

CUBES (Cassegrain U-Band Efficient Spectrograph) is the recently approved high-efficiency VLT spectrograph aimed to observe the sky in the UV ground-based region (305-400 nm) with a high-resolution mode (~20K) and a low-resolution mode (~5K). In this paper we will briefly describe the requirements and the design of the several software packages involved in the project, namely the instrument control software, the exposure time calculator, the end-to-end simulator, and the data reduction software suite. We will discuss how the above mentioned blocks cooperate to build up a "software ecosystem" for the CUBES instrument, and to support the users from the proposal preparation to the science-grade data products.

Gamma-ray telescopes in space are bombarded by large fluxes of charged particles, photons and secondary neutrons. These particles and radiation pose a threat to the nominal operation of satellites and limit the detection sensitivity of gamma-ray instruments. The background noise generated in gamma-ray space detectors by impinging particles is always much higher than the astrophysical signal to be detected. In this chapter, we present the different types of orbits suitable for gamma-ray missions, discussing their advantages and disadvantages, as well as the value of experiments embarked in stratospheric balloons. We then review the physical properties of all the background components in the different orbits and the stratosphere.

We report a promising candidate of a low-luminosity active galactic nucleus (AGN) at z=5 that was selected from the first near-infrared images of the JWST CEERS project. This source, named CEERS-AGN-z5-1 at absolute 1450 A magnitude M1450=-19.5+/-0.3, was found via a visual selection of compact sources from a catalog of Lyman break galaxies at z>4, taking advantage of the superb spatial resolution of the JWST/NIRCam images. The 21 photometric data available from CFHT, HST, Spitzer, and JWST suggest that the continuum shape of this source is reminiscent of that for an unobscured AGN, and there is a clear color excess in the filters where the redshifted Hbeta+[O III] and Halpha are covered. The estimated line luminosity is L(Hbeta+[OIII]) =10^43.0 erg s^-1 and L(Halpha) =10^42.9 erg s^-1 with the corresponding rest-frame equivalent width EW(Hbeta+[OIII]) =1100 A and EW(Halpha) =1600 A, respectively. Our SED fitting analysis favors the scenario that this object is either a strong broad-line emitter or even a super-Eddington accreting black hole (BH) emerging in a metal-poor host galaxy. The bolometric luminosity, L_bol=2.5+/-0.3 x 10^44 erg s^-1, is consistent with those of z<0.35 broad-line AGNs with M_BH ~ 10^6 Msun accreting at the Eddington limit. This new AGN population at the first 1.1 billion years of the universe may close the gap between the observed BH mass range at high redshift and that of BH seeds. Spectroscopic confirmation is awaited to secure the redshift and its AGN nature.

Galaxies are thought to reside inside of large gravitationally bound structures of dark matter, so-called haloes. While the smallest of these haloes host no or only a few stars, the biggest host entire clusters of galaxies. Over cosmic history, haloes often collided and merged, forming bigger and bigger structures. Merger trees, i.e., catalogues of haloes evolving and connections between them as they grow and merge, have become a vital tool in describing and understanding the history of cosmological objects such as our Galaxy. Semi-analytical models, built on top of such merger trees, are a common approach for theoretical studies in cosmology. The semi-analytical nature of such models is especially beneficial when the dynamic range in spatial and time scales that need to be considered becomes too large for numerical simulations. Ancient Stars and Local Observables by Tracing Halos (A-SLOTH) is such a semi-analytical model and it is designed to simulate star formation in the early Universe in a fast and accessible way. It uses merger trees, either from numerical simulations or generated by statistical algorithms to describe the history of galaxies. The processes of baryonic physics, in particular gas cooling, star formation and stellar feedback are described with approximations and statistical models. The range of applications for this model is extensive and we, therefore, make it available to the scientific community.

The strong equivalence principle is violated by gravity theories of Milgromian dynamics (MOND) through the action of the external field effect. We test two different Lagrangian theories AQUAL and QUMOND based on their numerical solutions of the external field effect, by comparing two independent estimates of the mean external field strength of the nearby universe: a theory-deduced value from fitting the outer rotation curves of 114 galaxies and an empirical value from the large-scale distribution of cosmic baryons. The AQUAL-deduced external field strength from rotation curves agrees with that from the large-scale cosmic environment, while QUMOND-deduced value is somewhat higher. This suggests that AQUAL is likely to be preferred over QUMOND as an effective non-relativistic limit of a potential relativistic modified gravity theory.

The current meteorite taxonomy, a result of two centuries of meteorite research and tradition, entangles textural and genetic terms in a less than consistent fashion, with some taxa (like shergottites) representing varied lithologies from a single putative parent body while others (like pallasites) subsume texturally similar objects of multifarious solar system origins. The familiar concept of group as representative of one primary parent body is also difficult to define empirically. It is proposed that the classification becomes explicitly binominal throughout the meteorite spectrum, with classes referring to petrographically defined primary rock types, whereas groups retain a genetic meaning, but no longer tied to any assumption on the number of represented parent bodies. The classification of a meteorite would thus involve both a class and a group, in a two-dimensional fashion analogous to the way Van Schmus and Wood decoupled primary and secondary properties in chondrites. Since groups would not substantially differ, at first, from those in current use de facto, the taxonomic treatment of normal meteorites, whose class would bring no new information, would hardly change. Yet classes combined with high- or low-level groups would provide a standardized grid to characterize petrographically and/or isotopically unusual or anomalous meteorites (which make up the majority of represented meteorite parent bodies) for example, in relation to the carbonaceous/noncarbonaceous dichotomy. In the longer term, the mergers of genetically related groups, a more systematic treatment of lithology mixtures, and the chondrite/achondrite transition can further simplify the nomenclature.

The Sch\"{o}nberg-Chandrasekhar limit in post main sequence evolution for stars of masses in the range $1.4\lesssim M/M_{\odot}\lesssim 6$ gives the maximum pressure that the stellar core can withstand, once the central hydrogen is exhausted. It is usually expressed as a quadratic function of $1/\alpha$, with $\alpha$ being the ratio of the mean molecular weight of the core to that of the envelope. Here, we revisit this limit in scenarios where the pressure balance equation in the stellar interior may be modified, and in the presence of small stellar pressure anisotropy, that might arise due to several physical phenomena. Using numerical analysis, we derive a three parameter dependent master formula for the limit, and discuss various physical consequences. As a byproduct, in a limiting case of our formula, we find that in the standard Newtonian framework, the Sch\"{o}nberg-Chandrasekhar limit is best fitted by a polynomial that is linear, rather than quadratic, to lowest order in $1/\alpha$.

We present a new Planck CMB lensing-CMB temperature cross-correlation likelihood that can be used to constrain cosmology via the Integrated Sachs-Wolfe (ISW) effect. CMB lensing is an excellent tracer of ISW, and we use the latest PR4 Planck data maps and lensing reconstruction to produce the first public Planck likelihood to constrain this signal. We demonstrate the likelihood by constraining the CMB background temperature from Planck data alone, where the ISW-lensing cross-correlation is a powerful way to break the geometric degeneracy, substantially improving constraints from the CMB and lensing power spectra alone.

Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observations of NSs, combined with terrestrial laboratory constraints and further theoretical investigations, hold the promise to provide important insight into the properties of matter in a region of the quantum chromodynamics phase space that is otherwise inaccessible. This multidisciplinary endeavor imposes the following requirements for facilities and resources in the upcoming decade and beyond: * A next generation of gravitational wave detectors to uncover more double NS and neutron star-black hole mergers; * Sensitive radio telescopes to find the most massive and fastest spinning NSs; * Large-area, high-time-resolution and/or high angular resolution X-ray telescopes to constrain the NS mass-radius relation; * Suitable laboratory facilities for nuclear physics experiments to constrain the dense matter equation of state; * Funding resources for theoretical studies of matter in this regime; * The availability of modern large-scale high performance computing infrastructure. The same facilities and resources would also enable significant advances in other high-profile fields of inquiry in modern physics such as the nature of dark matter, alternative theories of gravity, nucleon superfluidity and superconductivity, as well as an array of astrophysics, including but not limited to stellar evolution, nucleosynthesis, and primordial black holes.

The precise estimation of the mass of galaxy clusters is a major issue for cosmology. Large galaxy cluster surveys rely on scaling laws that relate cluster observables to their masses. From the high resolution observations of ~ 45 galaxy clusters with NIKA2 and XMM-Newton instruments, the NIKA2 SZ Large Program should provide an accurate scaling relation between the thermal Sunyaev-Zel'dovich effect and the hydrostatic mass. In this paper, we present an exhaustive analysis of the hydrostatic mass of the well known galaxy cluster CL J1226.9+3332, the highest-redshift cluster in the NIKA2 SZ Large Program at z = 0.89. We combine the NIKA2 observations with thermal Sunyaev-Zel'dovich data from NIKA, Bolocam and MUSTANG instruments and XMM-Newton X-ray observations and test the impact of the systematic effects on the mass reconstruction. We conclude that slight differences in the shape of the mass profile can be crucial when defining the integrated mass at R500, which demonstrates the importance of the modeling in the mass determination. We prove the robustness of our hydrostatic mass estimates by showing the agreement with all the results found in the literature. Another key information for cosmology is the bias of the masses estimated assuming hydrostatic equilibrium hypothesis. Based on the lensing convergence maps from the Cluster Lensing And Supernova survey with Hubble (CLASH) data, we obtain the lensing mass estimate for CL J1226.9+3332. From this we are able to measure the hydrostatic-to-lensing mass bias for this cluster, that spans from 1 - bHSE/lens ~ 0.7 to 1, presenting the impact of data-sets and mass reconstruction models on the bias.

With a new generation of observatories coming online this decade, the process of characterizing exoplanet atmospheres will need to be reinvented. Currently mostly on the instrumental side, characterization bottlenecks will soon stand by the models used to translate spectra into atmospheric properties. Limitations stemming from our stellar and atmospheric models have already been highlighted. Here, we show that the current limitations of the opacity models used to decode exoplanet spectra propagate into an accuracy wall at ~0.5-1.0 dex (i.e., 3 to 10x) on the atmospheric properties, which is an order of magnitude above the precision targeted by JWST Cycle 1 programs and needed for, e.g., meaningful C/O-ratio constraints and biosignatures identification. We perform a sensitivity analysis using nine different opacity models and find that most of the retrievals produce harmonious fits owing to compensations in the form of >5$\sigma$ biases on the derived atmospheric parameters translating in the aforementioned accuracy wall. We suggest a two-tier approach to alleviate this problem involving a new retrieval procedure and guided improvements in opacity data, their standardization and optimal dissemination.

The pathfinder Dragonfly Spectral Line Mapper is a distributed aperture telescope based off of the Dragonfly Telephoto Array with additional instrumentation (the Dragonfly "Filter-Tilter") to enable ultranarrow bandpass imaging. The pathfinder is composed of three redundant optical tube assemblies (OTAs) which are mounted together to form a single field of view imaging telescope (where the effective aperture diameter increases as the square-root of the number of OTAs). The pathfinder has been on sky from March 2020 to October 2021 equipped with narrowband filters to provide proof-of-concept imaging, surface brightness limit measurements, on sky testing, and observing software development. Here we describe the pathfinder telescope and the sensitivity limits reached along with observing methods. We outline the current limiting factors for reaching ultra-low surface brightnesses and present a comprehensive comparison of instrument sensitivities to low surface brightness line emission and other methods of observing the ultra-faint line emission from diffuse gas. Finally, we touch on plans for the upcoming 120-OTA Dragonfly Spectral Line Mapper, which is currently under construction.

Telescope arrays allow high-performance wide-field imaging systems to be built more quickly and at lower cost than conventional telescopes. Distributed aperture telescopes (the premier example of which is the Dragonfly Telephoto Array) are a special type of array in which all telescopes point at roughly the same position in the sky. In this configuration the array performs like a large and optically very fast single telescope with unusually good control over systematic errors. In a few key areas, such as low surface brightness imaging over wide fields of view, distributed aperture telescopes outperform conventional survey telescopes by a wide margin. In these Proceedings we outline the rationale for distributed aperture telescopes, and highlight the strengths and weaknesses of the concept. Areas of observational parameter space in which the design excels are identified. These correspond to areas of astrophysics that are both relatively unexplored and which have unusually strong breakthrough potential.

The Dragonfly Spectral Line Mapper (DSLM) is the latest evolution of the Dragonfly Telephoto Array, which turns it into the world's most powerful wide-field spectral line imager. The DSLM will be the equivalent of a 1.6m aperture $f$/0.26 refractor with a built-in Integral Field Spectrometer, covering a five square degree field of view. The new telescope is designed to carry out ultra-narrow bandpass imaging of the low surface brightness universe with exquisite control over systematic errors, including real-time calibration of atmospheric variations in airglow. The key to Dragonfly's transformation is the "Filter-Tilter", a mechanical assembly which holds ultra-narrow bandpass interference filters in front of each lens in the array and tilts them to smoothly shift their central wavelength. Here we describe our development process based on rapid prototyping, iterative design, and mass production. This process has resulted in numerous improvements to the design of the DSLM from the initial pathfinder instrument, including changes to narrower bandpass filters and the addition of a suite of calibration filters for continuum light subtraction and sky line monitoring. Improvements have also been made to the electronics and hardware of the array, which improve tilting accuracy, rigidity and light baffling. Here we present laboratory and on-sky measurements from the deployment of the first bank of lenses in May 2022, and a progress report on the completion of the full array in early 2023.

We explore the possible application of linear covariance-based (LCB) filtering to line-intensity mapping (LIM) signal reconstructions. Originally introduced for reconstruction of the integrated Sachs-Wolfe effect in the cosmic microwave background, the LCB filter is an optimal map estimator that extends the simple Wiener filter by leveraging external correlated data. Given a detectable strong LIM-galaxy or LIM-LIM cross power spectrum, we show recovery of high-redshift, large-scale line-intensity fluctuations -- even in the presence of bright interloper emission -- in simulations of a futuristic [C II] LIM survey as well as simulated future iterations of the CO Mapping Array Project (COMAP). With sufficient galaxy abundances or low LIM survey noise, normalised cross-correlation between the LCB reconstruction and the true signal reaches 70-90% on large, linear comoving scales corresponding to $k\sim0.1$ Mpc$^{-1}$. This suggests the possible use of such signal reconstructions in astrophysical or cosmological contexts that require identifying the locations of line emissivity peaks and voids, although clear shortcomings exist on smaller scales. The successful application of the LCB filter in simulated LIM contexts highlights the importance of cross-correlations to studies of the reionising and reionised high-redshift universe with LIM and other large-scale structure surveys.

The origin of short gamma-ray bursts (sGRBs) is associated with outflows powered by the remnant of a binary neutron star merger. This remnant can be either a black hole or a highly magnetized, fastly spinning neutron star, also known as a magnetar. Here, we present the results of two relativistic magnetohydrodynamical (RMHD) simulations aimed at investigating the large-scale dynamics and propagation of magnetar collimated outflows through the environment surrounding the remnant. The first simulation evolves a realistic jet by injecting external simulation data, while the second evolves an analytical model jet with similar properties for comparison. We find that both outflows remain collimated and successfully emerge through the environment. However, they fail to attain relativistic velocities and only reach a mean maximum speed of ~0.7c for the realistic jet, and ~0.6c for the analytical jet. We also find that the realistic jet has a much more complex structure. The lack of highly relativistic speeds, that makes these jets unsuitable as short GRB sources, appears to be due to the specific injected properties and not general to all possible magnetar outflows.

We consider the gravity assist maneuver, that is, a correction of spacecraft motion at its passing near a planet, as a tool for evaluating the Eddington post-Newtonian parameters $\beta$ and $\gamma$, characterizing vacuum spherically symmetric gravitation fields in metric theories of gravity. We estimate the effect of variation in $\beta$ and $\gamma$ on a particular trajectory of a probe launched from the Earth's orbit and passing closely near Venus, where relativistic corrections slightly change the impact parameter of probe scattering in Venus's gravitational field. It is shown, in particular, that a change of $10^{-4}$ in $\beta$ or $\gamma$ leads to a shift of about 50 km in the probe's aphelion position.

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process.

The dynamical tide can play an important role in the orbital motion of close eccentric double white dwarf binaries. As the launching of the space-based gravitational-wave detector, the Laser Interferometer Space Antenna (LISA), is just around the corner, detection of gravitational wave signals from such systems is anticipated. In this paper, we discuss the influence of the dynamical tide on eccentric orbits, focusing on the effect on orbital precession. We show that in orbits with a high eccentricity, resonance can cause a large precession when a harmonic of the orbital frequency matches the natural frequencies of the normal modes of the star. In contrast to the case with circular orbits, each mode can encounter multiple resonances with different harmonics and these resonant regions can cover about 10% of the frequency space for orbits with close separations. In this case, the tidal precession effect is distinct from the other contributions and can be identified with LISA if the signal-to-noise ratio is high enough. However, within the highly eccentric-small separation region, the dynamical tide causes chaotic motion and the gravitational wave signal becomes unpredictable. Even not at resonance, the dynamical tide can contribute up to 20% of the precession for orbits close to Roche-lobe filling separation with low eccentricities and LISA can resolve these off-resonant dynamical tide effects within the low eccentricity-small orbital separation region of the parameter space. For lower mass systems, the dynamical tide effect can degenerate with the uncertainties of the eccentricity, making it unmeasurable from the precession rate alone. For higher mass systems, the radiation reaction effect becomes significant enough to constrain the eccentricity, allowing the measurement of the dynamical tide.

The Standard Model predicts a long-range force between fermions due to the exchange of a pair of neutrinos. However, this quantum force, proportional to $G_F^2/r^5$ in the massless neutrino limit, is feeble and has not yet been observed in experiments. In this paper, we compute this force, including the background effect caused by cosmic, reactor, and solar neutrinos. We find that the neutrino force can be significantly stronger in the presence of neutrino backgrounds. Remarkably, reactor neutrino background can strengthen the neutrino force by more than 20 orders of magnitude compared with the one in vacuum. We discuss the experimental prospects of detecting the neutrino force in neutrino backgrounds and find that the effect is close to the available sensitivity of the current fifth force experiments. The results are encouraging and a detailed experimental study is called for to check if the effect can be probed.

Earth rotation is determined by polar motion (PM) and length of day (lod). The excitation sources of PM are torques linked to fluid circulations ("geophysical excitations"), and those of lod to luni-solar tides ("astronomical excitations"). We explore the links between the rotations and revolutions of planets, following Lagrange's (1853) presentation of mechanics. The energy of a planet in motion in a central field is the sum of kinetic, centrifugal (planet dependent) and centripetal (identical for all planets) energies. For each planet, one can calculate a "constant of gravitation" Gp . For the giant planets, Gp decreases as a function of aphelia. There is no such organized behavior for the terrestrial planets. The perturbing potential of other planets generates a small angular contribution to the displacement : this happens to be identical to Einstein's famous formula for precession. Delays in the planet's perihelia follow a (-5/2) power law of a. This is readily understood in the Lagrange formalism (the centrifugal term takes over for small distances). The telluric planets have lost energy, probably transferred to the planets rotations. The ratio of areal velocities to rotation obeys a -5/2 power law of a. The ratio of areal velocity to integrated period R also fits a -5/2 power dependence, implying linearity of the energy exchange between revolution and rotation. For Einstein deformation of space-time by the Sun is the origin of the field perturbation. For Lagrange the perturbation could only be due to the interactions of torques. The perihelion delays, the areal velocities and the planetary rotations display power laws of aphelia, whose behavior contrasts with that of the kinetic moment. The areal velocity being linearly linked to the kinetic moment of planets, this must be the level at which the transfer is achieved.

A large class of scalar-tensor theories of gravity exhibit a screening mechanism that dynamically suppresses fifth forces in the Solar system and local laboratory experiments. Technically, at the scalar field equation level, this usually translates into nonlinearities which strongly limit the scope of analytical approaches. This article presents $femtoscope$ $-$ a Python numerical tool based on the Finite Element Method (FEM) and Newton method for solving Klein-Gordon-like equations that arise in particular in the symmetron or chameleon models. Regarding the latter, the scalar field behavior is generally only known infinitely far away from the its sources. We thus investigate existing and new FEM-based techniques for dealing with asymptotic boundary conditions on finite-memory computers, whose convergence are assessed. Finally, $femtoscope$ is showcased with a study of the chameleon fifth force in Earth orbit.

Interest in cosmological singularities has remarkably grown in recent times, particularly on future singularities with the discovery of late-time acceleration of the universe and dark energy. Recent work has seen a proper classification of such singularities into strong and weak based on their strength, with weak singularities being the likes of sudden, w and big freeze singularities and strong singularities like the big rip. This has led to a classification of such singularities in various types like Big rip is Type 1, w-singularity is type V etc. While singularities of type I-type IV have been discussed vividly by taking into account inhomogeneous equations of state (EOS), the same has not been attempted for type V singularities. So in this work we have discussed in detail about the formation of type V singularities in various cosmologies after considering inhomogeneous equations of state. We consider two inhomogeneous forms of the EOS in the context of four different cosmological backgrounds ; standard general relativistic cosmology, an asymptotically safe cosmology, a cosmology inspired by modified area-entropy relations, generalized uncertainty principles, holographic renormalization and Chern-Simons gravity( all of which can be coincidentally described by the same form of the modified Friedmann equation) and an f(R) gravity cosmology. We show in detail that one sees some very big differences in the occurence conditions of type V singularities when one makes such considerations. In the particular case of the f(R) gravity cosmology, we see that the type V singularities get completely removed. This work goes to show that the creation and formation of type V singularities is influenced most strongly by the form of the equation of state that one considers, way more so than what background cosmology one chooses.

This paper studies the instability of two-dimensional magnetohydrodynamic (MHD) systems on a sphere using analytical methods. The underlying flow consists of a zonal differential rotation and a toroidal magnetic field is present. Semicircle rules that prescribe the possible domain of the wave velocity in the complex plane for general flow and field profiles are derived. The paper then sets out an analytical study of the `clamshell instability', which features field lines on the two hemispheres tilting in opposite directions (Cally 2001, Sol. Phys. vol. 199, pp. 231--249). An asymptotic solution for the instability problem is derived for the limit of weak shear of the zonal flow, via the method of matched asymptotic expansions. It is shown that when the zonal flow is solid body rotation, there exists a neutral mode that tilts the magnetic field lines, referred to as the `tilting mode'. A weak shear of the zonal flow excites the critical layer of the tilting mode, which reverses the tilting direction to form the clamshell pattern and induces the instability. The asymptotic solution provides insights into properties of the instability for a range of flow and field profiles. A remarkable feature is that the magnetic field affects the instability only through its local behaviour in the critical layer.

Solar flare prediction is a central problem in space weather forecasting and has captivated the attention of a wide spectrum of researchers due to recent advances in both remote sensing as well as machine learning and deep learning approaches. The experimental findings based on both machine and deep learning models reveal significant performance improvements for task specific datasets. Along with building models, the practice of deploying such models to production environments under operational settings is a more complex and often time-consuming process which is often not addressed directly in research settings. We present a set of new heuristic approaches to train and deploy an operational solar flare prediction system for $\geq$M1.0-class flares with two prediction modes: full-disk and active region-based. In full-disk mode, predictions are performed on full-disk line-of-sight magnetograms using deep learning models whereas in active region-based models, predictions are issued for each active region individually using multivariate time series data instances. The outputs from individual active region forecasts and full-disk predictors are combined to a final full-disk prediction result with a meta-model. We utilized an equal weighted average ensemble of two base learners' flare probabilities as our baseline meta learner and improved the capabilities of our two base learners by training a logistic regression model. The major findings of this study are: (i) We successfully coupled two heterogeneous flare prediction models trained with different datasets and model architecture to predict a full-disk flare probability for next 24 hours, (ii) Our proposed ensembling model, i.e., logistic regression, improves on the predictive performance of two base learners and the baseline meta learner measured in terms of two widely used metrics True Skill Statistic (TSS) and Heidke Skill core (HSS), and (iii) Our result analysis suggests that the logistic regression-based ensemble (Meta-FP) improves on the full-disk model (base learner) by $\sim9\%$ in terms TSS and $\sim10\%$ in terms of HSS. Similarly, it improves on the AR-based model (base learner) by $\sim17\%$ and $\sim20\%$ in terms of TSS and HSS respectively. Finally, when compared to the baseline meta model, it improves on TSS by $\sim10\%$ and HSS by $\sim15\%$.

Formaldehyde (HCHO) is one one of the most important trace gas in the atmosphere, as it is a pollutant causing respiratory and other diseases. It is also a precursor of tropospheric ozone which damages crops and deteriorates human health. Study of HCHO chemistry and long-term monitoring using satellite data is important from the perspective of human health, food security and air pollution. Dynamic atmospheric chemistry models struggle to simulate atmospheric formaldehyde and often overestimate by up to two times relative to satellite observations and reanalysis. Spatial distribution of modelled HCHO also fail to match satellite observations. Here, we present deep learning approach using a simple super-resolution based convolutional neural network towards simulating fast and reliable atmospheric HCHO. Our approach is an indirect method of HCHO estimation without the need to chemical equations. We find that deep learning outperforms dynamical model simulations which involves complicated atmospheric chemistry representation. Causality establishing the nonlinear relationships of different variables to target formaldehyde is established in our approach by using a variety of precursors from meteorology and chemical reanalysis to target OMI AURA satellite based HCHO predictions. We choose South Asia for testing our implementation as it doesnt have in situ measurements of formaldehyde and there is a need for improved quality data over the region. Moreover, there are spatial and temporal data gaps in the satellite product which can be removed by trustworthy modelling of atmospheric formaldehyde. This study is a novel attempt using computer vision for trustworthy modelling of formaldehyde from remote sensing can lead to cascading societal benefits.

This report summarizes the findings of the CF1 Topical Subgroup to Snowmass 2021, which was focused on particle dark matter. One of the most important scientific goals of the next decade is to reveal the nature of dark matter (DM). To accomplish this goal, we must delve deep, to cover high priority targets including weakly-interacting massive particles (WIMPs), and search wide, to explore as much motivated DM parameter space as possible. A diverse, continuous portfolio of experiments at large, medium, and small scales that includes both direct and indirect detection techniques maximizes the probability of discovering particle DM. Detailed calibrations and modeling of signal and background processes are required to make a convincing discovery. In the event that a candidate particle is found through different means, for example at a particle collider, the program described in this report is also essential to show that it is consistent with the actual cosmological DM. The US has a leading role in both direct and indirect detection dark matter experiments -- to maintain this leading role, it is imperative to continue funding major experiments and support a robust R\&D program.

Previous studies from the astrophysics and laser physics communities have identified an interesting phenomenon wherein ultrarelativistic charged particles experiencing strong radiation reaction tend to move along special directions fixed by the local electromagnetic field. In the relativity literature these are known as the "principal null directions" (PNDs) of the Maxwell field. A particle in this regime has "Aristotelian" dynamics in the sense that its velocity (rather than acceleration) is determined by the local field. We study this Aristotelian equilibrium in detail, starting from the Landau-Lifshitz equation describing charged particle motion including radiation reaction. Using a Frenet-Serret frame adapted to the PNDs, we derive the Lorentz factor describing motion along the local PND, together with drift velocities reflecting slower passage from one PND to another. We derive conditions on the field configuration that are necessary for such an equilibrium to occur. We demonstrate agreement of our analytic formulas with full numerical solutions of the Landau-Lifshitz equation in the appropriate regime.

Since general relativity is the unique theory of massless spin 2 particles at large distances, the most reasonable way to have significant modifications is to introduce one or more light scalars that mediate a new long-range force. Most existing studies of such scalars invoke models that exhibit some kind of "screening" at short distances to hide the force from solar system tests. However, as is well known, such modifications also exhibit superluminality, which can be interpreted as a form of acausality. In this work we explore explicitly subluminal and causal scalar field models. In particular, we study a conformally coupled scalar $\phi$, with a small coupling to matter to obey solar system bounds, and a non-canonical kinetic term $K(X)$ ($X=(\partial\phi)^2/2$) that obeys all subluminality constraints and is hyperbolic. We consider $K(X)$ that is canonical for small $X$, but beyond some nonlinear scale enters a new scaling regime of power $p$, with $1/2<p<1$ (the DBI kinetic term is the limit $p=1/2$ and a canonical scalar is $p=1$). As opposed to screening (and superluminality), this new force becomes more and more important in the regime of high densities (and subluminality). We then turn to the densest environments to put bounds on this new interaction. We compute constraints from precession in binary systems such as Hulse-Taylor, we compute corrections to neutron star hydrostatic equilibrium, and we compute power in radiation, both tensor mode corrections and the new scalar mode, which can be important during mergers.