Papers for Friday, May 24 2024

Variations of the magnetic field within solar coronal mass ejections (CMEs) in the heliosphere depend on the CME`s magnetic structure as it leaves the solar corona and its subsequent evolution through interplanetary space. To account for this evolution, we developed a new numerical model of the inner heliosphere that simulates the propagation of a CME through a realistic background solar wind and allows various CME magnetic topologies. To this end, we incorporate the Gibson-Low CME model within our global MHD model of the inner heliosphere, GAMERA-Helio. We apply the model to study the propagation of the geoeffective CME that erupted on 3 April, 2010 with the aim to reproduce the temporal variations of the magnetic field vector during the CME passage by Earth. Parameters of the Gibson-Low CME are informed by STEREO white-light observations near the Sun. The magnetic topology for this CME - the tethered flux rope - is informed by in-situ magnetic field observations near Earth. We performed two simulations testing different CME propagation directions. For an in-ecliptic direction, the simulation shows a rotation of all three magnetic field components within the CME, as seen at Earth, similar to that observed. With a southward propagation direction, suggested by coronal imaging observations, the modeled By and Bz components are consistent with the ACE data, but the Bx component lacks the observed change from negative to positive. In both cases, the model favors the East-West orientation of the CME flux rope, consistent with the orientation previously inferred from the STEREO/HI heliospheric images.

This study is devoted to dilaton generation in propagation of magnetic dipole waves of pulsar in a galactic magnetic field. It is shown that in this process it's necessary to take into account the presence of a reflective index of matter distributed in the galaxy. According to obtained estimates influence of the gravitational field of pulsars and magnetars on radiation of magnetic dipole waves is negligible and tends to a small change of their amplitude, near about 1 per cent. That's why in this study one can neglect the impact of gravitational field on process of dilaton generation. Exact solution of the dilaton equation is found and angular distribution of dilaton radiation is obtained in propagation of magnetic dipole waves of pulsar in galactic magnetic field.

Recent studies using four-point correlations suggest a parity violation in the galaxy distribution, though the significance of these detections is sensitive to the choice of simulation used to model the noise properties of the galaxy distribution. In a recent paper, we introduce an unsupervised learning approach which offers an alternative method that avoids the dependence on mock catalogs, by learning parity violation directly from observational data. However, the Convolutional Neural Network (CNN) model utilized by our previous unsupervised approach struggles to extend to more realistic scenarios where data is limited. We propose a novel method, the Neural Field Scattering Transform (NFST), which enhances the Wavelet Scattering Transform (WST) technique by adding trainable filters, parameterized as a neural field. We first tune the NFST model to detect parity violation in a simplified dataset, then compare its performance against WST and CNN benchmarks across varied training set sizes. We find the NFST can detect parity violation with $4\times$ less data than the CNN and $32\times$ less than the WST. Furthermore, in cases with limited data the NFST can detect parity violation with up to $6\sigma$ confidence, where the WST and CNN fail to make any detection. We identify that the added flexibility of the NFST, and particularly the ability to learn asymmetric filters, as well as the specific symmetries built into the NFST architecture, contribute to its improved performance over the benchmark models. We further demonstrate that the NFST is readily interpretable, which is valuable for physical applications such as the detection of parity violation.

Starlink satellites can become extremely bright when sunlight reflects specularly to an observer on the ground. The observed brightness of such flares is consistent with a bidirectional reflectance function of the Starlink satellite chassis. These findings are applied to the case of an extreme flare that was reported as an Unidentified Aerial Phenomena by the pilots of two commercial aircraft.

We present an analysis of X-ray observations of the Ultraluminous X-ray source (ULX) in IZw18 based on archival data taken with Chandra, XMM-Newton, and Suzaku. This ULX is considered to be an intermediate-mass black hole candidate simply because it is in the lowest metallicity environment among ULXs, where formation of heavy black holes is facilitated. However, actual study of the ULX based on observations spanning for a long period has been too limited to determine its nature. In this study, we investigate the spectral evolution of the ULX up to 2014, combining the previously-unpublished Suzaku data with those from the other two satellites. We derive a positive correlation of $L\propto T_{\rm in}^{2.1\pm0.4}$ between the bolometric luminosity $L$ and inner-disk temperature $T_{\rm in}$ on the basis of the multi-color disk-blackbody model, where we exclude the Chandra data, which has the lowest luminosity and systematic residuals in the fitting. The nominal relation $L\propto T_{\rm in}^{4}$ for the standard disk is rejected at a significance level of 1.5 %. These results suggest that the ULX was in the slim-disk state during these observations except at the time of the Chandra observation, in which the ULX was likely to be in a different state. The apparent inner-disk radius appears negatively correlated with the inner-disk temperature. Moreover, we find a radial dependence of the disk temperature of $T (r)\propto r^{-p}$ with $p<0.75$, which also supports the hypothesis that the ULX has a slim disk. In consequence, the IZw18 ULX is most likely to be powered by a stellar-mass compact object in supercritical accretion.

JWST is unveiling a surprising lack of evolution in the number densities of ultraviolet-selected (UV) galaxies at redshift $z\gtrsim 10$. At the same time, observations and simulations are providing evidence for highly bursty star formation in high-$z$ galaxies, resulting in significant scatter in their UV luminosities. Galaxies in low-mass dark matter halos are expected to experience most stochasticity due to their shallow potential wells. Here, we explore the impact of a mass-dependent stochasticity using a simple analytical model. We assume that scatter in the $M_\mathrm{UV}-M_h$ relation increases towards lower halo masses, following the decrease in halo escape velocity, $\sigma_\mathrm{UV} \sim M_h^{-1/3}$, independent of redshift. Since low-mass halos are more dominant in the early universe, this model naturally predicts an increase in UV luminosity functions (LFs) at high redshifts compared to models without scatter. We make predictions for additional observables which would be affected by stochasticity and could be used to constrain its amplitude, finding: (i) galaxies are less clustered compared to the no-scatter scenario, with the difference increasing at higher-$z$; (ii) assuming star-bursting galaxies dominate the ionizing photon budget implies reionization starts earlier and is more gradual compared to the no-scatter case, (iii) at fixed UV magnitude galaxies should exhibit wide ranges of UV slopes, nebular emission line strengths and Balmer breaks. Comparing to observations, the mass-dependent stochasticity model successfully reproduces the observed LFs up to $z\sim12$. However, the model cannot match the observed $z\sim14$ LFs, implying additional physical processes enhance star formation efficiency in the earliest galaxies.

Elucidating the processes that shape the circumgalactic medium (CGM) is crucial for understanding galaxy evolution. Absorption and emission diagnostics can be interpreted using photoionization calculations to obtain information about the phase and ionization structure of the CGM. For simplicity, typically only the metagalactic background is considered in photoionization calculations, and local sources are ignored. To test this simplification, we perform Monte Carlo radiation transfer on 12 cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project with halo masses $10^{10.5}-10^{13} \mathrm{M}_{\odot}$ in the redshift range $z = 0-3.5$ to determine the spatial extent over which local sources appreciably contribute to the ionizing radiation field in the CGM. We find that on average, the contribution of stars within the galaxy is small beyond one-tenth of the virial radius, $R_{\mathrm{vir}}$, for $z < 1$. For $1<z<2$ and $M_{\mathrm{vir}} \sim 10^{11.5}$, the radius at which the contribution to the ionizing radiation field from stars within the galaxy and that from the UV background are equal is roughly 0.2 $R_{\mathrm{vir}}$. For $M_{\mathrm{vir}} > 10^{12} \mathrm{M}_{\odot}$ at $z \sim 1.5-2.5$ and for all $M_{\mathrm{vir}}$ considered at $z>3$ , this transition radius can sometimes exceed 0.5 $R_{\mathrm{vir}}$. We also compute the escape fraction at $R_{\mathrm{vir}}$, finding typical values of less than $0.1$, except in higher-mass halos ($M_{\mathrm{halo}} \gtrsim 10^{12} \mathrm{M}_{\odot}$), which have consistently high values of $\sim 0.5-0.6$. Our results indicate that at low redshift, it is reasonable to ignore the ionizing radiation from host-galaxy stars outside of 0.2 $R_{\mathrm{vir}}$, while at Cosmic Noon, local stellar ionizing radiation likely extends further into the CGM and thus should be included in photoionization calculations.

Circumgalactic Lyman-alpha (Ly$\alpha$) nebulae are gaseous halos around galaxies exhibiting luminous extended Ly$\alpha$ emission. This work investigates Ly$\alpha$ nebulae from deep imaging of $\sim12~\mathrm{deg}^2$ sky, targeted by the MAMMOTH-Subaru survey. Utilizing the wide-field capability of Hyper Suprime-Cam (HSC), we present one of the largest blind Ly$\alpha$ nebula selections, including QSO nebulae, Ly$\alpha$ blobs, and radio galaxy nebulae down to typical $2\sigma$ Ly$\alpha$ surface brightness of $(5-10)\times10^{-18}\mathrm{~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$. The sample contains 117 nebulae with Ly$\alpha$ sizes of 40 - 400 kpc, and the most gigantic one spans about 365 kpc, referred to as the Ivory Nebula. Combining with multiwavelength data, we investigate diverse nebula populations and associated galaxies. We find a small fraction of Ly$\alpha$ nebulae have QSOs ($\sim7\%$), luminous infrared galaxies ($\sim1\%$), and radio galaxies ($\sim 2\%$). Remarkably, among the 28 enormous Ly$\alpha$ nebulae (ELANe) exceeding 100 kpc, about $80\%$ are associated with UV-faint galaxies ($M_\mathrm{UV} > -22$), categorized as Type II ELANe. We underscore that Type II ELANe constitute the majority but remain largely hidden in current galaxy and QSO surveys. Dusty starburst and obscured AGN activity are proposed to explain the nature of Type II ELANe. The SED of stacking all Ly$\alpha$ nebulae also reveals signs of massive dusty star-forming galaxies with obscured AGNs. We propose a model to explain the dusty nature where the diverse populations of Ly$\alpha$ nebula capture massive galaxies at different evolutionary stages undergoing violent assembling. Ly$\alpha$ nebulae provide critical insights into the formation and evolution of today's massive cluster galaxies at cosmic noon.

The origin and evolution of CO gas in debris disks has been debated since its initial detection. The gas could have a primordial origin, as a remnant of the protoplanetary disk or a secondary exocometary origin. This paper investigates the origin of gas in two debris disks, HD110058 and HD131488, using HST observations of CI and CO, which play critical roles in the gas evolution. We fitted several electronic transitions of CI and CO rovibronic bands to derive column densities and temperatures for each system, revealing high CO column densities ($\sim$3-4 orders of magnitude higher than $\beta$ Pictoris), and low CI/CO ratios in both. Using the exogas model, we simulated the radial evolution of the gas in the debris disk assuming a secondary gas origin. We explored a wide range of CO exocometary release rates and $\alpha$ viscosities, which are the key parameters of the model. Additionally, we incorporated photodissociation due to stellar UV to the exogas model and found that it is negligible for typical CO-rich disks and host stars, even at a few au due to the high radial optical depths in the EUV. We find that the current steady-state secondary release model cannot simultaneously reproduce the CO and CI HST-derived column densities, as it predicts larger CI/CO ratios than observed. Our direct UV measurement of low CI/CO ratios agrees with results derived from recent ALMA findings and may point to vertical layering of CI, additional CI removal, CO shielding processes, or different gas origin scenarios.

We report on the discovery of Gliese 12 b, the nearest transiting temperate, Earth-sized planet found to date. Gliese 12 is a bright ($V=12.6$ mag, $K=7.8$ mag) metal-poor M4V star only $12.162\pm0.005$ pc away from the Solar System with one of the lowest stellar activity levels known for an M-dwarf. A planet candidate was detected by TESS based on only 3 transits in sectors 42, 43, and 57, with an ambiguity in the orbital period due to observational gaps. We performed follow-up transit observations with CHEOPS and ground-based photometry with MINERVA-Australis, SPECULOOS, and Purple Mountain Observatory, as well as further TESS observations in sector 70. We statistically validate Gliese 12 b as a planet with an orbital period of $12.76144\pm0.00006$ days and a radius of $1.0\pm{0.1}$ R$_\oplus$, resulting in an equilibrium temperature of $\sim$315K. Gliese 12 b has excellent future prospects for precise mass measurement, which may inform how planetary internal structure is affected by the stellar compositional environment. Gliese 12 b also represents one of the best targets to study whether Earth-like planets orbiting cool stars can retain their atmospheres, a crucial step to advance our understanding of habitability on Earth and across the Galaxy.

We train a novel deep learning architecture to perform likelihood-free inference on the value of the cosmological parameters from halo catalogs of the Quijote N-body simulations. Our model takes as input a halo catalog where each halo is characterized by its position, mass, and velocity moduli. By construction, our model is E(3) invariant and is designed to extract information hierarchically. Unlike graph neural networks, it does not require the transformation of the input halo (or galaxy) catalog into a graph. Given its simplicity, our model can process point clouds with large numbers of points. We discuss the advantages of this class of methods but also point out its limitations and potential ways to improve them for cosmological data.

Galactic outflows influence the evolution of galaxies not only by expelling gas from their disks but also by injecting energy into the circumgalactic medium (CGM). This alters or even prevents the inflow of fresh gas onto the disk and thus reduces the star formation rate. Supernovae (SNe) are the engines of galactic winds as they release thermal and kinetic energy into the interstellar medium (ISM). Cosmic rays (CRs) are accelerated at the shocks of SN remnants and only constitute a small fraction of the overall SN energy budget. However, their long live-times allow them to act far away from the original injection site and thereby to participate in the galactic wind launching process. Using high-resolution simulations of an isolated Milky Way-type galaxy with the moving-mesh code Arepo and the new multi-phase ISM model Crisp (Cosmic Rays and InterStellar Physics), we investigate how SNe and CRs launch galactic outflows and how the inclusion of CR-mediated feedback boosts the energy and mass entrained in the galactic wind. We find that the majority of thermal SN energy and momentum is used for stirring turbulence either directly or indirectly by causing fountain flows, thereby self-regulating the ISM and not for efficiently driving outflows to large heights. A simulation without CRs only launches a weak galactic outflow at uniformly high temperatures and low densities by means of the thermal pressure gradient. By contrast, most of the CR energy accelerated at SN remnants ($\sim80\%$) escapes the ISM and moves into the CGM. In the inner CGM, CRs dominate the overall pressure and are able to accelerate a large mass fraction in a galactic wind. This wind is turbulent and multi-phase with cold cloudlets embedded in dilute gas at intermediate temperatures ($\sim10^5$ K) and the CGM shows enhanced OVI and CIV absorption in comparison to a simulation without CRs.

We present Keck Cosmic Web Imager IFU observations around extended Ly$\alpha$ halos of 27 typical star-forming galaxies with redshifts $2.0 < z < 3.2$ drawn from the MOSFIRE Deep Evolution Field survey. We examine the average Ly$\alpha$ surface-brightness profiles in bins of star-formation rate (SFR), stellar mass ($M_*$), age, stellar continuum reddening, SFR surface density ($\rm \Sigma_{SFR}$), and $\rm \Sigma_{SFR}$ normalized by stellar mass ($\rm \Sigma_{sSFR}$). The scale lengths of the halos correlate with stellar mass, age, and stellar continuum reddening; and anti-correlate with star-formation rate, $\rm \Sigma_{SFR}$, and $\rm \Sigma_{sSFR}$. These results are consistent with a scenario in which the down-the-barrel fraction of Ly$\alpha$ emission is modulated by the low-column-density channels in the ISM, and that the neutral gas covering fraction is related to the physical properties of the galaxies. Specifically, we find that this covering fraction increases with stellar mass, age, and $E(B-V)$; and decreases with SFR, $\rm \Sigma_{SFR}$ and $\rm \Sigma_{sSFR}$. We also find that the resonantly scattered Ly$\alpha$ emission suffers greater attenuation than the (non-resonant) stellar continuum emission, and that the difference in attenuation increases with stellar mass, age, and stellar continuum reddening, and decreases with $\rm \Sigma_{sSFR}$. These results imply that more reddened galaxies have more dust in their CGM.

The local gravitational instability of rotating discs is believed to be an important mechanism in different astrophysical processes, including the formation of gas and stellar clumps in galaxies. We aim to study in three dimensions the local gravitational instability of two-component thick discs. We take as starting point a recently proposed analytic three-dimensional (3D) instability criterion for discs with non-negligible thickness which takes the form $Q_{\rm 3D}<1$, where $Q_{\rm 3D}$ is a 3D version of the classical 2D Toomre $Q$ parameter for razor-thin discs. Here we extend the 3D stability analysis to two-component discs, considering first the influence on $Q_{\rm 3D}$ of a second unresponsive component, and then the case in which both components are responsive. We present the application to two-component discs with isothermal vertical distributions, which can represent, for instance, galactic discs with both stellar and gaseous components. Finally, we relax the assumption of vertical isothermal distribution, by studying one-component self-gravitating discs with polytropic vertical distributions for a range of values of the polytropic index corresponding to convectively stable configurations. We find that $Q_{\rm 3D}<1$, where $Q_{\rm 3D}$ can be computed from observationally inferred quantities, is a robust indicator of local gravitational instability, depending only weakly on the presence of a second component and on the vertical gradient of temperature or velocity dispersion. We derive a sufficient condition for local gravitational instability in the midplane of two-component discs, which can be employed when both components have $Q_{\rm 3D}>1$.

The Pristine-\textit{Gaia} synthetic catalogue provides reliable photometric metallicities for $\sim$30 million FGK stars using the Pristine survey model and Gaia XP spectra. We perform the first low-to-medium-resolution spectroscopic follow-up of bright (G<15) and distant (up to 35 kpc) very and extremely metal-poor (V/EMP, [Fe/H]<-2.5) red giant branch stars from this. We use Isaac Newton Telescope/Intermediate Dispersion Spectrograph (INT/IDS) observations centred around the calcium triplet region ideal for V/EMP stars. We find that 76\% of our stars indeed have [Fe/H]<-2.5 with these inferred spectroscopic metallicities and only 3\% are outliers with [Fe/H] > -2.0. We report a success rate of 77\% and 38\% in finding stars with [Fe/H]<-2.5 and -3.0 respectively. This will allow for 10,000-20,000 homogeneously analysed EMP stars using the WEAVE survey follow-up of Pristine EMP candidates. We associate 20\%, 46\%, and 34\% of the stars to be confined to the disc plane, or to have inner and outer halo orbits, respectively. We also associate these V/EMP stars to known accretion events such as Gaia-Enceladus-Sausage (GES), LMS-1/Wukong, Thamnos, Helmi streams, Sagittarius, Sequoia, etc. For the stars that orbit close to the disc plane, we find that the prograde region with low vertical action is overdense with a significance of 4$\sigma$ as compared to its retrograde counterpart. We also find three new (brightest) members of the most metal-poor stellar stream, C-19, one of which is 50$^\circ$ away from the main body of the stream. Our measured mean metallicity, velocity dispersion, and stream width are consistent with the literature, but our results favour a higher distance ($\sim$21.5 kpc) for the stream. We publish a catalogue (and 1D spectra) of 215 V/EMP stars from this spectroscopic follow-up and showcase the power of chemokinematic analysis of V/EMP end.

We used the Keck Planet Imager and Characterizer (KPIC) to obtain high-resolution (R$\sim$35,000) K-band spectra of kappa Andromedae b, a planetary-mass companion orbiting the B9V star, kappa Andromedae A. We characterized its spin, radial velocity, and bulk atmospheric parameters through use of a forward modeling framework to jointly fit planetary spectra and residual starlight speckles, obtaining likelihood-based posterior probabilities. We also detected H$_{2}$O and CO in its atmosphere via cross correlation. We measured a $v\sin(i)$ value for kappa And b of $38.42\pm{0.05}$ km/s, allowing us to extend our understanding of the population of close in bound companions at higher rotation rates. This rotation rate is one of the highest spins relative to breakup velocity measured to date, at close to $50\%$ of breakup velocity. We identify a radial velocity $-17.35_{-0.09}^{+0.05}$ km/s, which we use with existing astrometry and RV measurements to update the orbital fit. We also measure an effective temperature of $1700\pm{100}$ K and a $\log(g)$ of $4.7\pm{0.5}$ cgs dex.

Using Keck Planet Imager and Characterizer (KPIC) high-resolution ($R$~35000) spectroscopy from 2.29-2.49 $\mu$m, we present uniform atmospheric retrievals for eight young substellar companions with masses of ~10-30 $M_\textrm{Jup}$, orbital separations spanning ~50-360 au, and $T_\textrm{eff}$ between ~1500-2600 K. We find that all companions have solar C/O ratios, and metallicities, to within the 1-2$\sigma$ level, with the measurements clustered around solar composition. Stars in the same stellar associations as our systems have near-solar abundances, so these results indicate that this population of companions is consistent with formation via direct gravitational collapse. Alternatively, core accretion outside the CO snowline would be compatible with our measurements, though the high mass ratios of most systems would require rapid core assembly and gas accretion in massive disks. On a population level, our findings can be contrasted with abundance measurements for directly imaged planets with m<10 $M_\textrm{Jup}$, which show tentative atmospheric metal enrichment. In addition, the atmospheric compositions of our sample of companions are distinct from those of hot Jupiters, which most likely form via core accretion. For two companions with $T_\textrm{eff}$~1700-2000 K (kap And b and GSC 6214-210 b), our best-fit models prefer a non-gray cloud model with >3$\sigma$ significance. The cloudy models yield 2-3$\sigma$ lower $T_\textrm{eff}$ for these companions, though the C/O and [C/H] still agree between cloudy and clear models at the $1\sigma$ level. Finally, we constrain 12CO/13CO for three companions with the highest S/N data (GQ Lup b, HIP 79098 b, and DH Tau b), and report $v$sin($i$) and radial velocities for all companions.

The application of clustering algorithms to the Gaia astrometric catalog has revolutionized our census of stellar populations in the Milky Way, including the discovery of many new, dispersed structures. We focus on one such structure, Theia 456 (COIN-Gaia-13), a loosely bound collection of ~320 stars spanning ~120 pc that has previously been shown to exhibit kinematic, chemical, and gyrochronal coherency, indicating a common origin. We obtain follow-up radial velocities and supplement these with Gaia astrometry to perform an in-depth dynamical analysis of Theia 456. By integrating stellar orbits through a Milky Way potential, we find the currently dispersed structure coalesced into a small cluster in the past. Via Bayesian modeling, we derive a kinematic age of 245 +/- 3 Myr (statistical), a half-mass radius of 9 +/- 2 pc, and an initial velocity dispersion of 0.14 +/- 0.02 km/s. Our results are entirely independent of model isochrones, details of stellar evolution, and internal cluster dynamics, and the statistical precision in our age derivation rivals that of the most precise age-dating techniques known today, though our imperfect knowledge of the Milky Way potential and simple spherical model for Theia 456 at birth add additional uncertainties. Using posterior predictive checking, we confirm these results are robust under reasonable variations to the Milky Way potential. Such low density structures that are disrupted by the Galactic tides before virializing may be ubiquitous, signifying that Theia 456 is a valuable benchmark for studying the dynamical history of stellar populations in the Milky Way.

We present the results from a study of ~9,600 Broad-Line selected AGN with host galaxies detected from the Sloan Digital Sky Survey Data Release 17 (SDSS DR17). We compute ensemble variability statistics based on the comparison of the original SDSS photometric data with spectrophotometric measurements obtained days to decades later in the Sloan g-, r-, and i-bands. Galaxy and AGN templates have been fitted to the SDSS spectra to isolate the AGN component from the host galaxy. The sources have absolute magnitudes in the range -24<M_i<-18 and lie at redshifts less than z~0.9. A variability analysis reveals that the anti-correlation between luminosity and variability amplitude continues down to log(L_bol[erg/s]) = 43.5, demonstrating that the relationship extends by 4 orders of magnitude in AGN luminosity. To further explore the connection between AGN luminosity and variability, we determine the black hole mass and the accretion rate through measurement of the H-Beta line width and the monochromatic luminosity at rest-frame 5,100 angstrom. Our results suggest that the accretion rate is the dominant parameter impacting the amplitude of variability and that the anti-correlation between accretion rate and amplitude extends to rates as low as 1% Eddington. Moreover, we also identify an anti-correlation between variability amplitude and black hole mass, with the correlation appearing strongest among the AGN with low accretion rates.

The radio transient ASKAP J173608.2-321735, at the position (l,b)= (356.0872,-0.0390), was serendipitously observed by The HI/OH/Recombination Line Survey of the Galactic Center (THOR-GC) at three epochs in March 2020, April 2020 and February 2021. The source was detected only on 2020 April 11 with flux density 20.6 +/- 1.1 mJy at 1.23 GHz and in-band spectral index alpha = -3.1 +/- 0.2. The commensal VLA Low-band Ionsophere and Transient Experiment (VLITE) simultaneously detected the source at 339 MHz with a flux density 122.6 +/- 20.4 mJy, indicating a spectral break below 1 GHz. The rotation measure in April 2020 was 63.9 +/- 0.3rad/m2, which almost triples the range of the variable rotation measure observed by Wang et al. (2021) to ~130 rad/m2. The polarization angle, corrected for Faraday rotation, was 97 +/- 6 degrees. The 1.23 GHz linear polarization was 76.7% +/- 3.9% with wavelength-dependent depolarization indicating Faraday depth dispersion sigma_phi = 4.8^{+0.5}_{-0.7} rad/m2. We find an upper limit to circular polarization |V|/I < 10.1%. Interpretation of the data in terms of diffractive scattering of radio waves by a plasma near the source indicates electron density and line-of-sight magnetic field strength within a factor 3 of n_e ~2 cm^{-3} and B_par ~2 x 10^5 microgauss. Combined with causality limits to the size of the source, these parameters are consistent with the low-frequency spectral break resulting from synchrotron self-absorption, not free-free absorption. A possible interpretation of the source is a highly supersonic neutron star interacting with a changing environment.

Luminous, compact emission-line galaxies (LCGs) are the most abundant class of confirmed Lyman continuum (LyC) emitters. An optical integral field study of the nearby LCG NGC 2366 reveals an outflow originating at the star cluster 'knot B' thought to clear a channel via mechanical feedback that enables LyC escape. We observed NGC 2366 with the Chandra X-ray Observatory and detect X-ray emission from a point source coincident with the apex of the outflow at knot B. The point-like nature and variability of the X-ray emission suggests accretion onto a compact object. The accretion could produce sufficient kinetic energy to power the outflow.

The Simons Observatory (SO), a next-generation ground-based CMB experiment in its final stages of construction, will target primordial $B$-modes with unprecedented sensitivity to set tight bounds on the amplitude of inflationary gravitational waves. Aiming to infer the tensor-to-scalar ratio $r$ with precision $\sigma(r=0) \leq 0.003$, SO will rely on powerful component-separation algorithms to distinguish the faint primordial signal from stronger sources of large-scale $B$-modes such as Galactic foregrounds and weak gravitational lensing. We present an analysis pipeline that performs delensing and foreground cleaning simultaneously by including multifrequency CMB data and a lensing $B$-mode template in a power-spectrum-based likelihood. Here, we demonstrate this algorithm on masked SO-like simulations containing inhomogeneous noise and non-Gaussian foregrounds. The lensing convergence is reconstructed from high-resolution simulations of the CMB and external mass tracers. Using optimized pixel weights for power spectrum estimation, the target precision for SO's nominal design is achieved and delensing reduces $\sigma(r)$ by 27-37%, depending on foreground complexity.

We present an effective field theory (EFT) approach to extract fundamental cosmological parameters from the Lyman-alpha forest flux fluctuations as an alternative to the standard simulation-based techniques. As a first application, we re-analyze the publicly available one-dimensional Lyman-alpha flux power spectrum (FPS) data from the Sloan Digital Sky Survey (SDSS). Our analysis relies on informative priors on EFT parameters which we extract from a combination of public hydrodynamic simulation and emulator data. Assuming the concordance cosmological model, our one-parameter analysis yields a $2\%$ measurement of the late time mass fluctuation amplitude $\sigma_8 = 0.841\pm 0.017$, or equivalently,the structure growth parameter $S_8 = 0.852\pm 0.017$, consistent with the standard cosmology. Combining our EFT likelihood with Planck + baryon acoustic oscillation data, we find a new constraint on the total neutrino mass, $\sum m_\nu<$ 0.08 eV (at 95\% CL). Our study defines priorities for the development of EFT methods and sets the benchmark for cosmological analyses of the Lyman-alpha forest data from the Dark Energy Spectroscopic Instrument.

Following the milestone of detecting gravitational waves (GWs) from merging compact binaries, the next significant watershed moment in GW astronomy lies in detecting GWs from core-collapse supernovae (CCSNe). In this Letter, I describe the possibility of detecting the GW linear memory -- a phenomenon resulting from a combination of aspherical matter ejection and anisotropic neutrino emission during stellar collapse using GW detectors on the Moon. This would grant unprecedented access to the sub-Hz/Hz GW frequency range, which is inaccessible to current and future terrestrial GW detectors. I demonstrate that three-dimensional CCSNe model matter and neutrino GW waveforms may be detectable by both seismometer and interferometeric lunar GW detectors, with the latter design proposal extending 10 kpc to megaparsec detection distances.

We searched the Milky Way Plane along a 6-deg swath for pulses of monochromatic light as faint as 15th mag (V band) using a wide-field telescope equipped with a prism. Pulses with duration less than 1 second that occur more often than once every 10 minutes would be detected, and pulses arriving less frequently would be detected with proportionally lower probability. No unexplainable monochromatic emission, pulsed or continuous, was detected. The detection threshold corresponds to a 70 GW laser having a diffraction-limited 10-meter aperture located 1 kiloparsec away (depending on wavelength). Previous searches for laser emission from more than 5000 stars found none. Previous all-sky surveys at optical and radio wavelengths have revealed thousands of unexpected objects in the universe that exhibited extraordinary spectral emission, but none were technological. Hypotheses of our Milky Way Galaxy teeming with advanced life must be demoted.

Many novel methods have been proposed to mitigate stellar activity for exoplanet detection as the presence of stellar activity in radial velocity (RV) measurements is the current major limitation. Unlike traditional methods that model stellar activity in the RV domain, more methods are moving in the direction of disentangling stellar activity at the spectral level. The goal of this paper is to present a novel convolutional neural network-based algorithm that efficiently models stellar activity signals at the spectral level, enhancing the detection of Earth-like planets. We trained a convolutional neural network to build the correlation between the change in the spectral line profile and the corresponding RV, full width at half maximum (FWHM) and bisector span (BIS) values derived from the classical cross-correlation function. This algorithm has been tested on three intensively observed stars: Alpha Centauri B (HD128621), Tau ceti (HD10700), and the Sun. By injecting simulated planetary signals at the spectral level, we demonstrate that our machine learning algorithm can achieve, for HD128621 and HD10700, a detection threshold of 0.5 m/s in semi-amplitude for planets with periods ranging from 10 to 300 days. This threshold would correspond to the detection of a $\sim$4$\mathrm{M}_{\oplus}$ in the habitable zone of those stars. On the HARPS-N solar dataset, our algorithm is even more efficient at mitigating stellar activity signals and can reach a threshold of 0.2 m/s, which would correspond to a 2.2$\mathrm{M}_{\oplus}$ planet on the orbit of the Earth. To the best of our knowledge, it is the first time that such low detection thresholds are reported for the Sun, but also for other stars, and therefore this highlights the efficiency of our convolutional neural network-based algorithm at mitigating stellar activity in RV measurements.

An accurate calculation of their abundance is crucial for numerous aspects of cosmology related to primordial black holes (PBHs). For example, placing constraints on the primordial power spectrum from constraints on the abundance of PBHs (or vice-versa), calculating the mass function observable today, or predicting the merger rate of (primordial) black holes observable by gravitational wave observatories such as LIGO, Virgo and KAGRA. In this chapter, we will discuss the different methods used for the calculation of the abundance of PBHs forming from large-amplitude cosmological perturbations, assuming only a minimal understanding of modern cosmology. Different parameters to describe cosmological perturbations will be considered (including different choices for the window function), and it will be argued that the compaction is typically the most appropriate choice. Different methodologies for calculating the abundance and mass function are explained, including \emph{Press-Schechter}-type and peaks theory approaches.

We applied a Density-Based Clustering algorithm on samples of galaxies and galaxy systems belonging to 53 rich superclusters from the \textit{Main SuperCluster Catalogue} (MSCC) to identify the presence of ``central regions'', or \emph{cores}, in these large-scale structures. \emph{Cores} are defined here as large gravitationally bound galaxy structures, comprised of two or more clusters and groups, with sufficient matter density to survive cosmic expansion and virialize in the future. We identified a total of 105 galaxy structures classified as \emph{cores}, which exhibit a high density contrast of mass and galaxies. The Density-based \textit{Core} Catalogue (DCC), presented here, includes \emph{cores} that were previously reported in well-known superclusters of the Local Universe, and also several newly identified ones. We found that 83\% of the rich superclusters in our sample have at least one \emph{core}. While more than three \emph{cores} with different dynamical state are possible, the presence of a single \emph{core} in the superclusters is more common. Our work confirms the existence of nucleation regions in the internal structure of most rich superclusters and points to the fact that these \emph{cores} are the densest and most massive features that can be identified in the cosmic web with high probability for future virialization.

The hydrodynamic escape driven by external or internal energy sources sculpts the population of low mass close-in planets. However, distinguishing between the driving mechanisms responsible for the hydrodynamic escape of hydrogen-rich atmospheres is a complex task due to the involvement of many physical factors. My simulations show that the hydrodynamic escape can be driven solely by thermal energy deposited in the lower layers of the atmosphere due to the heat flux originating from the planetary core or bolometric heating from the star even in the absence of other energy sources, as long as the planet's Jeans parameter is below 3. Otherwise, stellar extreme ultraviolet irradiation or tidal forces are necessary in driving the escape, which means that the Jeans parameter is incapable of distinguishing the driving mechanisms, as it is only related to the properties of planet. Here, an upgraded Jeans parameter that takes into account tidal forces is introduced, which allows us to accurately categorize the driving mechanisms. The results show that when the upgraded Jeans parameter falls below 3 or exceeds 6, the atmospheric escape is primarily driven by tidal forces or extreme ultraviolet radiation from the host star, respectively. In the range of 3 to 6, both factors can trigger the escape of the atmosphere. I find that planets with high gravitational potential and low stellar irradiation are more likely to undergo subsonic escape, although transonic escape is prevalent among most planets. Moreover, the ionization status is significantly dependent on the gravitational potential. The upgraded Jeans parameter, which is closely related to the underlying physics, provides a concise method to categorize the driving mechanisms of hydrodynamic escape. The results can be applied to planetary evolution calculations.

Based on the newly acquired dense gas observations from the JCMT MALATANG survey and X-ray data from Chandra, we explore the correlation between hot gas and HCN $J=4 \rightarrow 3$, HCO$^+\ J=4 \rightarrow 3$ emission for the first time at sub-kiloparsec scale of five nearby star-forming galaxies, namely M82, M83, IC 342, NGC 253, and NGC 6946. We find that both HCN $J=4 \rightarrow 3$ and HCO$^+\ J=4 \rightarrow 3$ line luminosity show a statistically significant correlation with the 0.5${-}$2 keV X-ray emission of the diffuse hot gas ($L_{\rm 0.5 - 2\,keV}^{\rm gas}$). The Bayesian regression analysis gives the best fit of ${\rm log}(L_{\rm 0.5-2\,keV}^{\rm gas} /{\rm erg\,s^{-1}})=2.39\,{\rm log}(L'_{\rm HCN(4-3)} /{\rm K\,km\,s^{-1}\,pc^{2}})+24.83$ and ${\rm log}(L_{\rm 0.5-2\,keV}^{\rm gas} /{\rm erg\,s^{-1}})=2.48\,{\rm log}(L'_{\rm HCO^{+}(4-3)} /{\rm K\,km\,s^{-1}\,pc^{2}})+23.84$, with dispersion of $\thicksim$0.69 dex and 0.54 dex, respectively. At the sub-kiloparsec scale, we find that the power-law index of the $L_{\rm 0.5 - 2\,keV}^{\rm gas}$ ${-}$ star formation rate (SFR) relation is ${\rm log}(L_{\rm 0.5-2\,keV}^{\rm gas} /{\rm erg\,s^{-1}})=1.80\,{\rm log} ({\rm SFR} /M_\odot\,{\rm yr}^{-1})+39.16$, deviated from previous linear relations at global scale. This implies that the global property of hot gas significantly differs from individual resolved regions, which is influenced by the local physical conditions close to the sites of star formation.

Filament eruptions often result in flares and coronal mass ejections (CMEs). Most studies attribute the filament eruptions to their instabilities or magnetic reconnection. In this study, we report a unique observation of a filament eruption whose initiation process has not been reported before. This large-scale filament, with a length of about 360 Mm crossing an active region, is forced to erupted by two small-scale erupting filaments pushing out from below. This process of multi-filament eruption results in an M6.4 flare in the active region NOAA 13229 on 25th February 2023. The whole process can be divided into three stages: the eruptions of two active-region filaments F1 and F2; the interactions between the erupting F1, F2, and the large-scale filament F3; and the eruption of F3. Though this multi-filament eruption occurs near the northwest limb of the solar disk, it produces a strong halo CME that causes a significant geomagnetic disturbance. Our observations present a new filament eruption mechanism, in which the initial kinetic energy of the eruption is obtained from and transported to by other erupting structures. This event provides us a unique insight into the dynamics of multi-filament eruptions and their corresponding effects on the interplanetary space.

Theories of modified gravity suggest that the propagation speed of gravitational wave (GW) $v_g$ may deviate from the speed of light $c$. A constraint can be placed on the difference between $c$ and $v_g$ with a simple method that uses the arrival time delay between GW and electromagnetic (EM) wave simultaneously emitted from a burst event. We simulated the joint observation of GW and short Gamma-Ray burst (sGRB) signals from Binary Neutron Star (BNS) merger events in different observation campaigns, involving advanced LIGO (aLIGO) in design sensitivity and Einstein Telescope (ET) joint-detected with \textit{Fermi}/GBM. As a result, the relative precision of constraint on $v_g$ can reach $\sim 10^{-17}$ (aLIGO) and $\sim 10^{-18}$ (ET), which are one and two orders of magnitude better than that from GW170817, respectively. We continue to obtain the bound of graviton mass $m_g \leq 7.1(3.2)\times 10^{-20}\,$eV with aLIGO (ET). Applying the Standard-Model Extension (SME) test framework, the constraint on $v_g$ allows us to study the Lorentz violation in the nondispersive, nonbirefringent limit of the gravitational sector. We obtain the constraints of the dimensionless isotropic coefficients $\bar{s}_{00}^{(4)}$ at mass dimension $d = 4$, which are $-1\times 10^{-15}< \bar{s}_{00}^{(4)}<9\times 10^{-17}$ for aLIGO and $-4\times 10^{-16}< \bar{s}_{00}^{(4)}<8\times 10^{-18}$ for ET.

Molecular deuteration is commonly seen in starless cores and is expected to occur on a timescale comparable to that of the core contraction. Thus, the deuteration serves as a chemical clock, allowing us to investigate dynamical theories of core formation. We aim to provide a 3D cloud description for the starless core L 1498 located in the nearby low-mass star-forming region Taurus, and explore the possible core formation mechanism of L 1498. We carried out non-local thermal equilibrium radiative transfer with multi-transition observations of the high-density tracer N$_2$H$^+$ to derive the density and temperature profiles of the L 1498 core. Combining with the spectral observations of the deuterated species, ortho-H$_2$D$^+$, N$_2$D$^+$, and DCO$^+$, we derived the abundance profiles for observed species and performed chemical modeling of the deuteration profiles across L 1498 to constrain the contraction timescale. We present the first ortho-H$_2$D$^+$ (1$_{10}$-1$_{11}$) detection toward L 1498. We find a peak molecular hydrogen density of $1.6_{-0.3}^{+3.0}\times10^{5}$~cm$^{-3}$, a temperature of 7.5$_{-0.5}^{+0.7}$~K, and a N$_2$H$^+$ deuteration of 0.27$_{-0.15}^{+0.12}$ in the center. We derive a lower limit of the core age for L 1498 of 0.16~Ma which is compatible with the typical free-fall time, indicating that L 1498 likely formed rapidly.

(ABRIDGED) CONTEXT: The Villafranca project is combining Gaia data with ground-based surveys to analyze Galactic stellar groups with OB stars. AIMS: We want to analyze Stock 18 within the Villafranca project, a very young stellar cluster with a symmetrical and compact H II region around it. METHODS: We analyze the core, massive-star population, extinction, distance, membership, internal dynamics, density profile, age, IMF, total mass, stellar variability, and Galactic location of Stock 18 with Gaia data and ground-based spectroscopy. RESULTS: Stock 18 is a very young (~1.0 Ma) cluster located at a distance of 2.91+-0.10 kpc dominated by the GLS 13 370 system, whose primary is an O9 V star. We propose that Stock 18 was in a very compact state (~0.1 pc) about 1.0 Ma ago and that most massive stars were ejected at that time without significantly affecting the less massive stars as a result of multi-body dynamical interactions. Given its age close to 1.0 Ma, the dynamical interactions took place very soon after massive star formation. Well defined expanding stellar clusters have been observed before but none as young as this one. The IMF is top heavy but if we discard the ejected ones it becomes nearly canonical. Therefore, this is another example in addition to the one we previously found (the Bermuda cluster) of (a) a very young cluster with an already evolved PDMF (b) that has significantly contributed to the future population of free-floating compact objects. If confirmed in more clusters, the number of such compact objects may be higher in the Milky Way than previously thought. Stock 18 has a variable extinction with an average value of R_5495 higher than the canonical one of 3.1. The cluster is above our Galactic mid-plane and has a distinct motion with respect to its surrounding old population, which is possibly an influence of the Perseus spiral arm.

We propose a scenario that can explain the early-time inflation and the late-time dark energy within a unified framework. A scalar potential combining power-law and exponential type in a context of extended Jordan-Brans-Dicke gravity is critically important for this realization. A realistic scenario can be achieved in a two-field model in which one directional motion in field space realizes the slow-roll inflation. The inflaton ends up with oscillatory period and turns its direction to another direction that is identified as the quintessence field, giving rise to the dark energy at late times. The inflaton oscillation is expected to realize efficient heating if parametric amplification works. Along the quintessence direction, the present universe is on the way to reach the asymptotic fixed point. We search for successful parameter region, taking potential function in the form of low-order field powers times decreasing exponential in two dimensional field space.

We present a comprehensive photometric study of RR Lyrae stars in the M3 globular cluster, utilising a vast dataset of 3140 optical ($UBVRI$) CCD images spanning 35 years from astronomical data archives. We have successfully identified previously known 238 RR Lyrae stars from the photometric data, comprising 178 RRab, 49 RRc, and 11 RRd stars. Multi-band periodogram was used to significantly improve the long-term periods of $65\%$ of RR Lyrae stars in our sample, thanks to the unprecedentedly long temporal coverage of the observations. The light curve templates were used to obtain accurate and precise mean magnitudes and amplitudes of all RR Lyrae variables. We combined optical ($UBVRI$) and near-infrared (NIR, $JHK_{s}$) photometry of RR Lyrae variables to investigate their location in the colour-magnitude diagrams as well as the pulsation properties such as period distributions, Bailey diagrams and amplitude ratios. The Period-Luminosity relations in $R$ and $I$ bands and Period-Wesenheit relations were derived after excluding outliers identified in CMDs. The Period-Wesenheit relations calibrated via the theoretically predicted relations were used to determine a distance modulus of $\mu = 15.04 \pm 0.04 \,{\rm (stats)} \pm 0.19 \,{\rm {(syst.)}}$ mag (using metal-independent $W_{BV}$ Wesenheit) and $\mu = 15.03 \pm 0.04 \,{\rm (stats)} \pm 0.17 \,{\rm {(syst.)}}$ mag (using metal-dependent $W_{VI}$ Wesenheit). Our distance measurements are in excellent agreement with published distances to M3 in the literature. We also employed an artificial neural network based comparison of theoretical and observed light curves to determine physical parameters (mass, luminosity, and effective temperature) for $79$ non-Blazhko RRab stars that agree with limited literature measurements.

The advent of high-contrast imaging instruments combined with medium-resolution spectrographs allows spectral and temporal dimensions to be combined with spatial dimensions to detect and potentially characterize exoplanets with higher sensitivity. We develop a new method to effectively leverage the spectral and spatial dimensions in integral-field spectroscopy (IFS) datasets using a supervised deep-learning algorithm to improve the detection sensitivity to high-contrast exoplanets. We begin by applying a data transform whereby the IFS datasets are replaced by cross-correlation coefficient tensors obtained by cross-correlating our data with young gas giant spectral template spectra. This transformed data is then used to train machine learning (ML) algorithms. We train a 2D CNN and 3D LSTM with our data. We compare the ML models with a non-ML algorithm, based on the STIM map of arXiv:1810.06895. We test our algorithms on simulated young gas giants in a dataset that contains no known exoplanet, and explore the sensitivity of algorithms to detect these exoplanets at contrasts ranging from 1e-3 to 1e-4 at different radial separations. We quantify the sensitivity using modified receiver operating characteristic curves (mROC). We discover that the ML algorithms produce fewer false positives and have a higher true positive rate than the STIM-based algorithm, and the true positive rate of ML algorithms is less impacted by changing radial separation. We discover that the velocity dimension is an important differentiating factor. Through this paper, we demonstrate that ML techniques have the potential to improve the detection limits and reduce false positives for directly imaged planets in IFS datasets, after transforming the spectral dimension into a radial velocity dimension through a cross-correlation operation.

The new generation of observatories and instruments (VLT/ERIS, JWST, ELT) motivate the development of robust methods to detect and characterise faint and close-in exoplanets. Molecular mapping and cross-correlation for spectroscopy use molecular templates to isolate a planet's spectrum from its host star. However, reliance on signal-to-noise ratio (S/N) metrics can lead to missed discoveries, due to strong assumptions of Gaussian independent and identically distributed noise. We introduce machine learning for cross-correlation spectroscopy (MLCCS); the method aims to leverage weak assumptions on exoplanet characterisation, such as the presence of specific molecules in atmospheres, to improve detection sensitivity for exoplanets. MLCCS methods, including a perceptron and unidimensional convolutional neural networks, operate in the cross-correlated spectral dimension, in which patterns from molecules can be identified. We test on mock datasets of synthetic planets inserted into real noise from SINFONI at K-band. The results from MLCCS show outstanding improvements. The outcome on a grid of faint synthetic gas giants shows that for a false discovery rate up to 5%, a perceptron can detect about 26 times the amount of planets compared to an S/N metric. This factor increases up to 77 times with convolutional neural networks, with a statistical sensitivity shift from 0.7% to 55.5%. In addition, MLCCS methods show a drastic improvement in detection confidence and conspicuity on imaging spectroscopy. Once trained, MLCCS methods offer sensitive and rapid detection of exoplanets and their molecular species in the spectral dimension. They handle systematic noise and challenging seeing conditions, can adapt to many spectroscopic instruments and modes, and are versatile regarding atmospheric characteristics, which can enable identification of various planets in archival and future data.

The cosmic infrared background (CIB) traces star-forming galaxies throughout cosmic history, with emission peaking at $z\sim1-2$. CIB anisotropies are present at the far-infrared frequencies observed by cosmic microwave background (CMB) experiments such as $\textit{Planck}$. These contain a lot of astrophysical and cosmological information, but are hard to separate from the dust emission in our own Milky Way galaxy, especially on large scales where the Milky Way contamination severely dominates. This galactic component is often cleaned using information from other galactic tracers such as neutral hydrogen (HI). In this work we use HI data from the HI4PI survey to clean the 353, 545, and 857 GHz $\textit{Planck}$ NPIPE single-frequency maps using a needlet internal linear combination (NILC) method, with $\texttt{pyilc}$. This allows us to preserve the CIB anisotropy information on $\textit{all}$ scales, while reducing the variance sourced by the galactic contamination. We also create a NILC CMB map from the $\textit{Planck}$ NPIPE data, to subtract a CMB template from the 353 GHz map. Our resulting CIB maps are appropriate for cross-correlation studies with cosmological tracers such as CMB lensing maps down to very low $\ell$ ($\ell\sim10$), while achieving similar performance to previous works on intermediate scales. The use of the NPIPE data additionally allows us to achieve lower instrumental noise in the maps than in previous works. We use our maps, in combination with the $\textit{Planck}$ NPIPE CMB lensing reconstruction, to measure the CIB-CMB lensing cross correlation down to $\ell\sim10$. We make various versions of our maps publicly available to the community for further use in cross-correlation studies, along with a script (and the intermediate data products required) to produce dust-cleaned CIB maps on an arbitrary region of sky.

Contrastive learning (CL) has emerged as a potent tool for building meaningful latent representations of galaxy properties across a broad spectrum of wavelengths, ranging from optical and infrared to radio frequencies. These representations facilitate a variety of downstream tasks, including galaxy classification, similarity searches, and parameter estimation, which is why they are often referred to as foundation models. In this study, we employ CL on the latest extended DR from CALIFA survey, which encompasses 895 galaxies with enhanced spatial resolution. We demonstrate that CL can be applied to IFU surveys, even with small training sets, to meaningful embedding where galaxies are well-separated based on their physical properties. We discover that the strongest correlations in the embedding space are observed with the EW of Ha morphology, stellar metallicity, age, stellar surface mass density, the [NII]/Ha ratio, and stellar mass, in descending order of correlation strength. Additionally, we illustrate the feasibility of unsupervised separation of galaxy populations along the SFMS, successfully identifying the BC and the RS in a two-cluster scenario, and the GV population in a three-cluster scenario. Our findings indicate that galaxy luminosity profiles have minimal impact on the construction of the embedding space, suggesting that morphology and spectral features play a more significant role in distinguishing between galaxy populations. Moreover, we explore the use of CL for detecting variations in galaxy population distributions across different environments, including voids, clusters, filaments and walls. Nonetheless, we acknowledge the limitations of the CL and our specific training set in detecting subtle differences in galaxy properties, such as the presence of an AGN or other minor scale variations that exceed the scope of primary parameters like stellar mass or morphology.

Fast Radio Bursts (FRBs) are millisecond dispersed radio pulses of predominately extra-galactic origin. Although originally discovered at GHz frequencies, most FRBs have been detected between 400 to 800 MHz. Nevertheless, only a handful of FRBs were detected below 400 MHz. Searching for FRBs at low frequencies is computationally challenging due to increased dispersive delay that must be accounted for. However, the wide field of view (FoV) of low-frequency telescopes - such as the the Murchison Widefield Array (MWA), and prototype stations of the low-frequency Square Kilometre Array (SKA-Low) - makes them promising instruments to open a low-frequency window on FRB event rates, and constrain FRB emission models. The standard approach, inherited from high-frequencies, is to form multiple tied-array beams to tessellate the entire FoV and perform the search on the resulting time series. This approach, however, may not be optimal for low-frequency interferometers due to their large FoVs and high spatial resolutions leading to a large number of beams. Consequently, there are regions of parameter space in terms of number of antennas and resolution elements (pixels) where interferometric imaging is computationally more efficient. Here we present a new high-time resolution imager BLINK implemented on modern Graphical Processing Units (GPUs) and intended for radio astronomy. The main goal for this imager is to become part of a fully GPU-accelerated FRB search pipeline. We describe the imager and present its verification on real and simulated data processed to form all-sky and widefield images from the MWA and prototype SKA-Low stations. We also present and compare benchmarks of the GPU and CPU code executed on laptops, desktop computers, and Australian supercomputers. The code is publicly available at this https URL and can be applied to data from any radio telescope.

The enigmatic and intriguing phenomenon of the "soft excess" observed in the X-ray spectra of luminous quasars continues to be a subject of considerable interest and debate in the field of high-energy astrophysics. This study focuses on the quasar HE 1029-1401 ($z=0.086$, $\log(L_{\rm{bol}}/[\rm{erg\,s^{-1}}])= 46.0 \pm 0.2$), with a particular emphasis on investigating the properties of the hot corona and the physical origin of the soft excess. In this study, we present the results of a joint \textit{XMM-Newton}/\textit{NuSTAR} monitoring campaign of this quasar conducted in May 2022. The source exhibits a cold and narrow Fe $\rm{K}\alpha$ emission line at 6.4 keV, in addition to the detection of a broad component. Our findings suggest that the soft excess observed in HE 1029-1401 can be adequately explained by Comptonized emission originating from a warm corona. Specifically, fitting the spectra with two \nthcomp\, component we found that the warm corona is characterized by a photon index ($\Gamma^{w}$) of $2.75\pm0.05$ and by an electron temperature ($kT_{e}^{w}$) of $0.39^{+0.06}_{-0.04}$ keV, while the optical depth ($\tau^{w}$) is found to be $23\pm3$. We also test more physical models for the warm corona, corresponding to two scenarios: pure Comptonization and Comptonization plus reflection. Both models provide a good fit to the data, and are in agreement with a radially extended warm corona having a size of a few tens of gravitational radii.

Magnesium is one of the important elements in stellar physics as an electron donor and in Galactic Archaeology as a discriminator of different stellar populations. However, previous studies of Mg I and Mg II lines in metal-poor benchmark stars have flagged problems with magnesium abundances inferred from one-dimensional (1D), hydrostatic models of stellar atmospheres, both with or without the local thermodynamic equilibrium (LTE) approximation. We here present 3D non-LTE calculations for magnesium in FG-type dwarfs, and provide corrections for 1D LTE abundances. The 3D non-LTE corrections reduce the ionisation imbalances in the benchmark metal-poor stars HD84937 and HD140283 from $-0.16$ dex and $-0.27$ dex in 1D LTE, to just $-0.02$ dex and $-0.09$ dex respectively. We then applied our abundance corrections to 1D LTE literature results for stars in the thin disc, thick disc, $\alpha$-rich halo, and $\alpha$-poor halo. We find that the 3D non-LTE results show a richer substructure in [Mg/Fe]-[Fe/H] in the $\alpha$-poor halo, revealing two subpopulations at the metal-rich end. These two subpopulations are also separated in kinematics, supporting the astrophysical origin of the separation. While the more magnesium-poor subpopulation is likely to be debris from a massive accreted galaxy, Gaia-Enceladus, the other subpopulation may be related to a previous identified group of stars, called Eos. The presence of additional separation in [Mg/Fe] suggests that previous Mg abundance measurements may have been limited in the precision by the 1D and LTE approximations, highlighting the importance of 3D non-LTE modelling.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information.

The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated, and its technical potentials and limitations. Links to articles providing more details and technical background are included. NISP's 16 HAWAII-2RG (H2RG) detectors with a plate scale of 0.3" pix^-1 deliver a field-of-view of 0.57deg^2. In photo mode, NISP reaches a limiting magnitude of ~24.5AB mag in three photometric exposures of about 100s exposure time, for point sources and with a signal-to-noise ratio (SNR) of 5. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux ~2x10^-16erg/s/cm^2 integrated over two resolution elements of 13.4A, in 3x560s grism exposures at 1.6 mu (redshifted Ha). Our calibration includes on-ground and in-flight characterisation and monitoring of detector baseline, dark current, non-linearity, and sensitivity, to guarantee a relative photometric accuracy of better than 1.5%, and relative spectrophotometry to better than 0.7%. The wavelength calibration must be better than 5A. NISP is the state-of-the-art instrument in the NIR for all science beyond small areas available from HST and JWST - and an enormous advance due to its combination of field size and high throughput of telescope and instrument. During Euclid's 6-year survey covering 14000 deg^2 of extragalactic sky, NISP will be the backbone for determining distances of more than a billion galaxies. Its NIR data will become a rich reference imaging and spectroscopy data set for the coming decades.

The near-infrared calibration unit (NI-CU) onboard Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid is a mission in ESA's 'Cosmic Vision 2015-2025' framework, to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and dark energy. Calibrating photometric and spectrometric measurements of galaxies to better than 1.5% accuracy in a survey homogeneously mapping ~14000 deg^2 of extragalactic sky requires a very detailed characterisation of near-infrared (NIR) detector properties, as well their constant monitoring in flight. To cover two of the main contributions - relative pixel-to-pixel sensitivity and non-linearity characteristics - as well as support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous (<12% variations) and temporally stable illumination (0.1%-0.2% over 1200s) over the NISP detector plane, with minimal power consumption and energy dissipation. NI-CU is covers the spectral range ~[900,1900] nm - at cryo-operating temperature - at 5 fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of >=100 from ~15 ph s^-1 pixel^-1 to >1500 ph s^-1 pixel^-1. For this functionality, NI-CU is based on LEDs. We describe the rationale behind the decision and design process, describe the challenges in sourcing the right LEDs, as well as the qualification process and lessons learned. We also provide a description of the completed NI-CU, its capabilities and performance as well as its limits. NI-CU has been integrated into NISP and the Euclid satellite, and since Euclid's launch in July 2023 has started supporting survey operations.

We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from the combination of weak gravitational lensing and galaxy clustering data. The breath of Euclid's data will also foster a wide variety of scientific analyses. The Flagship simulation was developed to provide a realistic approximation to the galaxies that will be observed by Euclid and used in its scientific analyses. We ran a state-of-the-art N-body simulation with four trillion particles, producing a lightcone on the fly. From the dark matter particles, we produced a catalogue of 16 billion haloes in one octant of the sky in the lightcone up to redshift z=3. We then populated these haloes with mock galaxies using a halo occupation distribution and abundance matching approach, calibrating the free parameters of the galaxy mock against observed correlations and other basic galaxy properties. Modelled galaxy properties include luminosity and flux in several bands, redshifts, positions and velocities, spectral energy distributions, shapes and sizes, stellar masses, star formation rates, metallicities, emission line fluxes, and lensing properties. We selected a final sample of 3.4 billion galaxies with a magnitude cut of H_E<26, where we are complete. We have performed a comprehensive set of validation tests to check the similarity to observational data and theoretical models. In particular, our catalogue is able to closely reproduce the main characteristics of the weak lensing and galaxy clustering samples to be used in the mission's main cosmological analysis. (abridged)

The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level.

We provide an early assessment of the imaging capabilities of the Euclid space mission to probe deeply into nearby star-forming regions and associated very young open clusters, and in particular to check to what extent it can shed light on the new-born free-floating planet population. This paper focuses on a low-reddening region observed in just one Euclid pointing where the dust and gas has been cleared out by the hot sigma Orionis star. One late-M and six known spectroscopically confirmed L-type substellar members in the sigma Orionis cluster are used as benchmarks to provide a high-purity procedure to select new candidate members with Euclid. The exquisite angular resolution and depth delivered by the Euclid instruments allow us to focus on bona-fide point sources. A cleaned sample of sigma Orionis cluster substellar members has been produced and the initial mass function (IMF) has been estimated by combining Euclid and Gaia data. Our sigma Orionis substellar IMF is consistent with a power-law distribution with no significant steepening at the planetary-mass end. No evidence of a low-mass cutoff is found down to about 4 Jupiter masses at the young age (3 Myr) of the sigma Orionis open cluster.

As part of the Euclid Early Release Observations (ERO) programme, we analyse deep, wide-field imaging from the VIS and NISP instruments of two Milky Way globular clusters (GCs), namely NGC 6254 (M10) and NGC 6397, to look for observational evidence of their dynamical interaction with the Milky Way. We search for such an interaction in the form of structural and morphological features in the clusters' outermost regions, which are suggestive of the development of tidal tails on scales larger than those sampled by the ERO programme. Our multi-band photometric analysis results in deep and well-behaved colour-magnitude diagrams that, in turn, enable an accurate membership selection. The surface brightness profiles built from these samples of member stars are the deepest ever obtained for these two Milky Way GCs, reaching down to $\sim30.0$ mag~arcsec$^{-2}$, which is about $1.5$ mag arcsec$^{-2}$ below the current limit. The investigation of the two-dimensional density map of NGC 6254 reveals an elongated morphology of the cluster peripheries in the direction and with the amplitude predicted by $N$-body simulations of the cluster's dynamical evolution, at high statistical significance. We interpret this as strong evidence for the first detection of tidally induced morphological distortion around this cluster. The density map of NGC 6397 reveals a slightly elliptical morphology, in agreement with previous studies, which requires further investigation on larger scales to be properly interpreted. This ERO project thus demonstrates the power of Euclid in studying the outer regions of GCs at an unprecedented level of detail, thanks to the combination of large field of view, high spatial resolution, and depth enabled by the telescope. Our results highlight the future Euclid survey as the ideal data set to investigate GC tidal tails and stellar streams.

Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from about 0.5 Mpc to 8.8 Mpc. Our assessment of the surface brightness depths in the stacked Euclid images confirms previous estimates in 100 arcsec^2 regions of 1sigma=30.5 mag/arcsec^2 for VIS, but slightly deeper than previous estimates for NISP with 1sigma=29.2-29.4 mag/arcsec^2. By combining Euclid HE, YE, and IE into RGB images, we illustrate the large field-of-view covered by a single Reference Observing Sequence, together with exquisite detail on parsec scales in these nearby galaxies. Radial surface brightness and color profiles demonstrate galaxy colors in agreement with stellar population synthesis models. Standard stellar photometry selection techniques find approximately 1.3 million stars across the 6 galaxy fields. Euclid's resolved stellar photometry allows us to constrain the star-formation histories of these galaxies, by disentangling the distributions of young stars, as well as asymptotic giant branch and red giant branch stellar populations. We finally examine 2 galaxies individually for surrounding satellite systems. Our analysis of the ensemble of dwarf satellites around NGC6744 reveals a new galaxy, EDwC1, a nucleated dwarf spheroidal at the end of a spiral arm. Our new census of the globular clusters around NGC2403 yields 9 new star-cluster candidates, 8 of which with colors indicative of evolved stellar populations. In summary, our investigation of the 6 Showcase galaxies demonstrates that Euclid is a powerful probe of the anatomy of nearby galaxies [abridged].

We present an analysis of Euclid observations of a 0.5 deg$^2$ field in the central region of the Fornax galaxy cluster that were acquired during the performance verification phase. With these data, we investigate the potential of Euclid for identifying GCs at 20 Mpc, and validate the search methods using artificial GCs and known GCs within the field from the literature. Our analysis of artificial GCs injected into the data shows that Euclid's data in $I_{\rm E}$ band is 80% complete at about $I_{\rm E} \sim 26.0$ mag ($M_{V\rm } \sim -5.0$ mag), and resolves GCs as small as $r_{\rm h} = 2.5$ pc. In the $I_{\rm E}$ band, we detect more than 95% of the known GCs from previous spectroscopic surveys and GC candidates of the ACS Fornax Cluster Survey, of which more than 80% are resolved. We identify more than 5000 new GC candidates within the field of view down to $I_{\rm E}$ mag, about 1.5 mag fainter than the typical GC luminosity function turn-over magnitude, and investigate their spatial distribution within the intracluster field. We then focus on the GC candidates around dwarf galaxies and investigate their numbers, stacked luminosity distribution and stacked radial distribution. While the overall GC properties are consistent with those in the literature, an interesting over-representation of relatively bright candidates is found within a small number of relatively GC-rich dwarf galaxies. Our work confirms the capabilities of Euclid data in detecting GCs and separating them from foreground and background contaminants at a distance of 20 Mpc, particularly for low-GC count systems such as dwarf galaxies.

The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution.

We make use of the unprecedented depth, spatial resolution, and field of view of the Euclid Early Release Observations of the Perseus galaxy cluster to detect and characterise the dwarf galaxy population in this massive system. The Euclid high resolution VIS and combined VIS+NIR colour images were visually inspected and dwarf galaxy candidates were identified. Their morphologies, the presence of nuclei, and their globular cluster (GC) richness were visually assessed, complementing an automatic detection of the GC candidates. Structural and photometric parameters, including Euclid filter colours, were extracted from 2-dimensional fitting. Based on this analysis, a total of 1100 dwarf candidates were found across the image, with 638 appearing to be new identifications. The majority (96%) are classified as dwarf ellipticals, 53% are nucleated, 26% are GC-rich, and 6% show disturbed morphologies. A relatively high fraction of galaxies, 8%, are categorised as ultra-diffuse galaxies. The majority of the dwarfs follow the expected scaling relations. Globally, the GC specific frequency, S_N, of the Perseus dwarfs is intermediate between those measured in the Virgo and Coma clusters. While the dwarfs with the largest GC counts are found throughout the Euclid field of view, those located around the east-west strip, where most of the brightest cluster members are found, exhibit larger S_N values, on average. The spatial distribution of the dwarfs, GCs, and intracluster light show a main iso-density/isophotal centre displaced to the west of the bright galaxy light distribution. The ERO imaging of the Perseus cluster demonstrates the unique capability of Euclid to concurrently detect and characterise large samples of dwarfs, their nuclei, and their GC systems, allowing us to construct a detailed picture of the formation and evolution of galaxies over a wide range of mass scales and environments.

We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus cluster hosts 70000$\pm$2800 GCs and $1.6\times10^{12}$ L$_\odot$ of diffuse light from the BCG+ICL in the near-infrared H$_E$. This accounts for 37$\pm$6% of the cluster's total stellar luminosity within this radius. The ICL and ICGCs share a coherent spatial distribution, suggesting a common origin or that a common potential governs their distribution. Their contours on the largest scales (>200 kpc) are offset from the BCG's core westwards by 60 kpc towards several luminous cluster galaxies. This offset is opposite to the displacement observed in the gaseous intracluster medium. The radial surface brightness profile of the BCG+ICL is best described by a double Sérsic model, with 68$\pm$4% of the H$_E$ light in the extended, outer component. The transition between these components occurs at ~50 kpc, beyond which the isophotes become increasingly elliptical and off-centred. The radial ICGC number density profile closely follows the BCG+ICL profile only beyond this 50 kpc radius, where we find an average of 60 GCs per $10^9$ M$_\odot$ of diffuse stellar mass. The BCG+ICL colour becomes increasingly blue with radius, consistent with the stellar populations in the ICL having subsolar metallicities [Fe/H]~-0.6. The colour of the ICL, and the specific frequency and luminosity function of the ICGCs suggest that the ICL+ICGCs were tidally stripped from the outskirts of massive satellites with masses of a few $\times10^{10}$ M$_\odot$, with an increasing contribution from dwarf galaxies at large radii.

We present the first analysis of the Euclid Early Release Observations (ERO) program that targets fields around two lensing clusters, Abell 2390 and Abell 2764. We use VIS and NISP imaging to produce photometric catalogs for a total of $\sim 500\,000$ objects. The imaging data reach a $5\,\sigma$ typical depth in the range 25.1-25.4 AB in the NISP bands, and 27.1-27.3 AB in the VIS band. Using the Lyman-break method in combination with photometric redshifts, we identify $30$ Lyman-break galaxy (LBG) candidates at $z>6$ and 139 extremely red sources (ERSs), most likely at lower redshift. The deeper VIS imaging compared to NISP means we can routinely identify high-redshift Lyman breaks of the order of $3$ magnitudes, which reduces contamination by brown dwarf stars and low-redshift galaxies. Spectroscopic follow-up campaigns of such bright sources will help constrain both the bright end of the ultraviolet galaxy luminosity function and the quasar luminosity function at $z>6$, and constrain the physical nature of these objects. Additionally, we have performed a combined strong lensing and weak lensing analysis of A2390, and demonstrate how Euclid will contribute to better constraining the virial mass of galaxy clusters. From these data, we also identify optical and near-infrared counterparts of known $z>0.6$ clusters, which exhibit strong lensing features, establishing the ability of Euclid to characterize high-redshift clusters. Finally, we provide a glimpse of Euclid's ability to map the intracluster light out to larger radii than current facilities, enabling a better understanding of the cluster assembly history and mapping of the dark matter distribution. This initial dataset illustrates the diverse spectrum of legacy science that will be enabled by the Euclid survey.

This paper presents a search for high redshift galaxies from the Euclid Early Release Observations program "Magnifying Lens." The 1.5 deg$^2$ area covered by the twin Abell lensing cluster fields is comparable in size to the few other deep near-infrared surveys such as COSMOS, and so provides an opportunity to significantly increase known samples of rare UV-bright galaxies at $z\approx6-8$ ($M_{\rm UV}\lesssim-22$). Beyond their still uncertain role in reionisation, these UV-bright galaxies are ideal laboratories from which to study galaxy formation and constrain the bright-end of the UV luminosity function. Of the 501994 sources detected from a combined $Y_{\rm E}$, $J_{\rm E}$, and $H_{\rm E}$ NISP detection image, 168 do not have any appreciable VIS/$I_{\rm E}$ flux. These objects span a range in spectral colours, separated into two classes: 139 extremely red sources; and 29 Lyman-break galaxy candidates. Best-fit redshifts and spectral templates suggest the former is composed of both $z\gtrsim5$ dusty star-forming galaxies and $z\approx1-3$ quiescent systems. The latter is composed of more homogeneous Lyman break galaxies at $z\approx6-8$. In both cases, contamination by L- and T-type dwarfs cannot be ruled out with Euclid images alone. Additional contamination from instrumental persistence is investigated using a novel time series analysis. This work lays the foundation for future searches within the Euclid Deep Fields, where thousands more $z\gtrsim6$ Lyman break systems and extremely red sources will be identified.

We present the results of a study aimed at characterizing the kinematics of the inner regions of the halo globular cluster M75 (NGC 6864) based on data acquired as part of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters. Our analysis includes the first determination of the line-of-sight velocity dispersion profile in the core region of M75. By using MUSE/NFM observations, we obtained a sample of $\sim 1900$ radial velocity measurements from individual stars located within $16''$ (corresponding to about $r < 3 r_c$ where $r_c$ is the estimated core radius of the system) from the cluster center. After an appropriate selection of the most accurate velocity measures, we determined the innermost portion of the velocity dispersion profile, finding that it is characterized by a constant behavior and a central velocity dispersion of $\sigma_0\sim 9$ km s$^{-1}$. The simultaneous King model fitting to the projected velocity dispersion and density profiles allowed us to check and update previous determinations of the main structural parameters of the system. We also detected a mild hint of rotation in the central $\sim 7''$ from the center, with an amplitude of just $\sim 1.0$ km s$^{-1}$ and a position angle of the rotation axis of PA$_0 = 174°$. Intriguingly, the position angle is consistent with that previously quoted for the suspected rotation signal in the outer region of the cluster. Taking advantage of the high quality of the photometric catalog used for the analysis of the MUSE spectra, we also provide updated estimates of the cluster distance, age, and reddening.

Baryon acoustic oscillation data from the first year of the Dark Energy Spectroscopic Instrument (DESI) provide near percent-level precision of cosmic distances in seven bins over the redshift range $z=0.1$-$4.2$. We use this data, together with other distance probes, to constrain the cosmic expansion history using some well-motivated physical classes of dark energy. In particular, we explore three physics-focused behaviors of dark energy from the equation of state and energy density perspectives: the thawing class (matching many simple quintessence potentials), emergent class (where dark energy comes into being recently, as in phase transition models), and mirage class (where phenomenologically the distance to CMB last scattering is close to that from a cosmological constant $\Lambda$ despite dark energy dynamics). All three classes fit the data at least as well as $\Lambda$CDM, and indeed can improve on it by $\Delta\chi^2\approx -5$ to $-17$ for the combination of DESI BAO with CMB and supernova data, while having one more parameter. The mirage class does essentially as well as $w_0w_a$CDM while having one less parameter. These classes of dynamical behaviors highlight worthwhile avenues for further exploration into the nature of dark energy.

We present a proof-of-principle implementation of the first fully covariant filtering scheme applied to relativistic fluid turbulence. The filtering is performed with respect to special observers, identified dynamically as moving with the "bulk of the flow". This means that filtering does not depend on foliations of spacetime but rather on the intrinsic fibration traced out by the observers. The covariance of the approach means that the results may be lifted into an arbitrary, curved spacetime. This practical step follows theoretical work showing that the residuals introduced by filtering a fine-scale ideal fluid can be represented by a non-ideal fluid prescription at the coarse scale. We interpret such non-ideal terms using a simple first-order gradient model, which allows us to extract effective turbulent viscosities and conductivity. A statistical regression on these terms shows that the majority of their variation may be explained by the thermodynamic properties of the filtered fluid and invariants of its flow, such as the shear and vorticity. This serves as a validation of the method and enables us to fit a functional, power-law form for the non-ideal coefficients -- an approach that may be used practically to give a sub-grid closure model in large-eddy simulations.

Most stripped envelope supernova progenitors are formed through binary interaction, losing hydrogen and/or helium from their outer layers. An emerging class of supernovae with the highest degree of envelope-stripping are thought to be the product of stripping by a NS companion. However, relatively few examples are known and the outcomes of such systems can be diverse and are poorly understood at present. Here, we present spectroscopic observations and high cadence multi-band photometry of SN 2023zaw, a low ejecta mass and rapidly evolving supernova. SN 2023zaw was discovered in a nearby spiral galaxy at D = 39.7 Mpc, with significant Milky Way extinction, $E(B-V) = 0.21$, and significant (but uncertain) host extinction. Bayesian evidence comparison reveals that nickel is not the only power source and an additional energy source is required to explain our observations. Our models suggest an ejecta mass of $M_{\rm ej} \sim 0.07\,\rm M_\odot$ and a synthesised nickel mass of $M_{\rm ej} \sim 0.007\,\rm M_\odot$ is required to explain the explosion. However an additional heating from a magnetar or interaction with circumstellar material is required to power the early light curve.

We propose a novel control approach that combines offline supervised learning to address the challenges posed by non-linear phase reconstruction using unmodulated pyramid wavefront sensors (P-WFS) and online reinforcement learning for predictive control. The control approach uses a high-order P-WFS to drive a tip-tilt stage and a high-dimensional mirror concurrently. Simulation results demonstrate that our method outperforms traditional control techniques, showing significant improvements in performance under challenging conditions such as faint stars and poor seeing, and exhibits robustness against variations in atmospheric conditions.

To detect additional bodies in binary systems, we performed a potent approach of orbital period variation analysis. In this work, we present 90 new mid-eclipse times of a short-period eclipsing binary system. Observations were made using two telescopes from 2014 to 2024, extending the time span of the O-C diagram to 24 years. The data obtained in the last seven years indicate significant deviations in the O-C diagram from the models obtained in previous studies. We investigated whether this variation could be explained by mechanisms such as the LTT effect or Applegate. To investigate the cyclic behaviour observed in the system with the LTT effect, we modelled the updated O-C diagram using different models including linear/quadratic terms and additional bodies. The updated O-C diagram is statistically consistent with the most plausible solutions of models that include multiple brown dwarfs close to each other. However, it has been found that the orbit of the system is unstable on short time scales. Using three different theoretical definitions, we have found that the Applegate mechanism cannot explain the variation in the orbital period except for the model containing the fifth body. Therefore, due to the complex nature of the system, further mid-eclipse time is required before any conclusions can be drawn about the existence of additional bodies.

Over recent decades, astronomy has entered the era of massive data and real-time surveys. This is improving the study of transient objects - although they still contain some of the most poorly understood phenomena in astrophysics, as it is inherently more difficult to obtain data on them. In order to help detect these objects in their brightest state, we have built a quasi-real time transient detection system for the XMM-Newton pipeline: the Search for Transient Objects in New detections using Known Sources (STONKS) pipeline. STONKS detects long-term X-ray transients by automatically comparing new XMM-Newton detections to any available archival X-ray data at this position, sending out an alert if the amplitude of variability between observations is over 5. This required an initial careful cross-correlation and flux calibration of various X-ray catalogs from different observatories (XMM-Newton, Chandra, Swift, ROSAT, and eROSITA). We also systematically computed the XMM-Newton upper limits at the position of any X-ray source covered by the XMM-Newton observational footprint, even without any XMM-Newton counterpart. The behavior of STONKS was then tested on all 483 observations performed with imaging mode in 2021. Over the 2021 testing run, STONKS provided $0.7^{+0.7}_{-0.5}$ alerts per day, about 80% of them being serendipitous. STONKS also detected targeted tidal disruption events, ensuring its ability to detect other serendipitous events. As a byproduct of our method, the archival multi-instrument catalog contains about one million X-ray sources, with 15% of them involving several catalogs and 60% of them having XMM-Newton upper limits. STONKS demonstrates a great potential for revealing future serendipitous transient X-ray sources, providing the community with the ability to follow-up on these objects a few days after their detection.

The Gaia-Sausage-Enceladus merger was a major event in the history of the Milky Way. Studies on Milky Way satellite dwarf galaxies show that key elemental abundance patterns, which probe different nucleosynthetic channels, reflect the host galaxy's star formation history. We gather Mg, Fe, Ba, and Eu abundance measurements for Gaia-Sausage-Enceladus stars from the SAGA database and use [Fe/Mg], [Ba/Mg], [Eu/Mg], and [Eu/Ba], as a function of [Fe/H] to constrain the star formation history of Gaia-Sausage-Enceladus. We use the known star formation histories and elemental abundance patterns of the Sculptor and Fornax dwarf spheroidal galaxies as comparison. The elemental abundance ratios of [Fe/Mg], [Ba/Mg], [Eu/Mg], and [Eu/Ba] all increase with [Fe/H] in Gaia-Sausage- Enceladus. The [Eu/Mg] begins to increase at [Fe/H]= -2.0 and continues steadily, contrasting with the Sculptor dSph galaxy. The [Eu/Ba] increases and remains high across the [Fe/H] range, contrasting with that of the Sculptor dSph galaxy and deviating from the Fornax dSph galaxy at high [Fe/H]. The [Ba/Mg] is higher than those of the Sculptor dSph galaxy at the lowest [Fe/H] and gradually increases, similar to the Fornax dSph galaxy. We constrain three main properties of the Gaia-Sausage-Enceladus star formation history: 1) star formation started gradually, 2) it extended for over 2 Gyr, and 3) it was quenched around [Fe/H] of -0.5, likely when it fell into the Milky Way.

We aim to identify QSO2 candidates in the redshift desert using optical and infrared photometry. At this intermediate redshift range, most of the prominent optical emission lines in QSO2 sources (e.g. CIV1549; [OIII]4959,5008) fall either outside the wavelength range of the SDSS optical spectra or in particularly noisy wavelength ranges, making QSO2 identification challenging. Therefore, we adopted a semi-supervised machine learning approach to select candidates in the SDSS galaxy sample. Recent applications of machine learning in astronomy focus on problems involving large data sets, with small data sets often being overlooked. We developed a few-shot learning approach for the identification and classification of rare-object classes using limited training data (200 sources). The new AMELIA pipeline uses a transfer-learning based approach with decision trees, distance-based, and deep learning methods to build a classifier capable of identifying rare objects on the basis of an observational training data set. We validated the performance of AMELIA by addressing the problem of identifying QSO2s at 1 $\leq$ z $\leq$ 2 using SDSS and WISE photometry, obtaining an F1-score above 0.8 in a supervised approach. We then used AMELIA to select new QSO2 candidates in the redshift desert and examined the nature of the candidates using SDSS spectra, when available. In particular, we identified a sub-population of [NeV]3426 emitters at z $\sim$ 1.1, which are highly likely to contain obscured AGNs. We used X-ray and radio cross-matching to validate our classification and investigated the performance of photometric criteria from the literature showing that our candidates have an inherent dusty nature. Finally, we derived physical properties for our QSO2 sample using photoionisation models and verified the AGN classification using an SED fitting.

Following the widespread practice of exoplanetary transit simulations, various presumed components of an extrasolar system can be examined in numerically simulated transits, including exomoons, rings around planets, and the deformation of exoplanets. Template signals can then be used to efficiently search for light curve features that mark specific phenomena in the data, and they also provide a basis for feasibility studies of instruments and search programs. In this paper, we present a method for exocomet transit light curve calculations using arbitrary dust distributions in transit. The calculations, spanning four distinct materials (carbon, graphite, pyroxene, and olivine), dust grain sizes ($100$\,nm -- $300$\,nm, $300$\,nm -- $1000$\,nm, and $1000$\,nm -- $3000$\,nm) encompass light curves in VRJHKL bands. We also investigated the behavior of scattering colors. We show that multicolor photometric observations are highly effective tools in the detection and characterization of exocomet transits. They provide information on the dust distribution of the comet (encoded in the light curve shape), while the color information itself can reveal the particle size change and material composition of the transiting material, in relation to the surrounding environment. We also show that the typical cometary tail can result in the wavelength dependence of the transit timing. We demonstrate that multi-wavelength observations can yield compelling evidence for the presence of exocomets in real observations.

In this work, we obtain Hubble constant ($H_0$) estimates by using two galaxy cluster gas mass fraction measurement samples, Type Ia supernovae luminosity distances and the validity of the cosmic distance duality relation. Notably, the angular diameter distance (ADD) to each galaxy cluster in the samples is determined by combining its gas mass fraction measurement with galaxy clustering observations, more precisely, the $\Omega_b / \Omega_m$ ratio. Such a combination results in a $H_0$ estimate that is independent of a specific cosmological framework. In one of the samples, the gas fraction measurements were calculated in spherical shells at radii near $r_{\rm 2500}$ (44 data points), while in the other (103 data points) the measurements were obtained by integrating X-ray observations across all radii $< r_{\rm 500}$. The quantities $r_{2500}$ and $r_{\rm 500}$ define the radii within which the mean cluster density is 2500 or 500 times the critical density of the Universe at the cluster's redshift. We find $H_0=71.4^{+10.5}_{-8.9}$ km/s/Mpc and $H_0=75.5^{+9.9}_{-7.9}$ km/s/Mpc at 68\% CL, respectively. We also investigate the impact on the $H_0$ determination by exploring the precision and number of gas mass fraction data by performing a data Monte Carlo simulation. Our simulations show that future measurements could achieve a precision of up to 3.5\% for $H_0$.

NGC 777 provides an example of a phenomenon observed in some group-central ellipticals, in which the temperature profile shows a central peak, despite the short central cooling time of the intra-group medium. We use deep Chandra X-ray observations of the galaxy, supported by uGMRT 400 MHz radio imaging, to investigate the origin of this hot core. We confirm the centrally-peaked temperature profile and find that entropy and cooling time both monotonically decline to low values (2.62 [+0.19, -0.18] keV cm$^2$ and 71.3 [+12.8, -13.1] Myr) in the central ~700 pc. Faint diffuse radio emission surrounds the nuclear point source, with no clear jets or lobes but extending to ~10 kpc on a northwest-southeast axis. This alignment and extent agree well with a previously identified filamentary H$\alpha$+[NII] nebula. While cavities are not firmly detected, we see X-ray surface brightness decrements on the same axis at 10-20 kpc radius which are consistent with the intra-group medium having been pushed aside by expanding radio lobes. Any such outburst must have occurred long enough ago for lobe emission to have faded below detectability. Cavities on this scale would be capable of balancing radiative cooling for at least ~240 Myr. We consider possible causes of the centrally peaked temperature profile, including gravitational heating of gas as the halo relaxes after a period of AGN jet activity, and heating by particles leaking from the remnant relativistic plasma of the old radio jets.

Polycyclic Aromatic Hydrocarbons (PAHs) are organic molecules which comprise the smallest particles of dust in the interstellar medium (ISM). Due to their broad/complex emission profiles, obtaining kinematics is a challenge with traditional methods, especially before the advent of the JWST. In this work we employ Principal Component Analysis (PCA) tomography to analyse JWST/NIRSpec IFU data of 3 nearby luminous infrared galaxies (LIRGs), namely, NGC 3256 N, NGC 3256 S and NGC 7469. We detect the signature of rotation in the second principal component of the 3.3 $\mu$m PAH feature in all three targets. We construct velocity maps from the principal components for the 3.3 $\mu$m PAH feature, Br$\beta$ (2.625$\mu$m) and molecular hydrogen, H$_2$ 1-0 S(1) (2.12$\mu$m). We find that in each target, the PAHs qualitatively follow the rotation of the galaxy, consistent with the rotational signature derived from Br$\beta$ and H$_2$. There are however some differences, with the PAH rotation in NGC 3256 N appearing at a different position angle, which suggest differences in the motion of the dust as compared to the gas. This kind of analysis opens a new window into this key component of the ISM.

Filament eruptions often lead to coronal mass ejections (CMEs) on the Sun and are one of the most energetic eruptive phenomena in the atmospheres of other late-type stars. However, the detection of filament eruptions and CMEs on stars beyond the solar system is challenging. Here we present six filament eruption cases on the Sun and show that filament material obscuring part of the solar disk can cause detectable dimming signatures in sun-as-a-star flux curves of He II 304 A. Those filament eruptions have similar morphological features, originating from small filaments inside active regions and subsequently strongly expanding to obscure large areas of the solar disk or the bright flare regions. We have tracked the detailed evolution of six obscuration dimmings and estimated the dimming properties, such as dimming depths, dimming areas, and duration. The largest dimming depth among the six events under study is 6.2% accompanied by the largest dimming area of 5.6\% of the solar disk area. Other events have maximum dimming depths in a range of around 1% to 3% with maximum areas varying between about 3% to 4% of the solar disk area. The duration of the dimming spans from around 0.4 hours to 7.0 hours for the six events under study. A positive correlation was found between the dimming depth and area, which may help to set constraint on the filament sizes in stellar observations.

The standard cosmological model, with its six independent parameters, has proven remarkably successful in describing the evolution of the Universe. One of these parameters, the optical depth to reionization $\tau_{\mathrm{reio}}$, represents the scatterings that Cosmic Microwave Background (CMB) photons will experience after decoupling from the primordial plasma as the intergalactic medium transitions from neutral to ionized. $\tau_{\mathrm{reio}}$ depends on the neutral hydrogen fraction $x_{\mathrm{HI}}$, which, in turn, should theoretically depend on cosmology. We present a novel method to establish the missing link between cosmology and reionization timeline. We discover the timeline has a universal shape well described by the Gompertz mortality law, applicable to any cosmology within our simulated data. This enables us to map cosmology to reionization using symbolic regression and to treat $\tau_{\mathrm{reio}}$ as a derived parameter. Reanalyzing CMB data with our universal $x_{\mathrm{HI}}$ tightens the constraint on $\tau_{\mathrm {reio}}$ by more than one order of magnitude to $\approx 1\%$ and reduced the error on the amplitude of the primordial fluctuations by a factor of 2.5 compared to Planck's PR3 constraint. While our results rely on the astrophysical assumptions in our simulations, the methodology itself is independent of these assumptions; after all, $\tau_{\mathrm{reio}}$ is fundamentally a function of cosmology.

We investigate the effect of spin on equal and unequal mass binary neutron star mergers using finite-temperature, composition-dependent Hempel-Schaffner-Bielich (SFHo) equation of state, via 3+1 general relativistic hydrodynamics simulations which take into account neutrino emission and absorption. Equal mass cases that have a mass of $M_{1,2}$ =$1.27$, $1.52$ and $2.05M_{\odot}$, result in a supramassive neutron star, a delayed, and a prompt collapse to a black hole, respectively. For all cases, we analyse the effect of initial spin on dynamics, on the structure of the final remnant, its spin evolution, the amount and composition of the ejected matter, gravitational waves, neutrino energy and luminosities. We show that in equal mass binary neutron star mergers, the ejected mass could reach $\sim0.085M_{\odot}$ for highly aligned-spins ($\chi=0.67$). The black hole which results from such highly spinning, high-mass binary neutron star merger reaches a dimensionless spin of $0.92$; this is the highest spin reached in binary neutron star mergers, to date.

Low-loss deposited dielectrics are beneficial for the advancement of superconducting integrated circuits for astronomy. In the microwave band ($\mathrm{\sim}$1-10 GHz) the cryogenic and low-power dielectric loss is dominated by two-level systems. However, the origin of the loss in the millimeter-submillimeter band ($\mathrm{\sim}$0.1-1 THz) is not understood. We measured the loss of hydrogenated amorphous SiC (a-SiC:H) films in the 0.27-100 THz range using superconducting microstrip resonators and Fourier-transform spectroscopy. The agreement between the loss data and a Maxwell-Helmholtz-Drude dispersion model suggests that vibrational modes above 10 THz dominate the loss in the a-SiC:H above 200 GHz.

While much has been learned in recent decades about the X-ray emission of the extragalactic Ultraluminous X-ray sources (ULXs), their radiative output in the UV band remains poorly constrained. Understanding of the full ULX spectral energy distribution (SED) is imperative to constrain the accretion flow geometry powering them, as well as their radiative power. Here we present constraints on the UV emission of the pulsating ULX (PULX) NGC~1313~X--2 based on the absence of nebular {He{\sc ii} $\lambda$4686} emission in its immediate environment. To this end, we first perform multi-band spectroscopy of the ULX to derive three realistic extrapolations of the SED into the unaccessible UV, each predicting varying levels of UV luminosity. We then perform photo-ionization modelling of the bubble nebula and predict the {He{\sc ii} $\lambda$4686} fluxes that should have been observed based on each of the derived SEDs. We then compare these predictions with the derived upper limit on {\heii} from MUSE data, which allows us to infer a UV luminosity $L_\mathrm{UV} \lesssim 1 \times 10^{39}$ erg/s in the PULX NGC~1313~X--2. Comparing the UV luminosity inferred with other ULXs, our work suggests there may be an intrinsic difference between hard and soft ULXs, either related to different mass-transfer rates and/or the nature of the accretor. However, a statistical sample of ULXs with inferred UV luminosities is needed to fully determine the distinguishing features between hard and soft ULXs. Finally, we discuss ULXs ionising role in the context of the nebular {He{\sc ii} $\lambda$4686} line observed in star-forming, metal-poor galaxies.

The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, hence providing important anchors for planet formation models. We aim to improve the characterisation of the architecture of this key system, and in particular the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we aim to use it as a test bench for the robustness of the planetary mass determination with each technique. We perform a global analysis of all available photometry and radial velocity. We also try different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by each method. We show how stellar activity is preventing us from obtaining a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis resulted in a robust mass determination for planets c to g, with precision of 12% for the mass of planet c, and better than 10% for planets d to g. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable during the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain.

Irregularity in spectrum of the primary cosmic rays (PCR) mass composition at an energy of ~ 10 PeV is considered. To assess changes of the PCR mass composition, the X-ray emulsion chamber (XREC) method and the halo-method based on the XEC were used. The study of changes in the PCR mass composition was carried out based on the experimentally obtained characteristics of the extensive air showers (EAS) trunks, where fluctuations of these characteristics are minimal and information about the primary interaction of PCR nuclei with atmospheric atoms is maximally preserved. It is shown that at an energy of ~ 10 PeV there is a local maximum in the heavy nuclei fraction. This maximum is correlated with the PCR sources, - variable stars SR (red giants and super-giants) and WR (Wolf-Rayet).

We studied blue horizontal branch stars (BHBs), and calculated their radial velocities. Spectra of these stars have been obtained with moderate signal-to-noise ratio for five blue horizontal-branch stars using the 2 meter telescope and Echelle Spectrograph in Ondrejov observatory, Czech republic.

The chemical composition of a star's atmosphere reflects the chemical composition of its birth environment. Therefore, it should be feasible to recognize stars born together that have scattered throughout the galaxy, solely based on their chemistry. This concept, known as "strong chemical tagging", is a major objective of spectroscopic studies, but has yet to yield the anticipated results. We assess the existence and the robustness of the relation between chemical abundances and birth place using known member stars of open clusters. We followed a supervised machine learning approach, using chemical abundances obtained from APOGEE DR17, observed open clusters as labels and different data preprocessing techniques. We found that open clusters can be recovered with any classifier and on data whose features are not carefully selected. In the sample with no field stars, we obtain an average accuracy of $75.2\%$ and we find that the prediction accuracy depends mostly on the uncertainties of the chemical abundances. When field stars outnumber the cluster members, the performance degrades. Our results show the difficulty of recovering birth clusters using chemistry alone, even in a supervised scenario. This clearly challenges the feasibility of strong chemical tagging. Nevertheless, including information about ages could potentially enhance the possibility of recovering birth clusters.

Primordial black holes still represent a viable candidate for a significant fraction, if not for the totality, of dark matter. If these compact objects have masses of order tens of solar masses, their coalescence can be observed by current and future ground-based gravitational wave detectors. Therefore, finding new gravitational wave signatures associated with this dark matter candidate can either lead to their detection or help constraining their abundance. In this work we consider the phenomenology of primordial black holes in dense environments, in particular globular clusters. We model the internal structure of globular clusters in a semi-analytical fashion, and we derive the expected merger rate. We show that, if primordial black holes are present in globular clusters, their contribution to the GW background can be comparable to other well-known channels, such as early- and late-time binaries, thus enhancing the detectability prospects of primordial black holes and demonstrating that this contribution needs to be taken into account.

We demonstrate that "natural inflation", also known as "axion inflation", can be compatible with Planck 2018 measurements of the cosmic microwave background, while predicting an exponentially small tensor-to-scalar ratio, e.g., $r\sim 10^{-15}$. The strong suppression of $r$ arises from dynamics of the radial component of the complex scalar field, whose phase is the axion. Such tiny values of $r$ remain well below the threshold for detection by CMB-S4 or Simons Observatory B-mode searches. The model is testable with the running $\alpha_s$ of the spectral index, which is within reach of next-generation CMB and large-scale structure experiments, motivating the running as a primary science goal for future experiments.

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. We describe and test a "multi-zone" technique which is designed to tackle this difficult problem in 3D general relativistic magnetohydrodynamic (GRMHD) simulations. We simulate accretion on a non-spinning SMBH ($a_*=0$) using initial conditions from a large scale galaxy simulation, and achieve steady state over 8 decades in radius. The density scales with radius as $\rho \propto r^{-1}$ inside the Bondi radius $R_B$, which is located at $R_B=2\times 10^5 \,r_g$ ($\approx 60\,{\rm pc}$ for M87) where $r_g$ is the gravitational radius of the SMBH; the plasma-$\beta\sim$ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is $\approx 1\%$ of the analytical Bondi accretion rate $\dot{M}_B$; and there is continuous energy feedback out to $\approx 100R_B$ (or beyond $>\,{\rm kpc}$) at a rate $\approx 0.02 \dot{M}c^2$. Surprisingly, any ordered rotation in the external medium does not survive as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multi-zone method, we simulate GRMHD accretion for a wide range of Bondi radii, $R_{\rm B} \sim 10^2 - 10^7\,r_{\rm g}$, and find that $\dot{M}/\dot{M}_B\approx (R_B/6\, r_g)^{-0.5}$.

PSR J2215+5135 (J2215) is a `redback' spider pulsar, where the intrabinary shock (IBS) wraps around the pulsar rather than the stellar-mass companion. Spider orbital light curves are modulated, dominated by their binary companion thermal emission in the optical band and by IBS synchrotron emission in the X-rays. We report on new XMM-Newton X-ray and U-band observations of J2215. We produce orbital light curves and use them to model the system properties. From our optical light curve modeling, we find a neutron star mass $M_{NS}=2.15\pm0.10 M_\odot$, somewhat lower than previously reported. From the X-ray analysis, we find that the IBS still wraps around the pulsar, but with a wind/companion wind momenta ratio unusually close to unity, implying a flatter IBS geometry {than seen in other spiders. Estimating the companion wind momentum and speed from the X-ray light curve, we find a companion mass-loss rate of ${\dot M}_c\gtrsim10^{-10}$ M$_\odot$ yr$^{-1}$, so that J2215 may become an isolated millisecond pulsar in $\sim 1$ Gyr. Our X-ray analyses support models of magnetic reconnection and particle acceleration in the highly magnetized relativistic IBS.

