Papers for Friday, Nov 03 2023

We provide a detailed study of natural inflation with a periodic non-minimal coupling, which is a well-motivated inflationary model that admits an explicit UV completion. We demonstrate that this construction can satisfy the most recent observational constraints from Planck and the BICEP/Keck collaborations. We also compute the corresponding relic gravitational wave background due to tensor perturbations and show that future space-borne interferometers, such as DECIGO, BBO and ALIA, may be able to detect it. Next, we extend this analysis and establish the validity of these results in a multi-field model featuring an additional $R^2$ term in the action, which allows us to interpolate between natural and scalaron (a.k.a.~Starobinsky) inflation. We investigate the conditions under which the aforementioned future interferometers will have the capability to differentiate between pure natural inflation and natural-scalaron inflation. The latter analysis could open the door to distinguishing between single-field and multi-field inflation through gravitational wave observations in more general contexts.

The apparent tension between the luminosity functions of red supergiant (RSG) stars and of RSG progenitors of Type II supernovae (SNe) is often referred to as the RSG problem and it motivated some to suggest that many RSGs end their life without a SN explosion. However, the luminosity functions of RSG SN progenitors presented so far were biased to high luminosities, because the sensitivity of the search was not considered. Here, we use limiting magnitudes to calculate a bias-corrected RSG progenitor luminosity function. We find that only $(36\pm11)\%$ of all RSG progenitors are brighter than a bolometric magnitude of $-7\,\text{mag}$, a significantly smaller fraction than $(56\pm5)\%$ quoted by Davies & Beasor (2020). The larger uncertainty is due to the relatively small progenitor sample, while uncertainties on measured quantities such as magnitudes, bolometric corrections, extinction, or SN distances, only have a minor impact, as long as they fluctuate randomly for different objects in the sample. When comparing the luminosity function of RSG SN progenitors to Type M supergiants in the Large Magellanic cloud, we find that they are consistent, due to the flatter shape of the progenitor luminosity function. The RSG progenitor luminosity function, hence, does not imply the existence of failed SNe. The presented statistical method is not limited to progenitor searches, but applies to any situation in which a measurement is done for a sample of detected objects, but the probed quantity or property can only be determined for part of the sample.

Gravitational microlensing is one of the strongest observational techniques to observe non-luminous astrophysical bodies. Existing microlensing observations provide tantalizing evidence of a population of low-mass objects whose origin is unknown. These events may be caused by terrestrial-mass free-floating planets or by exotic objects such as primordial black holes. However, the nature of these objects cannot be resolved on an event-by-event basis, as the induced light curve is degenerate for lensing bodies of identical mass. One must instead statistically compare \textit{distributions} of lensing events to determine the nature of the lensing population. While existing surveys lack the statistics required to identify multiple subpopulations of lenses, this will change with the launch of the Nancy Grace Roman Space Telescope. Roman's Galactic Bulge Time Domain Survey is expected to observe hundreds of low-mass microlensing events, enabling a robust statistical characterization of this population. In this paper, we show that by exploiting features in the distribution of lensing event durations, Roman will be sensitive to a subpopulation of primordial black holes hidden amongst a background of free-floating planets. Roman's reach will extend to primordial black hole dark matter fractions as low as $f_\text{PBH} = 10^{-4}$ at peak sensitivity, and will be able to conclusively determine the origin of existing ultrashort-timescale microlensing events. A positive detection would provide evidence that a significant fraction of the cosmological dark matter consists of macroscopic, non-luminous objects.

Galaxies of different types are not equally distributed in the Local Universe. In particular, the supergalactic plane is prominent among the brightest ellipticals, but inconspicuous among the brightest disk galaxies. This striking difference provides a unique test for our understanding of galaxy and structure formation. Here we use the SIBELIUS DARK constrained simulation to confront the predictions of the standard Lambda Cold Dark Matter ($\Lambda$CDM) model and standard galaxy formation theory with these observations. We find that SIBELIUS DARK reproduces the spatial distributions of disks and ellipticals and, in particular, the observed excess of massive ellipticals near the supergalactic equator. We show that this follows directly from the local large-scale structure and from the standard galaxy formation paradigm, wherein disk galaxies evolve mostly in isolation, while giant ellipticals congregate in the massive clusters that define the supergalactic plane. Rather than being anomalous as earlier works have suggested, the distributions of giant ellipticals and disks in the Local Universe and in relation to the supergalactic plane are key predictions of the $\Lambda$CDM model.

We present 0.075 (~400 pc) resolution ALMA observations of the [CII] and dust continuum emission from the host galaxy of the z=6.5406 quasar, P036+03. We find that the emission arises from a thin, rotating disk with an effective radius of 0.21" (1.1 kpc). The velocity dispersion of the disk is consistent with a constant value of 66.4+-1.0 km/s, yielding a scale height of 80+-30 pc. The [CII] velocity field reveals a distortion that we attribute to a warp in the disk. Modeling this warped disk yields an inclination estimate of 40.4+-1.3 degrees and a rotational velocity of 116+-3 km/s. The resulting dynamical mass estimate of (1.96+-0.10) x 10^10 Msun is lower than previous estimates, which strengthens the conclusion that the host galaxy is less massive than expected based on local scaling relations between the black hole mass and the host galaxy mass. Using archival MUSE Ly-alpha observations, we argue that counter-rotating halo gas could provide the torque needed to warp the disk. We further detect a region with excess (15-sigma) dust continuum emission, which is located 1.3 kpc northwest of the galaxy's center and is gravitationally unstable (Toomre-Q < 0.04). We posit this is a star-forming region whose formation was triggered by the warp, because the region is located within a part of the warped disk where gas can efficiently lose angular momentum. The combined ALMA and MUSE imaging provides a unique view of how gas interactions within the disk-halo interface can influence the growth of massive galaxies within the first billion years of the universe.

Future space-based far-infrared astrophysical observatories will require exquis-itely sensitive detectors consistent with the low optical backgrounds. The PRobe far-Infrared Mission for Astrophysics (PRIMA) will deploy arrays of thousands of superconducting kinetic inductance detectors (KIDs) sensitive to radiation between 25 and 265 $\mu$m. Here, we present laboratory characterization of prototype, 25 -- 80 $\mu$m wavelength, low-volume, aluminum KIDs designed for the low-background environment expected with PRIMA. A compact parallel plate capacitor is used to minimize the detector footprint and suppress TLS noise. A novel resonant absorber is designed to enhance response in the band of interest. We present noise and optical efficiency measurements of these detectors taken with a low-background cryostat and a cryogenic blackbody. A microlens-hybridized KID array is found to be photon noise limited down to about 50 aW with a limiting detector NEP of about $6.5 \times 10^{-19}~\textrm{W/Hz}^{1/2}$. A fit to an NEP model shows that our optical system is well characterized and understood down to 50 aW. We discuss future plans for low-volume aluminum KID array development as well as the testbeds used for these measurements.

