Papers for Thursday, Jun 29 2023

Distance duality relation (DDR) marks a fundamental difference between expanding and nonexpanding Universes, as an expanding metric causes angular diameter distance smaller than luminosity distance by an extra factor of $(1+z)$. Here we report a test of this relation using two independent samples of ultracompact radio sources observed at 2.29 GHz and 5.0 GHz. The test with radio sources involves only geometry, so it is independent of cosmological models. Since the observed radio luminosities systematically increase with redshift, we do not assume a constant source size. Instead, we start with assuming the intensive property, luminosity density, does not evolve with redshift and then infer its evolution from the resultant DDR. We make the same assumption for both samples, and find it results in the same angular size-redshift relation. Interestingly, the resultant DDR is fully consistent with a nonexpanding Universe. Imposing the DDR predicted by the expanding Universe, we infer the radio luminosity density evolves as $\rho_L\propto(1+z)^3$. However, the perfect agreement with a nonexpanding Universe under the assumption of constant luminosity densities poses a conspiracy and fine-tuning problem: the size and luminosity density of ultracompact radio sources evolve in the way that precisely mimics a nonexpanding Universe.

Astrophysical explorations are underpinned by large-scale stellar spectroscopy surveys, necessitating a paradigm shift in spectral fitting techniques. Our study proposes three enhancements to transcend the limitations of the current spectral emulation models. We implement an attention-based emulator, adept at unveiling long-range information between wavelength pixels. We leverage a domain-specific fine-tuning strategy where the model is pre-trained on spectra with fixed stellar parameters and variable elemental abundances, followed by fine-tuning on the entire domain. Moreover, by treating wavelength as an autonomous model parameter, akin to neural radiance fields, the model can generate spectra on any wavelength grid. In the case with a training set of O(1000), our approach exceeds current leading methods by a factor of 5-10 across all metrics.

The existence of a cosmic neutrino background has been inferred indirectly from cosmological surveys through its effect on the linear-theory evolution of primordial density perturbations, as well as from measurements of the primordial abundances of light elements. Constraints on the masses of the three neutrino species imply that at least two of them move non-relativistically today. As a consequence, non-linear evolution of density perturbations results in the formation of neutrino halos around dark-matter halos. We study whether these neutrino halos can be detected in the foreseeable future through measurements of weak gravitational lensing of the cosmic microwave background, thus providing, possibly, the first beyond-linear-theory signature of cosmic neutrinos.

The Magellanic Stream (MS) - an enormous ribbon of gas spanning $140^\circ$ of the southern sky trailing the Magellanic Clouds - has been exquisitely mapped in the five decades since its discovery. However, despite concerted efforts, no stellar counterpart to the MS has been conclusively identified. This stellar stream would reveal the distance and 6D kinematics of the MS, constraining its formation and the past orbital history of the Clouds. We have been conducting a spectroscopic survey of the most distant and luminous red giant stars in the Galactic outskirts. From this dataset, we have discovered a prominent population of 13 stars matching the extreme angular momentum of the Clouds, spanning up to $100^\circ$ along the MS at distances of $60-120$ kpc. Furthermore, these kinemetically-selected stars lie along a [$\alpha$/Fe]-deficient track in chemical space from $-2.5 < \mathrm{[Fe/H]} < -0.5$, consistent with their formation in the Clouds themselves. We identify these stars as high-confidence members of the Magellanic Stellar Stream. Half of these stars are metal-rich and closely follow the gaseous MS, whereas the other half are more scattered and metal-poor. We argue that the metal-rich stream is the recently-formed tidal counterpart to the MS, and speculate that the metal-poor population was thrown out of the SMC outskirts during an earlier interaction between the Clouds. The Magellanic Stellar Stream provides a strong set of constraints - distances, 6D kinematics, and birth locations - that will guide future simulations towards unveiling the detailed history of the Clouds.

We propose a machine-learning-based technique to determine the number density of radio sources as a function of their flux density, for use in next-generation radio surveys. The method uses a convolutional neural network trained on simulations of the radio sky to predict the number of sources in several flux bins. To train the network, we adopt a supervised approach wherein we simulate training data stemming from a large domain of possible number count models going down to fluxes a factor of 100 below the threshold for source detection. We test the model reconstruction capabilities as well as benchmark the expected uncertainties in the model predictions, observing good performance for fluxes down to a factor of ten below the threshold. This work demonstrates that the capabilities of simple deep learning models for radio astronomy can be useful tools for future surveys.

The merger of neutron stars drives a relativistic jet which can be observed as a short gamma-ray burst. A strong large-scale magnetic field is necessary to launch the relativistic jet. However, the magnetohydrodynamical mechanism to build up this magnetic field remains uncertain. Here we show that the $\alpha\Omega$ dynamo mechanism driven by the magnetorotational instability builds up the large-scale magnetic field inside the long-lived binary neutron star merger remnant by performing an {\it ab initio} super-high resolution neutrino-radiation magnetohydrodynamics merger simulation in full general relativity. As a result, the magnetic field induces the Poynting-flux dominated relativistic outflow with the luminosity $\sim 10^{51}$\,erg/s and magnetically-driven post-merger mass ejection with the mass $\sim 0.1M_\odot$. Therefore, the magnetar scenario in binary neutron star mergers is possible. These can be the engines of short-hard gamma-ray bursts and very bright kilonovae. Therefore, this scenario is testable in future observation.

We introduce a novel method for gravitational wave detection that combines: 1) hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational wave signals; and 2) graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a detector network. These spatiotemporal-graph AI models are tested for signal detection of gravitational waves emitted by quasi-circular, non-spinning and quasi-circular, spinning, non-precessing binary black hole mergers. For the latter case, we needed a dataset of 1.2 million modeled waveforms to densely sample this signal manifold. Thus, we reduced time-to-solution by training several AI models in the Polaris supercomputer at the Argonne Leadership Supercomputing Facility within 1.7 hours by distributing the training over 256 NVIDIA A100 GPUs, achieving optimal classification performance. This approach also exhibits strong scaling up to 512 NVIDIA A100 GPUs. We then created ensembles of AI models to process data from a three detector network, namely, the advanced LIGO Hanford and Livingston detectors, and the advanced Virgo detector. An ensemble of 2 AI models achieves state-of-the-art performance for signal detection, and reports seven misclassifications per decade of searched data, whereas an ensemble of 4 AI models achieves optimal performance for signal detection with two misclassifications for every decade of searched data. Finally, when we distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, our AI ensemble is capable of processing a decade of gravitational wave data from a three detector network within 3.5 hours.

