The Union3 "Spline-Interpolated Distance Moduli" model posterior has been distributed for third-party cosmology analysis. The posterior prefers a large value of $\Omega_M$, a small absolute value of $w_0$, and a negative $w_a$, but still accommodates $\Lambda$CDM; the supernova data alone are not strongly constraining. The posterior is built assuming an underlying model and prior, both of which must be made to conform with any new model and prior being analyzed. The posterior is calculated for a prior that is not flat but rather has non-trivial structure in $\Omega_M$-$w_0$-$w_a$; the associated likelihood is slightly shifted relative to the posterior. The posterior for a prior that is flat in $\Omega_M$-$w_0$-$w_a$ is also shifted relative to the original, but not at a level that is statistically significant. The misapplication of Union3 results in the "DESI2024 VI" cosmology fits are inconsequential.

The evolution of the Hubble parameter $H(z)$ with redshift $z$ is estimated from the Pantheon+ data of Type Ia supernovae, for the $\Lambda$CDM model and the three special cases of the eternal coasting (EC) cosmological model with three different spatial geometries. The scatter associated with $H(z)$ is seen to grow markedly with redshift. This behaviour, which is deduced directly from the SNe Hubble diagram, raises the question of whether the universe is undergoing a stochastic expansion, which scenario can offer an explanation for the Hubble tension in cosmology. From the estimated $H(z)$ values, the present value of the Hubble parameter $H_0$ is evaluated for each of these models through regression, and the scatter using the Monte Carlo method. Bayesian comparison between these models is carried out using the data of 35 cosmic chronometers (CC). The comparative study favours the $\Lambda$CDM model, with some strong evidence. However, exclusion of four outlier CC data points with small errorbars leads to large reduction in the Bayes factor value. The unusually large value of Bayes factor obtained while using the full set of CC data raises some concerns about its tension with other data, such as that of the SNe Ia. While using the remaining 31 CC data points, it is observed that the resulting Bayes factor still favours the $\Lambda$CDM model, but with a much smaller value of the Bayes factor. When EC models are compared among themselves, the $\Omega = 2$ model has strong evidence than the $\Omega = 1$ (also known as $R_h = ct$) and the $\Omega = 0$ (Milne-type) models.

Recent observations by the Large High Altitude Air Shower Observatory (LHAASO) detected Ultra High Energy (UHE) photons in the range 100 TeV to 1.4 PeV from twelve sources including Crab nebula. The detection of these photons demands the presence of at least PeV energy particle in the source. It is important to understand particle acceleration and radiation emission processes in such source. One of those twelve sources, LHAASO J2108+5157 does not show any association or counterparts at any other wavelength. In search of counterpart, we surveyed the region with Giant Metrewave Radio Telescope (GMRT) at 650 MHz frequency. GMRT observation revel radio emission from an extended source within the PSF of LHAASO which shows disk-jet morphology. Considering the spatial association and extent of the source, it is plausible that particle acceleration to PeV energies originates from this source.

The Boltzmann equation relates the equilibrium phase space distribution of stars in the Milky Way to the Galaxy's gravitational potential. However, observations of stellar populations are biased by extinction from foreground dust, which complicates measurements of the potential in the disk and towards the Galactic center. Using the kinematics of Red Clump and Red Branch stars in Gaia DR3, we use machine learning to simultaneously estimate both the unbiased stellar phase space density and the gravitational potential. The unbiased phase space density is obtained through a learned "dust efficiency factor" -- an observational selection function that accounts for dust extinction. The potential and the dust efficiency are parameterized by fully connected neural networks and are completely data driven. We validate the dust efficiency using a recent three-dimensional dust map in this work, and examine the potential in a companion paper.

Gravitational Wave (GW) sources are standard sirens that provide an independent way to map the cosmic expansion history by combining with an independent redshift measurement either from an electromagnetic counterpart for a bright siren or using different statistical techniques for dark sirens. In this analysis, we perform the first Blinded Mock Data Challenge (Blinded-MDC) to test the robustness in inferring the value of Hubble constant $H_0$ for a dark siren technique which depends on astrophysical mass distribution of Binary Black Holes (BBHs). We have considered different analysis setups for the Blinded-MDC to test both statistical and systematic uncertainties and demonstrate the capabilities in inferring $H_0$ with detector sensitivity as per the fourth observation run of LIGO-Virgo-KAGRA. We find that when the astrophysical population of BBHs matches with the underlying assumption of the model, a cosmological pipeline can recover the injected parameters using the observed mass distribution. However, when the mock mass distribution of the astrophysical population depends slightly on redshift and one is ignorant about it in analyzing the data, it can cause a systematic discrepancy in the inferred value of $H_0$ by about $1.5\sigma$, above the statistical fluctuations due to noise and a limited number of events. The discrepancy in $H_0$ is arising due to astrophysical mis-modeling, and in the future, elaborate studies will be required to mitigate systematic uncertainties due to unknown astrophysical populations of BBHs. This MDC framework sets the road map for inspecting the precision and accuracy of standard siren cosmology and provides the first insight into the robustness of the population-dependent cosmology inference in a blinded analysis setup.

``Little Red Dots'' (LRDs) are an abundant high-redshift population newly discovered by the James Webb Space Telescope (JWST). They are characterized by a red color in the rest-frame optical band, compact morphology, and broad Balmer emission lines (${\rm FWHM} \gtrsim 1000~{\rm km\,s^{-1}}$) that suggest an AGN nature. Using a method of pixel-by-pixel color selection and relaxing the compactness criteria, we identify three of the first dual LRD candidates in the COSMOS-Web survey with projected separations of $0.\!\!^{\prime\prime}2-0.\!\!^{\prime\prime}4$ (1-2 pkpc at their photometric redshifts). A comparison between existing LRD samples and mock data reveals that the projected separations of these dual LRD candidates are unlikely to result from chance projections of objects at different redshifts. In one case (CW-B5-15958), the dual LRD includes two bright sources ($m_{\rm F444W}=24.3$ and $24.8$) with characteristic V-shape spectral energy distribution (SEDs) and photometric redshifts consistent with each other. We find that CW-B5-15958 has a faint off-centered component and a companion galaxy. In the other two dual systems, the brighter LRD exhibits a V-shape SED, while the fainter LRD ($m_{\rm F444W}\gtrsim26$) is undetected in both F115W and F150W. These discoveries suggest that the angular auto-correlation function (ACF) of LRDs exhibits a significant excess ($\sim3\times10^2$ times) on sub-arcsec (kilo-parsec) separations compared to the extrapolation of a power-law ACF of JWST-found AGNs measured over $10^{\prime\prime}-100^{\prime\prime}$. Follow-up spectroscopic confirmation of their redshifts and the construction of a larger sample are essential to advance our understanding of the evolution of supermassive black holes and the importance of mergers in the early universe.

Early JWST observations are providing growing evidence for a ubiquitous population of accreting supermassive black holes (BHs) at high redshift, many of which appear overmassive compared to the empirically-derived local scaling relation between black hole mass and host galaxy stellar mass. In this study, we leverage predictions from the semi-analytical Cosmic Archaeology Tool (CAT) to reconstruct the evolutionary pathways for this overmassive BH population, investigating how they assemble over cosmic time and interact with their host galaxies. We find that the large $M_{\rm BH}-M_{\rm star}$ ratios can be explained if light and heavy BH seeds grow by short, repeated episodes of super-Eddington accretion, triggered by major galaxy mergers. On average, we find that BH-galaxy co-evolution starts in earnest only at $z < 8$, when $\simeq 30\%$ of the final galaxy stellar mass has formed outside the massive black hole host. Our model suggests that super-Eddington bursts of accretion last between $0.5-3$ Myr, resulting in a duty cycle of $1-4 \%$ for the target BH sample. The boost in luminosity of BHs undergoing super-Eddington accretion helps explaining the luminosity function of Active Galactic Nuclei observed by JWST. At the same time, a large population of these overmassive BHs are predicted to be inactive, with Eddington ratio $\lambda_{\rm Edd} < 0.05$, in agreement with recent observations.

In this study, we explore the characteristics of `low-mass' ($\log(M_{\star}/M_{\odot}) \leq 8$) and `intermediate-mass' ($8 \lt \log(M_{\star}/M_{\odot}) \leq 10$) host galaxies of Type Ia supernovae (SNe Ia) from the second data release (DR2) of the Zwicky Transient Facility survey and investigate their correlations with different sub-types of SNe Ia. We use the photospheric velocities measured from the Si II $\lambda$6355 feature, SALT2 light-curve stretch ($x_1$) and host-galaxy properties of SNe Ia to re-investigate the existing relationship between host galaxy mass and Si II $\lambda$6355 velocities. We also investigate sub-type preferences for host populations and show that while the more energetic and brighter 91T-like SNe Ia tends to populate the younger host populations, 91bg-like SNe Ia populate in the older populations. Our findings suggest High Velocity SNe Ia (HV SNe Ia) not only comes from the older populations but they also come from young populations as well. Therefore, while our findings can partially provide support for HV SNe Ia relating to single degenerate progenitor models, they indicate that HV SNe Ia other than being a different population, might be a continued distribution with different explosion mechanisms. We lastly investigate the specific rate of SNe Ia in the volume-limited SN Ia sample of DR2 and compare with other surveys.

We propose a novel scenario in which scalar perturbations, that seed the large scale structure of the Universe, are generated without relying on a scalar field (the inflaton). In this framework, inflation is driven by a de Sitter space-time (dS), where tensor metric fluctuations (i.e., gravitational waves) naturally arise from quantum vacuum oscillations, and scalar fluctuations are generated via second-order tensor effects. We compute the power spectrum of such scalar fluctuations and show it to be consistent with near scale-invariance. We derive the necessary conditions under which scalar perturbations become significant and much larger than the tensor modes, and we identify a natural mechanism to end inflation via a transition to a radiation-dominated phase. Our proposed mechanism could remove the need for a model-dependent scenario: the choice of a scalar field, as the inflaton, to drive inflation.

We used the LOFAR telescope to monitor SN 2023ixf, a core-collapse supernova in M101, between 8 and 368 days post-explosion. We report non-detections down to ~80 {\mu}Jy sensitivity at 144 MHz. Our non-detections are consistent with published radio observations at higher frequencies. At the time, we are not able to constrain the properties of low-frequency absorption due to the progenitor star's circumstellar medium via these LOFAR observations.

In this work, I investigate the impact of Dark Energy Spectroscopic Instrument (DESI) Baryonic Acoustic Oscillations (BAO) data on cosmological parameters, focusing on the inflationary spectral index $n_s$, the amplitude of scalar perturbations $A_s$, and the matter density parameter $\omega_m$. By examining different models of late-time new physics, the inflationary parameters were revealed to be stable when compared with the baseline dataset that used the earlier BAO data from the SDSS collaboration. When combined with Cosmic Microwave Background (CMB) and type Ia supernovae (SNeIa), DESI BAO data leads to a slight reduction in $\omega_m$ (less than 2\%) and modest changes in $A_s$ and $n_s$, if compared with the same combination but using SDSS BAO data instead, suggesting a subtle shift in matter clustering. These effects may be attributed to a higher expansion rate from dynamical dark energy, changes in the recombination period, or modifications to the matter-radiation equality time. Further analyses of models with dynamical dark energy and free curvature show a consistent trend of reduced $\omega_m$, accompanied by slight increases in both $n_s$ and $H_0$. The results emphasize the importance of the DESI BAO data in refining cosmological parameter estimates and highlight the stability of inflationary parameters across different late-time cosmological models.

We used random forest algorithms to classify all objects in a large portion of the sky, using optical light curves obtained, or built from images provided, by the Zwicky Transient Facility (ZTF). We compare different selection sets based on alerts or complete light curves derived from different photometric selection algorithms. The AGN candidates thus selected are cross-matched with objects in the NASA-Sloan Atlas (NSA) of local galaxies, with $M_*<2\times10^{10}M_\odot$. The AGN nature of these candidates is verified and characterized using archival optical spectra from SDSS. We further establish the fraction of candidates with counterparts in the eROSITA data release 1 catalog of X-ray sources. From an initial sample of 506 candidates, 415 have good-quality spectra. Among these 415 objects, we found significant broad Balmer lines in the spectra for $86\%$ (357) of the candidates. When considering BPT classifications, an additional 5 candidates were confirmed, resulting in $87\%$ (362) confirmed candidates. Specifically, broad Balmer lines were detected in $94\%$-$98\%$ of the AGN candidates selected from complete light curves and in $80\%$ of those selected from the less frequent ZTF alerts. The black hole masses estimated from the spectra range from $2.2\times10^6M_\odot$ to $4.2\times10^7M_\odot$, reaching lower values for the candidates selected using the more sensitive light curves. The black hole masses obtained cluster around $0.1\%$ of the stellar mass of the host from the NSA catalog. Two-thirds of the AGN candidates are classified as Seyfert or Composite by their narrow emission line ratios (BPT diagnostics) while the rest are star-forming. Almost all the candidates classified as Seyfert and over $50\%$ of those classified as star-forming have significant BELs. We found X-ray counterparts for $67\%$ of the candidates that fall in the footprint of the eROSITA-DE DR1.

