Papers for Tuesday, Feb 13 2024

Context. Weak gravitational lensing offers a powerful method to investigate the projected matter density distribution within galaxy clusters, granting crucial insights into the broader landscape of dark matter on cluster scales. Aims. In this study, we make use of the large photometric galaxy cluster data set derived from the publicly available Third Data Release of the Kilo-Degree Survey, along with the associated shear signal. Our primary objective is to model the peculiar sharp transition in the cluster profile slope, i.e. what is commonly referred to as the splashback radius. The data set under scrutiny includes 6962 galaxy clusters, selected by AMICO on the KiDS-DR3 data, in the redshift range of 0.1 < z < 0.6, all observed at signal-to-noise ratio greater than 3.5. Methods. Employing a comprehensive Bayesian analysis, we model the stacked excess surface mass density distribution of the clusters. We adopt a model from recent results on numerical simulations, that captures the dynamics of both orbiting and infalling materials, separated by the region where the density profile slope undergoes a pronounced deepening. Results. We find that the adopted profile successfully characterizes the cluster masses, consistent with previous works, and models the deepening of the slope of the density profiles measured with weak-lensing data up to the outskirts. Moreover, we measure the splashback radius of galaxy clusters and show that its value is close to the radius within which the enclosed overdensity is 200 times the mean matter density of the Universe, while theoretical models predict a larger value consistent with a low accretion rate. This points to a potential bias of optically selected clusters that are preferentially characterized by a high density at small scales compared to a pure mass-selected cluster sample.

In the hot intracluster medium (ICM) in galaxy clusters, plasma microinstabilities may play an important role in the transport of heat and momentum on the large scales. In this paper, we continue our investigation of the effect of whistler suppression of thermal conductivity on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters and contribute to the observed turbulence. We use a closure for the heat flux inspired by kinetic simulations and show that MTI turbulence with whistler suppression exhibits a critical transition: for modest suppression of the conductivity, the MTI turbulent velocities decrease in agreement with previous MTI scaling laws. However, for suppression above a critical threshold, the MTI loses its ability to maintain equipartition-level magnetic fields through a small-scale dynamo, and the system enters a ``death-spiral''. We propose a model to explain this critical transition, and speculate that conditions in the hot ICM are favourable to the upkeep of the dynamo. Additionaly, with whistler suppression high-$\beta$ regions are brought out of thermal equilibrium while the efficiency of MTI turbulent driving is reduced. Finally, we show that external turbulence interferes with the MTI and leads to lower levels of turbulence. While individually both external turbulence and whistler suppression weaken the MTI, we find that they can exhibit a complex interplay when acting in conjunction, with external turbulence boosting the whistler-suppressed thermal conductivity and even reviving a ``dead'' MTI. Our study illustrates how extending magnetohydrodynamics with a simple prescription for microscale plasma physics can lead to the formation of a complicated dynamical system and demonstrates that further work is needed in order to bridge the gap between micro- and macro scales in galaxy clusters.

We report discovery and characterization of a main-sequence G star orbiting a dark object with mass $1.90\pm 0.04\,M_{\odot}$. The system was discovered via Gaia astrometry and has an orbital period of 731 days. We obtained multi-epoch RV follow-up over a period of 600 days, allowing us to refine the Gaia orbital solution and precisely constrain the masses of both components. The luminous star is a $\gtrsim 12$ Gyr-old, low-metallicity halo star near the main-sequence turnoff ($T_{\rm eff} \approx 6000$ K; $\log\left(g/\left[{\rm cm\,s^{-2}}\right]\right)\approx 4.0; \rm [Fe/H]\approx-1.25$; $M\approx0.79\,M_{\odot}$) with a highly enhanced lithium abundance. The RV mass function sets a minimum companion mass for an edge-on orbit of $M_2 > 1.67\,M_{\odot}$, well above the Chandrasekhar limit. The Gaia inclination constraint, $i=68.8\pm 1.4$ deg, then implies a companion mass of $M_2= 1.90\pm 0.04\,M_{\odot}$. The companion is most likely a massive neutron star: the only viable alternative is two massive white dwarfs in a close binary, but this scenario is disfavored on evolutionary grounds. The system's low eccentricity ($e=0.122\pm 0.003$) disfavors dynamical formation channels and implies that the neutron star formed with very little mass loss ($\lesssim 1\,M_{\odot}$) and with a weak natal kick ($v_{\rm kick}\lesssim 10\,\rm km\,s^{-1}$). The current orbit is too small to have accommodated the neutron star progenitor as a red supergiant or super-AGB star. The simplest formation scenario -- isolated binary evolution -- requires the system to have survived stable mass transfer or common envelope evolution with a donor-to-accretor mass ratio $>10$. The system, which we call Gaia NS1, is likely a progenitor of symbiotic X-ray binaries and long-period millisecond pulsars. Its discovery challenges binary evolution models and bodes well for Gaia's census of compact objects in wide binarie

We present a high-resolution spectral study of Fe L-shell extinction by the diffuse interstellar medium (ISM) in the direction of the X-ray binaries Cygnus X-1 and GX 339-4, using the XMM-Newton reflection grating spectrometer. The majority of interstellar Fe is suspected to condense into dust grains in the diffuse ISM, but the compounds formed from this process are unknown. Here, we use the laboratory cross sections from Kortright & Kim (2000) and Lee et al. (2009) to model the absorption and scattering profiles of metallic Fe, and the crystalline compounds fayalite (Fe$_2$SiO$_4$), ferrous sulfate (FeSO$_4$), hematite ($\alpha$-Fe$_2$O$_3$), and lepidocrocite ($\gamma$-FeOOH), which have oxidation states ranging from Fe$^{0}$ to Fe$^{3+}$. We find that the observed Fe L-shell features are systematically offset in energy from the laboratory measurements. An examination of over two dozen published measurements of Fe L-shell absorption finds a 1-2 eV scatter in energy positions of the L-shell features. Motivated by this, we fit for the best energy-scale shift simultaneously with the fine structure of the Fe L-shell extinction cross sections. Hematite and lepidocrocite provide the best fits ($\approx +1.1$ eV shift), followed by fayalite ($\approx +1.8$ eV shift). However, fayalite is disfavored, based on the implied abundances and knowledge of ISM silicates gained by infrared astronomical observations and meteoritic studies. We conclude that iron oxides in the Fe$^{3+}$ oxidation state are good candidates for Fe-bearing dust. To verify this, new absolute photoabsorption measurements are needed on an energy scale accurate to better than 0.2 eV.

We update the dust model present within the Simba galaxy simulations with a self-consistent framework for the co-evolution of dust and molecular hydrogen populations in the interstellar medium, and use this to explore $z \geq 6$ galaxy evolution. In addition to tracking the evolution of dust and molecular hydrogen abundances, our model fully integrates these species into the Simba simulation, explicitly modelling their impact on physical processes such as star formation and cooling through the inclusion of a novel two-phase sub-grid model for interstellar gas. In running two high-resolution simulations down to $z \sim 6$ we find that our Simba-EoR model displays a generally tighter concordance with observational data than fiducial Simba. Additionally we observe that our Simba-EoR models increase star formation activity at early epochs, producing larger dust-to-gas ratios consequently. Finally, we discover a significant population of hot dust at $\sim 100$ K, aligning with contemporaneous observations of high-redshift dusty galaxies, alongside the large $\sim 20$ K population typically identified.

In the recent Astro2020 Decadal Report, ''Pathways to Discovery in Astronomy and Astrophysics for the 2020s'' Adaptive Optics (AO) was identified as a crucial technology for a variety of reasons. These included an emphasis on high-contrast imaging and AO systems as being part of future technology development especially with application to the two US ELT projects. Instrument upgrades were also identified for existing 4m to 10m class telescopes which would incorporate upgrades to existing AO systems. As noted in the Report: (1) ''the central role of AO instrumentation and the importance of further development are rapidly growing, with novel concepts pushing toward wider area'', (2) ''Visible AO has high potential scientific return by opening up an entire wavelength regime to high angular resolution studies. The goal is to exploit the smaller diffraction limit of telescopes in the optical, yet both the coherence length and time decrease at shorter wavelengths requiring wavefront sensing at high spatial and temporal frequencies that are currently technologically challenging. This is an important developing area for the 2020s - 2030s.'', and (3) ''Such investments in AO systems development is a key risk mitigation strategy for ELTs, whose full resolution and sensitivity potential can only be realized with AO, and which is recognized as the most important technical risk for both GMT and TMT''. A workshop was held in May, 2023 to develop a Community Response document (this document) to provide feedback and suggested priorities to various funding agencies, such as NSF, NASA, and DoE, as to the AO Research and Development priorities to meet the technical and science objectives outlined in Astro2020 for ground-based AO, both stand-alone and in support of space missions.

We present the results of models of impulsively heated coronal loops using the 1-D hydrodynamic Adaptively Refined Godunov Solver (ARGOS) code. The impulsive heating events (which we refer to as "nanoflares") are modeled by discrete pulses of energy along the loop. We explore the occurrence of cold condensations due to the effective equivalent of thermal non-equilibrium (TNE) in loops with steady heating, and examine its dependence on nanoflare timing and intensity and also nanoflare location along the loop, including randomized distributions of nanoflares. We find that randomizing nanoflare distributions, both in time/intensity and location, diminishes the likelihood of condensations as compared to distributions with regularly occurring nanoflares with the same average properties. The usual criteria that condensations are favored for heating near loop footpoints and with high cadences are more strict for randomized (as opposed to regular) nanoflare distributions, and for randomized distributions the condensations stay in the loop for a shorter amount of time. That said, condensations can sometimes occur in cases where the average values of parameters (frequency or location) are beyond the critical limits above which condensations do not occur for corresponding steady, non-randomized values of those parameters. These properties of condensations occurring due to randomized heating can be used in the future to investigate diagnostics of coronal heating mechanisms.

The detection of planetary transits in the light curves of active stars, featuring correlated noise in the form of stellar variability, remains a challenge. Depending on the noise characteristics, we show that the traditional technique that consists of detrending a light curve before searching for transits alters their signal-to-noise ratio, and hinders our capability to discover exoplanets transiting rapidly-rotating active stars. We present nuance, an algorithm to search for transits in light curves while simultaneously accounting for the presence of correlated noise, such as stellar variability and instrumental signals. We assess the performance of nuance on simulated light curves as well as on the TESS light curves of 438 rapidly-rotating M dwarfs. For each dataset, we compare our method to 5 commonly-used detrending techniques followed by a search with the Box-Least-Square algorithm. Overall, we demonstrate that nuance is the most performant method in 93% of cases, leading to both the highest number of true positives and the lowest number of false positive detections. Although simultaneously searching for transits while modeling correlated noise is expected to be computationally expensive, we make our algorithm tractable and available as the JAX-powered Python package nuance, allowing its use on distributed environments and GPU devices. Finally, we explore the prospects offered by the nuance formalism, and its use to advance our knowledge of planetary systems around active stars, both using space-based surveys and sparse ground-based observations.