We compared magnetograms from the KPVT/SPM, SoHO/MDI, SOLIS/VSM, and SDO/HMI with the aim of probing the effect on measured solar magnetism of the variation in instrument response with time, magnetogram signal level, and position on the solar disc. Taking near-simultaneous observations from the various instruments, we examined the surface coverage by magnetic activity and the effect of cross-calibrating the various instruments under different assumptions. By comparing the surface coverage by magnetic activity in the observations from the various instruments, we traced the effect of the time variation in instrument response on the longitudinal magnetogram signal and disc-integrated unsigned magnetic flux. This yielded evidence of acute changes in the response of MDI and VSM with certain events such as the SoHO vacation in 1998 and the upgrade of the VSM CCD camera in 2009. Excluding these changes, the effect of instrument instability on the magnetogram signal and disc-integrated magnetic flux appears to be rather benign, with an associated uncertainty of less than 2%. We determined the magnetogram signal ratio between each instrument pairing as a function of magnetogram signal level and distance from disc centre and with it cross-calibrated the various instruments. We compared the result with that from repeating the cross-calibration with the overall magnetogram signal ratio. This allowed us to estimate the uncertainty in the magnetogram signal associated with the variation in instrument response with magnetogram signal level and distance from disc centre to be about 8% to 14%. The corresponding uncertainty in the disc-integrated magnetic flux is about 7% to 23 %. The results here will be useful to the interpretation of SPM, MDI, VSM, and HMI magnetograms. As examples, we applied our findings to selected results from earlier studies based on such data.

