Papers for Thursday, Nov 07 2024

Modeling strong gravitational lenses is computationally expensive for the complex data from modern and next-generation cosmic surveys. Deep learning has emerged as a promising approach for finding lenses and predicting lensing parameters, such as the Einstein radius. Mean-variance Estimators (MVEs) are a common approach for obtaining aleatoric (data) uncertainties from a neural network prediction. However, neural networks have not been demonstrated to perform well on out-of-domain target data successfully - e.g., when trained on simulated data and applied to real, observational data. In this work, we perform the first study of the efficacy of MVEs in combination with unsupervised domain adaptation (UDA) on strong lensing data. The source domain data is noiseless, and the target domain data has noise mimicking modern cosmology surveys. We find that adding UDA to MVE increases the accuracy on the target data by a factor of about two over an MVE model without UDA. Including UDA also permits much more well-calibrated aleatoric uncertainty predictions. Advancements in this approach may enable future applications of MVE models to real observational data.

The Rayleigh Taylor instability (RTI) occurs at the interface between two fluids of different densities, notably when a heavier fluid sits above a lighter one in an effective gravitational field. This instability is relevant to many astrophysical systems where relativistic effects are significant. We examine the linear theory of relativistic Rayleigh Taylor instability (RRTI) in a magnetized medium, allowing for relativistic fluid motion parallel to the interface. To simplify our derivations, we use an "intermediate frame" where fluids on both sides have the same Lorentz factor. Our analysis yields the dispersion relation for RRTI. We find that the instability occurs when the Atwood number $\mathcal{A}$ = $(\rho_1 h_1 - \rho_2 h_2) / (\rho_1 h_1 + \rho_2 h_2) >0$, without requiring relativistic correction. Relativistic motion increases the effective inertia ($\rho \rightarrow \gamma_*^2 \rho$), weakening the magnetic field's suppression of the instability. In the laboratory frame, the instability growth rate is reduced due to time dilation. These analytical results may inform studies of instabilities in systems such as microquasars, active galactic nuclei, gamma-ray bursts, and radio pulsars.

A revised catalogue of 310 Galactic supernova remnants (SNRs) is presented, along with some statistics of their properties. 21 SNRs have been added to the catalogue since the previous published version from 2019, and 5 entries have been removed, as they have been identified as HII regions. Also discussed are some basics statistics of the remnants in the catalogue, the selection effects that apply to the identification of Galactic SNRs and their consequences.

The detection of satellite thermal emission at millimeter wavelengths is presented using data from the 3rd-Generation receiver on the South Pole Telescope (SPT-3G). This represents the first reported detection of thermal emission from artificial satellites at millimeter wavelengths. Satellite thermal emission is shown to be detectable at high signal-to-noise on timescales as short as a few tens of milliseconds. An algorithm for downloading orbital information and tracking known satellites given observer constraints and time-ordered observatory pointing is described. Consequences for cosmological surveys and short-duration transient searches are discussed, revealing that the integrated thermal emission from all large satellites does not contribute significantly to the SPT-3G survey intensity map. Measured satellite positions are found to be discrepant from their two-line element (TLE) derived ephemerides up to several arcminutes which may present a difficulty in cross-checking or masking satellites from short-duration transient searches.

We present MARFA (Molecular atmospheric Absorption with Rapid and Flexible Analysis) -- a streamlined efficient tool for line-by-line calculation of atmospheric absorption signatures in the form of PT lookup tables, which may be used in radiative transfer codes. The tool is intended for computations in the IR and visible spectral regions. Core line-by-line scheme features nine-grid interpolation technique, which strikes good balance between speed and accuracy for calculations in the scenario of unknown continuum function and large lines cut-offs. The model features high flexibility, allowing fast recalculations with different line shapes, $\chi$-factors, line cut-offs conditions, and other parameters, making it valuable for planetary studies where atmospheric and spectroscopic data are sparse or uncertain. The MARFA tool is provided in two ways: through a web interface for onboarding and immediate usage, and as open-source code available in a public repository for advanced utilization, development and contributions.

The Daisy World model has long served as a foundational framework for understanding the self-regulation of planetary biospheres, providing insights into the feedback mechanisms that may govern inhabited exoplanets. In this study, we extend the classic Daisy World model through the lens of Semantic Information Theory (SIT), aiming to characterize the information flow between the biosphere and planetary environment -- what we term the \emph{information architecture} of Daisy World systems. Our objective is to develop novel methodologies for analyzing the evolution of coupled planetary systems, including biospheres and geospheres, with implications for astrobiological observations and the identification of agnostic biosignatures. To operationalize SIT in this context, we introduce a version of the Daisy World model tailored to reflect potential conditions on M-dwarf exoplanets, formulating a system of stochastic differential equations that describe the co-evolution of the daisies and their planetary environment. Analysis of this Exo-Daisy World model reveals how correlations between the biosphere and environment intensify with rising stellar luminosity, and how these correlations correspond to distinct phases of information exchange between the coupled systems. This \emph{rein control} provides a quantitative description of the informational feedback between the biosphere and its host planet. Finally, we discuss the broader implications of our approach for developing detailed ExoGaia models of inhabited exoplanetary systems, proposing new avenues for interpreting astrobiological data and exploring biosignature candidates.

The time until black hole formation in a binary neutron-star (NS) merger contains invaluable information about the nuclear equation of state (EoS) but has thus far been difficult to measure. We propose a new way to constrain the merger remnant's NS lifetime, which is based on the tendency of the NS remnant neutrino-driven winds to enrich the ejected material with helium. Based on the He I $\lambda 1083.3$ nm line, we show that the feature around 800-1200 nm in AT2017gfo at 4.4 days seems inconsistent with a helium mass fraction of $X_{\mathrm{He}} \gtrsim 0.05$ in the polar ejecta. Recent neutrino-hydrodynamic simulations of merger remnants are only compatible with this limit if the NS remnant collapses within 20-30 ms. Such a short lifetime implies that the total binary mass of GW170817, $M_\mathrm{\rm tot}$, lay close to the threshold binary mass for direct gravitational collapse, $M_\mathrm{thres}$, for which we estimate $M_{\mathrm{thres}}\lesssim 2.93 M_\odot$. This upper bound on $M_\mathrm{thres}$ yields upper limits on the radii and maximum mass of cold, non-rotating NSs, which rule out simultaneously large values for both quantities. In combination with causality arguments, this result implies a maximum NS mass of $M_\mathrm{max}\lesssim2.3 M_\odot$. The combination of all limits constrains the radii of 1.6 M$_\odot$ NSs to about 12$\pm$1 km for $M_\mathrm{max}$ = 2.0 M$_\odot$ and 11.5$\pm$1 km for $M_\mathrm{max}$ = 2.15 M$_\odot$. This $\sim2$ km allowable range then tightens significantly for $M_\mathrm{max}$ above $\approx2.15$ M$_\odot$. This rules out a significant number of current EoS models. The short NS lifetime also implies that a black-hole torus, not a highly magnetized NS, was the central engine powering the relativistic jet of GRB170817A. Our work motivates future developments... [abridged]

