Papers for Wednesday, Jan 10 2024

JWST has detected many overmassive galactic systems at $z > 4$, where the mass of the black hole, $M_\bullet$, is $10-100$ times larger than expected from local relations, given the host's stellar mass, $M_\star$. This Letter presents a model to describe these overmassive systems in the high-$z$ Universe. We suggest that the black hole mass is the main driver of high-$z$ star formation quenching. SMBHs globally impact their high-$z$ galaxies because their hosts are physically small, and the black holes have duty cycles close to unity at $z > 4$. In this regime, we assume that black hole mass growth is regulated by the quasar's output, while stellar mass growth is quenched by it and uncorrelated to the global properties of the host halo. We find that the ratio $M_\bullet/M_\star$ controls the average star formation efficiency: if $M_\bullet/M_\star > 8\times 10^{18} (n \Lambda/f_{edd})[(\Omega_b M_h)/(\Omega_m M_\star) - 1]$, then the galaxy is unable to form stars efficiently. Once this ratio exceeds the threshold, a runaway process brings the originally overmassive system towards the local $M_\bullet - M_\star$ relation. Furthermore, the $M_\bullet - M_\star$ relation evolves with redshift as $\propto (1+z)^{5/2}$. At $z \sim 5$, we find an overmassive factor of $\sim 55$, in excellent agreement with current JWST data and the high-$z$ relation inferred from those. Extending the black hole horizon farther in redshift and lower in mass will test this model and improve our understanding of the early co-evolution of black holes and galaxies.

Observations of ionised H$\alpha$ gas in disc galaxies with high star formation rates have ubiquitous and significant line broadening with widths $\sigma_{\rm H\alpha}\gtrsim 50-100\ {\rm km\ s^{-1}}$. To understand whether this broadening reflects gas turbulence within the interstellar medium (ISM) of galactic discs, or arises from off-the-plane emission in mass-loaded galactic winds, we perform radiation hydrodynamic (RHD) simulations of isolated Milky Way-mass disc galaxies in a gas-poor (low-redshift) and gas rich (high-redshift) condition and create mock H$\alpha$ emission line profiles. We find that the vast majority of the ${\rm H\alpha}$ emission is confined within the ISM, with extraplanar gas contributing mainly to the extended profile wings. This substantiates the \Halpha emission line as a tracer of mid-plane disc dynamics. We investigate the relative contribution of diffuse and dense ${\rm H\alpha}$ emitting gas, corresponding to DIG ($\rho \lesssim 0.1\ {\rm cm^{-3}}$, $T\sim 8\ 000\ {\rm K}$) and HII regions ($\rho \gtrsim 10\ {\rm cm^{-3}}$, $T\sim 10\ 000\ {\rm K}$), respectively, and find that DIG contributes $\lesssim 10 \%$ of the total ${\rm L}_{\rm H\alpha}$. However, the DIG can reach upwards of $\sigma_{\rm H\alpha} \sim 60-80\ {\rm km\ s^{-1}}$ while the HII regions are much less turbulent $\sigma_{\rm H\alpha}\sim10-40\ {\rm km\ s^{-1}}$. This implies that the $\sigma_{\rm H\alpha}$ observed using the full ${\rm H\alpha}$ emission line is dependent on the relative ${\rm H\alpha}$ contribution from DIG/HII regions and a larger $f_{\rm DIG}$ would shift $\sigma_{\rm H\alpha}$ to higher values. Finally, we show that $\sigma_{\rm H\alpha}$ evolves, in both the DIG and HII regions, with the galaxy gas fraction. Our high-redshift equivalent galaxy is roughly twice as turbulent, except for in the DIG which has a more shallow evolution.

We use the recent statistics of dual active galactic nuclei (AGN) in the $James \ Webb \ Space \ Telescope$ (JWST) data at $z \sim 3.4$ to address two aspects of the feedback and evolution scenarios of supermassive black hole binaries (SMBHB). We find that the JWST data provide evidence for the members of a binary black hole being 'lit' at the same time, rather than independently -- a scenario which is consistent with gas-rich mergers being responsible for concurrent AGN activity. This conclusion is supported by the recent NANOGrav Pulsar Timing Array (PTA) measurements, whose upper limits on the stochastic gravitational wave strain amplitude lie below those expected from extrapolating the dual AGN fraction. The results indicate either a 'stalling' of the binaries at the separations probed by NANOGrav, or rapid gas-driven inspirals.

Of order one in 10^3 quasars and high-redshift galaxies appears in the sky as multiple images as a result of gravitational lensing by unrelated galaxies and clusters that happen to be in the foreground. While the basic phenomenon is a straightforward consequence of general relativity, there are many non-obvious consequences that make multiple-image lensing systems (aka strong gravitational lenses) remarkable astrophysical probes in several different ways. This article is an introduction to the essential concepts and terminology in this area, emphasizing physical insight. The key construct is the Fermat potential or arrival-time surface: from it the standard lens equation, and the notions of image parities, magnification, critical curves, caustics, and degeneracies all follow. The advantages and limitations of the usual simplifying assumptions (geometrical optics, small angles, weak fields, thin lenses) are noted, and to the extent possible briefly, it is explained how to go beyond these. Some less well-known ideas are discussed at length: arguments using wavefronts show that much of the theory carries over unchanged to the regime of strong gravitational fields; saddle-point contours explain how even the most complicated image configurations are made up of just two ingredients. Orders of magnitude, and the question of why strong lensing is most common for objects at cosmological distance, are also discussed. The challenges of lens modeling, and diverse strategies developed to overcome them, are discussed in general terms, without many technical details.

Interpreting the observations of exoplanet atmospheres to constrain physical and chemical properties is typically done using Bayesian retrieval techniques. Because these methods require many model computations, a compromise is made between model complexity and run time. Reaching this compromise leads to the simplification of many physical and chemical processes (e.g. parameterised temperature structure). Here we implement and test sequential neural posterior estimation (SNPE), a machine learning inference algorithm, for exoplanet atmospheric retrievals. The goal is to speed up retrievals so they can be run with more computationally expensive atmospheric models, such as those computing the temperature structure using radiative transfer. We generate 100 synthetic observations using ARCiS (ARtful Modeling Code for exoplanet Science, an atmospheric modelling code with the flexibility to compute models in varying degrees of complexity) and perform retrievals on them to test the faithfulness of the SNPE posteriors. The faithfulness quantifies whether the posteriors contain the ground truth as often as we expect. We also generate a synthetic observation of a cool brown dwarf using the self-consistent capabilities of ARCiS and run a retrieval with self-consistent models to showcase the possibilities that SNPE opens. We find that SNPE provides faithful posteriors and is therefore a reliable tool for exoplanet atmospheric retrievals. We are able to run a self-consistent retrieval of a synthetic brown dwarf spectrum using only 50,000 forward model evaluations. We find that SNPE can speed up retrievals between $\sim2\times$ and $\geq10\times$ depending on the computational load of the forward model, the dimensionality of the observation, and the signal-to-noise ratio of the observation. We make the code publicly available for the community on Github.

The Chinese Space Station Telescope (CSST) slitless spectroscopic survey will observe objects to a limiting magnitude of $\sim23$\,mag (5$\sigma$, point sources) in $U$, $V$, and $I$ over 17,500 square degrees. The spectroscopic observations are expected to be highly efficient and complete for mapping galaxies over $0<z<1$ with secure redshift measurements at spectral resolutions of $R\sim200$, providing unprecedented datasets for cosmological studies. To quantitatively examine the survey potential, we develop a software tool, namely the CSST Emulator for Slitless Spectroscopy (\texttt{CESS}), to quickly generate simulated 1-D slitless spectra with limited computing resources. We introduce the architecture of \texttt{CESS} and the detailed process of creating simulated CSST slitless spectra. The extended light distribution of a galaxy induces the self-broadening effect on the 1-D slitless spectrum. We quantify the effect using morphological parameters: S\'ersic index, effective radius, position angle, and axis ratio. Moreover, we also develop a module for \texttt{CESS} to estimate the overlap contamination rate for CSST grating observations of galaxies in galaxy clusters. Applying \texttt{CESS} to the high-resolution model spectra of a sample of $\sim140$ million galaxies with $m_{z}<21$\,mag selected from the DESI LS DR9 catalogue, we obtain the simulated CSST slitless spectra. We examine the dependence of measurement errors on different types of galaxies due to instrumental and observational effects and quantitatively investigate the redshift completeness for different environments out to $z\sim1$. Our results show that the CSST spectroscopy is able to provide secure redshifts for about one-quarter of the sample galaxies.

The current orbital geometries of exoplanet systems offer a fossilized record of the systems' dynamical histories. A particularly rich set of dynamical mechanisms is available to exoplanets residing in multi-star systems, which may have their evolution shaped by the gravitational influence of bound stellar companions. In this work, we examine the joint distribution of stellar obliquities and orbital orientations for transiting exoplanets residing within astrometrically resolved binary and triple-star systems. We leverage existing constraints on stellar obliquities in exoplanet systems, together with astrometric measurements from Gaia DR3, to uncover a set of fully-aligned, "orderly" exoplanet systems that exhibit evidence of both spin-orbit and orbit-orbit alignment. We also find evidence that the observed distribution of orbit-orbit orientations in our sample is more strongly peaked toward alignment than an isotropic distribution. Our results may be indicative of efficient viscous dissipation by nodally recessing protoplanetary disks, demonstrating a regime in which stellar companions produce and maintain order in planetary systems, rather than enhancing misalignments.

Modern machine learning will allow for simulation-based inference from reionization-era 21cm observations at the Square Kilometre Array. Our framework combines a convolutional summary network and a conditional invertible network through a physics-inspired latent representation. It allows for an optimal and extremely fast determination of the posteriors of astrophysical and cosmological parameters. The sensitivity to non-Gaussian information makes our method a promising alternative to the established power spectra.

Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around $\approx 40-70\%$ of the Galactic CVs are post-period minimum systems referred to as "period bouncers", only a few dozen confirmed systems are known. We report the study and characterisation of a new eclipsing CV, SRGeJ041130.3+685350 (SRGeJ0411), discovered from a joint SRG/eROSITA and ZTF program. The optical spectrum of SRGeJ0411 shows prominent hydrogen and helium emission lines, typical for CVs. We obtained optical high-speed photometry to confirm the eclipse of SRGeJ0411 and determine the orbital period to be $P_\textrm{orb} \approx 97.530$ minutes. The spectral energy distribution suggests that the donor has an effective temperature of $\lesssim 1,800$ K. We constrain the donor mass with the period--density relationship for Roche-lobe-filling stars and find that $M_\textrm{donor} \lesssim 0.04\ M_\odot$. The binary parameters are consistent with evolutionary models for post-period minimum CVs, suggesting that SRGeJ0411 is a new period bouncer. The optical emission lines of SRGeJ0411 are single-peaked despite the system being eclipsing, which is typically only seen due to stream-fed accretion in polars. X-ray spectroscopy hints that the white dwarf in SRGeJ0411 could be magnetic, but verifying the magnetic nature of SRGeJ0411 requires further investigation. The lack of optical outbursts has made SRGeJ0411 elusive in previous surveys, and joint X-ray and optical surveys highlight the potential for discovering similar systems in the near future.