We derive the nitrogen and oxygen abundances in the Narrow Line Regions (NLRs) of a sample of 38 local ($z \: < \: 0.4$) Seyfert~2 nuclei. For that, we consider narrow optical emission line intensities and direct estimates of the electron temperatures ($T_{\rm e}$-method). We find nitrogen abundances in the range $7.6 \: < \: \rm 12+log(N/H) \: < \: 8.6$ (mean value $8.06\pm0.22$) or $\rm 0.4 \: < \: (N/N_{\odot}) \: < 4.7$, in the metallicity regime $8.3 \: < \: \rm 12+log(O/H) \: < \: 9.0$. Our results indicate that the dispersion in N/H abundance for a fixed O/H value in AGNs is in agreement with that for disc \ion{H}{ii} regions with similar metallicity. We show that Seyfert~2 nuclei follow a similar (N/O)-(O/H) relation to the one followed by star-forming objects. Finally, we find that active galaxies called as 'nitrogen-loud' observed at very high redshift ($z \: > \: 5$) show N/O values in consonance with those derived for local NLRs. This result indicates that the main star-formation event is completed in the early evolution stages of active galaxies.

Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics.

Imidogen (NH) is a reactive molecule whose presence in astrochemical environments is of interest due to its role in the formation of nitrogen-containing molecules and as a potential probe of nitrogen abundance. Spectroscopic NH monitoring is useful for Earth-based combustion and photolysis processes of ammonia and other nitrogen-containing species. NH is also relevant to ultracold molecular physics and plasma studies. To enable these diverse applications, high-quality molecular spectroscopic data is required. Here, a new line list with significant advantages over existing data is presented. Most notably, this line list models isotopologue spectroscopy and forbidden transitions (important for NH visible absorption), alongside some overall improvements to accuracy and completeness. This approach takes advantage of existing experimental data (from a previous MARVEL compilation) and perturbative line lists together with new MRCI ab initio electronic data. These are used to produce a novel variational spectroscopic model and trihybrid line list for the main 14N1H isotopologue, as well as isotopologue-extrapolated hybrid line lists for the 14N2H, 15N1H and 15N2H isotopologues. The new 14N1H ExoMol-style trihybrid line list, kNigHt, comprises 4,076 energy levels (1,078 experimental) and 327,014 transitions up to 47,500 cm-1 (211 nm) between five low-lying electronic states (X 3{\Sigma}-, a 1{\Delta}, b 1{\Sigma}+, A 3{\Pi} and c 1{\Pi}). For most anticipated applications aside from far-IR studies, this line list will be of sufficient quality; any improvements should focus on the b 1{\Sigma}+ energies, and the a 1{\Delta} - A 3{\Pi} and b 1{\Sigma}+ - A 3{\Pi} spin-orbit couplings.