Protoclusters are the densest regions in the distant universe ($z>2$) and are the progenitors of massive galaxy clusters ($M_{halo}>10^{14}{\rm M}_\odot$) in the local universe. They undoubtedly play a key role in early massive galaxy evolution and they may host the earliest sites of galaxy quenching or even induce extreme states of star formation. Studying protoclusters therefore not only gives us a window into distant galaxy formation but also provides an important link in our understanding of how dense structures grow over time and modify the galaxies within them. Current protocluster samples are completely unable to address these points because they are small and selected in a heterogeneous way. We propose the Roman-Cosmic Noon survey, whose centerpiece is an extremely deep (30ksec) and wide area (10 deg$^2$) prism slitless spectroscopy survey to identify the full range of galaxy structures at $2<z<3$. This survey will include 1500 uniformly selected protoclusters, their surrounding cosmic web environments, and at least 15,000 protocluster galaxies with $M_\star>10^{10.5} {\rm M}_\odot$ across the full range of star formation histories as well as many more lower mass star-forming galaxies. The survey will also contain field galaxies to much lower masses than in the High Latitude Wide Area Survey, but over an area dwarfing any current or planned deep spectroscopy probe at $z>2$. With the prism spectroscopy and some modest additional imaging this survey will measure precise stellar mass functions, quenched fractions, galaxy and protocluster morphologies, stellar ages, emission-line based SFRs, and metallicities. It will have extensive legacy value well beyond the key protocluster science goals.

We present an analysis of major-axis surface-brightness profiles of bars in a volume-limited sample of 182 barred spiral galaxies, using Spitzer 3.6 micron images. Unlike most previous studies, we use the entire bar profile, and we classify profiles into four categories. These are "Peak+Shoulders" (P+Sh) -- updating the classic "flat bar" profile -- and three subtypes of the classic "exponential" profile: (true) Exponential, "Two-Slope" (shallow inner slope + steeper outer slope), and "Flat-Top" (constant inner region, steep outer slope). P+Sh profiles are preferentially found in galaxies with high stellar masses, early Hubble types, red colours, and low gas fractions; the most significant factor is stellar mass, and previous correlations with Hubble type can be explained by the tendency of higher-mass galaxies to have earlier Hubble types. The most common type of non-P+Sh profile is Exponential, followed by Flat-Top profiles; all non-P+Sh profiles appear to have similar distributions of stellar mass, Hubble type, colour, and gas fraction. We also morphologically classify the bars of an inclined subsample into those with and without boxy/peanut-shaped (B/P) bulges; as previously reported, the presence of a B/P bulge is very strong function of stellar mass. Essentially all bars with B/P bulges have P+Sh profiles; we associate the profile shoulders with the outer, vertically thin part of the bar. We find a small number of P+Sh profiles in bars without clear B/P bulges, which may indicate that P+Sh formation precedes the formation of B/P bulges.

We investigate the role of magnetic field on the gas dynamics in the Galactic bulge region by three dimensional simulations with radiative cooling and heating. While high-temperature corona with $T>10^6\ {\rm K}$ is formed in the halo regions, the temperature near the Galactic plane is $\lesssim 10^4\ {\rm K}$ following the thermal equilibrium curve determined by the radiative cooling and heating. Although the thermal energy of the interstellar gas is lost by radiative cooling, the saturation level of the magnetic field strength does not significantly depend on the radiative cooling and heating. The magnetic field strength is amplified to $10\ {\rm \mu G}$ on average, and reaches several hundred ${\rm \mu G}$ locally. We find the formation of magnetically dominated regions at mid-latitudes in the case with the radiative cooling and heating, which is not seen in the case without radiative effect. The vertical thickness of the mid-latitude regions is $50-150\ {\rm pc}$ at the radial location of $0.4-0.8 \ {\rm kpc}$ from the Galactic center, which is comparable to the observed vertical distribution of neutral atomic gas. When we take the average of different components of energy density integrated over the Galactic bulge region, the magnetic energy is comparable to the thermal energy. We conclude that the magnetic field plays a substantial role in controlling the dynamical and thermal properties of the Galactic bulge region.

The semi-regular variable star V Canum Venaticorum (V CVn) is well-known for its unusual linear polarization position angle (PA). Decades of observing V CVn reveal a nearly constant PA spanning hundreds of pulsation cycles. This phenomenon has persisted through variability that has ranged by 2 magnitudes in optical brightness and through variability in the polarization amplitude over 0.3% and 6.9%. Additionally, the polarization fraction of V CVn varies inversely with brightness. This paper presents polarization measurements obtained over three pulsation cycles. We find that the polarization maximum does not always occur precisely at the same time as the brightness minimum. Instead, we observe a small lead or lag in relation to the brightness minimum, spanning a period of a few days up to three weeks. Furthermore, the PA sometimes exhibits a non-negligible rotation, especially at lower polarization levels. To elucidate the unusual optical behavior of V CVn, we present a list of literature sources that also exhibit polarization variability with a roughly fixed PA. We find this correlation occurs in stars with high tangential space velocities, i.e., "runaway" stars, suggesting that the long-term constant PA is related to how the circumstellar gas is shaped by the star's high-speed motion through the interstellar medium.

We aim at extracting the transmission spectrum of the HI Balmer lines of the ultra-hot Jupiter (UHJ) KELT-20b/MASCARA-2b from observations and to further compare the results with what obtained through forward modelling accounting for non-local thermodynamic equilibrium (NLTE) effects. We extract the line profiles from six transits obtained with the HARPS-N high-resolution spectrograph attached to the Telescopio Nazionale Galileo telescope. We compute the temperature-pressure (TP) profile employing the helios code in the lower atmosphere and the Cloudy NLTE code in the middle and upper atmosphere. We further use Cloudy to compute the theoretical planetary transmission spectrum in LTE and NLTE for comparison with observations. We detected the Halpha (0.79+/-0.03%; 1.25 Rp), Hbeta (0.52+/-0.03%; 1.17 Rp), and Hgamma (0.39+/-0.06%; 1.13 Rp) lines, while we detected the Hdelta line at almost 4 sigma (0.27+/-0.07%; 1.09 Rp). The models predict an isothermal temperature of about2200 K at pressures >10^-2 bar and of about 7700 K at pressures <10^-8 bar, with a roughly linear temperature rise in between. In the middle and upper atmosphere, the NLTE TP profile is up to about 3000 K hotter than in LTE. The synthetic transmission spectrum derived from the NLTE TP profile is in good agreement with the observed HI Balmer line profiles, validating our obtained atmospheric structure. Instead, the synthetic transmission spectrum derived from the LTE TP profile leads to significantly weaker absorption compared to the observations. Metals appear to be the primary agents leading to the temperature inversion in UHJs and the impact of NLTE effects on them increases the magnitude of the inversion. We find that the impact of NLTE effects on the TP profile of KELT-20b/MASCARA-2b is larger than for the hotter UHJ KELT-9b, and thus NLTE effects might be relevant also for planets cooler than KELT-20b/MASCARA-2b.

