Papers for Wednesday, Dec 13 2023

Earth-mass planets are expected to have atmospheres and experience thermal tides raised by the host star. These tides transfer energy to the planet that can counter the dissipation from bodily tides. Indeed, even a relatively thin atmosphere can drive the rotation of these planets away from the synchronous state. Here we revisit the dynamical evolution of planets undergoing thermal atmospheric tides. We use a novel approach based on a vectorial formalism, which is frame independent and valid for any configuration of the system, including any eccentricity and obliquity values. We provide the secular equations of motion after averaging over the mean anomaly and the argument of the pericenter, which are suitable to model the long-term spin and orbital evolution of the planet.

By June of 2020, the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spectrum Roentgen Gamma observatory had completed its first of the planned eight X-ray all-sky survey (eRASS1). The large effective area of the X-ray telescope makes it ideal for a survey of the faint X-ray diffuse emission over half of the sky with an unprecedented energy resolution and position accuracy. In this work, we produce the X-ray diffuse emission maps of the eRASS1 data with a current calibration, covering the energy range from 0.2 to 8.0 keV. We validated these maps by comparison with X-ray background maps derived from the ROSAT All Sky Survey (RASS). We generated X-ray images with a pixel area of 9 arcmin$^2$ using the observations available to the German eROSITA consortium. The contribution of the particle background to the photons was subtracted from the final maps. We also subtracted all the point sources above a flux threshold dependent on the goal of the subtraction, exploiting the eRASS1 catalog that will soon be available. The accuracy of the eRASS1 maps is shown by a flux match to the RASS X-ray maps, obtained by converting the eROSITA rates into equivalent ROSAT count rates in the standard ROSAT energy bands R4, R5, R6, and R7, within 1.25$\sigma$. We find small residual deviations in the R4, R5, and R6 bands, where eROSITA tends to observe lower flux than ROSAT (~11%), while a better agreement is achieved in the R7 band (~1%). The eRASS maps exhibit lower noise levels than RASS maps at the same resolution above 0.3 keV. We report the average surface brightness and total flux of different large sky regions as a reference. The detection of faint emission from diffuse hot gas in the Milky Way is corroborated by the consistency of the eRASS1 and RASS maps shown in this paper and by their comparable flux dynamic range.

The red damping wing from neutral hydrogen in the intergalactic medium is a smoking-gun signal of ongoing reionization. One potential contaminant of the intergalactic damping wing signal is dense gas associated with foreground galaxies, which can give rise to proximate damped Ly$\alpha$ absorbers. The Ly$\alpha$ imprint of such absorbers on background quasars is indistinguishable from the intergalactic medium within the uncertainty of the intrinsic quasar continuum, and their abundance at $z\gtrsim7$ is unknown. Here we show that the complex of low-ionization metal absorption systems recently discovered by deep JWST/NIRSpec observations in the foreground of the $z=7.54$ quasar ULAS~J1342$+$0928 can potentially reproduce the quasar's spectral profile close to rest-frame Ly$\alpha$ without invoking a substantial contribution from the intergalactic medium, but only if the absorbing gas is extremely metal-poor ($[{\rm O}/{\rm H}]\sim-3.5$). Such a low oxygen abundance has never been observed in a damped Ly$\alpha$ absorber at any redshift, but this possibility still complicates the interpretation of the spectrum. Our analysis highlights the need for deep spectroscopy of high-redshift quasars with JWST or ELT to "purify" damping wing quasar samples, an exercise which is impossible for much fainter objects like galaxies.

We simulate the formation of second-generation stars in young clusters with masses of $10^{5}$ and $10^6\ {\rm M_\odot}$ within $30-100\ {\rm Myr}$ after the formation of clusters. We assume the clusters move through a uniform interstellar medium with gas densities of $10^{-24}$ and $10^{-23}\ {\rm gcm}^{-3} $ and consider the stellar winds from asymptotic giant branch (AGB) stars, gas accretion onto the cluster, ram pressure, star formation, and photoionization feedback of our stellar systems including binary stars. We find that second-generation (SG) stars can be formed only within the $10^6\ {\rm M_\odot}$ cluster in the high-density simulation, where the cluster can accrete sufficient pristine gas from their surrounding medium, leading to efficient cooling required for the ignition of SG formation and sufficient dilution of the AGB ejecta. Hence, our results indicate that a denser environment is another requirement for the AGB scenario to explain the presence of multiple populations in globular clusters. On the other hand, the ionizing feedback becomes effective in heating the gas in our low-density simulations. As a result, the clusters cannot accumulate a considerable amount of pristine gas at their center. The gas mass within the clusters in these simulations is similar to that in young massive clusters (YMCs). Hence, our studies can provide a possible reason for the lack of gas, star formation, and SG stars in YMCs. Our results indicate that the ionizing stellar feedback is not a severe problem for SG formation; rather, it can help the AGB scenario to account for some observables.

Numerous sky background subtraction techniques have been developed since the first implementations of computer-based reduction of spectra. Kurtz & Mink (2000) described a singular value decomposition-based method which allowed them to subtract night sky background from multi-fiber spectroscopic observations without any additional sky observations. We hereby take this approach one step further with usage of non-negative matrix factorization instead of principal component analysis and generalize it to 2D spectra. This allows us to generate approximately 10 times as many valid eigenspectra because of non-negativity. We apply our method to short-slit spectra of low-mass galaxies originating from intermediate-resolution Echelle spectrographs (ESI at Keck, MagE at Magellan, X-Shooter at the VLT) when sources fill the entire slit. We demonstrate its efficiency even when no offset sky observations were collected.

With two central galaxies engaged in a major merger and a remarkable chain of 19 young stellar superclusters wound around them in projection, the galaxy cluster SDSS J1531+3414 ($z=0.335$) offers an excellent laboratory to study the interplay between mergers, AGN feedback, and star formation. New Chandra X-ray imaging reveals rapidly cooling hot ($T\sim 10^6$ K) intracluster gas, with two "wings" forming a concave density discontinuity near the edge of the cool core. LOFAR $144$ MHz observations uncover diffuse radio emission strikingly aligned with the "wings," suggesting that the "wings" are actually the opening to a giant X-ray supercavity. The steep radio emission is likely an ancient relic of one of the most energetic AGN outbursts observed, with $4pV > 10^{61}$ erg. To the north of the supercavity, GMOS detects warm ($T\sim 10^4$ K) ionized gas that enshrouds the stellar superclusters but is redshifted up to $+ 800$ km s$^{-1}$ with respect to the southern central galaxy. ALMA detects a similarly redshifted $\sim 10^{10}$ M$_\odot$ reservoir of cold ($T\sim 10^2$ K) molecular gas, but it is offset from the young stars by $\sim 1{-}3$ kpc. We propose that the multiphase gas originated from low-entropy gas entrained by the X-ray supercavity, attribute the offset between the young stars and the molecular gas to turbulent intracluster gas motions, and suggest that tidal interactions stimulated the "beads on a string" star formation morphology.

The small-scale structure of baryons in the intergalactic medium is intimately linked to their past thermal history. Prior to the $\gtrsim10^4$ K photoheating during the epoch of reionization, cold baryons may have closely traced the clumpy cosmic web of dark matter down to scales as low as $\lesssim1$ comoving kpc, depending on the degree of heating by the X-ray background. After the passage of the ionization front, this clumpy structure can persist for $\sim10^{8}$ years. The strong Ly$\alpha$ damping wings detected towards a few of the highest redshift quasars, in addition to their smaller-than-expected Ly$\alpha$-transmissive proximity zones, suggest that they have ionized and heated the foreground intergalactic medium less than $10^7$ years ago. Signatures of the pre-reionization small-scale structure should thus persist in their intergalactic surroundings. Here we explore how the persistence of this clumpy structure can affect the statistics of Ly$\alpha$ transmission inside the transparent proximity zones of $z\gtrsim7$ quasars by post-processing a suite of small-volume hydrodynamical simulations with 1D ionizing radiative transfer. We find that the Ly$\alpha$ flux power spectrum and flux PDF statistics of ten $z=7.5$ proximity zones, with realistic observational parameters, could distinguish the gaseous structure of a $T_{\rm IGM}\sim2$ K CDM model from warm dark matter models with particle masses $m_{\rm WDM}>10$ keV and X-ray heated models with $f_{\rm X}f_{\rm abs}>0.1$ ($T_{\rm IGM}(z=7.5)\gtrsim275$ K) at the $2\sigma$ level.

Halos with masses in excess of the atomic limit are believed to be ideal environments in which to form heavy black hole seeds with masses above 10^3 Msun. In cases where the H_2 fraction is suppressed this is expected to lead to reduced fragmentation of the gas and the generation of a top heavy initial mass function. In extreme cases this can result in the formation of massive black hole seeds. Resolving the initial fragmentation scale and the resulting protostellar masses has, until now, not been robustly tested. Cosmological simulations were performed with the moving mesh code Arepo using a primordial chemistry network until z = 11. Three haloes with masses in excess of the atomic cooling mass were then selected for detailed examination via zoom-ins. The highest resolution simulations resolve densities up to 10^-6 g cm^-3 (10^18 cm^-3) and capture a further 100 yr of fragmentation behaviour at the center of the halo. Our simulations show intense fragmentation in the central region of the halos, leading to a large number of near-solar mass protostars. Despite the increased fragmentation the halos produce a protostellar mass spectrum that peaks at higher masses relative to standard Population III star forming halos. The most massive protostars have accretion rates of 10^-3-10^-1 Msun yr^-1 after the first 100 years of evolution, while the total mass of the central region grows at 1 Msun yr^-1. Lower resolution zoom-ins show that the total mass of the system continues to accrete at 1 Msun yr^-1 for at least 10^4 yr, although how this mass is distributed amongst the rapidly growing number of protostars is unclear. However, assuming that a fraction of stars can continue to accrete rapidly the formation of a sub-population of stars with masses in excess of 10^3 Msun is likely in these halos.

Previously we demonstrated that the magnetorotational instability (MRI) grows vigorously in eccentric disks, much as it does in circular disks, and we investigated the nonlinear development of the eccentric MRI without vertical gravity. Here we explore how vertical gravity influences the magnetohydrodynamic (MHD) turbulence stirred by the eccentric MRI. Similar to eccentric disks without vertical gravity, the ratio of Maxwell stress to pressure, or the Shakura--Sunyaev alpha parameter, remains ~0.01, and the local sign flip in the Maxwell stress persists. Vertical gravity also introduces two new effects. Strong vertical compression near pericenter amplifies reconnection and dissipation, weakening the magnetic field. Angular momentum transport by MHD stresses broadens the mass distribution over eccentricity at much faster rates than without vertical gravity; as a result, spatial distributions of mass and eccentricity can be substantially modified in just ~5 to 10 orbits. MHD stresses in the eccentric debris of tidal disruption events may power emission $\gtrsim$1 yr after disruption.

Blue large-amplitude pulsators (BLAPs) make up a rare class of hot pulsating stars with effective temperatures of $\approx$30,000 K and surface gravities of 4.0 - 5.0 dex (cgs). The evolutionary origin and current status of BLAPs is not well understood, largely based on a lack of spectroscopic observations and no available mass constraints. However, several theoretical models have been proposed that reproduce their observed properties, including studies that identify them as pulsating helium-core pre-white dwarfs (He-core pre-WDs). We present here follow-up high-speed photometry and phase-resolved spectroscopy of one of the original 14 BLAPs, OGLE-BLAP-009, discovered during the Optical Gravitational Lensing Experiment. We aim to explore its pulsation characteristics and determine stellar properties such as mass and radius in order to test the consistency of these results with He-core pre-WD models. Using the mean atmospheric parameters found using spectroscopy, we fit a spectral energy distribution to obtain a preliminary estimate of the radius, luminosity and mass by making use of the Gaia parallax. We then compare the consistency of these results to He-core pre-WD models generated using MESA, with predicted pulsation periods implemented using GYRE. We find that our mass constraints are in agreement with a low-mass He-core pre-WD of $\approx$0.30 M$_{\odot}$.

Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here we compute the evolution of stars having such a BH at their center. We find that such objects can be surprisingly long-lived, with the lightest black holes having no influence over stellar evolution, while more massive ones consume the star over time to produce a range of observable consequences. Models of the Sun born about a BH whose mass has since grown to approximately $10^{-6}~\rm{M}_\odot$ are compatible with current observations. In this scenario, the Sun would first dim to half its current luminosity over a span of 100 Myr as the accretion starts to generate enough energy to quench nuclear reactions. The Sun would then expand into a fully-convective star, where it would shine luminously for potentially several Gyr with an enriched surface helium abundance, first as a sub-subgiant star, and later as a red straggler, before becoming a sub-solar-mass BH. We also present results for a range of stellar masses and metallicities. The unique internal structures of stars harboring BHs may make it possible for asteroseismology to discover them, should they exist. We conclude with a list of open problems and predictions.

