We present results from the search for an isotropic gravitational-wave background using Advanced LIGO and Advanced Virgo data from O1 through O4a, the first part of the fourth observing run. This background is the accumulated signal from unresolved sources throughout cosmic history and encodes information about the merger history of compact binaries throughout the Universe, as well as exotic physics and potentially primordial processes from the early cosmos. Our cross-correlation analysis reveals no statistically significant background signal, enabling us to constrain several theoretical scenarios. For compact binary coalescences which approximately follow a 2/3 power-law spectrum, we constrain the fractional energy density to $\Omega_{\rm GW}(25{\rm Hz})\leq 2.0\times 10^{-9}$ (95% cred.), a factor of 1.7 improvement over previous results. Scale-invariant backgrounds are constrained to $\Omega_{\rm GW}(25{\rm Hz})\leq 2.8\times 10^{-9}$, representing a 2.1x sensitivity gain. We also place new limits on gravity theories predicting non-standard polarization modes and confirm that terrestrial magnetic noise sources remain below detection threshold. Combining these spectral limits with population models for GWTC-4, the latest gravitational-wave event catalog, we find our constraints remain above predicted merger backgrounds but are approaching detectability. The joint analysis combining the background limits shown here with the GWTC-4 catalog enables improved inference of the binary black hole merger rate evolution across cosmic time. Employing GWTC-4 inference results and standard modeling choices, we estimate that the total background arising from compact binary coalescences is $\Omega_{\rm CBC}(25{\rm Hz})={0.9^{+1.1}_{-0.5}\times 10^{-9}}$ at 90% confidence, where the largest contribution is due to binary black holes only, $\Omega_{\rm BBH}(25{\rm Hz})=0.8^{+1.1}_{-0.5}\times 10^{-9}$.

We investigate the association between galaxies and neutral OI absorption systems at z~6, which trace metal-enriched gas during the epoch of reionization. We identify 40 galaxies across six quasar fields, residing in 15 overdensities within 300 kpc of the background sightlines. Five OI absorption systems are associated with five of these overdensities, yielding a covering fraction of $0.27^{+0.13}_{-0.10}$ within 300 kpc. The absorption occurs beyond typical virial radii, indicating that the gas traces extended overdensity environments rather than individual galaxy halos, unlike the z~0 CGM which is largely bound to halos. These galaxy-associated absorbers account for $\sim35\%$ of all OI systems seen in blind quasar surveys, implying the remainder arise in lower-mass galaxies below our detection threshold or in dense neutral IGM pockets. The CGM around these galaxies contains $\gtrsim 2\times10^6~M_{\odot}$ of oxygen, comparable to the ISM oxygen mass of the galaxies themselves, suggesting that the surrounding environment holds as much metal mass as the galaxies. All five galaxy-associated systems show significantly higher $\log(N_{\rm CII}/N_{\rm OI})$ ratios than absorbers lacking galaxy associations. Furthermore, relative abundance ratios ([Si/O], [C/O]) reveal that four of the five exhibit enrichment patterns consistent with Population III nucleosynthesis. These rare systems offer a unique window into the role of first-generation stars in shaping the early metal enrichment of galaxies and their environments.

