Papers for Wednesday, Nov 09 2022

The Square Kilometre Array Observatory (SKAO) is perhaps the most ambitious radio telescope envisaged yet. It will enable unprecedented studies of the Sun, the corona and the heliosphere and help to answer many of the outstanding questions in these areas. Its ability to make a vast previously unexplored phase space accessible, also promises a large discovery potential. The Indian solar and heliospheric physics community have been preparing for this science opportunity. A significant part of this effort has been towards playing a leading role in pursuing science with SKAO precursor instruments. This article briefly summarises the current status of the various aspects of work done as a part of this enterprise and our future goals.

Impulsive supernova feedback and non-standard dark matter models, such as self-interacting dark matter (SIDM), are the two main contenders for the role of the dominant core formation mechanism at the dwarf galaxy scale. Here we show that the impulsive supernova cycles that follow episodes of bursty star formation leave distinct features in the distribution function of stars: groups of stars with similar ages and metallicities develop overdense shells in phase space. If cores are formed through supernova feedback, we predict the presence of such features in star-forming dwarf galaxies with cored host halos. Their systematic absence would favor alternative dark matter models, such as SIDM, as the dominant core formation mechanism.

In recent years, deep learning has infiltrated every field it has touched, reducing the need for specialist knowledge and automating the process of knowledge discovery from data. This review argues that astronomy is no different, and that we are currently in the midst of a deep learning revolution that is transforming the way we do astronomy. We trace the history of astronomical connectionism from the early days of multilayer perceptrons, through the second wave of convolutional and recurrent neural networks, to the current third wave of self-supervised and unsupervised deep learning. We then predict that we will soon enter a fourth wave of astronomical connectionism, in which finetuned versions of an all-encompassing 'foundation' model will replace expertly crafted deep learning models. We argue that such a model can only be brought about through a symbiotic relationship between astronomy and connectionism, whereby astronomy provides high quality multimodal data to train the foundation model, and in turn the foundation model is used to advance astronomical research.

A complete census of Active Galactic Nuclei (AGN) is a prerequisite for understanding the growth of supermassive black holes across cosmic time. A significant challenge toward this goal is the whereabouts of heavily obscured AGN that remain uncertain. This paper sets new constraints on the demographics of this population by developing a methodology that combines X-ray spectral information with priors derived from multiwavelength observations. We select X-ray AGN in the Chandra COSMOS Legacy survey and fit their $2.2-500\mu m$ spectral energy distributions with galaxy and AGN templates to determine the mid-infrared ($6\mu m$) luminosity of the AGN component. Empirical correlations between X-ray and $6\mu m$ luminosities are then adopted to infer the intrinsic accretion luminosity at X-rays for individual AGN. This is used as prior information in our Bayesian X-ray spectral analysis to estimate physical properties, such as line-of-sight obscuration. Our approach breaks the degeneracies between accretion luminosity and obscuration that affect X-ray spectral analysis, particularly for the most heavily obscured (Compton-Thick) AGN with low photon counts X-ray spectra. The X-ray spectral results are then combined with the selection function of the Chandra COSMOS Legacy survey to derive the AGN space density and a Compton-Thick fraction of $21.0^{+16.1}_{-9.9}\%$ at redshifts $z<0.5$. At higher redshift, our analysis suggests upper limits to the Compton-Thick AGN fraction of $\le 40\%$. These estimates are at the low end of the range of values determined in the literature and underline the importance of multiwavelength approaches for tackling the challenge of heavily obscured AGN demographics.

NGC 5273 is a known optical and X-ray variable AGN. We analyze new and archival IR, optical, UV, and X-ray data in order to characterize its long-term variability from 2000 to 2022. At least one changing-look event occurred between 2011 and 2014, when the AGN changed from a Type 1.8/1.9 Seyfert to a Type 1. It then faded considerably at all wavelengths, followed by a dramatic but slow increase in UV/optical brightness between 2021 and 2022. We propose that NGC 5273 underwent multiple changing-look events between 2000 and 2022 -- starting as a Type 1.8/1.9, NGC 5273 changes-look to a Type 1 only temporarily in 2002 and again in 2014, reverting back to a Type 1.8/1.9 by 2005 and 2017, respectively. In 2022, it is again a Type 1 Seyfert with optical and NIR broad emission lines. We characterize the changing-look events and their connection to the dynamic accretion and radiative processes in NGC 5273.

The Sun emits copious amounts of photons and neutrinos in an approximately spatially isotropic distribution. Diffuse $\gamma$-rays and ultra-high energy (UHE) neutrinos from extragalactic sources may subsequently interact and annihilate with the emitted solar photons and neutrinos respectively. This will in turn induce an anisotropy in the cosmic ray background due to attenuation of the $\gamma$-ray and UHE neutrino flux by the solar radiation. Measuring this reduction, therefore, presents a simple and powerful astrophysical probe of electroweak interactions. In this letter we compute such anisotropies, which at the Earth (Sun) can be $\simeq 2\times 10^{-3}\,(0.5)\%$ and $\simeq 1\times 10^{-16}\,(2\times 10^{-14})\%$ for TeV scale $\gamma$-rays and PeV scale UHE neutrinos respectively. We briefly discuss exciting observational prospects for experiments such as the Fermi Gamma-Ray Space Telescope Large Area Telescope (Fermi LAT), High Energy Stereoscopic System (H.E.S.S), High-Altitude Water Cherenkov (HAWC) detector and IceCube. The potential for measuring $\gamma$-ray attenuation at orbital locations of other active satellites such as the Parker Solar Probe and James Webb Space Telescope is also explored.

We present early observations and analysis of the double-peaked Type IIb supernova (SN IIb) 2021zby. $TESS$ captured the prominent early shock cooling peak of SN 2021zby within the first $\sim$10 days after explosion with a 30-minute cadence. We present optical and near-infrared spectral series of SN 2021zby, including three spectra during the shock cooling phase. Using a multi-band model fit, we find that the inferred properties of its progenitor are consistent with a red supergiant or yellow supergiant, with an envelope mass of $\sim$0.3-3.0 M$_\odot$ and an envelope radius of $\sim$50-350$ R_\odot$. These inferred progenitor properties are similar to those of other SNe IIb with double-peak feature, such as SNe 1993J, 2011dh, 2016gkg and 2017jgh. This study further validates the importance of the high cadence and early coverage in resolving the shape of the shock cooling light curve, while the multi-band observations, especially UV, is also necessary to fully constrain the progenitor properties.

Inferring accurate posteriors for high-dimensional representations of the brightness of gravitationally-lensed sources is a major challenge, in part due to the difficulties of accurately quantifying the priors. Here, we report the use of a score-based model to encode the prior for the inference of undistorted images of background galaxies. This model is trained on a set of high-resolution images of undistorted galaxies. By adding the likelihood score to the prior score and using a reverse-time stochastic differential equation solver, we obtain samples from the posterior. Our method produces independent posterior samples and models the data almost down to the noise level. We show how the balance between the likelihood and the prior meet our expectations in an experiment with out-of-distribution data.

Minuscule primordial black holes before the end and after inflation can serve as "time capsules" bringing back energy from the past to a later epoch when they evaporate. As these black holes behave like matter, while the rest of the Universe content behaves like radiation, the mass fraction of these black holes, that is tiny at formation, becomes significant later. If sufficiently small, these black holes will evaporate while the Universe is still radiation dominated. We revisit this process and point out that gravitons produced during the evaporation behave as ``dark radiation". If the initial black holes are uniformly distributed so will be the gravitons and in this case they will be free of Silk damping and avoid current limits on "dark radiation" scenarios. Seeds for such black holes can arise during the last phases of inflation. We show here that with suitable parameters, this background graviton field can resolve the Hubble tension. We present current observational constraints on this scenario and suggest upcoming observational tests to prove or refute it. Finally, we also elaborate on the graviton background produced by particle annihilation during the Planck era or shortly after inflation.

We present the polarization profiles of 22 pulsars in the globular cluster 47 Tucanae using observations from the MeerKAT radio telescope at UHF-band (544-1088 MHz) and report precise values of dispersion measure (DM) and rotation measure (RM). We use these measurements to investigate the presence of turbulence in electron density and magnetic fields. The structure function of DM shows a break at $\sim 30$ arcsec ($\sim 0.6$ pc at the distance of 47 Tucanae) that suggests the presence of turbulence in the gas in the cluster driven by the motion of wind-shedding stars. On the other hand, the structure function of RM does not show evidence of a break. This non-detection could be explained either by the limited number of pulsars or by the effects of the intervening gas in the Galaxy along the line of sight. Future pulsar discoveries in the cluster could help confirm the presence and localise the turbulence.

The evolution of circumbinary discs and planets is often studied using two-dimensional (2D) numerical simulations, although recent work suggests that 3D effects may significantly alter the structure of the inner cavity created by the binary. In this study, we present the results of 3D hydrodynamical simulations of circumbinary discs that orbit around analogues of the Kepler-16 and Kepler-34 systems, including the effect of stellar heating and radiative cooling on the thermal disc structure. We find that compared to their 2D counterparts, the structures of the cavities in 3D circumbinary disc models appear to reach a quasi-stationary state more rapidly, and in a subset of our runs the evidence for this is unambiguous. Furthermore, the sizes and eccentricities of the inner cavity are smaller in 3D compared to 2D. We attribute this difference to enhanced spiral wave dissipation in disc regions above the midplane, where the cooling time is of the order of the dynamical timescale, resulting in smaller inner cavity sizes in 3D disc models. Our results suggest that migrating planets should park closer to the central binary in 3D models of circumbinary discs, and point to the importance of including the 3D structure when simulating circumbinary discs and planets.