The ubiquity of "peas-in-a-pod" architectural patterns and the existence of the radius valley each present a striking population-level trend for planets with $R_{p} \leq 4 R_{\oplus}$ that serves to place powerful constraints on the formation and evolution of these subgiant worlds. As it has yet to be determined whether the strength of this peas-in-a-pod uniformity differs on either side of the radius valley, we separately assess the architectures of systems containing only small ($R_{p} \leq 1.6 R_{\oplus}$), rocky planets from those harboring only intermediate-size ($1.6 R_{\oplus} < R_{p} \leq 4 R_{\oplus}$), volatile-rich worlds to perform a novel statistical comparison of intra-system planetary uniformity across compositionally distinct regimes. We find that, compared to their volatile-rich counterparts, rocky systems are less uniform in mass ($2.6\sigma$), but more uniform in size ($4.0\sigma$) and spacing ($3.0\sigma$). We provide further statistical validation for these results, demonstrating that they are not substantially influenced by the presence of mean motion resonances, low-mass host stars, alternative bulk compositional assumptions, sample size effects, or detection biases. We also obtain tentative evidence ($>2 \sigma$ significance) that the enhanced size uniformity of rocky systems is dominated by the presence of super-Earths ($1 R_{\oplus} \leq R_{p} \leq 1.6 R_{\oplus}$), while their enhanced mass diversity is driven by the presence of sub-Earth ($R_{p} < 1 R_{\oplus}$) worlds.