We address in this work the nature and evolution of the long-period compact star sources, which has recently added several unexpected members. The central hypothesis is that particle winds drive their evolution, being an important factor for these relatively old sources. We show the consistency of this picture and remark some unsolved problems and caveats within it.

High-redshift quasars are thought to live in the densest regions of space which should be made evident by an overdensity of galaxies around them. However, campaigns to identify these overdensities through the search of Lyman Break Galaxies (LBGs) and Lyman $\alpha$ emitters (LAEs) have had mixed results. These may be explained by either the small field of view of some of the experiments, the broad redshift ranges targeted by LBG searches, and by the inherent large uncertainty of quasar redshifts estimated from UV emission lines, which makes it difficult to place the Ly-$\alpha$ emission line within a narrowband filter. Here we present a three square degree search ($\sim 1000$ pMpc) for LAEs around the $z=6.9$ quasar VIKJ2348-3054 using the Dark Energy CAMera (DECam), housed on the 4m Blanco telescope, finding 38 LAEs. The systemic redshift of VIK J2348--3054 is known from ALMA [CII] observations and place the Ly-$\alpha$ emission line of companions within the NB964 narrowband of DECam. This is the largest field of view LAE search around a $z>6$ quasar conducted to date. We find that this field is $\sim$ 10 times more overdense when compared to the Chandra Deep-Field South, observed previously with the same instrumental setup as well as several combined blank fields. This is strong evidence that VIKJ2348-3054 resides in an overdensity of LAEs over several Mpc. Surprisingly, we find a lack of LAEs within 5 physical Mpc of the quasar and take this to most likely be evidence of the quasar suppressing star formation in its immediate vicinity. This result highlights the importance of performing overdensity searches over large areas to properly assess the density of those regions of the Universe.

We use the SAMI Galaxy Survey to examine the drivers of galaxy spin, $\lambda_{R_e}$, in a multi-dimensional parameter space including stellar mass, stellar population age (or specific star formation rate) and various environmental metrics (local density, halo mass, satellite vs. central). Using a partial correlation analysis we consistently find that age or specific star formation rate is the primary parameter correlating with spin. Light-weighted age and specific star formation rate are more strongly correlated with spin than mass-weighted age. In fact, across our sample, once the relation between light-weighted age and spin is accounted for, there is no significant residual correlation between spin and mass, or spin and environment. This result is strongly suggestive that present-day environment only indirectly influences spin, via the removal of gas and star formation quenching. That is, environment affects age, then age affects spin. Older galaxies then have lower spin, either due to stars being born dynamically hotter at high redshift, or due to secular heating. Our results appear to rule out environmentally dependent dynamical heating (e.g. galaxy-galaxy interactions) being important, at least within $1R_e$ where our kinematic measurements are made. The picture is more complex when we only consider high-mass galaxies ($M_*\gtrsim 10^{11}$M$_{\odot}$). While the age-spin relation is still strong for these high-mass galaxies, there is a residual environmental trend with central galaxies preferentially having lower spin, compared to satellites of the same age and mass. We argue that this trend is likely due to central galaxies being a preferred location for mergers.