Solar flares are known to leave imprints on the magnetic field at the photosphere, often manifested as an abrupt and permanent change in the downward-directed Lorentz force in localized areas inside the active region. Our study aims to differentiate eruptive and confined solar flares based on the vertical Lorentz force variations. We select 26 eruptive and 11 confined major solar flares (stronger than the GOES M5 class) observed during 2011-2017. We analyze these flaring regions using SHARP vector-magnetograms obtained from the NASA's Helioseismic and Magnetic Imager (HMI). We also compare data corresponding to 2 synthetic flares from a $\delta$--sunspot simulation reported in Chatterjee et al. [Phys. Rev. Lett. 116, 101101 (2016)]. We estimate the change in the horizontal magnetic field and the total Lorentz force integrated over an area around the polarity inversion line (PIL) that encompasses the location of the flare. Our results indicate a rapid increase of the horizontal magnetic field along the flaring PIL, accompanied by a significant change in the downward-directed Lorentz force in the same vicinity. Notably, we find that all the confined events under study exhibit a total change in Lorentz force of < $1.8 \times 10^{22}$ dyne. This threshold plays an important factor in effectively distinguishing eruptive and confined flares. Further, our analysis suggests that the change in total Lorentz force also depends on the reconnection height in the solar corona during the associated flare onset. The results provide significant implications for understanding the flare-related upward impulse transmission for the associated coronal mass ejection.

We present the evolution of Ly$\alpha$ emission derived from 53 galaxies at $z=6.6-13.2$ that are identified by multiple JWST/NIRSpec spectroscopy programs of ERO, ERS, GO, DDT, and GTO. These galaxies fall on the star-formation main sequence and are the typical star-forming galaxies with UV magnitudes of $-22.5\leq M_\mathrm{UV}\leq-17.0$. We find that 15 out of 53 galaxies show Ly$\alpha$ emission at the $>3\sigma$ levels, and obtain Ly$\alpha$ equivalent width (EW) measurements and stringent $3\sigma$ upper limits for the 15 and 38 galaxies, respectively. Confirming that Ly$\alpha$ velocity offsets and line widths of our galaxies are comparable with those of low-redshift Ly$\alpha$ emitters, we investigate the redshift evolution of the Ly$\alpha$ EW. We find that Ly$\alpha$ EWs statistically decrease towards high redshifts on the Ly$\alpha$ EW vs. $M_{\rm UV}$ plane for various probability distributions of the uncertainties. We then evaluate neutral hydrogen fractions $x_{\rm HI}$ with the Ly$\alpha$ EW redshift evolution and the cosmic reionization simulation results on the basis of a Bayesian inference framework, and obtain $x_{\rm HI}=0.59^{+0.15}_{-0.33}$, $0.81^{+0.09}_{-0.26}$, and $0.99^{+0.01}_{-0.05}$ at $z\sim7$, $8$, and $9-13$, respectively. These moderately large $x_{\rm HI}$ values are consistent with the Planck CMB optical depth measurement and previous $x_{\rm HI}$ constraints from galaxy and QSO Ly$\alpha$ damping wing absorptions, and strongly indicate a late reionization history. Such a late reionization history suggests that major sources of reionization would emerge late and be hosted by moderately massive halos in contrast with the widely-accepted picture of abundant low-mass objects for the sources of reionization.

Coronagraphs are highly sensitive to wavefront errors, with performance degrading rapidly in the presence of low-order aberrations. Correcting these aberrations at the coronagraphic focal plane is key to optimal performance. We present two new methods based on the sequential phase diversity approach of the "Fast and Furious" algorithm that can correct low-order aberrations through a coronagraph. The first, called "2 Fast 2 Furious," is an extension of Fast and Furious to all coronagraphs with even symmetry. The second, "Tokyo Drift," uses a deep learning approach and works with general coronagraphic systems, including those with complex phase masks. Both algorithms have 100% science uptime and require effectively no diversity frames or additional hardware beyond the deformable mirror and science camera, making them suitable for many high contrast imaging systems. We present theory, simulations, and preliminary lab results demonstrating their performance.

Jupiter-family Comets (JFCs) exhibit a wide range of activity levels and mass-loss over their orbits. We analyzed high-cadence observations of 42 active JFCs with the wide-field Asteroid Terrestrial-impact Last Alert System (ATLAS) survey in 2020-2021. We measured dust production rates of the JFCs using the Af\rho parameter and its variation as a function of heliocentric distance. There is a tendency for our JFC sample to exhibit a maximum Af\rho after perihelion, with 254P/McNaught and P/2020 WJ5 (Lemmon) having their maximum Af\rho over a year after perihelion. On average, the rate of change of activity post-perihelion was shallower than that pre-perihelion. We also estimated the mass maximum loss rate for 17 of the JFCs in our sample, finding 4P/Faye to be the most active. We present a subset of comets whose measured Af\rho have been interpolated and extrapolated to a common distance of 2 au pre-perihelion and post-perihelion. From these measurements we found no correlation of intrinsic activity with current perihelion distance. For three of the JFCs in our sample, 6P/d'Arrest, 156P/Russell-LINEAR and 254P/McNaught, there was no visible coma but a constant absolute magnitude which we attributed to a probable detection of the nucleus. We derived upper limits for the nuclear radii of \leq 2.1 +/- 0.3 km, \leq 2.0 +/- 0.2 km and \leq 4.0 +/- 0.8 km respectively. Finally, we found that 4P/Faye, 108P/Ciffreo, 132P/Helin-Roman-Alu 2, 141P/Machholz 2, and 398P/Boattini experienced outbursts between 2020 and 2022.

Non-linear force-free extrapolations are a common approach to estimate the 3D topology of coronal magnetic fields based on photospheric vector magnetograms. The force-free assumption is a valid approximation at coronal heights, but for the dense plasma conditions in the lower atmosphere, this assumption is not satisfied. In this study, we utilize multi-height magnetic field measurements in combination with physics-informed neural networks to advance solar magnetic field extrapolations. We include a flexible height-mapping, which allows us to account for the different formation heights of the observed magnetic field measurements. The comparison to analytical and simulated magnetic fields demonstrates that including chromospheric magnetic field measurements leads to a significant improvement of our magnetic field extrapolations. We also apply our method to chromospheric line-of-sight magnetograms, from the Vector Spectromagnetograph (VSM) on the Synoptic Optical Long-term Investigations of the Sun (SOLIS) observatory, in combination with photospheric vector magnetograms, from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO). The comparison to observations in extreme ultraviolet wavelengths shows that the additional chromospheric information leads to a better agreement with the observed coronal structures. In addition, our method intrinsically provides an estimate of the corrugation of the observed magnetograms. With this new approach, we make efficient use of multi-height magnetic field measurements and advance the realism of coronal magnetic field simulations.

The heaviest chemical elements are naturally produced by the rapid neutron-capture process (r-process) during neutron star mergers or supernovae. The r-process production of elements heavier than uranium (transuranic nuclei) is poorly understood and inaccessible to experiments, so must be extrapolated using nucleosynthesis models. We examine element abundances in a sample of stars that are enhanced in r-process elements. The abundances of elements Ru, Rh, Pd, and Ag (atomic numbers Z = 44 to 47, mass numbers A = 99 to 110) correlate with those of heavier elements (63 <= Z <= 78, A > 150). There is no correlation for neighboring elements (34 <= Z <= 42 and 48 <= Z <= 62). We interpret this as evidence that fission fragments of transuranic nuclei contribute to the abundances. Our results indicate that neutron-rich nuclei with mass numbers >260 are produced in r-process events.