The intermediate neutron capture process (i-process) operates at neutron densities between those of the slow and rapid neutron-capture processes. It can be triggered by the ingestion of protons in a convective helium-burning region. One possible astrophysical site is low-mass low-metallicity asymptotic giant branch (AGB) stars. We study here the possibility that actinides (particularly Th and U) may be significantly synthesized through i-process nucleosynthesis in AGB stars. We computed a 1 $M_{\odot}$ model at [Fe/H] $= -2.5$ with the stellar evolution code STAREVOL. We used a nuclear network of 1160 species from H to Cf coupled to the transport processes. During the proton ingestion event, the neutron density goes up to $\sim 10^{15}$ cm$^{-3}$. While most of the nuclear flow cycles in the neutron-rich Pb-Bi-Po region, a non-negligible fraction leaks towards heavier elements and eventually synthesizes actinides. The surface enrichment in Th and U is subject to nuclear and astrophysical model uncertainties that could be lowered in the future, in particular by a detailed analysis of the nuclear inputs that affect the neutron capture rates of neutron-rich isotopes between Pb and Pa. One stellar candidate that may confirm the production of actinides by the i-process is the carbon-enhanced metal-poor r/s star J0949-1617, which shows Th lines in its spectrum. Its surface abundance is shown to be reasonably well reproduced by our AGB model. Combined with cosmochronometry, this finding opens the way to dating the i-process event and thus obtaining a lower limit for the age of CEMP-r/s stars. Such a dating is expected to be accurate only if surface abundances of Th and U can be extracted simultaneously. This work shows that actinides can be synthesized in AGB stars through the i-process. As a consequence, the r-process may not be the sole mechanism for the production of U and Th.

Context. Sulfur (S) is one of the lesser-studied $\alpha$-elements. Published investigations of its behavior have so far focused on local stars, and only a few clusters of the Milky Way have been considered to study this topic. We aim to study the S content of the globular cluster Ruprecht 106 -- which has never before been studied for this purpose, but is known to present low levels of the [$\alpha$/Fe] abundance ratio -- and the open cluster Trumpler 5. The only star studied so far in Trumpler 5 shows an unexpectedly low abundance of S. Aims. With this work, we aim to provide the first S abundance in Ruprecht 106 and to investigate the S content of Trumpler 5 with a larger sample of stars. The open cluster Trumpler 20 is considered as a reference object. Methods. We performed a standard abundance analysis based on 1D model atmospheres in local thermodynamical equilibrium (LTE) and on high-resolution and high-signal-to-noise-ratio UVES-slit and UVES/FLAMES spectra. We also applied corrections for nonLTE. The metallicities of the targets were obtained by studying equivalent widths. Sulfur abundances were derived from multiplets 1, 6, and 8 by spectrosynthesis. Results. We find that the metallicities of Ruprecht 106 and Trumpler 5 are [Fe/H]= -1.37+/-0.11 and [Fe/H]= -0.49+/-0.14, respectively. Ruprecht 106 is less S-rich than the other Galactic clusters at similar metallicity. The low S content of Ruprecht 106, [S/Fe]NLTE= -0.52+/-0.13, is consistent with its shortage of $\alpha$-elements. This supports an extra-galactic origin of this cluster. We obtained a new and more robust S content value of Trumpler 5 of about [S/Fe]NLTE= 0.05+/-0.20. According to our results, Trumpler 5 follows the trend of the Galactic disk in the [S/Fe]LTE versus [Fe/H] diagram. Our results for Trumpler 20, of namely [Fe/H]= 0.06+/-0.15 and [S/Fe]NLTE= -0.28+/-0.21, are in agreement with those in the literature.

A model of geocoronal solar wind charge exchange (SWCX) emission was built and compared to five Suzaku detections of bright geocoronal SWCX events. An exospheric neutral hydrogen distribution model, charge exchange cross sections, solar wind ion data taken with the ACE and WIND satellites, and magnetic field models of the Earth's magnetosphere are all combined in order to predict time-variable geocoronal SWCX emission depending on line-of-sight directions of the Suzaku satellite. The modeled average intensities of O VII emission lines were consistent with the observed ones within a factor of three in four out of the five cases except for an event in which a line-of-sight direction was toward the night side of the high-latitude magnetosheath and a major geomagnetic storm was observed. Those of O VIII emission lines were underestimated by a factor of three or more in all the five cases. On the other hand, the modeled O VII and O VIII light curves reproduced the observed ones after being scaled by ratios between the observed and modeled average intensities. In particular, short-term variations due to line-of-sight directions traversing cusp regions during an orbital motion of the Suzaku satellite were reproduced. These results are discussed in the context of model uncertainties.

We aim to evaluate the behaviour of our 22 high-redshift (2.2 < z < 3.7) and high-luminosity (47.39 < Lbol < 48.36) quasars in the context of the 4-Dimensional Eigenvector 1. Our approach involves studying quasar physics through spectroscopic exploration of UV and optical emission line diagnostics. We are using new observations from ISAAC/VLT and mainly from the SDSS to cover the optical and the UV rest-frames, respectively. Emission lines are characterised both through a quantitative parametrisation of the line profiles, and by decomposing the emission line profiles using multicomponent fitting routines. We provide spectrophotometric properties and line profile measurements for Hb+[O iii], as well as for Si iv+O iv], C iv+He ii and the 1900 blend. Six out of the 22 objects present a significant blueshifted component on the Hb profile, and in 14/22 cases an Hb outflowing component associated to [O iii] is detected. The majority of [O iii] emission line profiles show blueshifted velocities larger than 250 km s^-1. [O iii] and C iv blueshifts show very high amplitudes and a high degree of correlation. Line width and shift are correlated for both [O iii] and C iv, suggesting that emission from outflowing gas is providing a substantial broadening to both lines. Otherwise, the links between C iv centroid velocity at half intensity (c(1/2)), Eddington ratio (L/LEdd), and bolometric luminosity are found to be in agreement with previous studies of high-luminosity quasars. Our analysis suggests that the behaviour of quasars of very high luminosity all along the main sequence is strongly affected by powerful outflows involving a broad range of spatial scales. The main sequence correlations remain valid at high redshift and high luminosity even if a systematic increase in line width is observed. Scaling laws based on UV Al iii and Hb emission lines are equally reliable estimators of MBH.

The outpouring of radiation during an X-ray burst can affect the properties of accretion discs around neutron stars: the corona can cool and collapse, the inner regions can be bled away due to enhanced accretion, and the additional heating will lead to changes in the disc height. In this paper, we investigate whether radiation from bursts can cause the disc to distort through a warping instability. Working in the limit of isotropic viscosity and linear growth, we find that bursts are more likely to drive disc warps when they have larger luminosities and longer durations. Therefore, warps will be most probable during intermediate duration bursts (IMDBs) and superbursts with evidence for photospheric radius expansion. Further, the development of warps depends on the disc viscosity with larger values of $\alpha$ increasing the likelihood of warp growth. We perform time-dependent evolution calculations of the development of warps during Type I bursts and IMDBs. Depending on the initial warp prior to the burst, we find the burst produces warps at $r$ <~ 50 $r_g$ that rapidly grow and decay on second-long timescales, or ones that grow more slowly and cover a large fraction of the disc. The pulsations of warp at small radii appear to have the properties needed to explain the achromatic fluctuations that have been observed during the tails of some IMDBs. The large scale, slowly growing warps could account for the large reflection strengths and absorbing column densities inferred late in the 4U 1820-30 and 4U 1636-53 superbursts.

We present the discovery of a gravitationally lensed dust-reddened QSO at z = 2.517, identified in a survey for QSOs by infrared selection. Hubble Space Telescope imaging reveals a quadruply lensed system in a cusp configuration, with a maximum image separation of ~1.8\arcsec. We find that compared to the central image of the cusp, the neighboring brightest image is anomalous by a factor of ~ 7 - 10, which is the largest flux anomaly measured to date in a lensed QSO. Incorporating high-resolution Jansky Very Large Array radio imaging and sub-mm imaging with the Atacama Large (sub-)Millimetre Array, we conclude that a low-mass perturber is the most likely explanation for the anomaly. The optical through near-infrared spectrum reveals that the QSO is moderately reddened with E(B - V) = 0.7 - 0.9. We see an upturn in the ultraviolet spectrum due to ~ 1% of the intrinsic emission being leaked back into the line of sight, which suggests that the reddening is intrinsic and not due to the lens. The QSO may have an Eddington ratio as high as L/L_Edd ~ 0.2. Consistent with previous red QSO samples, this source exhibits outflows in its spectrum as well as morphological properties suggestive of it being in a merger-driven transitional phase. We find a host-galaxy stellar mass of log M_*/M_Sun = 11.4, which is higher than the local M_BH vs. M_* relation, but consistent with other high redshift QSOs. When de-magnified, this QSO is at the knee of the luminosity function, allowing for the detailed study of a more typical moderate-luminosity infrared-selected QSO at high redshift.