We present a detailed analysis of the kinematics of SiO maser stars around the center of the Milky Way, Sagittarius A* (Sgr A*). We used the archive data in the SiO v=1, J=2-1 emission line obtained by the Atacama Large Millimeter/Submillimeter Array (ALMA) in 2017 and 2021 (#2016.1.00940.S, PI Darling, J. and #2019.1.00292.S, PI Paine, J.). We detected 37 SiO maser stars in the channel maps and derived their angular offsets relative to Sgr A* and LSR radial velocities. We derived the proper motions of 35 stars by comparing their angular offsets in the two epochs. The proper motions of Wolf-Rayet and O star in the Nuclear Star Cluster (NSC) are reported to be rather random, except for the co-moving clusters IRS13E and IRS13N (Tsuboi et al. 2022). However, the derived proper motions of SiO maser stars do not look completely random. The proper motions of the SiO maser stars show a tendency to lie along the Galactic plane. The proper motion amplitudes of SiO maser stars are larger than the LSR velocity amplitudes. We estimated the 3D motions from the proper motions and LSR velocities. Many 3D velocities are near to or larger than the upper limit velocities for Kepler orbits around Sgr A*, whose mass is assumed to be 4x10^6 Msun. These indicate that the SiO maser stars around Sgr A* are members of the Nuclear Stellar Disk rather than the NSC.

We explore the properties of interferometric data from high-redshift 21~cm measurements using the Murchison Widefield Array. These data contain redshifted 21~cm signal, contamination from continuum foreground sources, and radiometric noise. The 21~cm signal from the Epoch of Reionization is expected to be highly-Gaussian, which motivates the use of the power spectrum as an effective statistical tool for extracting astrophysical information. We find that foreground contamination introduces non-Gaussianity into the distribution of measurements, and then use this information to separate Gaussian from non-Gaussian signal. We present improved upper limits on the 21cm EoR power spectrum from the MWA using a Gaussian component of the data, based on the existing analysis from Nunhokee et al (2025). This is extracted as the best-fitting Gaussian to the measured data. Our best 2 sigma (thermal+sample variance) limit for 268 hours of data improves from (30.2~mK)^2 to (23.0~mK)^2 at z=6.5 for the EW polarisation, and from (39.2~mK)^2 to (21.7~mK)^2 = 470~mK^2 in NS. The best limits at z=6.8 (z=7.0) improve to P < (25.9~mK)^2 (P < (32.0~mK)^2), and k = 0.18h/Mpc (k = 0.21h/Mpc). Results are compared with realistic simulations, which indicate that leakage from foreground contamination is a source of the non-Gaussian behaviour.

We utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) Y1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude $S_8=\sigma_8\sqrt{\Omega_m/0.3}$. The large-scale structure's statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF -- a complementary higher-order statistic. By defining the PDF over the standard deviation-normalized convergence map we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find $S_8=0.860^{+0.066}_{-0.109}$ from the PDF alone and $S_8=0.798^{+0.029}_{-0.042}$ from the combination of PDF and power spectrum (68% CL). The PDF improves the power spectrum-only constraint by about 10%.

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

We report discovery of two CO clouds which are likely falling down to the Galactic plane at more than $35$ km s$^{-1}$. The clouds show head-tail distributions elongated perpendicular to the Galactic plane at $l=331.6^{\circ}$ and $b=0^{\circ}$ as revealed by an analysis of the Mopra CO $J=$1-0 survey data. We derived the distance of the clouds to be $2.46 \pm 0.18$ kpc based on the Gaia Data Release 3. The CO clouds have molecular masses of $4.8\times 10^3\ M_{\odot}$ and $3.5\times 10^3\ M_{\odot}$, respectively, and show kinetic temperature of 30-50 K as derived from the line intensities of the $^{13}$CO $J$=2-1, $^{12}$CO $J$=1-0, and $^{13}$CO $J$=1-0 emission. The temperature in the heads of the clouds is significantly higher than 10 K of the typical molecular clouds, although no radiative heat source is found inside or close to the clouds. Based on the results, we interpret that the present clouds are falling onto the Milky Way disk and are significantly heated up by the strong shock interaction with the disk HI gas. We suggest that the clouds represent part of the HI intermediate velocity clouds falling to the Galactic plane which were converted into molecular clouds by shock compression. This is the first case of falling CO clouds having direct observed signatures of the falling motion including clear directivity and shock heating. Possible implications of the CO clouds in the evolution of the Galactic interstellar medium are discussed.