We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $\tau_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the FRB with a Milky Way-like host galaxy at $z_{\mathrm{host}} = 0.554$ of stellar mass $\mathrm{log}_{10}(M_{\star, \mathrm{host}}) = 10.20^{+0.04}_{-0.03} ~M_\odot$. We identify two intervening galaxy halos at redshifts $z_{\mathrm{igh1}} = 0.492$ and $z_{\mathrm{igh2}} = 0.438$, with low impact parameters, $b_{\mathrm{igh1}} = 43.0_{-11.3}^{+11.3}$ kpc and $b_{\mathrm{igh2}} = 36.1_{-11.3}^{+11.3}$ kpc, and intermediate stellar masses, $\mathrm{log}_{10}(M_{\star, \mathrm{igh1}}) = 10.01^{+0.02}_{-0.02} ~M_\odot$ and $\mathrm{log}_{10}(M_{\star, \mathrm{igh2}}) = 10.60^{+0.02}_{-0.02} ~M_\odot$. The presence of two such galaxies suggests that the sightline is significantly overcrowded compared to the median sightline to this redshift, as inferred from the halo mass function. We perform a detailed analysis of the sightline toward FRB 20221219A, constructing both DM and scattering budgets. Our results suggest that, unlike most well-localized sources, the host galaxy does not dominate the observed scattering. Instead, we posit that an intersection with a single partially ionized cloudlet in the circumgalactic medium of an intervening galaxy could account for the substantial scattering in FRB 20221219A and remain in agreement with typical electron densities inferred for extra-planar dense cloud-like structures in the Galactic and extragalactic halos (e.g., high-velocity clouds).