This chapter aims at providing the tools and knowledge to understand and model the plasma environment surrounding comets in the innermost part near the nucleus. In particular, our goal is to give an updated post-Rosetta view of this ionised environment: what we knew, what we confirmed, what we overturned, and what we still do not understand.

The Colorado Ultraviolet Transit Experiment (CUTE) is a 6U NASA CubeSat carrying on-board a low-resolution, near-ultraviolet (2479-3306 A) spectrograph. It has a Cassegrain telescope with a rectangular primary to maximize the collecting area, given the shape of the satellite bus, and an aberration correcting grating to improve the image quality, and thus spectral resolution. CUTE, launched on the 27th of September 2021 to a Low Earth Orbit, is designed to monitor transiting extra-solar planets orbiting bright, nearby stars to improve our understanding of planet atmospheric escape and star-planet interaction processes. We present here the CUTE autONomous daTa ReductiOn pipeLine (CONTROL), developed for reducing CUTE data. The pipeline has been structured with a modular approach, which also considers scalability and adaptability to other missions carrying on-board a long-slit spectrograph. The CUTE data simulator has been used to generate synthetic observations used for developing and testing the pipeline functionalities. The pipeline has been tested and updated employing ight data obtained during commissioning and initial science operations of the mission.

Although weak lensing (WL) is a powerful method to estimate a galaxy cluster mass without any dynamical assumptions, a model bias can arise when the cluster density profile departs from the assumed model profile. In a merging system, the bias is expected to become most severe because the constituent halos undergo significant structural changes. In this study, we investigate WL mass bias in binary cluster mergers using a suite of idealized hydrodynamical simulations. Realistic WL shear catalogs are generated by matching the source galaxy properties, such as intrinsic shape dispersion, measurement noise, source densities, etc., to those from Subaru and {\it Hubble Space Telescope} observations. We find that, with the typical mass-concentration ($M$-$c$) relation and the Navarro-Frenk-White (NFW) profile, the halo mass bias depends on the time since the first pericenter passage and increases with the mass of the companion cluster. The time evolution of the mass bias is similar to that of the concentration, indicating that, to first order, the mass bias is modulated by the concentration change. For a collision between two $\sim10^{15}~M_{\odot}$ clusters, the maximum bias amounts to $\sim60\%$. This suggests that previous WL studies may have significantly overestimated the mass of the clusters in some of the most massive mergers. Finally, we apply our results to three merger cases: Abell 2034, MACS J1752.0+4440, and ZwCl 1856.8+6616, and report their mass biases at the observed epoch, as well as their pre-merger masses, utilizing their merger shock locations as tracers of the merger phases.

We present ALMA Band 7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866 $\mu$m, ruling out the possibility that S5-z17-1 is a low-$z$ dusty starburst with a star-formation rate (SFR) of $\gtrsim 30\,M_{\odot}$ yr$^{-1}$. We detect a 5.1$\sigma$ line feature at $338.726\pm0.007$ GHz exactly coinciding with the JWST source position, with a 2% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$\mu$m at $z=16.01$ or [CII]158$\mu$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11-13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11-17 is favored in every candidate as a very blue (UV continuum slope of $\approx-2.3$) and luminous ($M_{\rm UV}$ $\approx$ [$-$24:$-$21]) system. Still, we find in some candidates that reasonable SED fits ($\Delta$ $\chi^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$-5, where such populations predicted from their luminosity functions and EW([OIII]+H$\beta$) distributions are abundant in the survey volumes used for the $z\sim$11-17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in some of the $z\sim$11-17 candidates.

The Event Horizon Telescope (EHT) Collaboration has successfully produced images of two supermassive black holes, enabling novel tests of black holes and their accretion flows on horizon scales. The EHT has so far published total intensity and linear polarization images, while upcoming images may include circular polarization, rotation measure, and spectral index, each of which reveals different aspects of the plasma and space-time. The next-generation EHT (ngEHT) will greatly enhance these studies through wider recorded bandwidths and additional stations, leading to greater signal-to-noise, orders of magnitude improvement in dynamic range, multi-frequency observations, and horizon-scale movies. In this paper, we review how each of these different observables informs us about the underlying properties of the plasma and the spacetime, and we discuss why polarimetric studies are well-suited to measurements with sparse, long-baseline coverage.

While supermassive black hole masses have been cataloged across cosmic time, only a few dozen of them have robust spin measurements. By extending and improving the existing Event Horizon Telescope (EHT) array, the next-generation Event Horizon Telescope (ngEHT) will enable multifrequency, polarimetric movies on event horizon scales, which will place new constraints on the space-time and accretion flow. By combining this information, it is anticipated that the ngEHT may be able to measure tens of supermassive black hole masses and spins. In this white paper, we discuss existing spin measurements and many proposed techniques with which the ngEHT could potentially measure spins of target supermassive black holes. Spins measured by the ngEHT would represent a completely new sample of sources that, unlike pre-existing samples, would not be biased towards objects with high accretion rates. Such a sample would provide new insights into the accretion, feedback, and cosmic assembly of supermassive black holes.

Stellar evolutionary tracks for $0.12 \le M/M_\odot \le 1.0$ have been computed for each of several variations in the abundances of C, N, and O, assuming mass-fraction helium abundances $Y = 0.25$ and $0.29$, and 11 metallicities in the range $-1.5 \le$ [Fe/H] $\le -0.5$, in 0.2 dex increments. Such computations are provided for mixtures with [O/Fe] between $+0.4$ and $+0.8$, for different C:N:O ratios at a fixed value of [CNO/Fe], and for enhanced C. Computer codes are provided to interpolate within these grids to produce isochrones for ages $> 7$ Gyr and to generate magnitudes and colours for many broad-band filters using bolometric corrections based on MARCS model atmospheres and synthetic spectra. The models are compared with (i) similar computations produced by other workers, (ii) observed UV, optical, and IR colour-magnitude diagrams (CMDs), (iii) the effective temperatures, $(V-I_C)_0$ and $(V-K_S)_0$ colours of Pop. II stars in the solar neighbourhood, and (iv) empirical data for the absolute magnitude of the tip of the giant branch (TRGB). The isochrones are especially successful in reproducing the observed morphologies of observed CMDs and in satisfying the TRGB constraints. They also fare quite well in explaining the IR colours of low mass stars in globular clusters, indicating that they have [O/Fe] $\approx +0.6$, though some challenges remain.

We examine ultra-diffuse galaxies (UDGs) and their relation to non-UDGs in mass--radius--luminosity space. We begin by publishing Keck/KCWI spectroscopy for the Coma cluster UDG Y358, for which we measure both a recessional velocity and velocity dispersion. Our recessional velocity confirms association with the Coma cluster and Y358's status as a UDG. From our velocity dispersion (19 $\pm$ 3 km s$^{-1}$) we calculate a dynamical mass within the half-light radius which provides evidence for a core in Y358's dark matter halo. We compare this dynamical mass, along with those for globular cluster (GC)-rich/-poor UDGs in the literature, to mass profiles for isolated, gas-rich UDGs and UDGs in the NIHAO/FIRE simulations. We find GC-poor UDGs have dynamical masses similar to isolated, gas-rich UDGs, suggesting an evolutionary pathway may exist between the two. Conversely, GC-rich UDGs have dynamical masses too massive to be easily explained as the evolution of the isolated, gas-rich UDGs. The simulated UDGs match the dynamical masses of the GC-rich UDGs. However, once compared in stellar mass -- halo mass space, the FIRE/NIHAO simulated UDGs do not match the halo masses of either the isolated, gas-rich UDGs or the GC-rich UDGs at the same stellar mass. Finally, we supplement our data for Y358 with other UDGs that have measured velocity dispersions in the literature. We compare this sample to a wide range of non-UDGs in mass--radius--luminosity space, finding UDGs have a similar locus to non-UDGs of similar luminosity with the primary difference being their larger half-light radii.

We present the general relativistic electrodynamics and magnetohydrodynamics with a helically coupled scalar field. We consider three component system with the fluid, scalar field and electromagnetic fields with the helical coupling. We derive three exact formulations: the covariant formulation, the ADM formulation, and the fully nonlinear and exact perturbation formulation. We also derive the weak-gravity limit with fully relativistic fluid and fields. The latter two formulations are presented in the cosmological context.

There are increasing interests in binary supermassive black holes (SMBHs), but merging binaries with separations smaller than ~1 light days (~10^2 gravitational radii for 10^8 Msun), which are rapidly evolving under control of gravitational waves, are elusive in observations. In this paper, we discuss fates of mini-disks around component SMBHs for three regimes: 1) low rates (advection-dominated accretion flows: ADAFs); 2) intermediate rates; 3) super-Eddington accretion rates. Mini-disks with intermediate rates are undergoing evaporation through thermal conduction of hot corona forming a hybrid radial structure. When the binary orbital periods are shorter than sound propagation timescales of the evaporated mini-disks, a new instability, denoted as sound instability, arises because the disks will be highly twisted so that they are destroyed. We demonstrate a critical separation of A_{crit}~10^2 Rg from the sound instability of the mini-disks and the cavity is full of hot gas. For those binaries, component SMBHs are accreting with Bondi mode in the ADAF regime, showing periodic variations resulting from Doppler boosting effects in radio from the ADAFs due to orbital motion. In the mean while, the circumbinary disks (CBDs) are still not hot enough (ultraviolet deficit) to generate photons to ionize gas for broad emission lines. For slightly super-Eddington accretion of the CBDs, MgII line appears with decreases of UV deficit, and for intermediate super-Eddington Balmer lines appear, but CIV line never unless CBD accretion rates are extremely high. Moreover, if the CBDs are misaligned with the binary plane, it is then expected to have optical periodical variations with about ten times radio periods.