In models of warm dark matter, there is an appreciable population of high momentum particles in the early universe, which free stream out of primordial over/under densities, thereby prohibiting the growth of structure on small length scales. The distance that a dark matter particle travels without obstruction, known as the free streaming length, depends on the particle's mass and momentum, but also on the cosmological expansion rate. In this way, measurements of the linear matter power spectrum serve to probe warm dark matter as well as the cosmological expansion history. In this work, we focus on ultra-light wave wave dark matter (WWDM) characterized by a typical comoving momentum $q_\ast$ and mass $m$. We first derive constraints on the WWDM parameter space $(q_\ast, m)$ using Lyman-$\alpha$ forest observations due to a combination of the free-streaming effect and the white-noise effect. We next assess how the free streaming of WWDM is affected by three modified expansion histories: early matter domination, early dark energy, and very early dark energy.

The existence of planes of satellite galaxies has been identified as a long-standing challenge to $\Lambda$CDM cosmology, due to the rarity of satellite systems in cosmological simulations that are as extremely flattened and as strongly kinematically correlated as observed structures. Here we investigate a recently proposed new metric to measure the overall degree of ''planarity'' of a satellite system, which was used to claim consistency between the Milky Way satellite plane and $\Lambda$CDM. We study the behavior of the ''planarity'' metric under several features of anisotropy present in $\Lambda$CDM satellite systems but unrelated to satellite planes. Specifically, we consider the impact of oblate or prolate distributions, the number of satellites, clustering of satellites, and radial and asymmetric distributions ('lopsidedness'). We also investigate whether the metric is independent of the orientation of the studied satellite system. We find that all of these features of anisotropy result in the metric inferring an increased degree of ''planarity'', despite none of them having any direct relation to satellite planes. The metric is also highly sensitive to the orientation of the studied system (or chosen coordinate system): there is almost no correlation between the metric's reported degrees of ''planarity'' for identical random systems rotated by 90{\deg}. Our results demonstrate that the new proposed metric is unsuitable to measure overall ''planarity'' in satellite systems. Consequently, no consistency of the observed Milky Way satellite plane with $\Lambda$CDM can be inferred using this metric.

Models predict that the well-studied mass-radius relation of white dwarf stars depends on the temperature of the star, with hotter white dwarfs having larger masses at a given radius than cooler stars. In this paper, we use a catalog of 26,041 DA white dwarfs observed in Sloan Digital Sky Survey Data Releases 1-19. We measure the radial velocity, effective temperature, surface gravity, and radius for each object. By binning this catalog in radius or surface gravity, we average out the random motion component of the radial velocities for nearby white dwarfs to isolate the gravitational redshifts for these objects and use them to directly measure the mass-radius relation. For gravitational redshifts measured from binning in either radius or surface gravity, we find strong evidence for a temperature-dependent mass-radius relation, with warmer white dwarfs consistently having greater gravitational redshifts than cool objects at a fixed radius or surface gravity. For warm white dwarfs, we find that their mean radius is larger and mean surface gravity is smaller than those of cool white dwarfs at 5.2{\sigma} and 6.0{\sigma} significance, respectively. Selecting white dwarfs with similar radii or surface gravities, the significance of the difference in mean gravitational redshifts between the warm and cool samples is >6.1{\sigma} and >3.6{\sigma} for measurements binned in radius and surface gravity, respectively, in the direction predicted by theory. This is an improvement over previous implicit detections, and our technique can be expanded to precisely test the white dwarf mass-radius relation with future surveys.

We present Swift Ultraviolet Optical Telescope (UVOT) observations of the deep field GOODS-N in four near-UV filters. A catalog of detected galaxies is reported, which will be used to explore galaxy evolution using ultraviolet emission. Swift/UVOT observations probe galaxies at $z \lesssim 1.5$ and combine a wide field of view with moderate spatial resolution; these data complement the wide-field observations of GALEX and the deep, high angular resolution observations by HST. Using our catalog of detected galaxies, we calculate the UV galaxy number counts as a function of apparent magnitude and compute the UV luminosity function and its evolution with redshift. From the luminosity function fits in various redshift bins, we calculate the star formation rate density as a function of redshift and find evolution consistent with past works. We explore how different assumptions such as dust attenuation corrections can dramatically change how quickly the corrected star formation rate density changes with redshift. At these low redshifts, we find no trend between UV attenuation and redshift or absolute magnitude with significant scatter in the UV spectral slope $\beta$. This dataset will complement the extensive observations of GOODS-N already in the literature.

We present the initial results of MEAD (Measuring Extinction and Abundances of Dust), with a focus on the dust extinction features observed in our JWST near- and mid-infrared spectra of nine diffuse Milky Way sightlines ($1.2 \leq A(V) \leq 2.5$). For the first time, we find strong correlations between the 10 $\mu$m silicate feature strength and the column densities of Mg, Fe and O in dust. This is consistent with the well-established theory that Mg- and Fe-rich silicates are responsible for this feature. We obtained an average stoichiometry of the silicate grains in our sample of Mg:Fe:O = 1.1:1:11.2, constraining the grain composition. We find variations in the feature properties, indicating that different sightlines contain different types of silicates. In the average spectrum of our sample, we tentatively detect features around 3.4 and 6.2 $\mu$m, which are likely caused by aliphatic and aromatic/olefinic hydrocarbons, respectively. If real, to our knowledge, this is the first detection of hydrocarbons in purely diffuse sightlines with $A(V) \leq 2.5$, confirming the presence of these grains in diffuse environments. We detected a 3 $\mu$m feature toward HD073882, and tentatively in the sample average, likely caused by water ice (or solid-state water trapped on silicate grains). If confirmed, to our knowledge, this is the first detection of ice in purely diffuse sightlines with $A(V) \leq 2.5$, supporting previous findings that these molecules can exist in the diffuse ISM.

UHECR are evaluated in the frame role of different nuclei composition . Most of the past and present models are considering proton or iron as their main courier. Some attention has been paid to the role of the UHECR light nuclei in recent years. We update here the lightest nuclei UHECR model, able to explain the nearest AGN or Star Burst sources with the few observed Hot Spot clustering in AUGER and TA array data. Any additional components of the heaviest nuclei with the highest energy, more bent and smeared, may also fit recent AUGER and TA homogeneous records at those energy edges.

In the early stages of planetary system formation, young exoplanets gravitationally interact with their surrounding environments and leave observable signatures on protoplanetary disks. Among these structures, a pair of nearly symmetric spiral arms can be driven by a giant protoplanet. For the double-spiraled SAO 206462 protoplanetary disk, we obtained three epochs of observations spanning 7 yr using the Very Large Telescope's SPHERE instrument in near-infrared $J$-band polarized light. By jointly measuring the motion of the two spirals at three epochs, we obtained a rotation rate of $-0.85^\circ\pm0.05^\circ~{\rm yr}^{-1}$. This rate corresponds to a protoplanet at $66\pm3$ au on a circular orbit dynamically driving both spirals. The derived location agrees with the gap in ALMA dust-continuum observations, indicating that the spiral driver may also carve the observed gap. What is more, a dust filament at $\sim$63 au observed by ALMA coincides with the predicted orbit of the spiral-arm-driving protoplanet. This double-spiraled system is an ideal target for protoplanet imaging.

We present high resolution WIYN/NEID echelle spectroscopy (R $\approx70$,000) of the supernova (SN) 2023ixf in M101, obtained 1.51 to 18.51 days after explosion over nine epochs. Daily monitoring for the first four days after explosion shows narrow emission features ($\leq200$ km s$^{-1}$), exhibiting predominantly blueshifted velocities, that rapidly weaken, broaden, and vanish in a manner consistent with radiative acceleration and the SN shock eventually overrunning or enveloping the full extent of dense circumstellar medium (CSM). The most rapid evolution is in the He I emission, which is visible on day 1.51 but disappears by day 2.62. We measure the maximum pre-SN speed of He I to be 25 $^{+0}_{-5} \pm2$ km s$^{-1}$, where the error is attributable to the uncertainty in how much the He I had already been radiatively accelerated, and to measurement of the emission line profile. The radiative acceleration of material is likely driven by the shock-CSM interaction, and the CSM is accelerated to $\geq200$ km s$^{-1}$ before being completely swept up by the SN shock to $\sim 2000$ km s$^{-1}$. We compare the observed spectra with spherically-symmetric r16wb HERACLES/CMFGEN model spectra and find the line evolution to generally be consistent with radiative acceleration and optical depth effects. The progenitor of SN2023ixf underwent an enhanced mass loss phase $\gtrsim 4$ year prior to core-collapse, creating a dense, asymmetric CSM region extending out to approximately $r_{CSM} = 3.7 \times 10^{14}$ ($v_\textrm{shock}$/9500 km s$^{-1}$) cm.

We investigate the evolution of active galactic nucleus jets on kiloparsec-scales due to their interaction with the clumpy interstellar medium (ISM) of the host galaxy and, subsequently, the surrounding circumgalactic environment. Hydrodynamic simulations of this jet-environment interaction are presented for a range of jet kinetic powers, peak densities of the multiphase ISM, and scale radii of the larger-scale environment -- characteristic of either a galaxy cluster or poor group. Synthetic radio images are generated by considering the combination of synchrotron radiation from the jet plasma and free-free absorption from the multiphase ISM. We find that jet propagation is slowed by interactions with a few very dense clouds in the host galaxy ISM, producing asymmetries in lobe length and brightness which persist to scales of tens of kpc for poor group environments. The classification of kiloparsec-scale jets is highly dependent on surface brightness sensitivity and resolution. Our simulations of young active sources can appear as restarted sources, showing double-double lobe morphology, high core prominence (CP > 0.1), and the expected radio spectra for both the inner- and outer-lobe components. We qualitatively reproduce the observed inverse correlation between peak frequency and source size, and find that the peak frequency of the integrated radio spectrum depends on ISM density but not the jet power. Spectral turnover in resolved young radio sources therefore provides a new probe of the ISM.

With the low Earth orbit environment becoming increasingly populated with artificial satellites, rockets, and debris, it is important to understand the effects they have on radio astronomy. In this work, we undertake a multi-frequency, multi-epoch survey with two SKA-Low station prototypes located at the SKA-Low site, to identify and characterise radio frequency emission from orbiting objects and consider their impact on radio astronomy observations. We identified 152 unique satellites across multiple passes in low and medium Earth orbits from 1.6 million full-sky images across 13 selected ${\approx}1$ MHz frequency bands in the SKA-Low frequency range, acquired over almost 20 days of data collection. Our algorithms significantly reduce the rate of satellite misidentification, compared to previous work, validated through simulations to be $<1\%$. Notably, multiple satellites were detected transmitting unintended electromagnetic radiation, as well as several decommissioned satellites likely transmitting when the Sun illuminates their solar panels. We test alternative methods of processing data, which will be deployed for a larger, more systematic survey at SKA-Low frequencies in the near future. The current work establishes a baseline for monitoring satellite transmissions, which will be repeated in future years to assess their evolving impact on radio astronomy observations.

Radio galaxies can extend far beyond the stellar component of their originating host galaxies, and their radio emission can consist of multiple discrete components. Furthermore, the apparent source structure will depend on survey sensitivity, resolution and the observing frequency. Associated discrete radio components and their originating host galaxy are typically identified through a visual comparison of radio and mid-infrared survey images. We present the first data release of Radio Galaxy Zoo, an online citizen science project that enlists the help of citizen scientists to cross-match extended radio sources from the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) and the Australia Telescope Large Area Survey (ATLAS) surveys, often with complex structure, to host galaxies in 3.6 um infrared images from the Wide-field Infrared Survey Explorer (WISE) and the Spitzer Space Telescope. This first data release consists of 100,185 classifications for 99,146 radio sources from the FIRST survey and 583 radio sources from the ATLAS survey. We include two tables for each of the FIRST and ATLAS surveys: 1)the identification of all components making up each radio source; and 2) the cross-matched host galaxies. These classifications have an average reliability of 0.83 based on the weighted consensus levels of our citizen scientists. The reliability of the DR1 catalogue has been further demonstrated through several parallel studies which used the pre-release versions of this catalogue to train and prototype machine learning-based classifiers. We also include a brief description of the radio source populations catalogued by RGZ DR1.