Building on the success of the Starobinsky model in describing the inflationary period of the universe, we investigate two simple generalizations of this model and their constraints imposed by the reheating epoch. The first generalization takes the form $R^{2p}$, while the second is the $\alpha$-Starobinsky model. We first focus on the case where $p=1$ or equivalently, $\alpha=1$, which corresponds to the original Starobinsky model. We derive exact consistency relations between observables and cosmological quantities, without neglecting any terms, and impose the reheating condition $0 < \omega_{re} < 0.25$, where $\omega_{re}$ is the equation of state parameter at the end of reheating. This allows us to obtain new bounds for $n_s$ and $r$ that satisfy this condition and apply them to other observables and cosmological quantities. We repeat this process for the cases where $p \neq 1$ and $\alpha \neq 1$ and find that these generalizations only result in minor modifications of the Starobinsky model, including the potential and the bounds on observables and cosmological quantities.

Failed supernovae (FSN) are a possible channel for the formation of heavy stellar-mass black holes ($M_{ BH}>\sim 30$ M$_\odot$). However, the effects of metallicity, rotation and magnetic field on the islands of explodabilty of massive stars are not clear. Here, we simulate the stellar structure and evolution in the mass range between 6 and 55 $M_{\odot}$ with different initial rotational velocities, metallicities, and magnetic fields from zero-age main sequence (ZAMS) to pre-collapse. We find that the rapid rotating stars can remain lower $\rm ^{12}C$ mass fraction at the time of C ignition, which allows the transition, from convective carbon burning to radiative burning, to occur at lower $M_{\rm ZAMS}$ than those from stars without rotation. However, the rapid rotation is unfavorable for FSN occurring but is conducive to long gamma-ray bursts (lGRBs) because it results in the specific angular momentum in the CO core is greater than the last stable orbit at core collapse. The increasing metallicity does not affect FSN islands, but high metallicity inhibits rotational mixing and is unfavorable for producing lGRBs. A magnetic field can constrain the mass-loss rate even for rapid rotating stars, resulting in higher mass at pre-collapse. The magnetic braking triggered by the magnetic field can reduce the rotation velocity for high-metallicity models, which decreases the specific angular momentum in the CO core and is favorable for FSN occurring. We suggest that the heavy-mass black holes detected by LIGO may originate from rapidly rotating massive stars with strong magnetic fields, rather than those with very low metallicity.

We investigate observational signatures of explosive chromospheric evaporation during a white-light flare (WLF) that occurred on 2022 August 27. Using the moment analysis, bisector techniques, and the Gaussian fitting method, red-shifted velocities of less than 20 km/s are detected in low-temperature spectral lines of Ha, C I and Si IV at the conjugated flare kernels, which could be regarded as downflows caused by chromospheric condensation. Blue-shifted velocities of about 30-40 km/s are found in the high-temperature line of Fe XXI, which can be interpreted as upflows driven by chromospheric evaporation. A nonthermal hard X-ray (HXR) source is co-spatial with one of the flare kernels, and the Doppler velocities are temporally correlated with the HXR fluxes. The nonthermal energy flux is estimated to be at least (1.3+-0.2)*10^10 erg/s/cm^2. The radiation enhancement at Fe I 6569.2 A and 6173 A suggests that the flare is a WLF. Moreover, the while-light emission at Fe I 6569.2 A is temporally and spatially correlated with the blue shift of Fe XXI line, suggesting that both the white-light enhancement and the chromospheric evaporation are triggered and driven by nonthermal electrons. All our observations support the scenario of an electron-driven explosive chromospheric evaporation in the WLF.

We quantify the differences between stellar accelerations in disk galaxies formed in a MONDian universe relative to galaxies with the identical baryonic matter distributions and a fitted cold dark matter halo. In a Milky Way-like galaxy the maximal transverse acceleration is ${\cal {O}}(10^{-9})$ arcseconds per year per decade, well beyond even the most optimistic extrapolations of current capabilities. Conversely, the maximum difference in the line-of-sight acceleration is ${\cal {O}}(1)$ centimetre per second per decade at solar distances from the galactic centre. This level of precision is within reach of plausible future instruments.

This paper is concerned with the kinetic treatment of quasi-stationary axisymmetric collisionless accretion disc plasmas. The conditions of validity of the kinetic description for non-relativistic magnetized and gravitationally-bound plasmas of this type are discussed. A classification of the possible collisionless plasma regimes which can arise in these systems is proposed, which can apply to accretion discs around both stellar-mass compact objects and galactic-center black holes. Two different classifications are determined, which are referred to respectively as energy-based and magnetic field-based classifications. Different regimes are pointed out for each plasma species, depending both on the relative magnitudes of kinetic and potential energies and the magnitude of the magnetic field. It is shown that in all cases, there can be quasi-stationary Maxwellian-like solutions of the Vlasov equation. The perturbative approach outlined here permits unique analytical determination of the functional form for the distribution function consistent, in each kinetic regime, with the explicit inclusion of finite Larmor radius-diamagnetic and/or energy-correction effects.

We analyze all the available Atacama Large Millimeter / submillimeter Array archival data of the nearby Type-II Seyfert galaxy NGC 1068, including new 100 GHz data with the angular resolution of 0\farcs05, which was not included in previous continuum spectral analysis. By combining with the literature data based on the Very Large Array, we investigate the broadband radio continuum spectrum of the central $\lesssim7$ pc region of NGC 1068. We found that the flux density is between $\approx$10-20 mJy at 5-700 GHz. Due to the inability of the model in previous studies to account for the newly added 100 GHz data point, we proceeded to update the models and make the necessary adjustments to the parameters. One possible interpretation of this broadband radio spectrum is a combination of emission from the jet base, the dusty torus, and the compact X-raying corona with the magnetic field strength of $\approx20$ G on scales of $\approx30$ Schwarzschild radii from the central black hole. In order to firmly identify the compact corona by omitting any other possible extended components (e.g., free-free emission from ionized gas around), high-resolution/sensitivity observations achieved by next-generation interferometers will be necessary.

Observations suggested that the occurrence rate of hot Jupiters (HJs) in open clusters is largely consistent with the field ($\sim1\%$) but in the binary-rich cluster M67, the rate is $\sim5\%$. How does the cluster environment boost HJ formation via the high-eccentricity tidal migration initiated by the extreme-amplitude von Zeipel-Lidov-Kozai (XZKL) mechanism forced by a companion star? Our analytical treatment shows that the cluster's collective gravitational potential alters the companion's orbit slowly, which may render the star-planet-companion configuration XZKL-favourable, a phenomenon only possible for very wide binaries. We have also performed direct Gyr $N$-body simulations of the star cluster evolution and XZKL of planets' orbit around member stars. We find that an initially-single star may acquire a companion star via stellar scattering and the companion may enable XZKL in the planets' orbit. Planets around an initially-binary star may also be XZKL-activated by the companion. In both scenarios, the companion's orbit has likely been significantly changed by star scattering and the cluster potential before XZKL occurs in the planets' orbits. Across different cluster models, 0.8\%-3\% of the planets orbiting initially-single stars have experienced XZKL while the fraction is 2\%-26\% for initially-binary stars. Notably, the ejection fraction is similar to or appreciably smaller than XZKL. Around a star that is binary at 1 Gyr, 13\%-32\% of its planets have undergone XZKL, and combined with single stars, the overall XZKL fraction is 3\%-21\%, most affected by the cluster binarity. If 10\% of the stars in M67 host a giant planet, our model predicts an HJ occurrence rate of $\sim1\%$. We suggest that HJ surveys target old, high-binarity, not-too-dense open clusters and prioritise wide binaries to maximise HJ yield.

