Papers for Friday, Sep 08 2023

International LOFAR stations are powerful radio telescopes, however they are delivered without the tooling necessary to convert their raw data stream into standard data formats that can be used by common processing pipelines, or science-ready data products. udpPacketManager is a C and C++ library that was developed with the intent of providing a faster-than-realtime software package for converting raw data into arbitrary data formats based on the needs of observers working with the Irish LOFAR station (I-LOFAR), and stations across Europe. It currently offers an open-source solution for both offline and online (pre-)processing of telescope data into a wide variety of formats.

By combining the JWST/NIRCam JADES and CEERS extragalactic datasets, we have uncovered a sample of twenty-one T and Y brown dwarf candidates at best-fit distances between 0.1 - 4.2 kpc. These sources were selected by targeting the blue 1$\mu$m - 2.5$\mu$m colors and red 3$\mu$m - 4.5$\mu$m colors that arise from molecular absorption in the atmospheres of T$_{\mathrm{eff}} < $ 1300K brown dwarfs. We fit these sources using multiple models of low-mass stellar atmospheres and present the resulting fluxes, sizes, effective temperatures and other derived properties for the sample. If confirmed, these fits place the majority of the sources in the Milky Way thick disk and halo. We observe proper motion for seven of the candidate brown dwarfs with directions in agreement with the plane of our galaxy, providing evidence that they are not extragalactic in nature. We demonstrate how the colors of these sources differ from selected high-redshift galaxies, and explore the selection of these sources in planned large-area JWST NIRCam surveys. Deep imaging with JWST/NIRCam presents an an excellent opportunity for finding and understanding these very cold low-mass stars at kpc distances.

Post-merger galaxies are unique laboratories to study the triggering and interplay of star-formation and AGN activity. Combining new, high resolution, 10 GHz Jansky Very Large Array (VLA) observations with archival radio surveys, we have examined the radio properties of 28 spheroidal post-merger galaxies. We find a general lack of extended emission at (sub-)kiloparsec scales, indicating the prevalence of compact, nuclear radio emission in these post-merger galaxies, with the majority (16/18; 89\%) being radio-quiet at 10 GHz. Using multi-wavelength data, we determine the origin of the radio emission, discovering 14 new radio AGN and 4 post-mergers dominated by emission from a population of supernova remnants. Among the radio AGN, almost all are radio-quiet (12/14; 86\%). We discover a new dual AGN (DAGN) candidate, J1511+0417, and investigate the radio properties of the DAGN candidate J0843+3549. 4 of these radio AGN are hosted by SF emission-line galaxies, suggesting that radio AGN activity may be present during periods of SF activity in post-mergers. The low jet powers and compact morphologies of these radio AGN also point to a scenario in which AGN feedback may be efficient in this sample of post-mergers. Lastly, we present simulated, multi-frequency observations of the 14 radio AGN with the Very Long Baseline Array (VLBA) and the VLBI capabilities of the Next Generation Very Large Array (ngVLA) to assess the feasibility of these instruments in searches for supermassive black hole binaries (SMBHBs).

Observations with the Five-Hundred-Meter Aperture Spherical Telescope have revealed the presence of a marginally-resolved source of 21 cm emission from a location $\sim50'$ from the M94 galaxy, without a stellar counterpart down to the surface brightness limit of the DESI Imaging Legacy Survey ($\sim29.15$ mag arcsec$^{-2}$ in the $g$ band). The system (hereafter Cloud-9) has round column density isocontours and a line width consistent with thermal broadening from gas at $T\sim2\times10^4$ $K$. These properties are unlike those of previously detected dark HI clouds and similar to the expected properties of REionization-Limited-HI Cloud (RELHICs), namely, starless dark matter (DM) halos filled with gas in hydrostatic equilibrium and in thermal equilibrium with the cosmic ultraviolet background. At the distance of M94, $d\sim4.7$ Mpc, we find that Cloud-9 is consistent with being a RELHIC inhabiting a Navarro-Frenk-White (NFW) DM halo of mass, $M_{200}\sim5\times10^{9}$ $M_{\odot}$, and concentration, $c_{\rm NFW}\sim13$. Although the agreement between the model and observations is good, Cloud-9 appears to be slightly, but systematically, more extended than expected for $\Lambda$CDM RELHICs. This may imply either that Cloud-9 is much closer than implied by its recessional velocity, $v_{\rm CL9}\sim300$ km s$^{-1}$, or that its halo density profile is flatter than NFW, with a DM mass deficit greater than a factor of $10$ at radii $r\lesssim1$ kpc. Further observations may aid in constraining these scenarios better and help elucidate whether Cloud-9 is the first ever observed RELHIC, a cornerstone prediction of the $\Lambda$CDM model on the smallest scales.

The 3.5 keV line is a purported emission line observed in galaxies, galaxy clusters, and the Milky Way whose origin is inconsistent with known atomic transitions and has previously been suggested to arise from dark matter decay. We systematically re-examine the bulk of the evidence for the 3.5 keV line, attempting to reproduce six previous analyses that found evidence for the line. Surprisingly, we only reproduce one of the analyses; in the other five we find no significant evidence for a 3.5 keV line when following the described analysis procedures on the original data sets. For example, previous results claimed 4$\sigma$ evidence for a 3.5 keV line from the Perseus cluster; we dispute this claim, finding no evidence for a 3.5 keV line. We find evidence for background mismodeling in multiple analyses. We show that analyzing these data in narrower energy windows diminishes the effects of mismodeling but returns no evidence for a 3.5 keV line. We conclude that there is little robust evidence for the existence of the 3.5 keV line. Some of the discrepancy of our results from those of the original works may be due to the earlier reliance on local optimizers, which we demonstrate can lead to incorrect results. For ease of reproducibility, all code and data are publicly available.

Polaris is the nearest and brightest classical Cepheid, and pulsates with a period of about 4 days. It has long been known as a single-lined spectroscopic binary with an orbital period of 30 yr. Historical photometric and spectroscopic records indicate that, until recently, the pulsation period has been increasing at a rate of about 4.5 s/yr, and that the amplitude of the pulsation declined for most of the 20th century, but more recently halted its decline and began to increase. Here we report an analysis of the more than 3600 individual radial velocity measurements of Polaris available from the literature over the past 126 yr. We find that the pulsation period is now becoming shorter, and that the amplitude of the velocity variations has stopped increasing, and may be getting smaller again. We also find tantalising evidence that these changes in pulsation behaviour over the last century may be related to the binary nature of the system, as they seem to occur near each periastron passage, when the secondary comes within 29 stellar radii of the Cepheid in its eccentric orbit. This suggests the companion may be perturbing the atmosphere of the Cepheid and altering its pulsation properties at each encounter. After removal of the pulsation component of the velocities, we derive a much improved spectroscopic orbit for the binary that should serve as the basis for a more accurate determination of the dynamical masses, which are still rather uncertain.

With the rapid advance of wide-field surveys it is increasingly important to perform combined cosmological probe analyses. We present a new pipeline for simulation-based multi-probe analyses, which combines tomographic large-scale structure (LSS) probes (weak lensing and galaxy clustering) with cosmic microwave background (CMB) primary and lensing data. These are combined at the $C_\ell$-level, yielding 12 distinct auto- and cross-correlations. The pipeline is based on $\texttt{UFalconv2}$, a framework to generate fast, self-consistent map-level realizations of cosmological probes from input lightcones, which is applied to the $\texttt{CosmoGridV1}$ N-body simulation suite. It includes a non-Gaussian simulation-based covariance, several data compression schemes, and a neural network emulator for accelerated theoretical predictions. We validate our framework, apply it to a simulated $12\times2$pt tomographic analysis of KiDS, BOSS, and $\textit{Planck}$, and forecast constraints for a $\Lambda$CDM model with a variable neutrino mass. We find that, while the neutrino mass constraints are driven by the CMB data, the addition of LSS data helps to break degeneracies and improves the constraint by up to 35%. For a fiducial $M_\nu=0.15\mathrm{eV}$, a full combination of the above CMB+LSS data would enable a $3\sigma$ constraint on the neutrino mass. We explore data compression schemes and find that MOPED outperforms PCA. We also study the impact of an internal lensing tension, parametrized by $A_L$, on the neutrino mass constraint, finding that the addition of LSS to CMB data including all cross-correlations is able to mitigate the impact of this systematic. $\texttt{UFalconv2}$ and a MOPED compressed $\textit{Planck}$ CMB primary + CMB lensing likelihood are made publicly available. [abridged]

The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest halos as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyze the catalogs of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass halos with $M_{200} \sim 10^{11}$ M$_{\odot}$. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive halos is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intra-cluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than $10$ GC tracers are particularly prone to severe uncertainties.Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive halos suggested by GC kinematics around some UDGs

Milky Way Cepheid variables with accurate {\it Hubble Space Telescope} photometry have been established as standards for primary calibration of the cosmic distance ladder to achieve a percent-level determination of the Hubble constant ($H_0$). These 75 Cepheid standards are the fundamental sample for investigation of possible residual systematics in the local $H_0$ determination due to metallicity effects on their period-luminosity relations. We obtained new high-resolution ($R\sim81,000$), high signal-to-noise ($S/N\sim50-150$) multi-epoch spectra of 42 out of 75 Cepheid standards using ESPaDOnS instrument at the 3.6-m Canada-France-Hawaii Telescope. Our spectroscopic metallicity measurements are in good agreement with the literature values with systematic differences up to $0.1$ dex due to different metallicity scales. We homogenized and updated the spectroscopic metallicities of all 75 Milky Way Cepheid standards and derived their multiwavelength ($GVIJHK_s$) period-luminosity-metallicity and period-Wesenheit-metallicity relations using the latest {\it Gaia} parallaxes. The metallicity coefficients of these empirically calibrated relations exhibit large uncertainties due to low statistics and a narrow metallicity range ($\Delta\textrm{[Fe/H]}=0.6$~dex). These metallicity coefficients are up to three times better constrained if we include Cepheids in the Large Magellanic Cloud and range between $-0.21\pm0.07$ and $-0.43\pm0.06$ mag/dex. The updated spectroscopic metallicities of these Milky Way Cepheid standards were used in the Cepheid-Supernovae distance ladder formalism to determine $H_0=72.9~\pm 1.0$\textrm{~km~s$^{-1}$~Mpc$^{-1}$}, suggesting little variation ($\sim 0.1$ ~km~s$^{-1}$~Mpc$^{-1}$) in the local $H_0$ measurements due to different Cepheid metallicity scales.