The vertical evolution of galactic discs is governed by the sub-structures within them. We examine the diversity of kinematic sub-structure present in the first 12 galaxies observed from the GECKOS survey, a VLT/MUSE large programme providing a systematic study of 35 edge-on, Milky Way-mass disc galaxies. Employing the nGIST analysis pipeline, we derive the mean line-of-sight stellar velocity ($V_{\star}$), velocity dispersion ($\sigma_{\star}$), skew ($h_{3}$), and kurtosis ($h_{4}$) for the sample, and examine 2D maps and 1D line profiles. Visually, the majority of this sample (8/12) are found to possess boxy-peanut bulges and host the corresponding kinematic structure predicted for stellar bars viewed in projection. Four galaxies exhibit strong evidence for the presence of nuclear discs, including central $h_{3}$-$V_{\star}$ anti-correlations, `croissant'-shaped central depressions in $\sigma_{\star}$ maps, strong gradients in $h_{3}$, and positive $h_{4}$ plateaus over the expected nuclear disc extent. The strength of the $h_{3}$ feature corresponds to the size of the nuclear disc, measured from the $h_{3}$ turnover radius. We can explain the features within the kinematic maps of all sample galaxies via disc structure(s) alone. We do not find any need to invoke the existence of dispersion-dominated bulges. Obtaining the specialised data products for this paper and the broader GECKOS survey required significant development of existing integral field spectroscopic (IFS) analysis tools. Therefore, we also present the nGIST pipeline: a modern, sophisticated, and easy-to-use pipeline for the analysis of galaxy IFS data. We conclude that the variety of kinematic sub-structures seen in GECKOS galaxies requires a contemporary view of galaxy morphology, expanding on the traditional view of galaxy structure, and uniting the kinematic complexity observed in the Milky Way with the extragalactic.

We present $24$ cosmological dark matter (DM)--only zoom-in simulations of a Milky Way (MW) analog with initial conditions appropriate for scenarios where non-cold DM is a fraction of the total DM abundance (f-NCDM models), as the second installment of the COZMIC suite. We initialize our simulations using transfer functions, $T_{\mathrm{f-NCDM}}(k)\equiv\sqrt{P_{\mathrm{f-NCDM}}(k)/P_{\mathrm{CDM}}(k)}$ (where $P(k)$ is the linear matter power spectrum), with an initial suppression similar to thermal-relic warm DM (WDM) followed by a constant-amplitude plateau. We simulate suppression wave numbers $[22.8,~ 32.1,~ 41.8,~ 52.0,~ 57.1,~ 95.3]~\mathrm{Mpc}^{-1}$, corresponding to thermal-relic WDM masses $m_{\mathrm{WDM}}\in [3,~ 4,~ 5,~ 6,~ 6.5,~ 10]~\mathrm{keV}$, and plateau amplitudes $\delta\in [0.2,~ 0.4,~ 0.6,~ 0.8]$. We model the subhalo mass function in terms of the suppression wave number and $\delta$. Integrating these models into a forward model of the MW satellite galaxy population yields new limits on f-NCDM scenarios, with suppression wave numbers greater than $46$ and $ 40~\mathrm{Mpc}^{-1}$ for $\delta=0.2$, $0.4$, respectively, at $95\%$ confidence. The current data do not constrain $\delta>0.4$. We map these limits to scenarios where a fraction $f_{\mathrm{WDM}}$ of DM behaves as a thermal relic, which yields the following bounds on cosmologies with a mixture of WDM and CDM: $m_{\mathrm{WDM}}>3.6,~ 4.1,~ 4.6,~ 4.9,~ 5.4~\mathrm{keV}$ for $f_{\mathrm{WDM}}=0.5,~ 0.6,~ 0.7,~ 0.8,~ 0.9$, respectively, at $95\%$ confidence. The current data do not constrain WDM fractions $f_{\mathrm{WDM}}<0.5$. Our results affirm that low-mass halo abundances are sensitive to partial suppression in $P(k)$, indicating the possibility of using galactic substructure to reconstruct $P(k)$ on small scales.

The advancement of the Event Horizon Telescope has enabled the study of relativistic jets in active galactic nuclei down to sub-parsec linear scales even at high redshift. Quasi-simultaneous multifrequency observations provide insights into the physical conditions in compact regions and allow testing accretion theories. Initially we aimed at measuring the magnetic field strength close to the central supermassive black hole in NRAO 530 (1730-130) by studying frequency-dependent opacity of the jet matter, Faraday rotation and the spectral index in the mm-radio bands. NRAO 530 was observed quasi-simultaneously at 15, 22, 43, 86, and 227 GHz at four different very long baseline interferometer (VLBI) networks. By the means of imaging and model-fitting, we aligned the images, taken at different frequencies. We explored opacity along the jet and distribution of the linearly polarized emission in it. Our findings reveal that the jet of NRAO 530 at 86 and 227 GHz is transparent down to its origin, with 70 mJy emission detected at 227 GHz potentially originating from the accretion disk. The magnetic field strength near the black hole, estimated at $5r_\mathrm{g}$, is $3\times10^3-3\times10^4$ G (depending on the central black hole mass). These values represent some of the highest magnetic field strengths reported for active galaxies. We also report the first ever VLBI measurement of the Faraday rotation at 43-227 GHz, which reveals rotation measure values as high as -48000 rad/m2 consistent with higher particle density and stronger magnetic fields at the jet's outset. The complex shape of the jet in NRAO 530 is in line with the expected behavior of a precessing jet, with a period estimated to be around $6\pm4$~years.

The presence of supermassive black holes (SMBHs, $M_{\bullet}\sim 10^{6-10}~M_{\odot}$) in the first cosmic Gyr ($z\gtrsim 6$) challenges current models of BH formation and evolution. We propose a novel mechanism for the formation of early SMBH seeds based on primordial black holes (PBHs). We assume a non-Gaussian primordial power spectrum as expected in inflationary models; these scenarios predict that PBHs are initially clustered and preferentially formed in the high-$\sigma$ fluctuations of the large-scale density field, out of which dark matter (DM) halos are originated. Our model accounts for (i) PBH accretion and feedback, (ii) DM halo growth, and (iii) gas dynamical friction. PBHs lose angular momentum due to gas dynamical friction, sink into a dense core, where BH binaries form and undergo a runaway merger, eventually leading to the formation of a single, massive seed. This mechanism starts at $z\sim 20-40$ in rare halos ($M_h\sim 10^7\ M_\odot$ corresponding to $\sim 5-7\sigma$ fluctuations), and provides massive ($\sim 10^{4-5}~ M_{\odot}$) seeds by $z\sim 10-30$. We derive a physically-motivated seeding prescription that provides the mass of the seed, $ M_{\rm seed}(z)=3.1\times 10^{5}\ { M_{\odot}}[(1+z)/10]^{-1.2}$, and seeded halo, $ M_{h}(z)=2\times 10^{9}\ {M_{\odot}}[(1+z)/10]^{-2}e^{-0.05z}$ as a function of redshift. This seeding mechanism requires that only a small fraction of DM is constituted by PBHs, namely $f_{\rm PBH}\sim 3 \times 10^{-6}$. We find that $z\sim 6-7$ quasars can be explained with $6\times 10^4 M_{\rm \odot}$ seeds planted at $z\sim 32$, and growing at sub-Eddington rates, $\langle\lambda_{\rm E}\rangle\sim 0.55$. The same scenario reproduces the BH mass of GNz11 at $z=10.6$, while UHZ1 ($z=10.1$) and GHZ9 ($z=10$) data favour instead slightly later ($z\sim 20-25$), more massive ($10^5~M_{\rm \odot}$) seeds. [Abridged]

Our understanding of the process of terrestrial planet formation has grown markedly over the past 20 years, yet key questions remain. This review begins by first addressing the critical, earliest stage of dust coagulation and concentration. While classic studies revealed how objects that grow to $\sim$meter sizes are rapidly removed from protoplanetary disks via orbital decay (seemingly precluding growth to larger sizes), this chapter addresses how this is resolved in contemporary, streaming instability models that favor rapid planetesimal formation via gravitational collapse of solids in over-dense regions. Once formed, planetesimals grow into Mars-Earth-sized planetary embryos by a combination of pebble- and planetesimal accretion within the lifetime of the nebular disk. After the disk dissipates, these embryos typically experience a series of late giant impacts en route to attaining their final architectures. This review also highlights three different inner Solar System formation models that can match a number of empirical constraints, and also reviews ways that one or more might be ruled out in favor of another in the near future. These include (1) the Grand Tack, (2) the Early Instability and (3) Planet Formation from Rings. Additionally, this chapter discusses formation models for the closest known analogs to our own terrestrial planets: super-Earths and terrestrial exoplanets in systems also hosting gas giants. Finally, this review lays out a chain of events that may explain why the Solar System looks different than more than 99% of exoplanet systems.

Planets form from the same cloud of molecular gas and dust as their host stars. Confirming if planetary bodies acquire the same refractory element composition as their natal disc during formation, and how efficiently volatile elements are incorporated into growing planets, is key to linking the poorly constrained interior composition of rocky exoplanets to the observationally-constrained composition of their host star. Such comparisons also afford insight into the planet formation process. This work compares planetary composition with host-star composition using observations of a white dwarf that has accreted planetary material and its F-type star wide binary companion as a reference for the composition of the natal molecular gas and dust. Spectroscopic analysis reveals abundances of Fe, Mg, Si, Ca, and Ti in both stars. We use the white dwarf measurements to estimate the composition of the exoplanetary material and the F-type companion to constrain the composition of the material the planet formed from. Comparing planetary material to the composition of its natal cloud, our results reveal that the planetary material is depleted in moderate refractories (Mg, Si, Fe) relative to the refractory material (Ca, Ti). Grouping elements based on their condensation temperatures is key to linking stellar and planetary compositions. Fractionation during formation or subsequent planetary evolution leads to the depletion of moderate refractories from the planetary material accreted by the white dwarf. This signature, as seen for bulk Earth, will likely be present in the composition of many exoplanets relative to their host-stars.

Detecting exoplanets and other faint sources of emitted and reflected light near a bright star requires deeply suppressing the starlight while efficiently transmitting the dim light from its surroundings. This suppression can be carried out by coronagraphs, nulling interferometers, and starshades. This chapter provides a brief overview of these technologies, emphasizing coronagraphs.

Several missions following Planck are currently under development, which will provide high-precision measurements of the Cosmic Microwave Background (CMB) anisotropies. Specifically, measurements of the E modes will become nearly limited by cosmic variance, which, especially when considering the sharpness of the E-mode transfer functions, may allow for the ability to detect deviations from the concordance model in the CMB data. We investigate the capability of upcoming missions to scrutinize models that have been proposed to address large-scale anomalies observed in the temperature spectra from WMAP and Planck. To this purpose, we consider four benchmarks that modify the CMB angular power spectra at large scales: models producing suppression, a dip, and amplification in the primordial scalar power spectrum, as well as a beyond-Lambda CDM prescription of dark energy. Our analysis shows that large-scale measurements from LiteBIRD will be able to distinguish between various types of primordial and late-time models that predict modifications to the angular spectra at these scales. Moreover, if these deviations from the standard cosmological model are determined to be systematic and do not reflect the true universe model, future experiments could potentially dismiss these features as statistical fluctuations. We also show that additional measurements from CMB-S4 can impose more stringent constraints by probing correlated signals that these models predict at smaller scales (l>100). A byproduct of our analysis is that a recently proposed "Dark Dimension" scenario, featuring power amplification at large scales, is strongly bound by current data, pushing the deviation from the standard model to unobservable scales. Overall, our results demonstrate that future CMB measurements can provide valuable insights into large-scale anomalies that are present in the current CMB data.

A high resolution fourth-order Padé scheme is used to simulate locally isothermal 3D disk turbulence driven by the vertical shear instability (VSI) using 268.4 M points. In the early non-linear period of axisymmetric VSI, angular momentum transport by vertical jets creates correlated N-shaped radial profiles of perturbation vertical and azimuthal velocity. This implies dominance of positive perturbation vertical vorticity layers and a recently discovered angular momentum staircase with respect to radius ($r$). These features are present in 3D in a weaker form. The 3D flow consists of vertically and azimuthally coherent turbulent shear layers with small vortices with all three components are active. Previously observed large persistent vortices in the interior of the domain driven by the Rossby wave instability are absent. We speculate that this is due to a weaker angular momentum staircase in 3D in the present simulations than in a previous simulation. The turbulent viscosity parameter $\alpha(r)$ increases linearly with $r$. At intermediate resolution, the value of $\alpha(r)$ at mid-radius is close to that of a previous simulation. The specific kinetic energy spectrum with respect to radial wavenumber has a power law region with exponent $-1.84$, close to the value $-2$ expected for shear layers. The spectrum with respect to azimuthal wavenumber has a $-5/3$ region and lacks a $-5$ region reported in an earlier study. Finally, it is found axisymmetric VSI has artifacts at late times, including a very strong angular momentum staircase, which in 3D is present weakly in the upper layers of the disk.

JWST is opening many avenues for exploration. For cold brown dwarfs and exoplanets, JWST has opened the door to the mid-infrared wavelength region, where such objects emit significant energy. For the first time, astronomers have access to mid-infrared spectroscopy for objects colder than 600 K. The first spectra appear to validate the model suite known as ATMO 2020++: atmospheres which include disequilibrium chemistry and have a non-adiabatic pressure-temperature relationship. Preliminary fits to JWST spectroscopy of Y dwarfs show that the slope of the energy distribution from lambda = 4.5 um to lambda = 10 um is very sensitive to gravity. We explore this phenomenon using PH3-free ATMO 2020++ models and updated WISE W2 - W3 colors. We find that an absolute 4.5 um flux measurement constrains temperature, and the ratio of the 4.5 um flux to the 10 - 15 um flux is sensitive to gravity and less sensitive to metallicity. We identify 10 T dwarfs with red W2 - W3 colors which are likely to be very low gravity, young, few-Jupiter-mass objects; one of these is the previously known COCONUTS-2b. The unusual Y dwarf WISEPA J182831.08+265037.8 is blue in W2 - W3 and we find that the 4 to 18 um JWST spectrum is well reproduced if the system is a pair of high gravity 400 K dwarfs. Recently published JWST colors and luminosity-based effective temperatures for late-T and Y dwarfs further corroborate the ATMO 2020++ models, demonstrating the potential for significant improvement in our understanding of cold very low-mass bodies in the solar neighborhood.

We use the GALPROP cosmic ray (CR) framework to model the Galactic CR distributions and associated non-thermal diffuse emissions up to PeV energies. We consider ensembles of discrete, finite lifetime CR sources, e.g.\ supernova remnants (SNRs), for a range of creation rates and lifetimes. We find that global properties of the CR sources are likely not directly recoverable from the current `snapshot' of the historic injection and propagation of CRs within the Galaxy that are provided by the data. We show that models for the diffuse $\gamma$ rays based on the discrete/time-dependent scenarios we consider are able to explain LHAASO very-/ultra-high energy (VHE/UHE) $\gamma$-ray data with up to 50\% contribution by unresolved leptonic sources at the highest energies. Over the models that we consider, variations in the diffuse VHE emissions can be $\sim$25\%, which is comparable to those for steady-state models that we investigated in earlier work. Such variations due to the discrete/finite nature of the CR sources are an important factor that are necessary to construct accurate physical models of the diffuse emissions from the Galaxy at VHE/UHEs.