From the chemo-dynamical properties of tidal debris in the Milky Way, it has been inferred that the dwarf satellites that have been disrupted had different chemical abundances from their present-day counterparts of similar mass that survive today, specifically, they had lower [Fe/H] and higher [Mg/Fe]. Here we use the ARTEMIS simulations to study the relation between the chemical abundances of disrupted progenitors of MW-mass galaxies and their stellar mass, and the evolution of the stellar mass - metallicity relations (MZR) of this population with redshift. We find that these relations have significant scatter, which correlates with the accretion redshifts ($z_{\rm acc}$) of satellites, and with their cold gas fractions. We investigate the MZRs of dwarf populations accreted at different redshifts and find that they have similar slopes, and also similar with the slope of the MZR of the surviving population ($\approx 0.32$). However, the entire population of disrupted dwarfs displays a steeper MZR, with a slope of $\approx 0.48$, which can be explained by the changes in the mass spectrum of accreted dwarf galaxies with redshift. We find strong relations between the (mass-weighted) $\langle z_{\rm acc} \rangle$ of the disrupted populations and their global chemical abundances ($\langle$[Fe/H]$\rangle$ and $\langle$[Mg/Fe]$\rangle$), which suggests that chemical diagnostics of disrupted dwarfs can be used to infer the types of merger histories of their hosts. For the case of the MW, our simulations predict that the bulk of the disrupted population was accreted at $\langle z_{\rm acc} \rangle \approx 2$, in agreement with other findings. We also find that disrupted satellites form and evolve in denser environments, closer to their hosts, than their present-day counterparts.

We compare the performance of the popular Bondi-Hoyle-Lyttleton (BHL) accretion scheme with a simple mass-flux scheme applied to stellar-mass black holes (BHs) across six levels of increasing spatial resolution. Simulating the formation of black holes within cosmological mini-haloes at $z \sim 20$, we investigate both the core-collapse and direct-collapse channels, which result in BHs of initial mass $10.8 \, \text{M}_\odot$ and $270 \, \text{M}_\odot$ respectively. Our explicit focus on the stellar-mass range pushes the maximum resolution down to sub-$10^{-3} \, \text{pc}$ regimes, where more complicated gas dynamics are resolved. We observe efficient growth and rotationally supported, $\sim$$10^{-1} \, \text{pc}$-scale discs around all $270 \, \text{M}_\odot$ BHs independent of resolution and accretion scheme, though clumps, bars, and spiral arm structures impact stability at high resolution. We analyse the effect of these instabilities on the accretion cycle. In contrast, all bar one of the $10.8 \, \text{M}_\odot$ BHs fail to attract a disc and experience modest growth, even when characteristic scales of accretion and dynamical friction are reasonably resolved. While the two accretion schemes somewhat converge in mass growth for the $270 \, \text{M}_\odot$ case over $1 \, \text{Myr}$, the greater degree of gas fragmentation induces more randomness in the evolution of the $10.8 \, \text{M}_\odot$ BHs. We conclude that early universe black holes of $M_{\text{BH}} \sim 10^1 \, \text{M}_\odot$ struggle to grow even in gas-rich environments without feedback in comparison to seeds of $M_{\text{BH}} \sim 10^2 \, \text{M}_\odot$, and the latter exhibit convergent growth histories across accretion schemes below a spatial resolution of $1 \times 10^{-3} \, \text{pc}$.

Relics of ancient accretion events experienced by the Milky Way are predominantly located within the stellar halo of our Galaxy. However, debris from different objects display overlapping distributions in dynamical spaces, making it extremely challenging to properly disentangle their contribution to the build-up of the Galaxy. To shed light on this chaotic context, we started a program aimed at the homogeneous chemical tagging of the local halo of the Milky Way, focusing on the component in retrograde motion, since this is expected to host a large fraction of stars accreted from past mergers. The A Walk on the Retrograde Side (WRS) project targets retrograde halo stars in the Solar Neighborhood having accurate $6$-D phase space information available, measuring the precise chemical abundance of several chemical elements from high-resolution spectroscopy. In this first paper, we present the project and the analysis of high-resolution spectra obtained with UVES at VLT and PEPSI at LBT for $186$ stars. Accurate radial velocity and chemical abundance of several elements have been obtained for all the target stars. In particular we focus on the chemical composition of a specific subset of substructures identified dynamically in the literature. Our study reveals that two among the more recently discovered structures in the retrograde halo, namely Antaeus / L-RL$64$ and ED-$3$, have identical chemical patterns and similar integrals of motion, suggesting a common origin. In turn, the abundance patterns of this unified system differ from that of Gaia-Enceladus, confirming that it is an independent structure. Finally, Sequoia exhibits a different chemistry with respect to that of Gaia-Enceladus at $\mathrm{[Fe/H]} < -1.5$ dex, showcasing an excess of stars with lower Mg and Ca in the common metallicity range.

(Abridged) Polarized emission from masers is an excellent tool to study magnetic fields in maser sources. The linear polarization of most masers is understood as an interplay of maser saturation and anisotropic pumping. However, for the latter mechanism, no quantitative modeling has been presented yet. We present a comprehensive model of maser polarization, including quantitative modeling of both anisotropic pumping and the effects of maser saturation on the polarization of masers. We extend regular maser excitation modeling with a dimension that describes the molecular population alignments, as well as including the linear polarization dimension to the radiative transfer. The results of the excitation analysis yield the anisotropic pumping and decay parameters, that are subsequently used in one-dimensional proper maser polarization radiative transfer modeling. We present the anisotropic pumping parameters for a variety of transitions from class I CH$_3$OH masers, H$_2$O masers and SiO masers. SiO masers are highly anisotropically pumped due to them occurring in the vicinity of a late-type star, that irradiates the maser region with a strong directional radiation field. Class I CH$_3$OH masers and H$_2$O masers occur in association with shocks, and they are modestly anisotropically pumped due to the anisotropy of the excitation region. Our modeling constitutes the first quantitative constraints on the anisotropic pumping of masers. We find that anisotropic pumping can explain the high polarization yields of SiO masers, as well as the modest polarization of unsaturated class I CH$_3$OH masers. We predict that the $183$ GHz H$_2$O maser is strongly anisotropically pumped. Finally, we outline a mechanism through which non-Zeeman circular polarization is produced, when the magnetic field changes direction along the propagation through an anisotropically pumped maser.

GN-z11 is the highest redshift galaxy spectroscopically confirmed with the Hubble Space Telescope (HST). Previous measurements of the effective radius of GN-z11 utilized galfit, which is not optimized to measure structural parameters for such a faint, distant object. Using a new software program called forcepho on HST data for the first time, we derive a size from images in the F160W band obtained both from the complete CANDELS survey and additional midcycle observations in order to contribute to the knowledge base on the size evolution, size-luminosity, and size-mass relation of early galaxies. We find a half-light radius mean of 0''.036 \(\pm\) 0''.006 corresponding to a physical size of 0.15 \(\pm\) 0.025 kpc. This size, smaller than the point spread function, is dramatically smaller than previous estimates with shallower HST data using galfit but consistent with recent measurements using forcepho on new JWST data arXiv:2302.07234. Such a small size, combined with the JWST/NIRSpec spectroscopic observations arXiv:2305.12492, suggests that GN-z11's high luminosity is dominated by an AGN.

Fine-tuning studies whether some physical parameters, or relevant ratios between them, are located within so-called life-permitting intervals of small probability outside of which carbon-based life would not be possible. Recent developments have found estimates of these probabilities that circumvent previous concerns of measurability and selection bias. However, the question remains if fine-tuning can indeed be known. Using a mathematization of the epistemological concepts of learning and knowledge acquisition, we argue that most examples that have been touted as fine-tuned cannot be formally assessed as such. Nevertheless, fine-tuning can be known when the physical parameter is seen as a random variable and it is supported in the nonnegative real line, provided the size of the life-permitting interval is small in relation to the observed value of the parameter.

Bursty star formation -- a key prediction for high-redshift galaxies from cosmological simulations explicitly resolving stellar feedback in the interstellar medium -- has recently been observed to prevail among galaxies at redshift $z \gtrsim 6$. Line intensity mapping (LIM) of the 158 $\mu$m [C II] line as a star formation rate indicator offers unique opportunities to tomographically constrain cosmic star formation at high redshift, as an alternative to observations of individually detected galaxies. To understand effects of bursty star formation on [C II] LIM, which remain unexplored in previous studies, we present an analytic modeling framework for high-$z$ galaxy formation and [C II] LIM signals that accounts for bursty star formation histories induced by delayed supernova feedback. We use it to explore and characterize how bursty star formation can impact and thus complicate the interpretation of the [C II] luminosity function and power spectrum. Our simple analytic model indicates that bursty star formation is most important for low halo masses, and in the power spectrum it can create a substantial excess in the large-scale clustering term. This distortion results in a power spectrum shape which cannot be explained by invoking a mass-independent scatter. We conclude that burstiness must be accounted for when modeling and analyzing [C II] datasets from the early universe, and that in the extreme, the signature of burstiness may be detectable with first-generation experiments such as TIME, CONCERTO, and CCAT-DSS.

Acceleration and transport of solar energetic particles (SEPs) causes their abundances, measured at constant velocity, to be enhanced or suppressed as a function of each ion's magnetic rigidity, and hence its atomic mass-to-charge ratio A/Q. Ion charges, in turn, depend upon source electron temperature. In small "impulsive" SEP events, arising from solar jets, acceleration during magnetic reconnection causes steep power-law abundance enhancements. These impulsive SEP events can have 1000-fold enhancements of heavy elements from sources at ~2.5 MK, and similar enhancements of 3He/4He and of streaming electrons that drive type-III radio bursts. Gamma-ray lines show that solar flares also accelerate 3He-rich ions, but their electrons and ions remain trapped on magnetic loops so they dissipate their energy as X-rays, gamma-rays, heat, and light. "Gradual" SEPs accelerated at shock waves, driven by fast coronal mass ejections (CMEs), can show power-law abundance enhancements or depressions, even with see ions from the ambient solar corona. In addition, shocks can reaccelerate seed particles from residual impulsive SEPs with their pre-existing signature heavy-ion enhancements. Different patterns of abundance often show that heavy elements are dominated by a different source from that of H and He. Nevertheless, the SEP abundances averaged over many large events define the abundances of the corona itself, which is found to differ from the solar photosphere as a function of the first ionization potential (FIP) since ions, with FIP < 10 eV, are driven upward by forces of electromagnetic waves which neutral atoms, with FIP > 10 eV, cannot feel. Thus, SEPs provide a measurement of element abundances in the solar corona, distinct from the solar wind, and may even better define the photosphere for some elements.