We present the largest compilation to date of optical observations during and following fast radio bursts (FRBs). The data set includes our dedicated simultaneous and follow-up observations, as well as serendipitous archival survey observations, for a sample of 15 well-localized FRBs, including 8 repeating and 7 one-off sources. Our simultaneous (and nearly simultaneous with a $0.4$-sec delay) optical observations of 7 (1) bursts from the repeating FRB 20220912A provide the deepest such limits to date for any extragalactic FRB, reaching a luminosity limit of $\nu L_\nu\lesssim 10^{42}$ erg s$^{-1}$ ($\lesssim 2\times10^{41}$ erg s$^{-1}$); these observations are also the deepest to date in terms of optical flux to radio fluence ratio of $f_{\rm opt}/F_{\rm radio}\lesssim 10^{-7}$ ms$^{-1}$ ($\lesssim 10^{-8}$ ms$^{-1}$), and place a limit on the flux ratio of $f_{\rm opt}/f_{\rm radio}\lesssim 0.02$ on a msec timescale or $\lesssim 2\times 10^{-5}$ ($\lesssim 10^{-6}$) on a sec timescale. These simultaneous limits provide useful constraints in the context of FRB emission models, such as the pulsar magnetosphere and pulsar nebula models. Interpreting all available optical limits in the context of the FRB synchrotron maser model, we find that they constrain the flare energies to $\lesssim 10^{43}-10^{47}$ erg (depending on the distances of the various repeating FRBs, with $\lesssim 10^{39}$ erg for SGR 1935+2154). These limits are generally at least an order of magnitude larger than the energies inferred from the FRBs themselves, although in the case of FRB 20220912A our simultaneous and rapid follow-up observations severely restrict the model parameter space. We conclude by exploring the potential of future rapid response and simultaneous observations with large optical telescopes.

The origin of star-formation in customarily passively evolving early-type massive galaxies is poorly understood. We present a case study of a massive galaxy, I~Zw~81, inside the Bootes void. The void galaxy is known to host an active galactic nuclei (AGN). Our detailed 2D decomposition of the surface brightness distribution in the $Canada$ $France$ $Hawaii$ $Telescope$ ($CFHT$) g- and r-bands revealed multiple structural components such as a nuclear point source, a bar, a ring, and an inner exponential disk followed by an outer low surface brightness (LSB) disk. I~Zw~81 turns out to be a disk-dominated galaxy with lenticular morphology. The modelling of the multi-wavelength spectral energy distribution (SED) shows that the galaxy is star-forming (SF), and belongs to the blue cloud. We find that the optical (g$-$r) color of the bar is bluer than the disks, and the far- and near-ultraviolet emission inside the galaxy observed with Ultraviolet Imaging Telescope (UVIT) onboard {\em AstroSat} is concentrated in the central few kpc region enclosing the bar. The strong bar might be playing a pivotal role in driving the gas inflow and causing SF activity in tandem with the minor merger-like interactions as evident from the deep $CFHT$ data. The low-luminosity AGN is insufficient to quench the central SF. The results are peculiar from the standpoint of a massive barred lenticular galaxy.

With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors.

On 9 October, 2022, an extremely powerful gamma-ray burst, GRB 221009A, was detected by several instruments. Despite being obstructed by the Milky Way galaxy, its afterglow outburst outshone all other GRBs seen before. LHAASO detected several thousands very-high energy photons extending up to 18 TeV. Detection of such energetic photons are unexpected due to the large opacity of the Universe. It is possible that in the afterglow epoch the intrinsic very high-energy photon flux from the source might have increased manifolds, which could compensate the attenuation by pair-production with the extragalactic background light. We propose such a scenario and show that very high-energy photons can be observed on the Earth from the interaction of very-high energy protons with the seed synchrotron photons in the external forward shock region of the GRB jet.

The missing baryon problem is one of the major unsolved problems in astronomy. Fast radio bursts (FRBs) are bright millisecond pulses with unknown origins. The dispersion measure of FRBs is defined as the electron column density along the line of sight, and accounts for every ionized baryon. Here we measure the baryon content of the Universe using 22 localized FRBs. Unlike previous works that fixed the value of dispersion measure of FRB host galaxies and ignored the inhomogeneities of the intergalactic medium (IGM), we use the probability distributions of dispersion measures contributed by host galaxies and IGM from the state-of-the-art IllustrisTNG simulations. We derive the cosmic baryon density of $\Omega_b=0.0490^{+0.0036}_{-0.0033}$ (1$\sigma$), with a precision of 7.0%. This value is dramatically consistent with other measurements, such as the cosmic microwave background and Big Bang nucleosynthesis. Our work supports that the baryons are not missing, but residing in the IGM.

Recent measurement of the primordial $^4$He abundance $Y_p$ from EMPRESS suggests a cosmological scenario with the effective number of neutrino species deviated from the standard value and non-zero lepton asymmetry. We argue that the extension of the standard cosmological model would be more demanded when the Hubble tension is taken into account, in which derived baryon density could be somewhat higher than the standard $\Lambda$CDM framework. We also discuss the issue by assuming early dark energy whose energy density can have a sizable fraction at the epoch of big bang nucleosynthesis. We show that the existence of early dark energy can reduce some tension implied by the EMPRESS $Y_p$ results.

Evolutionary tracks of stars with masses on the main sequence $0.84M_\odot\le M_\textrm{ZAMS}\le 0.95M_\odot$ and initial metal abundances $Z=0.006$ and $Z=0.01$ were computed under various assumptions about the mass loss rate at the red giant stage as well as at the AGB and the post--AGB stages. Among 160 evolutionary sequences we selected nearly 30 sequences where the final thermal flash of the helium shell source occurs in the early post--AGB stage when the ratio of the hydrogen envelope mass to the stellar mass ranges from 0.01 to 0.08. Selected evolutionary sequences were used for calculation of initial and inner boundary conditions used in solution of the equations of radiation hydrodynamics and turbulent convection describing evolution of stellar pulsations after the helium flash. Among about three dozen hydrodynamic models we found the three ones demonstrating almost eightfold increase of the pulsation period observed in V725 Sgr during the last century as well as the gradual transformation of fairly regular pulsations with period $\Pi\approx 12$ day to semi--regular non--linear oscillations with period $80~\textrm{day} \lesssim\Pi\lesssim 90~\textrm{day}$. We conclude that the anomalous growth of the pulsation period in V725 Sgr is due to the final thermal flash of the helium shell source that occurred in the early post--AGB star with mass $M\approx 0.53M_\odot$ and the mass of the hydrogen envelope ranging from $0.013M_\odot$ to $0.019M_\odot$.

The Edelweiss collaboration performs light Dark Matter (DM) particles searches with germanium bolometer collecting charge and phonon signals. Thanks to the Neganov-Trofimov-Luke (NTL) effect, a RMS resolution of 4.46 electron-hole pairs was obtained on a massive (200g) germanium detector instrumented with a NbSi Transition Edge Sensor (TES) operated underground at the Laboratoire Souterrain de Modane (LSM). This sensitivity made possible a search for WIMP using the Migdal effect down to 32 MeV/C$^{2}$ and exclude cross-sections down to 10$^{-29}$ cm$^2$. It is the first measurement in cryogenic germanium with such thermal sensor, proving the high relevance of this technology. Furthermore, such TES have shown sensitivity to out-of-equilibrium phonons, paving the way for EDELWEISS new experience CRYOSEL. This is an important step in the development of Ge detectors with improved performance in the context of the EDELWEISS-SubGeV program.

Context: Gas is now detected in many extrasolar systems around mature stars aged between 10 Myr to $\sim$ 1 Gyr with planetesimal belts. Gas in these mature disks is thought to be released from planetesimals and has been modelled using a viscous disk approach. At low densities, this may not be a good assumption as the gas could be blown out by the stellar wind instead. Methods: We developed an analytical model for A to M stars that can follow the evolution of gas outflows and target when the transition occurs between a disk or a wind. The crucial criterion is the gas density for which gas particles stop being protected from stellar wind protons impacting at high velocities on radial trajectories. Results: We find that: 1) Belts of radial width $\Delta R$ with gas densities $< 7 \, (\Delta R/50 {\rm \, au})^{-1}$ cm$^{-3}$ would create a wind rather than a disk, which would explain the recent outflowing gas detection in NO Lup. 2) The properties of this belt wind can be used to measure stellar wind properties such as their densities and velocities. 3) Debris disks with low fractional luminosities $f$ are more likely to create gas winds, which could be observed with current facilities. Conclusions: The systems containing low gas masses such as Fomalhaut or TWA 7 or more generally, debris disks with fractional luminosities $f \lesssim 10^{-5} (L_\star/L_\odot)^{-0.37} $ or stellar luminosity $\gtrsim 20 \, L_\odot$ (A0V or earlier) would rather create gas outflows (or belt winds) than gas disks. Gas observed to be outflowing at high velocity in the young system NO Lup could be an example of such belt winds. The detection of these gas winds is possible with ALMA (CO and CO$^+$ could be good wind tracers) and would allow us to constrain the stellar wind properties of main-sequence stars, which are otherwise difficult to measure.