Image segmentation plays a critical role in unlocking the mysteries of the universe, providing astronomers with a clearer perspective on celestial objects within complex astronomical images and data cubes. Manual segmentation, while traditional, is not only time-consuming but also susceptible to biases introduced by human intervention. As a result, automated segmentation methods have become essential for achieving robust and consistent results in astronomical studies. This review begins by summarizing traditional and classical segmentation methods widely used in astronomical tasks. Despite the significant improvements these methods have brought to segmentation outcomes, they fail to meet astronomers' expectations, requiring additional human correction, further intensifying the labor-intensive nature of the segmentation process. The review then focuses on the transformative impact of machine learning, particularly deep learning, on segmentation tasks in astronomy. It introduces state-of-the-art machine learning approaches, highlighting their applications and the remarkable advancements they bring to segmentation accuracy in both astronomical images and data cubes. As the field of machine learning continues to evolve rapidly, it is anticipated that astronomers will increasingly leverage these sophisticated techniques to enhance segmentation tasks in their research projects. In essence, this review serves as a comprehensive guide to the evolution of segmentation methods in astronomy, emphasizing the transition from classical approaches to cutting-edge machine learning methodologies. We encourage astronomers to embrace these advancements, fostering a more streamlined and accurate segmentation process that aligns with the ever-expanding frontiers of astronomical exploration.

Here, we present the 1UVA catalog of time variable ultraviolet (UV) sources. We describe a new analysis pipeline, the VAriable Source Clustering Analysis (VASCA). We apply the pipeline to 10 years of data from the GALaxy Evolution eXplorer (GALEX) satellite. We analyse a sky area of 302 deg$^2$ , resulting in the detection of 4202 time-variable UV sources. We cross correlate these sources with multi-frequency data from the Gaia satellite and the SIMBAD database, finding an association for 3655 sources. The source sample is dominated by Active Galactic Nuclei ($\approx$73 %) and stars ($\approx$24 %). We look at UV and multi-frequency properties of these sources, focusing on the stellar population. We find UV variability for four White Dwarfs. One of them, WD J004917.14-252556.81, has recently been found to be the most massive pulsating White Dwarf. Its Spectral Energy Distribution shows no sign of a stellar companion. The observed flux variability is unexpected and difficult to explain.

The global magnetic field in the solar corona is known to contain free magnetic energy and magnetic helicity above that of a current-free (potential) state. But the strength of this non-potentiality and its evolution over the solar cycle remain uncertain. Here we model the corona over Solar Cycle 24 using a simplified magneto-frictional model that retains the magnetohydrodynamic induction equation but assumes relaxation towards force-free equilibrium, driven by solar surface motions and flux emergence. The model is relatively conservative compared to some others in the literature, with free energy approximately 20-25% of the potential field energy. We find that unsigned helicity is about a factor 10 higher at Maximum than Minimum, while free magnetic energy shows an even greater increase. The cycle averages of these two quantities are linearly correlated, extending a result found previously for active regions. Also, we propose a practical measure of eruptivity for these simulations, and show that this increases concurrently with the sunspot number, in accordance with observed coronal mass ejection rates. Whilst shearing by surface motions generates 50% or more of the free energy and helicity in the corona, we show that active regions must emerge with their own internal helicity otherwise the eruptivity is substantially reduced and follows the wrong pattern over time.

The spin-down law of pulsars is generally perturbed by two types of timing irregularities: glitches and timing noise. Glitches are sudden changes in the rotational frequency of pulsars, while timing noise is a discernible stochastic wandering in the phase, period, or spin-down rate of a pulsar. We present the timing results of a sample of glitching pulsars observed using the Ooty Radio Telescope (ORT) and the upgraded Giant Metrewave Radio Telescope (uGMRT). Our findings include timing noise analysis for 17 pulsars, with seven being reported for the first time. We detected five glitches in four pulsars and a glitch-like event in J1825-0935. The frequency evolution of glitch in pulsars, J0742-2822 and J1740-3015, is presented for the first time. Additionally, we report timing noise results for three glitching pulsars. The timing noise was analyzed separately in the pre-glitch region and post-glitch regions. We observed an increase in the red noise parameters in the post-glitch regions, where exponential recovery was considered in the noise analysis. Timing noise can introduce ambiguities in the correct evaluation of glitch observations. Hence, it is important to consider timing noise in glitch analysis. We propose an innovative glitch verification approach designed to discern between a glitch and strong timing noise. The novel glitch analysis technique is also demonstrated using the observed data.

In this paper, we compile a \emph{Fermi} sample of the \emph{long} GRBs from \emph{Fermi} observations with 15 years of the Fermi-GBM catalogue with identified redshift, in which the GOLD sample contains 123 long GRBs at $z\le5.6$ and the FULL sample contains 151 long GRBs with redshifts at $z\le8.2$. The Amati relation (the $E_{\rm p,i}$-$E_{\rm iso}$ correlation) are calibrated at $z<1.4$ by a Gaussian Process from the latest observational Hubble data (OHD) with the cosmic chronometers method to obtain GRBs at high-redshift $z\ge1.4$ which can be used to constrain cosmological models via the Markov chain Monte Carlo (MCMC) method. With the cosmology-independent GRBs with the GOLD sample at $z\ge1.4$ and the Pantheon sample of type Ia supernovae (SNe Ia) at $0.01<z\leq2.3$, we obtain $\Omega_{\rm m} = 0.354\pm0.018, H_0 = 73.05\pm0.002$ for the flat $\Lambda$CDM model; $w_0 = -1.22^{+0.18}_{-0.15}$ for the flat $w$CDM model; and $w_{a} = -1.12^{+0.45}_{-0.83}$ for the flat Chevallier-Polarski-Linder model at the 1$\sigma$ confidence level. Our results with the GOLD and FULL sample are almost identical, which are more stringent than the previous results with GRBs.

Planet formation takes place in protoplanetary discs around young T-Tauri stars. PDS 70 is one of the first confirmed examples of a system where the planets are currently forming in gaps in the disc, and can be directly imaged. One of the main early influences on planet formation is the lifetime of the protoplanetary disk, which is limited by the intense stellar X-ray and UV radiation. Stellar coronal activity and accretion of material onto the star are both potential sources of XUV radiation. Previous \textit{Swift} observations detected UV emission, which were consistent with a low rate of accretion. We present follow up observations with the XMM-Newton observatory, which observed PDS 70 simultaneously in X-ray and UV in order to determine intensity of XUV radiation in the system, and identify if the source is coronal, accretion, or both. We detect a strong source in both X-ray and UV, with an average X-ray 0.2-12 keV luminosity of $1.37\times10^{30}\ \mathrm{erg\ s}^{-1}$, and a possible flare which increased the luminosity to $2.8\times10^{30}\ \mathrm{erg\ s}^{-1}$. The UV flux density is in excess of what would be expected from chromospheric emission, and supports the interpretation that PDS 70 has continuing weak accretion less than $\sim10^{-10}\ \mathrm{M_{\odot}\ yr^{-1}}$. The implications of the detected X-ray and UV radiation are that the disc is likely to be in the final stages of dispersal, and will be completely evaporated in the next million years, bringing an end to the primary planet formation process.

Planetary rovers can use onboard data analysis to adapt their measurement plan on the fly, improving the science value of data collected between commands from Earth. This paper describes the implementation of an adaptive sampling algorithm used by PIXL, the X-ray fluorescence spectrometer of the Mars 2020 Perseverance rover. PIXL is deployed using the rover arm to measure X-ray spectra of rocks with a scan density of several thousand points over an area of typically 5 x 7 mm. The adaptive sampling algorithm is programmed to recognize points of interest and to increase the signal-to-noise ratio at those locations by performing longer integrations. Two approaches are used to formulate the sampling rules based on past quantification data: 1) Expressions that isolate particular regions within a ternary compositional diagram, and 2) Machine learning rules that threshold for a high weight percent of particular compounds. The design of the rulesets are outlined and the performance of the algorithm is quantified using measurements from the surface of Mars. To our knowledge, PIXL's adaptive sampling represents the first autonomous decision-making based on real-time compositional analysis by a spacecraft on the surface of another planet.