We leveraged low to high frequency ALMA observations to investigate the cold gas and dust in 10 QSOs at $z\gtrsim 6$. Our analysis of the CO(6-5) and CO(7-6) emission lines in the selected QSOs provided insights into their molecular gas masses, averaging around $10^{10}\ \rm M_\odot$, consistent with typical values for high-redshift QSOs. Proprietary and archival ALMA observations in bands 8 and 9 enabled, for the first time, precise constraints on the dust properties and SFR of 4 QSOs in our sample. Examination of the redshift distribution of dust temperatures revealed a general trend of increasing $T_{\rm dust}$ with redshift, which is in agreement with theoretical expectations. On the contrary, investigation of the dust emissivity index indicated a generally constant value with redshift, suggesting shared dust properties among sources. We computed a mean cold dust SED considering all 10 QSOs that offers a comprehensive view of high-$z$ QSO's dust properties. The QSOs marked by more intense supermassive black hole growth (HYPERION QSOs) showed -- on average -- smaller dust masses and higher gas-to-dust ratios, while having $\rm H_2$ gas reservoirs consistent with other QSOs at the same redshift. Beyond supporting the paradigm that high-$z$ QSOs reside in highly star-forming galaxies, our findings portrayed an interesting evolutionary path at $z>6$. Our study suggested that QSOs at $z\gtrsim 6$ are undergoing rapid galaxy growth, potentially regulated by strong outflows. In the $M_{\rm BH}-M_{\rm dyn}$ plane, our high-$z$ QSOs lie above the relation measured locally. Their inferred evolutionary path portends a convergence towards the massive end of the local relation, supporting their candidacy as progenitors of local massive galaxies. The observed pathway involves intense BH growth followed by substantial galaxy growth, in contrast with a symbiotic growth scenario.

This study explores the impact of a turbulent scattering mechanism, akin to those influencing solar and galactic cosmic rays propagating in the interplanetary medium, on the population of suprathermal electrons in the solar wind. We employ a Fokker-Planck equation to model the radial evolution of electron pitch angle distributions under the action of magnetic focusing, which moves the electrons away from isotropy, and of a diffusion process that tends to bring them back to it. We compare the steady-state solutions of this Fokker-Planck equation with data obtained from the Solar Orbiter and Parker Solar Probe missions and find a remarkable agreement, varying the turbulent mean free path as the sole free parameter in our model. The obtained mean free paths are of the order of the astronomical unit, and display weak dependence on electron energy within the $100$ eV to $1$ keV range. This value is notably lower than Coulomb collision estimates but aligns well with observed mean free paths of low-rigidity solar energetic particles events. The strong agreement between our model and observations leads us to conclude that the hypothesis of turbulent scattering at work on electrons at all heliospheric distances is justified. We discuss several implications, notably the existence of a low Knudsen number region at large distances from the Sun, which offers a natural explanation for the presence of an isotropic ``halo'' component at all distances from the Sun -- electrons being isotropized in this distant region before travelling back into the inner part of the interplanetary medium.

Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a synthetic UV to near-infrared broadband image atlas for a complete stellar-mass selected sample of 1154 galaxies extracted from the TNG50 cosmological simulation at $z=0$. The images have a high spatial resolution (100 pc) and a wide field of view (160 kpc). In addition to the dust-obscured images, we also release dust-free images and physical parameter property maps with matching characteristics. As a sanity check and preview application we discuss the UVJ diagram of the galaxy sample. We investigate the effect of dust attenuation on the UVJ diagram and find that it affects both the star-forming and the quiescent galaxy populations. The quiescent galaxy region is polluted by younger and star-forming highly inclined galaxies, while dust attenuation induces a separation in inclination of the star-forming galaxy population, with low-inclination galaxies remaining at the blue side of the diagram and high-inclination galaxies systematically moving towards the red side. This image atlas can be used for a variety of other applications, including galaxy morphology studies and the investigation of local scaling relations. We publicly release the images and parameter maps, and we invite the community to use them.

Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and near-infrared (NIR) wavelengths. We find that, on average, the effective radius in every band exceeds the stellar mass effective radius, and that this excess systematically decreases with increasing wavelength. The optical g-band (NIR Ks-band) effective radius is on average 58% (13%) larger than the stellar mass effective radius. Effective radii measured from dust-obscured images are systematically larger than those measured from dust-free images, although the effect is limited (8.7% in the g-band, 2.1% in the Ks-band). We find that stellar population gradients are the dominant factor (about 80%) in driving the wavelength dependence of the effective radius, and that differential dust attenuation is a secondary factor (20%). Comparing our results to recent observational data, we find offsets in the absolute values of the median effective radii, up to 50% for the population of blue galaxies. We find better agreement in the slope of the wavelength dependence of the effective radius, with red galaxies having a slightly steeper slope than green-blue galaxies. Comparing our effective radii with those of galaxies from the Siena Galaxy Atlas in separate bins in z-band absolute magnitude and g-z colour, we find excellent agreement for the reddest galaxies, but again significant offsets for the blue populations: up to 70% for galaxies around Mz=-21.5. This difference in median effective radius for the bluer galaxies is most probably due to (abridged...).

High-mass stars have an enormous influence on the evolution of the interstellar medium in galaxies, so it is important that we understand how they form. We examine the central clumps within a sample of seven infrared-dark clouds (IRDCs) with a range of masses and morphologies. We use 1 pc-scale observations from NOEMA and the IRAM 30-m telescope to trace dense cores with 2.8 mm continuum, and gas kinematics in C$^{18}$O, HCO$^+$, HNC, and N$_2$H$^+$ ($J$=1$-$0). We supplement our continuum sample with six IRDCs observed at 2.9 mm with ALMA, and examine the relationships between core- and clump-scale properties. We have developed a fully-automated multiple-velocity component hyperfine line-fitting code called mwydyn which we employ to trace the dense gas kinematics in N$_2$H$^+$ (1$-$0), revealing highly complex and dynamic clump interiors. We find that parsec-scale clump mass is the most important factor driving the evolution; more massive clumps are able to concentrate more mass into their most massive cores - with a log-normally distributed efficiency of around 9% - in addition to containing the most dynamic gas. Distributions of linewidths within the most massive cores are similar to the ambient gas, suggesting that they are not dynamically decoupled, but are similarly chaotic. A number of studies have previously suggested that clumps are globally collapsing; in such a scenario, the observed kinematics of clump centres would be the direct result of gravity-driven mass inflows that become ever more complex as the clumps evolve, which in turn leads to the chaotic mass growth of their core populations.

The tilt of the bipolar magnetic region (BMR) is crucial in the Babcock-Leighton process for the generation of the poloidal magnetic field in the Sun. We extend the work of Jha et al. (2020) and analyze the recently reported tracked BMR catalogue based on AutoTAB Sreedevi et al. (2023) from Michelson Doppler Imager (1996-2011) and Helioseismic and Magnetic Imager (2010-2018). Using the tracked information of BMRs based on AutoTAB, we confirm that the distribution of bmax reported by Jha et al. (2020) is not because of the BMRs are picked multiple times at the different phases of their evolution instead it is also present if we consider each BMRs only once. Moreover, we find that the slope of Joy's law initially increases slowly with the increase of bmax. However, when bmax>2.5 kG, gamma_0 decreases. The decrease of observed $gamma_0$ with bmax provides a hint to a nonlinear tilt quenching in the Babcock-Leighton process.

Following Shin et al. (2023b), which is a part of the Systematic KMTNet Planetary Anomaly Search series (i.e., a search for planets in the 2016 KMTNet prime fields), we conduct a systematic search of the 2016 KMTNet sub-prime fields using a semi-machine-based algorithm to identify hidden anomalous events missed by the conventional by-eye search. We find four new planets and seven planet candidates that were buried in the KMTNet archive. The new planets are OGLE-2016-BLG-1598Lb, OGLE-2016-BLG-1800Lb, MOA-2016-BLG-526Lb, and KMT-2016-BLG-2321Lb, which show typical properties of microlensing planets, i.e., giant planets orbit M dwarf host stars beyond their snow lines. For the planet candidates, we find planet/binary or 2L1S/1L2S degeneracies, which are an obstacle to firmly claiming planet detections. By combining the results of Shin et al. (2023b) and this work, we find a total of nine hidden planets, which is about half the number of planets discovered by eye in 2016. With this work, we have met the goal of the systematic search series for 2016, which is to build a complete microlensing planet sample. We also show that our systematic searches significantly contribute to completing the planet sample, especially for planet/host mass ratios smaller than $10^{-3}$, which were incomplete in previous by-eye searches of the KMTNet archive.

Outflows play a pivotal role in star formation as one of its most visible markers and a means of transporting mass, momentum, and angular momentum from the infalling gas into the surrounding molecular cloud. Their wide reach (at least $\sim 1000$s of au) is a contrast to typical disk sizes ($\sim 10-100$ au). We employ high-resolution three-dimensional nested-grid nonideal magnetohydrodynamic (MHD) simulations to study outflow properties in the Class 0 phase. We find that no disk wind is driven from the extended centrifugal disk that has weak magnetic coupling. The low-velocity winds emerge instead from the infalling magnetic pseudodisk. Much of the disk actually experiences an infall of matter rather than outflowing gas. Some of the pseudodisk wind (PD-wind) moves inward to regions above the disk and either falls onto the disk or proceeds upward. The upward flow gives the impression of a disk wind above a certain height even if the gas is originally emerging from the pseudodisk. The PD-wind has the strongest flow coming from a disk interaction zone that lies just outside the disk and is an interface between the inwardly advected magnetic field of the pseudodisk and the outwardly diffusing magnetic field of the disk. The low-velocity wind exhibits the features of a flow driven by the magnetic pressure gradient force in some regions and those of a magnetocentrifugal wind in other regions. We interpret the structure and dynamics of the outflow zone in terms of the basic physics of gravity, angular momentum, magnetic fields, and nonideal MHD.