The Indian Pulsar Timing Array (InPTA) employs unique features of the upgraded Giant Metrewave Radio Telescope (uGMRT) to monitor dozens of the International Pulsar Timing Array (IPTA) millisecond pulsars (MSPs), simultaneously in the 300-500 MHz and the 1260-1460 MHz bands. This dual-band approach ensures that any frequency-dependent delays are accurately characterized, significantly improving the timing precision for pulsar observations, which is crucial for pulsar timing arrays. We present details of InPTA's second data release that involves 7 yrs of data on 27 IPTA MSPs. This includes sub-banded Times of Arrival (ToAs), Dispersion Measures (DM), and initial timing ephemerides for our MSPs. A part of this dataset, originally released in InPTA's first data release, is being incorporated into IPTA's third data release which is expected to detect and characterize nanohertz gravitational waves in the coming years. The entire dataset is reprocessed in this second data release providing some of the highest precision DM estimates so far and interesting solar wind related DM variations in some pulsars. This is likely to characterize the noise introduced by the dynamic inter-stellar ionised medium much better than the previous release thereby increasing sensitivity to any future gravitational wave search.

We present a comprehensive analysis of the relationship between galaxies and the intergalactic medium (IGM) during the late stages of cosmic reionization, based on the complete JWST EIGER dataset. Using deep NIRCam $3.5\,\mathrm{\mu m}$ slitless spectroscopy, we construct a sample of 948 [\OIII]$\lambda5008$-emitting galaxies with $-21.4\lesssim M_\mathrm{UV}\lesssim -17.2$spanning$5.335.70$. Statistical tests using mock catalogs with realistic galaxy clustering but no correlation with the transmission field confirm that the observed correlations are unlikely to arise by chance. The evolving signals can be explained by stronger absorption in overdense regions, combined with the competing influences of local radiation fields and the rising background radiation. While local radiation dominates ionization of the surrounding IGM at earlier times, the background becomes increasingly important, eventually surpassing the impact of nearby galaxies. These results support an inside-out progression of reionization, with ionized regions originating around clustered, star-forming galaxies and gradually extending into underdense regions.

Atmospheres play a crucial role in planetary habitability. Around M dwarfs and young Sun-like stars, planets receiving the same insolation as the present-day Earth are exposed to intense stellar X-rays and extreme-ultraviolet (XUV) radiation. This study explores the fundamental question of whether the atmosphere of present-day Earth could survive in such harsh XUV environments. Previous theoretical studies suggest that stellar XUV irradiation is sufficiently intense to remove such atmospheres completely on short timescales. In this study, we develop a new upper-atmospheric model and re-examine the thermal and hydrodynamic responses of the thermospheric structure of an Earth-like N2-O2 atmosphere, on an Earth-mass planet, to an increase in the XUV irradiation. Our model includes the effects of radiative cooling via electronic transitions of atoms and ions, known as atomic line cooling, in addition to the processes accounted for by previous models. We demonstrate that atomic line cooling dominates over the hydrodynamic effect at XUV irradiation levels greater than several times the present level of the Earth. Consequentially, the atmosphere's structure is kept almost hydrostatic, and its escape remains sluggish even at XUV irradiation levels up to a thousand times that of the Earth at present. Our estimates for the Jeans escape rates of N2-O2 atmospheres suggest that these 1 bar atmospheres survive in early active phases of Sun-like stars. Even around active late M dwarfs, N2-O2 atmospheres could escape significant thermal loss on timescales of gigayears. These results give new insights into the habitability of terrestrial exoplanets and the Earth's climate history.

LiteBIRD is a next-generation satellite mission to measure the polarization of the cosmic microwave background (CMB) radiation. On large angular scales the B-mode polarization of the CMB carries the imprint of primordial gravitational waves, and its precise measurement would provide a powerful probe of the epoch of inflation. The goal of LiteBIRD is to achieve a measurement of the characterizing tensor to scalar ratio $r$ to an uncertainty of $\delta r=0.001$. In order to achieve this goal we will employ a kilo-pixel superconducting detector array on a cryogenically cooled sub-Kelvin focal plane with an optical system at a temperature of 4~K. We are currently considering two detector array options; transition edge sensor (TES) bolometers and microwave kinetic inductance detectors (MKID). In this paper we give an overview of LiteBIRD and describe a TES-based polarimeter designed to achieve the target sensitivity of 2~$\mu$K$\cdot$arcmin over the frequency range 50 to 320~GHz.

Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration.