To address questions about the physical nature and origin of spiral arms in galaxies, it is necessary to measure their dynamical properties, such as the angular speed $\Omega_p$ or the corotation radius. Observations suggest that galaxies may contain several independent spiral patterns simultaneously. It was shown that so-called non-linear resonance coupling plays an important role in such systems. We aim to identify cases of independent spiral patterns for galaxies with a flat rotation curve, and to investigate what relative pattern velocities $\Omega^{out}_{p}/\Omega^{in}_{p}$ they could have in principle for all possible cases of coupling between the main resonances. We solve equations for the main resonance positions (1:1, 2:1, 4:1) and estimate the ratio $\varpi$ of the corotation radii for two subsequent patterns. For six close galaxies with flat rotation curves, we collect the measurements of the corotation radii in the literature, using at least three different methods in each case for credibility, and measure the $\varpi$ ratio. We find $\varpi$ ratios for all possible cases for the main resonances. For three cases we get $\varpi>3$, meaning that it will be difficult to fit two or even more spiral patterns in the disc. These ratios have been used to derive the wind up time for spirals, estimated to be several galactic rotations. We find that three pairs of coupling cases, including vastly acknowledged in galaxies $OLR_{in}=CR_{out} \& CR_{in}=IUHR_{out}$, have very close $\varpi$ ratios, hence should be found simultaneously, as observed. We find strongly confirmed apparent resonance coupling for six galaxies, and show that the observed $\varpi$ is in agreement with theory. In two of them we identify a previously unreported form of simultaneous coupling, namely $OLR_{in}=OUHR_{out} \& OUHR_{in}=CR_{out}$, which was also predicted from the proximity of $\varpi$.

Given the complex nature of galaxies' interstellar medium (ISM), multi-wavelength data are required to probe the interplay among gas, dust, and stellar populations. Spiral galaxies are ideal laboratories for such a goal as they are rich in gas and dust. Using carbon monoxide (CO) along with GALEX far-ultraviolet (FUV) and Spitzer near-infrared (NIR) data we probe the correlations amongst the properties of stellar populations, gas, and dust over the disc of the spiral galaxy NGC~1055 at multiple angular resolutions, i.e. 2, 4, and 17 arcsec corresponding to a linear size of 144 pc, 288 pc, and 1.2 kpc, respectively. Our results indicate an asymmetry in the physical conditions along the galaxy's disc, i.e. the gas is slightly more extended and brighter, and molecular gas mass is higher on the disc's eastern side than the western side. All physical properties (i.e. molecular gas mass, CO line ratios, stellar mass, NIR emission) decrease from the centre going outwards in the disc with some exceptions (i.e. the extinction, FUV radiation, and the [3.6]-[4.5] colour). Our analysis indicates that the colour gets bluer (metallicity increases) halfway through the disc, then redder (metallicity decreases) going outwards further in the disc.

Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport varying the initial conditions of the disk and substructure to find the best match to the observed distributions. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure, namely a perturbation strong enough to be able to trap particles, and that is formed early in the evolution of the disk, that is within 0.4Myr. Agreement is reached by relatively high initial disk masses ($10^{-2.3}M_{\star}\leqslant M_{disk}\leqslant10^{-0.5}M_{\star}$) and moderate levels of turbulence ($10^{-3.5}\leqslant\alpha\leqslant 10^{-2.5}$). Other disk parameters play a weaker role. Only opacities with high absorption efficiency can reproduce the observed spectral indices. Our results extend to the whole population that substructure is likely ubiquitous, so far assessed only in individual disks and implies that most "smooth" disks hide unresolved substructure.