What is the numerical reproducibility of a stellar system (including its discs) when evolving only a sub-set of (partially-evolved) smoothed particle hydrodynamics (SPH) particles? To investigate this, we modelled the evolution of 29 star forming clumps that were extracted from our previous simulations that investigated the formation and early evolution of low-mass star clusters. These clumps were evolved using a three-dimensional smoothed particle radiation magnetohydrodynamics code, where we included or excluded non-ideal magnetohydrodynamics to match the cluster simulation. While star formation proceeded as expected, we were unable to identically reproduce any of the systems present at the end of the cluster simulations. However, the final distributions of stellar mass, stellar system mass, disc mass, and disc radii were reproduced statistically; unfortunately, the distribution of average magnetic field strengths in the discs was not reproduced statistically, but this may be a result of our updated algorithms governing the evolution of the magnetic field. Therefore, given that our clumps yield stellar masses that are statistically similar to those in the original low-mass star clusters, we have demonstrated that we can statistically reproduce systems (aside from their magnetic field strength) by evolving a subset of SPH particles. Therefore, clumps such as these can be used as initial conditions to investigate the formation of isolated stars from less-contrived initial environments.

Directly observing exoplanets with coronagraphs is impeded by the presence of speckles from aberrations in the optical path, which can be mitigated in hardware with wavefront control as well as in post-processing. This work explores using an instrument model in post-processing to separate astrophysical signals from residual aberrations in coronagraphic data. The effect of wavefront error (WFE) on the coronagraphic intensity consists of a linear contribution and a quadratic contribution. When either of the terms is much larger than the other, the instrument response can be approximated by a transfer matrix mapping WFE to detector plane intensity. From this transfer matrix, a useful projection onto instrumental modes that removes the dominant error modes can be derived. We apply this projection to synthetically generated Roman Space Telescope hybrid Lyot coronagraph (HLC) data to extract "robust observables," which can be used instead of raw data for applications such as detection testing. The projection improves planet flux ratio detection limits by about 28% in the linear regime and by over a factor of 2 in the quadratic regime, illustrating that robust observables can increase sensitivity to astrophysical signals and improve the scientific yield from coronagraphic data. While this approach does not require additional information such as observations of reference stars or modulations of a deformable mirror, it can and should be combined with these other techniques, acting as a model-informed prior in an overall post-processing strategy.

We report observational evidence of highly turbulent ionized gas kinematics in a sample of 20 Lyman continuum (LyC) emitters (LCEs) at low redshift ($z\sim 0.3$). Detailed Gaussian modeling of optical emission line profiles in high-dispersion spectra consistently shows that both bright recombination and collisionally excited lines can be fitted as one or two narrow components with intrinsic velocity dispersion of $\sigma \sim$40-100 km s$^{-1}$, on top of a broader component with $\sigma \sim$100-300 km s$^{-1}$, which contributes up to $\sim$40% of the total flux and is preferentially blue-shifted from the systemic velocity. We interpret the narrow emission as highly ionized gas close to the young massive star clusters and the broader emission as a signpost of unresolved ionized outflows, resulting from stellar winds and supernova feedback. We find a significant correlation between the width of the broad emission and the LyC escape fraction, with strong LCEs exhibiting more complex and broader line profiles than galaxies with weaker or undetected LyC emission. Our findings provide new observational evidence supporting predictions from models and simulations, which suggest that gas turbulence and outflows resulting from strong radiative and mechanical feedback play a key role in clearing channels through which LyC photons escape from galaxies. We propose that the detection of blue-shifted broad emission in the nebular lines of compact extreme emission-line galaxies can provide a new indirect diagnostic of Lyman photon escape, which could be useful to identify potential LyC leakers in the epoch of reionization with the JWST.

Flares are intense explosions on the solar and stellar surfaces, and solar flares are sometimes accompanied by filament or prominence eruptions. Recently, a large filament eruption associated with a superflare on a solar-type star EK Dra was discovered for the first time. The absorption of the H$\alpha$ spectrum initially exhibited a blueshift with the velocity of $510$ (km s$^{-1}$), and decelerated in time probably due to gravity. Stellar coronal mass ejections (CMEs) were thought to occur, although the filament eruption did not exceed the escape velocity under the surface gravity. To investigate how such filament eruption occur and whether CMEs are associated with the filament eruption or not, we perform one-dimensional hydrodynamic simulation of the flow along an expanding magnetic loop emulating a filament eruption under adiabatic and unsteady conditions. The loop configuration and expanding velocity normal to the loop are specified in the configuration parameters, and we calculate the line-of-sight velocity of the filament eruption using the velocities along and normal to the loop. We found that (i) the temporal variations of the H$\alpha$ spectrum for EK Dra can be explained by falling filament eruption in the loop with longer time and larger spatial scales than that of the Sun, and (ii) the stellar CMEs are also thought to be associated with the filament eruption from the superflare on EK Dra, because the rarefied loop unobserved in the H$\alpha$ spectrum needs to expand faster than the escape velocity, whereas the observed filament eruption does not exceed the escape velocity.

We present results from the ``Quasar hosts Unveiled by high Angular Resolution Techniques" (QUART) survey studying the Circumgalactic Medium (CGM) by observing rest-frame UV emission lines Ly$\alpha$, C IV and He II around two radio-loud quasars, 3C 9 (z=2.02) and 4C 05.84 (z=2.32), using Keck Cosmic Web Imager (KCWI). We detect large-scale Ly$\alpha$ nebulae around both quasars with projected diameters $\sim$ 100 kpc, with spatially resolved, embedded 15-30 kpc He II and C IV nebulae around both quasars as well as kinematically distinct He II and C IV nebulae at a physical separation of $\sim$ 15 kpc from both quasars. Observations of H$\alpha$, H$\beta$, and [O III] emission using Keck MOSFIRE spectroscopically confirm that the Ly$\alpha$ nebulae extend to companion galaxies and that these quasars are in a protogroup/protocluster environment. We confirm that the He II and C IV emission is kinematically and spatially coincident with the companion galaxies. We estimate the virial masses of the companion galaxies, their metallicities, and star formation rates, and investigate the sources of ionization. We measure the dynamical mass of the host dark matter halos and estimate that the dark matter halos of these systems will grow to a mass of 2 $\times 10^{14}$ M$_{\odot}$ (3C 9) and 2 $\times 10^{13}$ M$_{\odot}$ (4C 05.84) by z=0. The combined CGM and companion galaxies observations indicate Ly$\alpha$ substructure can indicate the presence of companion galaxies in the CGM.

We aim to understand the influence of evolutionary modified surface abundances on the strength of the stellar wind in the helium-dominated subdwarfs CD-46 8926 and CD-51 11879. We used our own optical spectroscopy combined with archival ultraviolet spectroscopy and photometry to derive basic parameters and surface abundances of selected stars. The resulting parameterst served as input for the METUJE stellar wind code, which predicts the wind structure of these stars. We compared the derived wind parameters with the predictions derived for solar abundances. The optical analysis showed that both subdwarfs have effective temperatures in excess of 60 kK and a strong overabundance of carbon in the case of CD-46 8926 and nitrogen in the case of CD-51 11879. We interpret the abundance patterns as being a result of enrichment by the products of nuclear reactions. The modified abundances reduce the wind mass-loss rate by tens of percent. The reduction improves the predicted wind line profiles in comparison to observations. The change in helium abundance does not have a strong effect on the wind parameters. As a result of a lower estimated bolometric luminosity and mass-loss rate and a larger distance, the expected X-ray luminosities become lower and agree with observationas. The nucleosynthesis does not significantly alter the strength of the wind of hot subdwarfs, but the inclusion of proper stellar parameters improves the agreement with observational wind characteristics. Our analysis indicates that subdwarfs overabundant in helium are typically able to launch wind. This conclusion is supported by data gathered for thousands of subdwarfs from the literature, which shows that subdwarfs overabundant in helium avoid the region in the Kiel diagram where the winds are predicted to be absent. This can be interpreted in terms of the gravitational settling of helium, which is suppressed by the winds.

The nature of dark energy is one of the fundamental problems in cosmology. Introduced to explain the apparent acceleration of the Universe's expansion, its origin remains to be determined. In this paper, we illustrate a result that may significantly impact understanding the relationship between dark energy and structure formation in the late-epoch Universe. Our analysis exploits a scale-dependent energy functional, initially developed for image visualization, to compare the physical and geometrical data that distinct cosmological observers register on their celestial spheres. In the presence of late-epoch gravitational structures, this functional provides a non-perturbative technique that allows the standard Friedmann-Lema\^itre-Robertson-Walker (FLRW) observer to evaluate a measurable, scale-dependent difference between the idealized FLRW past light cone and the physical light cone. From the point of view of the FLRW observer, this difference manifests itself as a redshift-dependent correction $\Lambda^{(corr)}(z)$ to the FLRW cosmological constant $\Lambda^{(FLRW)}$. At the scale where cosmological expansion couples with the local virialized dynamics of gravitational structures, we get $\Lambda^{(corr)}(z)\sim 10^{-52}\,m^{-2}$, indicating that the late-epoch structures induce an effective cosmological constant that is of the same order of magnitude as the assumed value of the FLRW cosmological constant, a result that may lead to an interpretative shift in the very role of dark energy.

Kuiper belt objects show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of $^{26}$Al. We aim to explain the density trend in the context of streaming instability and pebble accretion. Small pebbles experience lofting into the atmosphere of the disk, being exposed to UV and partially losing their ice via desorption. Conversely, larger pebbles are shielded and remain more icy. We use a shearing box model including gas and solids, the latter split into ices and silicate pebbles. Self-gravity is included, allowing dense clumps to collapse into planetesimals. We find that the streaming instability leads to the formation of mostly icy planetesimals, albeit with an unexpected trend that the lighter ones are more silicate-rich than the heavier ones. We feed the resulting planetesimals into a pebble accretion integrator with a continuous size distribution, finding that they undergo drastic changes in composition as they preferentially accrete silicate pebbles. The density and masses of large KBOs are best reproduced if they form between 15 and 22\,AU. Our solution avoids the timing problem because the first planetesimals are primarily icy, and $^{26}$Al is mostly incorporated in the slow phase of silicate pebble accretion. Our results lend further credibility to the streaming instability and pebble accretion as formation and growth mechanisms.

This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified.