We present a dynamical study of 39 active Centaurs and 17 high-perihelion (q$>$4.5 au) JFCs with a focus on investigating recent orbital changes as potential triggers for comet-like activity. We have identified a common feature in the recent dynamical histories of all active Centaurs and JFC in our sample that is not present in the history of the majority of inactive population members: a sharp decrease in semi-major axis and eccentricity occurring within the last several hundred years prior to observed activity. We define these rapid orbital changes as `a-jumps'. Our results indicate that these orbital reshaping events lead to shorter orbital periods and subsequently greater average per-orbit heating of Centaur nuclei. We suggest the a-jumps could therefore be a major trigger of cometary activity on Centaurs and JFCs. Our results further imply that analyses of the recent dynamical histories could be used to identify objects that are currently active or may become active soon, where we have identified three such Centaurs with recent a-jumps that should be considered high-priority targets for observational monitoring to search for activity.

We investigate two microflares of GOES classes A9 and C1 (after background subtraction) observed by STIX onboard Solar Orbiter with exceptionally strong nonthermal emission. We complement the hard X-ray imaging and spectral analysis by STIX with co-temporal observations in the (E)UV and visual range by AIA and HMI, in order to investigate what makes these microflares so sufficient in high-energy particle acceleration. We performed case studies of two microflares observed by STIX on October 11, 2021 and November 10, 2022 that showed unusually hard microflare X-ray spectra with power-law indices of the electron flux distributions delta = 2.98 and delta = 4.08 during their nonthermal peaks and photon energies up to 76 keV and 50 keV respectively. For both events under study, we found that one footpoint is located within a sunspot covering areas with mean magnetic flux densities in excess of 1500 G, suggesting that the hard electron spectra are caused by the strong magnetic fields in which the flare loops are rooted. In addition, we revisited the unusually hard RHESSI microflare initially published by Hannah et al. (2008b) and found that in this event also one flare kernel was located within a sunspot, which corroborates the result from the two hard STIX microflares under study. We conclude that the characteristics of the strong photospheric magnetic fields inside sunspot umbrae and penumbrae where the flare loops are rooted play an important role in the generation of exceptionally hard X-ray spectra in these microflares.

We present resolved images of the inner disk component around HD 141569 using the Magellan adaptive optics system with the Clio2 1 - 5 $\mu$m camera, offering a glimpse of a complex system thought to be in a short evolutionary phase between protoplanetary and debris disk stages. We use a reference star along with the KLIP algorithm for PSF subtraction to detect the disk inward to about 0.24" (~25 au assuming a distance of 111 pc) at high signal-to-noise ratios at $L'$ (3.8 $\mu$m), $Ls$ (3.3 $\mu$m), and narrowband $Ice$ (3.1 $\mu$m). We identify an arc or spiral arm structure at the southeast extremity, consistent with previous studies. We implement forward modeling with a simple disk model within the framework of an MCMC sampler to better constrain the geometrical attributes and photometry using our KLIP-reduced disk images. We then leverage these modeling results to facilitate a comparison of the measured brightness in each passband to find a reduction in scattered light from the disk in the $Ice$ filter, implying significant absorption due to water ice in the dust. Additionally, our best-fit disk models exhibit peak brightness in the southwestern, back-scattering region of the disk, which we suggest to be possible evidence of 3.3 $\mu$m PAH emission. However, we point out the need for additional observations with bluer filters and more complex modeling to confirm these hypotheses.

We investigate the viability of primordial black hole (PBH) formation in the Standard Model (SM) in a scenario that does not rely on specific inflationary features or any exotic physics such as phase transitions or non-minimal coupling to gravity. If the Brout-Englert-Higgs (BEH) field lies exactly at the transition between metastability and stability, its potential exhibits an inflexion point due to radiative corrections. The BEH can act like a stochastic curvaton field, leading to a non-Gaussian tail of large curvature fluctuations that later collapse into PBHs when they re-enter inside the horizon. This scenario would require a precise value of the top-quark mas to ensure the Higgs stability, which is disfavored but still consistent with the most recent measurements. However, we also find that large curvature fluctuations are also generated on cosmological scales that are inconsistent with cosmic microwave background (CMB) observations. We therefore conclude that the SM cannot have led to the formation of PBHs based on this mechanism. Nevertheless, a variation of the scenario based on the Palatini formulation of gravity may have provided the conditions to produce stellar-mass PBHs with an abundance comparable to dark matter, without producing too large curvature fluctuations on cosmological scales.

Using images collected with the WFC3 camera on board of the Hubble Space Telescope, we detect stellar clumps in continuum-subtracted $H\alpha$ and ultraviolet (F275W filter), such clumps are often embedded in larger regions (star-forming complexes) detected in the optical (F606W filter). We model the photometry of these objects using BAGPIPES to obtain their stellar population parameters. The median mass-weighted stellar ages are 27 Myr for $H\alpha$ clumps and 39 Myr for F275W clumps and star-forming complexes, the oldest stars in the complexes can be older than $\sim$300 Myr which indicates that star-formation is sustained for long periods of time. Stellar masses vary from 10$^{3.5}$ to 10$^{7.1}$ $M_\odot$, with star-forming complexes being more massive objects in the sample. Clumps and complexes found further away from the host galaxy are younger, less massive and less obscured by dust. We interpret these trends as due to the effect of ram-pressure in different phases of the interstellar medium. $H\alpha$ clumps form a well-defined sequence in the stellar mass--SFR plane with slope 0.73. Some F275W clumps and star-forming complexes follow the same sequence while others stray away from it and passively age. The difference in stellar age between a complex and its youngest embedded clump scales with the distance between the clump and the center of the complex, with the most displaced clumps being hosted by the most elongated complexes. This is consistent with a fireball-like morphology, where star-formation proceeds in a small portion of the complex while older stars are left behind producing a linear stellar population gradient. The stellar masses of star-forming complexes are consistent with the ones of globular clusters, but stellar mass surface densities are lower by 2 dex, and their properties are more consistent with the population of dwarf galaxies in clusters.

We probe the neutral circumgalactic medium (CGM) along the major axes of NGC891 and NGC4565 in 21-cm emission out to $\gtrsim 100$kpc using the Green Bank Telescope (GBT), extending our previous minor axes observations. We achieve an unprecedented $5\sigma$ sensitivity of $6.1\times 10^{16}$ cm$^{-2}$ per 20 km s$^{-1}$ velocity channel. We detect HI with diverse spectral shapes, velocity widths, and column densities. We compare our detections to the interferometric maps from the Westerbork Synthesis Radio Telescope (WSRT) obtained as part of the HALOGAS survey. At small impact parameters, $> 31-43\%$ of the emission detected by the GBT cannot be explained by emission seen in the WSRT maps, and it increases to $> 64-73\%$ at large impact parameters. This implies the presence of diffuse circumgalactic HI. The mass ratio between HI in the CGM and HI in the disk is an order of magnitude larger than previous estimates based on shallow GBT mapping. The diffuse HI along the major axes pointings is corotating with the HI disk. The velocity along the minor axes pointings is consistent with an inflow and/or fountain in NGC891 and an inflow/outflow in NGC4565. Including the circumgalactic HI, the depletion time and the accretion rate of NGC4565 are sufficient to sustain its star formation. In NGC891, most of the required accreting material is still missing.

We studied the linear and nonlinear evolution of the Vertical Shear Instability (VSI) in axisymmetric models of protoplanetary disks, focusing on the transport of angular momentum, the produced temperature perturbations, and the applicability of local stability conditions. We modeled the gas-dust mixture via high-resolution two-moment (M1) radiation-hydrodynamical simulations including stellar irradiation with frequency-dependent opacities. We found that, given sufficient depletion of small grains (with a dust-to-gas mass ratio of $10\%$ of our nominal value of $10^{-3}$ for $<0.25$ $\mu$m grains), the VSI can operate in surface disk layers while being inactive close to the midplane, resulting in a suppression of the VSI body modes. The VSI reduces the initial vertical shear in bands of approximately uniform specific angular momentum, whose formation is likely favored by the enforced axisymmetry. Similarities with Reynolds stresses and angular momentum distributions in 3D simulations suggest that the VSI-induced angular momentum mixing in the radial direction may be predominantly axisymmetric. The stability regions in our models are well explained by local stability criteria, while the employment of global criteria is still justifiable up to a few scale heights above the midplane, at least as long as VSI modes are radially optically thin. Turbulent heating produces only marginal temperature increases of at most $0.1\%$ and $0.01\%$ in the nominal and dust-depleted models, respectively, peaking at a few (approximately three) scale heights above the midplane. We conclude that it is unlikely that the VSI can, in general, lead to any significant temperature increase since that would either require it to efficiently operate in largely optically thick disk regions or to produce larger levels of turbulence than predicted by models of passive irradiated disks.

The vertical shear instability (VSI) is a hydrodynamical instability likely to produce turbulence in the dead zones of protoplanetary disks. Various aspects of this instability remain to be understood, including the disk regions where it can operate and the physical phenomena leading to its saturation. In this work, we studied the growth and evolution of secondary instabilities parasitic to the VSI, examining their relation with its saturation in axisymmetric radiation-hydrodynamical simulations of protoplanetary disks. We also constructed stability maps for our disk models, considering temperature stratifications enforced by stellar irradiation and radiative cooling and incorporating the effects of dust-gas collisions and molecular line emission. We found that the flow pattern produced by the interplay of the axisymmetric VSI modes and the baroclinic torque forms bands of nearly uniform specific angular momentum. In the high-shear regions in between these bands, the Kelvin-Helmholtz instability (KHI) is triggered. The significant transfer of kinetic energy to small-scale eddies produced by the KHI and possibly even the baroclinic acceleration of eddies limit the maximum energy of the VSI modes, likely leading to the saturation of the VSI. A third instability mechanism, consisting of an amplification of eddies by baroclinic torques, forms meridional vortices with Mach numbers up to $\sim 0.4$. Our stability analysis suggests that protoplanetary disks can be VSI-unstable in surface layers up to tens of au for reasonably high gas emissivities, even in regions where the midplane is stable. This picture is consistent with current observations of disks showing thin midplane millimeter-sized dust layers while appearing vertically extended in optical and near-infrared wavelengths.

Here we report an optical quasi-periodic oscillation (QPO) with a period of $\sim$134~day detected in g- and r-band light curves of a narrow-line Seyfert 1 galaxy TXS 1206+549 at redshift of 1.34 with the data from the Zwicky Transient Facility (ZTF) observations. After considering the trial factor, the significance levels in the two bands are 3.1 $\sigma$ and 2.6 $\sigma$, respectively. The QPO signal presents about 10 cycles ranging from 2018 March to 2021 December lasting $\sim$4 year. A near-sinusoidal profile also appears in the folded light curves by a phase-resolved analysis. Interestingly, in the simultaneous light curve as the time scale of ZTF observations, a potential periodic signal with similar period is detected in the o-band light curve from the Asteroid Terrestrial-Impact Last Alert System data, additionally a weak peak is also detected at the similar period in the gamma-ray light curve obtained from the \emph{Fermi} Gamma-ray Space Telescope data. Some potential origins of periodicities in active galactic nuclei are discussed for the QPO reported here.

Organic features lead to two distinct types of Class 0/I low-mass protostars: hot corino sources, and warm carbon-chain chemistry (WCCC) sources. Some observations suggest that the chemical variations between WCCC sources and hot corino sources are associated with local environments, as well as the luminosity of protostars. We conducted gas-grain chemical simulation in collapsing protostellar cores, and found that the fiducial model predicts abundant carbon-chain molecules and COMs, and reproduces WCCC and hot corino chemistry in the hybrid source L483. By changing values of some physical parameters, including the visual extinction of ambient clouds ($A_{\rm V}^{\rm amb}$), the cosmic-ray ionization rate ($\zeta$), the maximum temperature during the warm-up phase ($T_{\rm max}$), and the contraction timescale of protostars ($t_{\rm cont}$), we found that UV photons and cosmic rays can boost WCCC features by accelerating the dissociation of CO and CH$_4$ molecules. On the other hand, UV photons can weaken the hot corino chemistry by photodissociation reactions, while the dependence of hot corino chemistry on cosmic rays is relatively complex. The $T_{\rm max}$ does not affect WCCC features, while it can influence hot corino chemistry by changing the effective duration of two-body surface reactions for most COMs. The long $t_{\rm cont}$ can boost WCCC and hot corino chemistry, by prolonging the effective duration of WCCC reactions in the gas phase and surface formation reactions for COMs, respectively. Subsequently, we ran a model with different physical parameters to reproduce scarce COMs in prototypical WCCC sources. The scarcity of COMs in prototypical WCCC sources can be explained by insufficient dust temperature in the inner envelopes to activate hot corino chemistry. Meanwhile, the High $\zeta$ and the long $t_{\rm cont}$ favors the explanation for scarce COMs in these sources.

We report the discovery of a bow-shock pulsar wind nebula (PWN), named Potoroo, and the detection of a young pulsar J1638-4713 that powers the nebula. We present a radio continuum study of the PWN based on 20-cm observations obtained from the Australian Square Kilometre Array Pathfinder (ASKAP) and MeerKAT. PSR J1638-4713 was identified using Parkes radio telescope observations at frequencies above 3 GHz. The pulsar has the second-highest dispersion measure of all known radio pulsars (1553 pc/cm^3), a spin period of 65.74 ms and a spin-down luminosity of 6.1x10^36 erg/s. The PWN has a cometary morphology and one of the greatest projected lengths among all the observed pulsar radio tails, measuring over 21 pc for an assumed distance of 10 kpc. The remarkably long tail and atypically steep radio spectral index are attributed to the interplay of a supernova reverse shock and the PWN. The originating supernova remnant is not known so far. We estimated the pulsar kick velocity to be in the range of 1000-2000 km/s for ages between 23 and 10 kyr. The X-ray counterpart found in Chandra data, CXOU J163802.6-471358, shows the same tail morphology as the radio source but is shorter by a factor of 10. The peak of the X-ray emission is offset from the peak of the radio total intensity (Stokes I) emission by approximately 4.7", but coincides well with circularly polarised (Stokes V) emission. No infrared counterpart was found.

We show that massive young star clusters may be possible candidates that can accelerate Galactic cosmic rays (CRs) in the range of $10^7\hbox{--}10^9$ GeV (between the `knee' and `ankle'). Various plausible scenarios such as acceleration at the wind termination shock (WTS), supernova shocks inside these young star clusters, etc. have been proposed,since it is difficult to accelerate particles up to the $10^7\hbox{--}10^9$ GeV range in the standard paradigm of CR acceleration in supernova remnants. We consider a model for the production of different nuclei in CRs from massive stellar winds using the observed distribution of young star clusters in the Galactic plane. We present a detailed calculation of CR transport in the Galaxy, taking into account the effect of diffusion, interaction losses during propagation, and particle re-acceleration by old supernova remnants to determine the all-particle CR spectrum. Using the maximum energy estimate from the Hillas criterion, we argue that a young massive star cluster can accelerate protons up to a few tens of PeV. Upon comparison with the observed data, our model requires a CR source spectrum with an exponential cutoff of $5\times 10^7 Z$ GeV ($50\,Z$~PeV) from these clusters together with a cosmic-ray injection fraction of $\sim 5\%$ of the wind kinetic energy. We discuss the possibility of achieving these requirements in star clusters, and the associated uncertainties, in the context of considering star clusters as the natural accelerator of the `second component' of Galactic cosmic rays.

Several ongoing and upcoming large scale structure surveys aim to explore the nonlinear regime of structure formation with high precision. Making reliable cosmological inferences from these observations necessitates precise theoretical modeling of the mildly nonlinear regime. In this work we explore how the choice of nonlinear prescription would impact parameter estimation from cosmic shear measurements for a Euclid-like survey. Specifically, we employ two different nonlinear prescriptions of Haloft and the Effective Field Theory of the Large Scale Structure and compare their measurements for the three different cosmological scenarios of $\Lambda$CDM, $w$CDM and $(w_0,w_a)$CDM. We also investigate the impact of different nonlinear cutoff schemes on parameter estimation. We find that the predicted errors on most parameters shrink considerably as smaller scales are included in the analysis, with the amount depending on the nonlinear prescription and the cutoff scheme used. We use predictions from the Halofit model to analyze the mock data from DarkSky $N$-body simulations and quantify the parameter bias introduced in the measurements due to the choice of nonlinear prescription. We observe that $\sigma_8$ and $n_{\rm{s}}$ have the largest measurement bias induced by inaccuracies of the Halofit model.

Although thick disk is a structure prevalent in local disk galaxies and also present in our home Galaxy, its formation and evolution is still unclear. Whether the thick disk is born thick and/or gradually heated to be thick after formation is under debate. To disentangle these two scenarios, one effective approach is to inspect the thickness of young disk galaxies in the high redshift Universe. In this work we study the vertical structure of 191 edge-on galaxies spanning redshift from 0.2 to 5 using JWST NIRCAM imaging observations. For each galaxy, we retrieve the vertical surface brightness profile at 1 ${R_e}$ and fit a sech$^2$ function that has been convolved with the line spread function. The obtained scale height of galaxies at $z>1.5$ show no clear dependence on redshift, with a median value in remarkable agreement with that of the Milky Way's thick disk. This suggests that local thick disks are already thick when they were formed in the early times and secular heating is unlikely the main driver of thick disk formation. For galaxies at $z<1.5$, however, the disk scale height decreases systematically towards lower redshift, with low-redshift galaxies having comparable scale height with that of the Milky Way's thin disk. This cosmic evolution of disk thickness favors an upside-down formation scenario of galaxy disks.

The origin of apparently young $\alpha$-rich stars in the Galaxy is still a matter of debate in Galactic archaeology, whether they are genuinely young or might be products of binary evolution and merger/mass accretion. We aim to shed light on the nature of young $\alpha$-rich stars in the Milky Way by studying their distribution in the Galaxy thanks to an unprecedented sample of giant stars that cover different Galactic regions and have precise asteroseismic ages, chemical, and kinematic measurements. We analyze a new sample of $\sim$ 6000 stars with precise ages coming from asteroseismology. Our sample combines the global asteroseismic parameters measured from light curves obtained by the K2 mission with stellar parameters and chemical abundances obtained from APOGEE DR17 and GALAH DR3, then cross-matched with Gaia DR3. We define our sample of young $\alpha$-rich stars and study their chemical, kinematic, and age properties. We investigate young $\alpha$-rich stars in different parts of the Galaxy and we find that the fraction of young $\alpha$-rich stars remains constant with respect to the number of high-$\alpha$ stars at $\sim$ 10%. Furthermore, young $\alpha$-rich stars have kinematic and chemical properties similar to high-$\alpha$ stars, except for [C/N] ratios. This suggests that these stars are not genuinely young, but products of binary evolution and merger/mass accretion. Under that assumption, we find the fraction of these stars in the field to be similar to that found recently in clusters. This fact suggests that $\sim$ 10% of the low-$\alpha$ field stars could also have their ages underestimated by asteroseismology. This should be kept in mind when using asteroseismic ages to interpret results in Galactic archaeology.

The recent analysis on the central compact object in the HESS J1731-347 remnant suggests interestingly small values for its mass and radius. Such an observation favors soft nuclear models that may be challenged by the observation of massive compact stars. In contrast, the recent PREX-II experiment, concerning the neutron skin thickness of $^{208}$Pb, points towards stiff equations of state that favor larger compact star radii. In the present study, we aim to explore the compatibility between stiff hadronic equations of state (favored by PREX-II) and the HESS J1731-347 remnant in the context of hybrid stars. For the construction of hybrid equations of state we use three widely employed Skyrme models combined with the well-known vector MIT bag model. Furthermore we consider two different scenarios concerning the energy density of the bag. In the first case, that of a constant bag parameter, we find that the resulting hybrid equations of state are strongly disfavored by the observation of $\sim2 M_\odot$ pulsars. However, the introduction of a Gaussian density dependence yields results that are compatible with the conservative $2 M_\odot$ constraint. The utilization of recent data based on the observation of PSR J0030+0451, PSR J0952-0607 and GW190814 allows for the imposition of additional constraints on the relevant parameters and the stiffness of the two phases. Interestingly, we find that the derived hybrid equations of state do not satisfy the PSR J0030+0451 constraints in $1\sigma$ and only marginally agree with the $2\sigma$ estimations. In addition, we estimate that the observation of massive pulsars, like PSR~J0952-0607, in combination with the existence of HESS J1731-347, may require a strong phase transition below $\sim 1.7n_0$. Finally, we show that the supermassive compact object involved in GW190814 could potentially be explained as a rapidly rotating hybrid star.

SY Cha is a T Tauri star surrounded by a protoplanetary disk with a large cavity seen in the millimeter continuum but has the spectral energy distribution (SED) of a full disk. Here we report the first results from JWST-MIRI Medium Resolution Spectrometer (MRS) observations taken as part of the MIRI mid-INfrared Disk Survey (MINDS) GTO Program. The much improved resolution and sensitivity of MIRI-MRS compared to Spitzer enables a robust analysis of the previously detected H2O, CO, HCN, and CO2 emission as well as a marginal detection of C2H2. We also report the first robust detection of mid-infrared OH and ro-vibrational CO emission in this source. The derived molecular column densities reveal the inner disk of SY Cha to be rich in both oxygen and carbon bearing molecules. This is in contrast to PDS 70, another protoplanetary disk with a large cavity observed with JWST, which displays much weaker line emission. In the SY Cha disk, the continuum, and potentially the line, flux varies substantially between the new JWST observations and archival Spitzer observations, indicative of a highly dynamic inner disk.

The transit method is one of the most relevant exoplanet detection techniques, which consists of detecting periodic eclipses in the light curves of stars. This is not always easy due to the presence of noise in the light curves, which is induced, for example, by the response of a telescope to stellar flux. For this reason, we aimed to develop an artificial neural network model that is able to detect these transits in light curves obtained from different telescopes and surveys. We created artificial light curves with and without transits to try to mimic those expected for the extended mission of the Kepler telescope (K2) in order to train and validate a 1D convolutional neural network model, which was later tested, obtaining an accuracy of 99.02 % and an estimated error (loss function) of 0.03. These results, among others, helped to confirm that the 1D CNN is a good choice for working with non-phased-folded Mandel and Agol light curves with transits. It also reduces the number of light curves that have to be visually inspected to decide if they present transit-like signals and decreases the time needed for analyzing each (with respect to traditional analysis).

Radio observations of stars trace the plasma conditions and magnetic field properties of stellar magnetospheres and coronae. Depending on the plasma conditions at the emitter site, radio emission in the metre- and decimetre-wave bands is generated via different mechanisms such as gyrosynchrotron, electron cyclotron maser instability, and plasma radiation processes. The ongoing LOFAR Two-metre Sky Survey (LoTSS) and VLA Sky Survey (VLASS) are currently the most sensitive wide-field radio sky surveys ever conducted. Because these surveys are untargeted, they provide an opportunity to study the statistical properties of the radio-emitting stellar population in an unbiased manner. Here, we perform an untargeted search for stellar radio sources down to sub-mJy level using these radio surveys. We find that the population of radio-emitting stellar systems is mainly composed of two distinct categories: chromospherically active stellar (CAS) systems and M dwarfs. We also seek to identify signatures of a gradual transition within the M-dwarf population from chromospheric/coronal acceleration close to the stellar surface similar to that observed on the Sun, to magnetospheric acceleration occurring far from the stellar surface similar to that observed on Jupiter. We determine that radio detectability evolves with spectral type, and we identify a transition in radio detectability around spectral type M4, where stars become fully convective. Furthermore, we compare the radio detectability vs spectra type with X-ray and optical flare (observed by TESS) incidence statistics. We find that the radio efficiency of X-ray/optical flares, which is the fraction of flare energy channelled into radio-emitting charges, increases with spectral type. These results motivate us to conjecture that the emergence of large-scale magnetic fields in CAS systems and later M dwarfs leads to an increase in radio efficiency.

One of the most prominent sources for energetic particles in our solar system are huge eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs), which usually drive shocks that accelerate charged particles up to relativistic energies. In particular, energetic electron beams can generate radio bursts through the plasma emission mechanism, for example, type II and accompanying herringbone bursts. Here, we investigate the acceleration location, escape, and propagation directions of various electron beams in the solar corona and compare them to the arrival of electrons at spacecraft. To track energetic electron beams, we use a synthesis of remote and direct observations combined with coronal modelling. Remote observations include ground-based radio observations from the Nancay Radioheliograph (NRH) combined with space-based extreme-ultraviolet and white-light observations from the Solar Dynamics Observatory (SDO), the Solar Terrestrial Relations Observatory (STEREO) and Solar Orbiter (SolO). We also use direct observations of energetic electrons from the STEREO and Wind spacecraft. These observations are then combined with a three-dimensional (3D) representation of the electron acceleration locations that combined with results from magneto-hydrodynamic models of the solar corona is used to investigate the origin and link of electrons observed remotely at the Sun to in situ electrons. We observed a type II radio burst followed by herringbone bursts that show single-frequency movement through time in NRH images. The movement of the type II burst and herringbone radio sources seems to be influenced by the regions in the corona where the CME is more capable of driving a shock. We also found similar inferred injection times of near-relativistic electrons at spacecraft to the emission time of the type II and herringbone bursts.

The effect of accretion bursts on massive young stellar objects (MYSOs) represents a new research field in the study of young stars and their environment. We aim to investigate the impact of an accretion burst on massive disks with different types of envelopes and to study the effects of an accretion burst on the temperature structure and the chemistry of the disk. We focus on water and methanol as chemical species for this paper. The thermochemical code of PRODIMO (PROtoplanetary DIsk MOdel) is used to perform simulation of high mass protoplanetary disk models with different types of envelopes under the presence of an accretion burst. The models in question represent different evolutionary stages of protostellar objects. We calculate and show the chemical abundances in three phases of the simulation (pre-burst, burst, and post-burst). More heavily embedded disks show higher temperatures. The impact of the accretion burst is mainly characterized by the desorption of chemical species present in the disk and envelope from the dust grains to the gas phase. When the post-burst phase starts, the sublimated species freeze out again. The degree of sublimation depends strongly on the type of envelope the disk is embedded in. An accretion burst in more massive envelopes produces stronger desorption of the chemical species. However, our models show that the timescale for the chemistry to reach the pre-burst state is independent of the type of envelope. The study shows that the disk's temperature increases with a more massive envelope enclosing it. Thus, the chemistry of MYSOs in earlier stages of their evolution reacts stronger to an accretion burst than at later stages where the envelope has lost most of its mass or has been dissipated. The study of the impact of accretion bursts could also provide helpful theoretical context to the observation of methanol masers in massive disks.

Recently, the emergence of cosmological tension has raised doubts about the consistency of the $\Lambda$CDM model. In order to constrain the neutrino mass within a consistent cosmological framework, we investigate three massive neutrinos with normal hierarchy (NH) and inverted hierarchy (IH) in both the axion-like EDE (Axi-EDE) model and the AdS-EDE model. We use the joint datasets including cosmic microwave background (CMB) power spectrum from Planck 2018, Pantheon of type Ia supernova, baryon acoustic oscillation (BAO) and $H_0$ data from SH0ES. For the $\nu$Axi-EDE model, we obtain $\sum m_{\nu,\mathrm{NH}} < 0.152$ eV and $\sum m_{\nu,\mathrm{IH}} < 0.178$ eV, while for the $\nu$AdS-EDE model, we find $\sum m_{\nu,\mathrm{NH}} < 0.135$ eV and $\sum m_{\nu,\mathrm{IH}} < 0.167$ eV. Our results exhibit a preference for the normal hierarchy in both the $\nu$Axi-EDE model and the $\nu$AdS-EDE model.

In 2020, hydrogen chloride (HCl) in the gas phase was discovered in the atmosphere of Mars with the Atmospheric Chemistry Suite (ACS) onboard the Trace Gas Orbiter (TGO) mission (Korablev et al., 2021). Its volume mixing ratio (VMR) shows a seasonal increase of up to 5 ppbv during the perihelion season, followed by a sudden drop to undetectable levels, contradicting previous estimations of the HCl lifetime of several months. In the Earth's stratosphere, heterogeneous uptake of HCl onto water ice is known to be a major sink for this species. Modelling of associated chemistry involving heterogeneous reactions indicates that H2O ice becomes the most effective sink for HCl above 20 km with the characteristic time shorter than 12 hours. In this work, we use simultaneous measurements of water ice particles and HCl abundance obtained by the ACS instrument and show particular structures in the vertical profiles, forming detached layers of gas at the ice free altitudes ('ice-holes'). We demonstrate that the heterogeneous uptake of HCl onto water ice operates on Mars and is potentially a major mechanism regulating the HCl abundance in the atmosphere of Mars.

Axions, and axion-like particles, have come back into fashion in the last decades as a possible solution to the galactic-scale crisis suffered by the cold dark matter model. In the framework of the wave Dark Matter model, we have carried out a Jeans Analysis on eight dwarf spheroidal galaxies that are orbiting around the Milky Way, and we have constrained the boson mass. Differently to a previous analysis, we adopted an anisotropy parameter that varies with the distance from the centre of the galaxy to assess whether this assumption would help to resolve, or at least alleviate, the well-known tension with the value of the boson mass favoured by the cosmological analysis. Our results indicate that, differently to what happens in ultra-faint dwarf galaxies, such a tension cannot be lifted introducing a variable anisotropy parameter, leaving as a possible solution the existence of additional axion or axion-like particles with higher masses as naturally predicted in the Axiverse.

It is well known that stars move away from their birth location over time via radial migration. This dynamical process makes computing the correct chemical evolution, e.g., metallicity gradients, of galaxies very difficult. This dynamical process makes inferring the chemical evolution of observed galaxies from their measured abundance gradients very difficult. One way to account for radial migration is to infer stellar birth radii for individual stars. Many attempts to do so have been performed over the last years, but are limited to the Milky Way as computing the birth position of stars requires precise measurements of stellar metallicity and age for individual stars that cover large Galactic radii. Fortunately, recent and future surveys will provide numerous opportunities for inferring birth radii for external galaxies such as the Large Magellanic Cloud (LMC). In this paper, we investigate the possibility of doing so using the NIHAO cosmological zoom-in simulations. We find that it is theoretically possible to infer birth radii with a ~ 25% median uncertainty for individual stars in galaxies with i) orderliness of the orbits, $\langle v_\phi \rangle/\sigma_{v} >$ 2, ii) a dark matter halo mass greater or equal to approximately the LMC mass (~ 2 x 10$^{11} M_\odot$), and iii) after the average azimuthal velocity of the stellar disk reaches ~70% of its maximum. From our analysis, we conclude that it is possible and useful to infer birth radii for the LMC and other external galaxies that satisfy the above criteria.

In the homoeoidal expansion, a given ellipsoidally stratified density distribution, and its associated potential, are expanded in the (small) density flattening parameter $\eta$, and usually truncated at the linear order. The truncated density-potential pair obeys exactly the Poisson equation, and it can be interpreted as the first-order expansion of the original ellipsoidal density-potential pair, or as a new autonomous system. In the first interpretation, in the solutions of the Jeans equations the quadratic terms in $\eta$ must be discarded (``$\eta$-linear'' solutions), while in the second (``$\eta$-quadratic'') all terms are retained. In this work we study the importance of the quadratic terms by using the ellipsoidal Plummer model and the Perfect Ellipsoid, which allow for fully analytical $\eta$-quadratic solutions. These solutions are then compared with those obtained numerically for the original ellipsoidal models, finding that the $\eta$-linear models already provide an excellent approximation of the numerical solutions. As an application, the $\eta$-linear Plummer model (with a central black hole) is used for the phenomenological interpretation of the dynamics of the weakly flattened and rotating globular cluster NGC 4372, confirming that this system cannot be interpreted as an isotropic rotator, a conclusion reached previously with more sophisticated studies.

Observational breakthroughs in the exoplanet field of the last decade motivated the development of numerous theoretical models describing atmospheres and mass loss, which is believed to be one of the main drivers of planetary evolution. We aim to outline for which types of close-in planets in the Neptune-mass range the accurate treatment of photoionisation effects is most relevant concerning atmospheric escape and the parameters relevant for interpreting observations. We developed the CHAIN (Cloudy e Hydro Ancora INsieme) model combining 1D hydrodynamic upper atmosphere model with the non-LTE photoionisation and radiative transfer code Cloudy accounting for photochemistry, detailed atomic level populations, and chemical reactions for all elements up to zinc. We apply CHAIN to model the upper atmospheres of a range of Neptune-like planets with masses between 1 and 50 M$_{\oplus}$, varying also the orbital parameters. For the majority of warm and hot Neptunes, we find slower and denser outflows, with lower ion fractions, compared to the predictions of the hydrodynamic model alone. Furthermore, we find significantly different temperature profiles between CHAIN and the hydrodynamic model alone, though the peak values are similar for similar atmospheric compositions. The mass-loss rates predicted by CHAIN are higher for hot, strongly irradiated planets and lower for more moderate planets. All differences between the two models are strongly correlated with the amount of high-energy irradiation. Finally, we find that the hydrodynamic effects impact significantly ionisation and heating. The impact of the precise photoionisation treatment provided by Cloudy strongly depends on the system parameters. This suggests that some of the simplifications typically employed in hydrodynamic modelling might lead to systematic errors when studying planetary atmospheres, even at a population-wide level.

We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251\deg $\le l \le$ 358\deg and 2\deg $\le l \le$ 61\deg at $|b| \le 1.5\deg $). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $\mu$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance.

Ocean worlds, or icy bodies in the outer solar system that have or once had subsurface liquid water oceans, are among the most compelling topics of astrobiology. Typically, confirming the existence of a subsurface ocean requires close spacecraft observations. However, combining our understanding of the chemistry that takes place in a subsurface ocean with our knowledge of the building blocks that formed potential ocean worlds provides an opportunity to identify tracers of endogenic activity in the surface volatiles of Pluto and Triton. We show here that the current composition of the volatiles on the surfaces and in the atmospheres of Pluto and Triton are deficient in carbon, which can only be explained by the loss of CH4 through a combination of aqueous chemistry and atmospheric processes. Furthermore, we find that the relative nitrogen and water abundances are within the range observed in building block analogs, comets, and chondrites. A lower limit for N/Ar in Pluto's atmosphere also suggests source building blocks that have a cometary or chondritic composition, all pointing to an origin for their nitrogen as NH3 or organics. Triton's lower abundance of CH4 compared to Pluto, and the detection of CO2 at Triton but not at Pluto points to aqueous chemistry in a subsurface ocean that was more efficient at Triton than Pluto. These results have applications to other large Kuiper Belt objects as well as the assessment of formation locations and times for the four giant planets given future probe measurements of noble gas abundances and isotope ratios.

The possible existence of stellar halos in low-mass galaxies is being intensely discussed nowadays after some recent discoveries of stars located in the outskirts of dwarf galaxies of the Local Group. RR Lyrae stars can be used to identify the extent of these structures, taking advantage of the minimization of foreground contamination they provide. In this work we use RR Lyrae stars obtained from Gaia DR3, DES, ZTF, and Pan-STARRS1 to explore the outskirts of $45$ ultra-faint dwarf galaxies. We associate the stars with a host galaxy based on their angular separations, magnitudes and proper motions. We find a total of $120$ RR Lyrae stars that belong to $21$ different galaxies in our sample. We report seven new RR Lyrae stars in six ultra-faint dwarf galaxies (Hydrus I, Ursa Major I, Ursa Major II, Grus II, Eridanus II and Tucana II). We found a large number of new possible members in Bootes I and Bootes III as well, but some of them may actually belong to the nearby Sagittarius stream. Adding to our list of $120$ RR Lyrae stars the observations of other ultra-faint dwarf galaxies that were out of the reach of our search, we find that at least $10$ of these galaxies have RR Lyrae stars located at farther distances than $4$ times their respective half-light radius, which implies that at least $33\%$ of the 30 ultra-faint dwarfs with RR Lyrae star population have extended stellar populations.

Estimators for $n$-point clustering statistics in Fourier-space demand that modern surveys of large-scale structure be transformed to Cartesian coordinates to perform Fast Fourier Transforms (FFTs). In this work, we explore this transformation in the context of pixelised line intensity maps (LIM), highlighting potential biasing effects on power spectrum measurements. Current analyses often avoid a complete resampling of the data by approximating survey geometry as rectangular in Cartesian space, an increasingly inaccurate assumption for modern wide-sky surveys. Our simulations of a $20\,{\times}\,20\,\text{deg}^2$ 21cm LIM survey at $0.34\,{<}\,z\,{<}\,0.54$ show this assumption biases power spectrum measurements by ${>}\,20\%$ across all scales. We therefore present a more robust framework for regridding the voxel intensities onto a 3D FFT field by coordinate transforming large numbers of Monte-Carlo sampling particles. Whilst this unbiases power spectrum measurements on large scales, smaller-scale discrepancies remain, caused by structure smoothing and aliasing from separations unresolved by the grid. To correct these effects, we introduce modelling techniques, higher-order particle assignments, and interlaced FFT grids to suppress the aliased power. Using a Piecewise Cubic Spline (PCS) particle assignment and an interlaced FFT field, we achieve sub-percent accuracy up to 80% of the Nyquist frequency for our 21cm LIM simulations. We find a more subtle hierarchical improvement in results for higher-order assignment schemes, relative to the gains made for galaxy surveys, which we attribute to the extra complexity in LIM from additional discretising steps. Python code accompanying this paper is available at github.com/stevecunnington/gridimp.

We analyze time-dependent dark energy equations of state through linear and nonlinear structure formation and their quintessence potentials, characterized by fast, recent transitions, inspired by parameter space studies of Horndeski models. The influence of dark energy on structures comes from modifications to the background expansion rate and from perturbations as well. In order to compute the structures growth, we employ a generalization of the \emph{spherical collapse} formalism that includes perturbations of fluids with pressure. We numerically solve the equations of motion for the perturbations and the field. Our analysis suggests that a true Heaviside step transition is a good approximation for most of the considered models, since most of the quantities weakly depend on the transition speed. We find that transitions occurring at redshifts $z_{\rm t}\gtrsim 2$ cannot be distinguished from the $\Lambda$CDM model if dark energy is freezing, i.e., the corresponding equation of state tends to $-1$. For fast, recent transitions, the redshift at which the properties of dark energy have the most significant effect is $z=0.6\pm 0.2$. We also find that in the freezing regime, the $\sigma_8$ values can be lowered by about $8\%$, suggesting that those models could relieve the $\sigma_8$-tension. Additionally, freezing models generally predict faster late-time merging rates but a lower number of massive galaxies at $z=0$. Finally, the matter power spectrum for smooth dark energy shows a low-wavenumber peak which is absent in the clustering case.

We apply a novel model based on convolutional neural networks (CNNs) to identify gravitationally-lensed galaxies in multi-band imaging of the Hyper Suprime Cam Subaru Strategic Program (HSC-SSP) Survey. The trained model is applied to a parent sample of 2 350 061 galaxies selected from the $\sim$ 800 deg$^2$ Wide area of the HSC-SSP Public Data Release 2. The galaxies in HSC Wide are selected based on stringent pre-selection criteria, such as multiband magnitudes, stellar mass, star formation rate, extendedness limit, photometric redshift range, etc. Initially, the CNNs provide a total of 20 241 cutouts with a score greater than 0.9, but this number is subsequently reduced to 1 522 cutouts by removing definite non-lenses for further inspection by human eyes. We discover 43 definite and 269 probable lenses, of which 97 are completely new. In addition, out of 880 potential lenses, we recovered 289 known systems in the literature. We identify 143 candidates from the known systems that had higher confidence in previous searches. Our model can also recover 285 candidate galaxy-scale lenses from the Survey of Gravitationally lensed Objects in HSC Imaging (SuGOHI), where a single foreground galaxy acts as the deflector. Even though group-scale and cluster-scale lens systems were not included in the training, a sample of 32 SuGOHI-c (i.e., group/cluster-scale systems) lens candidates was retrieved. Our discoveries will be useful for ongoing and planned spectroscopic surveys, such as the Subaru Prime Focus Spectrograph project, to measure lens and source redshifts in order to enable detailed lens modelling.

A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any pre-existing gas disc to be removed, via fly-bys or ejective feedback from an active galactic nucleus. In this Letter, we study the evolution of a Milky Way mass galaxy in the FIREbox cosmological volume that displays a thin, counter-rotating gas disc with respect to its stellar component at low redshift. In contrast to scenarios involving gas ejection, we find that pre-existing gas is mainly removed via the conversion into stars in a central starburst, triggered by a merging satellite galaxy. The newly-accreted, counter-rotating gas eventually settles into a kinematically misaligned disc. About 4.4 (8 out of 182) of FIREbox galaxies with stellar masses larger than 5e9 Msun at z=0 exhibit gas-star kinematic misalignment. In all cases, we identify central starburst-driven depletion as the main reason for the removal of the pre-existing co-rotating gas component, with no need for feedback from, e.g., a central active black hole. However, during the starburst, the gas is funneled towards the central regions, likely enhancing black hole activity. By comparing the fraction of misaligned discs between FIREbox and other simulations and observations, we conclude that this channel might have a non-negligible role in inducing kinematic misalignment in galaxies.

We perform a data-constrained simulation with the zero-$\beta$ assumption to study the mechanisms of strong rotation and failed eruption of a filament in active region 11474 on 2012 May 5 observed by Solar Dynamics Observatory and Solar Terrestrial Relations Observatory. The initial magnetic field is provided by nonlinear force-free field extrapolation, which is reconstructed by the regularized Biot-Savart laws and magnetofrictional method. Our simulation reproduces most observational features very well, e.g., the filament large-angle rotation of about $130 ^{\circ}$, the confined eruption and the flare ribbons, allowing us to analyze the underlying physical processes behind observations. We discover two flux ropes in the sigmoid system, an upper flux rope (MFR1) and a lower flux rope (MFR2), which correspond to the filament and hot channel in observations, respectively. Both flux ropes undergo confined eruptions. MFR2 grows by tether-cutting reconnection during the eruption. The rotation of MFR1 is related to the shear-field component along the axis. The toroidal field tension force and the non-axisymmetry forces confine the eruption of MFR1. We also suggest that the mutual interaction between MFR1 and MFR2 contributes to the large-angle rotation and the eruption failure. In addition, we calculate the temporal evolution of the twist and writhe of MFR1, which is a hint of probable reversal rotation.

Till today, the nature of Dark Matter (DM) remains elusive despite all our efforts. This missing matter of the universe has not been observed by the already operating DM direct-detection experiments, but we can infer its gravitational effects. Galaxies and clusters of galaxies are most likely to contain DM trapped to their gravitational field. This leads us to the natural assumption that compact objects might contain DM too. Among the compact objects exist in galaxies, neutron stars considered as natural laboratories, where theories can be tested, and observational data can be received. Thus, many models of DM they have proposed it's presence in those stars. By employing the two fluid model, we discovered a stable area in the M-R diagram of a celestial formation consisting of neutron and DM that is substantial in size and vast in dimensions. This formation spans hundreds of kilometers in diameter and possesses a mass equivalent to 100 or more times that of our sun. To elucidate, this entity resembles an enormous celestial body of DM, with a neutron star at its core. This implies that a supramassive stellar compact entity can exist without encountering any issues of stability and without undergoing a collapse into a black hole. In any case, the present theoretical prediction can, if combined with corresponding observations, shed light on the existence of DM and even more on its basic properties.

Nova 1670 is a historical transient bearing strong similarities to a recently-recognized type of stellar eruptions known as red novae, which are thought to be powered by stellar mergers. The remnant of the transient, CK Vul, is observable today mainly through cool circumstellar gas and dust, and recombining plasma, but we have no direct view on the stellar object. Within the merger hypothesis, we aim to infer the most likely makeup of the progenitor system that resulted in Nova 1670. We collect and summarize the literature data on the physical properties of the outburst and the remnant, and on the chemical composition of the circumstellar material which resulted from optical and submillimeter observations of the circumstellar gas of CK Vul. Simple simulations yield the form and level of mixing of material associated with the merger. Products of nuclear burning are identified, among them ashes of H burning in the CNO cycles and the MgAl chain, as well as of partial He burning. Based on the luminosity and chemical composition of the remnant, we find that the progenitor primary had to be an evolutionarily advanced red-giant branch star of a mass of 1-2 M$_{\odot}$. The secondary was either a very similar giant, or a He white dwarf. While the eruption event was mainly powered by accretion, we estimate that about 12% of total energy might have come from He burning activated during the merger. The coalescence of a first-ascent giant with a He white dwarf created a star with a rather unique internal structure and composition, which resemble these of early R-type carbon stars. Nova 1670 was the result of a merger between a He white dwarf and a first-ascent red giant and is likely now evolving to become an early R-type carbon star.

We show how one can test the cosmological Poisson equation by requiring only the validity of three main assumptions: the energy-momentum conservation equations of matter, the equivalence principle, and the cosmological principle. We first point out that one can only measure the combination ${\mathcal M}\equiv \Omega_m^{(0)}\mu$, where $\mu$ quantifies the deviation of the Poisson equation from the standard one and $\Omega_m^{(0)}$ is the fraction of matter density at present. Then we employ a recent model-independent forecast for the growth rate $f(z)$ and the expansion rate $E(z)$ to obtain constraints on ${\mathcal M}$ for a survey that approximates a combination of the Dark Energy Spectroscopic Instrument (DESI) and Euclid. We conclude that a constant ${\mathcal M}$ can be measured with a relative error $\sigma_{\mathcal{M}}=4.1\%$, while if ${\mathcal M}$ is arbitrarily varying in redshift, it can be measured to within $19.3\%$ (1 $\sigma$ c.l.) at redshift $z=0.9$, and 20-30\% up to $z=1.5$. We also project our constraints on some parametrizations of ${\mathcal M}$ proposed in literature, while still maintaining model-independence for the background expansion, the power spectrum shape, and the non-linear corrections. Generally speaking, as expected, we find much weaker model-independent constraints than found so far for such models. This means that the cosmological Poisson equation remains quite open to various alternative gravity and dark energy models.

Non-linear structure formation for fermionic dark matter particles leads to dark matter density profiles with a degenerate compact core surrounded by a diluted halo. For a given fermion mass, the core has a critical mass that collapses into a supermassive black hole (SMBH). Galactic dynamics constraints suggest a $\sim 100$ keV/$c^2$ fermion, which leads to $\sim 10^7 M_\odot$ critical core mass. Here, we show that baryonic (ordinary) matter accretion drives an initially stable dark matter core to SMBH formation and determine the accreted mass threshold that induces it. Baryonic gas density $\rho_b$ and velocity $v_b$ inferred from cosmological hydro-simulations and observations produce sub-Eddington accretion rates triggering the baryon-induced collapse in less than a Gyr. This process produces active galactic nuclei in galaxy mergers and the high-redshift Universe. For TXS 2116-077, merging with a nearby galaxy, the observed $3\times 10^7 M_\odot$ SMBH, for $Q_b = \rho_b/v_b^3 = 0.125 M_\odot/(100 \text{km/s pc})^3$, forms in $\approx 0.6$ Gyr, consistent with the $0.5$-$2$ Gyr merger timescale and younger jet. For the farthest central SMBH detected by the \textit{Chandra} X-ray satellite in the $z=10.3$ UHZ1 galaxy observed by the James Webb Space Telescope (\textit{JWST}), the mechanism leads to a $4\times 10^7 M_\odot$ SMBH in $87$-$187$ Myr, starting the accretion at $z=12$-$15$. The baryon-induced collapse can also explain the $\approx 10^7$-$10^8 M_\odot$ SMBHs revealed by the JWST at $z\approx 4$-$6$. After its formation, the SMBH can grow to a few $10^9 M_\odot$ in timescales shorter than a Gyr via sub-Eddington baryonic mass accretion.

The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including 9 previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6)\Msun and 1.412(6)\Msun, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves: $\dot{P}_{\rm B} = -1.37(7)\times10^{-12}$, which is in agreement with the rate predicted by general relativity for these masses.

Venus may have had both an Earth-like climate as well as extensive water oceans and active (or incipient) plate tectonics for an extended interval of its history. The topographical power spectrum of Venus provides important clues to the planet's past evolution. By drawing detailed contrast with the strong low-order odd-$l$ dominated global topography of Earth, we demonstrate that the relatively flat Venusian topography can be interpreted to have arisen from the transition from active terrestrial-like plate tectonics to the current stagnant lid configuration at a time $\tau = 544^{+886}_{-193}$ million years before present. This scenario is plausible if loss of oceans and the attendant transition to a CO$_2$-dominated atmosphere were accompanied by rapid continental-scale erosion, followed by gradual lava resurfacing at an outflow rate $\sim$ 1 km$^{3}$ yr$^{-1}$. We study Venus' proposed topographical relaxation with a global diffusion-like model that adopts terrestrial erosion rates scaled to account for the increased rainfall and temperatures that would accompany a planet-wide transition from an Earth-like climate to the runaway greenhouse climate that could ultimately yield present-day Venus, with an estimate of $5.1^{+1.8}_{-1.1}$ Myr if the global erosion operated as efficiently as that of a typical bedrock river basin on Earth.

Measurements with widefield radio interferometers often include the near-infinite gradient between the sky and the horizon. This causes aliasing inherent to the measurement itself, and is purely a consequence of the Fourier basis. For this reason, the horizon is often attenuated by the instrumental beam down to levels deemed inconsequential. However, this effect is enhanced via our own Galactic plane as it sets over the course of a night. We show all-sky simulations of the Galactic plane setting in a radio interferometer in detail for the first time. We then apply these simulations to the Murchison Widefield Array to show that a beam attenuation of 0.1% is not sufficient in some precision science cases. We determine that the noise statistics of the residual data image are drastically more Gaussian with aliasing removal, and explore consequences in simulation for cataloging of extragalactic sources and 21-cm Epoch of Reionization detection via the power spectrum.

$\it{Context.}$ Various high-energy phenomena in the Universe are associated with blazars, powerful active galaxies with jets pointing to the observer. Novel results relating blazars to high-energy neutrinos, cosmic rays, and even possible manifestations of new particle physics, are often based on statistical analyses of blazar samples, and uniform sky coverage is important for many of these studies. $\it{Aims.}$ Here, we construct a uniform full-sky catalog of blazars selected by their optical emission. $\it{Methods.}$ We define criteria of isotropy, making a special effort to cover the Galactic plane region, and compile an isotropic sample of blazars with GAIA optical magnitudes $G<18^{\rm m}$, corrected for the Galactic absorption. The sources are taken from full-sky samples selected by parsec-scale radio emission or by high-energy gamma-ray flux, both being known to efficiently select blazar-like objects. $\it{Results.}$ We present a catalog of 651 optically bright blazars, uniformly distributed in the sky, together with their radio, optical, X-ray and gamma-ray fluxes, and an isotropic sample of 336 confirmed BL Lac type objects. $\it{Conclusions.}$ This catalog may be used in future statistical studies of energetic neutrinos, cosmic rays and gamma rays.

Cosmic rays with energies above $10^{19}$ eV, observed in 1999-2004 by the High Resolution Fly's Eye (HiRes) experiment in the stereoscopic mode, were found to correlate with directions to distant BL Lac type objects (BL Lacs), suggesting non-standard neutral particles travelling for cosmological distances without attenuation. This effect could not be tested by newer experiments because of their inferior angular resolution. The distribution in the sky of BL Lacs associated with cosmic rays was found to deviate from isotropy, which might give a clue to the interpretation of the observed anomaly. However, previous studies made use of a sample of BL Lacs which was anisotropic by itself, thus complicating these interpretations. Here, we use a recently compiled isotropic sample of BL Lacs and the same HiRes data to confirm the presence of correlations and to strengthen the case for the local large-scale structure pattern in the distribution of the correlated events in the sky. Further tests of the anomaly await new precise cosmic-ray data.

Cosmological simulations play a crucial role in elucidating the effect of physical parameters on the statistics of fields and on constraining parameters given information on density fields. We leverage diffusion generative models to address two tasks of importance to cosmology -- as an emulator for cold dark matter density fields conditional on input cosmological parameters $\Omega_m$ and $\sigma_8$, and as a parameter inference model that can return constraints on the cosmological parameters of an input field. We show that the model is able to generate fields with power spectra that are consistent with those of the simulated target distribution, and capture the subtle effect of each parameter on modulations in the power spectrum. We additionally explore their utility as parameter inference models and find that we can obtain tight constraints on cosmological parameters.

The magnetic field through the magnetic reconnection process affects the dynamics and structure of astrophysical systems. Numerical simulations are the tools to study the evolution of these systems. However, the resolution, dimensions, resistivity, and turbulence of the system are some important parameters to take into account in the simulations. In this paper, we investigate the evolution of magnetic energy in astrophysical simulations by performing a standard test problem for MHD codes, Orszag-Tang. We estimate the numerical dissipation in the simulations using state-of-the-art numerical simulation code in astrophysics, PLUTO. The estimated numerical resistivity in 2D simulations corresponds to the Lundquist number $\approx 10^{4}$ in the resolution of $512\times512$ grid cells. It is also shown that the plasmoid unstable reconnection layer can be resolved with sufficient resolutions. Our analysis demonstrates that in non-relativistic magnetohydrodynamics simulations, magnetic and kinetic energies undergo conversion into internal energy, resulting in plasma heating.

We present a novel multimodal dataset developed by expert astronomers to automate the detection and localisation of multi-component extended radio galaxies and their corresponding infrared hosts. The dataset comprises 4,155 instances of galaxies in 2,800 images with both radio and infrared modalities. Each instance contains information on the extended radio galaxy class, its corresponding bounding box that encompasses all of its components, pixel-level segmentation mask, and the position of its corresponding infrared host galaxy. Our dataset is the first publicly accessible dataset that includes images from a highly sensitive radio telescope, infrared satellite, and instance-level annotations for their identification. We benchmark several object detection algorithms on the dataset and propose a novel multimodal approach to identify radio galaxies and the positions of infrared hosts simultaneously.

The Jovian magnetic field, being the strongest and largest planetary one in the solar system, could offer us new insights into possible microscopic scale new physics, such as a non-zero mass of the Standard Model (SM) photon or a light dark photon kinetically mixing with the SM photon. We employ the immense data set from the latest Juno mission, which provides us unprecedented information about the magnetic field of the gas giant, together with a more rigorous statistical approach compared to the literature, to set strong constraints on the dark photon mass and kinetic mixing parameter, as well as the SM photon mass. The constraint on the dark photon parameters is independent of whether dark photon is (part of) dark matter or not, and serves as the most stringent one in a certain regime of the parameter space.

We revisit the electroweak phase transition in the scalar singlet extension of the standard model with a $\mathbb{Z}_2$ symmetry. In significant parts of the parameter space the phase transition occurs in two steps - including canonical benchmarks used in experimental projections for gravitational waves. Domain walls produced in the first step of the transition seed the final step to the electroweak vacuum, an effect which is typically neglected but leads to an exponentially enhanced tunnelling rate. We improve previous results obtained for the seeded transition, which made use of the thin-wall or high temperature approximations, by using the mountain pass algorithm that was recently proposed as a useful tool for seeded processes. We then determine the predictions of the seeded transition for the latent heat, bubble size and characteristic time scale of the transition. Differences compared to homogeneous transitions are most pronounced when there are relatively few domain walls per hubble patch, potentially leading to an enhanced gravitational wave signal. We also provide a derivation of the percolation criteria for a generic seeded transition, which applies to the domain wall seeds we consider as well as to strings and monopoles.

Quantum fluctuations in spacetime can, in some cases, lead to distortion in astronomical images of faraway objects. In particular, a stochastic model of quantum gravity predicts an accumulated fluctuation in the path length $\Delta L$ with variance $\langle \Delta L^2\rangle\sim l_pL$ over a distance $L$, similar to a random walk, and assuming no spatial correlation above length $l_p$; it has been argued that such an effect is ruled out by observation of sharp images from distant stars. However, in other theories, such as the pixellon (modeled on the Verlinde-Zurek (VZ) effect), quantum fluctuations can still accumulate as in the random walk model while simultaneously having large distance correlations in the fluctuations. Using renormalization by analytic continuation, we derive the correlation transverse to the light propagation, and show that image distortion effects in the pixellon model are strongly suppressed in comparison to the random walk model, thus evading all existing and future constraints. We also find that the diffraction of light rays does not lead to qualitative changes in the blurring effect.

We present a new search for dark matter using planetary atmospheres. We point out that annihilating dark matter in planets can produce ionizing radiation, which can lead to excess production of ionospheric $H_3^+$. We apply this search strategy to the night side of Jupiter near the equator. The night side has zero solar irradiation, and low latitudes are sufficiently far from ionizing auroras, leading to an effectively background-free search. We use Cassini data on ionospheric $H_3^+$ emission collected 3 hours either side of Jovian midnight, during its flyby in 2000, and set novel constraints on the dark matter-nucleon scattering cross section down to about $10^{-38}$ cm$^2$. We also highlight that dark matter atmospheric ionization may be detected in Jovian exoplanets using future high-precision measurements of planetary spectra.

The calculation of loop corrections to the correlation functions of quantum fields during inflation or in the de~Sitter background presents greater challenges than in flat space due to the more complicated form of the mode functions. While in flat space highly sophisticated approaches to Feynman integrals exist, similar tools still remain to be developed for cosmological correlators. However, usually only their late-time limit is of interest. We introduce the method-of-region expansion for cosmological correlators as a tool to extract the late-time limit, and illustrate it with several examples for the interacting, massless, minimally coupled scalar field in de~Sitter space. In particular, we consider the in-in correlator $\langle\phi^2(\eta,q)\phi(\eta,k_1)\phi(\eta,k_2)\rangle$, whose region structure is relevant to anomalous dimensions and matching coefficients in Soft de Sitter effective theory.

Previous works have argued that future gravitational-wave detectors will be able to probe the properties of astrophysical environments where binary coalesce, including accretion disks, but also dark matter structures. Most analyses have resorted to a Newtonian modelling of the environmental effects, which are not suited to study extreme-mass-ratio inspirals immersed in structures of ultra-light bosons. In this letter, we use relativistic perturbation theory to consistently study these systems in spherical symmetry. We compute the flux of scalar particles and the rate at which orbital energy (and angular momentum) is dissipated via gravitational radiation and depletion of scalars, i.e. dynamical friction. Our results suggest that the Laser Inteferometer Space Antenna will be able to probe ultra-light dark matter structures in the Galaxy by tracking the phase of extreme-mass-ratio inspirals.

We investigate how superradiance affects the generation of baryon asymmetry in a universe with rotating primordial black holes, considering a scenario where a scalar boson is coupled to the heavy right-handed neutrinos. We identify the regions of the parameter space where the scalar production is enhanced due to superradiance. This enhancement, coupled with the subsequent decay of the scalar into right handed neutrinos, results in the non-thermal creation of lepton asymmetry. We show that successful leptogenesis is achieved for masses of primordial black holes in the range of order $O(0.1~{\rm g}) - O(10~{\rm g})$ and the lightest of the heavy neutrino masses, $M_N \sim O(10^{12})~{\rm GeV}$. Consequently, regions of the parameter space, which in the case of Schwarzchild PBHs were incompatible with viable leptogenesis, can produce the observed matter-antimatter asymmetry.

In this note, we prove analytic bounds on the equation of state of a cosmological fluid composed of an arbitrary number of canonical scalars evolving in a negative multi-exponential potential. Because of the negative energy, the universe is contracting and eventually undergoes a big crunch. A contracting universe is a fundamental feature of models of ekpyrosis and cyclic universes, which are a proposed alternative to cosmic inflation to solve the flatness and horizon problems. Our analytic bounds set quantitative constraints on general effective theories of ekpyrosis. When applied to specific top-down constructions, our bounds can be used to determine whether ekpyrosis could in principle be realized. We point out some possible sources of tension in realizing the ekpyrotic universe in controlled constructions of string theory.

We propose a gauged $U(1)_{B-L}$ version of the light Dirac neutrino portal dark matter. The $U(1)_{B-L}$ symmetry provides a UV completion by naturally accommodating three right handed neutrinos from anomaly cancellation requirements which, in combination with the left handed neutrinos, form the sub-eV Dirac neutrinos after electroweak symmetry breaking. The particle content and the gauge charges are chosen in such a way that light neutrinos remain purely Dirac and dark matter, a gauge singlet Dirac fermion, remain stable. We consider both thermal and non-thermal production possibilities of dark matter and correlate the corresponding parameter space with the one within reach of future cosmic microwave background (CMB) experiments sensitive to enhanced relativistic degrees of freedom $\Delta N_{\rm eff}$. The interplay of dark matter, CMB, structure formation and other terrestrial constraints keep the scenario very predictive leading the $U(1)_{B-L}$ parameter space into tight corners.

We demonstrate an efficient magnetic dynamo due to precession-driven hydrodynamic turbulence in the local model. Dynamo growth rate increases with Poincar\'{e} and magnetic Prandtl numbers. Spectral analysis shows that the dynamo acts over a broad range of scales: at large (system size) and intermediate scales it is driven by 2D vortices and shear of the background precessional flow, while at smaller scales it is mainly driven by 3D inertial waves. These results are important for understanding magnetic field generation and amplification in precessing planets and stars.

Accurate estimation of thermospheric density is critical for precise modeling of satellite drag forces in low Earth orbit (LEO). Improving this estimation is crucial to tasks such as state estimation, collision avoidance, and re-entry calculations. The largest source of uncertainty in determining thermospheric density is modeling the effects of space weather driven by solar and geomagnetic activity. Current operational models rely on ground-based proxy indices which imperfectly correlate with the complexity of solar outputs and geomagnetic responses. In this work, we directly incorporate NASA's Solar Dynamics Observatory (SDO) extreme ultraviolet (EUV) spectral images into a neural thermospheric density model to determine whether the predictive performance of the model is increased by using space-based EUV imagery data instead of, or in addition to, the ground-based proxy indices. We demonstrate that EUV imagery can enable predictions with much higher temporal resolution and replace ground-based proxies while significantly increasing performance relative to current operational models. Our method paves the way for assimilating EUV image data into operational thermospheric density forecasting models for use in LEO satellite navigation processes.

We present a new mechanism to generate a coherently oscillating dark vector field from axion-SU(2) gauge field dynamics during inflation. The SU(2) gauge field acquires a nonzero background sourced by an axion during inflation, and it acquires a mass through spontaneous symmetry breaking after inflation. We find that the coherent oscillation of the dark vector field can account for dark matter in the mass range of $10^{-13}-1$ eV in a minimal setup. In a more involved scenario, the range can be wider down to the fuzzy dark matter region. One of the dark vector fields can be identified as the dark photon, in which case this mechanism evades the notorious constraints for isocurvature perturbation, statistical anisotropy, and the absence of ghosts that exist in the usual misalignment production scenarios. Phenomenological implications are discussed.

A transport-like framework for the study of magnetic reconnection mediated by self-driven turbulence is proposed, based on timescale separation between the reconnection time and the characteristic timescale of the turbulent fluctuations which arise in the reconnection layer. We argue that the mean fields remain on MHD scales even in collisionless cases. These observations provide theoretical justification for an efficient computational approach to this problem, which we discuss.

Motivated by precision computations of neutrino decoupling at MeV temperatures, we show how QED corrections to the thermal neutrino interaction rate can be related to the electron-positron spectral function as well as an effective $\bar{\nu}\nu\gamma$ vertex. The spectral function is needed both in a timelike and in a spacelike domain, and for both of its physical polarization states (transverse and longitudinal with respect to spatial momentum). Incorporating an NLO evaluation of this spectral function, an estimate of the $\bar{\nu}\nu\gamma$ vertex, and HTL resummation of scatterings mediated by soft Bose-enhanced $t$-channel photons, we compute the interaction rate as a function of the neutrino momentum and flavour. Effects on the $ -(0...2)\%$ level are found, noticeably smaller than a previous estimate of a related quantity.

The Brans-Dicke (BD) theory is one of the simplest scalar-tensor theories, which has potential relations with dark matter, dark energy, inflation, and primordial nucleosynthesis. The strongest constraint on the BD coupling constant is provided by the Shapiro time delay measurement in the solar system by Cassini. Constraints from gravitational wave (GW) events are subject to asymmetric binaries (binaries with different ``sensitivity'' such as neutron star-black hole (NSBH), white dwarf-neutron star, or white dwarf-black hole binary). The third Gravitational-Wave Transient Catalog (GWTC) reports an NSBH merger event, GW200115, making it possible to constrain BD by GW. With the aid of this source and Bayesian Markov-chain Monte Carlo (MCMC) analyses, we derive the 90\% credible lower bound on the modified parameter of BD as $\omega_{\rm BD}>4.75$ by using dominant (2, 2) mode correction. Extending to scalar-tensor theory, we have the constraint $\varphi_{-2}>-7.94\times10^{-4}$, which is consistent with the bound $7.3\times10^{-4}$ from the LIGO Laboratory. Usually, asymmetric binary systems have a significant mass ratio; in such cases, higher harmonic modes cannot be neglected. Our work considers higher harmonic corrections from BD and provides a tighter constraint of $\omega_{\rm BD}>5.06$. We also consider the suspected event GW190426\_152155 as an NSBH event. We find $\omega_{\rm BD}>1.25$ when employing the dominant (2, 2) mode only and $\omega_{\rm BD}>1.47$ when including the higher harmonic modes. Additionally, we take into account a BD-like theory, known as screened modified gravity (SMG), and provide the coupling constant constraints from both with and without high mode corrections, by using data from both GW200115 and GW190426\_152155.

We consider axion cloud decay due to the axion-photon conversion with multi-pole background magnetic fields. We focus on the $\ell=m=1$ and $n=2$ mode for the axion field configuration since it has the largest growth rate associated with superradiant instability. Under the existence of a background multi-pole magnetic field, the axion field can be converted into the electromagnetic field through the axion-photon coupling. Then the decay rate due to the dissipation of the converted photons is calculated in a successive approximation. We found that the decay rate is significantly dependent on the azimuthal quantum number characterizing the background magnetic field, and can be comparable to or larger than the growth rate of the superradiant instability.

Gas Electron Multiplier (GEM) detectors are crucial for enabling high-resolution X-ray polarisation of astrophysical sources when coupled to custom pixel readout ASIC in Gas Pixel Detectors (GPD), as in the Imaging X-ray Polarimetry Explorer (IXPE), the Polarlight cubesat pathfinder and the PFA telescope onboard the future large enhanced X-ray Timing and Polarimetry (eXTP) Chinese mission. The R&D efforts of the IXPE collaboration have resulted in mature GPD technology. However, limitations in the classical wet-etch or laser-drilled fabrication process of GEMs motivated our exploration of alternative methods. This work focuses on investigating a plasma-based etching approach for fabricating GEM patterns at Fondazione Bruno Kessler (FBK). The objective is to improve the aspect ratio of the GEM holes, to mitigate the charging of the GEM dielectric which generates rate-dependent gain changes. Unlike the traditional wet-etch process, Reactive Ion Etching (RIE) enables more vertical etching profiles and thus better aspect ratios. Moreover, the RIE process promises to overcome non-uniformities in the GEM hole patterns which are believed to cause systemic effects in the azimuthal response of GPDs equipped with either laser-drilled or wet-etch GEMs. We present a GEM geometry with 20 micrometers in diameter and 50 micrometers pitch, accompanied by extensive characterisation (SEM and PFIB) of the structural features and aspect ratios. The collaboration with INFN Pisa and Turin enabled us to compare the electrical properties of these detectors and test their performance in their use as electron multipliers in GPDs. Although this R&D work is in its initial stages, it holds promise for enhancing the sensitivity of the IXPE mission in X-ray polarimetry measurements through GEM pattern with more vertical hole profiles.

Plug-and-Play (PnP) algorithms are appealing alternatives to proximal algorithms when solving inverse imaging problems. By learning a Deep Neural Network (DNN) behaving as a proximal operator, one waives the computational complexity of optimisation algorithms induced by sophisticated image priors, and the sub-optimality of handcrafted priors compared to DNNs. At the same time, these methods inherit from the versatility of optimisation algorithms allowing to minimise a large class of objective functions. Such features are highly desirable in radio-interferometric (RI) imaging in astronomy, where the data size, the ill-posedness of the problem and the dynamic range of the target reconstruction are critical. In a previous work, we introduced a class of convergent PnP algorithms, dubbed AIRI, relying on a forward-backward backbone, with a differentiable data-fidelity term and dynamic range-specific denoisers trained on highly pre-processed unrelated optical astronomy images. Here, we show that AIRI algorithms can benefit from a constrained data fidelity term at the mere cost of transferring to a primal-dual forward-backward algorithmic backbone. Moreover, we show that AIRI algorithms are robust to strong variations in the nature of the training dataset: denoisers trained on MRI images yield similar reconstructions to those trained on astronomical data. We additionally quantify the model uncertainty introduced by the randomness in the training process and suggest that AIRI algorithms are robust to model uncertainty. Eventually, we propose an exhaustive comparison with methods from the radio-astronomical imaging literature and show the superiority of the proposed method over the state-of-the-art.

General relativity (GR) may be modified by adding an extra warped noncompact spatial dimension such that it is indistinguishable from GR as far as local tests of gravity are concerned. However, such a modified theory of gravity should have intriguing consequences on various aspects of black holes which can help us probe how the presence of an extra dimension affects gravity in the strong field regime, and help us distinguish it from GR. Therefore, we have studied massive scalar perturbations of four-dimensional rotating black holes in the Randall Sundrum II braneworld scenario. These black holes are endowed with a tidal charge that contains information pertaining to the extra spatial dimension in the braneworld model. Such black hole spacetimes are also noteworthy because they permit the black hole's rotation parameter to exceed unity, a possibility strictly forbidden by the general theory of relativity. Consequently, they offer valuable insights into exploring the repercussions of modifications to Einstein's theory through possible future observations. Our approach involves the numerical solution of the perturbed field equations using the continued fractions method. First we ascertain the quasinormal mode spectra of the rotating braneworld black hole. We then thoroughly investigate the existence of quasibound states and the associated superradiant instability. Such a superradiant cloud of bosons around the black hole are called gravitational atoms and are invaluable observational probes of ultralight bosonic particles predicted in various extensions of the Standard Model of particle physics. In comparison to four-dimensional Kerr black hole, we report distinctive signatures of the tidal charge and the rotation parameter, which manifest as signals of the extra dimension on both the quasinormal mode and the formation of gravitational atoms.

Templates modeling just the dominant mode of gravitational radiation are generally sufficient for the unbiased parameter inference of near-equal-mass compact binary mergers. However, neglecting the subdominant modes can bias the inference if the binary is significantly asymmetric, very massive, or has misaligned spins. In this work, we explore if neglecting these subdominant modes in the parameter estimation of non-spinning binary black hole mergers can bias the inference of their population-level properties such as mass and merger redshift distributions. Assuming the design sensitivity of advanced LIGO-Virgo detector network, we find that neglecting subdominant modes will not cause a significant bias in the population inference, although including them will provide more precise estimates. This is primarily due to the fact that asymmetric binaries are expected to be rarer in our detected sample, due to their intrinsic rareness and the observational selection effects. The increased precision in the measurement of the maximum black hole mass can help in better constraining the upper mass gap in the mass spectrum.

Motivated by recent achievements of a full general relativistic method in estimating the mass-to-distance ratio of supermassive black holes hosted at the core of active galactic nuclei, we introduce the new concept redshift rapidity in order to express the Schwarzschild black hole mass and its distance from the Earth just in terms of observational quantities. The redshift rapidity is also an observable relativistic invariant that represents the evolution of the frequency shift with respect to proper time in the Schwarzschild spacetime. We extract concise and elegant analytic formulas that allow us to disentangle mass and distance to black holes in the Schwarzschild background and estimate these parameters separately. This procedure is performed in a completely general relativistic way with the aim of improving the precision in measuring cosmic distances to astrophysical compact objects. Our exact formulas are valid on the midline and close to the line of sight, having direct astrophysical applications for megamaser systems, whereas the general relations can be employed in black hole parameter estimation studies.

It is known that the appearance of microscopic objects with distinct topology and different Euler characteristics, such as instatons and wormholes, at the spacetime-foam level in Euclidean quantum gravity approaches, leads to spacetime topology changes. Such changes, in principle, may affect the field equations that arise through the semiclassical variation procedure of gravitational actions. Although in the case of Einstein-Hilbert action the presence of microscopic wormholes does not lead to any non-trivial result, when the Gauss-Bonnet term is added in the gravitational action, the above effective topological variation procedure induces an effective cosmological term that depends on the Gauss-Bonnet coupling and the wormhole density. Since the later in a dynamical spacetime is in general time-dependent, one obtains an effective dark energy sector of topological origin.

The broadband solar K-corona is linearly polarized due to Thomson scattering. Various strategies have been used to represent coronal polarization. Here, we present a new way to visualize the polarized corona, using observations from the 2023 April 20 total solar eclipse in Australia in support of the Citizen CATE 2024 project. We convert observations in the common four-polarizer orthogonal basis (0{\deg}, 45{\deg}, 90{\deg}, & 135{\deg}) to -60{\deg}, 0{\deg}, and +60{\deg} (MZP) polarization, which is homologous to R, G, B color channels. The unique image generated provides some sense of how humans might visualize polarization if we could perceive it in the same way we perceive color.

It was proposed that acoustic perturbations generated by the QCD phase transition could create an inverse turbulent cascade \cite{Kalaydzhyan:2014wca}. An assumption was that propagation toward smaller momenta could reach the wavelength of a few km, the Universe's size at the time. Such acoustic waves were proposed to be the source of gravity waves. The kilometer wavelength corresponds to year-long gravity wave today, which were recently discovered using pulsar correlations. This paper argues further that an acoustic turbulence must be an ensemble of shocks. This brings two consequences: First, shocks generate gravity waves much more efficiently than sound waves due to the intermittency of matter distribution. We reconsider gravity wave radiation, using universal emission theory for soft radiation, and argue that soft-momenta plateau should reach wavelengths of the order of shock mean free paths. Second, collisions of shock waves create local density excess, which may create primordial black holes. A good tool to settle this issue can be an evaluation of the {\em trapped surfaces} like it was done in studies of heavy ion collisions using AdS/CFT correspondence