As the fields of stellar and exoplanetary study grow and revolutionary new detection instruments are created, it is imperative that a homogeneous, precise source of stellar parameters is available. This first work of the gr8stars collaboration presents the all-sky magnitude limited sample of 5645 bright FGKM dwarfs, along with homogeneously derived spectroscopic parameters of a subset of 1716 targets visible from the Northern hemisphere. We have collected high-resolution archival and new spectra from several instruments. Spectrosocpic parameters are determined using the PAWS pipeline, employing both the curve-of-growth equivalent width method, and the spectral synthesis method. We achieve median uncertainties of 106K in stellar effective temperature, 0.08 dex in surface gravity, and 0.03 dex in metallicity. This paper also presents photometric stellar parameters for these dwarfs, determined using SED fitting. The full gr8stars sample selection, including derived spectroscopic and photometric parameters, is made available through an interactive online database. We also perform a kinematic analysis to classify these stars according to their Galactic component.

We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $\lambda$3mm and $\lambda$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN. We detect high linear polarization fractions (2-15%) and large rotation measures (RM $>10^{3.3}-10^{5.5}$ rad m$^{-2}$). For Sgr A* we report a mean RM of $(-4.2\pm0.3) \times10^5$ rad m$^{-2}$ at 1.3 mm, consistent with measurements over the past decade, and, for the first time, an RM of $(-2.1\pm0.1) \times10^5$ rad m$^{-2}$ at 3 mm, suggesting that about half of the Faraday rotation at 1.3 mm may occur between the 3 mm photosphere and the 1.3 mm source. We also report the first unambiguous measurement of RM toward the M87 nucleus at mm wavelengths, which undergoes significant changes in magnitude and sign reversals on a one year time-scale, spanning the range from -1.2 to 0.3 $\times\,10^5$ rad m$^{-2}$ at 3 mm and -4.1 to 1.5 $\times\,10^5$ rad m$^{-2}$ at 1.3 mm. Given this time variability, we argue that, unlike the case of Sgr A*, the RM in M87 does not provide an accurate estimate of the mass accretion rate onto the black hole. We put forward a two-component model, comprised of a variable compact region and a static extended region, that can simultaneously explain the polarimetric properties observed by both the EHT and ALMA. These measurements provide critical constraints for the calibration, analysis, and interpretation of simultaneously obtained VLBI data with the EHT and GMVA.

Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modeling the impact of PSF systematics on the cosmic shear 2PCF, and mitigating its impact on cosmological analysis. Our framework can describe leakage and modeling error from all spin-2 quantities contributed by the PSF second and higher moments, rather than just the second moments, using the cross-correlations between galaxy shapes and PSF moments. We interpret null tests using the HSC Year 3 (Y3) catalogs with this formalism, and find that leakage from the spin-2 combination of PSF fourth moments is the leading contributor to additive shear systematics, with total contamination that is an order of magnitude higher than that contributed by PSF second moments alone. We conducted a mock cosmic shear analysis for HSC Y3, and find that, if uncorrected, PSF systematics can bias the cosmological parameters $\Omega_m$ and $S_8$ by $\sim$0.3$\sigma$. The traditional second moment-based model can only correct for a 0.1$\sigma$ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contamination for HSC Y3 cosmic shear analysis. We also reanalyze the HSC Y1 cosmic shear analysis with our updated systematics model, and identify a 0.07$\sigma$ bias on $\Omega_m$ when using the more restricted second moment model from the original analysis. We demonstrate how to self-consistently use the method in both real space and Fourier space, assess shear systematics in tomographic bins, and test for PSF model overfitting.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.