Only a handful of massive starless core candidates have been discovered so far, but none of them have been fully confirmed. Within the MM1 clump in the filamentary infrared dark cloud G34.43+0.24 that was covered by the ALMA-ATOMS survey at Band 3 ($\sim2\arcsec$, 6000\,au) and the ALMA-QUARKS survey at Band 6 ($\sim 0.3\arcsec$, 900\,au), two prestellar core candidates MM1-C and E1 with masses of 71 and 20 \solarmass~and radii of 2100--4400\,au were discovered. The two cores show no obvious sign of star-formation activities. In particular, MM1-C is a very promising massive prestellar core candidate with a total gas mass of 71\,\solarmass. Within MM1-C, we detected two extremely dense substructures, C1 and C2, as characterized by their high densities of $\rm n_{H_2}\sim 10^{8-9} cm^{-3}$. Moreover, evidence of further fragmentation in C2 was also revealed. We have detected the primordial fragmentation in the earliest stage of massive star formation, and we speculate that MM1-C would be the birthplace of a massive multiple system. However, we cannot fully rule out the possibility that the massive prestellar core MM1-C will just form a cluster of low-mass stars if it undergoes further fragmentation.

We present a demonstration version of a commensal pipeline for Fast Radio Burst (FRB) searches using a real-time incoherent beam from the Murchison Widefield Array (MWA). The main science target of the pipeline are bright nearby FRBs from the local Universe which are the best candidates to probe FRB progenitors and understand physical mechanisms powering these extremely energetic events. The new MWA beamformer, known as the "MWAX multibeam beamformer", can form multiple incoherent and coherent beams commensally to any on-going MWA observations. One of the beams is currently used for FRB searches (tested in 10 kHz frequency resolution and time resolutions between 0.1 and 100 ms). A second beam is used for a Search for Extraterrestrial Intelligence (SETI). This paper focuses on the FRB search pipeline and its verification on selected known bright pulsars. The pipeline uses the FREDDA implementation of the Fast Dispersion Measure Transform algorithm (FDMT) for single pulse searches. Initially, it was tested during standard MWA observations, and more recently using dedicated observations of selected 11 bright pulsars. The pulsar PSR J0835-4510 (aka Vela) has been routinely used as the primary probe of the data quality because its folded profile was always detected in the frequency band 200 - 230 MHz with typical SNR >10. Similarly, the low DM pulsar PSR B0950+08 was always detected in folded profile in the frequency band 140 - 170 MHz, and so far has been the only object for which single pulses were detected. We present the estimated sensitivity of the search in the currently limited observing bandwidth of a single MWA coarse channel (1.28 MHz) and for the upgraded, future system with 12.8 MHz (10 channels) of bandwidth. Based on expected sensitivity and existing FRB rate measurements, we estimate an expected number of FRB detections to be between a few and a few tens per year.

We have been developing silicon-on-insulator (SOI) pixel detectors with a pinned depleted diode (PDD) structure, named "XRPIX", for X-ray astronomy. The PDD structure is formed in a thick p-type substrate, to which high negative voltage is applied to make it fully depleted. A pinned p-well is introduced at the backside of the insulator layer to reduce a dark current generation at the Si-SiO$_{2}$ interface and to fix the back-gate voltage of the SOI transistors. An n-well is further introduced between the p-well and the substrate to make a potential barrier between them and suppress a leakage current. An optimization study on the n-well dopant concentration is necessary because a higher dopant concentration could result in a higher potential barrier but also in a larger sense-node capacitance leading to a lower spectroscopic performance, and vice versa. Based on a device simulation, we fabricated five candidate chips having different n-well dopant concentrations. We successfully found out the best n-well design, which suppressed a large leakage current and showed satisfactory X-ray spectroscopic performance. Too low and too high n-well dopant concentration chips showed a large leakage current and degraded X-ray spectroscopic performance, respectively. We also found that the dependency of X-ray spectroscopic performance on the n-well dopant concentration can be largely explained by the difference in sense-node capacitance.

The radius valley (or gap) in the observed distribution of exoplanet radii, which separates smaller super-Earths from larger sub-Neptunes, is a key feature that theoretical models must explain. Conventionally, it is interpreted as the result of the loss of primordial H/He envelopes atop rocky cores. However, planet formation models predict that water-rich planets migrate from regions outside the snowline toward the star. Assuming water to be in the form of solid ice in their interior, many of these planets would be located in the radius gap, in disagreement with observations. Here we use an advanced coupled formation and evolution model that describes the planets' origin and evolution starting from moon-sized, planetary seed embryos in the protoplanetary disk to mature Gyr-old planetary systems. Employing new equations of state and interior structure models to treat water as vapor mixed with H/He, we naturally reproduce the valley at the observed location. The model results indicate that the valley separates less massive, in-situ, rocky super-Earths from more massive, ex-situ, water-rich sub-Neptunes. Furthermore, the occurrence drop at larger radii, the so-called radius cliff, is also matched by planets with water-dominated envelopes. Owing to our statistical approach, we can assess that the synthetic distribution of radii quantitatively agrees with observations for the close-in population of planets; but only if atmospheric photoevaporation is also acting, populating the super-Earth peak with evaporated rocky cores. Therefore, we provide a hybrid theoretical explanation of the radius gap and cliff caused by both formation (orbital migration) as well as evolution (atmospheric escape).

Millisecond pulsars (MSPs) are a kind of radio pulsars with short spin periods, playing a key role in many aspects of stellar astrophysics. In recent years, some more MSPs with wide orbits ($>30\,\rm d$) have been discovered, but their origin is still highly unclear. In the present work, according to an adiabatic power-law assumption for the mass-transfer process, we carried out a large number of complete binary evolution computations for the formation of MSPs with wide orbits through the iron core-collapse supernova (CCSN) channel, in which a neutron star (NS) originating from a CCSN accretes matter from a red-giant (RG) star and spun up to millisecond periods. We found that this channel can form the observed MSPs with wide orbits in the range of $30-1200\,{\rm d}$, in which the WD companions have masses in the range of $0.28-0.55\,\rm M_{\odot}$. We also found that almost all the observed MSPs can be reproduced by this channel in the WD companion mass versus orbital period diagram. We estimate that the Galactic numbers of the resulting MSPs from the CCSN channel are in the range of $\sim 0.9-1.4\times10^{6}$. Compared with the accretion-induced collapse channel, the CCSN channel provides a dominant way to produce MSPs with wide orbits.

We identified a new ultra-distant comet C/2019 E3 (ATLAS) exhibiting preperihelion cometary activity at heliocentric distances $\gtrsim\!20$ au, making it the fourth member of this population after C/2010 U3 (Boattini), C/2014 UN$_{271}$ (Bernardinelli-Bernstein), and C/2017 K2 (PANSTARRS). From serendipitous archival data, we conducted analyses of the comet, finding that the activity was consistent with steady-state behaviour, suggestive of sublimation of supervolatiles, that the cross-section of dust increased gradually on the inbound leg of the orbit, varying with heliocentric distances as $r_{\rm H}^{-1.5 \pm 0.4}$, and that the dust was produced at a rate of $\gtrsim\!10^2$ kg s$^{-1}$ within the observed timespan. Our modelling of the largely symmetric morphology of the comet suggests that the dust environment was likely dominated by mm-scale dust grains ejected at speeds $\lesssim\!0.4$ m s$^{-1}$ from the sunlit hemisphere of the nucleus. Assuming a typical geometric albedo of 0.05 and adopting several simplistic thermophysical models, we estimated the nucleus to be at least $\sim\!3$ km across. We also measured the colour of the comet to be consistent with other long-period comets, except being slightly bluer in $g-r$. With our astrometric measurements, we determined an improved orbit of the comet, based upon which we derived that the comet is dynamically new and that its perihelion distance will further shrink due to the Galactic tide. We conclude the paper by comparing the known characteristics of the known ultra-distant comets.