To understand high precision observations of exoplanets and brown dwarfs, we need detailed and complex general circulation models (GCMs) that incorporate hydrodynamics, chemistry, and radiation. In this study, we specifically examine the coupling between chemistry and radiation in GCMs and compare different methods for mixing opacities of different chemical species in the correlated-k assumption, when equilibrium chemistry cannot be assumed. We propose a fast machine learning method based on DeepSets (DS), which effectively combines individual correlated-k opacities (k-tables). We evaluate the DS method alongside other published methods like adaptive equivalent extinction (AEE) and random overlap with rebinning and resorting (RORR). We integrate these mixing methods into our GCM (expeRT/MITgcm) and assess their accuracy and performance for the example of the hot Jupiter HD~209458 b. Our findings indicate that the DS method is both accurate and efficient for GCM usage, whereas RORR is too slow. Additionally, we observe that the accuracy of AEE depends on its specific implementation and may introduce numerical issues in achieving radiative transfer solution convergence. We then apply the DS mixing method in a simplified chemical disequilibrium situation, where we model the rainout of TiO and VO, and confirm that the rainout of TiO and VO would hinder the formation of a stratosphere. To further expedite the development of consistent disequilibrium chemistry calculations in GCMs, we provide documentation and code for coupling the DS mixing method with correlated-k radiative transfer solvers. The DS method has been extensively tested to be accurate enough for GCMs, however, other methods might be needed for accelerating atmospheric retrievals.

The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics using single-crystal silicon, and developments in the fabrication of large-format, small-pixel, high readout rate CCD detectors with good spectral resolution, allowing a robust and cost-effective design. Further, AXIS will be responsive to target-of-opportunity alerts and, with onboard transient detection, will be a powerful facility for studying the time-varying X-ray universe, following on from the legacy of the Neil Gehrels (Swift) X-ray observatory that revolutionized studies of the transient X-ray Universe. In this paper, we present an overview of AXIS, highlighting the prime science objectives driving the AXIS concept and how the observatory design will achieve these objectives.

We investigate the possibility of learning the representations of cosmological multifield dataset from the CAMELS project. We train a very deep variational encoder on images which comprise three channels, namely gas density (Mgas), neutral hydrogen density (HI), and magnetic field amplitudes (B). The clustering of the images in feature space with respect to some cosmological/astrophysical parameters (e.g. $\Omega_{\rm m}$) suggests that the generative model has learned latent space representations of the high dimensional inputs. We assess the quality of the latent codes by conducting a linear test on the extracted features, and find that a single dense layer is capable of recovering some of the parameters to a promising level of accuracy, especially the matter density whose prediction corresponds to a coefficient of determination $R^{2}$ = 0.93. Furthermore, results show that the generative model is able to produce images that exhibit statistical properties which are consistent with those of the training data, down to scales of $k\sim 4h/{\rm Mpc}.$