We analyzed a sample of $\sim$113,000 galaxies ($\rm z < 0.3$) from the Sloan Digital Sky Survey, divided into star-forming, composite, Seyfert, and LINER types, to explore the relationships between UV-to-optical colors ($\rm u-r$), star formation rates (SFRs), specific star formation rates (sSFRs), stellar velocity dispersions ($\rm \sigma_{*}$), mass accretion rates onto the black hole ($\rm L_{[OIII]}/\sigma_{*}^{4}$), and Eddington ratios. Star-forming galaxies predominantly feature young, blue stars along the main-sequence (MS) line, while composite, Seyfert, and LINER galaxies deviate from this line, displaying progressively older stellar populations and lower SFRs. $\rm L_{[OIII]}/\sigma_{*}^{4}$ and Eddington ratios are highest in Seyfert galaxies, moderate in composite galaxies, and lowest in LINERs, with higher ratios associated with bluer colors, indicating a younger stellar population and stronger active galactic nucleus (AGN) activity. These trends suggest a strong correlation between sSFRs and Eddington ratios, highlighting a close connection between AGN and star formation activities. These results may imply an evolutionary sequence where galaxies transition from blue star-forming galaxies to red LINERs, passing through composite and Seyfert phases, driven primarily by gas supply, with AGN feedback playing a secondary role. While both radio luminosities ($\rm L_{1.4GHz}$) and Eddington ratios correlate with SFRs, their trends differ on the SFR$-$stellar mass ($\rm M_{*}$) plane, with radio luminosities increasing with stellar mass along the MS line, and no direct connection between radio luminosities and Eddington ratios. These findings may provide new insights into the interplay between star formation, AGN activity, and radio emission in galaxies, shedding light on their evolutionary pathways.

As one type of blue early-type galaxies, the evolutionary history and fate of star-forming lenticular galaxies (S0s) remain elusive. We selected 134 star-forming S0s from the SDSS-IV MaNGA survey and found that they have steep and warped size-mass relations, similar to quiescent S0s and red spirals, indicating that they may have similar gas dissipation scenarios. These galaxies have a higher central stellar mass surface density than normal blue spirals. The radial profiles of $D_{\rm n}4000$ and [Mgb/Fe] show that red spirals and quiescent S0s have similar old central populations and high [Mgb/Fe] values, suggesting rapid bulge formation, though red spirals exhibit a steeper gradient possibly due to residual star formation (SF) in outer regions. In contrast, star-forming S0s exhibit profiles between quiescent S0s/red spirals and normal blue spirals, with relatively flat $D_{\rm n}4000$ and [Mgb/Fe] gradients. More long-term SF history causes normal blue spirals to have very flat $D_{\rm n}4000$ and [Mgb/Fe] profiles, and the majority of them (79 $\pm$ 5 $\%$) have S$\acute{\rm e}$rsic index $<$ 2. We also found that the halo mass of star-forming S0s resembles that of quiescent S0s/red spirals, with 82 $\pm$ 5 $\%$ exceeding the critical mass ($M_{\rm halo} = 10^{12}$$M_{\odot}$h$^{-1}$). To supplement previous H\,{\sc i} detection of star-forming S0s covered by H\,{\sc i}MaNGA, we obtained new observation for H\,{\sc i} emission from 41 star-forming S0s in our sample using the Five-hundred-meter Aperture Spherical Radio Telescope. We found that the H\,{\sc i} mass distribution of star-forming S0s matches that of normal blue spirals, although both star-forming S0s and red spirals are relatively gas-poor, resulting in varying atomic gas depletion times due to different SF levels. Based on these observational results, we discuss the possible evolutionary scenarios of star-forming S0s.

Accurate primary beam calibration is essential for precise brightness measurements in radio astronomy. The VLA Low-band Ionosphere and Transient Experiment (VLITE) faces challenges in calibration due to the offset Cassegrain optics used in its commensal observing system. This study aims to develop a novel calibration method to improve accuracy with no impact on the Very Large Array (VLA) primary science observations. We used the apparent brightness of standard candles identified in VLITE's commensal data to develop 1D and 2D primary beam response models. These models accounted for operational changes and asymmetries caused by the subreflector and were validated against holographic methods and compact source light curves. The models achieved calibration accuracy within 3 percent across the field of view, significantly improving the precision of brightness measurements. The results were consistent with holography-derived solutions and performed reliably under different operational conditions. This improved calibration technique expands VLITE's capabilities for studying active galactic nuclei, transients, and pulsars. It offers a cost-effective alternative to traditional holographic methods, facilitating broader use in commensal observing systems.

The detection and constraint of the orbits of long-period giant planets is essential for enabling their further study through direct imaging. Recently, Venner et al. (2024) highlighted discrepancies between the solutions presented by Feng et al. (2022) and those from other studies, which primarily use orvara. We address these concerns by reanalyzing the data for HD 28185, GJ 229, HD 211847, GJ 680, HD 111031, and eps Ind A, offering explanations for these discrepancies. Based on a comparison between the methods used by Feng et al. (2022) and orvara, we find the discrepancies are primarily data-related rather than methodology-related. Our re-analysis of HD 28185 highlights many of the data-related issues and particularly the importance of parallax modeling for year-long companions. The case of eps Ind A b is instructive to emphasize the value of an extended RV baseline for accurately determining orbits of long period companions. Our orbital solutions highlight other causes for discrepancies between solutions including the combination of absolute and relative astrometry, clear definitions of conventions, and efficient posterior sampling for the detection of wide-orbit giant planets.

The Tibet AS$\gamma$ experiment successfully detected sub-PeV $\gamma$-rays from the Crab nebula using a Surface Array and underground muon detector. Considering this, we are building in Bolivia a new experiment to explore the Southern Hemisphere, looking for the origins of cosmic rays in our Galaxy. The name of this project is Andes Large area PArticle detector for Cosmic ray physics and Astronomy (ALPACA). A prototype array called ALPAQUITA, with $1/4$ the total area of the full ALPACA, started observations in September $2022$. In this paper we introduce the status of ALPAQUITA and the plans to extend the array. We also report the results of the observation of the moon shadow in cosmic rays.

The observation of the kilonova AT2017gfo and investigations of its light curves and spectra confirmed that neutron star mergers are sites of r-process nucleosynthesis. However, the identification of elements responsible for the spectral features is still challenging, particularly at the near-infrared wavelengths. In this study, we systematically searched for all possible near-infrared transitions of heavy elements using experimentally calibrated energy levels. Our analysis reveals that most candidate elements with strong absorption lines are lanthanides (Z=57-71) and actinides (Z=89-103). This is due to their complex structures leading to many low-lying energy levels, which results in strong transitions in the near-infrared range. Domoto et al. (2022) have shown that La III and Ce III can explain the absorption features at $\lambda\sim$ 12,000 - 15,000 A. While our results confirm that these two elements show strong infrared features, we additionally identify Gd III as the next most promising species. Due to its unique atomic structure involving the half-filled 4f and the outer 5d orbitals, Gd III has one of the lowest-lying energy levels, between which relatively strong transitions occur. We also find absorption lines caused by Gd III in the near-infrared spectrum of a chemically peculiar star HR 465, which supports their emergence in kilonova spectra. By performing radiative transfer simulations, we confirm that Gd III lines affect the feature at $\sim$ 12,000 A previously attributed to La III. Future space-based time-series observations of kilonova spectra will allow the identification of Gd III lines.

The cosmological principle asserting the large-scale uniformity of the Universe is a testable assumption of the standard cosmological model. We explore the constraints on anisotropic expansion provided by measuring directional variation in the Hubble constant, $H_0$, derived from differential zeropoint measurements of the Tully-Fisher distance estimator. We fit various models for directional variation in $H_0$ using the Tully-Fisher dataset from the all-sky Cosmicflows-4 catalog. The best-fit dipole variation has an amplitude of 0.063 $\pm$ 0.016 mag in the direction ($\ell,b$) = (142 $\pm$ 30$^{\circ}$, 52 $\pm$ 10$^{\circ}$). If this were due to anisotropic expansion it would imply a 3% variation in $H_0$, corresponding to $\Delta H_0$ = 2.10 $\pm$ 0.53 km/s/Mpc if $H_0$ = 70 km/s/Mpc, with a significance of 3.9$\sigma$. A model that includes this $H_0$ dipole is only weakly favored relative to a model with a constant $H_0$ and a bulk motion of the volume sampled by Cosmicflows-4 that is consistent with the standard $\Lambda$CDM cosmology. However, we show that with the expected Tully-Fisher data from the WALLABY and DESI surveys it should be possible to detect a 1% $H_0$ dipole anisotropy at 5.8$\sigma$ confidence and to distinguish it from the typical bulk flow predicted by $\Lambda$CDM over the volume of these surveys.

Dark matter halos with self-interacting dark matter (SIDM) experience a unique evolutionary phenomenon, in that their central regions eventually collapse to high density through the runaway gravothermal process after initially forming a large and low-density core. When coupled with orbital evolution, this is expected to naturally produce a large diversity in dark-matter halos' inner mass distribution, potentially explaining the diversity problem of dwarf galaxies. However, it remains unknown how the diversity in SIDM dark-matter halos propagates to the more easily observed luminous matter at the center of the halo, especially the stellar component. In this work, we use idealized N-body simulations with two species of particles (dark matter and stars) to study the response of the stellar properties of field and satellite dwarf galaxies to SIDM evolution and orbital effects on their halos. Galaxies' stellar components, including galaxy size, mass-to-light ratio, and stellar velocity dispersion, display a much larger scatter in SIDM than the standard cold dark matter (CDM) model. Importantly, we find signs of universality in the evolution pathways, or ``tidal tracks'', of SIDM dwarf satellites, which are physically interpretable and potentially parameterizable. This type of tidal-track model can be layered onto larger-scale, cosmological simulations to reconstruct the evolution of populations of SIDM dwarfs in cases where high-resolution simulations of galaxies are otherwise prohibitively expensive.

We develop an analytic method of inverting the Tolman-Oppenheimer-Volkov (TOV) relations to high accuracy. This method is developed from the strong correlations that are shown to exist between the neutron star mass-radius curve and the equation of state (EOS) or pressure-energy density relation. Selecting points that have masses equal to fixed fractions of the maximum mass, we find a semi-universal power-law relation between the central energy densities, pressures, sound speeds, chemical potentials and number densities of those stars, with the maximum mass and the radii of one or more fractional maximum mass points. Root-mean-square fitting accuracies, for EOSs without large first-order phase transitions, are typically 0.5\% for all quantities at all mass points. The method also works well, although less accurately, in reconstructing the EOS of hybrid stars with first-order phase transitions. These results permit, in effect, an analytic method of inverting an arbitrary mass-radius curve to yield its underlying EOS. We discuss applications of this inversion technique to the inference of the dense matter EOS from measurements of neutron star masses and radii as a possible alternative to traditional Bayesian approaches.

Weak gravitational lensing requires precise measurements of galaxy shapes and therefore an accurate knowledge of the PSF model. The latter can be a source of systematics that affect the shear two-point correlation function. A key stake of weak lensing analysis is to forecast the systematics due to the PSF. Correlation functions of galaxies and the PSF, the so-called $\rho$- and $\tau$-statistics, are used to evaluate the level of systematics coming from the PSF model and PSF corrections, and contributing to the two-point correlation function used to perform cosmological inference. Our goal is to introduce a fast and simple method to estimate this level of systematics and assess its agreement with state-of-the-art approaches. We introduce a new way to estimate the covariance matrix of the $\tau$-statistics using analytical expressions. The covariance allows us to estimate parameters directly related to the level of systematics associated with the PSF and provides us with a tool to validate the PSF model used in a weak-lensing analysis. We apply those methods to data from the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS). We show that the semi-analytical covariance yields comparable results than using covariances obtained from simulations or jackknife resampling. It requires less computation time and is therefore well suited for rapid comparison of the systematic level obtained from different catalogs. We also show how one can break degeneracies between parameters with a redefinition of the $\tau$-statistics. The methods developed in this work will be useful tools in the analysis of current weak-lensing data but also of Stage IV surveys such as Euclid, LSST or Roman. They provide fast and accurate diagnostics on PSF systematics that are crucial to understand in the context of cosmic shear studies.