In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30m antenna, we mapped N$_2$H$^+$(1-0), N$_2$D$^+$(1-0), H$^{13}$CO$^+$ (1-0) and (2-1), DCO$^+$(2-1), and $p$-NH$_2$D(1$_{11}$-1$_{01}$) towards the chosen cores. The missing $p$-NH$_3$ and N$_2$H$^+$(1-0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the GAS survey data (ammonia), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO$^+$/H$^{13}$CO$^+$ is one order of magnitude lower ($\sim$0.08) than that of the core central parts traced by the nitrogen-bearing species ($\sim$0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of $R_D$ with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for $R_D$(N$_2$D$^+$/N$_2$H$^+$) and R$_D$(DCO$^+$/HCO$^+$) and underestimate the $R_D$(NH$_2$D/NH$_3$).

A meteor outburst consisting of at least 22 meteors above the Baltic sea and southern Scandinavia that occurred on 30 October 2022 was recorded using multiple cameras. A bright fireball was followed by fainter meteors over a 10 second period. All the meteors were travelling on parallel trajectories. The goal of this study is to determine the atmospheric trajectories and photometric masses of meteors and to use these data to determine the specifics of the progenitor meteoroid break-up and cluster formation. Double and triple-station observations using video cameras were used for the calculation of the atmospheric trajectories and photometric masses of the meteors. Their relative positions and mass distribution were then used to determine the time and cause of the meteoroid fragmentation. The relative position of the cluster particles in the atmosphere and the distribution of their masses best correspond to the separation of the smaller fragments from the mass-dominant fragment 10.6$\pm$1.7 days before the collision with Earth, assuming a meteoroid bulk density of 1000 kg.m$^{-3}$. The ejection velocities are in the range 0.16-0.61 m.s$^{-1}$. The directions of the ejection velocities are bounded by a cone with an apex angle of $43^\circ$. The axis of this cone has ecliptic coordinates of $l=154^\circ$ and $b=26^\circ$ and is $66^\circ$ away from the direction to the Sun. Thermal stresses appear to be the most likely cause of such meteor cluster formation.

Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray selected BL Lac-type blazar 1ES 1959+650 in 2022 May 3-4 showed a significant linear polarization degree of $\Pi_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $\psi_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $\Pi_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $\Pi_\mathrm{O} \sim 5\%$ is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day.

A fundamental assumption of the data model used in radio interferometric calibration is that the sky model only consists of compact and discrete radio sources. This assumption breaks down when there are large scale diffuse structure such as the Galaxy visible in the observed data. No straightforward method currently exists to include such large scale diffuse structure in calibration and only indirect techniques such as excluding short baselines or filtering are used in practice as a remedy. In this paper, we propose a novel mechanism to include large scale diffuse sky models into direction dependent calibration of radio interferometers. We extend distributed calibration of radio interferometric data with both spectral and spatial regularization to include models for diffuse emission. We use shapelet basis functions to model both the diffuse sky structure as well as the spatial variation of systematic errors. The application of the direction dependent errors onto the diffuse sky is done in closed form using specific properties of shapelet basis functions, avoiding the need for expensive operations such as convolution. We provide extensive simulations showing the efficacy of our proposed method. We are able to overcome a major problem faced by existing calibration techniques, i.e., the suppression of large scale diffuse structure by not properly modeling such structure in calibration.

Indirect detection of gamma rays with ground-based observatories is currently the most sensitive experimental approach to characterize the gamma-ray sky at energies $>0.1$\,TeV. Ground-based detection of gamma-rays relies on the electromagnetic showers that gamma rays initiate in the Earth's atmosphere. In this chapter we will review the properties of electromagnetic air showers as well as the differences with respect to cosmic-ray showers that enable the rejection of the cosmic ray background. The experimental techniques that have been developed for ground-based detection of gamma rays will be introduced. These fall onto three main categories: air shower particle detectors, sampling Cherenkov arrays, and imaging atmospheric Cherenkov telescopes. Hybrid concepts as well as other experimental approaches are also discussed.

Solar radio bursts are some of the brightest emissions at radio frequencies in the solar system. The emission mechanisms that generate these bursts offer a remote insight into physical processes in solar coronal plasma, while fine spectral features hint at its underlying turbulent nature. During radio noise storms many hundreds of solar radio bursts can occur over the course of a few hours. Identifying and classifying solar radio bursts is often done manually although a number of automatic algorithms have been produced for this purpose. The use of machine learning algorithms for image segmentation and classification is well established and has shown promising results in the case of identifying Type II and Type III solar radio bursts. Here we present the results of a convolutional neural network applied to dynamic spectra of NenuFAR solar observations. We highlight some initial success in segmenting radio bursts from the background spectra and outline the steps necessary for burst classification.

Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z=1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (\lambda_rest = 1.3mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_\odot, with star formation rates of 25-160 M_\odot yr^-1 for the host galaxies. These QSOs have lower $L'_\mathrm{CO}$ than star-forming galaxies with the same L_IR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90 Myr. We find molecular gas masses in the range 2-40 x 10^9(alpha_CO/4) M_\odot, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z~1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population.

Methods. We performed measurements of the \MgII, \MgI, \FeIIa, \FeIIb, and \FeIIc\ spectral lines present in the spectra of the selected sample to determine the EW and velocity of the flows observed in the star-forming galaxies. Subsequently, we conducted $10^7$ bootstrap simulations using Spearman's rank correlation coefficient ($\rho_s$) to explore correlations with galaxy properties. Furthermore, we calculated the covering factor, gas density, and optical depth for the measured \ion{Fe}{II} doublets. Results. Our analysis revealed strong correlations between the EW of \ion{Mg}{II} lines and both $M_{*}$ ($\rho_s=0.43$, 4.5$\sigma$) and SFR ($\rho_s=0.42$, 4.4$\sigma$). For the \ion{Fe}{II} lines, we observed strong correlations between the EW and SFR ($\rho_s\sim0.65$, $>3.9\sigma$), with a weaker correlation for $M_{*}$ ($\rho_s\sim0.35$, $>1.9\sigma$). No notable correlations were found between velocity measurements of \ion{Mg}{II} line and $M_{*}$, SFR, or sSFR of the objects ($\rho_s\sim0.1)$. However, a negative strong correlation was found between the velocity of the \ion{Fe}{II} lines and the SFR of the galaxies ($\rho_s\sim-0.45$, $\sim3\sigma$). Our results align with previous studies but studying FIR-selected objects. Finally, we detected a candidate \textit{loitering outflow}, a recently discovered subtype of FeLoBAL quasar, at redshift of $z=1.4399$, exhibiting emission in \ion{C}{III}] and low line velocities ($|v|\lesssim$ 200 km/s).

In this study, we investigate the dynamics of the Universe during the observed late-time acceleration phase within the framework of the Weyl-type $f(Q,T)$ theory. Specifically, we consider a well-motivated model with the functional form $f(Q,T)=\alpha Q+\frac{\beta }{6\kappa ^2}T$, where $Q$ represents the scalar of non-metricity and $T$ denotes the trace of the energy-momentum tensor. In this context, the non-metricity $Q_{\mu\alpha\beta}$ of the space-time is established by the vector field $w_\mu$. The parameters $\alpha$ and $\beta$ govern the gravitational field and its interaction with the matter content of the Universe. By considering the case of dust matter, we obtain exact solutions for the field equations and observe that the Hubble parameter $H(z)$ follows a power-law behavior with respect to redshift $z$. To constrain the model parameters, we analyze various datasets including the $Hubble$, $Pantheon$ datasets, and their combination. Our results indicate that the Weyl-type $f(Q,T)$ theory offers a viable alternative to explain the observed late-time acceleration of the Universe avoiding the use of dark energy.

Detections of transiting planets from the upcoming PLATO mission are expected to face significant contamination from contaminating eclipsing binaries, resulting in false positives. To counter this, a ground-based programme to acquire time-critical photometry is pursued. Its principal aim is to obtain time-series observations of the planet candidate and its surrounding stars at the times of expected transits. This programme is part of the PLATO Ground-based Observations Programme, which also covers spectroscopic and imaging observations. The current photometric follow-up programme is assembling the required observational resources, executing benchmark observations, and defining strategies for the observations and their reporting. Post-launch, it will focus on coordinating photometric data collection and analysis, and will update candidate statuses in the PLATO follow-up database. Its work packages are outlined, covering specific tools, citizen contributions, standard and multi-colour observations, secondary eclipses, and reprocessing of archival photometry. Ground-based follow-up photometry will likely concentrate on longer-period candidates, given that false positives of short-period candidates will likely become identifiable in timeseries available from GAIA in the near future. Geographical considerations for follow-up observations from the first PLATO long-observation field LOPS2 are outlined, which lies in the southern hemisphere, with later fields expected to be more suitable for northern observers.

Parker's mean-field model includes two processes generating large-scale oscillatory dynamo waves: stretching of magnetic field lines by small-scale helical flows, and by differential rotation. In this work, we investigate the capacity of data-driven modal analysis, Dynamic Mode Decomposition, to identify coherent magnetic field structures of this model. In its canonical form, the only existing field scale corresponds to the dynamo instability. To take into account multi-scale nature of the dynamo, the model was augmented with coherent in time flow field, forcing small-scale magnetic field with a faster temporal evolution. Two clusters of DMD modes were obtained: the ``slow" cluster, located near the dynamo wave frequency and associated with its nonlinear self-interaction, and the ``fast" cluster, centered around the forcing frequency and resulting from the interaction between the wave and the flow. Compared to other widely used methods of data analysis, such as Fourier transform, DMD provides a natural spatiotemporal basis for the dynamo, related to its nonlinear dynamics. We assess how the parameters of the DMD model, rank and delay, influence its accuracy, and finally discuss the limitations of this approach when applied to randomly forced, more complex dynamo flows.

Dust rings in protoplanetary discs are often observed in thermal dust emission and could be favourable environments for planet formation. While dust rings readily form in gas pressure maxima, their long-term stability is key to both their observability and potential to assist in planet formation. We investigate the stability of the dust ring generated by interactions of a protoplanetary disc with a Neptune sized planet and consider its possible long term evolution using the FARGO3D Multifluid code. We look at the onset of the Rossby Wave Instability (RWI) and compare how the addition of dust in a disc can alter the stability of the gas phase. We find that with the addition of dust, the rings generated by planet disc interactions are more prone to RWI and can cause the gas phase to become unstable. The instability is shown to occur more easily for higher Stokes number dust, as it accumulates into a more narrow ring which triggers the RWI, while the initial dust fraction plays a more minor role in the stability properties. We show that the dusty RWI generates vortices that collect dust in their cores, which could be sites for further planetesimal formation. We conclude that the addition of dust can cause a ring in a protoplanetary disc to become more prone to instability leading to a different long term evolution compared to gas only simulations of the RWI.

A significant impediment to solving the coronal heating problem is that we currently only observe active region (AR) loops in their cooling phase. Previous studies showed that the evolution of cooling loop densities and apex temperatures are insensitive to the magnitude, duration, and location of energy deposition. Still, potential clues to how energy is released are encoded in the cooling phase properties. The appearance of coronal rain, one of the most spectacular phenomena of the cooling phase, occurs when plasma has cooled below 1MK, which sets constraints on the heating frequency, for example. Most observations of coronal rain have been made by imaging instruments. Here we report rare Hinode/EUV Imaging Spectrometer (EIS) observations of a loop arcade where coronal rain forms following an X2.1 limb flare. A bifurcation in plasma composition measurements between photospheric at 1.5MK and coronal at 3.5MK suggests that we are observing post-flare driven coronal rain. Increases in non-thermal velocities and densities with decreasing temperature (2.7MK to 0.6MK) suggest that we are observing the formation and subsequent evolution of the condensations. Doppler velocity measurements imply that a 10% correction of apparent flows in imaging data is reasonable. Emission measure analysis at 0.7MK shows narrow temperature distributions, indicating coherent behaviour reminiscent of that observed in coronal loops. The space-time resolution limitations of EIS suggest that we are observing the largest features or rain showers. These observations provide insights into the heating rate, source, turbulence, and collective behaviour of coronal rain from observations of the loop cooling phase.

We measured the dipole of the diffuse $\gamma$-ray background (DGB) identifying a highly significant time-independent signal coincidental with that of the Pierre Auger UHECR. The DGB dipole is determined from flux maps in narrow energy bands constructed from 13 years of observations by the Large Area Telescope (LAT) of the {\it Fermi} satellite. The $\gamma$-ray maps were clipped iteratively of sources and foregrounds similar to that done for the cosmic infrared background. The clipped narrow energy band maps were then assembled into one broad energy map out to the given energy starting at $E=2.74$ Gev, where the LAT beam falls below the sky's pixel resolution. Next we consider cuts in Galactic latitude and longitude to probe residual foreground contaminations from the Galactic Plane and Center. In the broad energy range $2.74 < E\leq115.5$ GeV the measured dipoles are stable with respect to the various Galactic cuts, consistent with an extragalactic origin. The $\gamma$-ray sky's dipole/monopole ratio is much greater than that expected from the DGB clustering component and the Compton-Getting effect origin with reasonable velocities. At $\simeq (6.5-7)\%$ it is similar to the Pierre Auger UHECRs with $E_{\rm UHECR}\ge 8$ EeV pointing to a common origin of the two dipoles. However, the DGB flux associated with the found DGB dipole reaches parity with that of the UHECR around $E_{\rm UHECR}\le 1$ EeV, perhaps arguing for a non-cascading mechanism if the DGB dipole were to come from the higher energy UHECRs. The signal/noise of the DGB dipole is largest in the $5-30$ GeV range, possibly suggesting the $\gamma$-photons at these energies are the ones related to cosmic rays.

Excess of gamma-rays with a spherical morphology around the Galactic center (GC) observed in the Fermi Large Area Telescope (LAT) data is one of the most intriguing features in the gamma-ray sky. The spherical morphology and the spectral energy distribution with a peak around a few GeV are consistent with emission from annihilation of dark matter particles. Other possible explanations include a distribution of millisecond pulsars (MSPs). One of the caveats of the MSP hypothesis is the relatively small number of associated MSPs near the GC. In this paper we perform a multi-class classification of Fermi-LAT sources using machine learning and determine the contribution from MSP-like sources among unassociated Fermi-LAT sources near the GC. We find that the unassociated MSP-like sources have a spectral energy distribution that is comparable with the GC excess if the contribution from the Fermi bubbles is taken into account. The spatial morphology and the source count distribution are also consistent with expectations for a population of MSPs that can explain the gamma-ray excess. Possible caveats of the contribution from the unassociated MSP-like sources include uncertainties in the diffuse emission model that can affect the number and the properties of the point-like sources detected near the GC and the contribution of the Fermi bubbles.

We apply the vertical Jeans equation to the Milky Way disk, in order to study non-axisymmetric variations in the thin disk surface density. We divide the disk plane into area cells with a 100 pc grid spacing, and use four separate subsets of the Gaia DR3 stars, defined by cuts in absolute magnitude, reaching distances up to 3 kpc. The vertical Jeans equation is informed by the stellar number density field and the vertical velocity field; for the former, we use maps produced via Gaussian Process regression; for the latter, we use Bayesian Neural Network radial velocity predictions, allowing us to utilize the full power of the Gaia DR3 proper motion sample. For the first time, we find evidence of a spiral arm in the form of an over-density in the dynamically measured disk surface density, detected in all four data samples, which also agrees very well with the spiral arm as traced by stellar age and chemistry. We fit a simple spiral arm model to this feature, and infer a relative over-density of roughly 20 % and a width of roughly 400 pc. We also infer a thin disk surface density scale length of 3.3--4.2 kpc, when restricting the analysis to stars within a distance of 2 kpc.

The Andromeda Galaxy is home to the annually erupting recurrent nova (RN) M31N 2008-12a (12a); the first nova found to host a nova super-remnant (NSR). A NSR is an immense structure surrounding a RN, created from many millions of eruptions sweeping up material in the local environment to form a shell tens of parsecs across. Theory has demonstrated that NSRs should be found around all RNe, even those systems with long periods between eruptions. Befittingly, the second NSR was found around the Galactic classical (and long suspected recurrent) nova, KT Eridani. In this Paper, we aim to find more of these phenomena through conducting the first ever survey for NSRs in M31 and the Large Magellanic Cloud (LMC). We find that the surroundings of fourteen RNe in M31 as well as the surroundings of the four RNe in the LMC do not show any evidence of vast parsec-scale structures in narrowband (H${\alpha}$ and [S II]) images, unlike the one clearly seen around 12a, and therefore conclude that observable NSRs are either rare structures, or they are too faint (or small) to be detected in our existing datasets. Yet, the NSR surrounding 12a would also likely to have been overlooked in our study if it were approximately one magnitude fainter. Searches for NSRs around other RNe 'masquerading' as classical novae may prove to be fruitful as would whole surveys of other Local Group galaxies.

We analyse electron acceleration by a large-scale electric field $E$ in a collisional hydrogen plasma under the solar flare coronal conditions based on approaches proposed by Dreicer and Spitzer for the dynamic friction force of electrons. The Dreicer electric field $E_{Dr}$ is determined as a critical electric field at which the entire electron population runs away. Two regimes of strong ($E \lesssim E_{Dr}$) and weak ($E \ll E_{Dr}$) electric field are discussed. It is shown that the commonly used formal definition of the Dreicer field leads to an overestimation of its value by about five times. The critical velocity at which the electrons of the ``tail'' of the Maxwell distribution become runaway under the action of the sub-Dreiser electric fields turns out to be underestimated by $\sqrt{3}$ times in some works because the Coulomb collisions between runaway and thermal electrons are not taken into account. The electron acceleration by sub-Dreicer electric fields generated in the solar corona faces difficulties.

This work explores characteristics of the negative polarization branch (NPB), which occurs in scattered light from rough surfaces, with particular focus on the effects of fine particles. Factors such as albedo, compression, roughness, and the refractive index are considered to determine their influence on the NPB. This study compiles experimental data and lunar observations to derive insights from a wide array of literature. Employing our proposed methodology, we estimate the representative grain sizes on the lunar surface to be $D \sim 1 \mathrm{-} 2 \mathrm{\mu m}$, with $D \lesssim 2 \mathrm{-} 4 \mathrm{\mu m}$, consistent with observed grain size frequency distributions in laboratory settings for lunar fines. Considering Mars, we propose that the finest particles are likely lacking ($D\gg 10 \mathrm{\mu m}$), which matches previous estimations. This study highlights the potential of multiwavelength, particularly near-infrared, polarimetry for precisely gauging small particles on airless celestial bodies. The conclusions provided here extend to cross-validation with grain sizes derived from thermal modeling, asteroid taxonomic classification, and regolith evolution studies.

Policy Brief on "AstroInformatics, Recommendations for Global Collaboration", distilled from panel discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. The deliberations encompassed a wide array of topics, including broad astroinformatics, sky surveys, large-scale international initiatives, global data repositories, space-related data, regional and international collaborative efforts, as well as workforce development within the field. These discussions comprehensively addressed the current status, notable achievements, and the manifold challenges that the field of astroinformatics currently confronts. The G20 nations present a unique opportunity due to their abundant human and technological capabilities, coupled with their widespread geographical representation. Leveraging these strengths, significant strides can be made in various domains. These include, but are not limited to, the advancement of STEM education and workforce development, the promotion of equitable resource utilization, and contributions to fields such as Earth Science and Climate Science. We present a concise overview, followed by specific recommendations that pertain to both ground-based and space data initiatives. Our team remains readily available to furnish further elaboration on any of these proposals as required. Furthermore, we anticipate further engagement during the upcoming G20 presidencies in Brazil (2024) and South Africa (2025) to ensure the continued discussion and realization of these objectives. The policy webinar took place during the G20 presidency in India (2023). Notes based on the seven panels will be separately published.

From 31 Earth-based and three Voyager 2 occultations spanning 1977--2006, we determine the orbital elements of the nine main Uranian rings with typical RMS residuals of 0.2 -- 0.4 km and 1-$\sigma$ errors in $a, ae,$ and $a\sin i$ of order 0.1 km, registered on an absolute radius scale accurate to 0.2 km at the 2-$\sigma$ level. The $\lambda$ ring shows more substantial scatter. In addition to the free modes $m=0$ in the $\gamma$ ring and $m=2$ in the $\delta$ ring, we find two additional outer Lindblad resonance (OLR) modes ($m=-1$ and $-2$) and a possible $m=3$ inner Lindblad resonance (ILR) mode in the $\gamma$ ring. No normal modes are detected for rings 6, 5, 4, $\alpha$, or $\beta$. Five normal modes are forced by small moonlets: the 3:2 inner ILR of Cressida with the $\eta$ ring, the 6:5 ILR of Ophelia with the $\gamma$ ring, the 23:22 ILR of Cordelia with the $\delta$ ring, the 14:13 ILR of Ophelia with the outer edge of the $\epsilon$ ring, and the counterpart 25:24 OLR of Cordelia with the ring's inner edge. We determine the width-radius relations for nearly all of the detected mode. We find no convincing evidence for librations of any of the rings. The Uranus pole direction at epoch TDB 1986 Jan 19 12:00 is $\alpha_P=77.311327\pm 0.000141^\circ$ and $\delta_P=15.172795\pm0.000618^\circ$. We determine the zonal gravitational coefficients $J_2=(3509.291\pm0.412)\times 10^{-6}, J_4=(-35.522\pm0.466)\times10^{-6}$, and $J_6$ fixed at $0.5\times 10^{-6}$, with a correlation coefficient $\rho(J_2,J_4)=0.9861$, for a reference radius $R=$25559 km. From the amplitudes and resonance radii of normal modes forced by moonlets, we determine the masses of Cressida, Cordelia, and Ophelia. Their estimated densities decrease systematically with increasing orbital radius and generally follow the radial trend of the Roche critical density for a shape parameter $\gamma=1.6$.

Discrepancies between cosmological parameter estimates from cosmic shear surveys and from recent Planck cosmic microwave background measurements challenge the ability of the highly successful $\Lambda$CDM model to describe the nature of the Universe. To rule out systematic biases in cosmic shear survey analyses, accurate redshift calibration within tomographic bins is key. In this paper, we improve photo-$z$ calibration via Bayesian hierarchical modeling of full galaxy photo-$z$ conditional densities, by employing $\textit{StratLearn}$, a recently developed statistical methodology, which accounts for systematic differences in the distribution of the spectroscopic training/source set and the photometric target set. Using realistic simulations that were designed to resemble the KiDS+VIKING-450 dataset, we show that $\textit{StratLearn}$-estimated conditional densities improve the galaxy tomographic bin assignment, and that our $\textit{StratLearn}$-Bayesian framework leads to nearly unbiased estimates of the target population means. This leads to a factor of $\sim 2$ improvement upon the previously best photo-$z$ calibration method. Our approach delivers a maximum bias per tomographic bin of $\Delta \langle z \rangle = 0.0095 \pm 0.0089$, with an average absolute bias of $0.0052 \pm 0.0067$ across the five tomographic bins.

We select 1,052,469 (754,635) thin disk stars from {\it Gaia} eDR3 and LAMOST DR7 in the range of Galactocentric radius $R$ (guiding center radius $R_\mathrm{g}$) from 8 to 11\,kpc to investigate the asymmetries between the North and South of the disk midplane. More specifically we analyze the vertical velocity dispersion profiles ($\sigma_{v_{z}}(z$)) in different bins of $R$ ($R_\mathrm{g}$) and $[\mathrm{Fe/H}]$. We find troughs in the profiles of $\sigma_{v_{z}}(z)$ located in both the North ($z \sim 0.7$\,kpc) and South ($z \sim -0.5$\,kpc) of the disk at all radial and chemical bins studied. The difference between the Northern and Southern vertical velocity dispersion profiles ($\Delta\sigma_{v_{z}}(|z|)$) shows a shift between curves of different $R$ and $R_\mathrm{g}$. A similar shift exists in these NS asymmetry profiles further divided into different $[\mathrm{Fe/H}]$ ranges. The sample binned with $R_\mathrm{g}$ more clearly displays the features in the velocity dispersion profiles. The shift in the peaks of the $\Delta\sigma_{v_{z}}$ profiles and the variation in the phase spiral shape binned by metallicity indicate the variation of the vertical potential profiles and the radial metallicity gradient. The wave-like signal in NS asymmetry of $\sigma_{v_{z}}(z)$ largely originates from phase spiral; while the NS asymmetry profiles of [Fe/H] only display a weak wave-like feature near solar radius. We perform a test particle simulation to qualitatively reproduce the observed results. A quantitative explanation of the NS asymmetry in the metallicity profile needs careful consideration of the spiral shape and the perturbation model, and we leave this for future work.

We propose and perform a joint analysis of the two different mass estimates of galaxy clusters, namely the hydrostatic and caustic techniques. Firstly, we show comprehensively that the mass bias between these two techniques can be possibly alleviated when cluster-specific assumptions constrained using the hydrostatic technique are utilized within the caustic technique. While at face value this demotes the caustic technique from a completely independent method, this allows one to further tighten the constraints on the cluster mass and subsequently, allow us to test modifications to gravity. Implementing the aforementioned formalism for two well-observed massive galaxy clusters, A2029 and A2142, we highlight the proof of concept. In the current implementation, we use this method to constrain the Chameleon screening and Vainshtein screening. As anticipated, we show that the joint analysis can help improve the constraints on these modified gravity scenarios.

Long-standing unexplained Venus atmosphere observations and chemical anomalies point to unknown chemistry but also leave room for the possibility of life. The unexplained observations include several gases out of thermodynamic equilibrium (e.g. tens of ppm O2, the possible presence of PH3 and NH3, SO2 and H2O vertical abundance profiles), an unknown composition of large, lower cloud particles, and the "unknown absorber(s)". Here we first review relevant properties of the Venus atmosphere and then describe the atmospheric chemical anomalies and how they motivate future astrobiology missions to Venus.

The chemical enrichment of X-ray-emitting hot atmospheres has hitherto been primarily studied in galaxy clusters. These studies revealed relative abundances of heavy elements that are remarkably similar to Solar. Here, we present measurements of the metal content of M89 (NGC 4552), an elliptical galaxy infalling into the Virgo cluster with a $\sim$10 kpc ram-pressure stripped X-ray tail. We take advantage of deep Chandra and XMM-Newton observations, and with particular attention to carefully modelling the spectra, we measure the O/Fe, Ne/Fe, Mg/Fe, Si/Fe and S/Fe ratios. Contrary to previous measurements in galaxy clusters, our results for the hot atmosphere of M89 suggest super-Solar abundance ratios with respect to iron (i.e. $\alpha$/Fe > 1), similar to its stellar components. Our analysis of the active galactic nucleus (AGN) activity in this system indicates that the AGN-induced outflow could have facilitated the stripping of the original galactic atmosphere, which has been replaced with fresh stellar mass loss material with super-Solar $\alpha$/Fe abundance ratios. Additionally, we report a new fitting bias in the RGS data of low-temperature plasma. The measured O/Fe ratios are >1$\sigma$ lower in multi-temperature models than a single temperature fit, leading to discrepancies in the calculations of supernova fractions derived from the metal abundances.

Low-frequency radio observations show an increasing number of radio galaxies located in galaxy clusters that display peculiar morphologies and spectral profiles. This is the result of the dynamical interaction of the galaxy with the surrounding medium. Studying this phenomenon is key to understanding the evolution of low-energy relativistic particles in the intracluster medium. We present a multi-frequency study of the three head-tail (HT) radio galaxies and the radio halo in the galaxy cluster ZwCl0634.1+4747. We make use of observations at four frequencies performed with LOFAR LBA (53 MHz), HBA (144 MHz), GMRT (323 MHz) and VLA (1518 MHz) data. The use of extremely low radio frequency observations, such as LOFAR at 53 and 144 MHz, allowed us to detect the extension of the tails up to a distance of ~ 1 Mpc. We extracted spectral profiles along the tails in order to identify possible departures from a pure ageing model, such as the Jaffe-Perola (JP) model, which only involves synchrotron and inverse-Compton losses. We found clear evidence of departures from this simple ageing model, such as surface brightness enhancement and spectral flattening along all of the tails. This can be interpreted as the consequence of particle re-acceleration along the tails. Possible explanations for this behaviour include the interaction between a shock and the radio tails or a turbulence-driven re-acceleration mechanism. We show that the latter scenario is able to reproduce the characteristic features that we observed in our profiles.

We present a novel approach to the numerical computation of quasi-normal modes, based on the first-order (in radial derivative) formulation of the equations of motion and using a matrix version of the continued fraction method. This numerical method is particularly suited to the study of static black holes in modified gravity, where the traditional second-order, Schr\"odinger-like, form of the equations of motion is not always available. Our approach relies on the knowledge of the asymptotic behaviours of the perturbations near the black hole horizon and at spatial infinity, which can be obtained via the systematic algorithm that we have proposed recently. In this work, we first present our method for the perturbations of a Schwarzschild black hole and show that we recover the well-know frequencies of the QNMs to a very high precision. We then apply our method to the axial perturbations of an exact black hole solution in a particular scalar-tensor theory of gravity. We also cross-check the obtained QNM frequencies with other numerical methods.

Parity violation in the gravitational sector is a prediction of many theories beyond general relativity. In the propagation of gravitational waves, parity violation manifests by inducing amplitude and/or velocity birefringence between right- and left-circularly polarized modes. We study how the stochastic gravitational wave background can be used to place constraints on these birefringent effects. We consider two model scenarios, one in which we allow birefringent corrections to become arbitrarily large, and a second in which we impose stringent theory priors. In the former, we place constraints on a generic birefringent gravitational-wave signal due to the current non-detection of a stochastic background from compact binary events. We find a joint constraint on birefringent parameters, $\kappa_D$ and $\kappa_z$, of $\mathcal{O}(10^{-1})$. In the latter scenario, we forecast constraints on parity violating theories resulting from observations of the future upgraded LIGO-Virgo-KAGRA network as well as proposed third-generation detectors. We find that third-generation detectors will be able to improve the constraints by at least two orders of magnitude, yielding new stringent bounds on parity violating theories. This work introduces a novel and powerful probe of gravitational parity violation with gravitational-wave data.

We present an analytic frequency-domain gravitational waveform model for an inspiraling binary whose center-of-mass undergoes a small acceleration, assumed to be constant during the detection, such as when it orbits a distant tertiary mass. The center-of-mass acceleration along the line of sight is incorporated as a new parameter that perturbs the standard TaylorF2 model. We calculate the wave phase to 3rd post-Newtonian order and first order in the acceleration. It is shown that acceleration most significantly modifies the wave phase in the low frequency portion of the signal, so ground-based detectors with a good sensitivity at low frequencies are the most effective at detecting this effect. We present a Fisher information calculation to quantify the detectability of this effect at advanced LIGO A Plus, Cosmic Explorer, and Einstein Telescope over the mass range of neutron stars and stellar-mass black holes, and discuss degeneracy between acceleration and other parameters. We also determine the parameter space where the acceleration is large enough that the wave phase model would have to be extended to nonlinear orders in the acceleration.

We present our brief response and some comments to the article [Phys. Rev. C 108. 065801 (2023)]. In [Phys. Rev. C 108. 065801 (2023)] contains at least one unfounded and false statement, and our answer illuminates this explicitly. We made a graphical comparison of our results [Phys. Rev. Phys. Rev. C 105, 065806 (2022)] with the [Phys. Rev. C 108. 065801 (2023)]. Since in [Phys. Rev. C 108. 065801 (2023)] our potentials from [Phys. Rev. C 105, 065806 (2022)] with the slight modification were used, the results of these works are almost identical but in [Phys. Rev. C 108. 065801 (2023)] the opposite was stated.

With about a hundred binary black hole (BBH) mergers detected by LIGO-Virgo-KAGRA, and with several hundreds expected in the current O4 run, GWs are revolutionizing our understanding of the universe. Some BBH sources are too faint to be individually detected, but collectively they may give rise to a stochastic GW background (SGWB). In this paper, we calculate the SGWB associated with BBH mergers dynamically assembled in dense star clusters, using state-of-the-art numerical models. We discuss the role of modeling the evolution of the mass distribution of BBH mergers, which has significant implications for model selection and parameter estimation, and could be used to distinguish between different channels of BBH formation. We demonstrate how the birth properties of star clusters affect the amplitude and frequency spectrum of the SGWB, and show that upcoming observation runs of ground-based GW detectors may be sensitive enough to detect it. Even in the case of a non-detection, we find that GW data can be used to constrain the highly uncertain cluster birth properties, which can complement direct observations of young massive clusters and proto-star clusters in the early universe by JWST.

Gravitational waves are a unique probe of the early Universe, as the Universe is transparent to gravitational radiation right back to the beginning. In this article, we summarise detection prospects and the wide scope of primordial events that could lead to a detectable stochastic gravitational wave background. Any such background would shed light on what (if anything) lies beyond the Standard Model, sometimes at remarkably high scales. We overview the range of strategies for detecting a stochastic gravitational wave background before delving deep into three major primordial events that can source such a background. Finally, we summarize the landscape of other sources of primordial backgrounds.

Recent advances in far-infrared detector technology have led to increases in raw sensitivity of more than an order of magnitude over previous state-of-the-art detectors. With such sensitivity, photon noise becomes the dominant noise component, even when using cryogenically cooled optics, unless a method of restricting the spectral bandpass is employed. The leading instrument concept features reflecting diffraction gratings which post-disperse the light that has been modulated by a polarizing Fourier transform spectrometer (FTS) onto a detector array, thereby reducing the photon noise on each detector. This paper discusses the development of a cryogenic (4 K) diffraction grating spectrometer which operates over the wavelength range from \SIrange{285}{500}{\micro \meter} and was used to post-disperse the output from a room-temperature polarizing FTS. Measurements of the grating spectral response and diffraction efficiency are presented as a function of both wavelength and polarization to characterize the instrumental performance.

We discuss a novel mechanism for the proto-neutron star acceleration assisted by the chiral separation effect which induces an axial vector current in a dense medium. We consider the process of neutrinos scattering off the background axial vector current of electrons. We show that anisotropy of either magnetic field or density in momentum space is essential for nonzero recoil and we call this mechanism the chiral anisotropy conversion. Assuming a strong magnetic field $B \simeq 10^{12}$ T and anisotropy by $\sim 10\%$, we find that the chiral anisotropy conversion can yield the velocity of order of typical pulsar kicks, i.e., $v_{\mathrm{kick}} \gtrsim 1000$ km/s.