We examine spatial variations of the C$^0$/CO abundance ratio ($X_{\mathrm{C/CO}}$) in the vicinity of the $\gamma$-ray supernova remnant W51C, based on [CI] ($^3P_1$-$^3P_0$), $^{12}$CO(1-0), and $^{13}$CO(1-0) observations with the ASTE and Nobeyama 45-m telescopes. We find that $X_{\mathrm{C/CO}}$ varies in a range of 0.02-0.16 (0.05 in median) inside the molecular clouds of $A_V>$100 mag, where photodissociation of CO by the interstellar UV is negligible. Furthermore, $X_{\mathrm{C/CO}}$ is locally enhanced up to by a factor of four near the W51C center, depending on the projected distance from the W51C center. In high-$A_V$ molecular clouds, $X_{\mathrm{C/CO}}$ is determined by the ratio of the cosmic-ray (CR) ionization rate to the H$_2$ density, and we find no clear spatial variation of the H$_2$ density against the projected distance. Hence, the high CR ionization rate may locally enhance $X_{\mathrm{C/CO}}$ near the W51C center. We also find that the observed spatial extent of the enhanced $X_{\mathrm{C/CO}}$ ($\sim$17 pc) is consistent with the diffusion distance of CRs with the energy of 100 MeV. The fact suggests that the low-energy CRs accelerated in W51C enhance $X_{\mathrm{C/CO}}$. The CR ionization rate at the $X_{\mathrm{C/CO}}$-enhanced cloud is estimated to be 3$\times$10$^{-16}$ s$^{-1}$ on the basis of time-dependent PDR simulations of $X_{\mathrm{C/CO}}$, the value of which is 30 times higher than that in the standard Galactic environment. These results demonstrate that [CI] is a powerful probe to investigate the interaction between CRs and the interstellar medium for a wide area in the vicinity of supernova remnants.

Remote brightenings often appear at outskirts of source active regions of solar eruptive events, nevertheless, their origin remains to be ascertained. In this study, we report imaging and spectroscopic observations of sequential remote brightenings with a combination of H$\alpha$ Imaging Spectrograph (HIS) onboard the Chinese H$\alpha$ Solar Explorer (CHASE), which is the first space-based solar telescope of China, and the Solar Dynamics Observatory. It is found that, during two successive M-class flares occurring on 2022 August 17, multiple ribbon-like brightenings appeared in sequence away from the flaring active region. Meanwhile, abundant cool filament materials moved downward to the sequential remote brightenings as visible at the H$\alpha$ red wing with a line-of-sight speed up to 70 km s$^{-1}$. The extrapolated three-dimensional magnetic field configuration shows that the sequential remote brightenings correspond to the footpoints of closed ambient field lines whose conjugate footpoints are rooted in the main flare site. We suggest that the sequential remote brightenings are most likely caused by the heating of interchange reconnection between the erupting flux rope and closed ambient field, during which the rope-hosting filament materials are transferred to the periphery of flaring active region along the closed ambient field rather than to the interplanetary space like in the scenario of the slow solar wind formation.

We develop extreme data compression for use in Bayesian model comparison via the MOPED algorithm, as well as more general score compression. We find that Bayes factors from data compressed with the MOPED algorithm are identical to those from their uncompressed datasets when the models are linear and the errors Gaussian. In other nonlinear cases, whether nested or not, we find negligible differences in the Bayes factors, and show this explicitly for the Pantheon-SH0ES supernova dataset. We also investigate the sampling properties of the Bayesian Evidence as a frequentist statistic, and find that extreme data compression reduces the sampling variance of the Evidence, but has no impact on the sampling distribution of Bayes factors. Since model comparison can be a very computationally-intensive task, MOPED extreme data compression may present significant advantages in computational time.

We present a multi-wavelength analysis using SMA, JCMT, NOEMA, JWST, HST, and SST of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A1489. These SMA-located sources both lie at z=4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at <2um, even with deep JWST/NIRCam observations. We investigate their stellar, ISM, and dynamical properties, including a pixel-level SED analysis to derive sub-kpc-resolution stellar-mass and Av maps. We find that 850.1 is one of the most massive and highly obscured, Av~5, galaxies known at z>4 with M*~10^11.8 Mo (likely forming at z>6), and 850.2 is one of the least massive and least obscured, Av~1, members of the z>4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z>3-4 based on imaging at <2um, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar mass distribution, Re(mass)~1kpc, but its expected evolution to z~1.5 and then z~0 matches both the properties of massive, quiescent galaxies at z~1.5 and ultra-massive early-type galaxies at z~0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature which we link to the active dynamical environment it resides within.

We use the Galaxy and Mass Assembly (GAMA) and the Deep Extragalactic Visible Legacy Survey (DEVILS) observational data sets to calculate the cosmic star formation rate (SFR) and active galactic nuclei (AGN) bolometric luminosity history (CSFH/CAGNH) over the last 12.5 billion years. SFRs and AGN bolometric luminosities were derived using the spectral energy distribution fitting code ProSpect, which includes an AGN prescription to self consistently model the contribution from both AGN and stellar emission to the observed rest-frame ultra-violet to far-infrared photometry. We find that both the CSFH and CAGNH evolve similarly, rising in the early Universe up to a peak at look-back time $\approx 10$~Gyr ($z \approx 2$), before declining toward the present day. The key result of this work is that we find the ratio of CAGNH to CSFH has been flat ($\approx 10^{42.5}\mathrm{erg \, s^{-1}M_{\odot}^{-1}yr}$) for $11$~Gyr up to the present day, indicating that star formation and AGN activity have been coeval over this time period. We find that the stellar masses of the galaxies that contribute most to the CSFH and CAGNH are similar, implying a common cause, which is likely gas inflow. The depletion of the gas supply suppresses cosmic star formation and AGN activity equivalently to ensure that they have experienced similar declines over the last 10 Gyr. These results are an important milestone for reconciling the role of star formation and AGN activity in the life cycle of galaxies.