A major prediction of most super-Eddington accretion theories is the presence of anisotropic emission from supercritical disks, but the degree of anisotropy and its dependency with energy remain poorly constrained observationally. A key breakthrough allowing to test such predictions was the discovery of high-excitation photoionized nebulae around Ultraluminous X-ray sources (ULXs). We present efforts to tackle the degree of anisotropy of the UV/EUV emission in super-Eddington accretion flows by studying the emission-line nebula around the archetypical ULX NGC~1313~X--1. We first take advantage of the extensive wealth of optical/near-UV and X-ray data from \textit{Hubble Space Telescope}, \textit{XMM-Newton}, \textit{Swift}-XRT and \textit{NuSTAR} observatories to perform multi-band, state-resolved spectroscopy of the source to constrain the spectral energy distribution (SED) along the line of sight. We then compare spatially-resolved \texttt{Cloudy} predictions using the observed line-of-sight SED with the nebular line ratios to assess whether the nebula `sees' the same SED as observed along the line of sight. We show that to reproduce the line ratios in the surrounding nebula, the photo-ionizing SED must be a factor $\approx 4$ dimmer in ultraviolet emission than along the line-of-sight. Such nearly-iosotropic UV emission may be attributed to the quasi-spherical emission from the wind photosphere. We also discuss the apparent dichotomy in the observational properties of emission-line nebulae around soft and hard ULXs, and suggest only differences in mass-transfer rates can account for the EUV/X-ray spectral differences, as opposed to inclination effects. Finally, our multi-band spectroscopy suggest the optical/near-UV emission is not dominated by the companion star.

The anisotropy of the redshift space bispectrum depends upon the orientation of the triangles formed by three $\vec{k}$ modes with respect to the line of sight. For a triangle of fixed size ($k_1$) and shape ($\mu,t$), this orientation dependence can be quantified in terms of angular multipoles $B_l^m(k_1,\mu,t)$ which contain a wealth of cosmological information. We propose a fast and efficient FFT-based estimator that computes bispectrum multipole moments $B_l^m$ of a 3D cosmological field for all possible $l$ and $m$ (including $m\neq 0$). The time required by the estimator to compute all multipoles from a gridded data cube of volume $N_g^3$ scales as $\sim N_g^3 \log{(N_g)}$ in contrast to the direct computation technique which requires time $\sim N_g^6$. Here, we demonstrate the formalism and validate othe estimator using a simulated non-Gaussian field for which the analytical expressions for all bispectrum multipoles are known. The estimated results are found to be in good agreement with the analytical predictions for all $16$ non-zero multipoles (up to $\ell= 6, m=6$). We expect the $m \neq 0$ bispectrum multipoles to significantly enhance the information available from galaxy redshift surveys, and future redshifted 21-cm observations.

Diagnostics based on the polarization properties of the synchrotron emission can provide precious information on both the ordered structure and the random level of the magnetic field. While this issue has been already analysed in the radio, the polarization data recently obtained by the mission IXPE have shown the need to extend this analysis to the X-rays. While our immediate target are young SNRs, the scope of this analysis is wider. We aim at extending the analysis to particle energy distributions more complex than a power law, as well as to investigate a wider range of cases involving a composition of ordered and random magnetic fields. Since only in a limited number of cases an analytical approach is possible, we have devised to this purpose an optimised numerical scheme, and we have directly used it to investigate particle energy distributions in the form of a power law with an exponential, or super-exponential cutoff. We have also considered a general combination of a ordered field plus an anisotropic random component. We have shown that the previously derived analytic formulae, valid for power-law distributions, may be good approximations of the polarization degree also in the more general case with a cutoff, as typically seen in X-rays. We have explicitly analysed the young SNRs SN 1006, Tycho and Cas A. In particular, for SN 1006, we have proved the consistency between the radio and X-ray polarization degrees, favouring the case of predominantly random field, with an anisotropic distribution. In addition, for the power-law case, we have investigated the effect of a compression on both ordered and random magnetic field components, aimed at describing the mid-age radio SNRs. This work allows a more efficient exploitation of radio and X-ray measurements of the synchrotron polarization, and is addressed to present observations with IXPE as well as to future projects.

There is no consensus yet on whether the precursor and the main burst of gamma-ray bursts (GRBs) have the same origin, and their jet composition is still unclear. In order to further investigate this issue, we systematically search 21 Fermi GRBs with both precursor and main burst for spectral analysis. We first perform Bayesian time-resolved spectral analysis and find that almost all the precursors and the main bursts (94.4$\%$) exhibit thermal components, and the vast majority of them have low-energy spectral index ($\alpha$) (72.2$\%$) that exceed the limit of synchrotron radiation. We then analyse the evolution and correlation of the spectral parameters and find that approximately half of the $\alpha$ (50$\%$) of the precursors and the main bursts evolve in a similar pattern, while peak energy ($E_{p}$) (55.6$\%$) behave similarly, and their evolution is mainly characterized by flux tracking; for the $\alpha-F$ (the flux) relation, more than half of the precursors and the main bursts (61.1$\%$) exhibit roughly similar patterns; the $E_{p}-F$ relation in both the precursor and main burst (100$\%$) exhibits a positive correlation of at least moderate strength. Next, we constrain the outflow properties of the precursors and the main bursts and find that most of them exhibit typical properties of photosphere radiation. Finally, we compare the time-integrated spectra of the precursors and the main bursts and find that nearly all of them are located in similar regions of the Amati relation and follow the Yonetoku relation. Therefore, we conclude that main bursts are continuations of precursors and they may share a common physical origin.

Motivated by the practical interest in the third-body perturbation as a natural cleaning mechanism for high-altitude Earth orbits, we investigate the dynamics stemming from the secular Hamiltonian associated with the lunar perturbation, assuming that the Moon lies on the ecliptic plane. The secular Hamiltonian defined in that way is characterized by two timescales. We compare the location and stability of the fixed points associated with the secular Hamiltonian averaged with respect to the fast variable with the corresponding periodic orbits of the full system. Focusing on the orbit of the Galileo satellites, it turns out that the two dynamics cannot be confused, as the relative difference depends on the ratio between the semi-major axis of Galileo and the one of the Moon, that is not negligible. The result is relevant to construct rigorously the Arnold diffusion mechanism that can drive a natural growth in eccentricity that allows a satellite initially on a circular orbit in Medium Earth Orbit to reenter into the Earth's atmosphere.

While core-collapse supernovae (SNe) often show early and consistent signs of circumstellar (CSM) interaction, some exhibit delayed signatures due to interaction with distant material around the progenitor star. Here we present the discovery in NEOWISE data of WTP19aalnxx, a luminous mid-infrared (IR) transient in the outskirts of the galaxy KUG 0022-007 at $\approx 190$ Mpc. First detected in 2018, WTP19aalnxx reaches a peak absolute (Vega) magnitude of $\approx-22$ at $4.6 \, \mu$m in $\approx3$ yr, comparable to the most luminous interacting SNe. Archival data reveal a $\gtrsim 5\times$ fainter optical counterpart detected since 2015, while follow-up near-IR observations in 2022 reveal an extremely red ($Ks-W2 \approx 3.7$ mag) active transient. Deep optical spectroscopy confirm strong CSM interaction signatures via intermediate-width Balmer emission lines and coronal metal lines. Modeling the broadband spectral energy distribution, we estimate the presence of $\gtrsim 10^{-2}$ M$_\odot$ of warm dust, likely formed in the shock interaction region. Together with the lack of nebular Fe emission, we suggest that WTP19aalnxx is a missed, low (optical) luminosity SN in an emerging family of core-collapse SNe distinguished by their CSM-interaction-powered mid-IR emission that outshines the optical bands. Investigating the Zwicky Transient Facility sample of SNe in NEOWISE data, we find $17$ core-collapse SNe ($\gtrsim 3$% in a volume-limited sample) without early signs of CSM interaction that exhibit delayed IR brightening, suggestive of dense CSM shells at $\lesssim 10^{17}$cm. We suggest that synoptic IR surveys offer a new route to revealing late-time CSM interaction and the prevalence of intense terminal mass loss in massive stars.

We report the discovery of an S-shaped morphology of the radio galaxy J0644$+$1043 imaged with a 30 $\mu$Jy sensitive 525 MHz broadband (band 3 $+$ 4) uGMRT map. Dedicated spectroscopic observations of the host galaxy carried out with the 2-meter Rozhen telescope yielded a redshift of 0.0488, giving a projected linear size of the peculiar radio structure of over 0.7 Mpc. This giant radio galaxy is powered by a black hole of mass 4.1$^{+9.39}_{-2.87}\times 10^8$ \msun, from which vicinity emanate well-collimated and knotty jets, each $\sim$100 kpc long. The entire radio structure, presumably due to the effective jet precession, is less than 50 Myr old, has a power of $\sim$6 $\times 10^{24}$ W Hz$^{-1}$ at 1.4 GHz and the observed morphological characteristics do not strictly conform to the traditional FR I or FR II categories.

This paper highlights methods from geostatistics that are relevant to the interpretation, intercomparison, and synthesis of atmospheric model data, with a specific application to exoplanet atmospheric modeling. Climate models are increasingly used to study theoretical and observational properties of exoplanets, which include a hierarchy of models ranging from fast and idealized models to those that are slower but more comprehensive. Exploring large parameter spaces with computationally-expensive models can be accomplished with sparse sampling techniques, but analyzing such sparse samples can pose challenges for conventional interpolation functions. Ordinary kriging is a statistical method for describing the spatial distribution of a data set in terms of the variogram function, which can be used to interpolate sparse samples across any number of dimensions. Variograms themselves may also be useful diagnostic tools for describing the spatial distribution of model data in exoplanet atmospheric model intercomparison projects. Universal kriging is another method that can synthesize data calculated by models of different complexity, which can be used to combine sparse samples of data from slow models with larger samples of data from fast models. Ordinary and universal kriging can also provide a way to synthesize model predictions with sparse samples of exoplanet observations and may have other applications in exoplanet science.

Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We report the discovery of an Earth-sized planet transiting the nearby (12 pc) inactive M3.0 dwarf Gliese 12 (TOI-6251) with an orbital period ($P_{\rm{orb}}$) of 12.76 days. The planet, Gliese 12b, was initially identified as a candidate with an ambiguous $P_{\rm{orb}}$ from TESS data. We confirmed the transit signal and $P_{\rm{orb}}$ using ground-based photometry with MuSCAT2 and MuSCAT3, and validated the planetary nature of the signal using high-resolution images from Gemini/NIRI and Keck/NIRC2 as well as radial velocity (RV) measurements from the InfraRed Doppler instrument on the Subaru 8.2 m telescope and from CARMENES on the CAHA 3.5 m telescope. X-ray observations with XMM-Newton showed the host star is inactive, with an X-ray-to-bolometric luminosity ratio of $\log L_{\rm X}/L_{\rm bol} \approx -5.7$. Joint analysis of the light curves and RV measurements revealed that Gliese 12b has a radius of 0.96 $\pm$ 0.05 $R_\oplus$, a 3$\sigma$ mass upper limit of 3.9 $M_\oplus$, and an equilibrium temperature of 315 $\pm$ 6 K assuming zero albedo. The transmission spectroscopy metric (TSM) value of Gliese 12b is close to the TSM values of the TRAPPIST-1 planets, adding Gliese 12b to the small list of potentially terrestrial, temperate planets amenable to atmospheric characterization with JWST.

Solar flares, especially C, M, and X class, pose significant risks to satellite operations, communication systems, and power grids. We present a novel approach for predicting extreme solar flares using HMI intensitygrams and magnetograms. By detecting sunspots from intensitygrams and extracting magnetic field patches from magnetograms, we train a Residual Network (ResNet) to classify extreme class flares. Our model demonstrates high accuracy, offering a robust tool for predicting extreme solar flares and improving space weather forecasting. Additionally, we show that HMI magnetograms provide more useful data for deep learning compared to other SDO AIA images by better capturing features critical for predicting flare magnitudes. This study underscores the importance of identifying magnetic fields in solar flare prediction, marking a significant advancement in solar activity prediction with practical implications for mitigating space weather impacts.

Multi-Messenger observations and theory of astrophysical objects is fast becoming a critical research area in the astrophysics scientific community. In particular, point-like objects like that of BL Lac, flat spectrum radio quasars (FSRQ), and blazar candidates of uncertain type (BCU) are of distinct interest among those who look at the synchrotron, Compton, neutrino, and cosmic ray emissions sourced from compact objects. Notably, there is also much interest in the correlation between multi-frequency observations of blazars and neutrino surveys on source demographics. In this review we look at such multi-frequency and multi-physics correlations of the radio, X-ray, and $\gamma$-ray fluxes of different classes of blazars from a collection of survey catalogues. This multi-physics survey of blazars shows that there are characteristic cross-correlations in the spectra of blazars when considering their multi-frequency and multi-messenger emission. Accompanying this will be a review of cosmic ray and neutrino emissions from blazars and their characteristics.

The Transiting Exoplanet Survey Satellite (TESS) has discovered hundreds of new worlds, with TESS planet candidates now outnumbering the total number of confirmed planets from $\textit{Kepler}$. Owing to differences in survey design, TESS continues to provide planets that are better suited for subsequent follow-up studies, including mass measurement through radial velocity (RV) observations, compared to Kepler targets. In this work, we present the TESS-Keck Survey's (TKS) Mass Catalog: a uniform analysis of all TKS RV survey data which has resulted in mass constraints for 126 planets and candidate signals. This includes 58 mass measurements that have reached $\geq5\sigma$ precision. We confirm or validate 32 new planets from the TESS mission either by significant mass measurement (15) or statistical validation (17), and we find no evidence of likely false positives among our entire sample. This work also serves as a data release for all previously unpublished TKS survey data, including 9,204 RV measurements and associated activity indicators over our three year survey. We took the opportunity to assess the performance of our survey, and found that we achieved many of our goals including measuring the mass of 38 small ($<4R_{\oplus}$) planets, nearly achieving the TESS mission's basic science requirement. In addition, we evaluated the performance of the Automated Planet Finder (APF) as survey support and observed meaningful constraints on system parameters due to its more uniform phase coverage. Finally, we compared our measured masses to those predicted by commonly used mass-radius relations and investigated evidence of systematic bias.

In this work, we use 8 years of deep near-infrared imaging to select and study a new set of 601 active galaxies identified through long-term near-infrared (NIR) variability in the UKIDSS Ultra Deep Survey (UDS). These objects are compared to 710 X-ray bright AGN detected by the Chandra X-ray observatory. We show that infrared variability and X-ray emission select distinct sets of active galaxies, finding only a 37 per cent overlap of galaxies detected by both techniques and confirming NIR-variable AGN to be typically X-ray quiet. Examining the mass functions of the active galaxies shows that NIR variability detects AGN activity in galaxies over a significantly wider range of host stellar mass compared to X-ray detection. For example, at z $\sim$ 1, variable AGN are identified among approximately 1 per cent of galaxies in a roughly flat distribution above the stellar mass completeness limit (> 10$^{9}$M$_{\odot}$), while X-ray detection primarily identifies AGN in galaxies of higher mass (> 10$^{10}$M$_{\odot}$). We conclude that long-term near-infrared variability provides an important new tool for obtaining more complete samples of AGN in deep survey fields.

The discovery of GW170817 provided the first direct gravitational-wave measurement of the Hubble constant, $H_0$, demonstrating the potential power of standard-siren cosmology. The dark siren approach can be utilized for gravitational-wave sources in the absence of an electromagnetic counterpart: one considers all galaxies contained within the localization volume as potential hosts. When statistically averaging over the potential host galaxies, one can weight them by their luminosities to account for physically-motivated prescriptions (e.g., tracing star formation or stellar mass). Using mock galaxy catalogs, we explore the impact of these weightings on the measurement of $H_0$, focusing on the bias in $H_0$ inference that results from incorrectly-weighted prescriptions. We find that assuming an incorrect galaxy host probability can lead to significant biases in $H_0$, up to about five times off from typical values inferred by current experiments. These biases are due to inconsistent galaxy weighted redshift distributions as well as preferentially selecting the incorrect host during inference. The magnitudes of these biases are influenced by the galaxy number density along each line of sight, the uncertainty in the measurement of the gravitational-wave luminosity distance, and correlations in the parameter space of galaxies. These biases may be overcome with future GW detectors that contain better GW localization, using a strategic choice of weighting prescription, or with increasing the SNR cut. We propose using hierarchical inference as a diagnosis of incorrectly-weighted prescriptions, which can further be used to simultaneously infer the correct weighting scheme and cosmological parameters.

The low carbon content of Earth and primitive meteorites compared to the Sun and interstellar grains suggests that carbon-rich grains were destroyed in the inner few astronomical units of the young solar system. A promising mechanism to selectively destroy carbonaceous grains is thermal sublimation within the soot line at $\gtrsim$ 300 K. To address whether such hot conditions are common amongst low-mass protostars, we observe CH$_3$CN transitions at 1, 2 and 3 mm with the NOrthern Extended Millimeter Array (NOEMA) toward seven low-mass and one intermediate-mass protostar ($L_{\rm{bol}} \sim2-300 L_\odot$), as CH$_3$CN is an excellent temperature tracer. We find $>$ 300 K gas toward all sources, indicating that hot gas may be prevalent. Moreover, the excitation temperature for CH$_3$OH obtained with the same observations is always lower ($\sim$135-250 K), suggesting that CH$_3$CN and CH$_3$OH have a different spatial distribution. A comparison of the column densities at 1 and 3 mm shows a stronger increase at 3 mm for CH$_3$CN than for CH$_3$OH. Since the dust opacity is lower at longer wavelengths, this indicates that CH$_3$CN is enhanced in the hot gas compared to CH$_3$OH. If this CH$_3$CN enhancement is the result of carbon-grain sublimation, these results suggests that Earth's initial formation conditions may not be rare.

Since neutrinos have mass differences, they could decay into one another. But their lifetimes are likely long, even when shortened by new physics, so decay likely impacts neutrinos only during long trips. This makes high-energy astrophysical neutrinos, traveling for up to billions of light-years, sensitive probes of decay. However, their sensitivity must be tempered by reality. We derive from them thorough bounds on the neutrino lifetimes accounting for critical astrophysical unknowns and the nuances of neutrino detection. Using the diffuse neutrino flux, we disfavor lifetimes $\tau \lesssim 20$-450 s $(m/{\rm eV})$, based on present IceCube data, and forecast factor-of-10 improvements by upcoming detectors. Using, for the first time, neutrinos from the active galaxy NGC 1068, extant unknowns preclude placing lifetime bounds today, but upcoming detectors could disfavor $\tau \sim 100$-5000 s $(m/{\rm eV})$.

In this chapter, we give an overview of our current understanding of the physics of accreting massive black hole binaries (MBHBs), with a special focus on the latest developments in numerical simulations and General-Relativistic Magnetohydrodynamics (GRMHD) simulations in particular. We give a self-contained global picture of how to model accretion onto MBHBs, analyzing different aspects of the system such as the dynamics of the circumbinary disk, mini-disks, outflows, the role of magnetic fields, and electromagnetic signatures. We discuss important questions and open problems related to these systems, what are the advantages and disadvantages of the different numerical approaches, and what robust knowledge we have built from simulations.

A Planck scale inflationary era -- in a quantum gravity theory predicting discreteness of quantum geometry at the fundamental scale -- produces the scale invariant spectrum of inhomogeneities with very small tensor-to-scalar ratio of perturbations and a hot big bang leading to a natural dark matter genesis scenario. Here we evoke the possibility that some of the major puzzles in cosmology would have an explanation rooted in quantum gravity.

With precision pulsar timing, measured values of a large set of pulsar parameters are obtainable. For some of those parameters, such as the time-derivatives of spin or orbital periods (in the case of binary pulsars), the measured values are not the intrinsic values of the parameters as they contain contributions from the dynamical effects. In the case of orbital period derivatives, the intrinsic values are essentially the general relativistic results. Pulsar timing solution also provides measurement of higher time-derivatives of orbital frequency for some pulsars. We specifically focus on the second time-derivative of the orbital frequency to explore its application in testing general relativity. In this work, we have provided a formalism to estimate the general relativistic contribution to the second derivative of the orbital frequency. We have calculated the dynamical effect contributions as well as the general relativistic contributions to the second time-derivative of the orbital period for real as well as synthetic pulsars. We find that the general relativistic contribution to the second time-derivative of the orbital period is negligibly small compared to the observed values of the real pulsars.

Thermal higgsino dark matter (DM), with a mass near 1.1 TeV, is one of the most well-motivated and untested DM candidates. Leveraging recent hydrodynamic cosmological simulations that give DM density profiles in Milky Way analogue galaxies we show that the line-like gamma-ray signal predicted from higgsino annihilation in the Galactic Center could be detected at high significance with the upcoming Cherenkov Telescope Array (CTA) and Southern Wide-field Gamma-ray Observatory (SWGO) for all but the most pessimistic DM profiles. We perform the most sensitive search to-date for the line-like signal using 15 years of data from the Fermi Large Area Telescope, coming within an order one factor of the necessary sensitivity to detect the higgsino for some Milky Way analogue DM density profiles. We show that H.E.S.S. has sub-leading sensitivity relative to Fermi for the higgsino at present. In contrast, we analyze H.E.S.S. inner Galaxy data for the thermal wino model with a mass near 2.8 TeV; we find no evidence for a DM signal and exclude the wino by over a factor of two in cross-section for all DM profiles considered. In the process, we identify and attempt to correct what appears to be an inconsistency in previous H.E.S.S. inner Galaxy analyses for DM annihilation related to the analysis effective area, which may weaken the DM cross-section sensitivity claimed in those works by around an order of magnitude.

Galactic double white dwarf (DWD) binaries are among the guaranteed sources for the Laser Interferometer Space Antenna (LISA), an upcoming space-based gravitational wave (GW) detector. Most DWDs in the LISA band are far from merging and emit quasimonochromatic GWs. As these sources are distributed throughout the Milky Way, they experience different accelerations in the Galactic gravitational potential, and therefore each DWD exhibits an apparent GW frequency chirp due to differential acceleration between the source and LISA. We examine how Galactic acceleration influences parameter estimation for these sources; and investigate how LISA observations could provide insight into the distribution of matter in the Galaxy.

We consider a long-range force, mediated by an ultralight scalar, which can give rise to violation of baryon number. This would lead to very different lifetimes for nucleons in different astrophysical environments. Possible signals of this scenario include a flux of O(10 MeV) solar neutrinos or anomalous heating of old neutron stars; we find the latter to yield the strongest current bounds, which could be improved in the coming years. The ultralight scalar employed here can potentially be a good dark matter candidate.

There has been an attempt to revive the visible QCD axion at the 10 MeV scale assuming that it exclusively couples to the first-generation quarks and the electron. This variant of the QCD axion is claimed to remain phenomenologically viable, partly due to a clever model construction that induces tree-level pion-phobia and exploits uncertainties inherent in the chiral perturbation theory. We confront this model with the cosmological domain wall problem, the quality issue and constraints arising from the electron electric dipole moment. It is also pointed out that the gluon loop-generated axion-top coupling can provide a very large contribution to rare $B$-meson decays, such that the present LHCb data for $B^0 \to K^{*0} e^+ e^-$ rule out the model for the axion mass larger than 30 MeV. There is a strong motivation for pushing the experimental analysis of $B \to K^{(*)} e^+ e^-$ to a lower $e^+ e^-$ invariant mass window, which will conclusively determine the fate of the model, as its contribution to this branching ratio significantly exceeds the Standard Model prediction.

We present a pair of seismometers capable of measurement in all six axes of rigid motion. The vacuum-compatible devices implement compact interferometric displacement sensors to surpass the sensitivity of typical electrical readout schemes. Together with the capability to subtract the sensitivity-limiting coupling of ground tilt into horizontal motion, our seismometers can widen the sensing band towards mHz frequencies. This has notable applications across a range of fields requiring access to low-frequency signals, such as seismology and climate research. We particularly highlight their potential application in gravitational-wave observatories (LIGO) for observation of intermediate-mass black holes ($\sim 1000\,M_\odot$). The sensors are based on a near-monolithic fused-silica design consisting of a fused-silica mass and fibre, showing improved stability and robustness to tilt drifts, alignment, and control compared to all-metal or mixed metal-silica designs. We demonstrate tilt sensitivity that surpasses the best commercial alternatives in a significantly reduced footprint compared to our previous iterations of these sensors.

We examine the shadow cast by a Kerr black hole, focusing on constraints on photons corresponding to different shadow boundaries. The photons are related to different orbital ranges and impact parameter values, creating a map of the shadow boundaries. Our analysis fixes also the conditions under which it is possible to observe an "imprint" of the black hole (outer) ergosurface and (outer) ergoregion on the Kerr black hole shadow boundary. The counter-rotating case resulted strongly constrained with respect to the co-rotating case, constituting a remarkable and significant difference where the counter-rotating component associated with the shadow boundary is strongly distinct from the co-rotating one. However, in this framework, even the co-rotating photons imply restrictions on conditions on the spins and planes, which are bounded by limiting values. We believe the results found here, being a tracer for the central black hole, can constitute new templates for the ongoing observations.

We consider resonant wave-like dark matter conversion into low-frequency radio waves in the Earth's ionosphere. Resonant conversion occurs when the dark matter mass and the plasma frequency coincide, defining a range $m_{ \text{DM} } \sim 10^{-9} - 10^{-8}$ eV where this approach is best suited. Owing to the non-relativistic nature of dark matter and the typical variational scale of the Earth's ionosphere, the standard linearized approach to computing dark matter conversion is not suitable. We therefore solve a second-order boundary-value problem, effectively framing the ionosphere as a driven cavity filled with a positionally-varying plasma. An electrically-small dipole antenna targeting the generated radio waves can be orders of magnitude more sensitive to dark photon and axion-like particle dark matter in the relevant mass range. The present study opens up a promising way of testing hitherto unexplored parameter space which could be further improved with a dedicated instrument.

We investigate the gravitational production of a scalar field $\chi$ with a mass exceeding the Hubble scale during inflation $m_\chi \gtrsim H_I$, employing both analytical and numerical approaches. We demonstrate that the steepest descent method effectively captures the epochs and yields of gravitational production in a compact and simple analytical framework. These analytical results align with the numerical solutions of the field equation. Our study covers three spacetime backgrounds: de Sitter, power-law inflation, and the Starobinsky inflation model. Within these models, we identify two distinct phases of particle production: during and after inflation. During inflation, we derive an accurate analytic expression for the particle production rate, accounting for a varying Hubble rate. After inflation, the additional burst of particle production depends on the inflaton mass around its minimum. When this mass is smaller than the Hubble scale during inflation, $H_I$, there is no significant extra production. However, if the inflaton mass is larger, post-inflation production becomes the dominant contribution. Furthermore, we explore the implications of gravitationally produced heavy fields for dark matter abundance, assuming their cosmological stability.

In Einstein-Aether gravity, we revisit the issue of linear stabilities of black holes against odd-parity perturbations on a static and spherically symmetric background. In this theory, superluminal propagation is allowed and there is a preferred timelike direction along the unit Aether vector field. If we choose the usual spherically symmetric background coordinates with respect to the Killing time $t$ and the areal radius $r$, it may not be appropriate for unambiguously determining the black hole stability because the constant $t$ hypersurfaces are not necessarily always spacelike. Unlike past related works of black hole perturbations, we choose an Aether-orthogonal frame in which the timelike Aether field is orthogonal to spacelike hypersurfaces over the whole background spacetime. In the short wavelength limit, we show that no-ghost conditions as well as radial and angular propagation speeds coincide with those of vector and tensor perturbations on the Minkowski background. Thus, the odd-parity linear stability of black holes for large radial and angular momentum modes is solely determined by constant coefficients of the Aether derivative couplings.

The most pragmatic first step in the all-but-inevitable 3rd-millennium Völkerwanderung of humanity throughout the Solar System is the establishment of a permanent human presence on the Moon. This research examines: 1. the human, agricultural, and technical water needs of a 100-person, 500 m x 100 m x 6 m self-sustaining lunar colony; 2. choosing a strategic location for the moonbase; 3. a heat drill model by which the needed lunar water ice could be sublimated; and 4. the robust water treatment and recovery infrastructure and water management personnel that would be needed for a self-sustaining moonbase.

Space-borne gravitational wave detectors like TianQin might encounter data gaps due to factors like micro-meteoroid collisions or hardware failures. Such glitches will cause discontinuity in the data and have been observed in the LISA Pathfinder. The existence of such data gaps presents challenges to the data analysis for TianQin, especially for massive black hole binary mergers, since its signal-to-noise ratio (SNR) accumulates in a non-linear way, a gap near the merger could lead to significant loss of SNR. It could introduce bias in the estimate of noise properties, and furthermore the results of the parameter estimation. In this work, using simulated TianQin data with injected a massive black hole binary merger, we study the window function method, and for the first time, the inpainting method to cope with the data gap, and an iterative estimate scheme is designed to properly estimate the noise spectrum. We find that both methods can properly estimate noise and signal parameters. The easy-to-implement window function method can already perform well, except that it will sacrifice some SNR due to the adoption of the window. The inpainting method is slower, but it can minimize the impact of the data gap.

The impact of variable material properties, such as temperature-dependent thermal conductivity and dynamical viscosity, on the dynamics of a fully compressible turbulent convection flow beyond the anelastic limit are studied in the present work by two series of three-dimensional direct numerical simulations in a layer of aspect ratio 4 with periodic boundary conditions in both horizontal directions. One simulation series is for a weakly stratified adiabatic background, one for a strongly stratified one. The Rayleigh number is $10^5$ and the Prandtl number is 0.7 throughout this study. The temperature dependence of material parameters is imposed as a power law with an exponent $\beta$. It generates a superadiabaticity $\varepsilon(z)$ that varies across the convection layer. Central statistical quantities of the flow, such as the mean superadiabatic temperature, temperature and density fluctuations, or turbulent Mach numbers are compared in the form of horizontal plane-time averaged profiles. It is found that the additional material parameter dependence causes systematic quantitative changes of all these quantities, but no qualitative ones. A growing temperature power law exponent $\beta$ also enhances the turbulent momentum transfer in the weak stratification case by 40\%, it reduces the turbulent heat transfer by up to 50\% in the strong stratification case.

The last few decades have provided abundant evidence for physics beyond the two standard models of particle physics and cosmology. As is now known, the by far largest part of our universe's matter/energy content lies in the `dark' and consists of dark energy and dark matter. Despite intensive efforts on the experimental as well as the theoretical side, the origins of both are still completely unknown. Screened scalar fields have been hypothesized as potential candidates for dark energy or dark matter. Among these, some of the most prominent models are the chameleon, symmetron, and environment-dependent dilaton. In this article, we present a summary containing the most recent experimental constraints on the parameters of these three models. For this, experimental results have been employed from the qBOUNCE collaboration, neutron interferometry, and Lunar Laser Ranging (LLR), among others. In addition, constraints are forecast for the Casimir And Non Newtonian force EXperiment (CANNEX). Combining these results with previous ones, this article collects the most up-to-date constraints on the three considered screened scalar field models.

Effective and ethical mentorship practices are crucial to improving recruitment and retention especially for historically minoritized groups (HMGs). Spectrum is a diversity, inclusion, equity, and accessibility (DEIA) grassroots organization committed to empowering equitable excellence through sustainable change. By improving transparency and DEIA within the fields of physics and astronomy, we can empower the next generation of diverse scientists and increase field retention. Starting within our home department at George Mason University and moving outwards, we ensure our students leave as advocates for DEIA and AJEDI (access, justice, equity, diversity, and inclusion) through education and mentorship. Spectrum is providing professionally trained peer mentors to aid students in all facets of their academic and personal lives. Although the peer mentoring program existed since the creation of Spectrum in Spring 2020, we have recently developed and implemented a formal mentorship training for both student and faculty mentors thus increasing the quality, trustworthiness, and confidence of our mentors. Using the latest mentorship research available, this training is developed by Spectrum for George Mason University, with the ability to implement the training at any institution.

The detection of light of thermal origin is the principal means by which humanity has learned about our world and the cosmos. In optical astronomy, in particular, direct detection of thermal photons and the resolution of their spectra have enabled discoveries of the broadest scope and impact. Such measurements, however, do not capture the phase of the thermal fields--a parameter that has proven crucial to transformative techniques in radio astronomy such as synthetic aperture imaging. Over the last 25 years, tremendous progress has occurred in laser science, notably in the phase-sensitive, broad bandwidth, high resolution, and traceable spectroscopy enabled by the optical frequency comb. In this work, we directly connect the fields of frequency comb laser spectroscopy and passive optical sensing as applied to astronomy, remote sensing, and atmospheric science. We provide fundamental sensitivity analysis of dual-comb correlation spectroscopy (DCCS), whereby broadband thermal light is measured via interferometry with two optical frequency combs. We define and experimentally verify the sensitivity scaling of DCCS at black body temperatures relevant for astrophysical observations. Moreover, we provide comparison with direct detection techniques and more conventional laser heterodyne radiometry. Our work provides the foundation for future exploration of comb-based broadband synthetic aperture hyperspectral imaging across the infrared and optical spectrum.