The Ly$\alpha$ emission line is a characteristic feature found in high-$z$ galaxies, serving as a probe of cosmic reionization. While previous works present various correlations between Ly$\alpha$ emission and physical properties of host galaxies, it is still unclear which characteristics predominantly determine the Ly$\alpha$ emission. In this study, we introduce a neural network approach to simultaneously handle multiple properties of galaxies. The neural-network-based prediction model that identifies Ly$\alpha$ emitters (LAEs) from six physical properties: star formation rate (SFR), stellar mass, UV absolute magnitude $M_\mathrm{UV}$, age, UV slope $\beta$, and dust attenuation $E(B-V)$, obtained by the SED fitting. The network is trained with galaxy samples from the VANDELS and MUSE spectroscopic surveys and achieves the performance of 77% true positive rate and 14% false positive rate. The permutation feature importance method shows that $\beta$, $M_\mathrm{UV}$, and $M_*$ are important for the prediction of LAEs. As an independent validation, we find that 91% of LAEs spectroscopically confirmed by the James Webb Space Telescope (JWST) have a probability of LAE higher than 70% in this model. This prediction model enables the efficient construction of a large LAE sample in a wide and continuous redshift space using only photometric data. We apply the prediction model to the JWST photometric galaxy sample and obtain Ly$\alpha$ fraction consistent with previous studies. Moreover, we demonstrate that the difference between the distributions of LAEs predicted by the model and the spectroscopically identified LAEs provides a strong constraint on the HII bubble size.

Merger events are thought to be an important phase in the assembly of massive galaxies. At the same time, Active Galactic Nuclei (AGN) play a fundamental role in the evolution of their star formation histories. Both phenomena can be observed at work in NGC 6240, a local prototypical merger, classified as an UltraLuminous InfraRed Galaxy (ULIRG) thanks to its elevated infrared luminosity. Interestingly, NGC 6240 hosts two AGN separated by 1.5''(~ 735 pc), detected in both X-ray and radio band. Taking advantage of the unprecedented sensitivity and wavelength coverage provided by the Integral Field Unit (IFU) of the NIRSpec instrument onboard JWST, we observed the nuclear region of NGC 6240 in a FoV of 3.7'' x 3.7''(1.9 x 1.9 kpc^2), to investigate gas kinematics and InterStellar Medium (ISM) properties with a high spatial resolution of ~ 0.1'' (or ~ 50 pc). We separated the different gas kinematic components through multi-Gaussian fitting and studied the excitation properties of the ISM from the NIR diagnostic diagram based on the H_2 1-0 S(1)/BrGamma and [Fe II]1.257micron/PaBeta lines ratios. We isolated the ionization cones of the two nuclei, and detected coronal lines emission from both of them. Using H_2 line ratios, we found that the molecular hydrogen gas is excited mostly by thermal processes. We computed a hot molecular gas mass of 1.3 x 10^5 M_sun and an ionized gas mass in the range of 10^5 - 10^7 M_sun. We studied with unprecedented spatial resolution and sensitivity the kinematics of the molecular and ionized gas phases. We revealed the complex structure of the molecular gas and found a blueshifted outflow near the Southern nucleus, together with filaments connecting a highly redshifted H_2 cloud with the two nuclei. We speculate on the possible nature of this H_2 cloud and propose two possible scenarios: either outflowing gas, or a tidal cloud falling onto the nuclei.

Weak lensing surveys, along with most other late-Universe probes, have consistently measured a lower amplitude of the matter fluctuation spectrum, denoted by the parameter $S_8$, compared to predictions from early-Universe measurements in cosmic microwave background data. Improper modelling of nonlinear scales may partially explain these discrepancies in lensing surveys. This study investigates whether the conventional approach to addressing small-scale biases remains optimal for Stage-IV lensing surveys. We demonstrate that conventional weak lensing estimators are affected by scale leakage from theoretical biases at nonlinear scales, which influence all observed scales. Using the BNT transform, we propose an $\ell$-cut methodology that effectively controls this leakage. The BNT transform reorganizes weak lensing data in $\ell$ space, aligning it with $k$ space, thereby reducing the mixing of nonlinear scales and providing a more accurate interpretation of the data. We evaluate the BNT approach by comparing HMcode, Halofit, Baryon Correction Model and AxionHMcode mass power spectrum models using Euclid-like survey configurations. Additionally, we introduce a new estimator to quantify scale leakage in both the BNT and noBNT approaches. Our findings show that BNT outperforms traditional methods, preserving cosmological constraints while significantly mitigating theoretical biases.

As part of the guaranteed time observations program Mid-Infrared Characterization Of Nearby Iconic galaxy Centers (MICONIC), we used the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board of the JWST to study the nearby merger NGC6240. We aim to characterise the dual active galactic nuclei (AGN), the ionised gas outflows and the main properties of the interstellar medium over a mapped area of 6.6"x7.7". We obtained integral field spectroscopic mid-infrared data of NGC6240, resolving both nuclei for the first time in the full 5-28{\mu}m spectral range. We modelled the emission lines through a kinematic decomposition, finding that the fine-structure lines in the southern (S) nucleus are broader than for the northern (N) nucleus (full width at half maximum of $\geq$1500 vs ~700 km s^{-1} on average). High excitation lines, such as [NeV], [NeVI], and [MgV], are clearly detected in the N nucleus. In the S nucleus, the same lines can be detected but only after a decomposition of the PAH features in the integrated spectrum, due to a combination of a strong mid-IR continuum, the broad emission lines, and the intense star formation (SF). The SF is distributed all over the mapped FoV of 3.5kpc x 4.1kpc (projected), with the maximum located around the S nucleus. Both nuclear regions appear to be connected by a bridge region detected with all the emission lines. Based on the observed MRS line ratios and the high velocity dispersion ({\sigma}~600 km s^{-1}), shocks are also dominating the emission in this system. We detected the presence of outflows as a bubble north-west from the N nucleus and at the S nucleus. We estimated a ionised mass outflow rate of 1.4$\pm$0.3 M yr^{-1} and 1.8$\pm$0.2 M yr^{-1}, respectively. Given the derived kinetic power of these outflows, both the AGN and the starburst could have triggered them. [Abridged]

Measurements of weak gravitational lensing using the cosmic microwave background and the shapes of galaxies have refined our understanding of the late-time history of the Universe. While optical surveys have been the primary source for cosmic shear measurements, radio continuum surveys offer a promising avenue. Relevant radio sources, principally star-forming galaxies, have populations with higher mean redshifts and are less affected by dust extinction compared to optical sources. We focus on the future mid frequency SKA radio telescope and explore the cross-correlation between radio cosmic shear and CMB lensing convergence ($\gamma_\mathrm{R}\times \kappa_\mathrm{CMB}$). We investigate its potential in constraining the redshift distribution of radio galaxy samples and improving cosmological parameter constraints, including the neutrino sector. Using simulations of the first phase of the SKA and the Simons Observatory as a CMB experiment, we show how this $\gamma_\mathrm{R}\times \kappa_\mathrm{CMB}$ cross-correlation can provide $\sim1 - 10\%$ calibration of the overall radio source redshift distribution, which in turn can significantly tighten otherwise degenerate measurements of radio galaxy bias. For the case of the next-generation full SKA, we find that the cross-correlation becomes more powerful than the equivalent with a \textit{Euclid}-like survey, with constraints $30\%$ tighter on $\Lambda$CDM parameters and narrower bounds on sum of neutrino masses at the level of $\sim 24\%$. These constraints are also driven by higher redshifts and larger scales than other galaxy-CMB cross-correlations, potentially shedding light on different physical models. Our findings demonstrate the potential of radio weak lensing in improving constraints, and establish the groundwork for future joint analyses of CMB experiments and radio continuum surveys.

Understanding the structure of the Galactic disk is crucial for understanding the formation and evolutionary history of the Milky Way. This study examines the structure of the Galactic disk by analyzing a sample of 138,667 primary red clump (RC) stars from the LAMOST and Gaia datasets. We have categorized these RC stars into mono-age populations and investigated their spatial distributions within the R - Z plane, estimating scale heights and lengths through the fitting of their vertical and radial density profiles. Our analysis indicates that the vertical profiles of these mono-age populations fit a dual-component disk model, where both components exhibit significant flaring, particularly in the outer disk regions. Within a constant Galactocentric radius R, the scale heights of the first component, representing the morphologically thin disk, rise with age. In contrast, the scale heights of the second component, corresponding to the morphologically thick disk, remain comparatively stable across different age groups. Additionally, the radial density profiles of both disk components predominantly peak within a radial range of 7.5-8.5 kpc. These findings underscore the importance of age as a crucial factor in shaping the spatial distribution and structural evolution of the Galactic disk, offering valuable insights into its complex dynamics and history.

Conventional calibration of Baryon Acoustic Oscillations (BAO) data relies on estimation of the sound horizon at drag epoch $r_d$ from early universe observations by assuming a cosmological model. We present a recalibration of two independent BAO datasets, SDSS and DESI, by employing deep learning techniques for model-independent estimation of $r_d$, and explore the impacts on $\Lambda$CDM cosmological parameters. Significant reductions in both Hubble ($H_0$) and clustering ($S_8$) tensions are observed for both the recalibrated datasets. Moderate shifts in some other parameters hint towards further exploration of such data-driven approaches.

We search for an excess of electrons and positrons in the interplanetary space from the decays of heavy neutrinos produced in nuclear reactions in the Sun. Using measurements of the electron spectra in the MeV range from the Ulysses and SOHO satellites, we report the strongest direct upper bound to date on the mixing between heavy neutral leptons with MeV masses and electron neutrinos, reaching $U_e^2\simeq 10^{-6}$ at $M_N=10$ MeV. Our sensitivity is predominantly constrained by the uncertainties in the propagation of electrons and positrons, particularly the diffusion coefficient in the inner Solar System, as well as the uncertainties in the astrophysical background. Enhancing our understanding of either of these factors could lead to a significant improvement in sensitivity.