The massive semi-detached binary with a long-term decreasing orbital period may involve a rapid mass-transfer phase in Case A, and thus they are good astrophysical laboratories for investigating the evolution of massive binary stars. In this work, by using the long-term observational light curves from the OGLE project and other data in the low-metallicity LMC, four semi-detached massive binaries with long-term decreases in the orbital periods are detected from 165 EB-type close binaries. It is found that the more massive component in S07798 is filling its Roche lobe where the period decrease is caused by mass transfer from the primary to the secondary. However, the other three (S03065, S12631, S16873) are semi-detached binaries with a lobe-filling secondary where the mass transfer between the components should cause the period to increase if the angular momentum is conservative. The long-term period decreases in these three systems may be caused by the angular momentum loss. Additionally, the orbital periods of three systems (S03065, S07798, S16873) are detected to show cyclic variation with periods shorter than 11 years, which can be plausibly explained by the presence of close-in third bodies in these massive binaries. Based on all of these results, it is suggested that the detected four semi-detached binaries almost have multiplicity. The companion stars are crucial for the origin and evolution of these massive close binaries.

We present ALMA observations on the high-redshift galaxy GHZ2 and report a successful detection of the rest-frame 88um atomic transition from doubly-ionized Oxygen at z=12.3327+/-0.0005. Based on these observations, combined with additional constraints on the [OIII] 52um line luminosity and previous JWST data, we argue that GHZ2 is likely powered by compact and young star formation, and show that it follows well-established relationships found for giant HII regions and metal-poor star-forming dwarf galaxies that are known to host bright super star clusters. Additionally, these observations provide new constraints on the Oxygen electron density (100 < n_e[cm^-3] < 4,000) and dynamical mass (M_dyn=3-8x10^8M_sun). The existence of these massive starburst systems 13.3Gyr ago might explain the origin of today's globular clusters, a long-standing question in astronomy. To test this, we present observational probes to investigate whether sources like GHZ2 are linked to the formation of today's globular clusters or other more massive compact stellar systems.

Cosmology is entering a very exciting time in its history, when a wealth of cutting-edge experiments are all starting to collect data, or about to. These experiments aim at addressing some of the most intriguing questions in fundamental physics, such as what is the nature of dark matter, is dark energy a cosmological constant or a varying field, what are the masses of the neutrinos, and more. While Lambda-CDM has emerged as a simple model that is consistent with most of the current data sets, we are starting to see some interesting deviations that deserve further exploration. This contribution provides an overview of upcoming projects and the science opportunities they will allow. In particular, we recall and comment the DESI year-1 BAO constraints and their implications for dark energy. We put some of the most recent results and outstanding questions in the perspective of the forthcoming observational program.