We present results from the characterization and optimization of six Skipper CCDs for use in a prototype focal plane for the SOAR Integral Field Spectrograph (SIFS). We tested eight Skipper CCDs and selected six for SIFS based on performance results. The Skipper CCDs are 6k $\times$ 1k, 15 $\mu$m pixels, thick, fully-depleted, $p$-channel devices that have been thinned to $\sim 250 \mu$m, backside processed, and treated with an antireflective coating. We optimize readout time to achieve $<4.3$ e$^-$ rms/pixel in a single non-destructive readout and $0.5$ e$^-$ rms/pixel in $5 \%$ of the detector. We demonstrate single-photon counting with $N_{\rm samp}$ = 400 ($\sigma_{\rm 0e^-} \sim$ 0.18 e$^-$ rms/pixel) for all 24 amplifiers (four amplifiers per detector). We also perform conventional CCD characterization measurements such as cosmetic defects ($ <0.45 \%$ ``bad" pixels), dark current ($\sim 2 \times 10^{-4}$ e$^-$/pixel/sec.), charge transfer inefficiency ($3.44 \times 10^{-7}$ on average), and charge diffusion (PSF $< 7.5 \mu$m). We report on characterization and optimization measurements that are only enabled by photon-counting. Such results include voltage optimization to achieve full-well capacities $\sim 40,000-63,000$ e$^-$ while maintaining photon-counting capabilities, clock induced charge optimization, non-linearity measurements at low signals (few tens of electrons). Furthermore, we perform measurements of the brighter-fatter effect and absolute quantum efficiency ($\gtrsim\, 80 \%$ between 450 nm and 980 nm; $\gtrsim\,90 \%$ between 600 nm and 900 nm) using Skipper CCDs.

Future far-infrared astrophysics observatories will require focal plane arrays containing thousands of ultra-sensitive, superconducting detectors, each of which needs to be optically coupled to the telescope. At longer wavelengths, many approaches have been developed including feedhorn arrays and macroscopic arrays of lenslets. However, with wavelengths as short as 25 microns, optical coupling in the far-infrared remains challenging. In this paper, we present a novel approach for fabricating far-infrared monolithic silicon microlens arrays using grayscale lithography and deep reactive ion etching. The design, fabrication, and characterization of the microlens arrays are discussed. We compare the designed and fabricated lens profile, and calculate that the fabricated lenses will achieve 84% encircled power for the designed detector, which is only 3% less than the designed performance. We also present methods developed for anti-reflection coating microlens arrays and for a silicon-to-silicon die bonding process to hybridize microlens arrays with detector arrays.

We provide a measurement of the deficit in the Sunyaev-Zel'dovich Compton-$y$ signal towards cosmic voids, by stacking a catalogue of 97,090 voids constructed with BOSS-DR12 data, on the $y$ maps built on data from the Atacama Cosmology Telescope (ACT) DR4 and the Planck satellite. We detect the void signal with a significance of $7.3\,\sigma$ with ACT and $9.7\,\sigma$ with Planck, obtaining agreements in the associated void radial $y$ profiles extracted from both maps. The inner-void profile (for angular separations within the void angular radius) is reconstructed with significances of $4.7\sigma$ and $6.1\sigma$ with ACT and Planck, respectively; we model such profile using a simple model that assumes uniform gas (under)density and temperature, which enables us to place constraints on the product $(-\delta_{\rm v}T_{\rm e})$ of the void density contrast (negative) and the electron temperature. The best-fit values from the two data sets are $(-\delta_{\rm v}T_{\rm e})=(6.5\pm 2.3)\times 10^{5}\,\text{K}$ for ACT and $(8.6 \pm 2.1)\times 10^{5}\,\text{K}$ for Planck ($68\%$ C.L.), which are in good agreement under uncertainty. The data allow us to place lower limits on the expected void electron temperature at $2.7\times10^5\,\text{K}$ with ACT and $5.1\times10^5\,\text{K}$ with Planck ($95\%$ C.L.); these results can transform into upper limits for the ratio between the void electron density and the cosmic mean as $n^{\rm v}_{\rm e}/\bar{n}_{\rm e}\leqslant 0.73$ and $0.49$ ($95\%$ C.L.), respectively. Our findings prove the feasibility of using tSZ observations to constrain the gas properties inside cosmic voids, and confirm that voids are under-pressured regions compared to their surroundings.

Efficiently analyzing maps from upcoming large-scale surveys requires gaining direct access to a high-dimensional likelihood and generating large-scale fields with high fidelity, which both represent major challenges. Using CAMELS simulations, we employ the state-of-the-art score-based diffusion models to simultaneously achieve both tasks. We show that our model, HIDM, is able to efficiently generate high fidelity large scale HI maps that are in a good agreement with the CAMELS's power spectrum, probability distribution, and likelihood up to second moments. HIDM represents a step forward towards maximizing the scientific return of future large scale surveys.

Most of a galaxy's mass is located out to hundreds of kiloparsecs beyond its stellar component. This diffuse reservoir of gas, the circumgalactic medium (CGM), acts as the interface between a galaxy and the cosmic web that connects galaxies. We present kiloparsec-scale resolution integral field spectroscopy of emission lines that trace cool ionized gas from the center of a nearby galaxy to 30 kpc into its CGM. We find a smooth surface brightness profile with a break in slope at twice the 90% stellar radius. At this radius, the gas also transitions from being photoionized by HII star-forming regions in the disk to being ionized by shocks or the extragalactic UV background at larger distances. These changes represent the boundary between the interstellar medium (ISM) and the CGM, revealing the shape and extent of the dominant reservoir of baryonic matter in galaxies for the first time.

In 2020, the HAWC Collaboration presented the first catalog of gamma-ray sources emitting above 56 TeV and 100 TeV. With nine sources detected, this was the highest-energy source catalog to date. Here, we present the results of re-analysis of the old data, along with additional data acquired since then. We use a new version of the reconstruction software with better pointing accuracy and improved gamma/hadron separation. We now see more than 25 sources above 56 TeV, with most sources being located in the Galactic plane. The vast majority of these seem to be leptonic pulsar wind nebulae, but some have been shown to have hadronic emission. We will show spectra and discuss possible emission mechanisms of some of the most interesting sources, including the ones the HAWC Collaboration considers PeVatron candidates.

Line intensity mapping is a promising probe of the universe's large-scale structure. We explore the sensitivity of the DSA-2000, a forthcoming array consisting of over 2000 dishes, to the statistical power spectrum of neutral hydrogen's 21 cm emission line. These measurements would reveal the distribution of neutral hydrogen throughout the near-redshift universe without necessitating resolving individual sources. The success of these measurements relies on the instrument's sensitivity and resilience to systematics. We show that the DSA-2000 will have the sensitivity needed to detect the 21 cm power spectrum at z=0.5 and across power spectrum modes of 0.03-31.32 h/Mpc with 0.1 h/Mpc resolution. We find that supplementing the nominal array design with a dense core of 200 antennas will expand its sensitivity at low power spectrum modes and enable measurement of Baryon Acoustic Oscillations (BAOs). Finally, we present a qualitative discussion of the DSA-2000's unique resilience to sources of systematic error that can preclude 21 cm intensity mapping.

We investigate the feasibility of using COmoving Lagrangian Acceleration (COLA) technique to efficiently generate galaxy mock catalogues that can accurately reproduce the statistical properties of observed galaxies. Our proposed scheme combines the subhalo abundance matching (SHAM) procedure with COLA simulations, utilizing only three free parameters: the scatter magnitude ($\sigma_{\rm scat}$) in SHAM, the initial redshift ($z_{\rm init}$) of the COLA simulation, and the time stride ($da$) used by COLA. In this proof-of-concept study, we focus on a subset of BOSS CMASS NGC galaxies within the redshift range $z\in [0.45, 0.55]$. We perform $\mathtt{GADGET}$ simulation and low-resolution COLA simulations with various combinations of $(z_{\rm init}, da)$, each using $1024^{3}$ particles in an $800~h^{-1}{\rm Mpc}$ box. By minimizing the difference between COLA mock and CMASS NGC galaxies for the monopole of the two-point correlation function (2PCF), we obtain the optimal $\sigma_{\rm scat}$. We have found that by setting $z_{\rm init}=29$ and $da=1/30$, we achieve a good agreement between COLA mock and CMASS NGC galaxies within the range of 4 to $20~h^{-1}{\rm Mpc}$, with a computational cost two orders of magnitude lower than that of the N-body code. Moreover, a detailed verification is performed by comparing various statistical properties, such as anisotropic 2PCF, three-point clustering, and power spectrum multipoles, which shows similar performance between GADGET mock and COLA mock catalogues with the CMASS NGC galaxies. Furthermore, we assess the robustness of the COLA mock catalogues across different cosmological models, demonstrating consistent results in the resulting 2PCFs. Our findings suggest that COLA simulations are a promising tool for efficiently generating mock catalogues for emulators and machine learning analyses in exploring the large-scale structure of the Universe.

Using empty-field `Quick Look' images from the first two epochs of the VLA Sky Survey (VLASS) observations, centred on the positions of $\sim3700$ individually radio-non-detected active galactic nuclei (AGNs) at $z\ge4$, we performed image stacking analysis to examine the sub-mJy emission at $3$ GHz. We found characteristic monochromatic radio powers of $P_\mathrm{char}=(2-13) \times 10^{24}$ W Hz$^{-1}$, $P_\mathrm{char}=2\times10^{24}-1.3\times10^{25}$ W Hz$^{-1}$, indicating that AGN-related radio emission is widespread in the sample. The signal-to-noise ratios of the redshift-binned median stacked maps are between $4-6$, and we expect that with the inclusion of the yet to be completed third-epoch VLASS observations, the detection limit defined as signal-to-noise ratio $\mathrm{SNR}\ge6$ could be reached, and the redshift dependence can be determined. To obtain information on the general spectral properties of the faint radio emission in high-redshift AGNs, we confined the sample to $\sim3000$ objects covered by both the VLASS and the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey. We found that the flux densities from the median stacked maps show a characteristic spectral index of $\alpha^*=-0.30\pm0.15$, which is in agreement with the median spectral index of the radio-detected $z\ge4$ AGNs from our high-redshift AGN catalogue. The three-band mid-infrared colour--colour diagram based on Wide-field Infrared Survey Explorer observations provides further support regarding the AGN contribution to the radio emission in the sub-mJy sample.

We analyze TYPHOON long slit absorption line spectra of the starburst barred spiral galaxy NGC 1365 obtained with the Progressive Integral Step Method covering an area of 15 square kpc. Applying a population synthesis technique, we determine the spatial distribution of ages and metallicity of the young and old stellar population together with star formation rates, reddening, extinction and the ratio R$_V$ of extinction to reddening. We detect a clear indication of inside-out growth of the stellar disk beyond 3 kpc characterized by an outward increasing luminosity fraction of the young stellar population, a decreasing average age and a history of mass growth, which was finished 2 Gyrs later in the outermost disk. The metallicity of the young stellar population is clearly super solar but decreases towards larger galactocentric radii with a gradient of -0.02 dex/kpc. On the other hand, the metal content of the old population does not show a gradient and stays constant at a level roughly 0.4 dex lower than that of the young population. In the center of NGC 1365 we find a confined region where the metallicity of the young population drops dramatically and becomes lower than that of the old population. We attribute this to infall of metal poor gas and, additionally, to interrupted chemical evolution where star formation is stopped by AGN and supernova feedback and then after several Gyrs resumes with gas ejected by stellar winds from earlier generations of stars. We provide a simple model calculation as support for the latter.

Context. Blazar AO 0235+164, located at redshift z = 0.94, has undergone several sharp multi-spectral-range flaring episodes during the last decades. In particular, the episodes peaking in 2008 and 2015, that received extensive multi-wavelength coverage, exhibited interesting behavior. Aims. We study the actual origin of these two observed flares by constraining the properties of the observed photo-polarimetric variability, those of the broad-band spectral energy-distribution and the observed time-evolution behavior of the source as seen by ultra-high resolution total-flux and polarimetric Very-long-baseline interferometry (VLBI) imaging. Methods. The analysis of VLBI images allows us to constrain kinematic and geometrical parameters of the 7 mm jet. We use the Discrete Correlation Function to compute the statistical correlation and the delays between emission at different spectral ranges. Multi-epoch modeling of the spectral energy distributions allows us to propose specific models of emission; in particular for the unusual spectral features observed in this source in the X-ray region of the spectrum during strong multi spectral-range flares. Results. We find that these X-ray spectral features can be explained by an emission component originating in a separate particle distribution than the one responsible for the two standard blazar bumps. This is in agreement with the results of our correlation analysis that do not find a strong correlation between the X-rays and the remaining spectral ranges. We find that both external Compton dominated and synchrotron self-Compton dominated models can explain the observed spectral energy distributions. However, synchrotron self-Compton models are strongly favored by the delays and geometrical parameters inferred from the observations.

We have carried out HST snapshot observations at 1.1 $\mu$m of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel space observatory. Our candidates comprise systems with flux densities at $500\,\mu$m$ S_{500}\geq 80$ mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the foreground lens candidate. After combining visual inspection, archival high-resolution observations, and lens subtraction, we divide the systems into different classes according to their lensing likelihood. We confirm 65 systems to be lensed. Of these, 30 are new discoveries. We successfully perform lens modelling and source reconstruction on 23 systems, where the foreground lenses are isolated galaxies and the background sources are detected in the HST images. All the systems are successfully modelled as a singular isothermal ellipsoid. The Einstein radii of the lenses and the magnifications of the background sources are consistent with previous studies. However, the background source circularised radii (between 0.34 kpc and 1.30 kpc) are $\sim$3 times smaller than the ones measured in the sub-mm/mm for a similarly selected and partially overlapping sample. We compare our lenses with those in the SLACS survey, confirming that our lens-independent selection is more effective at picking up fainter and diffuse galaxies and group lenses. This sample represents the first step towards characterising the near-IR properties and stellar masses of the gravitationally lensed dusty star-forming galaxies.

[Abridged] JWST observations of the Orion Bar have shown the incredible richness of PAH bands and their variation on small scales. We aim to probe the photochemical evolution of PAHs across the key zones of the photodissociation region (PDR) that is the Orion Bar using unsupervised machine learning. We use NIRSpec and MIRI IFU data from the JWST ERS Program PDRs4All. We lever bisecting k-means clustering to generate detailed spatial maps of the spectral variability in several wavelength regions. We discuss the variations in the cluster profiles and connect them to the local physical conditions. We interpret these variations with respect to the key zones: the HII region, the atomic PDR zone, and the three dissociation fronts. The PAH emission exhibits spectral variation that depends strongly on spatial position in the PDR. We find the 8.6um band to behave differently than all other bands which vary systematically with one another. We find uniform variation in the 3.4-3.6um bands and 3.4/3.3 intensity ratio. We attribute the carrier of the 3.4-3.6um bands to a single side group attached to very similarly sized PAHs. Cluster profiles reveal a transition between characteristic profiles classes of the 11.2um feature from the atomic to the molecular PDR zone. We find the carriers of each of the profile classes to be independent, and reason the latter to be PAH clusters existing solely deep in the molecular PDR. Clustering also reveals a connection between the 11.2 and 6.2um bands; and that clusters generated from variation in the 10.9-11.63um region can be used to recover those in the 5.95-6.6um region. Clustering is a powerful tool for characterizing PAH variability on both spatial and spectral scales. For individual bands as well as global spectral behaviours, we find UV-processing to be the most important driver of the evolution of PAHs and their spectral signatures in the Orion Bar.

Urea is a prebiotic molecule that has been detected in few sources of the interstellar medium (ISM) and in Murchison meteorite. Being stable against ultraviolet radiation and high-energy electron bombardment, urea is expected to be present in interstellar ices. Theoretical and experimental studies suggest that isocyanic acid (HNCO) and formamide (NH$_2$CHO) are possible precursors of urea. However, uncertainties still exist regarding its formation routes. Previous computational works characterised urea formation in the gas phase or in presence of few water molecules by reaction of formamide with nitrogen-bearing species. In this work, we investigated the reaction of HNCO + NH$_3$ on an 18 water molecules ice cluster model mimicking interstellar ice mantles by means of quantum chemical computations. We characterised different mechanisms involving both closed-shell and open-shell species at B3LYP-D3(BJ)/ma-def2-TZVP level of theory, in which the radical-radical H$_2$NCO + NH$_2$ coupling has been found to be the most favourable one due to being almost barrierless. In this path, the presence of the icy surfaces is crucial for acting as reactant concentrators/suppliers, as well as third bodies able to dissipate the energy liberated during the urea formation.

We present the results from an analysis of data from an \textit{XMM-Newton} observation of the accreting high mass X-ray binary pulsar GX 301$-$2. Spectral analysis in the non-flaring segment of the observation revealed that the equivalent width of the iron fluorescence emission is correlated with the observed absorption column density and the ratio of the iron K$\beta$ and K$\alpha$ line strength varied with the flux of the source. Coherent pulsations were detected with the spin period of the pulsar of 687.9$\pm$0.1 s, and a secondary pulsation was also detected with a period of 671.8$\pm$0.2 s, most prominent in the energy band of the iron line. At the spin period of the neutron star, the pulsation of the iron line has a low amplitude and the profile is different from the continuum. Pulse phase-resolved spectroscopy also revealed pulsations of the iron emission line during the non-flaring segment of the light curve. At the secondary period, both the iron line and the continuum have nearly identical pulse fraction and pulse profile. The additional periodicity can be attributed to the beat frequency between the spin of the neutron star and the Keplerian frequency of a stellar wind clump in retrograde motion around the neutron star. Reprocessed X-ray emissions originating from the clump can produce the observed secondary pulsations both in the continuum and the iron fluorescence line. The clump rotating around the neutron star is estimated to be approximately five lt-s away from the neutron star.

Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the $\textit{Voyager}$ spacecrafts. Although the volcanic activity in Io is more complex, in theory volcanism can be driven by runaway melting in which the tidal heating increases as the mantle thickness decreases. We show that this runaway melting mechanism is generic for a composite planetary body with liquid core and solid mantle, provided that (i) the mantle rigidity, $\mu$, is comparable to the central pressure, i.e. $\mu/ (\rho g R_{\rm P})\gtrsim0.1$ for a body with density $\rho$, surface gravitational acceleration $g$, and radius $R_{\rm P}$, (ii) the surface is not molten, (iii) tides deposit sufficient energy, and (iv) the planet has nonzero eccentricity. We calculate the approximate liquid core radius as a function of $\mu/ (\rho g R_{\rm P})$, and find that more than $90\%$ of the core will melt due to this runaway for $\mu/ (\rho g R_{\rm P})\gtrsim1$. From all currently confirmed exoplanets, we find that the terrestrial planets in the L98-59 system are the most promising candidates for sustaining active volcanism. However, uncertainties regarding the quality factors and the details of tidal heating and cooling mechanisms prohibit definitive claims of volcanism on any of these planets. We generate synthetic transmission spectra of these planets assuming Venus-like atmospheric compositions with an additional 5, 50, and $98\%$ SO$_2$ component, which is a tracer of volcanic activity. We find a $\gtrsim 3 \sigma$ preference for a model with SO$_2$ with 5-10 transits with $\textit{JWST}$ for L98-59bcd.

The influence of intergalactic magnetic fields on the strong gravitational lensing of blazar secondary gamma radiation is discussed. Currently, two cases of strong gravitational lensing of blazar gamma-radiation are known, where radiation is deflected by galaxies on the line of sight between the blazars and the Earth. The magnetic field can affect the movements of electron-positron pairs generated by primary radiation, and thereby change the directions of secondary gamma radiation. It modifies the equation of the gravitational lens and leads to the dependence of the observed signal in the secondary gamma radiation on the energy of photons and on the magnetic field. Accordingly, it is possible in principle to estimate the intergalactic magnetic fields from the time delay of signals, from the angular position of images (for future high-resolution gamma-ray telescopes) or from the shape of the observed energy spectrum. This method is demonstrated by the example of the blazar B0218+357. In this case however it is not possible to obtain useful constraints due to the large distances to the blazar and the lens galaxy. The result is only a lower limit on the magnetic field $B>2\times10^{-17}$G, which is weaker than other existing constraints. But the future discoveries of lensed blazars may provide more favourable opportunities for measuring the magnetic fields, especially with the help of new generation of gamma-ray telescopes such as e-ASTROGAM, GECAM, and SVOM as well as future gamma-ray telescopes with high angular resolution, ~0.1''.

In this work, we focus on the period-luminosity relation (PLR) of $\delta$ Sct stars, in which mode excitation and selection mechanisms are still poorly constrained, and whose structure and oscillations are affected by rotation. We review the PLRs in the recent literature, and add a new inference from a large sample of $\delta$ Sct. We highlight the difficulty in identifying the fundamental mode and show that rotation-induced surface effects can impact the measured luminosities, explaining the broadening of the PLR. We derive a tight relation between the low-order large separation and the fundamental radial mode frequency (F0) that holds for rotating stars, thus paving the way towards mode identification. We show that the PLRs we obtain for different samples are compatible with each other and with the recent literature, and with most observed $\delta$ Sct stars when taking rotation effects into account. We also find that the highest-amplitude peak in the frequency spectrum corresponds to the fundamental mode in most $\delta$ Sct, thus shedding some light on their elusive mode selection mechanism.

AM CVn-type systems are ultracompact, helium-accreting binary systems which are evolutionarily linked to the progenitors of thermonuclear supernovae and are expected to be strong Galactic sources of gravitational waves detectable to upcoming space-based interferometers. AM CVn binaries with orbital periods $\lesssim$ 20--23 min exist in a constant high state with a permanently ionised accretion disc. We present the discovery of TIC 378898110, a bright ($G=14.3$ mag), nearby ($309.3 \pm 1.8$ pc), high-state AM CVn binary discovered in TESS two-minute-cadence photometry. At optical wavelengths this is the third-brightest AM CVn binary known. The photometry of the system shows a 23.07172(6) min periodicity, which is likely to be the `superhump' period and implies an orbital period in the range 22--23 min. There is no detectable spectroscopic variability. The system underwent an unusual, year-long brightening event during which the dominant photometric period changed to a shorter period (constrained to $20.5 \pm 2.0$ min), which we suggest may be evidence for the onset of disc-edge eclipses. The estimated mass transfer rate, $\log (\dot{M} / \mathrm{M_\odot} \mathrm{yr}^{-1}) = -6.8 \pm 1.0$, is unusually high and may suggest a high-mass or thermally inflated donor. The binary is detected as an X-ray source, with a flux of $9.2 ^{+4.2}_{-1.8} \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 0.3--10 keV range. TIC 378898110 is the shortest-period binary system discovered with TESS, and its large predicted gravitational-wave amplitude makes it a compelling verification binary for future space-based gravitational wave detectors.

We provide global models of line-driven winds of B supergiants for metallicities corresponding to the Large and Small Magellanic Clouds. The velocity and density structure of the models is determined consistently from hydrodynamical equations with radiative force derived in the comoving frame and level populations computed from kinetic equilibrium equations. We provide a formula expressing the predicted mass-loss rates in terms of stellar luminosity, effective temperature, and metallicity. Predicted wind mass-loss rates decrease with decreasing metallicity as $\dot M\sim Z^{0.60}$ and are proportional to the stellar luminosity. The mass-loss rates increase below the region of the bistability jump at about 20\,kK because of iron recombination. In agreement with previous theoretical and observational studies, we find a smooth change of wind properties in the region of the bistability jump. With decreasing metallicity, the bistability jump becomes weaker and shifts to lower effective temperatures. At lower metallicities above the bistability jump, our predictions provide similar rates to those used in current evolutionary models, but our rates are significantly lower than older predictions below the bistability jump. Our predicted mass-loss rates agree with observational estimates derived from H$\alpha$ line assuming that observations of stellar winds from Galaxy and the Magellanic Clouds are uniformly affected by clumping. The models nicely reproduce the dependence of terminal velocities on temperature derived from ultraviolet spectroscopy.

We use the first publically available data from the Daniel K. Inouye Solar Telescope (DKIST) to track magnetic connections from the solar photosphere into the corona. We scrutinize relationships between chromospheric magnetism and bright chromospheric, transition region and coronal plasmas. In June 2022, the Visible Spectro-Polarimeter (ViSP) instrument targeted unipolar network within a decaying active region. ViSP acquired rastered scans with longitudinal Zeeman sensitivities of 0.25 Mx/cm2 (Fe I 630.2 nm) and 0.5 Mx/cm2 (Ca II 854.2 nm). ViSP was operated in a "low" resolution mode (0.214" slit width, spectral resolution R ~ 70,000) to produce polarization maps over a common area of 105" x 50". Data from SDO and IRIS are combined to ask: Why is only a fraction of emerging flux filled with heated plasma? What is the elemental nature of the plasmas? No correlations were found between heated plasma and properties of chromospheric magnetic fields derived from the WFA, on scales below supergranules. Processes hidden from our observations control plasma heating. While improved magnetic measurements are needed, these data indicate that "the corona is a self-regulating forced system" (Einaudi et al. 2021). Heating depends on the state of the corona, not simply on boundary conditions. Heating models based upon identifiable bipolar fields, including cool loops, tectonics and observable magnetic reconnection, are refuted for these regions with unipolar chromospheric magnetic fields.

This paper describes the design of a 5.5:1 bandwidth feed antenna and reflector system, intended for use in hydrogen intensity mapping experiments. The system is optimized to reduce systematic effects that can arise in these experiments from scattering within the feed/reflector and cross-coupling between antennas. The proposed feed is an ultra wideband Vivaldi style design and was optimized to have a smooth frequency response, high gain, and minimal shadowing of the reflector dish. This feed can optionally include absorptive elements which reduce systematics but degrade sensitivity. The proposed reflector is a deep parabolic dish with $f/d = 0.216$ along with an elliptical collar to provide additional shielding. The procedure for optimizing these design choices is described.

This paper introduces a technique called NKL, which cleans both polarized and unpolarized foregrounds from HI intensity maps by applying a Karhunen-Lo\`eve transform on the needlet coefficients. In NKL, one takes advantage of correlations not only along the line of sight, but also between different angular regions, referred to as ``chunks". This provides a distinct advantage over many of the standard techniques applied to map-space that one finds in the literature, which do not consider such spatial correlations. Moreover, the NKL technique does not require any priors on the nature of the foregrounds, which is important when considering polarized foregrounds. We also introduce a modified version of GNILC, referred to as MGNILC, which incorporates an approximation of the foregrounds to improve performance. The NKL and MGNILC techniques are tested on simulated maps which include polarized foregrounds. Their performance is compared to the GNILC, GMCA, ICA and PCA techniques. Two separate tests were performed. One at $1.84 < z < 2.55$ and the other at $0.31 < z < 0.45$. NKL was found to provide the best performance in both tests, providing a factor of 10 to 50 improvement over GNILC at $k < 0.1\,{\rm hMpc^{-1}}$ in the higher redshift case and $k < 0.03 \,{\rm hMpc^{-1}}$ in the lower redshift case. However, none of the methods were found to recover the power spectrum satisfactorily at all BAO scales.

Atmospheric escape in exoplanets has traditionally been observed using hydrogen Lyman-$\alpha$ and H-$\alpha$ transmission spectroscopy, but more recent detections have utilised the metastable helium triplet at 1083$~$nm. Since this feature is accessible from the ground, it offers new possibilities for studying atmospheric escape. Our goal is to understand how the observability of escaping helium evolves during the lifetime of a highly irradiated gas giant. We extend our previous work on 1-D self-consistent hydrodynamic escape from hydrogen-only atmospheres as a function of planetary evolution to the first evolution-focused study of escaping hydrogen-helium atmospheres. Additionally, using these novel models we perform helium triplet transmission spectroscopy. We adapt our previous hydrodynamic escape model to now account for both hydrogen and helium heating and cooling processes and simultaneously solve for the population of helium in the triplet state. To account for the planetary evolution, we utilise evolving predictions of planetary radii for a close-in 0.3$~M_{\rm Jup}$ gas giant and its received stellar flux in X-ray, hard and soft EUV, and mid-UV wavelength bins assuming a K dwarf stellar host. We find that the helium triplet signature diminishes with evolution. Our models suggest that young ($\lesssim 150$~Myr), close-in gas giants ($\sim 1$ to $2~R_{\rm Jup}$) should produce helium 1083$~$nm transit absorptions of $\sim 4\%$ or $\sim 7\%$, for a slow or fast-rotating K dwarf, respectively, assuming a 2$\%$ helium abundance.

The reconstruction of the large scale velocity field from the grouped Cosmicflows-4 (CF4) database is presented. The lognormal bias of the inferred distances and velocities data is corrected by the Bias Gaussianization correction (BGc) scheme, and the linear density and velocity fields are reconstructed by means of the Wiener filter (WF) and constrained realizations (CRs) algorithm. These tools are tested against a suite of random and constrained Cosmicflows-3-like mock data. The CF4 data consists of 3 main subsamples - the 6dFGS and the SDSS data - and the `others'. The individual contributions of the subsamples have been studied. The quantitative analysis of the velocity field is done mostly by the mean overdensity ($\Delta_L(R)$) and the bulk velocity ($V_{\mathrm{bulk}}(R)$) profiles of the velocity field out to $300\, h^{-1}{\rm Mpc}$. The $V_{\mathrm{bulk}}(R)$ and $\Delta_{\mathrm L}(R)$ profiles of the CF4 data without its 6dFGS component are consistent with the cosmic variance to within $1\sigma$. The 6dFGS sample dominates the $V_{\mathrm{bulk}}$ ($\Delta_{\mathrm L}$) profile beyond $\sim120\, h^{-1}{\rm Mpc}$, and drives it to roughly a $3.4\sigma$ ($-1.9\sigma$) excess (deficiency) relative to the cosmic variance at $R\sim250\ (190)\ \, h^{-1}{\rm Mpc}$. The excess in the amplitude of $V_{\mathrm{bulk}}$ is dominated by its Supergalactic X component, roughly in the direction of the Shapley Concentration. The amplitude and alignment of the inferred velocity field from the CF4 data is at $\sim(2\,-\,3)\,\sigma$ discrepancy with respect to the $\Lambda$CDM model. Namely, it is somewhat atypical but yet there is no compelling tension with the model.

We propose polarization of scattered optical light from intermediate Galactic latitude infrared cirrus as a new diagnostic to constrain models of interstellar dust and the anisotropic interstellar radiation field (aISRF). For single scattering by a sphere, with Mie scattering phase functions for intensity and polarized intensity for a dust model at a given wavelength (Sloan $r$ and $g$ bands), and with models of anisotropic illumination from the entire sky (represented in HEALPix), we develop the formalism for calculating useful summary parameters for an integrated flux nebula (IFN): average of the phase function weighted by the illumination, polarization angle ($\psi$), and polarization fraction ($p$). To demonstrate the diagnostic discrimination of polarization from scattered light, we report on the effects of different anisotropic illumination models and different dust models on the summary parameters for the Spider IFN. The summary parameters are also sensitive to the IFN location, as we illustrate using FRaNKIE illumination models. For assessing the viability of dust and aISRF models, we find that observations of $\psi$ and $p$ of scattered light are indeed powerful new diagnostics to complement joint modeling of the intensity of scattered light (related to the average phase function) and the intensity of thermal dust emission. However, optically thin IFNs that can be modelled using single scattering are faint and $p$ is not large, as it could be with Rayleigh scattering, and so these observations need to be carried out with care and precision. Results for the Draco nebula compared to the Spider illustrate the challenge.

For the analysis of data taken by Imaging Air Cherenkov Telescopes (IACTs), a large number of air shower simulations are needed to derive the instrument response. The simulations are very complex, involving computational and memory-intensive calculations, and are usually performed repeatedly for different observation intervals to take into account the varying optical sensitivity of the instrument. The use of generative models based on deep neural networks offers the prospect for memory-efficient storing of huge simulation libraries and cost-effective generation of a large number of simulations in an extremely short time. In this work, we use Wasserstein Generative Adversarial Networks to generate photon showers for an IACT equipped with the FlashCam design, which has more than $1{,}500$ pixels. Using simulations of the H.E.S.S. experiment, we demonstrate the successful generation of high-quality IACT images. The analysis includes a comprehensive study of the generated image quality based on low-level observables and the well-known Hillas parameters that describe the shower shape. We demonstrate for the first time that the generated images have high fidelity with respect to low-level observables, the Hillas parameters, their physical properties, as well as their correlations. The found increase in generation speed in the order of $10^5$ yields promising prospects for fast and memory-efficient simulations of air showers for IACTs.

XO-7b is a hot Jupiter transiting a $V = 10.52$ mag G0V-type star. The planetary system is interesting because the linear slope in the discovery radial-velocity (RV) data indicated a wide-orbit massive companion. In 2020 we started an RV campaign for the system with the main scientific goal to follow-up this linear slope, and to put constraints on the orbital period of the companion. Furthermore, we aimed at refining the system parameters and we wanted to probe transit timing variations (TTVs) of XO-7b in order to search for long-term dynamical signs of the companion of XO-7b in the observed-minus-calculated (O-C) data of mid-transit times. Apart from the discovery RVs, we obtained and analyzed 20 follow-up RV observations and TESS photometric data. The previously observed significant linear RV slope was not confirmed with the follow-up RV data, where we detected only a marginal linear slope with the opposite trend. If the announced companion really exists, the most convincing explanation is that both RV datasets were collected near its quadrature position. Based on the RVs we estimated the minimum orbital period, which is $P_\mathrm{orb,min,3} \gtrsim 7900 \pm 1660$ d, and the 'minimum' minimum mass of the companion, which is $(M_3 \sin i)_\mathrm{min} = 16.7 \pm 3.5~\mathrm{M_{Jup}}$. We did not find significant evidence of the companion of XO-7b in the O-C dataset of mid-transit times. We can again conclude that if the announced companion really exists, this is in agreement with previous results that distant companions of exoplanets are only known by RV solutions.

Context. The corona of the Sun is the part of the solar atmosphere with temperatures of over one million Kelvin, which needs to be heated internally in order to exist. This heating mechanism remains a mystery; we see large magnetically active regions in the photosphere lead to strong extreme UV (EUV) emission in the corona. On much smaller scales (on the order of tens of Mm), there are bipolar and multipolar regions that can be associated with evenly sized coronal bright points (CBPs). Aims. Our aim was to study the properties of CBPs in a statistical sense and to use continuous data from the SDO spacecraft, which makes it possible to track CBPs over their whole lifetime. Furthermore, we tested various rotation-speed profiles for CBPs in order to find out if the lower corona is co-rotating with the photosphere. Then we compiled a database with about 346 CBPs together with information of their sizes, shapes, appearance and disappearance, and their visibility in the EUV channels of the AIA instrument. We want to verify our methods with similar previous studies. Methods. We used the high-cadence data of the largest continuous SDO observation interval in 2015 to employ an automated tracking algorithm for CBPs. Some of the information (e.g., the total lifetime, the characteristic shape, and the magnetic polarities below the CBPs) still requires human interaction. Results. In this work we present statistics on fundamental properties of CBPs along with some comparison tables that relate, for example, the CBP lifetime with their shape. CBPs that are visible in all AIA channels simultaneously seem to be brighter in total and also have a stronger heating, and hence a higher total radiation flux. We compared the EUV emission visibility in different AIA channels with the CBP's shape and lifetime. ... (full version see pdf)

Next generation high contrast imaging instruments face a challenging trade off: they will be required to deliver data with high spectral resolution at a relatively fast cadence (minutes) and across a wide field of view (arcseconds). For instruments that employ focal plane wavefront sensing and therefore require super-Nyquist sampling, these requirements cannot simultaneously be met with a traditional lenslet integral field spectrograph (IFU). For the SPIDERS pathfinder instrument, we are demonstrating an imaging Fourier transform spectrograph (IFTS) that offers a different set of tradeoffs than a lenslet IFU, delivering up to R20,000 spectral resolution across a dark hole. We present preliminary results from the SPIDERS IFTS including a chromaticity analysis of its dark hole and demonstrate a spectral differential imaging (SDI) improvement of up to 40 $\times$, and a first application of spectro-coherent differential imaging, combining both coherent differential imaging (CDI) and SDI.

Observations of 21\,cm line from neutral hydrogen promise to be an exciting new probe of astrophysics and cosmology during the Cosmic Dawn and through the Epoch of Reionization (EoR) to when dark energy accelerates the expansion of the Universe. At each of these epochs, separating bright foregrounds from the cosmological signal is a primary challenge that requires exquisite calibration. In this paper, we present a new calibration method called \textsc{nucal} that extends redundant-baseline calibration, allowing spectral variation in antenna responses to be solved for by using correlations between visibilities measuring the same angular Fourier modes at different frequencies. By modeling the chromaticity of the beam-weighted sky with a tunable set of discrete prolate spheroidal sequences (DPSS), we develop a calibration loop that optimizes for spectrally smooth calibrated visibilities. Crucially, this technique does not require explicit models of the sky or the primary beam. With simulations that incorporate realistic source and beam chromaticity, we show that this method solves for unsmooth bandpass features, exposes narrowband interference systematics, and suppresses smooth-spectrum foregrounds below the level of 21\,cm reionization models, even within much of the so-called "wedge" region where current foreground mitigation techniques struggle. We show that this foreground subtraction can be performed with minimal cosmological signal loss for certain well-sampled angular Fourier modes, making spectral-redundant calibration a promising technique for current and next-generation 21\,cm intensity mapping experiments.

A stable-frequency transmitter with relative radial acceleration to a receiver will show a change in received frequency over time, known as a "drift rate''. For a transmission from an exoplanet, we must account for multiple components of drift rate: the exoplanet's orbit and rotation, the Earth's orbit and rotation, and other contributions. Understanding the drift rate distribution produced by exoplanets relative to Earth, can a) help us constrain the range of drift rates to check in a Search for Extraterrestrial Intelligence (SETI) project to detect radio technosignatures and b) help us decide validity of signals-of-interest, as we can compare drifting signals with expected drift rates from the target star. In this paper, we modeled the drift rate distribution for $\sim$5300 confirmed exoplanets, using parameters from the NASA Exoplanet Archive (NEA). We find that confirmed exoplanets have drift rates such that 99\% of them fall within the $\pm$53 nHz range. This implies a distribution-informed maximum drift rate $\sim$4 times lower than previous work. To mitigate the observational biases inherent in the NEA, we also simulated an exoplanet population built to reduce these biases. The results suggest that, for a Kepler-like target star without known exoplanets, $\pm$0.44 nHz would be sufficient to account for 99\% of signals. This reduction in recommended maximum drift rate is partially due to inclination effects and bias towards short orbital periods in the NEA. These narrowed drift rate maxima will increase the efficiency of searches and save significant computational effort in future radio technosignature searches.

We reanalyse the Pantheon+ supernova catalogue to compare a cosmology with non-FLRW evolution, the "timescape cosmology", with the standard $\Lambda$CDM cosmology. To this end, we consider the Pantheon+ supernova catalogue, which is the largest available Type Ia supernova dataset for a geometric comparison between the two models. We construct a covariance matrix to be as independent of cosmology as possible, including independence from the FLRW geometry and peculiar velocity with respect to FLRW average evolution. Within this framework, which goes far beyond most other definitions of "model independence", we introduce new statistics to refine Type Ia supernova (SneIa) light-curve analysis. In addition to conventional galaxy correlation functions used to define the scale of statistical homogeneity we introduce empirical statistics which enables a refined analysis of the distribution biases of SneIa light-curve parameters $\beta c$ and $\alpha x_1$. For lower redshifts, the Bayesian analysis highlights important features attributable to the increased number of low-redshift supernovae, the artefacts of model-dependent light-curve fitting and the cosmic structure through which we observe supernovae. This indicates the need for cosmology-independent data reduction to conduct a stronger investigation of the emergence of statistical homogeneity and to compare alternative cosmologies in light of recent challenges to the standard model. "Dark energy" is generally invoked as a place-holder for "new physics". Our from-first-principles reanalysis of the Pantheon+ catalogue supports future deeper studies of the interplay of matter and nonlinear spacetime geometry, in a data-driven setting. For the first time in 25 years, we find evidence that the Pantheon+ catalogue already contains such a wealth of data that with further reanalysis, a genuine "paradigm shift" may soon emerge. [Abridged]

We present a novel method for identifying cosmic web filaments using the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. We compare the use of cosmic density field estimates from the Delaunay Tessellation Field Estimator (DTFE) and the Monte Carlo Physarum Machine (MCPM), which is inspired by the slime mold organism, in the DisPerSE structure identification framework. The MCPM-based reconstruction identifies filaments with higher fidelity, finding more low-prominence/diffuse filaments and better tracing the true underlying matter distribution than the DTFE-based reconstruction. Using our new filament catalogs, we find that most galaxies are located within 1.5-2.5 Mpc of a filamentary spine, with little change in the median specific star formation rate and the median galactic gas fraction with distance to the nearest filament. Instead, we introduce the filament line density, {\Sigma}fil(MCPM), as the total MCPM overdensity per unit length of a local filament segment, and find that this parameter is a superior predictor of galactic gas supply and quenching. Our results indicate that most galaxies are quenched and gas-poor near high-line density filaments at z<=1. At z=0, quenching in log(M*/Msun)>10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. In high-line density filaments, satellites are strongly quenched, whereas centrals have reduced star formation, but not gas fraction, at z<=0.5. We discuss the prospect of applying our new filament identification method to galaxy surveys with SDSS, DESI, Subaru PFS, etc. to elucidate the effect of large-scale structure on galaxy formation.

We investigate how the photon polarization is affected by the interaction with axion-like particles (ALPs) in the rotating magnetic field of a neutron star (NS). Using quantum Boltzmann equations the study demonstrates that the periodic magnetic field of millisecond NSs enhances the interaction of photons with ALPs and creates a circular polarization on them. A binary system including an NS and a companion star could serve as a probe. When the NS is in front of the companion star with respect to the earth observer, there is a circular polarization on the previously linearly polarized photons as a result of the interaction with ALPs there. After a half-binary period, the companion star passes in front of the NS, and the circular polarization of photons disappears and changes to linear. The excluded parameter space for a millisecond NS with 300~Hz rotating frequency, highlights the coupling constant of $1.7\times10^{-11}~\text{GeV}^{-1}\leq g_{a\gamma\gamma}\leq1.6\times10^{-3}~\text{GeV}^{-1}$ for the ALP masses in the range of $7\times10^{-12}~\text{eV}\leq m_a\leq1.5\times 10^{3}~\text{eV}$.

This paper compiles the model parameters and zero-temperature properties of an extensive collection of published theoretical nuclear interactions, including 251 non-relativistic (Skyrme-like), 252 relativistic mean field (RMF) and point-coupling (RMF-PC), and 13 Gogny-like forces. This forms the most exhaustive tabulation of model parameters to date. The properties of uniform symmetric matter and pure neutron matter at the saturation density are determined. Symmetry properties found from the second-order term of a Taylor expansion in neutron excess are compared to the energy difference of pure neutron and symmetric nuclear matter at the saturation density. Selected liquid-droplet model parameters, including the surface tension and surface symmetry energy, are determined for semi-infinite surfaces. Liquid droplet model neutron skin thicknesses and dipole polarizabilities of the neutron-rich closed-shell nuclei $^{48}$Ca and $^{208}$Pb are compared to published theoretical Hartree-Fock and experimental results. In addition, the radii, binding energies, moments of inertia and tidal deformabilities of 1.2, 1.4 and 1.6 M$_\odot$ neutron stars are computed. An extensive correlation analysis of bulk matter, nuclear structure, and low-mass neutron star properties is performed and compared to nuclear experiments and astrophysical observations.

Wash-in leptogenesis is a powerful mechanism to generate the baryon asymmetry of the Universe that treats right-handed-neutrino interactions on the same footing as electroweak sphaleron processes: as mere spectator processes acting on the background of chemical potentials in the Standard Model plasma. Successful wash-in leptogenesis requires this chemical background to be CP-violating, which can be achieved by violating any of the more than ten global charges that are conserved in the Standard Model at very high temperatures. In this paper, we demonstrate that the primordial charge asymmetries required for wash-in leptogenesis can be readily produced by evaporating primordial black holes (PBHs). Our argument is based on the fact that the Hawking radiation emitted by PBHs contains more or less any state in the particle spectrum. Therefore, if heavy states with CP-violating decays are present in the ultraviolet, PBH evaporation will unavoidably lead to the production of these states. We illustrate this scenario by means of a simple toy model where PBH evaporation leads to the production of heavy particles that we call asymmetrons and whose decay results in a primordial charge asymmetry for right-handed electrons, which in turn sets the initial conditions for wash-in leptogenesis. We focus on the parameter region where the decay of the initial thermal asymmetron abundance occurs long before PBH evaporation and only results in a negligible primordial charge asymmetry. PBH evaporation at later times then serves as a mechanism to resurrect the asymmetron abundance and ensure the successful generation of the baryon asymmetry after all. We conclude that PBHs can act as asymmetry-producing machines that grant access to whatever CP-violating physics may be present in the ultraviolet, rekindling it at lower energies where it can be reprocessed into a baryon asymmetry by right-handed neutrinos.

Given the urgency to reduce fossil fuel energy production to make climate tipping points less likely, we call for resource-aware knowledge gain in the research areas on Universe and Matter with emphasis on the digital transformation. A portfolio of measures is described in detail and then summarized according to the timescales required for their implementation. The measures will both contribute to sustainable research and accelerate scientific progress through increased awareness of resource usage. This work is based on a three-days workshop on sustainability in digital transformation held in May 2023.

We study the effect of temperature on the global properties of static and slowly rotating self-gravitating Bose-Einstein condensate (BEC) stars within general relativity. We employ a recently developed temperature dependent BEC equation of state (EoS) to describe the stellar matter by assuming that the condensate can be described by a non-relativistic EoS. Stellar profiles are obtained using general relativistic Hartle-Thorne slow rotation approximation equations. We find that with increasing temperatures mass-radius values are found to be decreasing for the static and rotating cases; though presence of temperature supports high mass values at lower central densities. Countering effects of rotation and temperature on the BEC stellar structure have been analysed and quantified. We report that inclusion of temperature has significant effect on the rotating stellar profiles but negligible effect on the maximum mass, as in the case of static system. We have also studied the effect of EoS parameters -- boson mass and strength of the self-interaction -- on global properties of static and rotating BEC stars, in presence of temperature.

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.

Recent progress has revealed a number of constraints that cosmological correlators and the closely related field-theoretic wavefunction must obey as a consequence of unitarity, locality, causality and the choice of initial state. When combined with symmetries, namely homogeneity, isotropy and scale invariance, these constraints enable one to compute large classes of simple observables, an approach known as (boostless) cosmological bootstrap. Here we show that it is possible to relax the restriction of scale invariance, if one retains a discrete scaling subgroup. We find an infinite class of solutions to the weaker bootstrap constraints and show that they reproduce and extend resonant non-Gaussianity, which arises in well-motivated models such as axion monodromy inflation. We find no evidence of the new non-Gaussian shapes in the Planck data. Intriguingly, our results can be re-interpreted as a deformation of the scale-invariant case to include a complex order of the total energy pole, or more evocatively interactions with a complex number of derivatives. We also discuss for the first time IR-divergent resonant contributions and highlight an inconsequential inconsistency in the previous literature.