Sgr A East is the supernova remnant closest to the centre of the Milky Way. Its age has been estimated to be either very young, around 1-2 kyr, or about 10 kyr, and its exact origin remains unclear. We aspire to create a simple model of a supernova explosion that reproduces the shape, size, and location of Sgr A East. Using a simplified hydrodynamical code, we simulated the evolution of a supernova remnant in the medium around the Galactic centre. The latter consists of a nearby massive molecular cloud with which Sgr A East is known to be interacting and a wind from the nuclear star cluster. Our preferred models of the Sgr A East remnant are compatible with an age of around 10 kyr. We also find suitable solutions for older ages, but not for ages younger than 5 kyr. Our simulations predict that the supernova exploded at a distance of about 3.5 pc from the Galactic centre, below the Galactic plane, slightly eastwards from the centre and 3 pc behind it.

We present a catalogue of 4499 groups and clusters of galaxies from the first data release of the multi-filter (5 broad, 7 narrow) Southern Photometric Local Universe Survey (S-PLUS). These groups and clusters are distributed over 273 deg$^2$ in the Stripe 82 region. They are found using the PzWav algorithm, which identifies peaks in galaxy density maps that have been smoothed by a cluster scale difference-of-Gaussians kernel to isolate clusters and groups. Using a simulation-based mock catalogue, we estimate the purity and completeness of cluster detections: at $\mathrm{S/N>3.3}$ we define a catalogue that is 80% pure and complete in the redshift range $\mathrm{0.1<z<0.4}$, for clusters with $\mathrm{M_{200} > 10^{14}}$ M$_\odot$. We also assessed the accuracy of the catalogue in terms of central positions and redshifts, finding scatter of $\mathrm{\sigma_R=12}$ kpc and $\mathrm{\sigma_z=8.8 \times 10^{-3}}$, respectively. Moreover, less than 1% of the sample suffers from fragmentation or overmerging. The S-PLUS cluster catalogue recovers $\sim$80% of all known X-ray and Sunyaev-Zel'dovich selected clusters in this field. This fraction is very close to the estimated completeness, thus validating the mock data analysis and paving an efficient way to find new groups and clusters of galaxies using data from the ongoing S-PLUS project. When complete, S-PLUS will have surveyed 9300 deg$^{2}$ of the sky, representing the widest uninterrupted areas with narrow-through-broad multi-band photometry for cluster follow-up studies.

XDCP0044.0-2033 is the most massive galaxy cluster known at z>1.5 and its core shows a high density of galaxies which are experiencing mergers and hosting nuclear activity. We present a multi-wavelength study of a region located 157 kpc from the center of this galaxy cluster, for which we have photometric and spectroscopic multi-wavelength observations (high resolution HST images in F105W, F140W and F160W bands, NIR KMOS data in H and YJ bands and Chandra ACIS-S X-ray data). Our main goal is to investigate the environmental effects acting on the galaxies inhabiting this high density region. We find that the analyzed region hosts at least nine different sources, six of them confirmed to be cluster members within a narrow redshift range 1.5728<z<1.5762. These sources form two different complexes at a projected distance of $\sim$13 kpc, which are undergoing merging on an estimated timescale off 10-30 Myr. One of the sources shows the presence of a broad H alpha emission line and is classified as Type 1 AGN. This AGN is associated to an X-ray point-like source, whose emission appears moderately obscured (with intrinsic absorption $N_{H} \sim 10^{22} cm^{-2}$) and hosts a relatively massive black hole with mass $M_{BH} \sim 10^{7} M_{\odot}$, which is accreting with an Eddington ratio of $\sim$0.2. We conclude that the analyzed region is consistent with being the formation site of a secondary BCG. These findings, together with an in-depth analysis the X-ray morphology of the cluster, suggest a merging scenario of the entire cluster, with two massive halos both harbouring two rapidly evolving BCGs on the verge of being assembled. Our results are also consistent with the scenario in which the AGN phase in member galaxies is triggered by gas-rich mergers, playing a relevant role in the formation of the red sequence of elliptical galaxies observed in the center of local galaxy clusters.

Context. Radial velocity (RV) measurements induced by the presence of planets around late-type stars are contaminated by stellar signals that are of the order of a few meters per second in amplitude, even for the quietest stars. Those signals are induced by acoustic oscillations, convective granulation patterns, active regions co-rotating with the stellar surface, and magnetic activity cycles. Aims. This study investigates the properties of all coherent stellar signals seen on the Sun on timescales up to its sidereal rotational period. By combining HARPS and HARPS-N solar data spanning several years, we are able to clearly resolve signals on timescales from minutes to several months. Methods. We use a Markov Chain Monte Carlo (MCMC) mixture model to determine the quality of the solar data based on the expected airmass-magnitude extinction law. We then fit the velocity power spectrum of the cleaned and heliocentric RVs with all known variability sources, to recreate the RV contribution of each component. Results. After rejecting variations caused by poor weather conditions, we are able to improve the average intra-day root mean square (RMS) value by a factor of ~1.8. On sub-rotational timescales, we are able to fully recreate the observed RMS of the RV variations. In order to also include rotational components and their strong alias peaks introduced by nightly sampling gaps, the alias powers are accounted for by being redistributed to the central frequencies of the rotational harmonics. Conclusions. In order to enable a better understanding and mitigation of stellar activity sources, their respective impact on the total RV must be well-measured and characterized. We are able to recreate RV components up to rotational timescales, which can be further used to analyse the impact of each individual source of stellar signals on the detectability of exoplanets.

We study the performance of three pairs of tilted $\Lambda$CDM cosmological models, two pairs allowing for non-flat spatial hypersurfaces with CMB temperature and polarization power spectrum data (P18) from Planck, P18 lensing (lensing), and non-CMB data (non-CMB). For the six models, we measure cosmological parameters and study whether or not pairs of the data sets are mutually consistent in these models. Half of these models allow the lensing consistency parameter $A_L$ to be an additional free parameter, while the other three have $A_L = 1$. The tilted spatially-flat models assume the usual primordial spatial inhomogeneity power spectrum. The tilted non-flat models assume either the primordial power spectrum used in the Planck group analyses [Planck $P(q)$] or a recently computed power spectrum [new $P(q)$]. In the tilted non-flat models with $A_L=1$ we find differences between P18 data and non-CMB data cosmological parameter constraints, which are large enough to rule out the Planck $P(q)$ model at 3$\sigma$ but not the new $P(q)$ model. While both P18 data and non-CMB data separately favor a closed geometry when P18+non-CMB data are jointly analyzed the evidence in favor of non-flat hypersurfaces subsides. Differences between P18 data and non-CMB data cosmological constraints subside when $A_L$ is allowed to vary. From the most restrictive P18+lensing+non-CMB data combination we get almost model-independent constraints and find that the $A_L>1$ option is preferred over the $\Omega_k<0$ one, with the $A_L$ parameter, for all models, being larger than unity by $\sim 2.5\sigma$. According to the deviance information criterion, in the P18+lensing+non-CMB analysis, the varying $A_L$ option is on the verge of being strongly favored over the $A_L=1$ one, which could indicate a problem for the standard tilted flat $\Lambda$CDM model (Abridged abstract).

Statistical inference of population parameters of astrophysical sources is challenging. It requires accounting for selection effects, which stem from the artificial separation between bright detected and dim undetected sources that is introduced by the analysis pipeline itself. We show that these effects can be modeled self-consistently in the context of sequential simulation-based inference. Our approach couples source detection and catalog-based inference in a principled framework that derives from the truncated marginal neural ratio estimation (TMNRE) algorithm. It relies on the realization that detection can be interpreted as prior truncation. We outline the algorithm, and show first promising results.

During the last decades the industry has seen the number of Earth orbiting satellites rise, mostly due to the need to monitor Earth as well as to establish global communication networks. Nano, micro, and small satellites have been a prime tool for answering these needs, with large and mega constellations planned, leading to a potential launch gap. An effective and commercially appealing solution is the development of small launchers, as these can complement the current available launch opportunity offer, serving a large pool of different types of clients, with a flexible and custom service that large conventional launchers cannot adequately assure. Rocket Factory Augsburg has partnered with CEiiA for the development of several structures for the RFA One rocket. The objective has been the design of solutions that are low-cost, light, and custom-made, applying design and manufacturing concepts as well as technologies from other industries, like the aeronautical and automotive, to the aerospace one. This allows for the implementation of a New Space approach to the launcher segment, while also building a supply chain and a set of solutions that enables the industrialisation of such structures for this and future small launchers. The two main systems under development have been a versatile Kick-Stage, for payload carrying and orbit insertion, and a sturdy Payload Fairing. Even though the use of components off-the-shelf have been widely accepted in the space industry for satellites, these two systems pose different challenges as they must be: highly reliable during the most extreme conditions imposed by the launch, so that they can be considered safe to launch all types of payloads. This paper thus dives deep on the solutions developed in the last few years, presenting also lessons learned during the manufacturing and testing of these structures.