We investigate the dynamic evolution of gaseous region around FS~CMa post-mergers. Due to the slow rotation of a central B-type star, the dynamics is driven mainly by the magnetic field of the central star. Recent observations have allowed us to set a realistic initial conditions such as, the magnetic field value ($B_\star\approx6\times10^{3}G$), the mass of the central star ($M_\star=6M_\odot$), and the initial disc density $\rho_{d0}\in[10^{-13}\mathrm{g\,cm^{-3}},10^{-11}\mathrm{g \, cm^{-3}}] $. We use the PLUTO code to perform 2.5D-MHD simulations of thin and thick discs models. Especially relevant for the interpretation of the observed properties of FS~CMa post-mergers are the results for low-density discs, in which we find formation of a jet emerging from inner edge of the disc, as well as the formation of the so called "hot plasmoid" in the corona region. Jets are probably detected as discrete absorption components in the resonance lines of FS~CMa stars. Moreover, the magnetic field configuration in the low-density plasma region, favors the appearance of magnetocentrifugal winds from the disc. The currents toward the star created by the magnetic field may explain accidentally observed material infall. The disc structure is significantly changed due to the presence of the magnetic field. The magnetic field is also responsible for the formation of a hot corona as observed in several FS~CMa stars through the Raman lines. Our results are valid for all magnetic stars surrounded by a low density plasma, i.e., some of stars showing the B[e] phenomenon.

Core-collapse supernovae are among the most powerful explosions in the universe, emitting thermal neutrinos that carry away the majority of the gravitational binding energy released. These neutrinos create a diffuse supernova neutrino background (DSNB), one of the largest energy budgets among all radiation backgrounds. Detecting the DSNB is a crucial goal of modern high-energy astrophysics and particle physics, providing valuable insights in both core-collapse modeling, neutrino physics, and cosmic supernova rate history. In this review, we discuss the key ingredients of DSNB calculation and what we can learn from future detections, including black-hole formation and non-standard neutrino interactions. Additionally, we provide an overview of the latest updates in neutrino experiments, which could lead to the detection of the DSNB in the next decade. With the promise of this breakthrough discovery on the horizon, the study of DSNB holds enormous potential for advancing our understanding of the Universe.

We extend the study of the radio emission in early-type galaxies (ETGs) in the nearby Universe (recession velocity <7,500 km/s) as seen by the 150 MHz Low-Frequency ARray (LOFAR) observations and extend the sample from giant ETGs to massive (~6x10^10 - 3x10^11 solar masses) ETGs (mETGS) with -25 < MK < -23.5. Images from the second data release of the LOFAR Two-metre Sky Survey were available for 432 mETGs, 48% of which are detected above a typical luminosity of ~3x10^20 W/Hz. Most (85%) of the detected sources are compact, with sizes <4 kpc. The radio emission of 31 mETGs is extended on scales ranging from 2 to 180 kpc (median 12 kpc). In several cases, it is aligned with the host galaxy. We set a limit of ~1% to the fraction of remnant or restarted objects, which is ~16% of the extended sources. We found that the properties of the radio sources are connected with the stellar mass of the ETGs (the median radio power, the fraction of extended radio sources, and the link with the large-scale environment). However, these results only describe statistical trends because the radio properties of sources of similar stellar mass and environment show a large spread of radio properties. These trends break at the lowest host luminosities (MK>-24.5). This effect is strengthened by the analysis of even less massive ETGs, with -23.5 < MK < -21.5. This suggests that at a mass of ~2x10^11 solar masses, a general transition occurs from radio emission produced from radio-loud active galactic nuclei (AGN) to processes related to the host galaxy and (or) radio quiet AGN. At this luminosity, a transition in the stellar surface brightness profile also occurs from Sersic galaxies to those with a depleted stellar core, the so-called core galaxies. This finding is in line with previous results that indicated that only core galaxies host radio-loud AGN.

As an observer from within the Milky Way, it is difficult to determine its global structure. Despite extensive observational data from surveys at different wavelengths, we have no conclusive description of the structure of our own Galaxy. For very-high-energy (VHE) $\gamma$ rays, the most comprehensive catalogue of Galactic sources resulting from the H.E.S.S. Galactic Plane Survey (HGPS) shows a striking asymmetry in the distribution of the sources in the latitudinal direction. This could be the result of a local feature in the spatial distribution of the sources or it could be due to the position of the Sun above the Galactic plane. In this contribution, we estimate the position of the Sun based on the latitudinal flux profile of VHE $\gamma$-ray sources, assuming three mirror-symmetric models for the spatial distribution of the sources in three-dimensional space and taking into account the observational bias of the HGPS. We verify our method using simulations and find values for $z_{\odot}$ between $-6\,\mathrm{pc}$ and $94\,\mathrm{pc}$ depending on the considered model. Our results show that the position of the Sun has a significant impact on the observed source distribution and must therefore be taken into account when modelling the population of Galactic VHE $\gamma$ sources. However, it is not conclusive whether the Sun's offset from the Galactic plane is the only factor leading to the asymmetry in the latitudinal profile.

Cosmic Microwave Background (CMB) independent approaches are frequently used in the literature to provide estimates of Hubble constant ($H_0$). In this work, we report CMB independent constraints on $H_0$ in an anisotropic extension of $\Lambda$CDM model using the Big Bang Nucleosynthesis (BBN), Baryonic Acoustic Oscillations (BAO), Cosmic Chronometer (CC), and Pantheon+ compilation of Type Ia supernovae data. The analyses of the anisotropic model with different combinations of data sets reveal that anisotropy is strongly and positively correlated with $H_0$, and an anisotropy of the order $10^{-14}$ in the anisotropic model yields values of $H_0$ consistent with the direct local measurements.

We investigate the linear and nonlinear properties of the Goldreich-Schubert-Fricke (GSF) instability in stellar radiative zones with arbitrary local (radial and latitudinal) differential rotation. This instability may lead to turbulence that contributes to redistribution of angular momentum and chemical composition in stars. In our local Boussinesq model, we investigate varying the orientation of the shear with respect to the 'effective gravity', which we describe using the angle $\phi$. We first perform an axisymmetric linear analysis to explore the effects of varying $\phi$ on the local stability of arbitrary differential rotations. We then explore the nonlinear hydrodynamical evolution in three dimensions using a modified shearing box. The model exhibits both the diffusive GSF instability, and a non-diffusive instability that occurs when the Solberg-H\{o}iland criteria are violated. We observe the nonlinear development of strong zonal jets ("layering" in the angular momentum) with a preferred orientation in both cases, which can considerably enhance turbulent transport. By varying $\phi$ we find the instability with mixed radial and latitudinal shears transports angular momentum more efficiently (particularly if adiabatically unstable) than cases with purely radial shear $(\phi = 0)$. By exploring the dependence on box size, we find the transport properties of the GSF instability to be largely insensitive to this, implying we can meaningfully extrapolate our results to stars. However, there is no preferred length-scale for adiabatic instability, which therefore exhibits strong box-size dependence. These instabilities may contribute to the missing angular momentum transport required in red giant and subgiant stars and drive turbulence in the solar tachocline.