We performed shot analyses of X-ray and optical sub-second flares observed during the low/hard state of the 2018 outburst in MAXI J1820$+$070. Optical shots were less spread than X-ray shots. The amplitude of X-ray shots was the highest at the onset of the outburst, and they faded at the transition to the intermediate state. The timescale of shots was $\sim$0.2 s, and we detected the abrupt spectral hardening synchronized with this steep flaring event. The time evolution of optical shots was not similar to that of X-ray shots. These results suggest that accreting gas blobs triggered a series of magnetic reconnections at the hot inner accretion flow in the vicinity of the black hole, which enhanced X-ray emission and generated flaring events. The rapid X-ray spectral hardening would be caused by this kind of magnetic activity. Also, the synchrotron emission not only at the hot flow but also at the jet plasma would contribute to the optical rapid variability. We also found that the low/hard state exhibited six different phases in the hardness-intensity diagram and the correlation plot between the optical flux and the X-ray hardness. The amplitude and duration of X-ray shots varied in synchrony with these phases. This time variation may provide key information about the evolution of the hot flow, the low-temperature outer disk, and the jet-emitting plasma.

Most earlier studies have been limited to estimating the kinematic evolution of coronal mass ejections (CMEs), and only limited efforts have been made to investigate their thermodynamic evolution. We focus on the interplay of the thermal properties of CMEs with their observed global kinematics. We implement the Flux rope Internal State (FRIS) model to estimate variations in the polytropic index, heating rate per unit mass, temperature, pressure, and various internal forces. The model incorporates inputs of 3D kinematics obtained from the Graduated Cylindrical Shell (GCS) model. In our study, we chose nine fast-speed CMEs from 2010 to 2012. Our investigation elucidates that the selected fast-speed CMEs show a heat-release phase at the beginning, followed by a heat-absorption phase with a near-isothermal state in their later propagation phase. The thermal state transition, from heat release to heat absorption, occurs at around 3($\pm$0.3) to 7($\pm$0.7) $R_\odot$ for different CMEs. We found that the CMEs with higher expansion speeds experience a less pronounced sharp temperature decrease before gaining a near-isothermal state. The differential emission measurement (DEM) analysis findings, using multi-wavelength observation from SDO/AIA, also show a heat release state of CMEs at lower coronal heights. We also find the dominant internal forces influencing CME radial expansion at varying distances from the Sun. Our study shows the need to characterize the internal thermodynamic properties of CMEs better in both observational and modeling studies, offering insights for refining assumptions of a constant value of the polytropic index during the evolution of CMEs.