Detection of the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and the Epoch of Reionization is one of the final frontiers of modern observational cosmology. The inherently faint signal makes it susceptible to contamination by several sources like astrophysical foregrounds and instrumental systematics. Nevertheless, developments achieved in the recent times will combine to make signal detection possible with the upcoming Square Kilometer Array (SKA), both statistically and via tomography. This review describes an indigenously developed end-to-end pipeline that simulates sensitive interferometric observations. It mainly focuses on the requirements for \hi detection in interferometers. In its present form, it can mimic the effects of realistic point source foregrounds and systematics- calibration error and position error on 21-cm observations. The performance of the pipeline is demonstrated for test cases with 0.01\% calibration error and position error. Its performance is consistent across telescope, foreground, and signal models. The focus of the simulation pipeline during the initial stages was for EoR science. But since this is a general interferometric simulation pipeline, it will be helpful to the entire SKA user community, irrespective of the science goals.

Baryonic Acoustic Oscillations (BAOs) studies based on the clustering of voids and matter tracers provide important constraints on cosmological parameters related to the expansion of the Universe. However, modelling the void exclusion effect is an important challenge for fully exploiting the potential of these kind of analyses. We thus develop two numerical methods to describe the clustering of cosmic voids. Neither model requires additional cosmological information beyond that assumed within the galaxy de-wiggled model. The models consist in power spectra whose performance we assess in comparison to a parabolic model on both Patchy boxes and light-cones. Moreover, we test their robustness against systematic effects and the reconstruction technique. The void model power spectra and the parabolic model with a fixed parameter provide strongly correlated values for the Alcock-Paczynski ($\alpha$) parameter, for boxes and light-cones likewise. The resulting $\alpha$ values - for all three models - are unbiased and their uncertainties are correctly estimated. However, the numerical models show less variation with the fitting range compared to the parabolic one. The Bayesian evidence suggests that the numerical techniques are often favoured compared to the parabolic model. Moreover, the void model power spectra computed on boxes can describe the void clustering from light-cones as well as from boxes. The same void model power spectra can be used for the study of pre- and post-reconstructed data-sets. Lastly, the two numerical techniques are resilient against the studied systematic effects. Consequently, using either of the two new void models, one can more robustly measure cosmological parameters.

Nearby galaxies provide populations of stellar and non-stellar sources at a common distance and in quantifiable environments. All are observed through the Milky Way foreground, with varying degrees of contamination that depend on observed Galactic latitude and the distance and size of the target galaxy. This work uses Gaia Data Release 3 (DR3) to identify foreground sources via astrometric measurements and thus quantify foreground contamination for a large sample of nearby galaxies. There are approximately half a million Gaia sources in the directions of 1401 galaxies listed in the Local Volume Galaxy catalogue (D<11 Mpc), excluding the largest Local Group galaxies. About two thirds of the Gaia sources have astrometric properties consistent with foreground sources; these sources are brighter, redder, and less centrally-concentrated than non-foreground sources. Averaged over galaxies, foreground sources make up 50 per cent of Gaia sources at projected radius r50=1.06*a26, where a26 is the angular diameter at the B=26.5 isophote. Foreground sources make up 50 per cent of Gaia sources at apparent magnitude m(G,50)=20.50. This limit corresponds to the tip of the red giant branch absolute magnitude at D = 450 kpc, and to the globular cluster luminosity function peak absolute magnitude at 5 Mpc. Gaia data provide a powerful tool for removing foreground contamination in stellar population studies of nearby galaxies, although Gaia foreground removal will be incomplete beyond distances of 5 Mpc.

We report on results from GBT/ARGUS $^{12}$CO(1-0) observations for the giant low surface brightness galaxy Malin 1, which allow us to determine an upper limit for its CO mass, and hence its molecular gas mass and molecular gas mass surface density $\Sigma_{H_2}$. Although we performed very deep observations through 17 hours on source integration time, reaching a noise level of $\sim 0.2$ mK (T$^{*}_{A}$) with a corresponding extended source CO limit (3$\sigma$) of 0.09 K km s$^{-1}$, 19 times more sensitive than previous works, we do not detect the $^{12}$CO(1-0) emission line. However, the observations allow us to estimate an upper limit (3$\sigma$) for the CO mass of about $7.4 \times 10^9$ M$_\odot$ for the extended emission, and $1.4 \times 10^8$ M$_\odot$ for the central part of the galaxy. With these figures we conclude that the molecular gas surface density is lower than 0.3 M$_\odot$ pc$^{-2}$, and the corresponding molecular to atomic gas mass ratio is lower than 0.13. The evidence suggests a quite different physical conditions for the interstellar medium in Malin 1 compared to that of normal, high surface brightness spirals. This, in one way to another, keeps an usual molecular gas tracer as CO hidden from our observations, in spite of the diverse stellar and structural properties of Malin 1 observed by several authors since more than 30 years.

We argue that the core degenerate (CD) scenario of type Ia supernovae (SNe Ia) can explain the compact helium-rich circumstellar material (CSM) of SN 2020eyj. In the new channel of the CD scenario that we propose there are two major common envelope evolution (CEE) phases. After the white dwarf (WD) companion removes the hydrogen-rich envelope of the asymptotic giant branch star its spiralling-in halts at few solar radii from the core, rather than continuing to the carbon-oxygen (CO) core as in the hydrogen-rich SNe Ia-CSM CD scenario. Only hundreds to tens of thousands of years later, after the helium-rich core expands, does the WD enters a CEE with the helium-rich layer. By that time the hydrogen-rich envelope is at a large distance from the center. The WD merges with the CO core during the second CEE phase, and only after a merger to explosion delay (MED) time of weeks to tens of years the merger remnant explodes. The SN Ia ejecta collides with a helium-rich CSM at tens to hundreds of AU. We follow the evolution of two stellar models with initial masses of 5Mo and 7Mo to their asymptotic giant branch phase when they are supposed to engulf the WD companion. We find that there is a sufficiently massive CO core to merge with the WD in the frame of the CD scenario as well as a massive helium-rich layer, 0.3-1Mo, to account for the helium-rich CSM of SN 2020eyj.

The electromagnetic emission resulting from neutron star mergers have been shown to encode properties of the ejected material in their light curves. The ejecta properties inferred from the kilonova emission has been in tension with those calculated based on the gravitational wave signal and numerical relativity models. Motivated by this tension, we construct a broad set of surrogate light curve models derived for kilonova ejecta. The four-parameter family of two-dimensional anisotropic simulations and its associated surrogate explore different assumptions about the wind outflow morphology and outflow composition, keeping the dynamical ejecta component consistent. We present the capabilities of these surrogate models in interpolating kilonova light curves across various ejecta parameters and perform parameter estimation for AT2017gfo both without any assumptions on the outflow and under the assumption that the outflow must be representative of solar \emph{r}-process abundance patterns. Our parameter estimation for AT2017gfo shows these new surrogate models help alleviate the ejecta property discrepancy while also illustrating the impact of systematic modeling uncertainties on these properties, urging further investigation.

Analysis of galaxy--galaxy strong lensing systems is strongly dependent on any prior assumptions made about the appearance of the source. Here we present a method of imposing a data-driven prior / regularisation for source galaxies based on denoising diffusion probabilistic models (DDPMs). We use a pre-trained model for galaxy images, AstroDDPM, and a chain of conditional reconstruction steps called denoising diffusion reconstruction model (DDRM) to obtain samples consistent both with the noisy observation and with the distribution of training data for AstroDDPM. We show that these samples have the qualitative properties associated with the posterior for the source model: in a low-to-medium noise scenario they closely resemble the observation, while reconstructions from uncertain data show greater variability, consistent with the distribution encoded in the generative model used as prior.

Betelgeuse is a pulsating red supergiant whose brightness is semi periodically variable and in February 2020 reached a historical minimum, the Great Dimming. The aims of this study are to characterize Betelgeuse's variability based on available archival data and to study possible causes of light variability. Many spectra, from ultraviolet and optical regions, were evaluated for spectral analysis. The spectra were used primarily to determine radial velocities from different layers of atmosphere and their long{-}term evolution. Additionally, photometric data were analyzed in different filters as well, to construct light curves and to determine periods of the variability. Spectroscopic and photometric variability are compared to each other and given into a context with the Great Dimming. The two most dominant photometric periods are $ P_{1} = 2190 \pm 270 \: \rm d $ and $ P_{2} = 417 \pm 17 \: \rm d $, while the dominant optical radial velocity periods are $ P_{1, v_{\rm r}} = 2510 \pm 440 \: \rm d $ and $ P_{2, v_{\rm r}} = 415 \pm 11 \: \rm d $. In the same time, the radial velocity determined from ultraviolet spectra also shows variability and is distinctively different from the variability of photospheric velocity, undergoing longer periods of variability. We attribute these velocities to the velocities at the base of outflowing wind. We also report a maximum of stellar wind velocity during the Great Dimming, accompanied by the previously reported minimum of brightness and the maximum of photospheric radial velocity. After the Dimming, Betelgeuse mode of variability has fundamentally changed and is now instead following a shorter period of $ \sim 200 \: \rm d $.