We present WarmSPy, a numerical code in Python designed to solve for the perturbations' equations in warm inflation models and compute the corresponding scalar power spectrum at CMB horizon crossing. In models of warm inflation, a radiation bath of temperature $T$ during inflation induces a dissipation (friction) rate of strength $Q \propto T^c/\phi^m$ in the equation of motion for the inflaton field $\phi$. While for a temperature-independent dissipation rate ($c=0$) an analytic expression for the scalar power spectrum exists, in the case of a non-zero value for $c$ the set of equations can only be solved numerically. For $c>0$ ($c<0$), the coupling between the perturbations in the inflaton field and radiation induces a growing (decaying) mode in the scalar perturbations, generally parameterized by a multiplicative function $G(Q)$ which we refer to as the scalar dissipation function. Using WarmSPy, we provide an analytic fit for $G(Q)$ for the cases of $c=\{3,1,-1\}$, corresponding to three cases that have been realized in physical models. Compared to previous literature results, our fits are more robust and valid over a broader range of dissipation strengths $Q\in[10^{-7},10^{4}]$. Additionally, for the first time, we numerically assess the stability of the scalar dissipation function against various model parameters, inflationary histories as well as the effects of metric perturbations. As a whole, the results do not depend appreciably on most of the parameters in the analysis, except for the dissipation index $c$, providing evidence for the universal behaviour of the scalar dissipation function $G(Q)$.

We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3\sigma$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4\sigma$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal.

We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3\sigma$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years.

Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^{\alpha}$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $\alpha=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $\alpha=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $\alpha=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2\sigma$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array.