Relativistic jets manifest some of the most intriguing activities in the nuclear regions of active galaxies. Identifying the most powerful relativistic jets permits us to probe the most luminous accretion systems and, in turn, the most massive black holes. This paper reports the identification of one such object, PMN J1310$-$5552 ($z=1.56$), a blazar candidate of uncertain type detected with the Fermi Large Area Telescope (LAT) and Swift Burst Alert Telescope. The detection of broad emission lines in its optical spectra taken with the X-Shooter and Goodman spectrographs classifies it to be a flat-spectrum radio quasar. The analysis of the Goodman optical spectrum has revealed PMN J1310$-$5552 harbors a massive black hole (log scale $M_{\rm BH}=9.90\pm0.07$, in $M_{\odot}$) and luminous accretion disk (log scale $L_{\rm disk}=46.86\pm0.03$, in erg s$^{-1}$). The fitting of the observed big blue bump with the standard accretion disk model resulted in the log scale $M_{\rm BH}=9.81^{+0.19}_{-0.20}$ (in $M_{\odot}$) and $L_{\rm disk}=46.86^{+0.09}_{-0.09}$ (in erg s$^{-1}$), respectively. These parameters suggest PMN J1310$-$5552 hosts one of the most massive black holes and the most luminous accretion disks among the blazar population. The physical properties of this enigmatic blazar were studied by modeling the broadband spectral energy distribution considering the data from NuSTAR, Swift, Fermi-LAT, and archival observations. Overall, PMN J1310$-$5552 is a powerful `MeV' blazar with physical parameters similar to other members of this unique class of blazars. These results provide glimpses of monsters lurking among the unknown high-energy emitters and demonstrate the importance of ongoing wide-field sky surveys to discover them.

The NEO Coordination Centre (NEOCC) of the European Space Agency is an operational centre that, among other activities, computes the orbits of near-Earth objects and their probabilities of impact with the Earth. The NEOCC started providing information about near-Earth objects in 2012 on a dedicated web portal, accessible at this https URL. Since the beginning of the operational phase, many developments and improvements have been implemented regarding the software, the data provided, and the portal. One of the most important upgrades is that the NEOCC is now independently providing data through a newly developed Orbit Determination and Impact Monitoring system, named Aegis. All the data computed by Aegis is publicly available on the NEOCC web portal, and Aegis is also used to maintain all the major services offered. The most important services comprise an orbital catalogue of all known asteroids, a list of possible future impacts with the Earth (also called Risk List), a list of forthcoming close approaches, a set of graphical toolkits, and an on-demand ephemerides service. Many of the services are also available through dedicated APIs, which can be used to automatically retrieve data. Here we give an overview of the algorithms implemented in the Aegis software, and provide a summary of the services offered by the NEOCC that are supported by Aegis.

The $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmological model has been highly successful in predicting cosmic structure and evolution, yet recent precision measurements have highlighted discrepancies, especially in the Hubble constant inferred from local and early-Universe data. Gamma-ray bursts (GRBs) present a promising alternative for cosmological measurements, capable of reaching higher redshifts than traditional distance indicators. This work leverages GRBs to refine cosmological parameters independently of the $\Lambda$CDM framework. Using the Platinum compilation of long GRBs, we calibrate the Dainotti relations-empirical correlations among GRB luminosity properties-as standard candles through artificial neural networks (ANNs). We analyze both the 2D and 3D Dainotti calibration relations, leveraging an ANN-driven Markov Chain Monte Carlo approach to minimize scatter in the calibration parameters, thereby achieving a stable Hubble diagram. This ANN-based calibration approach offers advantages over Gaussian processes, avoiding issues such as kernel function dependence and overfitting. Our results emphasize the need for model-independent calibration approaches to address systematic challenges in GRB luminosity variability, ultimately extending the cosmic distance ladder in a robust way. By addressing redshift evolution and reducing systematic uncertainties, GRBs can serve as reliable high-redshift distance indicators, offering critical insights into current cosmological tensions.

This work considers the impact of resolution in the construction of mock observations of simulated galaxies. In particular, when building mock integral field spectroscopic observations from galaxy formation models in cosmological simulations, we investigate the possible systematics that may arise given the assumption that all galaxies above some stellar mass limit will provide unbiased and meaningful observable stellar kinematics. We build a catalogue of N-body simulations to sample the range of stellar particle resolutions within the EAGLE Ref0050N0752 simulation box and examine how their observable kinematics vary relative to a higher-resolution N-body control. We use these models to compile a table of the minimum number of particles-per-pixel to reach a given uncertainty in the fitted line-of-sight velocity distribution parameters. Further, we introduce a Voronoi-binning module to the mock observation code, SimSpin, in order to meet these minimum numbers. Using EAGLE, we show the impact of this shot noise on the observed spin-ellipticity plane and the recovery of this space when observations are binned with increasing numbers of particles. In conclusion, we advise binning mock images to meet at least 200 particles-per-pixel to avoid systematically under-estimating the velocity dispersion along a given line-of-sight. We demonstrate that this is important for comparing galaxies extracted from the same simulation, as well as between simulations of varying mass resolution and observations of real galaxies.

We present interferometric observations at 1 and 3\,mm with the Atacama Large Millimeter Array (ALMA) of the free-free continuum and mm-wavelength recombination line (mRRL) emission of the ionized core (within $\lsim$130\,au) of the young Planetary Nebula (PN) candidate M\,this http URL inner regions are concealed in the vast majority of similar objects. A spectral index for the mm-to-cm continuum of $\sim$0.9 indicates predominantly free-free emission from an ionized wind, with a minor contribution from warm dust. The mm-continuum emission in M\,2-9 reveals an elongated structure along the main symmetry axis of the large-scale bipolar nebula with a C-shaped curvature surrounded by a broad-waisted component. This structure is consistent with an ionized bent jet and a perpendicular compact dusty disk. The presence of a compact equatorial disk (of radius $\sim$50\,au) is also supported by red-shifted CO and \trecem\ absorption profiles observed from the base of the receding north lobe against the compact background continuum. The redshift observed in the CO absorption profiles likely signifies gas infall movements from the disk toward a central source. The mRRLs exhibit velocity gradients along the axis, implying systematic expansion in the C-shaped bipolar outflow. The highest expansion velocities ($\sim$80\,\kms) are found in two diagonally opposed compact regions along the axis, referred to as the high-velocity spots/shells (HVS), indicating either rapid wind acceleration or shocks at radial distances of $\sim$0\farc02-0\farc04 ($\sim$15-25\,au) from the center. A subtle velocity gradient perpendicular to the lobes is also found, suggestive of rotation. Our ALMA observations detect increased brightness and broadness in the mRRLs... (abridged).

We present an investigation of Complex Organic Molecules (COMs) in the spatially resolved Keplerian disk around V883 Ori, an eruptive young star, based on a spectral survey carried out with ALMA in Band 6 (220.7$-$274.9 GHz). We identified about 3,700 molecular emission lines and discovered 23 COMs in the disk. We estimated the column densities of COMs detected through the iterative LTE line fitting method. According to our analyses, using only optically thin lines is critical to deriving the reliable column densities of COMs. Therefore, covering a large frequency range is important for the studies of COMs. The most distinct phenomenon found from the spectra of the V883 Ori disk is that nitrogen-bearing COMs other than CH$_{3}$CN are missing, whereas various oxygen-bearing COMs, except for the CH$_2$OH-bearing molecules, are detected. The missing CH$_2$OH-bearing COMs may indicate the warm water-ice dominant environment for forming COMs. We compared our results with various objects in different evolutionary stages, from Class 0 hot corinos to a Solar System comet 67P/Churyumov-Gerasimenko, to examine the effect of evolution on the COM compositions. In general, the COMs abundances relative to methanol in V883 Ori are higher than in the hot corinos and hot cores, while they are comparable to the cometary values. This may indicate the planet-forming material chemically evolves in the disk midplane after being accreted from the envelope. In addition, as found in the comet 67P/Churyumov-Gerasimenko, nitrogen might also be trapped as ammonium salt within the dust grains in the V883 Ori disk.

Galaxy cluster mergers are excellent laboratories for studying a wide variety of different physical phenomena. Such a unique system is the distant SPT-CLJ2228-5828 cluster merger located at $z\approx 0.77$. Previous analyses via Sunyaev-Zeldovich and weak lensing data suggested that the system potentially was a dissociative cluster post-merger. In this work, we use new, deep XMM-Newton data to study the hot gas in X-rays, spectroscopic Gemini data to precisely determine the redshift of the two mass concentrations, and new HST data to improve the total mass estimates of the two components. We find that SPT-CLJ2228-5828 constitutes a pre-merging, double cluster system, instead of a post-merger. The merging process of the two clusters has started with their outskirt gas colliding with a $\sim 22^{\circ}-27^{\circ}$ on the plane of the sky. We fully characterize the surface brightness, gas density, temperature, pressure, and entropy profiles of the two merging clusters. The two systems have very similar X-ray properties with a moderate cluster mass of $M_{\text{tot}}\sim (2.1-2.4)\times 10^{14}\ M_{\odot}$. A $\approx 333$ kpc long gas bridge connecting the two clusters is detected at a $5.8\sigma$ level. The baryon overdensity of the excess bridge gas is $\delta_{\text{b}}\sim (75-320)$ across the length of the bridge and its gas mass is $M_{\text{gas}}\sim 1.4\times 10^{12}\ M_{\odot}$. Gas density and temperature jumps are also found across the gas bridge, revealing the existence of a weak shock front with a Mach number $\mathcal{M}\sim 1.1$. The gas pressure and entropy are also increased at the position of the shock front. We estimate the age of the shock front to be $\lesssim 100$ Myr and its kinetic energy $\sim 2.4\times 10^{44}$ erg s$^{-1}$. SPT-CLJ2228-5828 is the first such high-$z$ pre-merger with a gas bridge and a shock front to be studied in X-rays.

We study the electrostatic energy of binary ionic mixtures (BIMs) in the form of Coulomb crystals with the main focus on ordered crystals. We consider 15 different binary bcc-like lattices, accurately calculate their electrostatic energies and approximate them by a unified equation. These results extend those available in the literature. A detailed comparison with selected previous results is made, particularly, using previous calculations in the linear mixing rule approximation. The case of disordered BIMs is also outlined. The results are expected to be useful for exploring multi-component Coulomb systems in compact stars, laboratory plasmas, and technological applications.

The radiative cooling time of hot gas in the cool cores of many galaxy clusters and massive elliptical galaxies drops in the centre to below 100 million years. The mass cooling rates inferred from simple modelling of X-ray observations of these objects are very low, indicating that either AGN feedback is tightly balanced or that soft X-rays from cooling gas are somehow hidden from view. An intrinsic absorption model developed for application to galaxy clusters is used here to search for hidden cooling flows (HCFs) in seven nearby elliptical galaxies. Mass cooling rates of 0.5-8 solar masses per year are found in each galaxy. The absorbed cooling flow luminosity is in agreement with the observed Far Infrared (FIR) luminosity in each case, indicating absorbed emission is energetically capable of emerging in the FIR band. An observed lack of agreement between HCF rates and normal star formation rates suggests the cooled material must have an alternative fate, with low-mass star formation considered as the primary outcome.

It is essential to examine the physical or chemical properties of molecular gas in starburst galaxies to reveal the underlying mechanisms characterizing starbursts. We used non-negative matrix factorization (NMF) to extract individual molecular or physical components involved in the star formation process in NGC\,253. We used images of 148 transitions from 44 different species of the ALMA large program ALCHEMI. Additionally, we included the continuum images at ALMA Bands 3 and 7 from the same dataset. For the five NMF components (NF1--5), we obtained that their distributions correspond to various basic phenomena related to star formation: i) low-density gas extended through the galactic central molecular zone (NF2), ii) shocks (NF3), iii) starburst regions (NF4), and iv) young star-forming regions (NF5). The other component (NF1) is related to excitation; three components obtained by NMF (NF3, 1, and 5) show a strong dependence upon the upper state energies of transitions, and represent low-, intermediate-, and high-excitation, respectively. We also compared our results using principal component analysis (PCA) previously applied to the same dataset. Molecular components extracted from NMF are similar to the ones obtained from PCA. However, NMF is better at extracting components associated with a single physical component, while a single component in PCA usually contains information on multiple physical components. This is especially true for features with weak intensities like emission from outflows. Our results suggest that NMF can be one of promising methods interpreting molecular line survey data, especially in the upcoming era of wide-band receivers.

The paper presents a theoretical model describing the full power spectra of synchrotron radiation generated by relativistic electrons in a turbulent magnetic field. Using the theoretical model, numerical calculations of the complete power spectra of synchrotron radiation were performed for a turbulent field generated by the harmonic method. Additionally, a model sky map was constructed, demonstrating the structure of the spatial inhomogeneity of the synchrotron radiation power distribution as seen by an observer.

This paper presents a 610 MHz radio survey covering 1.94 square degrees around the North Ecliptic Pole (NEP), which includes parts of the AKARI (ADF-N) and Euclid, Deep Fields North. The median 5-sigma sensitivity is 28 microJy beam per beam, reaching as low as 19 microJy per beam, with a synthesised beam of 3.6 x 4.1 arcsec. The catalogue contains 1675 radio components, with 339 grouped into multi-component sources and 284 isolated components likely part of double radio sources. Imaging, cataloguing, and source identification are presented, along with preliminary scientific results. From a non-statistical sub-set of 169 objects with multi-wavelength AKARI and other detections, luminous infrared galaxies (LIRGs) represent 66 percent of the sample, ultra-luminous infrared galaxies (ULIRGs) 4 percent, and sources with L_IR < 1011 L_sun 30 percent. In total, 56 percent of sources show some AGN presence, though only seven are AGN-dominated. ULIRGs require three times higher AGN contribution to produce high-quality SED fits compared to lower luminosity galaxies, and AGN presence increases with AGN fraction. The PAH mass fraction is insignificant, although ULIRGs have about half the PAH strength of lower IR-luminosity galaxies. Higher luminosity galaxies show gas and stellar masses an order of magnitude larger, suggesting higher star formation rates. For LIRGs, AGN presence increases with redshift, indicating that part of the total luminosity could be contributed by AGN activity rather than star formation. Simple cross-matching revealed 13 ROSAT QSOs, 45 X-ray sources, and 61 sub-mm galaxies coincident with GMRT radio sources.

We present the results from a series of analyses on two parametric tests of gravity that modify the growth of linear, sub-horizon matter perturbations in the $\Lambda$CDM model. The first test, known as the $(\mu,\Sigma)$ framework, modifies the Poisson and lensing equations from General Relativity (GR). The second test introduces the growth index $\gamma$, which directly affects the time evolution of matter density perturbations. Our study is motivated by results from the analysis of the Planck-PR3 2018 spectra, which indicate a preference for $\Sigma_0 \neq 0$ and $\gamma_0 > 0.55$, both of which deviate from the $\Lambda$CDM predictions at a significance level of $\sim 2.5\sigma$. To clarify the nature of these anomalous results and understand how the lensing anomaly fits into the picture, we analyze the most recent Planck-PR4 spectra extracted from the updated \texttt{NPIPE} maps. Overall, the Planck-PR4 data show better consistency with GR. The updated likelihood \texttt{Camspec} provides constraints on $\Sigma_0$ and $\gamma_0$ that are consistent with GR within $1.5\sigma$ and $2\sigma$, respectively. The updated likelihoods \texttt{HiLLiPoP} and \texttt{LoLLiPoP} show even closer agreement, with all parameter values consistent with a $\Lambda$CDM cosmology within $1\sigma$. This enhanced consistency is closely correlated with the lensing anomaly. Across the different likelihoods, the tendency of $\Sigma_0$ and $\gamma_0$ to drift towards non-standard values matches the observed preference for $A_L > 1$, both of which are significantly reduced or disappear within the Planck-PR4 data.

The detection of Type I X-ray bursts is attributed to those seen by the Astronomical Netherlands Satellite (ANS) in September 1975 from the globular cluster NGC6624 containing the X-ray source 4U1820-303. I revisit these X-ray bursts, by re-analysing data from the Soft X-ray Experiment (SXX) onboard ANS, which were stored on microfiche. Earlier accounts of X-ray bursts had been reported; the first Type I X-ray burst recorded is the one observed by Vela 5B from Cen X-4 in July 1969.

Star formation in galaxies is regulated by dynamical and thermal processes. The photoelectric effect on small dust grains usually dominates the heating of the star-forming neutral atomic gas reservoir in metal-rich galaxies, while the lower dust-to-gas mass ratio and the higher luminosity of X-ray sources in metal-poor galaxies suggest that other heating mechanisms may be at play. We calculate the relative contributions of the photoelectric effect, photoionization by UV and X-ray photons, and ionization by cosmic rays to the total heating in a sample of 37 nearby galaxies reaching down to 3% the Milky Way metallicity. We use the statistical code MULTIGRIS together with a grid of Cloudy models propagating radiation from stellar clusters and X-ray sources to the ionized and neutral gas, each galaxy being described as a statistical distribution of many 1D components. Infrared cooling lines from the interstellar medium (ISM) are used as constraints to evaluate the most likely distributions and parameters. We show that the photoelectric effect heating dominates in high-metallicity galaxies (>1/18 the Milky Way value) while cosmic rays and especially photoionization from X-rays become predominant in low-metallicity galaxies. Our models predict reasonably well the X-ray source fluxes in the 0.3-8 keV band using indirect ISM tracers, illustrating that the adopted strategy makes it possible to recover the global intrinsic radiation field properties when X-ray observations are unavailable, for instance in early universe galaxies. Finally, we show that the photoelectric effect heating efficiency on PAHs may be recovered through the [CII]+[OI] / PAH observational proxy only if the other heating mechanisms are accounted for (abridged).

Building the radio sky template are crucial for detecting the 21 cm emission line signal from the Epoch of Reionization (EoR), as well as for other cosmological research endeavors. Utilizing data from the LOFAR Two-meter Sky Survey (LoTSS) at 150 MHz, we recalibrated the luminosity function for various types of radio sources, including High Excitation Radio Galaxies (HERGs), Low Excitation Radio Galaxies (LERGs), Radio-Quiet Active Galactic Nuclei (RQ-AGNs), and Star-Forming Galaxies (SFGs). We subsequently updated the Tiered Radio Extragalactic Continuum Simulation (T-RECS) code to generate refined mock radio source catalogues. The simulated source counts from this work align more closely with observed data at redshifts greater than $z>4$. Additionally, the differential source counts in total intensity within the flux density range of $0.1-1~\mathrm{mJy}$ closely mirror actual observations. Due to our model incorporating a lower number of faint sources compared to T-RECS, it predicts a reduced power spectrum for point sources, suggesting a potential advantage in studies in low frequency band.

It has been proposed that the flat rotation curves observed at large radii in disk galaxies can be interpreted as an effect of General Relativity (GR) instead of the presence of dark matter (DM) halos in Newtonian gravity. In Ciotti (2022) the problem is rigorously explored in the special setting of the weak-field, low-velocity gravitomagnetic limit of GR. The rotation curves are obtained for purely baryonic disk models with realistic density profiles, and compared with the predictions of Newtonian gravity for the same disks, in absence of DM. The rotation curves are indistinguishable, with percentual GR corrections at all radii of the order of $\approx 10^{-6}$ or less, so that DM halos are required in gravitomagnetism as in Newtonian gravity. From a more general point of view, a list of the most urgent problems that must be addressed by any proposed GR-based alternative to the existence of DM, is given.

Filaments stand as pivotal structures within the cosmic web. However, direct detection of the cold gas content of the filaments remains challenging due to its inherent low brightness temperature. With the TNG hydrodynamical simulations, we demonstrate the effectiveness of isolating faint filament HI signal from the FAST HI intensity mapping (IM) survey through pairwise stacking of galaxies, which yields an average HI filament signal amplitude of $\sim 0.28\ {\mu{\rm K}}$ at $z\simeq 0.1$. However, our simulations reveal a non-negligible contribution from HI-rich galaxies within or near the filaments. Particularly, the faint galaxies dominantly contribute to the extra filament HI signal. Our simulation also shows that the measurement uncertainty is produced by both thermal noise and background variation caused by brightness leakage from surrounding random galaxies. Given a fixed total observation time, a wide-field HI IM survey, which includes a large number of galaxy pairs, can simultaneously reduce thermal noise to below the filament signal level and minimize background variation to a negligible level. Through the end-to-end simulation, this work demonstrates the critical role of the galaxy pairwise stacking method in future filament HI detection, outlining a road map for filament HI detection in the next-generation HI IM surveys.

Binary stars are prevalent yet challenging to detect. We present a novel approach using convolutional neural networks (CNNs) to identify binary stars from low-resolution spectra obtained by the LAMOST survey. The CNN is trained on a dataset that distinguishes binaries from single main sequence stars based on their positions on the Hertzsprung-Russell diagram. The network achieves high accuracy with an area under the receiver operating characteristic curve of 0.949 on the test set. Its performance is further validated against known eclipsing binaries (97% detection rate) and binary stars identified by radial velocity variations (92% detection rate). Applying the trained CNN to a sample of one million main sequence stars from LAMOST DR10 and Gaia DR3 yields a catalog of 468,634 binary stars. This catalog includes 115 binary stars located beyond 10 kpc from the Sun and 128 cross-matched with known exoplanet hosts from the NASA Exoplanet Archive. This new catalog provides a valuable resource for future research on the properties, formation, and evolution of binary systems, particularly for statistically characterizing large populations.

We investigate whether Neural Networks (NN) can accurately differentiate between growth-rate data of the Large Scale Structure (LSS) of the Universe, simulated via two models: a cosmological constant and cold dark matter ({\Lambda}CDM) model and a tomographic Coupled Dark Energy (CDE) model. We build an NN classifier and test its accuracy in distinguishing between cosmological models. For our dataset, we generate f{\sigma_8}(z) growth-rate observables simulating a realistic Stage IV galaxy survey-like setup, for both {\Lambda}CDM and a tomographic CDE model, for various values of the model parameters. We then optimise and train our NN with Optuna, aiming to avoid overfitting and maximising the accuracy of the trained model. We conduct our analysis for both a binary classification, comparing between {\Lambda}CDM and a CDE model where only one tomographic coupling bin is activated, and a multiclass classification scenario where all the models are combined. For the case of binary classification, we find that our NN can confidently (with > 86% accuracy) detect non-zero values of the tomographic coupling regardless of the redshift range at which coupling is activated, and at a 100% confidence level, detect the {\Lambda}CDM model. For the multiclass classification task, we find that the NN performs adequately well at distinguishing between {\Lambda}CDM, a CDE model with low redshift coupling, and a model with high redshift coupling, with 99%, 79% and 84% accuracy respectively. By leveraging the power of machine learning, our pipeline can be a useful tool to analyse growth-rate data and maximise the potential of current surveys to probe for deviations from General Relativity.

When bright solar-system objects are observed by GHz-THz regime telescopes, off-axis signals bounce around locally and re-enter the signal path with a time delay, causing sinusoidal ripples in output spectra. Ripples that are unstable over time are challenging to remove. A typical detection limit for planetary spectral lines is a fraction of order 0.001 of continuum signal, restricting searches for minor atmospheric trace-gases. Modern wideband spectra of Venus demonstrate a plethora of effects, at three example telescopes spanning nearly a factor of 50 in frequency. Characterisation of instrumental effects as families of pure sine-waves via Fourier-analysis is shown to improve dynamic range by factors of a few. An example upper limit on sulphuric acid vapour in the Venusian mesosphere, from fully-automated data-cleaning of a 3.5 GHz band containing 10 line components, goes as deep as the best previously-published limit. The most challenging cases are searches for single lines of width comparable to ripple periods. Traditional polynomial-fitting approaches can be deployed to test for false positives, to demonstrate robustness at a level of zero 'fake lines' in more than 1000 comparisons. Fourier-based data-cleaning avoids subjectivity and can be fully automated, and synthetic spectra can be injected before processing to test to what degree signals are lost in cleaning. An ideal robustness strategy is mitigation at the data-acquisition stage, e.g. using slow drifts in target-velocity with respect to the telescope to isolate a planetary line from a quasi-static instrumental ripple pattern.

We explore synergies between the Nancy Grace Roman Space Telescope High Latitude Wide Area Survey (HLWAS) and CMB experiments, specifically Simons Observatory (SO) and CMB-Stage4 (S4). Our simulated analyses include weak lensing, photometric galaxy clustering, CMB lensing, thermal SZ, and cross-correlations between these probes. While we assume the nominal 16,500 square degree area for SO and S4, we consider multiple survey designs for Roman that overlap with Rubin Observatory's Legacy Survey of Space and Time (LSST): the 2000 square degree reference survey using four photometric bands, and two shallower single-band surveys that cover 10,000 and 18,000 square degree, respectively. We find a ~2x increase in the dark energy figure of merit when including CMB-S4 data for all Roman survey designs. We further find a strong increase in constraining power for the Roman wide survey scenario cases, despite the reduction in galaxy number density, and the increased systematic uncertainties assumed due to the single band coverage. Even when tripling the already worse systematic uncertainties in the Roman wide scenarios, which reduces the 10,000 square degree FoM from 269 to 178, we find that the larger survey area is still significantly preferred over the reference survey (FoM 64). We conclude that for the specific analysis choices and metrics of this paper, a Roman wide survey is unlikely to be systematics-limited (in the sense that one saturates the improvement that can be obtained by increasing survey area). We outline several specific implementations of a two-tier Roman survey (1000 square degree with 4 bands, and a second wide tier in one band) that can further mitigate the risk of systematics for Roman wide concepts.

PDS 70 is a unique system as it hosts a protoplanetary disk with two confirmed forming planets, making it an ideal target for characterizing dust in such disks. We present new high-contrast polarimetric differential imaging of PDS 70 using the $N\_R$ filter on SPHERE/ZIMPOL, combined with archival VLT/SPHERE data across five wavelengths ($N\_R$, $VBB$, $J$, $H$, and $Ks$) spanning seven epochs over eight years. For each epoch, we corrected smearing effects from instrument resolution, analyzed azimuthal brightness profiles, and derived intrinsic disk-integrated polarized reflectivity and brightness contrasts. Our analysis reveals significant temporal variability in both integrated polarized reflectivity and azimuthal brightness profiles, suggesting variable shadowing on the outer disk from unresolved inner disk structures. Nonetheless, we observe a systematic wavelength-dependent contrast between the near and far sides of the inclined disk, highlighting the need to consider shadowing from the inner disk and surface geometry of the reflecting disk in data interpretation.

Solar energetic particle (SEP) events are associated with coronal mass ejections (CMEs) and/or solar flares. SEPs travel through the corona and interplanetary space to reach Earth, posing a radiation hazard to spacecraft and astronauts working in space and the electronics on spacecraft. Due to the distinct magnetic field configuration and solar eruption kinematic properties associated with each event, the utilization of a data-driven model becomes essential for predicting SEP hazards. In this study, we use a developed model that utilizes photospheric magnetic field measurements and CME shock observations as inputs to simulate several historical SEP events associated with fast CME speeds (>700 km/s). The model includes an SEP source term aligned with the theory of diffusive shock acceleration by the CME shock. The performance of the model is accessed by comparing simulations and observations of SEP intensity time profiles at SOHO, ACE, STEREO-A and STEREO-B. The results generally matched observations well, particularly for protons below 40.0 MeV. However, discrepancies arose for higher-energy protons, notably for the events on 2011 March 7 and 2014 February 25, where the simulation tended to overestimate the proton flux . At STEREO-A, the modeled proton intensities for the SEP events on 2013 April 11 and 2011 March 7 display a very different behavior compared to observations because of the efficient transport in longitude caused by the weak magnetic field.

Gas metallicity, ionization parameter, and gas pressure can affect the observed ratios of specific strong emission lines within galaxies. While the theoretical strong lines diagnostics for gas metallicity, ionization parameters, and gas pressure in star-forming regions are well-established, theoretical diagnostics for active galactic nuclei (AGNs) narrow line regions are still lacking. In Zhu et al. (2023), we presented a new AGN model that provides the best predictions for observations spanning the UV, optical, and infrared wavelengths. This paper presents a suite of theoretical diagnostics for the gas metallicity, ionization parameter, gas pressure, and the peak energy in AGN ionizing radiation field $E_{peak}$ for AGN narrow-line regions spanning the UV and optical wavelengths. We investigate the model dependency on the ionization parameter, gas pressure, $E_{peak}$, and the nitrogen scaling relation and make recommendations on metallicity diagnostics that are most robust against these parameters. We test our new AGN metallicity diagnostics using optical galaxy spectra from Sloan Digital Sky Survey DR16. These tests show that the metallicities measured from different diagnostics in this paper are consistent within $\sim0.3$ dex. We compare consistent HII and AGN diagnostics and demonstrate that HII and AGN diagnostics should not be used interchangeably. With a wide wavelength coverage, we anticipate that these AGN diagnostics will enable new metallicity studies of galaxies dominated by AGN.

Interacting supernovae provide key insights into the mass-loss processes of massive stars and their circumstellar environments. By analyzing their photometric and spectroscopic properties, we can study the complex interactions between ejected material and circumstellar material (CSM). This paper highlights the diversity of interacting SNe, including Types IIn, Ibn, and Icn, and explores the challenges in understanding progenitor systems, CSM structures, and late-time evolution. Advances in high-cadence observations and modeling are crucial for improving our knowledge of these stellar explosions.

We study the evolution of eccentric, equatorial extreme-mass-ratio inspirals (EMRIs) immersed in the accretion disks of active galactic nuclei. We find that single gravitational-wave observations from these systems could provide measurements with ~ 10 % relative precision of, simultaneously, the disk viscosity and mass accretion rate of the central supermassive black hole. This is possible when the EMRI transitions, within the observation time, from supersonic to subsonic motion relative to the disk gas, for eccentricities e > ~ 0.025-0.1. The estimate of the accretion rate would assist in the identification of the EMRI's host galaxy, or the observation of a direct electromagnetic counterpart, improving the chances of using these sources as cosmological sirens. Our work highlights the rich phenomenology of binary evolution in astrophysical environments and the need to improve the modelling and analysis of these systems for future gravitational-wave astronomy.

We investigate a dark matter model that couples to the standard model through a one-loop interaction with neutrinos, where the mediator particles also generate neutrino masses. We perform a global fit that incorporates dark matter relic abundance, primordial nucleosynthesis, neutrino mass, collider and indirect detection constraints. Thanks to the loop suppression, large couplings are allowed, and we find that the model parameters are constrained on all sides. Dark matter masses from 10 MeV to a few TeV are allowed, but sub-GeV masses are preferred for the model to also account for the heaviest neutrino mass. Though our results are valid for a single neutrino mass eigenstate at a time, the model and methods are generalizable to the full 3-flavor case.

Driven by recent laboratory experiments and astronomical observations, significant advances have deepened our understanding of neutron-star physics. NICER's Pulse Profile Modeling has refined our knowledge of neutron star masses and radii, while gravitational-wave detections have revealed key insights into the structure of neutron stars. Particularly relevant is the extraction of the tidal deformability by the LIGO-Virgo collaboration and the most recent determination of stellar radii by NICER, both suggesting a relatively soft equation of state (EOS) at intermediate densities. Additionally, measurements from the PREX collaboration and from pulsar timing suggest instead that the EOS is stiff in the vicinity of saturation density and at the highest densities accessible to date. But how stiff can the EOS be at these very high densities? Recent events featuring compact objects near the "lower mass gap" have raised questions about the existence of very massive neutron stars. Motivated by this finding and in light of new refinements to theoretical models, we explore the possibility that these massive objects may indeed be rapidly rotating neutron stars. We explore how rotation affects both the maximum neutron star mass and their associated radii, and discuss the implications they may have on the equation of state.

The triple-$\alpha$ reaction from ternary continuum states at off-resonant energies, $\alpha+\alpha+\alpha \rightarrow ^{12}$C, still remains in unsolved problems. This direct process is estimated with a non-adiabatic Faddeev hyper-spherical harmonics and $R$-matrix expansion method. After reviewing this model, the resultant photo-disintegration of $^{12}$C($2^+_1 \rightarrow 0^+$) is shown to be in 10$^{-15}$--10$^{-3}$ pico-barn order for $0.15 \le E \le 0.35$ MeV. This is far below the values predicted by the recent adiabatic calculations. In spite of the large difference, the derived reaction rates are illustrated to be concordant with the current evaluated rates for $0.08 \le T_9 \le 3$ including helium burning temperatures. The difference below $E = 0.20$ MeV can be seen in the rates for $T_9 \le 0.07$. In comparison with the calculations, the triple-$\alpha$ reaction rates are found to be reduced by about 10$^{-4}$ at $T_9 = 0.05$, because of an accurate description for $^8$Be break-up. Uncertainties of the rates are also estimated by examining sensitivity to 3$\alpha$ potentials. With introducing three-body $S$-factors and a resonant term, the present rates are expressed in an analytic form, and they are provided in a tabular form for astrophysical applications. To update the evaluated rates, non-resonant sequential process between $\alpha+^8$Be may be eliminated by hand.

Time-delay interferometry (TDI) is essential in space-based gravitational wave (GW) detectors, effectively reducing laser noise and improving detection precision. As one of the most promising GW detectors, the space-based detectors are able to observe the effects from GW polarizations. The detection of GW additional polarizations carries significant implications, potentially revealing deviations from general relativity and opening avenues to explore alternative gravity theories. In this study, we examine the impacts of second-generation TDI combinations on GW polarization detection by simulating LISA, Taiji, and TianQin, including realistic orbital effects such as link length and angle variations. The detector performance is assessed through sensitivity curves derived from averaged response functions, as well as signal-to-noise ratio (SNR) of binary black holes (BBHs). For massive BBHs, the $\mathcal{A}$ and $\mathcal{E}$ channels typically offer the best sensitivity, while the $X$ channel in TianQin is most effective for detecting additional polarizations. For stellar-mass BBHs, the $\alpha$ channel provides the highest SNR for vector modes in LISA and Taiji specifically for lower-mass systems, while the $\mathcal{A}$ and $\mathcal{E}$ channels are optimal for higher masses or other polarizations. TianQin consistently favors the $X$ channel for additional polarizations. Our findings emphasize the importance of selecting high-sensitivity TDI combinations to enhance detection capabilities across different polarizations, deepening our insight into GW sources and the fundamental nature of spacetime

We provide a detailed derivation of the spectral density of the stochastic background generated by the superposition of coalescing compact binaries. We show how the expression often used in the literature emerges from an average over the extrinsic parameters of the binaries (times of arrival, polarization angles, arrival directions and orbit inclinations) and how the Stokes parameters related to circular and linear polarization are set to zero by such averaging procedure. We then consider the effect of shot noise, i.e. the fact that for the superposition of a finite number of sources these averages are only approximate, and we show how it generates circular and linear polarizations (even for isotropic backgrounds) as well as spatial anisotropies, and we compute them explicitly for a realistic population of binary black holes and binary neutron stars.

At third-generation (3G) gravitational-wave detector networks, compact binaries coalescences produce a ``confusion noise'' due to unresolved sources and to the error in the reconstruction of resolved sources, that can degrade the sensitivity to cosmological backgrounds. We show how to characterize from first-principles this astrophysical confusion noise by reconstructing the resolved sources at a detector network, subtracting them from the data stream of each detector of the network, and then computing the correlation among detector pairs of these ``partially cleaned'' data streams that, beside instrumental noise, contain the unresolved sources and the error on the resolved sources. In a two-detector correlation, this residual astrophysical background then acts as an effective correlated noise. We point out that its effect, in the search for a cosmological background, must be evaluated by correlating two detectors with the filter function that optimizes the given cosmological search (and not the search for the astrophysical background itself) and therefore depends on the specific cosmological signal that we want to extract from the data. We show how to obtain search-independent upper bounds on the effect of this astrophysical confusion noise, and how to characterize its actual effect on a given cosmological search. We then apply this methodology to an example of a 3G network and, for a realistic population of mergers, we evaluate explicitly the upper bound on the effect of the astrophysical confusion noise, as well as its actual effect in the search of selected examples of power-law cosmological backgrounds and of backgrounds featuring broad peaks, as in phase transitions.

We determine parameters of the renormalization group-consistent (RG-consistent) three-flavor color-superconducting Nambu-Jona-Lasinio (NJL) model that are suited to investigate possible compact-star configurations. Our goal is to provide viable quark-matter equation of state (EoS) that can generally be used for hybrid-star constructions. To that end, we mainly focus on quark-star properties in this work. By varying the vector and diquark coupling constants, we analyze their impact on the EoS, speed of sound (SoS), the maximum diquark gap, and the mass-radius relation. In almost all configurations, a stable color-flavor-locked (CFL) phase appears in the core of the maximum-mass configurations, typically spanning several kms in radius. In other cases, the star's two-flavor color-superconducting (2SC) branch of the EoS becomes unstable before reaching the CFL transition density. At neutron-star densities, the SoS squared reaches up to 0.6 and the CFL diquark gap up to 250 MeV. We argue that adding a hadronic EoS at lower densities by performing a Maxwell construction, does not increase the maximum mass substantially, thus we use the 2 solar mass constraint to constrain the NJL model parameters that are suited for the construction of hybrid-star EoS. We construct three examples of the hybrid star model, demonstrating that there is room for different CSC compositions. The hybrid EoS obtained in this way can have no 2SC matter or different ratios of 2SC and CFL quark matter in the core. We show that early hadron-quark transitions are possible that can modify the tidal deformability at 1.4 solar mass. We will provide tabulated EoS of the RG-consistent NJL model for these three parameter sets. We find that these EoS are consistent with the imposed constraints from astrophysics and perturbative QCD. They allow for different hybrid-star scenarios with a hadronic EoS that is soft at low densities.

In this work, we analyze the Einstein-scalar-Gauss-Bonnet (EsGB) theory of gravity in a cosmological context using the formalism of dynamical systems. We obtain the equations of motion of the theory and introduce an appropriate set of dynamical variables to allow for a direct comparison with the results from General Relativity (GR). We observe that the cosmological phase space features the same set of fixed points as in standard GR, i.e., radiation-dominated, matter-dominated, curvature-dominated, and exponentially-accelerated solutions independently of the values of the coupling function and the scalar field. Furthermore, the radiation-dominated fixed points are repellers and the exponentially accelerated fixed points are attractors in the phase space, thus allowing for cosmological solutions behaving qualitatively similar to the $\Lambda$CDM model, i.e., transitioning from a radiation-dominated phase into a matter-dominated phase, and later into a late-time cosmic acceleration phase supported by the scalar field potential. Following a reconstruction method through which we produce the cosmological solutions in the GR limit of the theory and introduce them into the general EsGB dynamical system, a numerical integration of the dynamical system shows that the EsGB theory provides cosmological solutions indistinguishable from those of the standard $\Lambda$CDM model, compatible with the current observations from the Planck satellite and weak-field solar system dynamics, while maintaining the scalar field and the coupling function finite and regular throughout the entire time evolution.