Recent Imaging X-ray Polarimetry Explorer (IXPE) observations of blazars tend to support the shock model for the X-ray emission, but report a low polarization degree ($\Pi\sim 10\%$) in X-rays compared with the previous theoretical expectations in the shock model. In order to reconcile the theoretical expectations with observations, we revisit the polarization of the shock emission by considering different kind of direction distribution for the shock-generated magnetic fields (sgMFs). Here, $w'_{\rm sg}\propto(\sin\theta')^{\zeta_{\rm sg}}$ with $\theta'=0$ along the shock normal direction is used to describe the direction distribution of the sgMFs in the shock co-moving frame. It is found that the polarization in the X-ray and radio emission for a general jet in blazars can be described as $\Pi\sim 44.5[1-\exp(-\zeta_{\rm sg}/2.6)]\%$ and $\Pi\sim 20[1-\exp(-\zeta_{\rm sg}/2.4)]\%$, respectively. Correspondingly, one can have $\zeta_{\rm sg}\sim 1-1.5$ according to the IXPE observations. Besides the sgMFs, the magnetic fields generated by the Richmyer-Meshkov instability (rmMFs) is supposed to present in the jets. The direction of the rmMFs is mainly distributed along the shock normal in the simulations and thus $w'_{\rm rm}\propto(\cos \theta')^{\zeta_{\rm rm}}$ is adopted to describe the direction distribution of rmMFs. We find that the rmMFs is likely to significantly affect the polarization properties at the low-frequency emission, especially when the sgMFs decay rapidly. Based on the contemporaneous radio and X-ray observations, we find the the emission of the electrons in the rmMFs make a significant contribution in the low-frequency emission and the the ordered background magnetic fields (obMFs) can be neglected.

Photoevaporation is thought to play an important role in the early planetary evolution. In this study, we investigate the diffusion limit of X-ray and ultraviolet induced photoevaporation in primordial atmospheres. We find that compositional fractionation resulting from mass loss is more significant than currently recognized because it is controlled by the conditions at the top of the atmosphere, where particle collisions are less frequent. Such fractionation at the top of the atmosphere develops a compositional gradient that extends downward. Mass outflow eventually reaches a steady state in which hydrogen loss is diffusion limited. We derive new analytic expressions for the diffusion-limited mass loss rate and the crossover mass.

We study how close passages of interstellar objects of planetary and substellar masses may affect the immediate and long-term dynamics of the Solar system. We consider two nominal approach orbits, namely, the orbits of actual interstellar objects 1I/'Oumuamua and 2I/Borisov, assuming them to be typical or representative for interstellar swarms of matter. Thus, the nominal orbits of the interloper in our models cross the inner part of the Solar system. Series of massive numerical experiments are performed, in which the interloper's mass is varied with a small step over a broad range. We find that, even if a Jovian-mass interloper does not experience close encounters with the Solar system planets (and this holds for our nominal orbits), our planetary system can be destabilised on timescales as short as several million years. In what concerns substellar-mass interlopers (free-floating brown dwarfs), an immediate (on a timescale of $\sim 10 - 100$ yr) consequence of such a MISO flyby is a sharp increase in the orbital eccentricities and inclinations of the outer planets. On an intermediate timescale ($\sim 10^3 - 10^5$ yr after the MISO flyby), Uranus or Neptune can be ejected from the system, as a result of their mutual close encounters and encounters with Saturn. On a secular timescale ($\sim 10^6 - 10^7$ yr after the MISO flyby), the perturbation wave formed by secular planetary interactions propagates from the outer Solar system to its inner zone.

The star formation histories (SFHs) of galaxies provide valuable insights into galaxy evolution and stellar physics. Understanding the SFHs enables the study of chemical enrichment of galaxies, star formation triggered by interactions, and the behavior of various stellar populations. This work investigates the SFHs of ten dwarf galaxies in the Local Group (LG), which spans a wide range of types, masses, luminosities, and metallicities. The analysis is based on our new sample of the member stars in the LG after removing the foreground dwarf stars by the near-infrared color-color diagram and the Gaia astrometric information. The samples include the most complete and pure red supergiants and asymptotic giant branch stars to gain valuable insights into the recent SFHs of the galaxies. The CMD fitting method is introduced to measure the SFH. The Padova isochrones are used to generate initial model CMDs, accounting for photometric errors and completeness through star field simulations to match the completeness and error distributions of the observed CMDs. Subsequently, the SFHs, distance modulus, and metallicity of the ten dwarf galaxies are determined by fitting the CMDs. The results indicate that the star formation rates (SFRs) of dwarf irregulars show a gradual increase, while those of dwarf ellipticals exhibit a gradual decrease from the past to the present. Furthermore, this work shows that the star formation activity in dwarf ellipticals persisted up to 30 Myr ago. A significant increasing feature in the SFH of NGC 6822 reveals star formation activity triggered by an interaction event.

The results of a Fourier analysis of high-precision Kepler photometry of 75 double-mode RR~Lyrae (RRd) stars observed during NASA's K2 Mission (2014-18) are presented. Seventy-two of the stars are `classical' RRd (cRRd) stars lying along a well-defined curve in the Petersen diagram and showing no evidence of Blazhko modulations. The remaining three stars are `anomalous' RRd (aRRd) stars that lie well below the cRRd curve in the Petersen diagram. These stars have larger fundamental-mode amplitudes than first-overtone amplitudes and exhibit Blazhko variations. Period-amplitude relations for the individual pulsation components of the cRRd stars are examined, as well as correlations involving Fourier phase-difference and amplitude-ratio parameters that characterize the light curves for the two radial modes. A simple statistical model relating the fundamental (P0) and first-overtone (P1) periods to [Fe/H] provides insight into the functional form of the Petersen diagram. A calibration equation for estimating [Fe/H]phot abundances of `classical' RRd stars is derived by inverting the model and using 211 field and 57 globular cluster cRRd stars with spectroscopic metallicities to estimate the model coefficients. The equation is used to obtain [Fe/H]phot for the full sample of 72 K2 cRRd stars and for 2130 cRRd stars observed by the ESA Gaia Mission. Of the 49 K2 cRRd stars that are in the Gaia DR3 catalogue only five were found to be correctly classified, the remainder having been misclassified `RRc' or `RRab'.

One of the most important objectives of solar physics is the physical understanding of the solar atmosphere, the structure of which is also described in terms of the density (N) and temperature (T) distributions of the atmospheric matter. Several multi-frequency analyses show that the characteristics of these distributions are still debated, especially for the outer coronal emission. We aim to constrain the T and N distributions of the solar atmosphere through observations in the centimetric radio domain. We employ single-dish observations from two of the INAF radio telescopes at the K-band frequencies (18 - 26 GHz). We investigate the origin of the significant brightness temperature ($T_B$) level that we detected up to the upper corona ($\sim 800$ Mm of altitude with respect to the photospheric solar surface). To probe the physical origin of the atmospheric emission and to constrain instrumental biases, we reproduced the solar signal by convolving specific 2D antenna beam models. The analysis of the solar atmosphere is performed by adopting a physical model that assumes the thermal bremsstrahlung as the emission mechanism, with specific T and N distributions. The modelled $T_B$ profiles are compared with those observed by averaging solar maps obtained during the minimum of solar activity (2018 - 2020). The T and N distributions are compatible (within $25\%$ of uncertainty) with the model up to $\sim 60$ Mm and $\sim 100$ Mm of altitude, respectively. The analysis of the role of the antenna beam pattern on our solar maps proves the physical nature of the atmospheric emission in our images up to the coronal tails seen in our $T_B$ profiles. The challenging analysis of the coronal radio emission at higher altitudes, together with the data from satellite instruments will require further multi-frequency measurements.

The Milky Way is a spiral galaxy comprising three main components: the Bulge, the Disk, and the Halo. Of particular interest is the Galactic disk, which holds a significant portion of the baryonic matter angular momentum and harbors at least two primary stellar populations: the thin and thick disks. Understanding the formation and evolution of the Galactic disk is crucial for comprehending the origins and development of our Galaxy. Stellar archaeology offers a means to probe the disk's evolution by listening to the cosmological narratives of its oldest and most pristine stars, specifically the metal-poor stars. In this study, we employed accurate photometric metallicity estimates and Gaia Early Data Release 3 astrometry to curate a pure sample of the oldest Galactic stars. This proceeding presents a summary of our primary findings.

We report on the discovery of the first ultra metal-poor (UMP) star 2MASS~J20500194$-$6613298 (J2050$-$6613; \mbox{[Fe/H] = $-4.05$}) selected from the Gaia BP/RP spectral catalog that belongs to the ancient Atari disk component. We obtained a high-resolution spectrum for the star with the MIKE spectrograph on the Magellan-Clay telescope. J2050$-$6613 displays a typical chemical abundance pattern for UMP stars, including carbon and zinc enhancements. In contrast, J2050$-$6613 shows extremely high [Sr/Fe] and [Sr/Ba] ratios compared to other stars in the [Fe/H] $<-4.0$ regime. J2050$-$6613 is most likely an early Population\,II star that formed from a gas cloud that was chemically enriched by a massive Population\,III hypernova (E $> 10^{52}$\,erg). Such a Population\,III core-collapse hypernova could simultaneously explain the origin of the abundance pattern of light and heavy elements of 2MASS~J2050$-$6613 if a large amount of Sr of $\sim10^{-5}$\,M$_{\odot}$ was produced, possibly by neutrino-driven \textbf{(wind)} ejecta. Therefore, the abundance pattern of 2MASS~J2050$-$6613 places important constraints on Sr-producing nucleosynthesis sources operating in the Atari progenitor at the earliest times.

The absorption and emission of light by exoplanet atmospheres encode details of atmospheric composition, temperature, and dynamics. Fundamentally, simulating these processes requires detailed knowledge of the opacity of gases within an atmosphere. When modeling broad wavelength ranges at high resolution, such opacity data for even a single gas can take up multiple gigabytes of system random-access memory (RAM). This aspect can be a limiting factor when considering the number of gases to consider in a simulation, the sampling strategy used for inference, or even the architecture of the system used for calculations. Here, we present cortecs, a Python tool for compressing opacity data. cortecs provides flexible methods for fitting the temperature, pressure, and wavelength dependencies of opacity data and for evaluating the opacity with accelerated, GPU-friendly methods. The package is actively developed on GitHub (https://github.com/arjunsavel/cortecs), and it is available for download with pip.

The precise measurement of neutron star (NS) spins can provide important insight into the formation and evolution of compact binaries containing NS. While traditional methods of NS spin measurement rely on pulsar observations, gravitational wave detections offer a complementary avenue. However, determining component spins with gravitational waves is hindered by the small dimensionless spins of the NS and the degeneracy in the mass and spin parameters. This degeneracy can be addressed by the inclusion of higher-order modes in the waveform, which are important for systems with unequal masses. This study shows the suitability of neutron star-black hole mergers, which are naturally mass-asymmetric, for precise NS spin measurements. We explore the effects of the black hole masses and spins, higher-mode content, inclination angle, and detector sensitivity on the measurement of NS spin. We find that networks with next-generation observatories like the Cosmic Explorer and the Einstein Telescope can distinguish NS dimensionless spin of 0.04 (0.1) from zero at $1-\sigma$ confidence for events within $\sim 350$ $(\sim 1000)$ Mpc. Networks with A+ and A$^{\sharp}$ detectors achieve similar distinction within $\sim 30$ $(\sim 70)$ Mpc and $\sim 50$ $(\sim 110)$ Mpc, respectively.

It is shown that the foundational axioms of MOND alone predict a strong correlation between a bulk measure of the baryonic surface density, $\Sigma_B$, and the corresponding dynamical one, $\Sigma_D$, of an isolated object, such as a galaxy. The correlation is encapsulated by its high- and low-$\Sigma_B$ behaviors. For $\Sigma_B\gg\Sigma_M\equiv a_0/2\pi G$ ($\Sigma_M$ is the critical MOND surface density) one has $\Sigma_D\approx\Sigma_B$. Their difference -- which would be interpreted as the contribution of dark matter -- is $\Sigma_P=\Sigma_D-\Sigma_B\sim\Sigma_M\ll\Sigma_B$. In the deep-MOND limit, $\Sigma_B\ll\Sigma_M$, one has $\Sigma_D\sim (\Sigma_M\Sigma_B)^{1/2}$. This is a primary prediction of MOND, shared by all theories that embody its basic tenets. Sharper correlations, even strict algebraic relations, $\Sigma_D(\Sigma_B)$, are predicted in specific MOND theories, for specific classes of mass distribution -- e.g., pure discs, or spherical systems -- and for specific definitions of the surface densities. I proceed to discuss such tighter correlations for the central surface densities of axisymmetric galactic systems, $\Sigma^0_B$ and $\Sigma^0_D$. Past work has demonstrated such relations for pure discs in the AQUAL and QUMOND theories. Here I consider them in broader classes of MOND theories. For most observed systems, $\Sigma^0_D$ can not be determined directly at present, but, in many cases, a good proxy for it is the acceleration integral $\mathcal{G}\equiv\int_0^\infty g_r d\ln~r$, where $g_r$ is the radial acceleration along a reflection symmetry of a system, such as a disc galaxy. $\mathcal{G}$ can be determined directly from the rotation curve. I discuss the extent to which $\mathcal{G}$ is a good proxy for $\Sigma^0_D$, and how the relation between them depends on system geometry, from pure discs, through disc-plus-bulge ones, to quasi-spherical systems.

Ellerman bombs (EBs) with significant H$\alpha$ wing emissions and ultraviolet bursts (UV bursts) with strong Si IV emissions are two kinds of small transient brightening events that occur in the low solar atmosphere.We numerically investigated the magnetic reconnection process between the emerging arch magnetic field and the lower atmospheric background magnetic field. We aim to find out if the hot UV emissions and much colder H$\alpha$ wing emissions can both appear in the same reconnection process and how they are located in the reconnection region. The open-source code NIRVANA was applied to perform the 2.5D magnetohydrodynamic (MHD) simulation. We developed the related sub-codes to include the more realistic radiative cooling process for the photosphere and chromosphere and the time-dependent ionization degree of hydrogen. The initial background magnetic field is 600 G, and the emerged magnetic field in the solar atmosphere is of the same magnitude, meaning that it results in a low- $\beta$ magnetic reconnection environment. We also used the radiative transfer code RH1.5D to synthesize the Si IV and H$\alpha$ spectral line profiles based on the MHD simulation results. Magnetic reconnection between emerged and background magnetic fields creates a thin, curved current sheet, which then leads to the formation of plasmoid instability and the nonuniform density distributions. The mix of hot tenuous and much cooler dense plasmas in the turbulent reconnection region can appear at about the same height, or even in the same plasmoid. The turbulent current sheet is always in a dense plasma environment with an optical depth larger than 6.5$\times$10$^{-5}$ due to the emerged magnetic field pushing high-density plasmas upward.

The origin of the quasi-periodic eruptions (QPEs) is possibly mass loss at the periastron of a body moving around the supermassive black hole (SMBH) in a high eccentric orbit. Such a tidally stripped star is expected to radiate gravitational wave thereby leading to shrinkage of the periastron distance, and thus will eventually be disrupted by the SMBH in the previous studies. This scenario predicts a gradually increasing mass transfer, contradicting the observed decay in the intensity of the QPEs in GSN 069. In this paper, we first revisited the orbital evolution of the stripped star. The effect of the mass transfer finally dominates the orbital evolution, resulting in the stripped star finally escaping the SMBH rather than being disrupted by it. Then we suggested a model of a tidal stripped WD moving inside an accretion disk for QPEs in GSN 069. The drag force by the disk can effectively reduce mass transfer and thus explain the observed decay in the intensity of the QPEs in GSN 069. The disk is likely a fallback disk of the tidal disruption event in GSN 069. Considering the evolution of its accretion rate, the increase in the intensity of the latest eruption can also be explained.

The strict periodicity of pulsars is the primary source of information we have to learn about their nature and environment, it allows us to challenge general relativity and measure gravitational waves. Identifying such a periodicity from a discrete set of arrival times is a difficult algorithmic problem, particularly when the pulsar is in a binary system. This challenge is especially acute in $\gamma$-ray pulsar astronomy as there are hundreds of unassociated Fermi-LAT sources awaiting a timing solution that will reveal their nature, and may allow adding them to pulsar timing arrays. The same issue arises when attempting to recover a strict periodicity for repeating fast radio bursts (FRBs). Such a detection would be a major breakthrough, providing us with the FRB source's age, magnetic field, and binary orbit. The problem of recovering a timing solution from sparse time-of-arrival (TOA) data is currently unsolvable for pulsars in binary systems and incredibly hard even for single pulsars. In a series of papers, we will develop an algorithmic set of tools that will allow us to solve the timing recovery problem under different regimes. In this paper, we frame the timing recovery problem as the problem of finding a short vector in a lattice and obtain the solution using off-the-shelf lattice reduction and sieving techniques. As a proof of concept, we solve PSR J0318+0253, a millisecond $\gamma$-ray pulsar discovered by FAST in a $\gamma$-ray directed search, in a few CPU-minutes. We discuss the assumptions of the standard lattice techniques and quantify their performance and limitations.

We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, suggesting that grain growth appears to be significant in disks at earlier protostellar phases. For the rest of the protostars, the dust polarization traces the magnetic field, whose morphology can be approximately classified into three categories: standard-hourglass, rotated-hourglass (with its axis perpendicular to outflow), and spiral-like morphology. 40.0% (+-3.0%) of the protostars exhibit a mean magnetic field direction approximately perpendicular to the outflow on several 100--1000 au scales. However, in the remaining sample, this relative orientation appears to be random, probably due to the complex set of morphologies observed. Furthermore, we classify the protostars into three types based on the C17O (3--2) velocity envelope's gradient: perpendicular to outflow, non-perpendicular to outflow, and unresolved gradient (<1.0~km/s/arcsec). In protostars with a velocity gradient perpendicular to outflow, the magnetic field lines are preferentially perpendicular to outflow, most of them exhibit a rotated hourglass morphology, suggesting that the magnetic field has been overwhelmed by gravity and angular momentum. Spiral-like magnetic fields are associated with envelopes having large velocity gradients, indicating that the rotation motions are strong enough to twist the field lines. All of the protostars with a standard-hourglass field morphology show no significant velocity gradient due to the strong magnetic braking.

We present a uniform and sensitive X-ray census of active galactic nuclei (AGNs) in the two nearest galaxy clusters, Virgo and Fornax, utilizing the newly released X-ray source catalogs from the first all-sky scan of SRG/eROSITA. A total of 50 and 10 X-ray sources are found positionally coincident with the nuclei of member galaxies in Virgo and Fornax, respectively, down to a 0.2--2.3 keV luminosity of $\sim10^{39}\rm~erg~s^{-1}$ and reaching out to a projected distance well beyond the virial radius of both clusters. The majority of the nuclear X-ray sources are newly identified. There is weak evidence that the nuclear X-ray sources are preferentially found in late-type hosts. Several hosts are dwarf galaxies with a stellar mass below $\sim10^{9}\rm~M_\odot$. We find that contamination by non-nuclear X-ray emission can be neglected in most cases, indicating the dominance of a genuine AGN. In the meantime, no nuclear X-ray source exhibits a luminosity higher than a few times $10^{41}\rm~erg~s^{-1}$. The X-ray AGN occupation rate is only $\sim$3\% in both clusters, apparently much lower than that in field galaxies inferred from previous X-ray studies. Both aspects suggest that the cluster environment effectively suppresses AGN activity. The findings of this census have important implications on the interplay between galaxies and their central massive black holes in cluster environments.

We present a catalog of sources detected by the Mikhail Pavlinsky ART-XC telescope onboard the SRG space observatory during the observations of the Galactic plane region near a longitude $l\simeq20$\deg\ (L20 field) in October 2019. The L20 field was observed four times in the scanning mode, which provided a uniform coverage of the sky region with a total area of $\simeq24$ sq. deg with a median sensitivity of $8\times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ (at 50% detection completeness) in the 4$-$12 keV. As a result, we have detected 29 X-ray sources at a statistically significant level, 11 of which have not been detected previously by other X-ray observatories. Preliminary estimates show that four of them can presumably be extragalactic in nature. We also show that the source SRGA J183220.1$-$103508 (CXOGSG J183220.8$-$103510), is most likely a galaxy cluster containing a bright radio galaxy at redshift $z\simeq0.121$.

We present deep total intensity and polarization observations of the Coma cluster at 1.4 and 6.6 GHz performed with the Sardinia Radio Telescope. By combining the single-dish 1.4 GHz data with archival Very Large Array observations we obtain new images of the central radio halo and of the peripheral radio relic where we properly recover the brightness from the large scale structures. At 6.6 GHz we detect both the relic and the central part of the halo in total intensity and polarization. These are the highest frequency images available to date for these radio sources in this galaxy cluster. In the halo, we find a localized spot of polarized signal, with fractional polarization of about 45%. The polarized emission possibly extends along the north-east side of the diffuse emission. The relic is highly polarized, up to 55%, as usually found for these sources. We confirm the halo spectrum is curved, in agreement with previous single-dish results. The spectral index is alpha=1.48 +/- 0.07 at a reference frequency of 1 GHz and varies from alpha ~1.1, at 0.1 GHz, up to alpha ~ 1.8, at 10 GHz. We compare the Coma radio halo surface brightness profile at 1.4 GHz (central brightness and e-folding radius) with the same properties of the other halos, and we find that it has one of the lowest emissivities observed so far. Reanalyzing the relic's spectrum in the light of the new data, we obtain a refined radio Mach number of M=2.9 +/- 0.1.

TOI-677 b is part of an emerging class of ``tidally-detached'' gas giants ($a/R_\star \gtrsim 11$) that exhibit large orbital eccentricities and yet low stellar obliquities. Such sources pose a challenge for models of giant planet formation, which must account for the excitation of high eccentricities without large changes in the orbital inclination. In this work, we present a new Rossiter-McLaughlin (RM) measurement for the tidally-detached warm Jupiter TOI-677 b, obtained using high-precision radial velocity observations from the PFS/Magellan spectrograph. Combined with previously published observations from the ESPRESSO/VLT spectrograph, we derive one of the most precisely constrained sky-projected spin-orbit angle measurements to date for an exoplanet. The combined fit offers a refined set of self-consistent parameters, including a low sky-projected stellar obliquity of $\lambda=3.2^{+1.6}_{-1.5}$ deg and a moderately high eccentricity of $e=0.460^{+0.019}_{-0.018}$, that further constrains the puzzling architecture of this system. We examine several potential scenarios that may have produced the current TOI-677 orbital configuration, ultimately concluding that TOI-677 b most likely had its eccentricity excited through disk-planet interactions. This system adds to a growing population of aligned warm Jupiters on eccentric orbits around hot ($T_{\rm eff}>6100$ K) stars.

The white-light continuum emissions in solar flares (i.e., white-light flares) are usually observed on the solar disk but, in a few cases, off the limb. Here we present on-disk as well as off-limb continuum emissions at 3600 {\AA} (in the Balmer continuum) in an X2.1 flare (SOL2023-03-03T17:52) and an X1.5 flare (SOL2023-08-07T20:46), respectively, observed by the White-light Solar Telescope (WST) on the Advanced Space-based Solar Observatory (ASO-S). These continuum emissions are seen at the ribbons for the X2.1 flare and on loops during the X1.5 event, in which the latter also appears in the decay phase. These emissions also show up in the pseudo-continuum images at Fe I {\lambda}6173 from the Helioseismic and Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO). In addition, the ribbon sources in the X2.1 flare exhibit significant enhancements in the Fe I line at 6569.2 {\AA} and the nearby continuum observed by the Chinese H{\alpha} Solar Explorer (CHASE). It is found that the on-disk continuum emissions in the X2.1 flare are related to a nonthermal electron-beam heating either directly or indirectly, while the off-limb emissions in the X1.5 flare are associated with thermal plasma cooling or due to Thomson scattering. These comprehensive continuum observations can provide good constraints on flare energy deposition models, which helps well understand the physical mechanism of white-light flares.

When low- and intermediate-mass stars evolve off the main sequence, they expand and cool into the red giant stages of evolution, which include those associated with shell H burning (the red giant branch), core He burning (the red clump), and shell He burning (the asymptotic giant branch). The majority of red giants have masses $< 2 M_\odot$, and red giants more massive than this are often excluded from major studies. Here we present a study of the highest-mass stars ($M > 3.0 M_\odot$) in the Kepler sample of 16,000 red giants. We begin by re-estimating their global seismic properties with new light curves, highlighting the differences between using the SAP and PDCSAP light curves provided by Kepler. We use the re-estimated properties to derive new mass estimates for the stars, ending with a final sample of 48 confirmed high-mass stars. We explore their oscillation envelopes, confirming the trends found in recent works such as low mean mode amplitude and wide envelopes. We find, through probabilistic means, that our sample is likely all core He burning stars. We measure their dipole and quadrupole mode visibilities and confirm that the dipole mode visibility tends to decrease with mass.

We study the magnetic and tidal interactions of a gas-giant exoplanet with its host star and with its exomoons, and focus on their retention. We briefly revisit the scaling law for planetary dynamo in terms of its mass, radius and luminosity. Based on the virial theorem, we construct an evolution law for planetary magnetic field and find that its initial entropy is important for the field evolution of a high-mass planet. We estimate the magnetic torques on orbit arising from the star-planet and planet-moon magnetic interactions, and find that it can compensate tidal torques and bypass frequency valleys where dynamical-tidal response is ineffective. For exomoon's retention we consider two situations. In the presence of a circumplanetary disk (CPD), by comparison between CPD's inner and outer radii, we find that planets with too strong magnetic fields or too small distance from its host star tend not to host exomoons. During the subsequent CPD-free evolution, we find, by comparison between planet's spindown and moon's migration timescales, that hot Jupiters with periods of several days are unlikely to retain large exomoons, albeit they could be surrounded by rings from the debris of tidally disrupted moons. In contrast, moons, if formed around warm or cold Jupiters, could have migration timescale is longer than the age of planetary system and be preserved. Finally, we estimate the radio power and flux density due to the star-planet and planet-moon magnetic interactions and give the upper limit of detection distance by FAST.

We present a detailed structural and gas kinematic study of the star-forming complex W5-NW. A cloud-cloud collision scenario unravels with evidences of collision induced star and cluster formation. Various signatures of cloud-cloud collision such as "complementary distribution" and "bridging-features" are explored. At the colliding region, the two clouds have complementary morphologies, where W5-NWb has a filamentary key-like shape which fits into the U-shaped cavity in W5-NWa that behaves like a keyhole. The interaction region between the two clouds is characterised by bridging features with intermediate velocities connecting the two clouds. A skewed V-shaped bridging feature is also detected at the site of collision. A robust picture of the molecular gas distribution highlighting the bridges is seen in the position-position-velocity diagram obtained using the SCOUSEPY algorithm. Star cluster formation with an over-density of Class I and Class II young stellar objects is also seen towards this cloud complex, likely triggered by the cloud collision event.

Star formation quenching is one of the key processes that shape the evolution of galaxies. In this study, we investigate the changes in molecular gas and star formation properties as galaxies transit from the star-forming main sequence to the passive regime. Our analysis reveals that as galaxies move away from the main sequence towards the green valley the radial profile of specific star formation rate surface density ($\Sigma_\mathrm{sSFR}$) is suppressed compared with main sequence galaxies out to a galactocentric radius of 1.5 $R_{e}$ ($\sim$ 7 kpc for our sample). By combining radial profiles of gas fraction ($f_\mathrm{gas}$) and star formation efficiency (SFE), we can discern the underlying mechanism that determines $\Sigma_\mathrm{sSFR}$ at different galactocentric radii. Analysis of relative contributions of $f_\mathrm{gas}$ and SFE to $\Sigma_\mathrm{sSFR}$ uncovers a diverse range of quenching modes. Star formation in approximately half of our quenching galaxies is primarily driven by a single mode (i.e. either $f_\mathrm{gas}$ or SFE), or a combination of both. A collective analysis of all galaxies reveals that the reduction in star formation within the central regions ($R$ $<$ 0.5 $R_{e}$) is primarily attributable to a decrease in SFE. Conversely, in the disk regions ($R$ $>$ 0.5 $R_{e}$), both $f_\mathrm{gas}$ and SFE contribute to the suppression of star formation. Our findings suggest that multiple quenching mechanisms may be at play in our sample galaxies, and even within a single galaxy. We also compare our observational outcomes with those from galaxy simulations and discuss the implications of our data.

Helioseismology has revealed an increase in the rotation rate with depth in a thin ($\sim$30 Mm) near-surface layer. The normalized rotational shear in this layer is independent of latitude. This rotational state is shown to be a consequence of the short characteristic time of near-surface convection compared to the rotation period, and the radial anisotropy of the convective turbulence. Analytical derivations within mean-field hydrodynamics reproduce the observed normalized rotational shear and are in agreement with numerical experiments on the radiative hydrodynamics of solar convection. The near-surface shear layer is the source of the global meridional flow important for the solar dynamo.

We analyze data retrieved by the Imaging Science System onboard the Cassini spacecraft to study the horizontal velocity and vorticity fields of Saturn's Polar Regions (latitudes 60-90$^\circ$N in June-December 2013 and 60-90$^\circ$S in October 2006 and July-December 2008), including the Northern region where the hexagonal wave is prominent. With the aid of an automated two dimensional correlation algorithm we determine two-dimensional maps of zonal and meridional winds, and deduce vorticity maps. We extract zonal averages of zonal winds, providing wind profiles that reach latitudes as high 89.5$^\circ$ in the south and 89.9$^\circ$ in the north. Wind measurements cover the intense polar cyclonic vortices that reach similar peak velocities of 150 ms-1 at 88.5$^\circ$. The hexagonal wave lies in the core of an intense eastward jet at planetocentric latitude 75.8$^\circ$N with motions that become non-zonal at the hexagonal feature. In the south hemisphere the peak of the eastward jet is located at planetocentric latitude 70.4$^\circ$S. A large anticyclone (the South Polar Spot, SPS), similar to the North Polar Spot (NPS) observed at the Voyager times (1980-81), has been observed in images from April 2008 to January 2009 in the South Polar Region at latitude -66.1$^\circ$ close to the eastward jet. The SPS does not apparently excite a wave on the jet. We analyze the stability of the zonal jets, finding potential instabilities at the flanks of the eastward jets around 70$^\circ$ and we measure the eddy wind components, suggesting momentum transfer from eddy motion to the westward jets closer to the poles.

Searching the data of gravitational-wave detectors for signals from compact binary mergers is a computationally demanding task. Recently, machine learning algorithms have been proposed to address current and future challenges. However, the results of these publications often differ greatly due to differing choices in the evaluation procedure. The Machine Learning Gravitational-Wave Search Challenge was organized to resolve these issues and produce a unified framework for machine-learning search evaluation. Six teams submitted contributions, four of which are based on machine learning methods and two are state-of-the-art production analyses. This paper describes the submission from the team TPI FSU Jena and its updated variant. We also apply our algorithm to real O3b data and recover the relevant events of the GWTC-3 catalog.

In the context of modeling spectral energy distributions (SEDs) for blazars, we extend the method that uses a convolutional neural network (CNN) to include external inverse Compton processes. The model assumes that relativistic electrons within the emitting region can interact and up-scatter external photon originating from the accretion disk, the broad-line region, and the torus, to produce the observed high-energy emission. We trained the CNN on a numerical model that accounts for the injection of electrons, their self-consistent cooling, and pair creation-annihilation processes, considering both internal and all external photon fields. Despite the larger number of parameters compared to the synchrotron self-Compton model and the greater diversity in spectral shapes, the CNN enables an accurate computation of the SED for a specified set of parameters. The performance of the CNN is demonstrated by fitting the SED of two flat-spectrum radio quasars, namely 3C 454.3 and CTA 102, and obtaining their parameter posterior distributions. For the first source, the available data in the low-energy band allowed us to constrain the minimum Lorentz factor of the electrons, $\gamma_{\rm min}$, while for the second source, due to the lack of these data, $\gamma_{\rm min} = 10^2$ was set. We used the obtained parameters to investigate the energetics of the system. The model developed here, along with one from B\'egu\'e et al. (2023), enables self-consistent, in-depth modeling of blazar broadband emissions within leptonic scenario.

The latest results of the most detailed analysis of multi-epoch polarization-sensitive observations of active galactic nuclei (AGN) jets at parsecs scales by very long baseline interferometry (VLBI) reveal several characteristic patterns of linear polarization distribution and its variability (Pushkarev et al., 2023; Zobnina et al., 2023). Some of the observed profiles can be reproduced by a simple model of a jet threaded by a helical magnetic field. However, none of the models presented to date can explain the observed polarization profiles with an increase in its degree towards the edges of the jet, and accompanied by a 'fountain' type electric vector pattern and its high temporal variability in the center. Based on simulations of the VLBI observations of relativistic jets, we show here that the observed transverse linear polarization profiles, atypical for the simple magnetic field models can be naturally produced assuming the finite resolution of VLBI arrays and precession of a jet on ten-years scales, observational indications of which are found in an increasing number of AGN. In our simulations, we qualitatively reproduce the distribution of the electric vector and its variability, though the polarization images are characterized by a bright spine due to weak smearing, which is poorly consistent with observations. More effective depolarization can be obtained in models with the suppressed emission of the jet spine.

Bernard Lyot invented the monochromatic birefringent filter in 1933 in order to investigate the coronal emissions of solar structures above the limb with the coronagraph installed at the Pic du Midi observatory. The filter was improved later and he made the first observations of the chromosphere above the solar disk in 1948, at Meudon. After his death, Grenat and Laborde continued the development in the frame of the coming International Geophysical Year (IGY 1957-1958). A modern H$\alpha$ heliograph was completed soon and the flare patrol started in 1956. This instrument was reproduced by two companies (SECASI and OPL) and disseminated around the world in order to contribute to the IGY common effort dedicated to the solar activity survey. We describe in this short paper the capabilities of one of these copies operating at Haute Provence station from 1958 to 1994.

Though efforts to detect them have been made with a variety of methods, no technique can claim a successful, confirmed detection of a moon outside the Solar System yet. Moon detection methods are restricted in capability to detecting moons of masses beyond what formation models would suggest, or they require surface temperatures exceeding what tidal heating simulations allow. We expand upon spectroastrometry, a method that makes use of the variation of the centre of light with wavelength as the result of an unresolved companion, which has previously been shown to be capable of detecting Earth-analogue moons around nearby exo-Jupiters, with the aim to place bounds on the types of moons detectable using this method. We derived a general, analytic expression for the spectroastrometric signal of a moon in any closed Keplerian orbit, as well as a new set of estimates on the noise due to photon noise, pointing inaccuracies, background and instrument noise, and a pixelated detector. This framework was consequently used to derive bounds on the temperature required for Solar System-like moons to be observable around super-Jupiters in nearby systems, with $\epsilon$ Indi Ab as an archetype. We show that such a detection is possible with the ELT for Solar System-like moons of moderate temperatures (150-300 K) in line with existing literature on tidal heating, and that the detection of large (Mars-sized or greater) icy moons of temperatures such as those observed in our Solar System in the very nearest systems may be feasible.

GW150914 marked the start of the gravitational wave (GW) era with the direct detection of binary black hole (BBH) merger by the LIGO-Virgo GW detectors. The event was temporally coincident with a weak signal detected by Fermi-GBM, which hinted towards the possibility of electromagnetic emission associated with the compact object coalescence. The detection of a short Gamma-Ray Burst (GRB) associated with GW1708017, along with several multi-wavelength detections, truly established that compact object mergers are indeed multi-messenger events. The Cadmium Zinc Telluride Imager onboard AstroSat can search for X-ray counterparts of the GW events and has detected over 600 GRBs since launch. Here we present results from our searches for counterparts coincident with GW triggers from the first three LIGO-Virgo-KAGRA (LVK) GW Transient Catalogs. For 72 out of 90 GW events for which AstroSat-CZTI data was available, we undertook a systematic search for temporally coincident transients and we detected no X-ray counterparts. We evaluate the upper limits on the maximum possible flux from the source in a 100 s window centered around each trigger, consistent with the GW localization of the event. Thanks to the high sensitivity of CZTI, these upper limits are highly competitive with those from other spacecraft. We use these upper limits to constrain the theoretical models that predict high-energy counterparts to the BBH mergers. We also discuss the probability of non-detections of BBH mergers at different luminosities and the implications of such non-detections from the ongoing fourth observing run of the LVK detectors.

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase.

We present the discovery of `A Big Ring on the Sky' (BR), the second ultra-large-scale structure (uLSS) found in MgII-absorber catalogues, following the previously reported Giant Arc (GA). In cosmological terms the BR is close to the GA - at the same redshift $z \sim 0.8$ and with a separation on the sky of only $\sim 12^\circ$. Two extraordinary uLSSs in such close configuration raises the possibility that together they form an even more extraordinary cosmological system. The BR is a striking circular, annulus-like, structure of diameter $\sim 400$ Mpc (proper size, present epoch). The method of discovery is as described in the GA paper, but here using the new MgII-absorber catalogues restricted to DR16Q quasars. Using the Convex Hull of Member Spheres (CHMS) algorithm, we estimate that the annulus and inner absorbers of the BR have departures from random expectations, at the density of the control field, of up to $5.2\sigma$. We present the discovery of the BR, assess its significance using the CHMS, Minimal Spanning Tree (MST), FilFinder and Cuzick & Edwards (CE) methods, show it in the context of the GA+BR system, and suggest some implications for the origins of uLSS and for our understanding of cosmology. For example, it may be that unusual geometric patterns, such as these uLSSs, have an origin in cosmic strings.

The large-scale structure in cosmology is highly non-Gaussian at late times and small length scales, making it difficult to describe analytically. Parameter inference, data reconstruction, and data generation tasks in cosmology are greatly aided by various machine learning models. In order to retain as much information as possible while solving these problems, this work operates at the field level, rather than at the level of summary statistics. The probability density function of the large-scale structure is learned with normalizing flows, a class of probabilistic generative models. Normalizing flows learn the transformation from a simple base distribution to a more complicated distribution, much like the matter content evolved to its present day complexities from a Gaussian field at early times. While the normalizing flows have accurately modelled 2-dimensional projections of the matter content, we find that denoising diffusion models are well-suited for volumetric data. A super-resolution emulator is developed for cosmological simulation volumes, generating high-resolution baryonic simulation volumes conditional on low-resolution dark matter simulations. The super-resolution emulator is trained to perform outpainting, and can thus upgrade very large cosmological volumes from low-resolution to high-resolution using an iterative outpainting procedure.

Any astrophysical object can, in principle, serve as a probe of the interaction between Dark Matter and regular, baryonic matter. This method is based on the potential observable consequences annihilations of captured Dark Matter has on the surface temperature of the object itself. In a series of previous papers we developed and validated simple analytic approximations for the total capture rates of Dark Matter (DM) valid in four distinct regions of the DM-nucleon scattering cross section ($\sigma$) vs. DM particle mass ($m_X$) parameter space. In this work we summarize those previous results and extend them significantly, by deriving a completely general, closed form solution for the total capture rate of Dark Matter in the multiscatter regime. Moreover, we demonstrate the existence of a region in the $\sigma$ vs. $m_X$ parameter space where the constraining power of any astrophysical object heated by annihilations of captured DM is lost. This corresponds to a maximal temperature ($T_{crit}$) any astrophysical object can have, such that it can still serve as a DM probe. Any object with observed temperature $T_{obs}>T_{crit}$ loses its DM constraining power. We provide analytic formulae that can be used to estimate $T_{crit}$ for any object.

Analytical models of galaxy-halo connection such as the Halo Occupation Distribution (HOD) model have been widely used over the past decades as a means to intensively test perturbative models on quasi-linear scales. However, these models fail to reproduce the galaxy-galaxy lensing signal on non-linear scales, over-predicting the observed signal up to 40%. With ongoing Stage-IV galaxy surveys such as DESI and EUCLID, it is now crucial to accurately model the galaxy-halo connection up to intra-halo scales to accurately estimate theoretical uncertainties of perturbative models. This paper compares the standard HOD model to an extended HOD framework that incorporates as additional features galaxy assembly bias and local environmental dependencies on halo occupation. These models have been calibrated against the observed clustering and galaxy-galaxy lensing signal of eBOSS Luminous Red Galaxies (LRG) and Emission Lines Galaxies (ELG) in the range 0.6 < z < 1.1. A combined clustering-lensing cosmological analysis is then performed on the simulated galaxy samples of both standard and extended HOD frameworks to quantify the systematic budget of perturbative models. The extended HOD model offers a more comprehensive understanding of the connection between galaxies and their surroundings. In particular, we found that the LRGs preferentially occupy denser and more anisotropic environments. Our results highlight the importance of considering environmental factors in galaxy formation models, with an extended HOD framework that reproduces the observed signal within 20% on scales below 10 Mpc/h. Our cosmological analysis reveals that our perturbative model yields similar constraints regardless of the galaxy population, with a better goodness of fit for the extended HOD. These results suggest that the extended HOD framework should be used to quantify modeling systematics.

We have successfully integrated $\Lambda_{\rm s}$CDM, a promising model for alleviating cosmological tensions, into a theoretical framework by endowing it with a specific Lagrangian from the VCDM model, a type-II minimally modified gravity. In this theory, we demonstrate that an auxiliary scalar field with a linear potential induces an effective cosmological constant, enabling the realization of an abrupt mirror AdS-dS transition in the late universe through a piecewise linear potential. To eliminate the sudden singularity in this setup and ensure stable transitions, we smooth out this potential. Realized within the VCDM theory, the $\Lambda_{\rm s}$CDM model facilitates two types of rapid smooth mirror AdS-dS transitions: (i) the agitated transition, associated with a smooth jump in the potential, where $\Lambda_{\rm s}$, and consequently $H$, exhibits a bump around the transition's midpoint; and (ii) the quiescent transition, associated with a smooth change in the slope of the potential, where $\Lambda_{\rm s}$ transitions gracefully. These transitions are likely to leave distinct imprints on the background and perturbation dynamics, potentially allowing the observational data to distinguish between them. This novel theoretical framework propels $\Lambda_{\rm s}$CDM into a fully predictive model capable of exploring the evolution of the Universe including the late-time AdS-dS transition epoch, and extends the applicability of the model. We believe further research is crucial in establishing $\Lambda_{\rm s}$CDM as a leading candidate or guide for a new concordance cosmological model.

The high-mass X-ray binary LS I +61{\deg}303 is composed of a Be-type star and a compact object. The emission is variable and periodic across the electromagnetic spectrum, from radio to very high-energy gamma rays. The orbital period is ~26.5 d, and the source also features a super-orbital period with a value of ~4.6 years. Long-term monitoring of the binary by the Owens Valley Radio Observatory (OVRO) at 15 GHz has now completed 13.8 years, which corresponds to three full cycles of the super-orbital period. We performed a timing analysis of the OVRO radio light curve, and we also combined the OVRO data with the full archive of previous radio observations and computed the discrete autocorrelation function. The most powerful features in the periodogram of the OVRO data are two peaks at P1 = 26.49 +/- 0.05 d and P2 = 26.93 +/- 0.05 d. Our measurement of the long-term period is P_long = 1698 +/- 196 d. Dividing the OVRO data into three segments of equal length showed that the two periods, P1 and P2, are present in the periodogram of each of the consecutive long-term cycles. Our analysis of the full radio archive resulted in the detection of the same three periods, and the autocorrelation function showed a regular pattern, proving the stability of the super-orbital modulation. We report a possible systematic modulation of the radio flux density with a timescale of approximately 40 years that has so far remained unnoticed. The physical model of a relativistic jet whose mass loading is modulated with the orbital period P1 and is precessing with the slightly larger period P2, giving rise to a beating with period P_long, had previously been able to reproduce the radio and gigaelectron volt emission. The ongoing presence and the stability of the periodic signals imply that this model is still the most plausible explanation for the physical processes at work in this source. (abridged)

The ability to measure precise and accurate stellar effective temperatures ($T_{\rm{eff}}$) and surface gravities ($\log(g)$) is essential in determining accurate and precise abundances of chemical elements in stars. Measuring $\log(g)$ from isochrones fitted to colour-magnitude diagrams of open clusters is significantly more accurate and precise compared to spectroscopic $\log(g)$. By determining the ranges of ages, metallicity, and extinction of isochrones that fit the colour-magnitude diagram, we constructed a joint probability distribution of $T_{\rm{eff}}$ and $\log(g)$. The joint photometric probability shows the complex correlations between $T_{\rm{eff}}$ and $\log(g)$, which depend on the evolutionary stage of the star. We show that by using this photometric prior while fitting spectra, we can acquire more precise spectroscopic stellar parameters and abundances of chemical elements. This reveals higher-order abundance trends in open clusters like traces of atomic diffusion. We used photometry and astrometry provided by the \textit{Gaia} DR3 catalogue, Padova isochrones, and Galactic Archaeology with HERMES (GALAH) DR4 spectra. We analysed the spectra of 1979 stars in nine open clusters, using MCMC to fit the spectroscopic abundances of 26 elements, $T_{\rm{eff}}$, $\log(g)$, $v_{\rm{mic}}$, and $v_{\rm{broad}}$. We found that using photometric priors improves the accuracy of abundances and $\log(g)$, which enables us to view higher-order trends of abundances caused by atomic diffusion in M67 and Ruprecht 147.

We present a statistical and multiwavelength photometric studies of young open cluster IC 1590. We identified 91 cluster members using $Gaia$ DR3 astrometry data using ensemble-based unsupervised machine learning techniques. From $Gaia$ EDR3 data, we estimate the best-fitted parameters for IC 1590 using the Automated Stellar Cluster Analysis package (ASteCA) yielding the distance $d$ $\sim$ 2.87 $\pm$ 0.02 Kpc, age $\sim$ 3.54 $\pm$ 0.05 Myr, metallicity $z$ $\sim$ 0.0212 $\pm$ 0.003, binarity value of $\sim$ 0.558, and extinction $A_v$ $\sim$ 1.252 $\pm$ 0.4 mag for an $R_v$ value of $\sim$ 3.322 $\pm$ 0.23. We estimate the initial mass function slope of the cluster to be $\alpha$ = 1.081 $\pm$ 0.112 for single stars and $\alpha$ = 1.490 $\pm$ 0.051 for a binary fraction of $\sim$ 0.558 in the mass range 1 M$_{\odot}$ $\leq$ m(M$_{\odot}$) $\leq$ 100 M$_{\odot}$. The $G$-band luminosity function slope is estimated to be $\sim$ 0.33 $\pm$ 0.09. We use $(J-H)$ versus $(H-K_s)$ color-color diagram to identify young stellar objects (YSOs). We found that all the identified YSOs have ages $\leq$ 2 Myr and masses $\sim$ 0.35 - 5.5 M$_{\odot}$. We also fit the radial surface density profile. Using the galpy we performed orbit analysis of the cluster. The extinction map for the cluster region has been generated using the PNICER technique and it is almost similar to the dust structure obtained the dust structure from the 500 $\mu$$m$ dust continuum emissions map of $Herschel$ SPIRE. We finally at the end discussed the star formation scenario in the cluster region.

While contact binary objects are common in the solar system, their formation mechanism is unclear. In this work we examine several contact binaries and calculate the necessary strength parameters that allow the two lobes to merge without the smaller of the two being gravitationally destroyed by the larger. We find a small but non-zero amount of cohesion or a large friction angle is required for the smaller lobe to survive the merging process, consistent with observations. This means it is possible for two previously separated rubble piles to experience a collapse of their mutual orbit and form a contact binary. The necessary strength required to survive this merger depends on the relative size, shape, and density of the body, with prolate shapes requiring more cohesion than oblate shapes.

A new meteor shower $\lambda$-Sculptorids produced by the comet 46P/Wirtanen was forecast for December 12, 2023. The predicted activity was highly uncertain, but generally considered to be low. Observations in Australia, New Zealand, and Oceania were solicited to help constrain the size distribution of meteoroids in the shower. This work aims to characterize the new meteor shower, by comparing the observed and predicted radiants and orbits, and to provide a calibration for future predictions. Global Meteor Network video cameras were used to observe the meteor shower. Multi-station observations were used to compute trajectories and orbits, while single-station observations were used to measure the flux profile. A total of 23 $\lambda$-Sculptorid orbits have been measured. The shower peaked at a zenithal hourly rate (ZHR) of $0.65^{+0.24}_{-0.20}$ meteors per hour at $\lambda_{\odot} = 259.988^{\circ} \pm 0.042^{\circ}$. Due to the low in-atmosphere speed of 15~km s$^{-1}$, the mean mass of observed meteoroids was 0.5~g ($\sim10$~mm diameter), an order of magnitude higher than predicted. The dynamical simulations of the meteoroid stream can only produce such large meteoroids arriving at Earth in 2023 with correct radiants when a very low meteoroid density of $\sim 100$~kg~m$^{-3}$ is assumed. However, this assumption cannot reproduce the activity profile. It may be reproduced by considering higher density meteoroids in a larger ecliptic plane-crossing time window ($\Delta T$ = 20 days) and trails ejected prior to 1908, but then the observed radiant structure is not reproduced.

In this paper, we focus on the scientific case of Cygnus OB2, a northern sky young massive stellar cluster (YMSC) located towards the Cygnus X star-forming complex. We consider a model that assumes cosmic ray acceleration occurring only at the termination shock of the collective wind of the YMSC and address the question of whether, and under what hypotheses, hadronic emission by the accelerated particles can account for the observations of Cygnus OB2 obtained by Fermi-LAT, HAWC and LHAASO. In order to do so, we carefully review the available information on this source, also confronting different estimates of the relevant parameters with ad hoc developed simulations. Once other model parameters are fixed, the spectral and spatial properties of the emission are found to be very sensitive to the unknown properties of the turbulent magnetic field. Comparison with the data shows that our suggested scenario is incompatible with Kolmogorov turbulence. Assuming Kraichnan or Bohm type turbulence spectra, the model accounts well for the Very High Energy (VHE) data, but fails to reproduce the centrally peaked morphology observed by Fermi-LAT, suggesting that additional effects might be important for lower energy $\gamma$-ray emission. We discuss how additional progress can be made with a more detailed and extended knowledge of the spectral and morphological properties of the emission.

We present a study of the detection efficiency for the TESS mission, focusing on the yield of longer-period transiting exoplanets ($P > 25$ days). We created the Transit Investigation and Recoverability Application (TIaRA) pipeline to use real TESS data with injected transits to create sensitivity maps which we combine with occurrence rates derived from Kepler. This allows us to predict longer-period exoplanet yields, which will help design follow-up photometric and spectroscopic programs, such as the NGTS Monotransit Program. For the TESS Year 1 and Year 3 SPOC FFI lightucurves, we find $2271^{+241}_{-138}$ exoplanets should be detectable around AFGKM dwarf host stars. We find $215^{+37}_{-23}$ exoplanets should be detected from single-transit events or "monotransits". An additional $113^{+22}_{-13}$ detections should result from "biennial duotransit" events with one transit in Year 1 and a second in Year 3. We also find that K dwarf stars yield the most detections by TESS per star observed. When comparing our results to the TOI catalogue we find our predictions agree within $1\sigma$ of the number of discovered systems with periods between 0.78 and 6.25\,days and agree to $2\sigma$ for periods between 6.25 and 25\,days. Beyond periods of 25 days we predict $403^{+64}_{-38}$ detections, which is 3 times as many detections as there are in the TOI catalogue with $>3\sigma$ confidence. This indicates a significant number of long-period planets yet to be discovered from \tess\ data as monotransits or biennial duotransits.

One common approach for solving collisions between protoplanets in simulations of planet formation is to employ analytical scaling laws. The most widely used one was developed by Leinhardt & Stewart (2012) from a catalog of ~ 180 N-body simulations of rubble-pile collisions. In this work, we use a new catalogue of more than 20,000 SPH simulations to test the validity and the prediction capability of Leinhardt & Stewart (2012) scaling laws. We find that these laws overestimate the fragmentation efficiency in the merging regime and they are not able to properly reproduce the collision outcomes in the super-catastrophic regime. In the merging regime, we also notice a significant dependence between the collision outcome, in terms of the largest remnant mass, and the relative mass of the colliding protoplanets. Here, we present a new set of scaling laws that are able to better predict the collision outcome in all regimes and it is also able to reproduce the observed dependence on the mass ratio. We compare our new scaling laws against a machine learning approach and obtain similar prediction efficiency.

Neptune's external mean motion resonances play an important role in sculpting the observed population of transneptunian objects (TNOs). The population of scattering TNOs are known to 'stick' to Neptune's resonances while evolving in semimajor axis ($a$), though simulations show that resonance sticking is less prevalent at $a\gtrsim200-250$ au. Here we present an extensive numerical exploration of the strengths of Neptune's resonances for scattering TNOs with perihelion distances $q=33$ au. We show that the drop-off in resonance sticking for the large $a$ scattering TNOs is not a generic feature of scattering dynamics, but can instead be attributed to the specific configuration of Neptune and Uranus in our solar system. In simulations with just Uranus removed from the giant planet system, Neptune's resonances are strong in the scattering population out to at least $\sim300$ au. Uranus and Neptune are near a 2:1 period ratio, and the variations in Neptune's orbit resulting from this near resonance are responsible for destabilizing Neptune's resonances for high-$e$ TNO orbits beyond the $\sim20$:1 resonance at $a\approx220$ au. Direct interactions between Uranus and the scattering population are responsible for slightly weakening Neptune's closer-in resonances. In simulations where Neptune and Uranus are placed in their mutual 2:1 resonance, we see almost no stable libration of scattering particles in Neptune's external resonances. Our results have important implications for how the strengths of Neptune's distant resonances varied during the epoch of planet migration when the Neptune-Uranus period ratio was evolving. These strength variations likely affected the distant resonant and detached TNO populations.

We present the optical-near infrared spectral energy distributions (SED) and near infrared variability properties of 30 low-redshift iron low-ionization Broad Absorption Line quasars (FeLoBALQs) and matched samples of LoBALQs and unabsorbed quasars. Significant correlations between the SED properties and accretion rate indicators found among the unabsorbed comparison sample objects suggest an intrinsic origin for SED differences. A range of reddening likely mutes these correlations among the FeLoBAL quasars. The restframe optical-band reddening is correlated with the location of the outflow, suggesting a link between the outflows and the presence of dust. We analyzed WISE variability and provide a correction for photometry uncertainties in an appendix. We found an anticorrelation between the variability amplitude and inferred continuum emission region size, and suggest that as the origin of the anticorrelation between variability amplitude and luminosity typically observed in quasars. We found that the LoBALQ optical emission line and other parameters are more similar to those of the unabsorbed continuum sample objects than the FeLoBALQs. Thus, FeLoBAL quasars are a special population of objects. We interpret the results using an accretion-rate scenario for FeLoBAL quasars. The high accretion rate FeLoBAL quasars are radiating powerfully enough to drive a thick, high-velocity outflow. Quasars with intermediate accretion rates may have an outflow, but it is not sufficiently thick to include FeII absorption. Low accretion rate FeLoBAL outflows originate in absorption in a failing torus, no longer optically thick enough to reprocess radiation into the near-IR.

We report the characterization of two planet candidates detected by the Transiting Exoplanet Survey Satellite (TESS), TOI-1199\:b and TOI-1273\:b, with periods of 3.7 and 4.6\,days, respectively. Follow-up observations for both targets, which include several ground-based light curves, confirmed the transit events. High-precision radial velocities from the SOPHIE spectrograph revealed signals at the expected frequencies and phases of the transiting candidates and allowed mass determinations with a precision of 8.4\% and 6.7\% for TOI-1199\:b and TOI-1273\:b, respectively. The planetary and orbital parameters were derived from a joint analysis of the radial velocities and photometric data. We find that the planets have masses of 0.239$\,\pm\,$0.020\,M$_{\mathrm{J}}$ \ and 0.222$\,\pm\,$0.015\,M$_{\mathrm{J}}$ \ and radii of 0.938$\,\pm\,$0.025\,R$_{\mathrm{J}}$ \ and 0.99$\,\pm\,$0.22\,R$_{\mathrm{J}}$,\ respectively. The grazing transit of TOI-1273\:b translates to a larger uncertainty in its radius, and hence also in its bulk density, compared to TOI-1199\:b. The inferred bulk densities of 0.358$\,\pm\,$0.041\,g\,cm$^{-3}$ \ and 0.28$\,\pm\,$0.11\,g\,cm$^{-3}$ \ are among the lowest known for exoplanets in this mass range, which, considering the brightness of the host stars ($V$$\approx$11\,mag), render them particularly amenable to atmospheric characterization via the transit spectroscopy technique. The better constraints on the parameters of TOI-1199\:b provide a transmission spectroscopy metric of 134\,$\pm$\,17, making it the better suited of the two planets for atmospheric studies.

We present the analysis of archival Very Large Telescope (VLT) MUSE and multi-epoch Hubble Space Telescope (HST) WFPC2 and WFC3/UVIS narrow-band observations of the remnant associated with the ejecta of the mid-nineteenth century outburst of the recurrent nova T Pyx. These data sets are used to investigate its true 3D physical structure and the nebular expansion patterns along the line of sight and on the plane of the sky. The VLT MUSE emission line maps and 3D visualisations based on position-position-velocity diagrams reveal T Pyx as a bipolar nebula, with a knotty toroidal structure as its waist best seen in H$\beta$ and two open bowl-shaped bipolar lobes (a diabolo) best revealed by the [OIII] emission lines. The comparison of multi-epoch HST WFPC2 and WFC3/UVIS narrow-band images and VLT MUSE emission line maps of T Pyx reveals the angular expansion of the remnant through the proper motion of individual knots and nebular features. The angular expansion is confirmed to be homologous in the period from 1994.2 to 2007.4 before the recent 2011 outburst, but there is suggestive evidence that the inner knots have experienced a higher expansion rate since then.

Magnetic storms on stars manifest as remarkable, randomly occurring changes of the luminosity over durations that are tiny in comparison to the normal evolution of stars. These stellar flares are bursts of electromagnetic radiation from X-ray to radio wavelengths, and they occur on most stars with outer convection zones. They are analogous to the events on the Sun known as solar flares, which impact our everyday life and modern technological society. Stellar flares, however, can attain much greater energies than those on the Sun. Despite this, we think that these phenomena are rather similar in origin to solar flares, which result from a catastrophic conversion of latent magnetic field energy into atmospheric heating within a region that is relatively small in comparison to normal stellar sizes. We review the last several decades of stellar flare research. We summarize multi-wavelength observational results and the associated thermal and nonthermal processes in flaring stellar atmospheres. Static and hydrodynamic models are reviewed with an emphasis on recent progress in radiation-hydrodynamics and the physical diagnostics in flare spectra. Thanks to their effects on the space weather of exoplanetary systems (and thus in our search for life elsewhere in the universe) and their preponderance in \emph{Kepler} mission data, white-light stellar flares have re-emerged in the last decade as a widely-impactful area of study within astrophysics. Yet, there is still much we do not understand, both empirically and theoretically, about the spectrum of flare radiation, its origin, and its time evolution. We conclude with several big-picture questions that are fundamental in our pursuit toward a greater understanding of these enigmatic stellar phonemena and, by extension, those on the Sun.

Detecting imprints of orbital eccentricity in gravitational wave signals promises to shed light on the formation mechanisms of binary black holes. To constrain formation mechanisms, distributions of eccentricity derived from numerical simulations of astrophysical formation channels are compared to the estimates of eccentricity inferred from GW signals. We report that the definition of eccentricity typically used in astrophysical simulations is inconsistent with the one used while modeling GW signals, with the differences mainly arising due to the choice of reference frequency used in both cases. We also posit a prescription for calculating eccentricity from astrophysical simulations by evolving ordinary differential equations obtained from post-Newtonian theory, and using the dominant ($\ell = m =2$) mode's frequency as the reference frequency; this ensures consistency in the definitions. On comparing the existing eccentricities of binaries present in the Cluster Monte Carlo catalog of globular cluster simulations with the eccentricities calculated using the prescription presented here, we find a significant discrepancy at $e \gtrsim 0.2$; this discrepancy becomes worse with increasing eccentricity. We note the implications this discrepancy has for existing studies, and recommend that care be taken when comparing data-driven constraints on eccentricity to expectations from astrophysical formation channels.

Localizing fast radio bursts (FRBs) to their host galaxies is an essential step to better understanding their origins and using them as cosmic probes. The CHIME/FRB Outrigger program aims to add VLBI-localization capabilities to CHIME, such that FRBs may be localized to tens of milliarcsecond precision at the time of their discovery, more than sufficient for host galaxy identification. The first-built outrigger telescope is KKO, located 66 kilometers west of CHIME. Cross-correlating KKO with CHIME can achieve arcsecond-scale localization in right ascension while avoiding the worst effects of the ionosphere. This paper presents measurements of KKO's performance throughout its commissioning phase, as well as a summary of its design and function. We demonstrate KKO's capabilities as a standalone instrument by producing full-sky images, mapping the angular and frequency structure of the primary beam, and measuring feed positions. To demonstrate the localization capabilities of the CHIME -- KKO baseline, we collected five separate observations each for a set of twenty bright pulsars, and aimed to measure their positions to within 5~arcseconds. All of these pulses were successfully localized to within this specification. The next two outriggers are expected to be commissioned in 2024, and will enable subarcsecond localizations for approximately hundreds of FRBs each year.

Data from the Kepler space telescope have led to the discovery of thousands of planet candidates. Most of these candidates are likely to be real exoplanets, but a significant number of false positives still contaminate the sample, especially in candidate lists from the K2 mission. Identifying and rejecting the false positives lurking in the planet candidates sample is important for prioritizing follow-up resources and measuring the most accurate population statistics. Here, we identify false positives in the K2 planet candidate sample using a technique called "ephemeris matching," in which we compare the period and transit time of different signals. When signals from different stars show the same period and time of transit, we can conclude that at least one of the two signals is contamination. We identify 42 false positives among published K2 planet candidates (nearly 2% of the complete list), one of which (K2-256 b) was previously validated as genuine exoplanet. This work increases the reliability of the K2 planet sample and helps boost confidence in the surviving planet candidates.

We report simulation studies of six low-energy electron-antineutrino detector designs, with the goal of determining their ability to resolve the direction to an antineutrino source. Such detectors with target masses on the one-ton scale are well-suited to reactor monitoring at distances of 5--25 meters from the core. They can provide accurate measurements of reactor operating power, fuel mix, and burnup, as well as unsurpassed nuclear non-proliferation information in a non-contact cooperating reactor scenario such as those used by IAEA. A number of groups around the world are working on programs to develop detectors similar to some of those in this study. Here, we examine and compare several approaches to detector geometry for their ability not only to detect the inverse beta decay (IBD) reaction, but also to determine the source direction of incident antineutrinos. The information from these detectors provides insight into reactor power and burning profile, which is especially useful in constraining the clandestine production of weapons material. In a live deployment, a non-proliferation detector must be able to isolate the subject reactor, possibly from a field of much-larger power reactors; directional sensitivity can help greatly with this task. We also discuss implications for using such detectors in longer-distance observation of reactors, from a few km to hundreds of km. We have modeled six abstracted detector designs, including two for which we have operational data for validating our computer modeling and analytical processes. We have found that the most promising options, regardless of scale and range, have angular resolutions on the order of a few degrees, which is better than any yet achieved in practice by a factor of at least two.

As the authors of the original paper "Detailed study of the astrophysical direct capture reaction $^{6}{\rm Li}(p, \gamma)^{7}{\rm Be}$ in a potential model approach" published in Ref. [1] we confirm that there are several points in the Comments which should be taken into the consideration.

Ultralight dark matter (ULDM) particles of mass $m_\phi \lesssim 1~{\rm eV}$ can form boson stars in DM halos. Collapse of boson stars leads to explosive bosenova emission of copious relativistic ULDM particles. In this work, we analyze sensitivity of terrestrial and space-based experiments to detect such relativistic scalar ULDM particles interacting through quadratic couplings with Standard Model constituents, including electrons, photons and gluons. We highlight key differences with searches for linear ULDM couplings. Screening of ULDM with quadratic couplings near the surface of the Earth can significantly impact observations in terrestrial experiments, motivating future space-based experiments. We demonstrate excellent ULDM discovery prospects, especially for quantum sensors, which can probe quadratic couplings orders below existing constraints by detecting bosenova events in the ULDM mass range $10^{-23}\,{\rm eV} \lesssim m_\phi \lesssim 10^{-5}\,{\rm eV}$. We also report updated constraints on quadratic couplings of ULDM in case it comprises cold DM.

Angular momentum and spin precession are expected to be generic features of a significant fraction of binary black hole systems. As such, it is essential to have waveform models that faithfully incorporate the effects of precession. Here, we assess how well the current state of the art models achieve this for waveform strains constructed only from the $\ell=2$ multipoles.Specifically, we conduct a survey on the faithfulness of the waveform models %(approximants) \texttt{SEOBNRv5PHM}, \texttt{TEOBResumS}, \texttt{IMRPhenomTPHM}, \texttt{IMRPhenomXPHM} to the numerical relativity (NR) surrogate \texttt{NRSur7dq4} and to NR waveforms from the \texttt{SXS} catalog. The former assessment involves systems with mass ratios up to six and dimensionless spins up to 0.8. The latter employs $317$ short and $23$ long \texttt{SXS} waveforms. For all cases, we use reference inclinations of zero and $90^\circ$. We find that all four models become more faithful as the mass ratio approaches unity and when the merger-ringdown portion of the waveforms are excluded. We also uncover a correlation between the co-precessing $(2,\pm2)$ multipole mismatches and the overall strain mismatch. We additionally find that for high inclinations, precessing $(2,\pm 1)$ multipoles that are more faithful than their $(2,\pm2)$ counterparts, and comparable in magnitude, improve waveform faithfulness. As a side note, we show that use of uniformly-filled parameter spaces may lead to an overestimation of precessing model faithfulness. We conclude our survey with a parameter estimation study in which we inject two precessing \texttt{SXS} waveforms (at low and high masses) and recover the signal with \texttt{SEOBNRv5PHM}, \texttt{IMRPhenomTPHM} and \texttt{IMRPhenomXPHM}. As a bonus, we present preliminary multidimensional fits to model unfaithfulness for Bayesian model selection in parameter estimation studies.

Future gravitational-wave detectors, especially the Laser Interferometer Space Antenna (LISA), will be sensitive to black hole binaries formed in astrophysical environments that promote large eccentricities and spin precession. Gravitational-wave templates that include both effects have only recently begun to be developed. The Efficient Fully Precessing Eccentric (EFPE) family is one such model, covering the inspiral stage with small-eccentricity-expanded gravitational-wave amplitudes accurate for eccentricities $e < 0.3$. In this work, we extend this model to cover a larger range of eccentricities. The new EFPE_ME model is able to accurately represent the leading-order gravitational-wave amplitudes to $e \leq 0.8$. Comparing the EFPE and the EFPE_ME models in the LISA band, however, reveals that there is no significant difference when $e_0 \leq 0.5$ for binaries at 4 years before merger, as radiation reaction circularizes supermassive black hole binaries too quickly. This suggests that the EFPE model may have a larger regime of validity in eccentricity space than previously thought, making it suitable for some inspiral parameter estimation with LISA data. On the other hand, for systems with $e_0 > 0.5$, the deviations between the models are significant, particularly for binaries with total masses below $10^5\, \mathrm{M}_{\odot}$. This suggests that the EFPE_ME model will be crucial to avoid systematic bias in parameter estimation with LISA in the future, once this model has been hybridized to include the merger and ringdown.

We combine the non-relativistic effective theory of dark matter (DM) - electron interactions with linear response theory to obtain a formalism that fully accounts for screening and collective excitations in DM-induced electronic transition rate calculations for general DM-electron interactions. In the same way that the response of a dielectric material to an external electric field in electrodynamics is described by the dielectric function, so in our formalism the response of a detector material to a DM perturbation is described by a set of generalised susceptibilities which can be directly related to densities and currents arising from the non-relativistic expansion of the Dirac Hamiltonian. We apply our formalism to assess the sensitivity of non-spin-polarised detectors, and find that in-medium effects significantly affect the experimental sensitivity if DM couples to the detector's electron density, while being decoupled from other densities and currents. Our formalism can be straightforwardly extended to the case of spin-polarised materials.

Upgrades beyond the current second generation of ground-based gravitational wave detectors will allow them to observe tens of thousands neutron star and black hole binaries. Given the typical minute-to-hour duration of neutron star signals in the detector frequency band, a number of them will overlap in the time-frequency plane resulting in a nonzero cross-correlation. We examine source confusion arising from overlapping signals whose time-frequency tracks cross. Adopting the median observed merger rate of $100$ Gpc$^{-3}$yr$^{-1}$, each neutron star binary signal overlaps with an average of 42(4)[0.5] other signals when observed from 2(5)[10] Hz. The vast majority of overlaps occur at low frequencies where the inspiral evolution is slow: 91% of time-frequency overlaps occur in band below 5 Hz. The combined effect of overlapping signals does not satisfy the central limit theorem and source confusion cannot be treated as stationary, Gaussian noise: on average 0.91(0.17)[0.05] signals are present in a single adaptive time-frequency bin centered at 2(5)[10] Hz. We quantify source confusion under a realistic neutron star binary population and find that parameter uncertainty typically increases by less than 1% unless there are overlapping signals whose detector-frame chirp mass difference is $\lesssim 0.01 M_{\odot}$ and the overlap frequency is $\gtrsim$ 40 Hz. Out of $1\times10^6$ simulated signals, 0.14% fall within this region of detector-frame chirp mass differences, but their overlap frequencies are typically lower than 40 Hz. Source confusion for ground-based detectors is significantly milder than the equivalent LISA problem.

In neutron stars, flavor-changing weak interactions determine the equilibrium fraction of protons over neutrons. In binary neutron-star mergers, violent changes in density modify this equilibrium value at timescales of milliseconds, comparable to those required for weak interactions to take place. As a result, the fraction of protons evolves out of phase with the density oscillations, giving rise to irreversible processes. The corresponding shift in pressure leads to dissipative work that can be modeled as an effective bulk-viscous correction. In this work, we derive the relevant equations of motion of Israel-Stewart hydrodynamics within this context. Using a toy model, we compute the second-order transport coefficients. Finally, we comment on the use of a realistic equation of state. Our results are expected to be useful for the study of viscous effects in numerical simulations of binary mergers.

An essay on Horndeski gravity, how it was formulated in the early 1970s and how it was 're-discovered' and widely adopted by Cosmologists more than thirty years later.

We study the fluctuation-dissipation relation for sound waves in the cosmic microwave background (CMB), employing effective field theory (EFT) for fluctuating hydrodynamics. Treating sound waves as the linear response to thermal radiation, we establish the fluctuation-dissipation relation within a cosmological framework. While dissipation is elucidated in established linear cosmological perturbation theory, the standard Boltzmann theory overlooks the associated noise, possibly contributing to inconsistencies in Lambda Cold Dark Matter ($\Lambda$CDM) cosmology. This paper employs EFT for fluctuating hydrodynamics in cosmological perturbation theory, deriving sound wave noise. Notably, the long-time limit of the noise spectrum is independent of viscosity details, resembling a Brownian motion bounded in a harmonic potential. The net energy transfer between the sound wave system and the radiation environment reaches a balance within Hubble time, suggesting the thermal equilibrium of the sound waves themselves. The induced density power spectrum is characterized as white noise dependent on the inverse of the entropy density, which is negligibly small on the CMB scale. The energy density of the entire sound wave system scales as $a^{-4}$, akin to radiation. While the numerical factor is not determined in the present calculation, the back reaction of the sound wave system to the background radiation may not be negligible, serving as a potential source for various fitting issues in $\Lambda$CDM cosmology.