Observing and timing a group of millisecond pulsars (MSPs) with high rotational stability enables the direct detection of gravitational waves (GWs). The GW signals can be identified from the spatial correlations encoded in the times-of-arrival of widely spaced pulsar-pairs. The Chinese Pulsar Timing Array (CPTA) is a collaboration aiming at the direct GW detection with observations carried out using Chinese radio telescopes. This short article serves as a `table of contents' for a forthcoming series of papers related to the CPTA Data Release 1 (CPTA DR1) which uses observations from the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Here, after summarizing the time span and accuracy of CPTA DR1, we report the key results of our statistical inference finding a correlated signal with amplitude $\log A_{\rm c}= -14.4 \,^{+1.0}_{-2.8}$ for spectral index in the range of $\alpha\in [-1.8, 1.5]$ assuming a GW background (GWB) induced quadrupolar correlation. The search for the Hellings-Downs (HD) correlation curve is also presented, where some evidence for the HD correlation has been found that a 4.6-$\sigma$ statistical significance is achieved using the discrete frequency method around the frequency of 14 nHz. We expect that the future International Pulsar Timing Array data analysis and the next CPTA data release will be more sensitive to the nHz GWB, which could verify the current results.

We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background.

Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.

The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.

The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data.

Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems.

This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data.

Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 years. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning of 2021, with individual pulsar datasets with timespans ranging from 14 to 25 years. These lead to improved constraints on annual parallaxes, secular variation of the orbital period, and Shapiro delay for a number of sources. Based on these results, we derived astrophysical parameters that include distances, transverse velocities, binary pulsar masses, and annual orbital parallaxes.

The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts.

We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature.

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy.

Pulsar Timing Array experiments probe the presence of possible scalar/pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results show that ultralight particles with masses $10^{-24.0}~\text{eV} \lesssim m \lesssim 10^{-23.2}~\text{eV}$ cannot constitute $100\%$ of the measured local dark matter density, but can have at most local density $\rho\lesssim 0.15$ GeV/cm$^3$.

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $\gamma=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission.

We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes ``Murriyang'' radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release.

The chromospheric Lyman-alpha line of neutral hydrogen (Ly$\alpha$; 1216 \r{A}) is the most intense emission line in the solar spectrum, yet until recently observations of flare-related Ly$\alpha$ emission have been scarce. Here, we examine the relationship between nonthermal electrons accelerated during the impulsive phase of three M3 flares that were co-observed by RHESSI, GOES, and SDO, and the corresponding response of the chromosphere in Ly$\alpha$. Despite having identical X-ray magnitudes, these flares show significantly different Ly$\alpha$ responses. The peak Ly$\alpha$ enhancements above quiescent background for these flares were 1.5%, 3.3%, and 6.4%. However, the predicted Ly$\alpha$ enhancements from FISM2 were consistently <2.5%. By comparing the properties of the nonthermal electrons derived from spectral analysis of hard X-ray observations, flares with a harder spectral index were found to produce a greater Ly$\alpha$ enhancement. The percentage of nonthermal energy radiated by the Ly$\alpha$ line during the impulsive phase was found to range from 2.0-7.9%. Comparatively, the radiative losses in He II (304 \r{A}) were found to range from 0.6-1.4% of the nonthermal energy while displaying enhancements above the background of 7.3-10.8%. FISM2 was also found to underestimate the level of He II emission in two out of the three flares. These results may have implications for space weather studies and modelling the response of the terrestrial atmosphere to changes in the solar irradiance, and will guide the interpretation of flare-related Ly$\alpha$ observations that will become available during Solar Cycle 25.

Numerous observations have shown that almost all galaxies in our Universe host supermassive black holes (SMBHs), but there is still much debate about their formation and evolutionary processes. Recently, gravitational waves (GWs) have been expected to be a new and important informative observation, in particular, in the low-frequency region by making use of the Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays (PTAs). As an evolutionary process of the SMBHs, we revisit a dark matter (DM) halo-SMBH coevolution model based on the halo merger tree employing an ansatz for the mass relation between the DM halos and the SMBHs at $z=6$. In this model, the mass of SMBHs grows through their mergers associated with the halo mergers, and hence the evolutionary information must be stored in the GWs emitted at the mergers. We investigate the stochastic gravitational background from the coalescing SMBH binaries, which the PTAs can detect, and also the GW bursts emitted at the mergers, which can be detected by the mHz band observations such as LISA. We also discuss the possibility of probing the mass relation between the DM halos and the SMBHs at high redshift by future GW observations.

We set an upper limit for the total TeV neutrino flux expected from the disk of our Galaxy in the region $|l|<30^{\circ}$ and $|b|<2^{\circ}$ probed by the ANTARES experiment. We include both the diffuse emission, due to the interaction of cosmic rays with the interstellar medium, and the possible contribution produced by gamma-ray Galactic sources. The neutrino diffuse emission is calculated under different assumptions for the cosmic ray spatial and energy distribution in our Galaxy. In particular, we assume that the total gamma-ray flux produced by all the sources, resolved and unresolved by H.E.S.S., is produced via hadronic interaction and, hence, is coupled with neutrino emission. We compare our total neutrino flux with the recent ANTARES measurement of the neutrino from the Galactic Ridge. We show that the ANTARES best-fit flux requires the existence of a large source component, close to or even larger than the most optimistic predictions obtained with our approach.

Through the use of axisymmetric 2D hydrodynamic simulations, we further investigate laterally propagating flames in X-ray bursts (XRBs). Our aim is to understand the sensitivity of a propagating helium flame to different nuclear physics. Using the Castro simulation code, we confirm the phenomenon of enhanced energy generation shortly after a flame is established after by adding ${}^{12}$C(p, ${\gamma}$)${}^{13}$N(${\alpha}$, p)${}^{16}$O to the network, in agreement with the past literature. This sudden outburst of energy leads to a short accelerating phase, causing a drastic alteration in the overall dynamics of the flame in XRBs. Furthermore, we investigate the influence of different plasma screening routines on the propagation of the XRB flame. We finally examine the performance of simplified-SDC, a novel approach to hydrodynamics and reaction coupling incorporated in Castro, as an alternative to operator-splitting.

We present an approach for accurate estimation of the covariance of 2-point correlation functions that requires fewer mocks than the standard mock-based covariance. This can be achieved by dividing a set of mocks into jackknife regions and fitting the correction term first introduced in Mohammad & Percival (2022), such that the mean of the jackknife covariances corresponds to the one from the mocks. This extends the model beyond the shot-noise limited regime, allowing it to be used for denser samples of galaxies. We test the performance of our fitted jackknife approach, both in terms of accuracy and precision, using lognormal mocks with varying densities and approximate EZmocks mimicking the DESI LRG and ELG samples in the redshift range of z = [0.8, 1.2]. We find that the Mohammad-Percival correction produces a bias in the 2-point correlation function covariance matrix that grows with number density and that our fitted jackknife approach does not. We also study the effect of the covariance on the uncertainty of cosmological parameters by performing a full-shape analysis. We find that our fitted jackknife approach based on 25 mocks is able to recover unbiased and as precise cosmological parameters as the ones obtained from a covariance matrix based on 1000 or 1500 mocks, while the Mohammad-Percival correction produces uncertainties that are twice as large. The number of mocks required to obtain an accurate estimation of the covariance for 2-point correlation function is therefore reduced by a factor of 40-60.

We have performed the first systematic search of the full GALEX data archive for astrophysical variability on timescales of seconds to minutes by rebinning data across the whole mission to 30-second time resolution. The result is the GALEX Flare Catalog (GFCAT) which describes 1426 ultraviolet variable sources, including stellar flares, eclipsing binaries, $\delta$ Scuti and RR Lyrae variables, and Active Galactic Nuclei (AGN). Many of these sources have never previously been identified as variable. We have also assembled a table of observations of ultraviolet flares and accompanying statistics and measurements, including energies, and of candidate eclipsing stars. This effort was enabled by a significantly-enhanced version of the gPhoton software for analyzing time-domain GALEX data; this gPhoton2 package is available to support follow-on efforts.

We revisit the evidence for CMB birefringence in the context of a rich Axiverse. Using probability density functions (PDFs) for various axion parameters, such as the mass and axion decay constant, we construct the PDF for the cosmic birefringence angle and investigate its properties. By relating the observed value of the birefringence angle to the mean or standard deviation of the constructed PDF, we constrain the shape of the input PDFs, providing insights into the statistical distribution of the Axiverse. We focus on three different types of axion potentials: cosine, quadratic, and asymptotically linear axion monodromy. Our analysis showcases the potential of cosmic birefringence in constraining the distribution of axion parameters and uncovering possible correlations among them. We additionally offer predictions for "birefringence tomography," anticipating future measurements of birefringence from lower multipoles, and show how it can be used to rule out simpler versions of the Axiverse. Our findings contribute to the ongoing exploration of the Axiverse and its implications for cosmic birefringence.

Transient surveys are a vital tool in exploring the dynamic universe, with radio transients acting as beacons for explosive and highly energetic astrophysical phenomena. However, performing commensal transient surveys using radio imaging can require a significant amount of computing power, data storage and time. With the instrumentation available to us, and with new and exciting radio interferometers in development, it is essential that we develop efficient methods to probe the radio transient sky. In this paper, we present results from an commensal short duration transient survey, on time scales of 8 seconds, 128 seconds and 1 hour, using data from the MeerKAT radio telescope. The dataset used was obtained as part of a galaxy observing campaign, and we focus on the field of NGC 5068. We present a quick, wide field imaging strategy to enable fast imaging of large datasets, and develop methods to efficiently filter detected transient candidates. No transient candidates were identified on the time scales of 8 seconds, 128 seconds and 1 hour, leading to competitive limits on the transient surface densities of $6.7{\times}10^{-5}$ deg$^{-1}$, $1.1{\times}10^{-3}$ deg$^{-1}$, and $3.2{\times}10^{-2}$ deg$^{-1}$ at sensitivities of 56.4 mJy, 19.2 mJy, and 3.9 mJy for the respective time scales. We find one possible candidate that could be associated with a stellar flare, that was rejected due to strict image quality control. Further short time-scale radio observations of this candidate could give definite results to its origin.

Weak Impulsive Narrowband Quiet Sun Emissions (WINQSEs) are a newly discovered class of radio emission from the solar corona. These emissions are characterized by their extremely impulsive, narrowband and ubiquitous nature. We have systematically been working on their detailed characterization, including their strengths, morphologies, temporal characteristics, energies, etc. This work is the next step in this series and focuses on the spectral nature of WINQSEs. Given that their strength is only a few percent of the background solar emission, we have adopted an extremely conservative approach to reliably identify WINQSES. Only a handful of WINQSEs meet all of our stringent criteria. Their flux densities lie in the 20 $-$ 50 Jy range and they have compact morphologies. For the first time, we estimate their bandwidths and find them to be less than 700 kHz, consistent with expectations based on earlier observations. Interestingly, we also find similarities between the spectral nature of WINQSEs and the solar radio spikes. This is consistent with our hypothesis that the WINQSEs are the weaker cousins of the type-III radio bursts and are likely to be the low-frequency radio counterparts of the nanoflares, originally hypothesized as a possible explanation for coronal heating.

Although deuterium (D) on Mars has received substantial attention, the deuterated ionosphere remains relatively unstudied. This means that we also know very little about non-thermal D escape from Mars, since it is primarily driven by excess energy imparted to atoms produced in ion-neutral reactions. Most D escape from Mars is expected to be non-thermal, highlighting a gap in our understanding of water loss from Mars. In this work, we set out to fill this knowledge gap. To accomplish our goals, we use an upgraded 1D photochemical model that fully couples ions and neutrals and does not assume photochemical equilibrium. To our knowledge, such a model has not been applied to Mars previously. We model the atmosphere during solar minimum, mean, and maximum, and find that the deuterated ionosphere behaves similarly to the H-bearing ionosphere, but that non-thermal escape on the order of 8000-9000 cm$^{-2}$s$^{-1}$ dominates atomic D loss under all solar conditions. The total fractionation factor, $f$, is $f=0.04$--0.07, and integrated water loss is 147--158 m GEL. This is still less than geomorphological estimates. Deuterated ions at Mars are likely difficult to measure with current techniques due to low densities and mass degeneracies with more abundant H ions. Future missions wishing to measure the deuterated ionosphere in situ will need to develop innovative techniques to do so.

The Fermi-LAT observations of SN 2023ixf, a Type II supernova in the nearby Pinwheel Galaxy, Messier 101 (M101), presents us with an excellent opportunity to constrain MeV-scale Axion-Like Particles (ALPs). By examining the photon decay signature from heavy ALPs that could be produced in the explosion, we improve the existing constraints on the ALP-photon coupling by up to a factor of $ \sim 2 $ for masses $ m_a \lesssim 3 $ MeV, with the exact value depending mostly on plasma properties of the collapsing core. This study demonstrates the relevance of core-collapse supernovae, also beyond the Magellanic Clouds, as probes of fundamental physics.

Dust coagulation in protoplanetary disks is not straightforward and is subject to several slow-down mechanisms, such as bouncing, fragmentation and radial drift to the star. Furthermore, dust grains in UV-shielded disk regions are negatively charged due to collisions with the surrounding electrons and ions, which leads to their electrostatic repulsion. For typical disk conditions, the relative velocities between micron-size grains are small and their collisions are strongly affected by the repulsion. On the other hand, collisions between pebble-size grains can be too energetic, leading to grain fragmentation. The aim of the present paper is to study a combined effect of the electrostatic and fragmentation barriers on dust evolution. We numerically solve the Smoluchowski coagulation-fragmentation equation for grains whose charging occurs under conditions typical for the inner disk regions, where thermal ionization operates. We find that dust fragmentation efficiently resupplies the population of small grains under the electrostatic barrier. As a result, the equilibrium abundance of sub-micron grains is enhanced by several orders of magnitude compared to the case of neutral dust. For some conditions with fragmentation velocities $\sim 1$ m s$^{-1}$, macroscopic grains are completely destroyed.

We consider in detail a possibility that the observed neutrino oscillations are due to refraction on ultralight scalar boson dark matter. We introduce the refractive mass squared, $\widetilde{m}^2$, and study its properties: dependence on neutrino energy, state of the background, etc. If the background is in a state of cold gas of particles, $\widetilde{m}^2$ shows a resonance dependence on energy. Above the resonance ($E \gg E_R $), we find that $\widetilde{m}^2$ has the same properties as the usual vacuum mass. Below the resonance, $\widetilde{m}^2$ decreases with energy, which allows to avoid the cosmological bound on the sum of neutrino masses. We consider the validity of the results: effects of multiple interactions with scalars, and modification of the dispersion relation. We show that for values of parameters of the system required to reproduce the observed neutrino masses, perturbativity is broken at low energies, which border above the resonance. If the background is in the state of coherent classical field, the refractive mass does not depend on energy but may show time dependence. It coincides with the refractive mass in a cold gas at high energies. Refractive nature of neutrino mass can be tested by searches of its dependence on energy.

Metamaterials based on regular two-dimensional arrays of thin wires have attracted renewed attention in light of a recently proposed strategy to search for dark matter axions. When placed in the external magnetic field, such metamaterials facilitate resonant conversion of axions into plasmons near their plasma frequency. Since the axion mass is not known a priori, a practical way to tune the plasma frequency of metamaterial is required. In this work, we have studied a system of two interpenetrating rectangular wire lattices where their relative position is varied. The plasma frequency as a function of their relative position in two dimensions has been mapped out experimentally, and compared with both a semi-analytic theory of wire-array metamaterials and numerical simulations. Theory and simulation yield essentially identical results, which in turn are in excellent agreement with experimental data. Over the range of translations studied, the plasma frequency can be tuned over a range of 16%.

In this article, on the basis of the Langevin equation applied to velocity fluctuations, we numerically model the Partial Variance of Increments, which is a useful tool to measure time and spatial correlations in space plasmas. We consider a Coupled Map Lattice model to relate the spatial scale of fluctuations, $k$, to some macro parameters of the systems, as the Reynolds number, $R_\lambda$, the $\kappa$ parameter of Kappa distributions, and a skewness parameter, $\delta$. To do so, we compute the Velocity Probability Density Function (PDF) for each spatial scale and different values of Reynolds number in the simulations. We fit the PDF with a Skew-Kappa distribution, and we obtain a numerical relationship between the level of turbulence of the plasma and the skewness of obtained distributions; namely $\langle \delta \rangle \sim R_\lambda^{-1/2}$. We expect the results exposed in this paper to be useful as a tool to characterize the turbulence in the context of space plasma and other environments.

Applying TQRPA calculations of Gamow--Teller strength functions in hot nuclei, we compute the (anti)neutrino spectra and energy loss rates arising from weak processes on hot $^{56}$Fe under pre-supernova conditions. We use a realistic pre-supernova model calculated by the stellar evolution code MESA. Taking into account both charged and neutral current processes, we demonstrate that weak reactions with hot nuclei can produce high-energy (anti)neutrinos. We also show that, for hot nuclei, the energy loss via (anti)neutrino emission is significantly larger than that for nuclei in their ground state. It is found that the neutral current de-excitation via the $\nu\bar\nu$-pair emission is presumably a dominant source of antineutrinos. In accordance with other studies, we confirm that the so-called single-state approximation for neutrino spectra might fail under certain pre-supernova conditions. }

It is well known that, during inflation, the conformal invariance of the electromagnetic action has to be broken in order to produce magnetic fields of observed strengths today. Often, to further enhance the strengths of the magnetic fields, parity is also assumed to be violated when the fields are being generated. In this work, we examine the evolution of the quantum state of the Fourier modes of the non-conformally coupled and parity violating electromagnetic field during inflation. We utilize tools such as the Wigner ellipse, squeezing parameters and quantum discord to understand the evolution of the field. We show that the violation of parity leads to an enhancement of the squeezing amplitude and the quantum discord (or, equivalently, in this context, the entanglement entropy) associated with a pair of opposite wave vectors for one of the two states of polarization (and a suppression for the other state of polarization), when compared to the case wherein parity is conserved. We highlight the similarities between the evolution of the Fourier modes of the electromagnetic field when parity is violated during inflation and the behavior of the modes of a charged, quantum, scalar field in the presence of a constant electric field in a de Sitter universe. We briefly discuss the implications of the results we obtain.

The low-frequency gravitons correspond to typical wavelengths that left the Hubble radius during the early inflationary stages of expansion and reentered after matter radiation equality. Consequently the temperature and the polarization anisotropies of the cosmic microwave background constrain the tensor-to-scalar-ratio in the aHz region but since the audio band and the MHz domain are sensitive to the post-inflationary expansion rate, the low-frequency determinations of the tensor-to-scalar-ratio can be combined with the high-frequency constraints. In this framework we examine the possibility that the low-frequency gravitons remain invisible in the aHz region but are still potentially detectable at much higher frequencies. Because the number of $e$-folds associated with the exit of the cosmic microwave background wavelengths depends both on the slow-roll parameters and on the total expansion rate after inflation, this approach leads to a set of lower bounds on the tensor-to-scalar-ratio.