The SPT-3G 2018 TT/TE/EE cosmic microwave background (CMB) data set (temperature and polarization) is used to place constraints on an axion-like model of early dark energy (EDE). These data do not favor axion-like EDE and place an upper limit on the maximum fraction of the total energy density $f_{\rm EDE}< 0.172$ (at the 95% confidence level, CL). This is in contrast with ACT DR4 which gives $f_{\rm EDE}=0.150^{+0.050}_{-0.078}$. When combining CMB measurements with measurements of the baryon acoustic oscillations and luminosity distance to Type Ia supernovae, we show that the tension with the S$H_0$ES measurement of the Hubble parameter goes up from 2.6$\sigma$ with Planck to 2.9$\sigma$ with Planck+SPT-3G 2018. The additional inclusion of ACT DR4 data leads to a reduction of the tension to $1.6\sigma$, but the discrepancy between ACT DR4 and Planck+SPT-3G 2018 casts some doubt on the statistical consistency of this joint analysis. The importance of improved measurements of the CMB at both intermediate and small scales (in particular the shape of the damping tail) as well as the interplay between temperature and polarization measurements in constraining EDE are discussed. Upcoming ground-based measurements of the CMB will play a crucial role in determining whether EDE remains a viable model to address the Hubble tension.

We present detailed and comprehensive data reduction and point-spread-function (PSF) model construction for all public JWST NIRCam imaging data from the COSMOS-Web treasury program (up to June 2023, totaling 0.28 ${\rm deg}^2$). We show that the NIRCam PSF has significant short-timescale temporal variations and random spatial variations in all four filters (F115W, F150W, F277W, and F444W). Combining NIRCam with archival HST imaging, we perform multiwavelength AGN+host image decomposition to study the properties of 143 X-ray-selected ($L_{\rm bol}=10^{43.6-47.2}$ erg s$^{-1}$) broad-line AGNs at $0.35\lesssim z \lesssim 3.5$. Leveraging the superb resolution, wavelength coverage, and sensitivity of NIRCam, we successfully detect host stellar emission after decomposing the central AGN point source in 142 objects. $\sim 2/3$ AGNs are in star-forming galaxies based on the UVJ diagram, suggesting no instantaneous negative AGN feedback. X-ray-selected broad-line AGN hosts follow a similar stellar mass-size relation as inactive galaxies, albeit with slightly smaller galaxy sizes. We find that although major mergers are rare ($\sim$7-22%) among the sample, more subtle non-axisymmetric features from stellar bars, spiral arms, and minor mergers are ubiquitous, highlighting the importance of secular processes and minor mergers in triggering AGN activity. For a subsample of 30 AGNs at $1<z<2.5$ with black hole mass measurements from single epoch spectra, they follow a similar black hole mass-stellar mass relation as local inactive early-type galaxies but reside preferentially near the upper envelope of nearby AGNs. We caution that selection biases and intrinsic differences of AGN populations at different redshifts may significantly affect their location on the black hole mass-stellar mass plane.

Magnetic reconnection is a ubiquitous phenomenon for magnetized plasmas and leads to the rapid reconfiguration of magnetic field lines. During reconnection events, plasma is heated and accelerated until the magnetic field lines enclose and capture the plasma within a circular configuration. These plasmoids could therefore observationally manifest themselves as hot spots that are associated with flaring behavior in supermassive black hole systems, such as Sagittarius A$^\ast$. We have developed a novel algorithm for identifying plasmoid structures, which incorporates watershed and custom closed contouring steps. From the identified plasmoids, we determine the plasma characteristics and energetics in magnetohydrodynamical simulations. The algorithm's performance is showcased for a high-resolution suite of axisymmetric ideal and resistive magnetohydrodynamical simulations of turbulent accretion discs surrounding a supermassive black hole. For validation purposes, we also evaluate several Harris current sheets that are well-investigated in the literature. Interestingly, we recover the characteristic power-law distribution of plasmoid sizes for both the black hole and Harris sheet simulations. This indicates that while the dynamics are vastly different, with different dominant plasma instabilities, the plasmoid creation behavior is similar. Plasmoid occurrence rates for resistive general relativistic magnetohydrodynamical simulations are significantly higher than for the ideal counterpart. Moreover, the largest identified plasmoids are consistent with sizes typically assumed for semi-analytical interpretation of observations. We recover a positive correlation between the plasmoid formation rate and a decrease in black-hole-horizon-penetrating magnetic flux. The developed algorithm has enabled an extensive quantitative analysis of plasmoid formation in black hole accretion simulations.

The recent launch of JWST has ushered in a new era of high-redshift astronomy by providing detailed insights into the gas and stellar populations of galaxies in the epoch of reionization. Interpreting these observations and translating them into constraints on the physics of early galaxy formation is a complex challenge that requires sophisticated models of star formation and the interstellar medium (ISM) in high-redshift galaxies. To this end, we present Version 1 of the Sphinx$^{20}$ public data release. Sphinx$^{20}$ is a full box cosmological radiation hydrodynamics simulation that simultaneously models the large-scale process of cosmic reionization and the detailed physics of a multiphase ISM, providing a statistical sample of galaxies akin to those currently being observed by JWST. The data set contains $\sim14,000$ mock images and spectra of the stellar continuum, nebular continuum, and 52 nebular emission lines, including Ly$\alpha$, for each galaxy in Sphinx$^{20}$ with a star formation rate $\geq0.3\ {\rm M_{\odot}\ yr^{-1}}$. All galaxy emission has been processed with dust radiative transfer and/or resonant line radiative transfer, and data is provided for ten viewing angles for each galaxy. Additionally, we provide a comprehensive set of intrinsic galaxy properties, including halo masses, stellar masses, star formation histories, and ISM characteristics (e.g., metallicity, ISM gas densities, LyC escape fractions). This paper outlines the data generation methods, presents a comparative analysis with JWST ERS and Cycle 1 observations, and addresses data set limitations. The Sphinx$^{20}$ data release can be downloaded at the following URL: https://github.com/HarleyKatz/SPHINX-20-data

In November 2020, the Swift team announced a major update to the calibration of the UltraViolet and Optical Telescope (UVOT) data to correct for the gradual loss of sensitivity over time. Beginning in roughly 2015, the correction affected observations in the three near ultraviolet (UV) filters, reaching levels of up to 0.3 mag immediately prior to the correction. Over the same time period, an increasing number of Type I superluminous supernovae (SLSNe-I) were discovered and studied. Many SLSNe-I are hot (T$_\textrm{eff}$ $\approx 10,000$ K) near peak, and therefore accurate UV data are imperative towards properly understanding their physical properties and energetics. We re-compute Swift UVOT photometry for SLSNe-I discovered between 2014 and 2021 with at least 5 Swift observations in 2015 or later. We calculate host-subtracted magnitudes for each SLSN and fit their spectral energy distributions with modified blackbodies to obtain the radius and temperature evolution. We also fit multi-band photometry using the Modular Open Source Fitter for Transients (MOSFiT) to obtain key parameters such as the spin period (P), magnetic field strength (B), ejecta mass (M$_\textrm{ej}$), and kinetic energy (E$_\textrm{kin}$). From our MOSFiT modeling, we also estimate the peak UV/optical luminosity (L$_\textrm{peak}$) and total radiative energy (E$_\textrm{rad}$). Under the assumption of magnetar-powered SLSNe we find several strong trends, including anti-correlations between P and both L$_\textrm{peak}$ and E$_\textrm{rad}$, a correlation between E$_\textrm{kin}$ and E$_\textrm{rad}$, and an anti-correlation between B and E$_\textrm{rad}$.

We report the spectroscopic analysis of 20 halo ab-type RR Lyrae stars with heliocentric distances between 15 and 165 kpc, conducted using medium-resolution spectra from the Magellan Inamori Kyocera Echelle (MIKE) spectrograph. We obtain the systemic line-of-sight velocities of our targets with typical uncertainties of 5-10 km s$^{-1}$, and compute orbital parameters for a subsample out to 50 kpc from the Galactic centre, including proper motion data from Gaia DR3. The orientation of our stars' orbits, determined for an isolated Milky Way and for a model perturbed by the Large Magellanic Cloud, appears to suggest an accreted origin for at least half of the sample. In addition, we derive atmospheric parameters and chemical abundance ratios for seven stars beyond 20 kpc. The derived $\alpha$-abundances of five of these stars follow a Milky Way halo-like trend, while the other two display an underabundance of $\alpha$-elements for their [Fe/H], indicating an association with accretion events. Furthermore, based on the [Sr/Ba] ratio, we can speculate about the conditions for the formation of a potential chemically peculiar carbon-enhanced metal-poor (CEMP) RR Lyrae star. By analysing the stars' orbital parameters and abundance ratios, we find hints of association of two of our stars with two massive satellites, namely the Large Magellanic Cloud and Sagittarius. Overall, our results are in line with the suggestion that the accretion of sub-haloes largely contributes to the outer halo stellar populations.

We present a study aimed at understanding the physical phenomena underlying the formation and evolution of galaxies following a data-driven analysis of spectroscopic data based on the variance in a carefully selected sample. We apply Principal Component Analysis (PCA) independently to three subsets of continuum-subtracted optical spectra, segregated into their nebular emission activity as quiescent, star-forming, and Active Galactic Nuclei (AGN). We emphasize that the variance of the input data in this work only relates to the absorption lines in the photospheres of the stellar populations. The sample is taken from the Sloan Digital Sky Survey (SDSS) in the stellar velocity dispersion range 100-150 km/s, to minimise the ``blurring'' effect of the stellar motion. We restrict the analysis to the first three principal components (PCs), and find that PCA segregates the three types with the highest variance mapping SSP-equivalent age, along with an inextricable degeneracy with metallicity, even when all three PCs are included. Spectral fitting shows that stellar age dominates PC1, whereas PC2 and PC3 have a mixed dependence of age and metallicity. The trends support - independently of any model fitting - the hypothesis of an evolutionary sequence from star-formation to AGN to quiescence. As a further test of the consistency of the analysis, we apply the same methodology in different spectral windows, finding similar trends, but the variance is maximal in the blue wavelength range, roughly around the 4000A break.