We present High-Level Science Products (HLSPs) containing light curves from MIT's Quick-Look Pipeline (QLP) from the second year of TESS' first Extended Mission (Sectors 40 - 55; 2021 July - 2022 September). In total, 12.2 million per-sector light curves for 5.7 million unique stars were extracted from 10-minute cadence Full-Frame Images (FFIs) and are made available to the community. As in previous deliveries, QLP HLSPs include both raw and detrended flux time series for all observed stars brighter than TESS magnitude T = 13.5 mag. Starting in Sector 41, QLP also produces light curves for select fainter M dwarfs. QLP has provided the community with one of the largest sources of FFI-extracted light curves to date since the start of the TESS mission.

Light curves feature many kinds of variability, including instrumental systematics, intrinsic stellar variability such as pulsations, and flux changes caused by transiting exoplanets or eclipsing binary stars. Detrending is a key pre-planet-search data processing step that aims to remove variability not due to transits. This data release note describes improvements to the Quick-Look Pipeline's detrending algorithm via the inclusion of quaternion data to remove short-timescale systematics. We describe updates to our procedure, intermediate data products outputted by the algorithm, and improvements to light curve precision.

An essential component of planetary climatology is knowledge of the tropospheric temperature field and its variability. Previous studies of Jupiter hinted at periodic behavior that was non-seasonal, as well as dynamical relationships between tropospheric and stratospheric temperatures. However, these observations were made over time frames shorter than Jupiter's orbit or they used sparse sampling. We derived upper-tropospheric (300-mbar) temperatures over 40 years, extending those studies to cover several orbits of Jupiter, revealing unexpected results. Periodicities of 4, 7 8-9 and 10-14 years were discovered that involved different latitude bands and seem disconnected from seasonal changes in solar heating. Anti-correlations of variability in opposite hemispheres were particularly striking at 16, 22 and 30 degrees from the equator. Equatorial temperature variations are also anticorrelated with those 60-70 km above. Such behavior suggests a top-down control of equatorial tropospheric temperatures from stratospheric dynamics. Realistic future global climate models must address the origins of these variations in preparation for their extension to a wider array of gas-giant exoplanets.

Alerts of potentially hazardous coronal mass ejections (CME) are based on the detection of rapid changes in remote observations of the solar atmosphere. This paper presents a method that detects and estimates the central coordinates of CME eruptions in Extreme Ultraviolet (EUV) data, with the dual aim of providing an early alert, and giving an initial estimate of the CME direction of propagation to a CME geometrical model. In particular, we plan to link the ALMANAC method to the CME detection and characterisation module of the Space Weather Empirical Ensemble Package (SWEEP), which is a fully automated modular software package for operational space weather capability currently being developed for the UK Meteorological Office. In this work, ALMANAC is applied to observations by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). As well as presenting the method, a proof of concept test is made on a limited set of data associated with twenty halo CMEs recorded by the Coordinated Data Analysis Workshop (CDAW) catalogue near the activity maximum of solar cycle 24. SDO/AIA data for each event is processed at 6 minute cadence to identify the on-disk location and time of each CME. The absolute mean deviance between the ALMANAC and CDAW source event coordinates are within 37.05 +- 29.71 minutes and 11.01 +- 10.39 degrees. These promising results give a solid foundation for future work, and will provide initial constraints to an automated CME alert and forecasting system.

During primordial star formation, the main cooling channel is provided by H$_{2}$ and super-molecules, such as H$_{2}$ or H$_{2}$, at sufficiently high densities. When the latter form at $n_{\rm H}$ $\geq$ $10^{14}$~cm$^{-3}$, collision-induced emission (CIE) provides efficient gas cooling. We investigate how CIE cooling affects the formation of metal-free binaries comparing simulations with and without this process. Irrespective of the cooling mechanism, we find a typical protostellar mass range between 0.01 to 100 M$_{\odot}$. However, models with only H$_{2}$ line cooling produce a greater number of low-mass protostars which exhibit stronger variations in their radial velocities than the high-mass protostars. Similarly, in models with both H$_{2}$ cooling and CIE cooling, significant variations in the radial velocities are found for protostars in the intermediate mass range. The initial number of fragments $N_{\rm max}$ decreases with increasing strength of turbulence. Cooling via super-molecules lets the most massive protobinaries (MMPBs) efficiently accrete mass. The maximum mass accretion rate $\dot M_{\rm max}$ for the MMPBs is more than an order of magnitude higher in the presence of CIE cooling than for pure H$_{2}$ line cooling. As a result, compact binaries with a semi-major axis as small as 3.57 au may form through the H$_{2}$ $-$ H$_{2}$ cooling channel. Our results indicate that in addition to the MMPBs most population III (Pop. III) binaries should be in eccentric i.e. non-circular orbits. This provides an important connection to the eccentric binaries reported in previous studies, which were found to exhibit rich temporal accretion signals during their evolution.

As part of the Portuguese Space Surveillance and Tracking (SST) program, a tracking radar and a double Wide Field of View Telescope system (4.3{\deg} x 2.3{\deg}) are being installed at the Pampilhosa da Serra Space Observatory (PASO) in the centre of continental Portugal, complementing an already installed deployable optical sensor for MEO and GEO surveillance. The tracking radar will track space debris in Low Earth Orbit (LEO) up to 1000 km and at the same time the telescope will also have LEO tracking capabilities. This article intends to discuss possible ways to take advantage of having these two sensors at the same location. Using both types of sensors takes advantage of the radar measurements which give precise radial velocity and distance to the objects, while the telescope gives better sky coordinates measurements. With the installation of radar and optical sensors, PASO can extend observation time of space debris and correlate information from optical and radar provenances in real time. During twilight periods both sensors can be used simultaneously to rapidly compute new TLEs for LEO objects, eliminating the time delays involved in data exchange between sites in a large SST network. This concept will not replace the need for a SST network with sensors in multiple locations around the globe, but will provide a more complete set of measurements from a given object passage, and therefore increase the added value for initial orbit determination, or monitoring of reentry campaigns of a given location. PASO will contribute to the development of new solutions to better characterize the objects improving the overall SST capabilities and constitute a perfect site for the development and testing of new radar and optical data fusion algorithms and techniques for space debris monitoring.

Recent multimessenger studies have provided evidence for high-energy neutrino sources that are opaque to GeV-TeV gamma-rays. We investigate the connection between high-energy neutrinos and gamma-rays in the active galaxy NGC 1068, and find that the neutrinos most likely come from regions within about 100 Schwarzschild radii. This is especially the case if neutrinos are produced via the photomeson production process, although the constraints could be alleviated if hadronuclear interactions are dominant. We consider the most favorable neutrino production regions, and discuss coronae, jets, winds, and their interactions with dense material. The results strengthen the importance of understanding dissipation mechanisms near the coronal region and the base of outflows. There could be a connection between active galactic nuclei with near-Eddington accretion and tidal disruptions events, in that neutrinos are produced in the obscured vicinity of supermassive black holes.

Type Ia supernovae (SNe Ia) are thought to be the result of thermonuclear explosions in white dwarfs (WDs). Commonly considered formation pathways include two merging WDs (the double degenerate channel), and a single WD accreting material from a H or He donor (the single degenerate channel). Since the predicted SNe Ia rates from WD in binaries are thought to be insufficient to explain the observed SNe Ia rate, it is important to study similar interactions in higher-order multiple star systems such as triple systems. We use the evolutionary population synthesis code Multiple Stellar Evolution (MSE) to study stellar evolution, binary interactions and gravitational dynamics of the triple-star systems. Also, unlike previous studies, prescriptions are included to simultaneously take into account the single and double degenerate channels, and we consider triples across the entire parameter space (including those with tight inner binaries). We explore the impact of typically ignored or uncertain physics such as fly-bys and CE prescription parameters on our results. The majority of systems undergo circular mergers to explode as SNe Ia, while eccentric collisions contribute to $0.4-4$ per cent of systems. The time-integrated SNe Ia rate from the triple channel is found to be $(3.60 \pm 0.04) \times {10^{-4}}\,\mathrm{M_{\odot}}^{-1}$ which is, surprisingly, similar to that of the isolated binary channel where the SNe Ia rate is $(3.2 \pm 0.1) \times {10^{-4}}\,\mathrm{M_{\odot}}^{-1}$. This implies that triples, when considering their entire parameter space, yield an important contribution to the overall SNe Ia rate.

We study the effects of light QCD axions on the stellar configuration of white dwarfs. At finite baryon density, the non-derivative coupling of the axion to nucleons displaces the axion from its in-vacuum minimum which implies a reduction of the nucleon mass. This dramatically alters the composition of stellar remnants. In particular, the modifications of the mass-radius relationship of white dwarfs allow us to probe large regions of unexplored axion parameter space without requiring it to be a significant fraction of dark matter.