We present new radio images of the supernova remnant (SNR) CTA 1 at 1420 and 408 MHz, and in the 21 cm line of H I observed with the Dominion Radio Astrophysical Observatory Synthesis Telescope and at 1420 MHz observed with the Effelsberg 100 m telescope. We confirm previously described continuum features and elaborate further on filamentary features identified using the high-resolution (1') maps from these new observations. We investigate the abrupt change in sign of rotation measure (RM) across the SNR, using the linear polarization observations in the four bands around 1420 MHz. Following X. H. Sun et al.'s (2011) investigation, we both confirm that the distribution of signs of the RMs for extragalactic sources in the area appears to match that of the shell, as well as combine the data from the four bands to estimate the relative depolarization and the intrinsic rotation measure of the SNR. We do not conclusively reject X. H. Sun et al.'s (2011) claim of a Faraday screen in the foreground causing the distribution of RMs that we observe; however, we do suggest an alternative explanation of a swept-up stellar wind from the progenitor star with a toroidal magnetic field. Finally, we expand on the analysis of the H I observations by applying the Rolling Hough Transform to isolate filamentary structure and better identify H I emission with the SNR. Further constraining the H I velocity channels associated with CTA 1, we use more recent Galactic rotation curves to calculate an updated kinematic distance of 1.09 +/- 0.2 kpc.

Arp 220 is the prototypical Ultraluminous Infrared Galaxy (ULIRG), and one of the brightest objects in the extragalactic far-infrared sky. It is the result of a merger between two gas rich spiral galaxies which has triggered starbursting activity in the merger nuclear regions. Observations with the Submillimeter Array centred at a frequency of 345 GHz and with a synthesised beamsize of 0.77 x 0.45 arcseconds were used to search for polarized dust emission from the nuclear regions of Arp 220. Polarized dust emission was clearly detected at 6 sigma significance associated with the brighter, western nucleus, with a peak polarization fraction of 2.7 +/- 0.35 per cent somewhat offset from the western nucleus. A suggestive 2.6 sigma signal is seen from the fainter eastern nucleus. The dust emission polarization is oriented roughly perpendicular to the molecular disk in the western nucleus suggesting that the magnetic field responsible is orientated broadly in the plane of the disk, but may be being reordered by the interaction between the two nuclei. Unlike more evolved interacting systems, we see no indication that the magnetic field is being reordered by the outflow from the western nucleus. These observations are the first detection of dust polarization, and thus of magnetic fields, in the core of a ULIRG.

Radio frequency interference (RFI) is a persistent contaminant in terrestrial radio astronomy. While new radio interferometers are becoming operational, novel sources of RFI are also emerging. In order to strengthen the mitigation of RFI in modern radio interferometers, we propose an on-line RFI mitigation scheme that can be run in the correlator of such interferometers. We combine statistics based on the energy as well as the polarization alignment of the correlated signal to develop an on-line RFI mitigation scheme that can be applied to a data stream produced by the correlator in real-time, especially targeted at low duty-cycle or transient RFI detection. In order to improve the computational efficiency, we explore the use of both single precision and half precision floating point operations in implementing the RFI mitigation algorithm. This ideally suits its deployment in accelerator computing devices such as graphics processing units (GPUs) as used by the LOFAR correlator. We provide results based on real data to demonstrate the efficacy of the proposed method.

The driving and excitation mechanisms of decay-less kink oscillations in coronal loops remain under debate. We aim to quantify and provide simple observational constraints on the photospheric driving of oscillating coronal loops in a few typical active region configurations: sunspot, plage, pores and enhanced-network regions. We then aim to investigate the possible interplay between photospheric driving and properties of kink oscillations in connected coronal loops. We analyse two unique datasets of the corona and photosphere taken at a high resolution during the first coordinated observation campaign between Solar Orbiter and the Swedish 1-m Solar Telescope (SST). A local correlation tracking method is applied on the SST/CRISP data to quantify the photospheric motions at the base of coronal loops. The same loops are then analysed in the corona by exploiting data from the Extreme Ultraviolet Imager on Solar Orbiter, and by using a wavelet analysis to characterize the kink oscillations. Each photospheric region shows dynamics with an overall increase in strength going from pore, plage, enhanced-network to sunspot regions. Differences are also seen in the kink-mode amplitudes of the corresponding coronal loops. This suggests the photosphere is involved in the driving of coronal kink oscillations. However, the few samples available does not allow to further establish the excitation mechanism yet. Despite oscillating coronal loops being anchored in seemingly "static" strong magnetic field regions as seen from coronal EUV observations, photospheric observations provide evidence for a continuous and significant driving at their base. The precise connection between photospheric driving and coronal kink oscillations remains to be further investigated. This study finally provides critical constraints on photospheric driving that can be tested in existing numerical models of coronal loops.

Magnetic fields present in the Large Scale Structure (LSS) of the Universe change polarization of radio waves arriving from distant extragalactic sources through the effect of Faraday rotation. This effect has been recently used to detect magnetic field in the LSS filaments based on the Rotation Measure data of the LOFAR Two-Meter Sky Survey (LoTSS). We notice that the same data also constrain the strength of the volume-filling magnetic field in the voids of the LSS. We use the LoTSS data to to derive an improved upper bound on the volume-filling field. The new upper bound provides an order of magnitude improvement on the previous Faraday rotation bounds. The new Faraday Rotation bound on the scale-invariant field that may originate from the epoch of inflation is also an order of magnitude lower than the bound on such field derived from the anisotropy analysis of the Cosmic Microwave Background.

The {\it Javalambre Photometric Local Universe Survey} (J-PLUS) is a {\it spectro-photometric} survey covering about 3,000~deg$^2$ in its third data release (DR3), and containing about 300,000 galaxies with high quality ({\it odds}$>0.8$) photometric redshifts (hereafter photo-$z$s). We use this galaxy sample to conduct a tomographic study of the counts and redshift angular fluctuations under Gaussian shells sampling the redshift range $z\in[0.05,0.25]$. We confront the angular power spectra of these observables measured under shells centered on 11 different redshifts with theoretical expectations derived from a linear Boltzmann code ({\tt ARFCAMB}). Overall we find that J-PLUS DR3 data are well reproduced by our linear, simplistic model. We obtain that counts (or density) angular fluctuations (hereafter ADF) are very sensitive to the linear galaxy bias $b_g(z)$, although weakly sensitive to radial peculiar velocities of the galaxy field, while suffering from systematics residuals for $z>0.15$. Angular redshift fluctuations (ARF), instead, show higher sensitivity to radial peculiar velocities and also higher sensitivity to the average uncertainty in photo-$z$s ($\sigma_{\rm Err}$), with no obvious impact from systematics. For $z<0.15$ both ADF and ARF agree on measuring a monotonically increasing linear bias varying from $b_g(z=0.05)\simeq 0.9\pm 0.06$ up to $b_g(z=0.15)\simeq 1.5\pm 0.05$, while, by first time, providing consistent measurements of $\sigma_{\rm Err}(z)\sim 0.014$ that are $\sim 40~\%$ higher than estimates from the photo-$z$ code {\tt LePhare}, ($\sigma_{\rm Err}^{\rm LePhare}=0.010$). As expected, this photo-$z$ uncertainty level prevents the detection of radial peculiar velocities in the modest volume sampled by J-PLUS DR3, although prospects for larger galaxy surveys of similar (and higher) photo-$z$ precision are promising.

Observational and/or astrophysical systematics modulating the observed number of luminous tracers can constitute a major limitation in the cosmological exploitation of surveys of the large scale structure of the universe. Part of this limitation arises on top of our ignorance on how such systematics actually impact the observed galaxy/quasar fields. In this work we develop a generic, hybrid model for an arbitrary number of systematics that may modulate observations in both an additive and a multiplicative way. This model allows us devising a novel algorithm that addresses the identification and correction for either additive and/or multiplicative contaminants. We test this model on galaxy mocks and systematics templates inspired from data of the third data release of the {\it Javalambre Photometric Local Universe Survey} (J-PLUS). We find that our method clearly outperforms standard methods that assume either an additive or multiplicative character for all contaminants in scenarios where both characters are actually acting on the observed data. In simpler scenarios where only an additive or multiplicative imprint on observations is considered, our hybrid method does not lie far behind the corresponding simplified, additive/multiplicative methods. Nonetheless, in scenarios of mild/low impact of systematics, we find that our hybrid approach converges towards the standard method that assumes additive contamination, as predicted by our model describing systematics. Our methodology also allows for the estimation of biases induced by systematics residuals on different angular scales and under different observational configurations, although these predictions necessarily restrict to the subset of {\em known/identified} potential systematics, and say nothing about ``unknown unknowns" possibly impacting the data.

Astrometric microlensing events occur when a massive object passes between a distant source and the observer, causing a shift of the light centroid. The precise astrometric measurements of the Gaia mission provide an unprecedented opportunity to detect and analyze these events, revealing properties of lensing objects such as their mass and distance. We develop and test the Gaia Astrometric Microlensing Events (GAME) Filter, a software tool to identify astrometric microlensing events and derive lensing object properties. We generated mock Gaia observations for different magnitudes, number of Gaia visits, and events extending beyond Gaia's observational run. We applied GAME Filter to these datasets and validated its performance. We also assessed the rate of false positives where binary astrometric systems are misidentified as microlensing events. GAME Filter successfully recovers microlensing parameters for strong events. Parameters are more difficult to recover for short events and those extending beyond Gaia's run, where only a fraction of the events is observed. The astrometric effect breaks the degeneracy in the microlensing parallax present in photometric microlensing. For fainter sources, the observed signal weakens, reducing recovered events and increasing parameter errors. However, even for Gaia G-band magnitude 19, parameters can be recovered for Einstein radii above two mas. Observing regions with varying numbers of Gaia visits has minimal impact on filter accuracy when the number of visits exceeds 90. Additionally, even if the peak of a microlensing event lies outside Gaia's run, microlensing parameters can still be recovered. We also tested the sensitivity to contamination and found that 5 percent of binary systems were misclassified. GAME Filter characterizes lenses with astrometry-only data for lens masses from approximately 1 to 20 solar masses and distances up to 6 kpc.

In recent years, complementary metal-oxide-semiconductor (CMOS) sensors have been demonstrated to have significant potential in X-ray astronomy, where long-term reliability is crucial for space X-ray telescopes. This study examines the long-term stability of a scientific CMOS sensor, focusing on its bias, dark current, readout noise, and X-ray spectral performance. The sensor was initially tested at -30 $^\circ$C for 16 months, followed by accelerated aging at 20 $^\circ$C. After a total aging period of 610 days, the bias map, dark current, readout noise, gain, and energy resolution exhibited no observable degradation. There are less than 50 pixels within the 4 k $\times$ 4 k array which show a decrease of the bias under 50 ms integration time by over 10 digital numbers (DNs). First-order kinetic fitting of the gain evolution predicts a gain degeneration of 0.73% over 3 years and 2.41% over 10 years. These results underscore the long-term reliability of CMOS sensors for application in space missions.

Time-delay interferometry (TDI) suppresses laser frequency noise by forming linear combinations of time-shifted interferometric measurements. The time-shift operation is implemented by interpolating discretely sampled data. To enable in-band laser noise reduction by eight to nine orders of magnitude, interpolation has to be performed with high accuracy. Optimizing the design of those interpolation methods is the focus of this work. Previous research that studied constant time-shifts suggested Lagrange interpolation as the interpolation method for TDI. Its transfer function performs well at low frequency but requires a high number of coefficients. Furthermore, when applied in TDI we observed prominent time-domain features when a time-varying shift scanned over a pure integer sample shift. To limit this effect we identify an additional requirement for the interpolation kernel: when considering time-varying shifts the interpolation kernel must be sufficiently smooth to avoid unwanted time-domain transitions that produce glitch-like features in power spectral density estimates. The Lagrange interpolation kernel exhibits a discontinuous first derivative by construction, which is insufficient for the application to LISA or other space-based GW observatories. As a solution we propose a novel design method for interpolation kernels that respect a predefined requirement on in-band interpolation residuals and that possess continuous derivatives up to a prescribed order. Using this method we show that an interpolation kernel with 22 coefficients is sufficient to respect LISA's picometre-requirement and to allow for a continuous first derivative which suppresses the magnitude of the time-domain transition adequately. The reduction from 42 (Lagrange interpolation) to 22 coefficients enables us to save computational cost and increases robustness against artefacts in the data.

The direct collapse scenario, which predicts the formation of supermassive stars (SMSs) as precursors to supermassive black holes (SMBHs), has been explored primarily under the assumption of metal-free conditions. However, environments exposed to strong far-ultraviolet (FUV) radiation, which is another requirement for the direct collapse, are often chemically enriched to varying degrees. In this study, we perform radiation hydrodynamic simulations of star-cluster formation in clouds with finite metallicities, $Z=10^{-6}$ to $10^{-2} Z_{\odot}$, incorporating detailed thermal and chemical processes and radiative feedback from forming stars. Extending the simulations to approximately two million years, we demonstrate that SMSs with masses exceeding $10^4~M_\odot$ can form even in metal-enriched clouds with $Z \lesssim 10^{-3} Z_{\odot}$. The accretion process in these cases, driven by "super-competitive accretion," preferentially channels gas into central massive stars in spite of small (sub-pc) scale fragmentation. At $Z \simeq 10^{-2} Z_{\odot}$, however, enhanced cooling leads to intense fragmentation on larger scales, resulting in the formation of dense star clusters dominated by very massive stars with $10^3 M_{\odot}$ rather than SMSs. These clusters resemble young massive or globular clusters observed in the distant and local universe, exhibiting compact morphologies and high stellar surface densities. Our findings suggest that SMS formation is viable below a metallicity threshold of approximately $10^{-3} Z_{\odot}$, significantly increasing the number density of massive seed black holes to levels sufficient to account for the ubiquitous SMBHs observed in the local universe. Moreover, above this metallicity, this scenario naturally explains the transition from SMS formation to dense stellar cluster formation.

In 2017, as the solar cycle approached solar minimum, an unusually long and large depression was observed in galactic cosmic ray (GCR) protons, detected with the Alpha Magnetic Spectrometer (AMS-02), lasting for the second half of that year. The depression, as seen in the Bartel rotation-averaged proton flux, has the form of a Forbush decrease (FD). Despite this resemblance, however, the cause of the observed depression does not have such a simple explanation as FDs, due to coronal mass ejections (CMEs), typically last for a few days at 1 AU rather than half a year. In this work, we seek the cause of the observed depression and investigate two main possibilities. First, we consider a mini-cycle - a temporary change in the solar dynamo that changes the behavior of the global solar magnetic field and, by this, the modulation of GCRs. Secondly, we investigate the behavior of solar activity, both CMEs and co-rotating/stream interactions regions (C/SIRs), during this period. Our findings show that, although there is some evidence for mini-cycle behavior prior to the depression, the depression is ultimately due to a combination of recurrent CMEs, SIRs and CIRs. A particular characteristic of the depression is that the largest impacts that help to create and maintain it are due to four CMEs from the same, highly active, magnetic source that persists for several solar rotations. This active magnetic source is unusual given the closeness of the solar cycle to solar minimum, which also helps to make the depression more evident.

Recent observational results from asteroseismic studies show that an important fraction of solar-like stars do not present detectable stochastically excited acoustic oscillations. This non-detectability seems to correlate with a high rotation rate in the convective envelope and a high surface magnetic activity. At the same time, the properties of stellar convection are affected by rotation and magnetism. We investigate the role of rotation in the excitation of acoustic modes in the convective envelope of solar-like stars, to evaluate its impact on the energy injected in the oscillations. We derive theoretical prescriptions for the excitation of acoustic waves in the convective envelope of rotating solar-like stars. We adopt the Rotating Mixing-Length Theory to model the influence of rotation on convection. We use the MESA stellar evolution code and the GYRE stellar oscillation code to estimate the power injected in the oscillations from our theoretical prescriptions. We demonstrate that the power injected in the acoustic modes is insensitive to the rotation if a Gaussian time-correlation function is assumed, while it can decrease by up to 60 % for a Lorentzian time-correlation function, for a $20 \Omega_{\odot}$ rotation rate. This result can allow us to better constrain the properties of stellar convection by studying observationally acoustic modes excitation. These results demonstrate how important it is to take into account the modification of stellar convection by rotation when evaluating the amplitude of the stellar oscillations it stochastically excites. They open the path for understanding the large variety of observed acoustic-mode amplitudes at the surface of solar-like stars as a function of surface rotation rates.

We revisit the origin of the observed misaligned rings in the circumtriple disk around GW Ori. Previous studies appeared to disagree on whether disk breaking is caused by the differential precession driven in the disk by the triple star system. In this letter, we show that the previous studies are in agreement with each other when using the same set of parameters. But for observationally motivated parameters of a typical protoplanetary disk, the disk is unlikely to break due to interactions with the triple star system. We run 3-dimensional hydrodynamical simulations of a circumtriple disk around GW Ori with different disk aspect ratios. For a disk aspect ratio typical of protoplanetary disks, $H/r \gtrsim 0.05$, the disk does not break. An alternative scenario for the gap's origin consistent with the expected disk aspect ratio involves the presence of giant circumtriple planets orbiting GW Ori.

We present a structural analysis of 521 galaxy candidates at 6.5 < z < 12.5, with $SNR > 10\sigma$ in the F444W filter, taken from the EPOCHS v1 sample, consisting of uniformly reduced deep JWST NIRCam data, covering the CEERS, JADES GOOD-S, NGDEEP, SMACS0723, GLASS and PEARLS surveys. We use standard software to fit single S\'ersic models to each galaxy in the rest-frame optical and extract their parametric structural parameters (S\'ersic index, half-light radius and axis-ratio), and \texttt{Morfometryka} to measure their non-parametric concentration and asymmetry parameters. We find a wide range of sizes for these early galaxies, but with a strong galaxy-size mass correlation up to $z \sim 12$ such that galaxy sizes continue to get progressively smaller in the high-redshift regime, following $R_{e} = 2.74 \pm 0.49 \left( 1 + z \right) ^{-0.79 \pm 0.08}$ kpc. Using non-parametric methods we find that galaxy merger fractions, classified through asymmetry parameters, at these redshifts remain consistent with those in literature, maintaining a value of $f_{m} \sim 0.12 \pm 0.07$ showing little dependence with redshift when combined with literature at $z > 4$. We find that galaxies which are smaller in size also appear rounder, with an excess of high axis-ratio objects. Finally, we artificially redshift a subsample of our objects to determine how robust the observational trends we see are, determining that observed trends are due to real evolutionary effects, rather than being a consequence of redshift effects.

Stellar spin is one of the fundamental quantities that characterize a star itself and its planetary system. Nevertheless, stellar spin-down mechanisms in protostellar and pre-main-sequence stellar phases have been a long-standing problem in the star formation theory. To realize the spin-down, previous axisymmetric models based on the conventional magnetospheric paradigm have to assume massive stellar winds or produce highly time-variable magnetospheric ejections. However, this picture has been challenged by both numerical simulations and observations. With a particular focus on the propeller regime for solar-mass stars, we propose a new picture of stellar spin-down based on our recent three-dimensional (3D) magnetohydrodynamic simulation and stellar evolution calculation. We show that failed magnetospheric winds, unique to 3D models, significantly reduce the spin-up accretion torque, which make it easier for the star to spin-down. Additionally, the amplitude of time variability associated with magnetospheric ejections is reduced by 3D effects. Our simulation demonstrates that the star spins down by generating a conical disk wind, driven by a rotating stellar magnetosphere. Our theoretical estimates, inspired by the numerical model, suggest that the conical disk wind is likely to play a crucial role in extracting stellar angular momentum during the protostellar phase. As magnetospheric accretion is expected to occur in other accreting objects such as proto-giant planets, this study will also contribute to the understanding of the angular momentum of such objects.

We present new radio observations of the galaxy cluster merger CIZA J0107.7+5408 (CIZA0107), a large, roughly equal mass, post-core passage, dissociative binary system at z = 0.1066. CIZA0107 is an elongated, disturbed system, hosting two subclusters with optical galaxy number density peaks offset from their associated X-ray density peaks and double-peaked diffuse radio structure. We present new 240-470 MHz and 2.0-4.0 GHz Very Large Array observations of CIZA0107. We image the diffuse emission at high resolution, constrain its integrated spectrum, and map the spectral index distribution. We confirm the presence of steep-spectrum ($\alpha$ $\sim$ -1.3) emission on a scale of about 0.5 Mpc in both subclusters. We identify two smaller ultra-steep spectrum ($\alpha$ $>$ -2) regions, superimposed on larger-scale radio emission associated with the southwestern subcluster. At 340 MHz, we detect a radio edge bounding the emission to the south and show that it is coincident with a weak (M $\sim$ 1.2) shock identified in the Chandra image. At 3 GHz, the emission does not show any corresponding edge-like feature, and in fact it extends beyond the shock. We investigate the nature of the emission in CIZA0107 and find that, while the system may host a double halo structure, we cannot rule out a scenario in which the emission arises from two relics projected on the central cluster regions.

We present JWST/NIRSpec IFU observations of the z=7.152 galaxy system B14-65666, as part of the GA-NIFS survey. Line and continuum emission in this massive system (log10(M*/Msol)=9.8+/-0.2) is resolved into two strong cores, two weaker clumps, and a faint arc, as seen in recent JWST/NIRCam imaging. Our dataset contains detections of [OII]3727,3729, [NeIII]3869,3968, Balmer lines (HBeta, HGamma, HDelta, HEpsilon, HZeta), [OIII]5007, and weak [OIII]4363. Each spectrum is fit with a model that consistently incorporates interstellar medium conditions (i.e., electron temperature, T_e, electron density, n_e, and colour excess, E(B-V)). The resulting line fluxes are used to constrain the gas-phase metallicity (~0.3-0.4 solar) and HBeta-based SFR (310+/-40 Msol/yr) for each region. Common line ratio diagrams (O32-R23, R3-R2, Ne3O2-R23) reveal that each line-emitting region lies at the intersection of local and high-redshift galaxies, suggesting low ionisation and higher metallicity compared to the predominantly lower-mass galaxies studied with the JWST/NIRSpec IFU so far at z>5.5. Spaxel-by-spaxel fits reveal evidence for both narrow (FWHM<400 km/s) and broad (FWHM >500 km/s) line emission, the latter of which likely represents tidal interaction or outflows. Comparison to ALMA [CII]158um and [OIII]88um data shows a similar velocity structure, and optical-far infrared diagnostics suggest regions of high Lyman continuum escape fraction and n_e. This source lies on the mass-metallicity relation at z>4, suggesting an evolved nature. The two core galaxies show contrasting properties (e.g., SFR, M*, gas-phase metallicity), suggesting distinct evolutionary pathways. Combining the NIRSpec IFU and ALMA data sets, our analysis opens new windows into the merging system B14-65666.

This paper investigates the linear response of a series of spheroidal stellar clusters, the Kuzmin-Kutuzov St\"ackel family, which exhibit a continuous range of flattening and rotation, extending from an isochrone sphere to a Toomre disk. The method successfully replicates the growing modes previously identified in published $N$-body simulations. It relies on the efficiency of the matrix method to quantify systematically the effects of rotation and flattening on the eigenmodes of the galaxy. We identify two types of bi-symmetric instabilities for the flatter models - the so-called bending and bar-growing modes - the latter of which persists even for very round models. As anticipated, in its least unstable configurations, the system becomes flatter as its rotational speed increases. More realistic equilibria will be required to achieve a better match to the main sequence of fast-slow rotators. The corresponding code is made public.

We develop a method for constraining the potential of the Milky Way using stellar streams with a known progenitor. The method expresses the stream in angle-action coordinates and integrates the orbits of the stars backwards in time to obtain the stripping point positions of the stream stars relative to the cluster. In the potential that generated the stream, the stars return approximately to the cluster centre. In a different potential, they are redirected to different locations. The free parameters of the model are estimated by maximising the degree of clustering of the stripping point distribution. We test this method with the stellar stream of the globular cluster M68 (NGC 4590). We use an N-body code to simulate the stream and generate a realistic star sample using a model of the Gaia selection function. We also simulate the expected observational uncertainties, and estimate the heliocentric distances and radial velocities of the stream stars from the cluster orbit. Using this sample of stars, we recover the values of four free parameters characterising the potential of the disc and the dark halo to an accuracy of about 10 per cent of the correct values. We show that this accuracy is improved up to about 2 per cent using the expected end-of-mission Gaia data. In addition, we obtain a strong correlation between the mass of the disc and the dark halo axis ratio, which is explained by the fact that the stream flows close and parallel to the disc plane.

The field of air shower physics, dedicated to understanding the development of cosmic-ray interactions with the Earth's atmosphere, faces a significant challenge regarding the muon content of air showers observed by the Pierre Auger Observatory, and numerous other observatories. Thorough comparisons between extensive air shower (EAS) measurements and simulations are imperative for determining the primary energy and mass of ultra-high energy cosmic rays. Current simulations employing state-of-the-art hadronic interaction models reveal a muon deficit compared to experimental measurements, commonly known as the "Muon Puzzle". The primary cause of this deficit lies in the uncertainties surrounding high-energy hadronic interactions. In this contribution, we discuss the integration of a new hadronic interaction model, Pythia 8, into the effort to resolve the Muon Puzzle. While the Pythia 8 model is well-tailored in the context of Large Hadron Collider (LHC) experiments, its application in air shower studies remained limited until now. However, recent advancements, particularly in the Angantyr model of Pythia 8, offer promising enhancements in describing hadron-nucleus interactions, thereby motivating its potential application in air shower simulations. We present results from EAS simulations conducted using CORSIKA 8, wherein Pythia is employed to model hadronic interactions.

For over two decades, CORSIKA 7 and its previous versions have been the leading Monte Carlo code for simulating extensive air showers. However, its monolithic Fortran-based software design and hand-optimized code has created challenges for maintenance, adaptation to new computing paradigms, and extensions for more complex simulations. Addressing these limitations, the CORSIKA 8 project represents a comprehensive rewrite of CORSIKA 7, seeing its core functionality re-envisioned in a modern and modular C++ framework. CORSIKA 8 has now reached a "physics-complete" state, offering a robust platform that encourages expert development for specialized applications. It supports high-energy hadronic interactions using models such as Sibyll 2.3d, QGSJet-II.04, EPOS-LHC, and Pythia 8.3, alongside the treatment of electromagnetic cascades with PROPOSAL 7.6.2. Key highlights are the support for multiple interaction media, including cross-media particle showers, and an enhanced calculation of radio emissions from particle showers. This contribution provides an overview of the current functionalities, showcases validation results of its simulations, and discusses future development plans.

The glitch signatures in $r_{010}$ for F-type stars (higher amplitude and period of the oscillatory component) are very different from those of G-type stars. The aim of this work is to analyse the signatures of these glitches and understand the origin of the differences in these signatures between G-type and F-type stars. We fit the glitch signatures in the frequencies, second differences, and $r_{010}$ ratios while assuming either a sinusoidal variation or a more complex expression. The fit provides the acoustic depth, and hence the position, of the bottom of the convective envelope for nine \textit{Kepler} stars and the Sun. We find that for F-type stars, the most commonly used fitting expressions for the glitch of the bottom of the convective envelope provide different measurements of the position of the bottom of the convective envelope for the three seismic indicators, while it is not the case for G-type stars. When adding an additional term in the fitting expression with twice the acoustic depth of the standard term (a contribution that accounts for the highly non-sinusoidal shape of the signature in the $r_{010}$ ratios), we find better agreement between the three seismic indicators and with the prediction of stellar evolution models. While the origin of this additional term is not yet understood, this may be an indication that the transition between the convective envelope and the underlying radiative zone is different for G- and F-type stars. This outcome brings new insight into the physics in these regions.

Model misspecification analysis strategies, such as anomaly detection, model validation, and model comparison are a key component of scientific model development. Over the last few years, there has been a rapid rise in the use of simulation-based inference (SBI) techniques for Bayesian parameter estimation, applied to increasingly complex forward models. To move towards fully simulation-based analysis pipelines, however, there is an urgent need for a comprehensive simulation-based framework for model misspecification analysis. In this work, we provide a solid and flexible foundation for a wide range of model discrepancy analysis tasks, using distortion-driven model misspecification tests. From a theoretical perspective, we introduce the statistical framework built around performing many hypothesis tests for distortions of the simulation model. We also make explicit analytic connections to classical techniques: anomaly detection, model validation, and goodness-of-fit residual analysis. Furthermore, we introduce an efficient self-calibrating training algorithm that is useful for practitioners. We demonstrate the performance of the framework in multiple scenarios, making the connection to classical results where they are valid. Finally, we show how to conduct such a distortion-driven model misspecification test for real gravitational wave data, specifically on the event GW150914.

We present a systematic study of the $BVRI$ colours of novae over the course of their eruptions. Where possible, interstellar reddening was measured using the equivalent widths of Diffuse Interstellar Bands (DIBs). Some novae lack spectra with sufficient resolution and signal-to-noise ratios; therefore, we supplement as necessary with 3D and 2D dust maps. Utilising only novae with DIB- or 3D-map-based $E(B-V)$, we find an average intrinsic $(B-V)_0$ colour of novae at $V$-band light curve peak of 0.18 with a standard deviation of 0.31, based on a sample of 23 novae. When the light curve has declined by 2 magnitudes ($t_2$), we find an average $(B-V)_0 = -0.02$ with a standard deviation of 0.19. These average colours are consistent with previous findings, although the spreads are larger than previously found due to more accurate reddening estimates. We also examined the intrinsic $(R-I)_0$ and $(V-R)_0$ colours across our sample. These colours behave similarly to $(B-V)_0$, except that the $(V-R)_0$ colour gets redder after peak, likely due to the contributions of emission line flux. We searched for correlations between nova colours and $t_2$, peak $V$-band absolute magnitude, and GeV $\gamma$-ray luminosity, but find no statistically significant correlations. Nova colours can therefore be used as standard "crayons" to estimate interstellar reddening from photometry alone, with 0.2--0.3 mag uncertainty. We present a novel Bayesian strategy for estimating distances to Galactic novae based on these $E(B-V)$ measurements, independent of assumptions about luminosity, built using 3D dust maps and a stellar mass model of the Milky Way.

The standard cosmological model currently in force, aka $\Lambda$CDM, has been plagued with a variety of tensions in the last decade or so, which puts it against the wall. At the core of the $\Lambda$CDM we have a rigid cosmological term, $\Lambda$, for the entire cosmic history. Recently, the results from the DESI collaboration suggested the possibility that dark energy (DE) should be dynamical rather than just a cosmological constant. Using a generic $w_0w_a$CDM parameterization, DESI favors quintessence behavior at $\sim 3\sigma$ c.l. However, to alleviate the tensions the DE needs more features. In the proposed $w$XCDM model of [43] the DE is actually a composite cosmic fluid with two components $(X,Y)$ acting sequentially: first $X$ (above a transition redshift $z_t$) and second $Y$ (below $z_t$). Fitting the model to the data, we find that the late component $Y$ behaves as quintessence, like DESI. However, to cure the $H_0$ and growth tensions, $X$ must behave as `phantom matter' (PM), which in contrast to phantom DE provides positive pressure at the expense of negative energy density. Using the SNIa (considering separately Pantheon$+$ and DESY5), cosmic chronometers, transversal BAO, LSS data, and the full CMB likelihood from Planck 2018, we find that both tensions can be completely cut down. We also compare the $w$XCDM with our own results using the standard $w$CDM and $w_0w_a$CDM parameterizations of the DE. In all cases, model $w$XCDM performs much better. Finally, we have repeated our analysis with BAO 3D data (replacing BAO 2D), and we still find that the dynamical DE models (including composite ones) provide a much better fit quality compared to $\Lambda$CDM. The growth tension is alleviated again, but in contrast, the $H_0$-tension remains significant, which is most likely reminiscent of the internal conflict in the BAO sector.

We present a new, cosmologically model-independent, statistical analysis of the Pantheon+ type Ia supernovae spectroscopic dataset, improving a standard methodology adopted by Lane et al. We use the Tripp equation for supernova standardisation alone, thereby avoiding any potential correlation in the stretch and colour distributions. We compare the standard homogeneous cosmological model, i.e., $\Lambda$CDM, and the timescape cosmology which invokes backreaction of inhomogeneities. Timescape, while statistically homogeneous and isotropic, departs from average Friedmann-Lema\^{\i}tre-Robertson-Walker evolution, and replaces dark energy by kinetic gravitational energy and its gradients, in explaining independent cosmological observations. When considering the entire Pantheon+ sample, we find very strong evidence ($\ln B> 5$) in favour of timescape over $\Lambda$CDM. Furthermore, even restricting the sample to redshifts beyond any conventional scale of statistical homogeneity, $z > 0.075$, timescape is preferred over $\Lambda$CDM with $\ln B> 1$. These results provide evidence for a need to revisit the foundations of theoretical and observational cosmology.

Primordial gravitational waves propagate almost unimpeded from the moment they are generated to the present epoch. Nevertheless, they are subject to convolution with a non-trivial transfer function. Within the standard thermal history, shifts in the temperature-redshift relation combine with damping effects by free streaming neutrinos to non-trivially process different wavelengths during radiation domination, with subsequently negligible effects at later times. Presuming a nearly scale invariant primordial spectrum, one obtains a characteristic late time spectrum, deviations from which would indicate departures from the standard thermal history. Given the paucity of probes of the early universe physics before nucleosynthesis, it is useful to classify how deviations from the standard thermal history of the early universe can be constrained from observations of the late time stochastic background. The late time spectral density has a plateau at high frequencies that can in principle be significantly enhanced or suppressed relative to the standard thermal history depending on the equation of state of the epoch intervening reheating and the terminal phase of radiation domination, imprinting additional features from bursts of entropy production, and additional damping at intermediate scales via anisotropic stress production. In this paper, we survey phenomenologically motivated scenarios of early matter domination, kination, and late time decaying particles as representative non-standard thermal histories, elaborate on their late time stochastic background, and discuss constraints on different model scenarios.

We present the discovery of 14 new (and recovery of 4 known) low accretion rate magnetic white dwarfs in post-common envelope binaries that emit strong cyclotron emission using the Zwicky Transient Facility (ZTF) light curves, doubling the known sample size. In addition, we discovered a candidate magnetic period bouncer and recovered three known ones. We confirmed the presence of cyclotron emission using low-resolution spectra in 19 objects. Using the ZTF light curves, follow-up spectra, and the spectral energy distribution, we measured the orbital period, magnetic field strength, and white dwarf temperature of each system. Although the phase-folded light curves have diverse shapes and show a much larger variability amplitude, we show that their intrinsic properties (e.g. period distribution, magnetic field strength) are similar to those of previously known systems. The diversity in light curve shapes can be explained by differences in the optical depth of the accretion spot and geometric differences, the inclination angle and the magnetic spot latitude. The evolutionary states of the longer period binaries are somewhat uncertain but are vary; we found systems consistent with being pre-polars, detached polars, or low-state polars. In addition, we discovered two new low-state polars that likely have brown dwarf companions and could be magnetic period bouncers.

We present improvements to the pointing accuracy of the South Pole Telescope (SPT) using machine learning. The ability of the SPT to point accurately at the sky is limited by its structural imperfections, which are impacted by the extreme weather at the South Pole. Pointing accuracy is particularly important during SPT participation in observing campaigns with the Event Horizon Telescope (EHT), which requires stricter accuracy than typical observations with the SPT. We compile a training dataset of historical observations of astronomical sources made with the SPT-3G and EHT receivers on the SPT. We train two XGBoost models to learn a mapping from current weather conditions to two telescope drive control arguments -- one which corrects for errors in azimuth and the other for errors in elevation. Our trained models achieve root mean squared errors on withheld test data of $2.14''$ in cross-elevation and $3.57''$ in elevation, well below our goal of $5''$ along each axis. We deploy our models on the telescope control system and perform further in situ test observations during the EHT observing campaign in 2024 April. Our models result in significantly improved pointing accuracy: for sources within the range of input variables where the models are best trained, average combined pointing error improved 33%, from $15.9''$ to $10.6''$. These improvements, while significant, fall shy of our ultimate goal, but they serve as a proof of concept for the development of future models. Planned upgrades to the EHT receiver on the SPT will necessitate even stricter pointing accuracy which will be achievable with our methods.

Ultralight axions (ULAs) with masses $10^{-33} \lesssim m/{\rm eV} \lesssim 10^{-12}$ are well motivated in string-inspired models and can be part or all of the dark energy or the dark matter in this range. Since the ULA field oscillates at a frequency $m$ that can be much larger than the expansion rate $H$, accurate and efficient calculation of cosmological observables requires an effective time averaged treatment. While these are well established for $m\gg 10 H_{\rm eq}$, the Hubble rate at matter radiation equality, here we extend and develop these techniques to cover the mass range $10^{-33} \lesssim m/{\rm eV} \lesssim 10^{-18}$. We implement this technique in a full cosmological Boltzmann code ($\mathrm{AxiECAMB}$) with numerical precision sufficiently accurate for current and next-generation cosmic microwave background as well as large-scale structure data analysis. New effects including the time averaging of metric perturbations and hydrostatic equilibrium of the effective fluid result in many orders of magnitude improvements for power spectra accuracy over some previous treatments such as $\mathrm{AxionCAMB}$ in some extreme regions of parameter space and order unity changes of the ULA effects near $\Lambda$CDM models. These improvements may impact the model parameters that resolve various tensions in $\Lambda$CDM at a comparable level.

Here, we consider a classically scale-invariant extension of the Standard Model (SM) with two-component dark matter (DM) candidates, including a Dirac spinor and a scalar DM. We probe the parameter space of the model, constrained by relic density and direct detection, and investigate the generation of gravitational waves (GWs) produced by an electroweak first-order phase transition. The analysis demonstrates that there are points in the parameter space, leading to a detectable GW spectrum arising from the first-order phase transition, which is also consistent with the DM relic abundance and direct detection bounds. These GWs could be observed by forthcoming space-based interferometers such as the Big Bang Observer, Decihertz Interferometer Gravitational-wave Observatory, and Ultimate-Decihertz Interferometer Gravitational-wave Observatory.

This paper describes the different steps to include the adiabatic tidal effects to the gravitational waveform amplitude for quasi-circular non-spinning compact binaries up to the second-and-a-half post-Newtonian (PN) order. The amplitude, that relates the two gravitational wave polarizations, is decomposed onto the basis of spin-weighted spherical harmonics of spin -2, parametrized by the two numbers $(\ell,m)$, where the modes of the waveform correspond to the coefficients of the decomposition. These modes are readily computed from the radiative multipole moments. They can be expressed in a PN-expanded form as well as in a factorized form, suitable to be directly included in effective-one-body models to describe more accurately the waveform of binary neutron stars. We also provide the energy flux and phasing evolution in time and frequency domain. The results presented in this article are collected in an ancillary file.

Cosmologically stable, light particles that came into thermal contact with the Standard Model in the early universe may persist today as a form of hot dark matter. For relics with masses in the eV range, their role in structure formation depends critically on their mass. We trace the evolution of such hot relics and derive their density profiles around cold dark matter halos, introducing a framework for their indirect detection. Applying this framework to axions -- a natural candidate for a particle that can reach thermal equilibrium with the Standard Model in the early universe and capable of decaying into two photons -- we establish stringent limits on the axion-photon coupling $g_{a \gamma} $ using current observations of dwarf galaxies, the Milky Way halo, and galaxy clusters. Our results set new bounds on hot axions in the $\mathcal{O}(1-10)\,$eV range.

Atom interferometers offer exceptional sensitivity to ultra-light dark matter (ULDM) through their precise measurement of phenomena acting on atoms. While previous work has established their capability to detect scalar and vector ULDM, their potential for detecting spin-2 ULDM remains unexplored. This work investigates the sensitivity of atom interferometers to spin-2 ULDM by considering several frameworks for massive gravity: a Lorentz-invariant Fierz-Pauli case and two Lorentz-violating scenarios. We find that coherent oscillations of the spin-2 ULDM field induce a measurable phase shift through three distinct channels: coupling of the scalar mode to atomic energy levels, and vector and tensor effects that modify the propagation of atoms and light. Atom interferometers uniquely probe all of these effects, while providing sensitivity to a different mass range from laser interferometers. Our results demonstrate the potential of atom interferometers to advance the search for spin-2 dark matter through accessing unexplored parameter space and uncovering new interactions between ULDM and atoms.

Dark matter captured in stars can act as an additional heat transport mechanism, modifying fusion rates and asteroseismoloigcal observables. Calculations of heat transport rates rely on approximate solutions to the Boltzmann equation, which have never been verified in realistic stars. Here, we simulate heat transport in the Sun, the Earth, and a brown dwarf model, using realistic radial temperature, density, composition and gravitational potential profiles. We show that the formalism developed in arXiv:2111.06895 remains accurate across all celestial objects considered, across a wide range of kinematic regimes, for both spin-dependent and spin-independent interactions where scattering with multiple species becomes important. We further investigate evaporation rates of dark matter from the Sun, finding that previous calculations appear robust. Our Monte Carlo simulation software Cosmion is publicly available.

We consider the gravitational collapse of a homogeneous pressureless ellipsoid. We have shown that the minimal size $r$ that the ellipsoid can reach during collapse depends on its initial eccentricity $e_0$ as $r\propto e_0^\nu$, where $\nu \approx 15/8$, and this dependence is very universal. We have estimated the parameters (in particular, the initial eccentricity) of a homogeneous pressureless ellipsoid, whereat it collapses directly into a black hole.

This work shows an anomalously enhanced response of the low-latitude ionosphere over the Indian sector under weak geomagnetic conditions (October 31, 2021) in comparison to a stronger event (November 04, 2021) under the influence of an Interplanetary Coronal Mass Ejection (ICME)-driven Magnetic Cloud (MC)-like and sheath regions respectively. The investigation is based on measurements of the Total Electron Content (TEC) from Ahmedabad (23.06$^\circ$N, 72.54$^\circ$E, geographic; dip angle: 35.20$^\circ$), a location near the northern crest of the Equatorial Ionization Anomaly (EIA) over the Indian region. During the weaker event, the observed TEC from the Geostationary Earth Orbit (GEO) satellites of Navigation with Indian Constellation (NavIC), showed diurnal maximum enhancements of about 20 TECU over quiet-time variations, as compared to the stronger event where no such enhancements are present. It is shown that storm intensity (SYM-H) or magnitude of the southward Interplanetary Magnetic Field (IMF) alone is unable to determine the ionospheric impacts of this space weather event. However, it is the non-fluctuating southward IMF and the corresponding penetration electric fields, for a sufficient interval of time, in tandem with the poleward neutral wind variations, that determines the strengthening of low-latitude electrodynamics of this anomalous event of October 31, 2021. Therefore, the present investigation highlights a case for further investigations of the important roles played by non-fluctuating penetration electric fields in determining a higher response of the low-latitude ionosphere even if the geomagnetic storm intensities are significantly low.

Neutrinos being massive could undergo non-radiative decay, a property for which the diffuse supernova neutrino background has a unique sensitivity. We extend previous analyses to explore our ability to disentangle predictions for the diffuse supernova neutrino background in presence or absence of neutrino non-radiative two-body decay. In a three-neutrino framework, we give predictions of the corresponding neutrino fluxes and the expected number of events in the Super-Kamiokande+Gadolinium, the Hyper-Kamiokande, the JUNO and the DUNE experiments. In our analysis, we employ supernova simulations from different groups and include current uncertainties from both the evolving core-collapse supernova rate and the fraction of failed supernovae. We perform the first Bayesian analysis to see our ability to disentangle the cases in presence and absence of neutrino decay. To this aim we combine the expected events in inverse beta-decay and the neutrino-argon detection channels. We also discuss neutrino-electron, neutrino-proton and of neutrino-oxygen scattering. Our investigation covers the different possible decay patterns for normal mass ordering, both strongly-hierarchical and quasi-degenerate as well as the inverted neutrino mass ordering.

We propose a revised formulation of General Relativity for cosmological contexts, in which the Einstein constant varies with the energy density of the Universe. We demonstrate that this modification has no direct phenomenological impact on the Universe's dynamics or on particle motion within the expanding cosmos. Assuming a state close to vacuum, here defined by the vanishing product of the Einstein coupling constant and the Universe's energy density, we perform a Taylor expansion of the theory. In this framework, the vacuum energy problem is addressed, and an additional constant pressure term, which induces a Chaplygin-like contribution to the dark energy equation of state, arises in the late-time dynamics. The correction to the late-time Hubble parameter is investigated by comparing theoretical predictions with the late Universe observational data. Our findings indicate that the current value of the vacuum energy is consistent with zero. Alternatively, the expansion used in our formulation would no longer be valid if the current state significantly deviates from the assumed near-vacuum condition. Implications of the modified $\Lambda$CDM model with respect to the Hubble tension are also discussed.

Living organisms have some common structures, chemical reactions and molecular structures. The organisms consist of cells with cell division, they have homochirality of protein and carbohydrate units, and metabolism, and genetics, and they are mortal. The molecular structures and chemical reactions underlying these features are common from the simplest bacteria to human beings. The origin of life is evolutionary with the emergence of a network of spontaneous biochemical reactions, and the evolution has taken place over a very long time. The evolution contains, however some "landmarks" and bottlenecks, which in a revolutionary manner directed the evolution, and the article tries to establish the order of these events. The article advocates that a possible order in the emergence of life is that the first milestone in prebiotic evolution is at the emergence of homochirality in proteins. The homochirality of peptides is, however, with instability and racemization which causes aging of the peptides and mortality. The metabolism and genetics are established through homochiral enzymes in the Earth's crust for $\approx$ 4 Gyr ago. Finally, the cells with cell division are established in the Hot Springs environment at the interface between the crust and the Hadean Ocean.

We derive the model of the Schwarzschild black hole immersed into a dark matter halo with a relativistic Hernquist profile, the Schwarzschild-Hernquist black hole, and obtain its tidal Love numbers and quasi-normal modes. We thoroughly compare our odd and even parity perturbation equations with the literature and point out that two distinct choices of matter perturbations lead to distinct spectra. We establish that the quasi-normal modes admit qualitatively distinct scaling laws in terms of dark matter densities for non-relativistic and relativistic halos. We develop a stable numerical scheme for computing tidal Love numbers based on asymptotic series expansions. We further comment upon the existence of matter configurations obeying the dominant energy condition that lead to multiple light rings.

This summary of the second Terrestrial Very-Long-Baseline Atom Interferometry (TVLBAI) Workshop provides a comprehensive overview of our meeting held in London in April 2024, building on the initial discussions during the inaugural workshop held at CERN in March 2023. Like the summary of the first workshop, this document records a critical milestone for the international atom interferometry community. It documents our concerted efforts to evaluate progress, address emerging challenges, and refine strategic directions for future large-scale atom interferometry projects. Our commitment to collaboration is manifested by the integration of diverse expertise and the coordination of international resources, all aimed at advancing the frontiers of atom interferometry physics and technology, as set out in a Memorandum of Understanding signed by over 50 institutions.

This article provides an overview of the current state of machine learning in gravitational-wave research with interferometric detectors. Such applications are often still in their early days, but have reached sufficient popularity to warrant an assessment of their impact across various domains, including detector studies, noise and signal simulations, and the detection and interpretation of astrophysical signals. In detector studies, machine learning could be useful to optimize instruments like LIGO, Virgo, KAGRA, and future detectors. Algorithms could predict and help in mitigating environmental disturbances in real time, ensuring detectors operate at peak performance. Furthermore, machine-learning tools for characterizing and cleaning data after it is taken have already become crucial tools for achieving the best sensitivity of the LIGO--Virgo--KAGRA network. In data analysis, machine learning has already been applied as an alternative to traditional methods for signal detection, source localization, noise reduction, and parameter estimation. For some signal types, it can already yield improved efficiency and robustness, though in many other areas traditional methods remain dominant. As the field evolves, the role of machine learning in advancing gravitational-wave research is expected to become increasingly prominent. This report highlights recent advancements, challenges, and perspectives for the current detector generation, with a brief outlook to the next generation of gravitational-wave detectors.

We investigate the gravitational wave emission from extreme mass ratio inspirals, key targets for the upcoming space-based detector LISA, considering the scenario where the lighter black hole in the binary is endowed by a long-lived, time-dependent scalar field configuration, known as a scalar wig. We develop a formalism to compute scalar perturbations for extreme mass ratio inspirals on circular orbits around Schwarzschild and Kerr black holes, and apply this framework to compute additional fluxes induced by the scalar wig as well as their dependence on the scalar field properties. Our computation provides strong indications that in this scenario the presence of the scalar field does not significantly affect the orbital motion and the gravitational waveform.

We present a novel mechanism for the irreducible production of magnetic monopoles from interactions of cosmic rays and interstellar medium (ISM). Resulting monopoles drain energy from galactic magnetic fields, disrupting their formation and sustainability. We generalize conventional Parker bounds to monopoles with extended energy spectrum and, considering cosmic ray ISM monopole production, set novel constraints from disruption of Milky Way Galactic magnetic fields and their seeds. Further, we set first constraints on disruption of galactic magnetic fields and their seeds of Andromeda galaxy, with results being competitive with distinct existing bounds. Unlike Parker limits of previous works that relied on cosmological monopoles, our constraints are independent of cosmological monopole production or their primordial abundance. Besides, we estimate new constraints on dipole magnetic moments generated from cosmic ray ISM interactions. We discuss implications for monopoles with generalized magnetic charges.

Based on the Newtonian mechanics, in this article, we present a heuristic derivation of the Friedmann equations, providing an intuitive foundation for these fundamental relations in cosmology. Additionally, using the first law of thermodynamics and Euler's equation, we derive a set of equations that, at linear order, coincide with those obtained from the conservation of the stress-energy tensor in General Relativity. This approach not only highlights the consistency between Newtonian and relativistic frameworks in certain limits but also serves as a pedagogical bridge, offering insights into the physical principles underlying the dynamics of the universe.

For gravitationally lensed type II signals, the phase of the dominant (2, 2) mode and the higher order (3, 3) mode is offset by $-\pi/12$. Using this, we develop a test for type II imagery by allowing the phases of the (2,2) and (3,3) modes to vary separately and introducing a new waveform parameter to represent the phase offset between the two. We use simulated, low mass ratio, precessing signals to show that the test can reproduce the $-\pi/12$ phase offset when detected by three detectors for H-L optimal SNR $\gtrsim$ 40 and $\mathcal{M} \leq 30$. We analyze GW190412 and GW190814 using this parameterization, measuring the offset to be $0.13^{+0.22}_{-0.17}$ for GW190412 and $-0.05^{+0.20}_{-0.22}$ for GW190814. We also measure the Bayes Factor in support of no phase offset, $\log_{10} \mathcal{B}_{\Delta \varphi = 0}$, to be $-0.14$ for GW190412 and $0.21$ for GW190814.

We investigate the QCD phase transition and its phase structure within Einstein-Maxwell-Dilaton-scalar system and compare the results with those obtained from the Einstein-Maxwell-Dilaton system. It is shown that both models reproduce behavior consistent with lattice QCD. In particular, the Einstein-Maxwell-Dilaton-scalar system exhibits a first-order phase transition in the pure gauge sector, aligning with predictions from Yang-Mills theory. Based on these models, we construct a holographic model for neutron stars, incorporating leptons to satisfy electric charge neutrality, and examine the cold equation of state, the mass-radius relation, and tidal deformability of neutron stars. It is demonstrated that the Einstein-Maxwell-Dilaton-scalar system enables us to describe neutron star properties that meet current astrophysical constraints.

This paper calculates the stochastic gravitational wave background from dark binaries with finite-range attractive dark forces, complementing previous works which consider long-range dark forces. The finiteness of the dark force range can dramatically modify both the initial distributions and evolution histories of the binaries. The generated gravitational wave spectrum is enhanced in the intermediate frequency regime and exhibits interesting "knee" and "ankle" features, the most common of which is related to the turn on of the dark force mediator radiation. Other such spectral features are related to changes in the binary merger lifetime and the probability distribution for the initial binary separation. The stochastic gravitational wave background from sub-solar-mass dark binaries is detectable by both space- and ground-based gravitational wave observatories.