Accreting white dwarfs are known to show signatures of wind-type outflows in the ultraviolet. At optical wavelengths, however, wind detections have only been reported for a few sources. We present GTC-10.4m optical spectroscopy of four accreting white dwarfs (BZ Cam, V751 Cyg, MV Lyr, and V425 Cas) observed during luminous epochs, when their optical emission is expected to be dominated by the accretion disc. We focused the analysis on four emission lines: H$\alpha$ and He I $\lambda$5876, $\lambda$6678, $\lambda$7065. Line profiles are complex and variable on short (minutes) and long (days to weeks) time scales, with transient absorption and emission components. Among them, we detect strong blue-shifted absorptions at $\gtrsim 1000$ km s$^{-1}$. These high-velocity components, present only in the blue wing of the emission lines, are observed in all four sources and could be associated with accretion disc winds. For MV Lyr and V425 Cas, these would represent the first detection of optical outflows in these objects, while in the case of BZ Cam and V751 Cyg, the presence of outflows has been previously reported. This study suggests that, in addition to ultraviolet winds, optical outflows might be also common in accreting white dwarfs. We discuss the observational properties of these winds and their possible similarity to those detected in accreting black holes and neutrons stars.

Despite the importance of Active Galactic Nuclei (AGN) in galaxy evolution, the mechanisms that fuel AGN activity remain poorly understood. Theoretical models suggest that major mergers of galaxies contribute strongly to AGN fuelling, particularly at high AGN luminosities. The connection between mergers and AGN activity has therefore been widely studied, although with contradictory results. Some studies find a strong connection between mergers and AGN, while others find merger fractions in AGN hosts to match those in the inactive galaxy population. To address these apparent contradictions, I present a complete and systematic analysis of detected merger fractions in AGN hosts from the literature. I assess if discrepancies between studies are indicative of systematic uncertainties and biases and analyse the detected merger fraction as a function of luminosity, redshift, and AGN selection method. X-ray selected AGN samples show comparable detected merger fractions across studies and major mergers do not dominate triggering in this AGN population. On the other hand, signatures of significant merger contribution to the AGN population are observed in a small fraction of primarily radio selected and reddened AGN samples. It is unclear if this is due to observational biases or physical differences in the host galaxies. There is no correlation between the detected merger fraction and AGN luminosity. This lack of correlation between detected merger fraction and AGN luminosity, which has previously been reported in the literature, cannot be explained by systematic uncertainties and observational biases.

We studied seven nearby narrow-line Seyfert 1 (NLS1) galaxies in $J$- and $Ks$-bands with redshifts varying from 0.019 to 0.092. This is the first multi-source study targeting the hosts of southern NLS1 galaxies. Our data was obtained with the FourStar instrument of the 6.5~m Magellan Baade telescope at the Las Campanas Observatory (Chile). The aim of our study is to determine the host galaxy morphologies of these sources by using GALFIT. We were able to model six out of the seven sources reliably. Our conclusion is that all of the reliably modelled sources are disk-like galaxies, either spirals or lenticulars. None of these sources present elliptical morphology. Our findings are in agreement with the hypothesis that disk-like galaxies are the main host of jetted NLS1 galaxies. Taking advantage of observations in two bands, we also produced a $J - Ks$ colour map of each source. Five of the six colour maps show significant dust extinction near the core of the galaxy -- a feature often seen in gamma-ray-detected jetted NLS1 galaxies, and interpreted to be a consequence of a past minor merger.

We investigate the X-ray variability properties of Seyfert1 Galaxies belonging to the BAT AGN Spectroscopic Survey (BASS). The sample includes 151 unobscured (N$_{\rm H}<10^{22}$ cm$^{-2}$) AGNs observed with XMM-Newton for a total exposure time of ~27 Ms, representing the deepest variability study done so far with high signal-to-noise XMM-Newton observations, almost doubling the number of observations analysed in previous works. We constrain the relation between the normalised excess variance and the 2-10 keV AGN luminosities, black hole masses and Eddington ratios. We find a highly significant correlation between $\sigma^{2}_{NXS}$ and $M_{\rm BH}$, with a scatter of ~0.85 dex. For sources with high $L_{2-10}$ this correlation has a lower normalization, confirming that more luminous (higher mass) AGNs show less variability. We explored the $\sigma^{2}_{NXS}$ vs $M_{\rm BH}$ relation for the sub-sample of sources with $M_{\rm BH}$ estimated via the "reverberation mapping" technique, finding a tighter anti-correlation, with a scatter of ~ 0.65 dex. We examine how the $\sigma^{2}_{NXS}$ changes with energy by studying the relation between the variability in the hard (3-10 keV) and the soft (0.2-1 keV)/medium (1-3 keV) energy bands, finding that the spectral components dominating the hard energy band are more variable than the spectral components dominating in softer energy bands, on timescales shorter than 10 ks.

We explore the dependence of global galaxy properties in the SIMBA simulation as a function of distance from filaments identified using DisPerSE. We exclude halos with mass $M_h>10^{13}M_\odot$ to mitigate the impact of group and cluster environments. Galaxies near filaments are more massive and have more satellites, which we control for by examining deviations from best-fit scaling relations. At $z=0$, star formation (SF) is significantly suppressed within $\lesssim 100$ kpc of filaments, more strongly for satellites, indicating substantial pre-processing in filaments. By $z=2$, the trend is weak and if anything indicates an increase in SF activity close to filaments. The suppression at $z\lesssim 1$ is accompanied by lowered HI fractions, and increased metallicities, quenched fractions, and dispersion-dominated systems. H2 fractions are not strongly suppressed when controlling for stellar mass, suggesting that star formation efficiency drives the drop in SF. By comparing amongst different SIMBA feedback variant runs, we show that the majority of SF suppression owes to filamentary shock-heating, but there is a non-trivial additional effect from AGN feedback. When looking around massive ($M_h>10^{13}M_\odot$) halos, those galaxies near filaments behave somewhat differently, indicating that filaments provide an additional environmental effect relative to halos. Finally, we compare SIMBA results to EAGLE and IllustrisTNG at $z=0$, showing that all models predict SF suppression within $\lesssim 100$ kpc of filaments, nonetheless, detailed differences may be observationally testable.