Current sheets are spatially localized almost-1D structures with intense plasma currents. They play a key role in storing the magnetic field energy and they separate different plasma populations in planetary magnetospheres, the solar wind, and the solar corona. Current sheets are primary regions for the magnetic field line reconnection responsible for plasma heating and charged particle acceleration. One of the most interesting and widely observed type of 1D current sheets is the rotational discontinuity, that can be force-free or include plasma compression. Theoretical models of such 1D current sheets are based on the assumption of adiabatic motion of ions, i.e. ion adiabatic invariants are conserved. We focus on three current sheet configurations, widely observed in the Earth magnetopause and magnetotail and in the near-Earth solar wind. Magnetic field in such current sheets is supported by currents carried by transient ions, which exist only when there is a sufficient number of invariants. In this paper, we apply a novel machine learning approach, AI Poincar'e, to determine parametrical domains where adiabatic invariants are conserved. For all three current sheet configurations, these domains are quite narrow and do not cover the entire parametrical range of observed current sheets. We discuss possible interpretation of obtained results indicating that 1D current sheets are dynamical rather than static plasma equilibria.

We strengthen the cosmological bound on the axion mass, by solving the momentum-dependent Boltzmann equations for axion-pion scatterings and by using a phenomenological production rate derived from pion-pion scattering data, overcoming the breakdown of chiral perturbation theory. Using present cosmological datasets we obtain $m_a\leq 0.24~\text{eV}$. To further improve the bound and exploit the reach of upcoming cosmological surveys, reliable non-perturbative calculations above the QCD crossover are needed.

Phenomenological studies of cosmological collider physics in recent years have identified many 1-loop inflation correlators as leading channels for discovering heavy new particles around or above the inflation scale. However, complete analytical results for these massive 1-loop correlators are currently unavailable. In this work, we embark on a program of bootstrapping inflation correlators with massive exchanges at 1-loop order, with the input of tree-level inflation correlators and the techniques of spectral decomposition in dS. As a first step, we present for the first time the complete and analytical results for a class of 4-point and 3-point inflation correlators mediated by massive scalar fields at the 1-loop order. Using the full result, we provide simple and reliable analytical approximations for the signals and the background in the squeezed limit. We also identify configurations of the scalar trispectrum where the oscillatory signal from the loop is dominant over the background.

A major threat to satellites is space debris with their low mass and high rotational speed. Accordingly, the short observation time of these objects is a major limitation in space research for appropriate detection and decision. As a result, these objects do not fully illuminated, leading to their incomplete images at any snapshot. In this paper, we propose a method to decrease the number of snapshots in a given observation time and using a limited number of spot beams per snapshot called the encoded aperture. To recover the space debris images, an inverse problem is defined based on compressive sensing methods. Also, we show that for satellite imaging the T V norm is more appropriate. We develop a procedure to recover space debris and satellites using L1 and T V norms. Using simulation results, we compare the results with the well-known SBL and SL0 norm in terms of the number of snapshots, MSE, SNR, and running time. It is shown that our proposed method can successfully recover the space objects images using a fewer number of snapshots.

Gravitational signatures of black hole superradiance are a unique probe of ultralight particles that are weakly-coupled to ordinary matter. The existence of an ultralight boson would lead spinning black holes with size comparable to the Compton wavelength of the boson to become superradiantly unstable to forming an oscillating cloud, spinning down the black hole, and radiating gravitational waves in the process. However, maximizing the chance of observing such signals or, in their absence, placing the strongest constraints on the existence of such particles, requires accurate theoretical predictions. In this work, we introduce a new gravitational waveform model, SuperRad, that models the dynamics, oscillation frequency, and gravitational wave signals of these clouds by combining numerical results in the relativistic regime with fits calibrated to analytical estimates, covering the entire parameter space of ultralight scalar and vector clouds with the lowest two azimuthal numbers ($m = 1$ and $2$). We present new calculations of the gravitational wave frequency evolution as the boson cloud dissipates, including using fully general-relativistic methods to quantify the error in more approximate treatments. Finally, as a first application, we assess the viability of conducting follow-up gravitational wave searches for ultralight vector clouds around massive black hole binary merger remnants. We show that LISA may be able to probe vector masses in the range from $1\times 10^{-16}$ eV to $6\times 10^{-16}$ eV using follow-up gravitational wave searches.

Space weather observations and modeling at Mars have begun but they must be significantly increased to support the future of Human Exploration on the Red Planet. A comprehensive space weather understanding of a planet without a global magnetosphere and a thin atmosphere is very different from our situation at Earth so there is substantial fundamental research remaining. It is expected that the development of suitable models will lead to a comprehensive operational Mars space weather alert (MSWA) system that would provide rapid dissemination of information to Earth controllers, astronauts in transit, and those in the exploration zone (EZ) on the surface by producing alerts that are delivered rapidly and are actionable. To illustrate the importance of such a system, we use a magnetohydrodynamic code to model an extreme Carrington-type coronal mass ejection (CME) event at Mars. The results show a significant induced surface field of nearly 3000 nT on the dayside that could radically affect unprotected electrical systems that would dramatically impact human survival on Mars. Other associated problems include coronal mass ejection (CME) shock-driven acceleration of solar energetic particles producing large doses of ionizing radiation at the Martian surface. In summary, along with working more closely with international partners, the next Heliophysics Decadal Survey must include a new initiative to meet expected demands for space weather forecasting in support of humans living and working on the surface of Mars. It will require significant effort to coordinate NASA and the international community contributions.

We compute the leading and sub-leading spin effects through the second post-Newtonian order (2PN) in spherical harmonic modes of gravitational waveforms from inspiralling compact binaries in non-circular orbits with non-precessing components. The two spin couplings, linear-in-spin (spin-orbit; SO) and quadratic-in-spin (spin-spin; SS), that appear in 2PN waveforms are computed with desired accuracy and explicit expressions for relevant modes are derived. The modes that have spin corrections through 2PN include $(\ell, |m|)$=$((2,2),\,(2,1),\,(3,3),\,(3,2),\,(3,1),\,(4,3),\,(4,1))$ modes. Closed form expressions for these modes for compact binaries in general orbits as well as in elliptical orbits are being provided. While the general orbit results can be used to study signals from binaries in orbits of arbitrary shape and nature, elliptical orbit results are applicable to systems with arbitrary eccentricities. We also express the elliptical orbit results as leading eccentric corrections to the circular results. Our prescription represents, the first, fully analytical treatment that combines spins, eccentricity and higher modes together and completes computation of spin effects through 2PN order. These should find immediate applications in inspiral-merger-ringdown modelling for eccentric mergers including the effect of non-precessing spins and higher modes as well as in parameter estimation analyses employing inspiral waveform.

We investigate the impact of stochastic quantum noise due to trans--Planckian effects on the primordial power spectrum for gravity waves during inflation. Given an energy scale Lambda, expected to be close to the Planck scale m_Pl and larger than the Hubble scale H, this noise is described in terms of a source term in the evolution equation for comoving modes k which changes its amplitude growth from early times as long as the mode physical wavelength is smaller than Lambda^-1. We model the source term as due to a gas of black holes in the trans--Planckian regime and the corresponding Hawking radiation. In fact, for energy scales larger than, or of the order of Lambda, it is expected that trapped surfaces may form due to large energy densities. At later times the evolution then follows the standard sourceless evolution. We find that this mechanism still leads to a scale-invariant power spectrum of tensor perturbations, with an amplitude that depends upon the ratio Lambda/m_Pl.

Astrochemical surface reactions are thought to be responsible for the formation of complex organic molecules, which are of potential importance for the origin of life. In a situation, when the chemical composition of dust surfaces is not precisely known, the fundamental knowledge concerning such reactions gains significance. We describe an experimental technique, which can be used to measure the energy released in reactions of a single pair of reactants. These data can be directly compared with the results of quantum chemical computations leading to unequivocal conclusions regarding the reaction pathways and the presence of energy barriers. It allows for predicting the outcomes of astrochemical surface reactions with higher accuracy compared to that achieved based on gas-phase studies. However, for the highest accuracy, some understanding of the catalytic influence of specific surfaces on the reactions is required. The new method was applied to study the reactions of C atoms with H2, O2, and C2H2. The formation of HCH, CO + O, and triplet cyclic-C3H2 products has been revealed, correspondingly.

Key to the simplicity of supergravity alpha-attractor models of inflation are Volkov-Akulov fermions, often in the form of nilpotent superfields. Here we explore the possibility of using the double-copy to construct theories of Dirac-Born-Infeld-Volkov-Akulov (DBIVA) coupled to supergravity. A color-dual bootstrap admits scattering amplitudes involving pions and vectors through five-point tree-level order by order in mass-dimension, but requires the introduction of a tr(F^3) operator. Gauge theories with this operator were recently found to require a tower of higher-derivative operators to be compatible with the duality between color and kinematics. Adjoint-type double-copy construction at its most conservative seems to require the UV completion of DBVIA + pure Poincare supergravity scattering amplitudes to a family of theories involving DBVIA-like particles coupled to Weyl-Einstein supergravity. We also point out an alternative solution to color-dual gauged pions that allows adjoint double-copy without a tower of higher derivative corrections but at the cost of exchange symmetry between scalars.