Can new cosmic physics be uncovered through tensions amongst datasets? Tensions in parameter determinations amongst different types of cosmological observation, especially the `Hubble tension' between probes of the expansion rate, have been invoked as possible indicators of new physics, requiring extension of the $\Lambda$CDM paradigm to resolve. Within a fully Bayesian framework, we show that the standard tension metric gives only part of the updating of model probabilities, supplying a data co-dependence term that must be combined with the Bayes factors of individual datasets. This shows that, on its own, a reduction of dataset tension under an extension to $\Lambda$CDM is insufficient to demonstrate that the extended model is favoured. Any analysis that claims evidence for new physics {\it solely} on the basis of alleviating dataset tensions should be considered incomplete and suspect. We describe the implications of our results for the interpretation of the Hubble tension.

Within the self-gravitating Bose-Einstein condensate (BEC) model of dark matter (DM), we argue that the axion-like self-interaction of ultralight bosons ensures the existence of both rarefied and dense phases in the DM halo core of dwarf galaxies. In fact, this stems from two independent solutions of the Gross-Pitaevskii equation corresponding to the same model parameters. For a small number of particles, this structure disappears, and the Gross-Pitaevskii equation reduces to stationary sine-Gordon equation, the one-dimensional anti-kink solution of which mimics a single-phase DM radial distribution in the halo core. Quantum mechanically, this solution corresponds to a zero-energy bound state of two particles in a closed scattering channel formed by the domain-wall potential with a finite asymptotics. To produce a two-particle composite with low positive energy and a finite lifetime, we appeal to the resonant transition of one asymptotically free particle of a pair from an open channel (with a model scattering potential) to the closed channel. Using the Feshbach resonance concept, the problem of two-channel quantum mechanics is solved in the presence of a small external influence which couples the two channels, and an analytical solution is obtained in the first approximation. Analyzing the dependence of scattering data on interaction parameters, we reveal a long-lived two-particle composite (dimer) state possessing a lifetime of millions of years. This result is rather surprising and supposes important implications of dimers' being involved in forming large DM structures. It is shown that the dimers' appearance is related with the regime of infinite scattering length due to resonance. The revealed dependence of the DM scattering length $a$ on the parameters of interactions can theoretically justify variation of $a$ in the DM dominated galaxies and its role for large DM structures.

We present a novel AI approach for high-resolution high-dynamic range synthesis imaging by radio interferometry (RI) in astronomy. R2D2, standing for "{R}esidual-to-{R}esidual {D}NN series for high-{D}ynamic range imaging", is a model-based data-driven approach relying on hybrid deep neural networks (DNNs) and data-consistency updates. Its reconstruction is built as a series of residual images estimated as the outputs of DNNs, each taking the residual dirty image of the previous iteration as an input. The approach can be interpreted as a learned version of a matching pursuit approach, whereby model components are iteratively identified from residual dirty images, and of which CLEAN is a well-known example. We propose two variants of the R2D2 model, built upon two distinctive DNN architectures: a standard U-Net, and a novel unrolled architecture. We demonstrate their use for monochromatic intensity imaging on highly-sensitive observations of the radio galaxy Cygnus~A at S band, from the Very Large Array (VLA). R2D2 is validated against CLEAN and the recent RI algorithms AIRI and uSARA, which respectively inject a learned implicit regularization and an advanced handcrafted sparsity-based regularization into the RI data. With only few terms in its series, the R2D2 model is able to deliver high-precision imaging, significantly superior to CLEAN and matching the precision of AIRI and uSARA. In terms of computational efficiency, R2D2 runs at a fraction of the cost of AIRI and uSARA, and is also faster than CLEAN, opening the door to real-time precision imaging in RI.

The physical collisions of two white dwarfs (WDs) (i.e. not slow mergers) have been shown to produce type-Ia-like supernovae (SNe) explosions. Most studies of WD-collisions have focused on zero impact-parameter (direct) collisions, which can also be studied in 2D. However, the vast majority of WD collisions arising from any evolutionary channels suggested to date, are expected to be indirect, i.e. have a non-negligible impact parameter upon collision. Here we use the highest resolution 3D simulations to date (making use of the AREPO code), in order to explore both direct and indirect collisions and the conditions in which they give rise to a detonation and the production of a luminous SNe. Using our simulations, we find a detonation criterion that can provide the critical impact parameter for an explosion to occur, depending on the density profile of the colliding WDs, their composition, and their collision velocities. We find that the initial velocity has a significant impact on the amount of 56Ni production from the explosion. Furthermore, the 56Ni production is also strongly dependent on the numerical resolution. While our results from the head-on collision with large initial velocities produce more 56Ni than previous simulations, those with small and comparable velocities produce significantly less 56Ni than previous works.

The effect of dust attenuation on a galaxy's light depends on a number of physical properties, such as geometry and dust composition, both of which can vary across the faces of galaxies. To investigate this variation, we continue analysis on star-forming regions in 29 galaxies studied previously. We analyse these regions using Swift/UVOT and WISE images, as well as SDSS/MaNGA emission line maps to constrain the relationship between the infrared excess (IRX) and the UV spectral index (beta) for each star forming region. This relationship can be used to constrain which dust attenuation law is appropriate for the region. We find that the value of Dn(4000) for a region is correlated with both IRX and beta, and that the gas-phase metallicity is strongly correlated with the IRX. This correlation between metallicity and IRX suggests that regardless of aperture, metal rich regions have steeper attenuation curves. We also find that integrated galactic light follows nearly the same IRX-beta relationship as that found for kiloparsec-sized star forming regions. This similarity may suggest that the attenuation law followed by the galaxy is essentially the same as that followed by the regions, although the relatively large size of our star forming regions complicates this interpretation because optical opacity and attenuation curves have been observed to vary within individual galaxies.

We investigate several panchromatic scaling relations (SRs) for the dwarf irregular galaxy NGC 1569 using IFU data from the Metal-THINGS Survey. Among the spatially resolved properties analyzed, we explore SRs between the stellar mass, SFR, molecular gas, total gas, baryonic mass, gas metallicity, gas fraction, SFE and effective oxygen yields. Such multiwavelength SRs are analyzed at a spatial resolution of 180 pc, by combining our IFU observations with data from the surveys THINGS, CARMA, and archival data from DustPedia. Although we recover several known relations, our slopes are different to previously reported ones. Our star formation main sequence, Kennicutt-Schmidt (KS) and molecular KS relations show higher SFRs, lower scatter, and higher correlations, with steeper (1.21), and flatter slopes (0.96, 0.58) respectively. The shape of the SRs including metallicity, stellar mass, and gas fraction are flat, with an average value of 12+log(O/H) $\sim$ 8.12 dex. The baryonic mass vs effective oxygen yields, and the stellar, gas and baryonic mass vs SFE show higher dispersions and lower correlations. Since we use the dust mass as a tracer of gas mass, we derive the Dust-to-Gas Ratio and the CO luminosity-to-molecular gas mass conversion factors, showing differences of 0.16 and 0.95 dex for the total and molecular gas surface density, respectively, in comparison to previously reported values. We use a self regulated feedback model to conclude that stellar feedback plays an important role generating outflows in NGC 1569.

With neutron star applications in mind, we developed a theory of diffusion in mixtures of superfluid, strongly interacting Fermi liquids. By employing the Landau theory of Fermi liquids, we determined matrices that relate the currents of different particle species, their momentum densities, and the partial entropy currents to each other. Using these results, and applying the quasiclassical kinetic equation for the Bogoliubov excitations, we derived general expressions for the diffusion coefficients, which properly incorporate all the Fermi liquid effects and depend on the momentum transfer rates between different particle species. The developed framework can be used as a starting point for systematic calculations of the diffusion coefficients (as well as other kinetic coefficients) in superfluid Fermi mixtures, particularly, in superfluid neutron stars.

We present the results of our analysis of Gaia19dke, an extraordinary microlensing event in the Cygnus constellation that was first spotted by the {\gaia} satellite. This event featured a strong microlensing parallax effect, which resulted in multiple peaks in the light curve. We conducted extensive photometric, spectroscopic, and high-resolution imaging follow-up observations to determine the mass and the nature of the invisible lensing object. Using the Milky Way priors on density and velocity of lenses, we found that the dark lens is likely to be located at a distance of $D_L =(3.05^{+4.10}_{-2.42})$kpc, and has a mass of $M_L =(0.51^{+3.07}_{-0.40}) M_\odot$. Based on its low luminosity and mass, we propose that the lens in Gaia19dke event is an isolated white dwarf.

Studies of young clusters have shown that a large fraction of O-/early B-type stars are in binary systems, where the binary fraction increases with mass. These massive stars are present in clusters of a few Myrs, but gradually disappear for older clusters. The lack of detailed studies of intermediate-age clusters has meant that almost no information is available on the multiplicity properties of stars with M $<$ 4 $M_{\odot}$. In this study we present the first characterization of the binary content of NGC 1850, a 100 Myr-old massive star cluster in the Large Magellanic Cloud, relying on a VLT/MUSE multi-epoch spectroscopic campaign. By sampling stars down to M = 2.5 $M_{\odot}$, we derive a close binary fraction of 24 $\pm$ 5 \% in NGC 1850, in good agreement with the multiplicity frequency predicted for stars of this mass range. We also find a trend with stellar mass (magnitude), with higher mass (brighter) stars having higher binary fractions. We modeled the radial velocity curves of individual binaries using The Joker and constrained the orbital properties of 27 systems, $\sim$17\% of all binaries with reliable radial velocities in NGC 1850. This study has brought to light a number of interesting objects, such as four binaries showing mass functions f(M) $>$ 1.25 $M_{\odot}$. One of these, star #47, has a peculiar spectrum, explainable with the presence of two disks in the system, around the visible star and the dark companion, which is a black hole candidate. These results confirm the importance and urgency of studying the binary content of clusters of any age.

Recently, the Lyman-$\alpha$ (Ly$\alpha$) forest flux auto-correlation function has been shown to be sensitive to the mean free path of hydrogen-ionizing photons, $\lambda_{\text{mfp}}$, for simulations at $z \geq 5.4$. Measuring $\lambda_{\text{mfp}}$ at these redshifts will give vital information on the ending of reionization. Here we present the first observational measurements of the Ly$\alpha$ forest flux auto-correlation functions in ten redshift bins from $5.1 \leq z \leq 6.0$. We use a sample of 35 quasar sightlines at $z > 5.7$ from the extended XQR-30 data set, this data has signal-to-noise ratios of $> 20$ per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly$\alpha$ systems. With these measurements, we introduce software tools to generate auto-correlation function measurements from any simulation. For an initial comparison, we show our auto-correlation measurements with simulation models for recently measured $\lambda_{\text{mfp}}$ values and find good agreements. Further work in modeling and understanding the covariance matrices of the data is necessary to get robust measurements of $\lambda_{\text{mfp}}$ from this data.

Clusters of galaxies possess the capability to accelerate cosmic rays (CRs) to very high energy up to $\sim10^{18}$~eV due to their large size and magnetic field strength which favor CR confinement for cosmological times. During their confinement, they can produce neutrinos and $\gamma-$rays out of interactions with the background gas and photon fields. In recent work, \cite{hussain2021high, hussain2023diffuse} have conducted three-dimensional cosmological magnetohydrodynamical (MHD) simulations of the turbulent intracluster medium (ICM) combined with multi-dimensional Monte Carlo simulations of CR propagation for redshifts ranging from $z \sim 5$ to $z = 0$ to study the multi-messenger emission from these sources. They found that when CRs with a spectral index in the range $1.5 - 2.5$ and cutoff energy $E_\mathrm{max} = 10^{16} - 10^{17}$~eV are injected into the system, they make significant contributions to the diffuse background emission of both neutrinos and gamma-rays. In this work, we have revisited this model and undertaken further constraints on the parametric space. This was achieved by incorporating the recently established upper limits on neutrino emission from galaxy clusters, as obtained by the IceCube experiment. We find that for CRs injected with spectral indices in the range $2.0 - 2.5$, cutoff energy $E_\mathrm{max} = 10^{16} - 10^{17}$~eV, and power corresponding to $(0.1-1)\%$ of the cluster luminosity, our neutrino flux aligns with the upper limits estimated by IceCube. Additionally, the resulting contribution from clusters to the diffuse $\gamma$-ray background (DGRB) remains significant with values of the order of $ \sim 10^{-5}\, \mathrm{MeV} \, \mathrm{cm}^{-2} \,\mathrm{s}^{-1} \, \mathrm{sr}^{-1}$ at energies above $500$ GeV.

Stellar flares cannot be spatially resolved, which complicates ascertaining the physical processes behind particular spectral signatures. Due to their proximity to Earth, solar flares can serve as a stepping stone for understanding their stellar counterparts, especially when using a Sun-as-a-star instrument and in combination with spatially resolved observations. We aim to understand the disk-integrated spectral behaviors of a confined X2.2 solar flare and its eruptive X9.3 successor as measured by HARPS-N. The behavior of multiple photospheric and chromospheric spectral lines are investigated by means of activity indices and contrast profiles. A number of different photospheric lines were also investigated by means of equivalent widths, and radial velocity measures, which are then related to physical processes directly observed in high-resolution observations made with the Swedish 1-meter Solar Telescope and SDO Our findings suggest a relationship between the evolving shapes of contrast profile time and the flare locations, which assists in constraining flare locations in disk-integrated observations. In addition, an upward bias was found in flare statistics based on activity indices derived from the Ca II H & K lines. In this case, much smaller flares cause a similar increase in the activity index as that produced by larger flares. H$\alpha$-based activity indices do not show this bias and are therefore less susceptible to activity jitter. Sodium line profiles show a strongly asymmetric response during flare activity, which is best captured with a newly defined asymmetrical sodium activity index. A strong flare response was detected in Mn I line profiles, which is unexpected and calls for further exploration. Intensity increases in H$\alpha$, H$\beta$, and certain spectral windows of AIA before the flare onset suggest their potential use as short-term flare predictors.

Context. Most stars from in groups which with time disperse, building the field population of their host galaxy. In the Milky Way, open clusters have been continuously forming in the disk up to the present time, providing it with stars spanning a broad range of ages and masses. Observations of the details of cluster dissolution are, however, scarce. One of the main difficulties is obtaining a detailed characterisation of the internal cluster kinematics, which requires very high quality proper motions. For open clusters, which are typically loose groups with some tens to hundreds of members, there is the additional difficulty of inferring kinematic structures from sparse and irregular distributions of stars. Aims. Here, we aim to analyse internal stellar kinematics of open clusters, and identify rotation, expansion or contraction patterns. Methods. We use Gaia Early Data Release 3 (EDR3) astrometry and Integrated Nested Laplace Approximations to perform vector-field inference and create spatio-kinematic maps of 1237 open clusters. The sample is composed of clusters for which individual stellar memberships were known, thus minimising contamination from field stars in the velocity maps. Projection effects were corrected using EDR3 data complemented with radial velocities from Gaia Data Release 2 and other surveys. Results. We report the detection of rotation patterns in 8 open clusters. Nine additional clusters display possible rotation signs. We also observe 14 expanding clusters, with 15 other objects showing possible expansion patterns. Contraction is evident in two clusters, with one additional cluster presenting a more uncertain detection. In total, 53 clusters are found to display kinematic structures. Within these, elongated spatial distributions suggesting tidal tails are found in 5 clusters. [abridged]

We present here the DECam Ecliptic Exploration Project (DEEP), a three year NOAO/NOIRLab Survey that was allocated 46.5 nights to discover and measure the properties of thousands of trans-Neptunian objects (TNOs) to magnitudes as faint as VR~27, corresponding to sizes as small as 20 km diameter. In this paper we present the science goals of this project, the experimental design of our survey, and a technical demonstration of our approach. The core of our project is "digital tracking," in which all collected images are combined at a range of motion vectors to detect unknown TNOs that are fainter than the single exposure depth of VR~23 mag. Through this approach we reach a depth that is approximately 2.5 magnitudes fainter than the standard LSST "wide fast deep" nominal survey depth of 24.5 mag. DEEP will more than double the number of known TNOs with observational arcs of 24 hours or more, and increase by a factor of 10 or more the number of known small (<50 km) TNOs. We also describe our ancillary science goals, including measuring the mean shape distribution of very small main belt asteroids, and briefly outline a set of forthcoming papers that present further aspects of and preliminary results from the DEEP program.

Recent independent observations using several different telescope systems an analysis methods have provided evidence of parity violation between the number of galaxies that spin in opposite directions. On the other hand, other studies argued that no parity violation can be identified. This paper provides detailed analysis, statistical inference, and reproduction of previous reports that show no preferred spin direction. Code and data used for the reproduction are publicly available. The results show that the data used in all of these studies agrees with the observation of a preferred direction as observed from Earth. In some of these studies the datasets were too small, or the statistical analysis was incomplete. In other papers the results were impacted by experimental design decisions that lead directly to show non-preferred direction. In some of these cases these decisions are not stated in the papers, but were revealed after further investigation in cases where the reproduction of the work did not match the results reported in the papers. These results show that the data used in all of these previous studies in fact agree with the contention that galaxies as observed from Earth have a preferred spin direction, and the distribution of galaxy spin directions as observed from Earth form a cosmological-scale dipole axis. This study also shows that the reason for the observations is not necessarily an anomaly in the large-scale structure, and can also be related to internal structure of galaxies.

Broad absorption line (BAL) quasars are characterized by gas clouds that absorb flux at the wavelength of common quasar spectral features, although blueshifted by velocities that can exceed 0.1c. BAL features are interesting as signatures of significant feedback, yet they can also compromise cosmological studies with quasars through their impact on accurate redshifts and measurements of the matter density distribution traced by the Lyman-alpha forest. The presence of BALs can also significantly contaminate the shape of the most prominent quasar emission lines and introduce systematic shifts in quasar redshifts. We present a catalog of BAL quasars discovered in the Dark Energy Spectroscopic Instrument (DESI) survey Early Data Release, which were observed as part of DESI Survey Validation, as well as the first two months of the main survey. We describe our method to automatically identify BAL quasars in DESI data, the quantities we measure for each BAL, and investigate the completeness and purity of this method with mock DESI observations. We mask the wavelengths of the BAL features and recompute the quasar redshifts, and find the new redshifts differ by 243 km/s on average for the BAL quasar sample. These new, more accurate redshifts are important to obtain the best measurements of quasar clustering, especially at small scales. Finally, we present some spectra of rarer classes of BALs that illustrate the potential of DESI data to identify such populations for further study.

The surface temperature of an asteroid is fundamental information for the design of an exploration mission and the interpretation of scientific observations. In addition, the thermal radiation of the asteroid causes a non-gravitational acceleration that induces secular changes in its orbit and spin. We have been developing a numerical calculation library for simulating the dynamics and thermophysics of asteroids. The asteroid dynamical simulator, \texttt{Astroshaper}, can calculate the temperature distribution based on a 3-dimensional shape model of an asteroid and predict the non-gravitational acceleration. In recent years, asteroid exploration missions such as Hayabusa2 and Hera have been equipped with thermal infrared imagers. The asteroid thermography can provide the thermal properties of the surface material of the target bodies. The functionality of thermophysical modeling in \texttt{Astroshaper} contributes to simulating the thermal environment on the asteroids, estimating the thermal properties, and predicting the dynamical evolution controlled by the non-gravitational effects.

We introduce \texttt{Weakhub}, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino-matter, neutrino-photon, and neutrino-neutrino interactions, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of $\beta$-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino-antineutrino annihilation, and nuclear de-excitation are included. We also present the detailed derivations of weak interactions and the coupling of them to the two-moment based general-relativistic multi-group radiation transport in the \texttt{G}eneral-relativistic \texttt{mu}ltigrid \texttt{nu}merical (\texttt{Gmunu}) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernova and binary neutron star postmerger remnant, and predict the effects of improved opacities in comparison to conventional ones for a binary neutron star postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino library \texttt{NuLib} in a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than $\sim 10\%$ in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modelling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations.

Context. The study of planet-crossing asteroids is of both practical and fundamental importance. As they are closer than asteroids in the Main Belt, we have access to a smaller size range, and this population frequently impacts planetary surfaces and can pose a threat to life. Aims. We aim to characterize the compositions of a large corpus of planet-crossing asteroids and to study how these compositions are related to orbital and physical parameters. Methods. We gathered publicly available visible colors of near-Earth objects (NEOs) from the Sloan Digital Sky Survey (SDSS) and SkyMapper surveys. We also computed SDSS-compatible colors from reflectance spectra of the Gaia mission and a compilation of ground-based observations. We determined the taxonomy of each NEO from its colors and studied the distribution of the taxonomic classes and spectral slope against the orbital parameters and diameter. Results. We provide updated photometry for 470 NEOs from the SDSS, and taxonomic classification of 7,401 NEOs. We classify 42 NEOs that are mission-accessible, including six of the seven flyby candidates of the ESA Hera mission. We confirm the perihelion dependance of spectral slope among S-type NEOs, likely related to a rejuvenation mechanism linked with thermal fatigue. We also confirm the clustering of A-type NEOs around 1.5-2 AU, and predict the taxonomic distribution of small asteroids in the NEO source regions in the Main Belt.

The Rosette Nebula is a young stellar cluster and molecular cloud complex, located at the edge of the southern shell of a middle-aged SNR Monoceros Loop (G205.5+0.5). We revisited the GeV gamma-ray emission towards the Rosette Nebula using more than 13 years of Fermi-LAT data. We tested several spatial models and found that compared to the result using the CO gas template only, the inclusion of the HII gas template can significantly improve the likelihood fit. We performed spectral analysis using the new spatial template. With both the gamma-ray observation and CO+HII gas data, we derived the cosmic ray spectrum of different components in the vicinity of the Rosette Nebula. We found the gamma-ray emissions from Rosette Nebula are substantially harder than previously reported, which may imply that Rosette Nebula is another example of a gamma-ray emitting young massive star cluster.

We present measurements of the relation between X-ray luminosity and star formation activity for a sample of normal galaxies spanning the redshift range between 0 and 0.25. We use data acquired by SRG/eROSITA for the performance and verification phase program called eROSITA Final Equatorial Depth Survey (eFEDS). The eFEDS galaxies are observed in the 0.2-2.3 keV band. Making use of a wide range of ancillary data, spanning from the ultraviolet (UV) to mid-infrared wavelengths (MIR), we estimated the star formation rate (SFR) and stellar mass ($M_{star}$) of 888 galaxies, using Code Investigating GALaxy Emission (CIGALE). We divided our sample of normal galaxies in star-forming (SFGs) and quiescent galaxies according to their position on the main sequence. We confirm a linear correlation between the X-ray luminosity and the SFR for our sample of SFGs, as shown previously in the literature. However, we find this relation to be strongly biased by the completeness limit of the eFEDS survey. Correcting for completeness, we find the fitted relation to be consistent with the literature. We also investigated the relation between X-ray emission from both LMXBs and HMXBs populations with $M_{star}$ and SFR, respectively. Correcting for completeness, we find our fitted relation to considerably scatter from the literature relation at high specific SFR ($SFR/M_{star}$). We conclude that without accounting for X-ray non-detections, it is not possible to employ eFEDS data to study the redshift evolution of the LMXBs and HMXBs contributions due to completeness issues. Furthermore, we find our sources to largely scatter from the expected Lx/SFR vs specific SFR relation at high redshift. We discuss the dependence of the scatter on the stellar mass, metallicity, or the globular cluster content of the galaxy.

With the advent of multi-wavelength electromagnetic observations of neutron stars, spanning many decades in photon energies, from radio wavelengths up to X-rays and $\gamma$-rays, it becomes possible to significantly constrain the geometry and the location of the associated emission regions. In this work, we use results from the modelling of thermal X-ray observations of PSR~J0030+0451 from the NICER mission and phase-aligned radio and $\gamma$-ray pulse profiles to constrain the geometry of an off-centred dipole able to reproduce the light-curves in these respective bands simultaneously. To this aim, we deduce a configuration with a simple dipole off-centred from the location of the centre of the thermal X-ray hot spots and show that the geometry is compatible with independent constraints from radio and $\gamma$-ray pulsations only, leading to a fixed magnetic obliquity of $\alpha \approx 75\deg$ and a line of sight inclination angle of $\zeta \approx 54\deg$. We demonstrate that an off-centred dipole cannot be rejected by accounting for the thermal X-ray pulse profiles. Moreover, the crescent shape of one spot is interpreted as the consequence of a small scale surface dipole on top of the large scale off-centred dipole.

The BL Lacertae object OJ 287 is a very unusual quasar producing a wobbling radio jet and some double-peaked optical outbursts with a possible period of about 12 yr for more than one century. This variability is widely explained by models of binary supermassive black hole (SMBH) or precessing jet/disk from a single SMBH. To enable an independent and nearly bias-free investigation on these possible scenarios, we explored the feasibility of extremely high-precision differential astrometry on its innermost restless jet at mm-wavelengths. Through re-visiting some existing radio surveys and very long baseline interferometry (VLBI) data at frequencies from 1.4 to 15.4 GHz and performing new Very Long Baseline Array (VLBA) observations at 43.2 GHz, we find that the radio source J0854$+$1959, 7.1 arcmin apart from OJ 287 and no clearly-seen optical and infrared counterparts, could provide a nearly ideal reference point to track the complicated jet activity of OJ 287. The source J0854$+$1959 has a stable GHz-peaked radio spectrum and shows a jet structure consisting of two discrete, mas-scale-compact and steep-spectrum components and showing no proper motion over about 8 yr. The stable VLBI structure can be interpreted by an episodic, optically thin and one-sided jet. With respect to its 4.1-mJy peak feature at 43.2 GHz, we have achieved an astrometric precision at the state-of-art level, about 10 $\mu$as. These results indicate that future VLBI astrometry on OJ 287 could allow us to accurately locate its jet apex and activity boundary, align its restless jet structure over decades without significant systematic bias, and probe various astrophysical scenarios.

Study of the tilt angles of solar bipolar magnetic regions is important because the tilt angles have an important role in the solar dynamo. We analysed the data on tilt angles of sunspot groups measured at the Mt. Wilson Observatory (MWOB) during the period 1917-1986 and Kodaikanal Observatory (KOB) during the period 1906-1986. We binned the daily tilt-angle data during each of the Solar Cycles 15-21 into different 5-deg. latitude intervals and calculated the mean value of the tilt angles in each latitude interval and the corresponding standard error. We fitted these binned data to Joy's Law (increase of the tilt angle with latitude), i.e. the linear relationship between tilt angle and latitude of an active region. The linear-least-square fit calculations were done by taking into account the uncertainties in both the abscissa (latitude) and ordinate (mean tilt angle). The calculations were carried out by using both the tilt-angle and area weighted tilt-angle data in the whole sphere, northern hemisphere, and southern hemisphere during the whole period and during each individual solar cycle. We find a significant difference (absolute north--south asymmetry) between the slopes of Joy's Law recovered from northern and southern hemispheres' whole period MWOB data of area-weighted tilt angles. Only the slope obtained from the southern hemisphere's MWOB data of a solar cycle is found to be reasonably well anti-correlated to the amplitude of the solar cycle. In the case of area weighted tilt-angle data, a good correlation is found between the absolute north--south asymmetry in the slope of a solar cycle and the amplitude of the solar cycle. The corresponding best-fit linear equations are found to be statistically significant.

We present the timing analysis of the intermediate polar TX Col in the X-ray band using the observations made by Chandra, Swift, and Suzaku during the years 2000, 2007, and 2009, respectively. The spin, orbital, and beat periods derived from these data are consistent with the earlier findings. We found that the spin modulation was dominant during the Chandra observation, whereas both orbital and beat modulations were dominant during the Swift and Suzaku observations. These findings and past X-ray observations indicate that TX Col is changing its accretion geometry from disc dominance to stream dominance and vice versa.

The open problems related to cosmological tensions in current times have opened new paths to study new probes to constrain cosmological parameters in standard and extended cosmologies, in particular, to determine at a local level the value of the Hubble constant $H_0$, through independent techniques. However, while standard Cosmological Constant Cold Dark Matter ($\Lambda$CDM) model has been well constrained and parts of extended cosmology have been intensively studied, the physics behind them aspects restrains our possibilities of selecting the best cosmological model that can show a significant difference from the first model. Therefore, to explore a possible deviation from a such model that can explain the current discrepancy on the $H_0$ value, in this work we consider adding the current local observables, e.g. Supernovae Type Ia (SNIa), $H(z)$ measurements, and Baryon Acoustic Observations (BAO) combined with two new calibrated Quasars (QSO) datasets using ultraviolet, x-ray and optical plane techniques. While these can be identified as part of the high-redshift standard candle objects, the main characteristics of these are based on fluxes distributions calibrated up to $z \sim 7 $. We consider five $H_0$ prior scenarios to develop these calibrations. Furthermore, we found that our estimations provide the possibility to relax the $H_0$ tension at 2$\sigma$ using a QSO ultraviolet sample in combination with late measurements showing higher values of $H_0$. Our results can be an initial start for more serious treatments in the quasars physics from ultraviolet, x-ray, and optical plane techniques behind the local observations as cosmological probes to relax the cosmological tensions problems.

The S-index is a measure of emission in the CaII H & K lines and is a widely used proxy of stellar magnetic activity. It has been assumed until now that the S-index is mainly affected by bright plage regions in the chromosphere. In particular, the effect of starspots on the S-index has been neglected. In this study we revisit this assumption. For this we analyze high-resolution observations of sunspots recorded in the CaII H spectral line at the Swedish 1-m Solar Telescope and determine the contrast of spots with respect to the quiet surroundings. We find that the CaII H line core averaged over whole sunspots (including superpenumbrae) is brighter than in the quiet surroundings and that the spot contrast in the line core is comparable to the facular contrast. This allows us to get a first estimate of the influence of spots on the S-index. We show that spots increase the S-index. While this increase is quite small for the Sun, it becomes significantly larger for more active stars. Further, we show that the inclusion of the contribution of spots to the S-index strongly affects the relationship between the S-index and stellar disk area coverages by spots and faculae, and present the new relations.

Gravitational wave (GW) events, particularly those connected to the merger of compact objects such as neutron stars, are believed to be the primary source of short gamma-ray bursts. To explore the very high energy (VHE) component of the emission from these events, the H.E.S.S. collaboration has dedicated a substantial effort and observing time to follow up on these events. During the second and third GW observing runs, H.E.S.S. was the first ground-based instrument to observe the GW170817 binary neutron star merger. In addition, H.E.S.S. followed four binary black hole mergers. The data acquired by H.E.S.S. was used to constrain the VHE emission from these events for the first time. H.E.S.S. also monitored the GW170817 source for approximately 50 hours and obtained limits that constrained the magnetic field in the merger remnant to $> 24 \mu G$. As the fourth GW observing run (O4) approaches, the H.E.S.S. collaboration has allocated significant observation time to the follow-up of GW events. This contribution provides an overview of the science results derived from the H.E.S.S. follow-up of GW events, a technical overview of the GW follow-up strategies for O4, and an update on H.E.S.S. activities during O4.

This contribution describes some recent advances in the parallelization of the generation and processing of radio signals emitted by particle showers in CORSIKA 8. CORSIKA 8 is a Monte Carlo simulation framework for modeling ultra-high energy particle cascades in astroparticle physics. The aspects associated with the generation and processing of radio signals in antennas arrays are reviewed, focusing on the key design opportunities and constraints for deployment of multiple threads on such calculations. The audience is also introduced to Gyges, a lightweight, header-only and flexible multithread self-adaptive scheduler written compliant with C++17 and C++20, which is used to distribute and manage the worker computer threads during the parallel calculations. Finally, performance and scalability measurements are provided and the integration into CORSIKA 8 is commented.

The binary stellar system $\eta$ Carinae is one of very few established astrophysical hadron accelerators. It seems likely that at least some fraction of the accelerated particles escape from the system. Copious target material for hadronic interactions and associated $\gamma$-ray emission exists on a wide range of spatial scales outside the binary system. This material creates a unique opportunity to trace the propagation of particles into the interstellar medium. Here we analyse $\gamma$-ray data from Fermi-LAT of $\eta$ Carinae and surrounding molecular clouds and investigate the many different scales on which escaping particles may interact and produce $\gamma$-rays. We find that interactions of escaping cosmic rays from $\eta$ Carinae in the wind region and the Homunculus Nebula could produce a significant contribution to the $\gamma$-ray emission associated with the system. Furthermore, we detect excess emission from the surrounding molecular clouds. The derived radial cosmic-ray excess profile is consistent with a steady injection of cosmic rays by a central source. However, this would require a higher flux of escaping cosmic rays from $\eta$ Carinae than provided by our model. Therefore it is likely that additional cosmic ray sources contribute to the hadronic $\gamma$-ray emission from the clouds.

Many important details on the mechanisms underlying the ejection of material during a (classical) nova eruption are still not understood. Here we present optical spectroscopy and narrow-band images of the nova V1425 Aql, 23 years after the nova eruption. We find that the ejecta consist of two significantly different components. The first resembles what is commonly seen in novae: a symmetric distribution, centred on the position of the underlying cataclysmic binary, and presenting both allowed (hydrogen and helium) and forbidden ([OIII] and [NII]) transitions. The second one, however, consists of approximately three times higher velocity material that is not visible in the allowed transitions, presents a significantly different [NII] - [OIII] ratio and is located at approximately 2.3 arcsec to the south-west of the position of the binary. Comparing the velocities and spatial extensions of the two ejecta, we find that both originated in the same nova eruption. We explore possible extrinsic and intrinsic mechanisms for the asymmetry of the high-velocity material in the form of asymmetrically distributed interstellar material and magnetic accretion, respectively, but find the available data to be inconclusive. From the expansion parallax, we derive the distance of the nova to 3.3(3) kpc.

In the environments where the abundance of heavy elements is low, rare events are expected to impact the chemical enrichment. Dwarf galaxies have small masses, low average metallicities and in general low star formation rates, and thus investigating the chemical enrichment provides understanding on the impact of each source of elements on the chemical abundance. Using a chemical evolution model in which the rarity is introduced, we investigate the impact of rare events on the chemical enrichment for Local Group dwarf galaxies. In the model, the occurrence of individual sources of elements is estimated with the star formation history derived by the colour-magnitude diagram. The abundance ratios of trans-iron elements to iron predicted by the model show the oscillation at the lowest metallicities because of the r-process events. In the case of a galaxy of a lower mass, the oscillation caused by neutron star mergers is also seen at higher metallicities, which suggests that the rarity can be important in lower-mass systems. Regarding the source of the chemical enrichment, we observe that the r-process sites seem to contribute more to the production of trans-iron elements at low metallicities, but massive stars of different rotating velocities also contribute to create part of the dispersion of the abundance ratios through the s-process. Both observational and theoretical data, including nucleosynthesis calculations and the chemical abundance of metal-poor stars, are needed to obtain deeper insights into the sources of the chemical enrichment at low metallicities.

Galactic conformity is the phenomenon in which a galaxy of a certain physical property is correlated with its neighbors of the same property, implying a possible causal relationship. The observed auto correlations of emission line galaxies (ELGs) from the highly complete DESI One-Percent survey exhibit a strong clustering signal on small scales, providing clear evidence for the conformity effect of ELGs. Building upon the original subhalo abundance matching (SHAM) method developed by Gao et al. (2022, 2023), we propose a concise conformity model to improve the ELG-halo connection. In this model, the number of satellite ELGs is boosted by a factor of $\sim 5$ in the halos whose central galaxies are ELGs. We show that the mean ELG satellite number in such central halos is still smaller than 1, and the model does not significantly increase the overall satellite fraction. With this model, we can well recover the ELG auto correlations to the smallest scales explored with the current data (i.e. $r_{\mathrm{p}} > 0.03$ $\mathrm{Mpc}\,h^{-1}$ in real space and at $s > 0.3$ $\mathrm{Mpc}\,h^{-1}$ in redshift space), while the cross correlations between luminous red galaxies (LRGs) and ELGs are nearly unchanged. Although our SHAM model has only 8 parameters, we further verify that it can accurately describe the ELG clustering in the entire redshift range from $z = 0.8$ to $1.6$. We therefore expect that this method can be used to generate high-quality ELG lightcone mocks for DESI.

Planet formation is directly linked to the birthing environment that protoplanetary disks provide. The disk properties determine whether a giant planet will form and how it evolves. The number of exoplanet and disk observations is consistently rising, however, it is not yet possible to directly link these two populations. Therefore, a deep theoretical understanding of how planets form is crucial. We performed numerical simulations of planet formation via pebble and gas accretion, while including migration, in a viscously evolving protoplanetary disk, with dust growing, drifting, and evaporating at the ice lines. In our investigation of the most favorable conditions for giant planet formation, we find that these are high disk masses, early formation, and a large enough disk to host a long-lasting pebble flux. However, small disks with the same mass allow more efficient gas accretion onto planetary cores, leading to more massive gas giants. Given the right conditions, high viscosity leads to more massive cores and it enhances gas accretion. It also causes faster type II migration rates, so the giants have a decreasing final position for increasing viscosity. Intermediate dust fragmentation velocities provide the necessary pebble sizes and radial drift velocities for maximized pebble accretion and pebble flux. An enhanced dust-to-gas ratio can compensate for lower disk masses, but early formation is still crucial. We conclude that there is no specific initial parameter that leads to giant planet formation; rather, it is the outcome of a combination of complementary factors. This also implies that the diversity of the exoplanet systems is the product of the intrinsic diversity of the protoplanetary disks and it is crucial to take advantage of the increasing number and quality of observations to constrain the disk population properties and ultimately devise planet formation theories.

We report the observations of single pulse emission from the pulsar B0809+74 at 150 MHz using the Polish LOFAR station, PL-611. The three major phenomena of subpulse drifting, nulling and mode changing associated with single pulse variations are prominently seen in these observations. The pulsar has a single component conal profile and the single pulses are primarily in the ``normal'' drift mode with periodicity ($P_3$) 11.1$\pm$0.5 $P$ for 96\% of the observing duration, while the shorter duration ``slow-drift'' mode has $P_3$ = 15.7$\pm$1.2 $P$. We were able to measure the phase behaviour associated with drifting from the fluctuation spectral analysis that showed identical linear phase variations across the pulse window for both modes despite their different periodic behaviour. Earlier studies reported that the transitions from the normal state to the slow-drift mode were preceded by the presence of nulling with typical durations of 5 to 10 periods. Our observations however seem to suggest that the transition to nulling follows shortly after the pulsar switches to the slow-drift mode and not at the boundary between the modes, with one instance of complete absence of nulling between mode switching. In addition we also detected a second type of short duration nulls not associated with the mode changing that showed quasi-periodic behaviour with periodicity, $P_N\sim44\pm7$. The variety of features revealed in the single pulse sequence makes PSR B0809+74 an ideal candidate to understand the physical processes in the Partially Screened Gap dominated by non-dipolar magnetic fields.

Intrinsic alignment (IA) of source galaxies represents an important contaminant for upcoming cosmic shear surveys. In particular, it is expected on general grounds that IA contains a B-mode while the weak lensing signal does not. Thus, a detection of B-modes offers the possibility to study directly the IA signal of the sources. Galaxy clusters exhibit strong IA and are therefore a natural candidate to look for a B-mode signal. We forecast the signal-to-noise ratio (SNR) for B-modes from IA of galaxy clusters in the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We use a perturbative model for the IA multipoles based on the Effective Field Theory of Intrinsic Alignments (EFT of IA), which has recently been validated against N-body simulations. For LSST, we forecast SNR $\approx 12$. We find that this detectability is not significantly impacted by different analysis choices. Lastly, we also apply our forecast to clusters in the SDSS redMaPPer and DESY1 samples, finding SNR $\approx 5$ and SNR $\approx 3$, respectively. This implies that a detection may already be within reach of current cluster samples.

Are Saturn's regular satellites young or old? And how old are Enceladus' cratered plains? To answer these questions we computed model surface ages of the most heavily cratered terrains on Saturn's regular icy satellites using new high-resolution outer Solar System evolution simulations, and coupled with improved estimates of the trans-Neptunian objects populations. The output of the simulations allowed us to construct a model impact chronology onto Saturn which automatically applies to its regular satellites. We used crater densities and our impact chronology onto Saturn to construct model impact-crater isochrons, i.e., the scaling of the satellite crater production function with time. The surface ages derived for the cratered plains on Mimas, Enceladus, Tethys, Dione and Rhea range from 4.1 Ga to 4.4 Ga, with the surfaces of Mimas and Enceladus roughly 200 Myr younger than those of the outer three satellites. Uncertainties in these ages are less than 300 Myr. The calculated model surface ages of these satellites are consistent over as much as two orders of magnitude in the observed crater diameter. The similarity of the crater production function amongst all satellites suggests that they were bombarded by a single impactor source. This work supports the idea that Saturn's regular satellites are ancient, and has implications for their formation and their tidal evolution.

Some solar-type stars are known to present faint, time-variable radio continuum emission whose nature is not clearly established. We report on Jansky Very Large Array observations of the nearby star $\epsilon$ Eridani at 10.0 and 33.0 GHz. We find that this star has flux density variations on scales down to days, hours and minutes. On 2020 Apr 15 it exhibited a radio pulse at 10.0 GHz with a total duration of about 20 minutes and a peak four times larger than the plateau of 40 $\mu$Jy present in that epoch. We were able to model the time behavior of this radio pulse in terms of the radiation from shocks ramming into the stellar wind. Such shocks can be produced by the wind interaction of violently expanding gas heated suddenly by energetic electrons from a stellar flare, similar to the observed solar flares. Because of the large temperature needed in the working surface to produce the observed emission, this has to be non thermal. It could be gyrosynchrotron or synchrotron emission. Unfortunately, the spectral index or polarization measurements from the radio pulse do not have enough signal-to-noise ratio to determine its nature.

Although it is well known that the bulk of dark matter (DM) has to be cold, the existence of an additional sub-dominant, hot species remains a valid possibility. In this paper we investigate the potential of the cosmic shear power spectrum to constrain such a mixed (hot plus cold) DM scenario with two additional free parameters, the hot-to-total DM fraction ($f_{\rm hdm}$) and the thermal mass of the hot component ($m_{\rm hdm}$). Running a Bayesian inference analysis for both the Kilo-Degree Survey cosmic shear data (KiDS) as well as the Cosmic Microwave Background (CMB) temperature and polarisation data from Planck, we derive new constraints for the mixed DM scenario. We find a 95 per cent confidence limit of $f_{\rm hdm}<0.08$ for a very hot species of $m_{\rm hdm}\leq20$ eV. This constraint is weakened to $f_{\rm hdm}<0.25$ for $m_{\rm hdm}\leq80$ eV. Scenarios with masses above $m_{\rm hdm}\sim200$ eV remain unconstrained by the data. Next to providing limits, we investigate the potential of mixed DM to address the clustering (or $S_8$) tension between lensing and the CMB. We find a reduction of the 2D ($\Omega_m - S_8$) tension from 2.9$\sigma$ to 1.6$\sigma$ when going from a pure cold DM to a mixed DM scenario. When computing the 1D gaussian tension on $S_8$ the improvement is milder, from 2.4$\sigma$ to 2.0$\sigma$.

To test theoretical models of massive star formation it is important to compare their predictions with observed systems. To this end, we conduct CO molecular line radiative transfer post-processing of 3D magneto-hydrodynamic (MHD) simulations of various stages in the evolutionary sequence of a massive protostellar core, including its infall envelope and disk wind outflow. Synthetic position-position-velocity (PPV) cubes of various transitions of CO, 13CO, and C18O emission are generated. We also carry out simulated Atacama Large Millimeter/submillimeter Array (ALMA) observations of this emission. We compare the mass, momentum and kinetic energy estimates obtained from molecular lines to the true values, finding that the mass and momentum estimates can have uncertainties of up to a factor of four. However, the kinetic energy estimated from molecular lines is more significantly underestimated. Additionally, we compare the mass outflow rate and momentum outflow rate obtained from the synthetic spectra with the true values. Finally, we compare the synthetic spectra with real examples of ALMA-observed protostars and determine the best fitting protostellar masses and outflow inclination angles. We then calculate the mass outflow rate and momentum outflow rate for these sources, finding that both rates closely match theoretical protostellar evolutionary tracks.

Although the present-day orbital distribution of minor bodies that go around the Sun between the orbit of Neptune and the Kuiper Cliff is well understood, past ~50 au from the Sun, our vision gets blurred as objects become fainter and fainter and their orbital periods span several centuries. Deep imaging using the largest telescopes can overcome the first issue but the problems derived from the second one are better addressed using data analysis techniques. Here, we make use of the heliocentric range and range-rate of the known Kuiper belt objects and their uncertainties to identify structures in orbital parameter space beyond the Kuiper Cliff. The distribution in heliocentric range there closely resembles that of the outer main asteroid belt with a gap at 70 au that may signal the existence of a dynamical analogue of the Jupiter family comets. Outliers in the distribution of mutual nodal distances suggest that a massive perturber is present beyond the heliopause.

Massive protostars launch accretion-powered, magnetically-collimated outflows, which play crucial roles in the dynamics and diagnostics of the star formation process. Here we calculate the shock heating and resulting free-free radio emission in numerical models of outflows of massive star formation within the framework of the Turbulent Core Accretion model. We post-process 3D magneto-hydrodynamic simulation snapshots of a protostellar disk wind interacts with an infalling core envelope, and calculate shock temperatures, ionization fractions, and radio free-free emission. We find heating up to ~10^7 K and near complete ionization in shocks at the interface between the outflow cavity and infalling envelope. However, line-of-sight averaged ionization fractions peak around ~10%, in agreement with values reported from observations of massive protostar G35.20-0.74N. By calculating radio continuum fluxes and spectra, we compare our models with observed samples of massive protostars. We find our fiducial models produce radio luminosities similar to those seen from low and intermediate-mass protostars that are thought to be powered by shock ionization. Comparing to more massive protostars, we find our model radio luminosities are ~10 to 100 times less luminous. We discuss how this apparent discrepancy either reflects aspects of our modeling related to the treatment of cooling of the post-shock gas or a dominant contribution in the observed systems from photoionization. Finally, our models exhibit 10-year radio flux variability of ~5%, especially in the inner 1000 au region, comparable to observed levels in some hyper-compact HII regions.

The Cryogenic IR echelle Spectrometer (CRIRES) instrument at the Very Large Telescope (VLT) was in operation from 2006 to 2014. Great strides in characterizing the inner regions of protoplanetary disks were made using CRIRES observations in the L- and M-band at this time. The upgraded instrument, CRIRES+, became available in 2021 and covers a larger wavelength range simultaneously. Here we present new CRIRES+ Science Verification data of the binary system S Coronae Australis (S CrA). We aim to characterize the upgraded CRIRES+ instrument for disk studies and provide new insight into the gas in the inner disk of the S CrA N and S systems. We analyze the CRIRES+ data taken in all available L- and M-band settings, providing spectral coverage from 2.9 to 5.5 $\mu$m. We detect emission from $^{12}$CO (v=1-0, v=2-1, and v=3-2), $^{13}$CO (v=1-0), hydrogen recombination lines, OH, and H$_2$O in the S CrA N disk. In the fainter S CrA S system, only the $^{12}$CO v=1-0 and the hydrogen recombination lines are detected. The $^{12}$CO v=1-0 emission in S CrA N and S shows two velocity components, a broad component coming from $\sim$0.1 au in S CrA N and $\sim$0.03 au in S CrA S and a narrow component coming from $\sim$3 au in S CrA N and $\sim$5 au in S CrA S. We fit local thermodynamic equilibrium slab models to the rotation diagrams of the two S CrA N velocity components and find that they have similar column densities ($\sim$1-7$\times$10$^{17}$ cm$^{-2}$), but that the broad component is coming from a hotter and narrower region. Two filter settings, M4211 and M4368, provide sufficient wavelength coverage for characterizing CO and H$_2$O at $\sim$5 $\mu$m, in particular covering low- and high-$J$ lines. CRIRES+ provides spectral coverage and resolution that are crucial complements to low-resolution observations, such as those with JWST, where multiple velocity components cannot be distinguished.

We study the effect of asymmetric fermionic dark matter (DM) on the thermal evolution of neutron stars (NSs). No interaction between DM and baryonic matter is assumed, except the gravitational one. Using the two-fluid formalism, we show that DM accumulated in the core of a star pulls inwards the outer baryonic layers of the star, increasing the baryonic density in the NS core. As a result, it significantly affects the star's thermal evolution by triggering an early onset of the direct Urca process and modifying the photon emission from the surface caused by the decrease of the radius. Thus, due to the gravitational pull of DM, the direct Urca process becomes kinematically allowed for stars with lower masses. Based on these results, we discuss the importance of NS observations at different distances from the Galactic center. Since the DM distribution peaks towards the Galactic center, NSs in this region are expected to contain higher DM fractions that could lead to a different cooling behavior.

Many cosmological models assume or imply that the total size of the universe is very large, perhaps even infinite. Here we argue instead that the universe might be comparatively small, in fact not much larger than the currently observed size. A concrete implementation of this idea is provided by the no-boundary proposal, in combination with a plateau-shaped inflationary potential. In this model, opposing effects of the weighting of the wave function and of the criterion of allowability of the geometries conspire to favour small universes. We point out that a small size of the universe also fits well with swampland conjectures, and we comment on the relation with the dark dimension scenario.

Axions and axion-like particles are well-motivated dark matter candidates. We propose an experiment that uses single photon detection interferometry to search for axions and axion-like particles in the galactic halo. We show that photon counting with a dark rate of 6E-6 Hz can improve the quantum sensitivity of axion interferometry by a factor of 15 compared to the quantum-enhanced heterodyne readout for 5-m long optical resonators. The proposed experimental method has the potential to be scaled up to kilometer-long facilities, enabling the detection or setting of constraints on the axion-photon coupling coefficient of 1E-17 - 1E-16 GeV-1 for axion masses ranging from 0.1 to 1 neV.

Modeling environments that are not in local thermal equilibrium, such as protoplanetary disks or planetary atmospheres, with molecular spectroscopic data from space telescopes requires knowledge of the rate coefficients of rovibrationally inelastic molecular collisions. Here, we present such rate coefficients in a temperature range from 10 to 500 K for collisions of CO$_2$ with He atoms in which CO$_2$ is (de)excited in the bend mode. They are obtained from numerically exact coupled-channel (CC) calculations as well as from calculations with the less demanding coupled-states approximation (CSA) and the vibrational close-coupling rotational infinite-order sudden (VCC-IOS) method. All of the calculations are based on a newly calculated accurate \textit{ab initio} four-dimensional CO$_2$-He potential surface including the CO$_2$ bend ($\nu_2$) mode. We find that the rovibrationally inelastic collision cross sections and rate coefficients from the CSA and VCC-IOS calculations agree to within 50% with the CC results at the rotational state-to-state level, except for the smaller ones and in the low energy resonance region, and to within 20% for the overall vibrational quenching rates except for temperatures below 50 K where resonances provide a substantial contribution. Our CC quenching rates agree with the most recent experimental data within the error bars. We also compared our results with data from Clary et al. calculated in the 1980's with the CSA and VCC-IOS methods and a simple atom-atom model potential based on ab initio Hartree-Fock calculations and found that their cross sections agree fairly well with ours for collision energies above 500 cm$^{-1}$, but that the inclusion of long range attractive dispersion interactions is crucial to obtain reliable cross sections at lower energies and rate coefficients at lower temperatures.

We study the tidal problem and the resulting $I$-Love-$Q$ approximate universal relations for rotating superfluid neutron stars in the Hartle-Thorne formalism. Superfluid stars are described in this work by means of a two-fluid model consisting of superfluid neutrons and all other charged constituents. We employ a stationary and axisymmetric perturbation scheme to second order around a static and spherically symmetric background. Recently, we used this scheme to study isolated rotating superfluid stars. In this paper it is applied to analyze the axially symmetric sector of the tidal problem in a binary system. We show that a consistent use of perturbative matching theory amends the original two-fluid formalism for the tidal problem to account for the possible non-zero value of the energy density at the boundary of the star. This is exemplified by building numerically different stellar models spanning three equations of state. Significant departures from universality are found when the correct matching relations are not taken into account. We also present an augmented set of universal relations for superfluid neutron stars which includes the contribution to the total mass of the star at second order, $\delta M$. Therefore, our results complete the set of universal relations for rotating superfluid stars, generalizing our previous findings in the perfect fluid case.

The dark sirens method enables us to use gravitational wave events without electromagnetic counterparts as tools for cosmology and tests of gravity. Furthermore, the dark sirens analysis code gwcosmo can now robustly account for information coming from both galaxy catalogues and the compact object mass distribution. We present here an extension of the gwcosmo code and methodology to constrain parameterized deviations from General Relativity that affect the propagation of gravitational waves. We show results of our analysis using data from the GWTC-3 gravitational wave catalogues, in preparation for application to the O4 observing run. After testing our pipelines using the First Two Years mock data set, we reanalyse 46 events from GWTC-3, and combine the posterior for BBH and NSBH sampling results for the first time. We obtain joint constraints on H0 and parameterized deviations from General Relativity in the Power Law + Peak BBH population model. With increased galaxy catalogue support in the future, our work sets the stage for dark sirens to become a powerful tool for testing gravity.

The cores of dense stars are a powerful laboratory for studying feebly-coupled particles such as axions. Some of the strongest constraints on axionlike particles and their couplings to ordinary matter derive from considerations of stellar axion emission. In this work we study the radiation of axionlike particles from degenerate neutron star matter via a lepton-flavor-violating (LFV) coupling that leads to muon-electron conversion when an axion is emitted. We calculate the axion emission rate per unit volume (emissivity) and by comparing with the rate of neutrino emission, we infer upper limits on the LFV coupling that are at the level of $|g_{ae\mu}| \lesssim 10^{-6}$. For the hotter environment of a supernova, such as SN 1987A, the axion emission rate is enhanced and the limit is stronger, at the level of $|g_{ae\mu}| \lesssim 10^{-11}$, competitive with laboratory limits. Interestingly, our derivation of the axion emissivity reveals that axion emission via the LFV coupling is suppressed relative to the familiar lepton-flavor-preserving channels by a factor of $T^2 E_{F,e}^2 / (m_\mu^2 - m_e^2)^2 \sim T^2/m_\mu^2$, which is responsible for the relatively weaker limits.