Papers for Thursday, May 30 2024

The presence of kiloparsec-sized bubble structures in both sides of the Galactic plan suggests active phases of Sgr A$^\star$, the central supermassive black hole of the Milky-Way in the last 1-6 Myr. The contribution of such event on the cosmic-ray flux measured in the solar neighborhood is investigated with numerical simulations. We evaluate whether the population of high-energy charged particles emitted by the Galactic Center could be sufficient to significantly impact the CR flux measured in the solar neighborhood. We present a set of 3D magnetohydrodynamical simulations, following the anisotropic propagation of CR in a Milky - Way like Galaxy. Independent populations of cosmic-ray are followed through time, originating from two different sources types, namely Supernovae and the Galactic Center. To assess the evolution of the CR flux spectrum properties, we split these populations into two independent energy groups of 100 GeV and 10 TeV. We find that the anisotropic nature of cosmic-ray diffusion dramatically affects the amount of cosmic-ray energy received in the solar neighborhood. Typical timescale to observe measurable changes in the CR spectrum slope is of the order 10 Myr, largely surpassing estimated ages of the Fermi bubbles in the AGN jet-driven scenario. We conclude that a cosmic-ray outburst from the Galactic center in the last few Myr is unlikely to produce any observable feature in the local CR spectrum in the TeV regime within times consistent with current estimates of the Fermi bubbles ages.

Context. Airless planetary objects have their surfaces covered by craters, and these can be used to study the characteristics of asteroid populations. Planetary surfaces present binary craters that are associated with the synchronous impact of binary asteroids. Aims. We identify binary craters on asteroids (1) Ceres and (4) Vesta, and aim to characterize the properties (size ratio and orbital plane) of the binary asteroids that might have formed them. Methods. We used global crater databases developed in previous studies and mosaics of images from the NASA DAWN mission high-altitude and low-altitude mapping orbits. We established selection criteria to identify craters that were most likely a product of the impact of a binary asteroid. We performed numerical simulations to predict the orientation of the binary craters assuming the population of impactors has mutual orbits coplanar with heliocentric orbits, as the current census of binary asteroids suggests. We compared our simulations with our survey of binary craters on Ceres and Vesta through a Kolmogorov-Smirnov test. Results. We find geomorphological evidence of 39 and 18 synchronous impacts on the surfaces of Ceres and Vesta, respectively. The associated binary asteroids are widely separated and similar in diameter. The distributions of the orientation of these binary craters on both bodies are statistically different from numerical impact simulations that assume binary asteroids with coplanar mutual and heliocentric orbits. Conclusions. Although the identification of binary craters on both bodies and the sample size are limited, these findings are consistent with a population of well-separated and similarly sized binary asteroids with nonzero obliquity that remains to be observed, in agreement with the population of binary craters identified on Mars.

The James Webb Space Telescope (JWST) has spectroscopically confirmed numerous galaxies at $z > 10$. While weak rest-ultraviolet emission lines have only been seen in a handful of sources, the stronger rest-optical emission lines are highly diagnostic and accessible at mid-infrared wavelengths with the Mid-Infrared Instrument (MIRI) of JWST. We report the photometric detection of the most distant spectroscopically confirmed galaxy JADES-GS-z14-0 at $z = 14.32^{+0.08}_{-0.20}$ with MIRI at $7.7\ \mu\mathrm{m}$. The most plausible solution for the stellar population properties is that this galaxy contains half a billion solar masses in stars with a strong burst of star formation in the most recent few million years. For this model, at least one-third of the flux at $7.7\ \mu\mathrm{m}$ comes from the rest-optical emission lines $\mathrm{H}\beta$ and/or $\mathrm{[OIII]}\lambda\lambda4959,5007$. The inferred properties of JADES-GS-z14-0 suggest rapid mass assembly and metal enrichment during the earliest phases of galaxy formation.

The splashback radius is one popular method to constrain the size of galaxy clusters. It is typically measured through the logarithmic derivative of the galaxy number density profile, since doing so is more observationally viable and computationally inexpensive compared to other methods. However, measuring the splashback radius through the galaxy number density has consistently produced smaller values of the splashback radius than those measured with dark matter density or other processes. Dynamical friction has been posited as one possible reason that splashback radii measured through galaxy number densities are reduced, since it decays the orbits of subhaloes within the halo, however, the effects of dynamical friction cannot be isolated within cosmological simulations. Here, we present idealized simulations starting with isolated galaxy clusters drawn from the IllustrisTNG cosmological simulation, where we isolate dynamical friction. We show that although dynamical friction can reduce measurements of the splashback radius, it does not have a significant effect on clusters with $M_\mathrm{200,mean} > 10^{14} \mathrm{M_\odot}$, and thus cannot completely account for previously measured discrepancies.

Galaxy number counts probe the evolution of galaxies over cosmic time, and serve as a valuable comparison point to theoretical models of galaxy formation. We present new galaxy number counts in eight photometric bands between 5 and 25 $\mu$m from the Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES) and the JWST Advanced Deep Extragalactic Survey (JADES) deep MIRI parallel, extending to unprecedented depth. By combining our new MIRI counts with existing data from Spitzer and AKARI, we achieve counts across 3-5 orders of magnitude in flux in all MIRI bands. Our counts diverge from predictions from recent semi-analytical models of galaxy formation, likely owing to their treatment of mid-infrared aromatic features. Finally, we integrate our combined JWST-Spitzer counts at 8 and 24 $\mu$m to measure the cosmic infrared background (CIB) light at these wavelengths; our measured CIB fluxes are consistent with those from previous mid-infrared surveys, but larger than predicted by models based on TeV blazar data.

There is now abundant observational evidence that star formation is a highly dynamical process that connects filament hierarchies and supernova feedback from galaxy scale kpc filaments and superbubbles, to giant molecular clouds (GMCs) on 100 pc scales and star clusters (1 pc). Here we present galactic multi-scale MHD simulations that track the formation of structure from galactic down to sub pc scales in a magnetized, Milky Way like galaxy undergoing supernova driven feedback processes. We do this by adopting a novel zoom-in technique that follows the evolution of typical 3-kpc sub regions without cutting out the surrounding galactic environment, allowing us to reach 0.28 pc resolution in the individual zoom-in regions. We find a wide range of morphologies and hierarchical structure including superbubbles, turbulence, kpc atomic gas filaments hosting multiple GMC condensations that are often associated with superbubble compression; down to smaller scale filamentary GMCs and star cluster regions within them. Gas accretion and compression ultimately drive filaments over a critical, scale - dependent, line mass leading to gravitational instabilities that produce GMCs and clusters. In quieter regions, galactic shear can produce filamentary GMCs within flattened, rotating disk-like structures on 100 pc scales. Strikingly, our simulations demonstrate the formation of helical magnetic fields associated with the formation of these disk like structures.

Accurate chemical compositions of star-forming regions are a critical diagnostic tool to characterize the star formation history and gas flows which regulate galaxy formation. However, the abundance discrepancy factor (ADF) between measurements from the "direct" optical electron temperature ($T_e$) method and from the recombination lines (RL) represents $\sim0.2$ dex systematic uncertainty in oxygen abundance. The degree of uncertainty for other elements is unknown. We conduct a comprehensive analysis of O$^{++}$ and N$^+$ ion abundances using optical and far-infrared spectra of a star-forming region within the nearby dwarf galaxy Haro 3, which exhibits a typical ADF. Assuming homogeneous conditions, the far-IR emission indicates an O abundance which is higher than the $T_e$ method and consistent with the RL value, as would be expected from temperature fluctuations, whereas the N abundance is too large to be explained by temperature fluctuations. Instead a component of highly obscured gas is likely required to explain the high far-IR to optical flux ratios. Accounting for this obscured component reduces both the IR-based metallicities and the inferred magnitude of temperature fluctuations, such that they cannot fully explain the ADF in Haro 3. Additionally, we find potential issues when predicting the RL fluxes from current atomic data. Our findings underscore the critical importance of resolving the cause of abundance discrepancies and understanding the biases between different metallicity methods. This work represents a promising methodology, and we identify further approaches to address the current dominant uncertainties.

Dynamical interactions in globular clusters (GCs) significantly impact the formation and evolution of binary sources, including cataclysmic variables (CVs). This study investigates the connection between dynamical states of GCs and X-ray luminosity ($L_{x}$) distributions of CV populations through both simulations and actual observations. Utilizing a Monte Carlo simulation tool, MOCCA, we categorize the simulated GCs into three different evolutionary stages which are referred to as Classes I/II/III. Significant differences are found in the $L_{x}$ distributions of the CVs among these three Classes. In observational aspects, we have analyzed 179 CV candidates in 18 GCs observed by the {\it Chandra} X-ray Observatory. By dividing these GCs into three Families of different dynamical ages, namely Families I/II/III, the $L_{x}$ distributions of the CV candidates also show significant differences among these three Families. Both simulations and observational results suggest that CVs in more dynamically evolved clusters (Class/Family III) exhibit brighter X-ray emission. This highlights the influence of the dynamical status of a GC on the properties of its hosted compact binaries. Similar to blue stragglers, CV populations can serve as tracers of a GC's dynamical history. Our findings provide insights for understanding the interplay between intracluster dynamics and the evolution of compact binaries in GCs.

The discovery by JWST of an abundance of luminous galaxies in the very early Universe suggests that galaxies developed rapidly, in apparent tension with many standard models. However, most of these galaxies lack spectroscopic confirmation, so their distances and properties are uncertain. We present JADES JWST/NIRSpec spectroscopic confirmation of two luminous galaxies at redshifts of $z=14.32^{+0.08}_{-0.20}$ and $z=13.90\pm0.17$. The spectra reveal ultraviolet continua with prominent Lyman-$\alpha$ breaks but no detected emission lines. This discovery proves that luminous galaxies were already in place 300 million years after the Big Bang and are more common than what was expected before JWST. The most distant of the two galaxies is unexpectedly luminous (M$_{\rm uv}=-20.81\pm0.16$) and is spatially resolved with a radius of 260 parsecs. Considering also the steep ultraviolet slope of the second galaxy ($\beta=-2.71\pm0.19$), we conclude that both are dominated by stellar continuum emission, showing that the excess of luminous galaxies in the early Universe cannot be entirely explained by accretion onto black holes. Galaxy formation models will need to address the existence of such large and luminous galaxies so early in cosmic history.

We present high-cadence photometric and spectroscopic observations of supernova (SN) 2024ggi, a Type II SN with flash spectroscopy features which exploded in the nearby galaxy NGC 3621 at $\sim$7 Mpc. The light-curve evolution over the first 30 hours can be fit by two power law indices with a break after 22 hours, rising from $M_V \approx -12.95$ mag at +0.66 days to $M_V \approx -17.91$ mag after 7 days. In addition, the densely sampled color curve shows a strong blueward evolution over the first few days and then behaves as a normal SN II with a redward evolution as the ejecta cool. Such deviations could be due to interaction with circumstellar material (CSM). Early high- and low-resolution spectra clearly show high-ionization flash features from the first spectrum to +3.42 days after the explosion. From the high-resolution spectra, we calculate the CSM velocity to be 37 $\pm~4~\mathrm{km\,s^{-1}} $. We also see the line strength evolve rapidly from 1.22 to 1.49 days in the earliest high-resolution spectra. Comparison of the low-resolution spectra with CMFGEN models suggests that the pre-explosion mass-loss rate of SN 2024ggi falls in a range of $10^{-3}$ to $10^{-2}$ M$_{\odot}$ yr$^{-1}$, which is similar to that derived for SN 2023ixf. However, the rapid temporal evolution of the narrow lines in the spectra of SN 2024ggi ($R_\mathrm{CSM} \sim 2.7 \times 10^{14} \mathrm{cm}$) could indicate a smaller spatial extent of the CSM than in SN 2023ixf ($R_\mathrm{CSM} \sim 5.4 \times 10^{14} \mathrm{cm}$) which in turn implies lower total CSM mass for SN 2024ggi.

We present an analysis of the number density of galaxies as a function of stellar mass (i.e., the stellar mass function, SMF) in the COSMOS field at z~3.3, making a comparison between the SMF in overdense environments and the SMF in the coeval field. In particular, this region contains the Elentári proto-supercluster, a system of 6 extended overdensities spanning ~70 cMpc on a side. A clear difference is seen in the high-mass slope of these SMFs, with overdense regions showing an increase in the ratio of high-mass galaxies to low-mass galaxies relative to the field, indicating a more rapid build-up of stellar mass in overdense environments. This result qualitatively agrees with analyses of clusters at z~1, though the differences between protocluster and field SMFs at z~3.3 are smaller. While this is consistent with overdensities enhancing the evolution of their member galaxies, potentially through increased merger rates, whether this enhancement begins in protocluster environments or even earlier in group environments is still unclear. Though the measured fractions of quiescent galaxies between the field and overdense environments do not vary significantly, implying that this stellar mass enhancement is ongoing and any starbursts triggered by merger activity have not yet quenched, we note that spectroscopic observations are biased towards star-forming populations, particularly for low-mass galaxies. If mergers are indeed responsible, high resolution imaging of Elentári and similar structures at these early epochs should then reveal increased merger rates relative to the field. Larger samples of well-characterized overdensities are necessary to draw broader conclusions in these areas.

We investigate the scenario in which primordial black holes (PBHs) with masses Mpbh < 10^9 g undergo Hawking evaporation, around the Big-Bang nucleosynthesis (BBN) epoch. The evaporation process modifies the Universe's expansion rate and the baryon-to-photon ratio, leading to an alteration of the primordial abundance of light nuclei. We present numerical solutions for the set of equations describing this physics, considering different values of PBH masses and abundances at their formation, showing how their evaporation impacts the abundances of light nuclei, obtained by incorporating the non-standard Hubble rate and baryon-to-photon ratio into the BBN code PArthENoPE. The results are then used to place upper bounds for the PBH relative abundance at formation in the range 10^8 g < Mpbh < 10^9 g, providing the strongest constraints existing to-date in this mass range.

Advances in time domain astronomy have produced a growing population of flares from galactic nuclei, including both tidal disruption events (TDEs) and flares in active galactic nuclei (AGN). Because TDEs are uncommon and AGN variability is abundant, large-amplitude AGN flares are usually not categorized as TDEs. While TDEs are normally channelled by the collisional process of two-body scatterings over relaxation timescale, the quadrupole moment of a gas disk alters the stellar orbits, allowing them to collisionlessly approach the central massive black hole (MBH). This leads to an effectively enlarged loss cone, the \emph{loss wedge}. Earlier studies found a moderate enhancement, up to a factor $\sim 2-3$, of TDE rates $\dot{N}_{\rm 2b} $ for a static axisymmetric perturbation. Here we study the loss wedge problem for an evolving AGN disk, which can capture large number of stars into the growing loss wedge over much shorter times. The rates $\dot{N}_{\rm cl}$ of collisionless TDEs produced by these time-evolving disks are much higher than the collisional rates $\dot{N}_{\rm 2b}$ in a static loss wedge. We calculate the response of a stellar population to the axisymmetric potential of an adiabatically growing AGN disk and find that the highest rates of collisionless TDEs are achieved for the largest (i) MBH masses $M_{\bullet}$ and (ii) disk masses $M_{\rm d}$. For $M_{\bullet}\sim 10^7 M_\odot$ and $M_{\rm d} \sim 0.1 M_{\bullet}$, the rate enhancement can be up to a factor $\dot{N}_{\rm cl}/\dot{N}_{\rm 2b} \sim 10$. The orbits of collisionless TDEs sometimes have a preferred orientation in apses, carrying implications for observational signatures of resulting flares.

In Patil et. al 2024a, we developed a multitaper power spectrum estimation method, mtNUFFT, for analyzing time-series with quasi-regular spacing, and showed that it not only improves upon the statistical issues of the Lomb-Scargle periodogram, but also provides a factor of three speed up in some applications. In this paper, we combine mtNUFFT with the harmonic F-test to test the hypothesis that a strictly periodic signal or its harmonic (as opposed to e.g. a quasi-periodic signal) is present at a given frequency. This mtNUFFT/F-test combination shows that multitapering allows detection of periodic signals and precise estimation of their frequencies, thereby improving both power spectrum estimation and harmonic analysis. Using asteroseismic time-series data for the Kepler-91 red giant, we show that the F-test automatically picks up the harmonics of its transiting exoplanet as well as certain dipole ($l=1$) mixed modes. We use this example to highlight that we can distinguish between different types of stellar oscillations, e.g., transient (damped, stochastically-excited) and strictly periodic (undamped, heat-driven). We also illustrate the technique of dividing a time-series into chunks to further examine the transient versus periodic nature of stellar oscillations. The harmonic F-test combined with mtNUFFT is implemented in the public Python package tapify (this https URL), which opens opportunities to perform detailed investigations of periodic signals in time-domain astronomy.

The dark matter subhalos orbiting in a galactic halo perturb the orbits of stars in thin stellar streams. Over time the random velocities in the streams develop non-Gaussian wings. The rate of velocity increase is approximately a random walk at a rate proportional to the number of subhalos, primarily those in the mass range $\approx 10^{6-7} M_\odot$. The distribution of random velocities in long, thin, streams is measured in simulated Milky Way-like halos that develop in representative WDM and CDM cosmologies. The radial velocity distributions are well modeled as the sum of a Gaussian and an exponential. The resulting MCMC fits find Gaussian cores of 1-2 km/sec and exponential wings that increase from 3 km/sec for 5.5 keV WDM, 4 km/sec for 7 keV WDM, to 6 km/sec for a CDM halo. The observational prospects to use stream measurements to constrain the nature of galactic dark matter are discussed.

The Single Aperture Large Telescope for Universe Studies (SALTUS) is a mission concept for a far-infrared observatory developed under the recent Astrophysics Probe Explorer opportunity from NASA. The enabling element of the program is a 14 m diameter inflatable primary mirror, M1. Due to its importance to SALTUS and potentially other space observatories, this paper focuses entirely on M1. We present a historical overview of inflatable systems, illustrating that M1 is the logical next step in the evolution of such systems. The process of design and manufacture is addressed. We examine how M1 performs in its environment in terms of operating temperature, interaction with the solar wind, and shape change due to non-penetrating particles. We investigate the longevity of the inflatant in detail and show it meets mission lifetime requirements with ample margin and discuss the development and testing to realize the flight M1.

High-number-density tracers of large-scale structure, such as the HI-rich galaxies measured by 21 cm intensity mapping, have low sampling noise, making them particularly promising as cosmological probes. At large scales, this sampling noise can be subdominant to other scale-independent contributions to the power spectrum, arising from nonlinear bias. This has important consequences for cosmological constraints obtained from such tracers, since it indicates that using the power spectrum does not lead to optimal constraints even in the linear regime. In this paper, we provide a conservative estimate of the possible improvement in constraining power of a 21cm survey if one were to use an optimal analysis strategy (such as field-level analysis), where only the true sampling noise enters the error budget. We find that improvements in uncertainties on some cosmological parameters can be as large as 50%, depending on redshift, foreground cleaning efficiency, scales used in the analysis, and instrumental noise. One byproduct of our work is measurements of bias parameters and stochasticity for neutral hydrogen in the IllustrisTNG simulation over a wide range of redshifts; we provide simple fitting formulas for these measurements. Our results motivate further exploration of new optimal analysis techniques and provide important insights into the constraining power of current and future 21 cm surveys.

Modern high-sensitivity radio telescopes are discovering an increased number of resolved sources with intricate radio structures and fainter radio emissions. These sources often present a challenge because source detectors might identify them as separate radio sources rather than components belonging to the same physically connected radio source. Currently, there are no reliable automatic methods to determine which radio components are single radio sources or part of multi-component sources. We propose a deep learning classifier to identify those sources that are part of a multi-component system and require component association on data from the LOFAR Two-Metre Sky Survey (LoTSS). We combine different types of input data using multi-modal deep learning to extract spatial and local information about the radio source components: a convolutional neural network component that processes radio images is combined with a neural network component that uses parameters measured from the radio sources and their nearest neighbours. Our model retrieves 94 per cent of the sources with multiple components on a balanced test set with 2,683 sources and achieves almost 97 per cent accuracy in the real imbalanced data (323,103 sources). The approach holds potential for integration into pipelines for automatic radio component association and cross-identification. Our work demonstrates how deep learning can be used to integrate different types of data and create an effective solution for managing modern radio surveys.

Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existence of a fainter, undetected population. We present a systematic search in the ZTF archive of alerts from 1 June 2019 to 1 September 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. Candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. While superluminous supernovae (SLSNe) cannot be fully rejected, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 $\pm$ 3 and 34 $\pm$ 1 days for the two events, respectively, which suggests that we have uncovered a population with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, like the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

Magnetic holes (MHs) are coherent magnetic field dips whose size ranges from fluid to kinetic scale, ubiquitously observed in the heliosphere and in planetary environments. Despite the longstanding effort in interpreting the abundance of observations, the origin and properties of MHs are still debated. In this letter, we investigate the interplay between plasma turbulence and MHs, using a 2D hybrid simulation initialized with solar wind parameters. We show that fully developed turbulence exhibits localized elongated magnetic depressions, whose properties are consistent with linear MHs frequently encountered in space. The observed MHs develop self-consistently from the initial magnetic field perturbations, by trapping hot ions with large pitch angles. Ion trapping produces an enhanced perpendicular temperature anysotropy that makes MHs stable for hundreds of ion gyroperiods, despite the surrounding turbulence. We introduce a new quantity, based on local magnetic field and ion temperature values, to measure the efficiency of ion trapping, with potential applications to the detection of MHs in satellite measurements. We complement this method by analyzing the ion velocity distribution functions inside MHs. Our diagnostics reveal the presence of trapped gyrotropic ion populations, whose velocity distribution is consistent with a loss cone, as expected for the motion of particles inside a magnetic mirror. Our results have potential implications for the theoretical and numerical modelling of MHs.

Jupiter's moon Io hosts extensive volcanism driven by tidal heating. The isotopic composition of Io's inventory of volatile elements, including sulfur and chlorine, reflects its outgassing and mass loss history and provides an avenue for exploring its evolution. We used millimeter observations of Io's atmosphere to measure sulfur isotopes in gaseous SO2 and SO, and chlorine isotopes in gaseous NaCl and KCl. We find $^{34}$S/$^{32}$S=0.0595$\pm$0.0038 ($\delta^{34}$S=+347$\pm$86 per mille), which is highly enriched compared to average Solar System values and indicates that Io has lost 94 to 99% of its available sulfur. Our measurement of $^{37}$Cl/$^{35}$Cl=0.403$\pm$0.028 ($\delta^{37}$Cl=+263$\pm$88 per mille) shows chlorine is similarly enriched. These measurements indicate that Io has been volcanically active for most or all of its history, with potentially higher outgassing and mass-loss rates at earlier times.

Reanalyzing contact binaries with space-based photometric data and investigating possible parameter changes can yield accurate samples for theoretical studies. We investigated light curve solutions and fundamental parameters for twenty contact binary systems. The most recent Transiting Exoplanet Survey Satellite (TESS) data is used to analyze. The target systems in the investigation have an orbital period of less than 0.58 days. Light curve solutions were performed using the PHysics Of Eclipsing BinariEs (PHOEBE) Python code version 2.4.9. The results show that systems had various mass ratios from q=0.149 to q=3.915, fillout factors from f=0.072 to f=0.566, and inclinations from i=52.8deg to i=87.3deg. The effective temperature of the stars was less than 7016 K, which was expected given the features of most contact binary stars. Twelve of the target systems' light curves were asymmetrical in the maxima, showing the O'Connell effect, and a starspot was required for light curve solutions. The estimation of the absolute parameters of the binary systems was presented using the a-P empirical relationship and discussed. The orbital angular momentum (J_0) of the systems was calculated. The positions of the systems were also depicted on the M-L, M-R, q-L_ratio, M_tot-J_0, and T-M diagrams.

Context. The YORP effect is the thermal torque generated by radiation from the surface of an asteroid. The effect is sensitive to surface topology, including small-scale roughness, boulders, and craters. Aims: The aim of this paper is to develop a computationally efficient semi-analytical model for the crater-induced YORP (CYORP) effect that can be used to investigate the functional dependence of this effect. Methods. This study obtains the temperature field in a crater over a rotational period in the form of a Fourier series, accounting for the effects of self-sheltering, self-radiation, and self-scattering. Results. We obtain the temperature field of a crater, accounting for the thermal inertia, crater shape, and crater location. We then find that the CYORP effect is negligible when the depth-to-diameter ratio is smaller than 0.05. In this case, it is reasonable to assume a convex shape for YORP calculations. Varying the thermal conductivity yields a consistent value of approximately 0.01 for the spin component of the CYORP coefficient, while the obliquity component is inversely related to thermal inertia, declining from 0.004 in basalt to 0.001 in metal. For a z-axis symmetric shape, the CYORP spin component vanishes, while the obliquity component persists. Our model confirms that the total YORP torque is damped by a few tens of percent by uniformly distributed small-scale surface roughness. Furthermore, for the first time, we calculate the change in the YORP torque at each impact on the surface of an asteroid explicitly and compute the resulting stochastic spin evolution more precisely. Conclusions. The semi-analytical method that we developed, which benefits from fast computation, offers new perspectives for future investigations of the YORP modeling of real asteroids and for the complete rotational and orbital evolution of asteroids accounting for collisions.

Solar magnetic fields are closely related to various physical phenomena on the sun, which can be extrapolated with different models from photospheric magnetograms. However, the Open Flux Problem (OFP), the underestimation of the magnetic field derived from the extrapolated model, is still unsolved. To minimize the impact of the OFP, we propose three evaluation parameters to quantitatively evaluate magnetic field models and determine the optimal free parameters in the models by constraining the coronal magnetic fields (CMFs) and the interplanetary magnetic fields (IMFs) with real observations. Although the OFP still exists, we find that magnetic field lines traced from the coronal models effectively capture the intricate topological configurations observed in the corona, including streamers and plumes. The OFP is lessened by using the HMI synoptic map instead of the GONG daily synoptic maps, and the PFSS+PFCS model instead of the CSSS model. For Carrington Rotation (CR) 2231 at the solar minimum, we suggest that the optimal parameters for the PFSS+PFCS model are $R_{\mathrm{ss}} = 2.2-2.5\ R_{\mathrm{sun}}$ and $R_{\mathrm{scs}} = 10.5-14.0\ R_{\mathrm{sun}}$, as well as for the CSSS model are $R_{\mathrm{cs}} = 2.0 - 2.4\ R_{\mathrm{sun}}$, $R_{\mathrm{ss}} = 11.0 - 14.7\ R_{\mathrm{sun}}$ and $a = 1.0\ R_{\mathrm{sun}}$. Despite the IMFs at 1 AU being consistent with the measurements by artificially increasing the polar magnetic fields, the IMFs near the sun are still underestimated. The OFP might be advanced by improving the accuracy of both the weak magnetic fields and polar magnetic fields, especially considering magnetic activities arising from interplanetary physical processes.

The existence of high-$z$ over-massive supermassive black holes represents a major conundrum in our understanding of black hole evolution. In this paper, we probe from the observational point of view how early Universe environmental conditions could have acted as an evolutionary mechanism for the accelerated growth of the first black holes. Under the assumption that the early Universe is dominated by dwarf galaxies, we investigate the hypothesis that dwarf-dwarf galaxy interactions trigger black hole accretion. We present the discovery of 82 dwarf-dwarf galaxy pairs and 11 dwarf galaxy groups using the Hubble Space Telescope, doubling existing samples. The dwarf systems span a redshift range of 0.13$<$z$<$1.5, and a stellar mass range of 7.24$<$log(M$_*$/\(M_\odot\))$<$9.73. We performed an X-ray study of a subset of these dwarf systems with Chandra and detected six new AGN, increasing the number of known dwarf-dwarf-merger-related AGN from one to seven. We then compared the frequency of these AGN in grouped/paired dwarfs to that of isolated dwarfs and found a statistically significant enhancement (4$\sigma$-6$\sigma$) in the interacting sample. This study, the first of its kind at the lowest mass scales, implies that the presence of a nearby dwarf neighbor is efficient in triggering black hole accretion. These results open new avenues for indirect studies of the emergence of the first supermassive black holes.

The magnetic field inside the sunspot umbra, as observed by the Full-disk MagnetoGraph (FMG) onboard the Advanced Space based Solar Observatory (ASO-S), was found to be experiencing a weakening. To address this issue, we employed a method developed by Xu et al. (2021) to correct the weakening in the data of 20 active regions observed by FMG during the period spanning December 29, 2022, to July 23, 2023. Research has revealed that the onset of magnetic field weakening occurs at a minimum magnetic field strength of 705 G, with the peak strength reaching up to 1931 G. We computed the change ratio (R1) of the unsigned magnetic flux within the sunspot umbra, considering measurements both before and after correction. The change ratio (R1) spans from 26% to 124%, indicating a significant increase in the unsigned magnetic flux within sunspot umbrae observed by FMG after correction. To illustrate this, we selected four active regions for comparison with data from the Helioseismic and Magnetic Imager (HMI). After correction, it is found that the unsigned magnetic flux in sunspot umbrae measured by FMG aligns more closely with that of HMI. This supports the effectiveness of the corrective method for FMG, despite imperfections, particularly at the umbra-penumbra boundary.

We develop a model of globular cluster (GC) formation within the cosmological hierarchy of structure formation. The model is rooted in the `two-phase' scenario of galaxy formation developed in Paper-I, where the fast accretion of dark matter halos at high redshift leads to the formation of self-gravitating, turbulent gas clouds that subsequently fragment into dynamically hot systems of dense sub-clouds with masses $\sim 10^6$-$10^7 M_\odot$. Here we elaborate on the formation, evolution, and fate of these sub-clouds, and show that some of the sub-clouds can be compactified via two distinctive channels into a 'supernova-free' regime to form two distinct populations of GCs. The model is simple, characterized by a small number of free parameters underpinned by physical considerations, and can be efficiently implemented into cosmological N-body simulations to generate a coherent sample of halos, galaxies, and GCs. Our model can reproduce a range of observations on GCs, including the mass function, the size-mass relation, the frequency per unit host galaxy/halo mass, the bimodal metallicity distribution, and the spatial profile. Predictions for GCs are made for both the local Universe and for redshift up to $z \approx 10$, and can be tested by upcoming observations.

The DESI collaboration have recently analyzed their first year of data, finding a preference for thawing dark energy scenarios when using parameterized equations of state for dark energy. We investigate whether this preference persists when the data is analyzed within the context of a well-studied field theory model of thawing dark energy, exponential quintessence. No preference for this model over $\Lambda$CDM is found, and both models are poorer fits to the data than the Chevallier-Polarski-Linder $w_0$--$w_a$ parameterization. We demonstrate that the worse fit is due to a lack of sharp features in the potential that results in a slowly-evolving dark energy equation of state that does not have enough freedom to simultaneously fit the combination of the supernovae, DESI, and cosmic microwave background data. Our analysis provides guidance for constructing dynamical dark energy models that are able to better accommodate the data.

The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense gas was derived at a density of $n$(H$_{2}$) = 10$^{5}$ cm$^{-3}$. The para-H$_2$CO lines reveal significantly higher temperatures than NH$_3$ (1,1)/(2,2) and FIR wavelengths. The dense clumps appear to correlate with the notable kinetic temperature. Among the four dense cores (N44, N46, N48, and N54), temperature gradients are observed on a scale of $\sim$0.1-0.3 pc. This suggests that the warm dense gas is influenced by internal star formation activity. With the exception of N54, the temperature profiles of these cores were fitted with power-law indices ranging from $-$0.3 to $-$0.5. This indicates that the warm dense gas is heated by radiation emitted from internally embedded protostar(s) and/or clusters. While there is no direct evidence supporting the idea that the dense gas is heated by shocks resulting from a past explosive event in the DR21 region, our measurements toward the DR21W1 region provide compelling evidence that the dense gas is indeed heated by shocks originating from the western DR21 flow. Higher temperatures appear to be associated with turbulence. The physical parameters of the dense gas in the DR21 filament exhibit a remarkable similarity to the results obtained in OMC-1 and N113. This may imply that the physical mechanisms governing the dynamics and thermodynamics of dense gas traced by H$_{2}$CO in diverse star formation regions may be dominated by common underlying principles despite variations in specific environmental conditions. (abbreviated)

This paper provides a comprehensive overview of the subsystems of the NIGHT instrument. NIGHT (the Near Infrared Gatherer of Helium Transits) is a narrowband, high-resolution spectrograph, marking the first dedicated survey instrument for exoplanetary atmosphere observations. Developed through a collaboration between the Observatory of Geneva and the Universite de Montreal, NIGHT aims to conduct an extensive statistical survey of helium atmospheres around 100+ exoplanets over several years. The instrument will report new detections of helium in exoplanet atmospheres and perform temporal monitoring of a subset of these. NIGHT measures absorption from the metastable helium state during exoplanet transits, observable in a triplet of lines around 1083nm. The instrument comprises a vacuum enclosure housing the spectrograph, a front end unit for fiber injection at the telescope's focal plane, and a calibration and control rack containing calibration light sources and control hardware. The spectrograph is optimized for efficiency, achieving a uniform throughput of approximately 71%. The primary disperser employs a VPH grating in a unique double-pass configuration, enabling a spectral resolution of 75,000 while maintaining high throughput. The detector is a HAWAII-1 infrared array, cooled to 85K, with the spectrograph operating at room temperature. Thanks to its relatively high throughput, NIGHT on a 2m class telescope is predicted to be as sensitive as existing instruments on 4m class telescopes. The front end unit injects starlight and sky background into two separate fibers leading to the spectrograph. It also performs near-infrared guiding and includes a mechanism for injecting calibration light. The assembly and optical alignment of NIGHT's spectrograph and front end unit are scheduled for July to September 2024, with the first light anticipated before early 2025.

Context: Discrete symmetries have found numerous applications in photonics and quantum mechanics, but remain little studied in fluid mechanics, particularly in astrophysics. Aims: We aim to show how PT and anti-PT symmetries determine the behaviour of linear perturbations in a wide class of astrophysical problems. They set the location of Exceptional Points in the parameter space and the associated transitions to instability, and are associated to the conservation of quadratic quantities that can be determined explicitly. Methods: We study several classical local problems: the gravitational instability of isothermal spheres and thin discs, the Schwarzschild instability, the Rayleigh-Bénard instability and acoustic waves in dust-gas mixtures. We calculate the locations and the order of the Exceptional Points with a method of resultant, as well as the conserved quantities in the different regions of the parameter space using Krein theory. Results: All problems studied here exhibit discrete symmetries, even though Hermiticity is broken by different physical processes (self-gravity, buoyancy, diffusion, drag). This analysis provides genuine explanations for certain instabilities, and for the existence of regions in the parameter space where waves do not propagate. Those correspond to breaking of PT and anti-PT symmetries respectively. Not all instabilities are associated to symmetry breaking (e.g. the Rayleigh-Benard instability).

We are developing an imaging polarimeter by combining a fine-pixel CMOS image sensor with a coded aperture mask as part of the cipher project, aiming to achieve X-ray polarimetry in the energy range of $10$$\unicode{x2013}$$30~\mathrm{keV}$. A successful proof-of-concept experiment was conducted using a fine-pixel CMOS sensor with a $2.5~\mathrm{\mu m}$ pixel size. In this study, we conducted beam experiments to assess the modulation factor (MF) of the CMOS sensor with a $1.5~\mathrm{\mu m}$ pixel size manufactured by Canon and to determine if there was any improvement in the MF. The measured MF was $8.32\% \pm 0.34\%$ at $10~\mathrm{keV}$ and $16.10\% \pm 0.68\%$ at $22~\mathrm{keV}$, exceeding those of the $2.5~\mathrm{\mu m}$ sensor in the $6$$\unicode{x2013}$$22~\mathrm{keV}$ range. We also evaluated the quantum efficiency of the sensor, inferring a detection layer thickness of $2.67 \pm 0.48~{\rm \mu m}$. To develop a more sensitive polarimeter, a sensor with a thicker detection layer, smaller pixel size, and reduced thermal diffusion effect is desirable.

Low-thrust trajectories play a crucial role in optimizing scientific output and cost efficiency in asteroid belt missions. Unlike high-thrust transfers, low-thrust trajectories require solving complex optimal control problems. This complexity grows exponentially with the number of asteroids visited due to orbital mechanics intricacies. In the literature, methods for approximating low-thrust transfers without full optimization have been proposed, including analytical and machine learning techniques. In this work, we propose new analytical approximations and compare their accuracy and performance to machine learning methods. While analytical approximations leverage orbit theory to estimate trajectory costs, machine learning employs a more black-box approach, utilizing neural networks to predict optimal transfers based on various attributes. We build a dataset of about 3 million transfers, found by solving the time and fuel optimal control problems, for different time of flights, which we also release open-source. Comparison between the two methods on this database reveals the superiority of machine learning, especially for longer transfers. Despite challenges such as multi revolution transfers, both approaches maintain accuracy within a few percent in the final mass errors, on a database of trajectories involving numerous asteroids. This work contributes to the efficient exploration of mission opportunities in the asteroid belt, providing insights into the strengths and limitations of different approximation strategies.

The main science aim of the BlackGEM array is to detect optical counterparts to gravitational wave mergers. Additionally, the array will perform a set of synoptic surveys to detect Local Universe transients and short time-scale variability in stars and binaries, as well as a six-filter all-sky survey down to ~22nd mag. The BlackGEM Phase-I array consists of three optical wide-field unit telescopes. Each unit uses an f/5.5 modified Dall-Kirkham (Harmer-Wynne) design with a triplet corrector lens, and a 65cm primary mirror, coupled with a 110Mpix CCD detector, that provides an instantaneous field-of-view of 2.7~square degrees, sampled at 0.564\arcsec/pixel. The total field-of-view for the array is 8.2 square degrees. Each telescope is equipped with a six-slot filter wheel containing an optimised Sloan set (BG-u, BG-g, BG-r, BG-i, BG-z) and a wider-band 440-720 nm (BG-q) filter. Each unit telescope is independent from the others. Cloud-based data processing is done in real time, and includes a transient-detection routine as well as a full-source optimal-photometry module. BlackGEM has been installed at the ESO La Silla observatory as of October 2019. After a prolonged COVID-19 hiatus, science operations started on April 1, 2023 and will run for five years. Aside from its core scientific program, BlackGEM will give rise to a multitude of additional science cases in multi-colour time-domain astronomy, to the benefit of a variety of topics in astrophysics, such as infant supernovae, luminous red novae, asteroseismology of post-main-sequence objects, (ultracompact) binary stars, and the relation between gravitational wave counterparts and other classes of transients

The divide between weak and strong lensing is of course artificial, in that both regimes are manifestations of the same physical phenomenon: gravity bending the path of light. Nevertheless, these two regimes have to a large extent been treated separately, with different methods. This review traces the development of methods combining weak-lensing and strong-lensing data for joint lens-mass reconstruction, with a particular emphasis on cluster lenses, where both effects occur. We conclude that so-called inverse methods have been successful in merging the two regimes insofar data analysis is concerned. However, a number of improvements seem to be needed. First, not many studies include weak lensing data beyond shear. In light of the unprecedented quality of the data of JWST and future surveys, this is a clear point of improvement. Especially so, since flexion terms have proven useful in determining sub-structures. Second, considering the amount of data available, and the complexity of non-parametric lenses, automating the processes of lens-mass reconstruction would be beneficial. Towards this end, invoking machine learning seems like a promising way forward. The silence of the literature on this latter point is in fact somewhat surprising.

Aims. We report the confirmation of a new transiting exoplanet orbiting the star TOI-5076. Methods. We present our vetting procedure and follow-up observations which led to the confirmation of the exoplanet TOI-5076b. In particular, we employed high-precision {\it TESS} photometry, high-angular-resolution imaging from several telescopes, and high-precision radial velocities from HARPS-N. Results. From the HARPS-N spectroscopy, we determined the spectroscopic parameters of the host star: T$\rm_{eff}$=(5070$\pm$143) K, log~g=(4.6$\pm$0.3), [Fe/H]=(+0.20$\pm$0.08), and [$\alpha$/Fe]=0.05$\pm$0.06. The transiting planet is a warm sub-Neptune with a mass m$\rm_p=$(16$\pm$2) M$\rm_{\oplus}$, a radius r$\rm_p=$(3.2$\pm$0.1)~R$\rm_{\oplus}$ yielding a density $\rho_p$=(2.8$\pm$0.5) g cm$^{-3}$. It revolves around its star approximately every 23.445 days. Conclusions. The host star is a metal-rich, K2V dwarf, located at about 82 pc from the Sun with a radius of R$_{\star}$=(0.78$\pm$0.01) R$_{\odot}$ and a mass of M$_{\star}$=(0.80$\pm$0.07) M$_{\odot}$. It forms a common proper motion pair with an M-dwarf companion star located at a projected separation of 2178 au. The chemical analysis of the host-star and the Galactic-space velocities indicate that TOI-5076 belongs to the old population of thin-to-thick-disk transition stars. The density of TOI-5076b suggests the presence of a large fraction by volume of volatiles overlying a massive core. We found that a circular orbit solution is marginally favored with respect to an eccentric orbit solution for TOI-5076b.

Low radio frequency spectral index measurements are a powerful tool to distinguish between different emission mechanisms and, in turn, to understand the nature of the sources. Besides the standard method of estimating the ``broadband" spectral index of sources from observations in two different frequency ``bands", if the observations were made with large instantaneous bandwidth, the ``in-band" spectral index can be determined, either using images of emission at multiple frequency ranges within a band or using the novel Multi Term-Multi Frequency Synthesis (MT-MFS) imaging algorithm. Here, using simulated upgraded Giant Metrewave Radio Telescope (uGMRT) data, we have systematically studied the reliability of various methods of spectral index estimation for sources with a wide range of signal-to-noise ratio (SNR). It is found that, for synthetic uGMRT point source data, the MT-MFS imaging algorithm produces in-band spectral indices for SNR~$\lesssim100$ that have errors $\gtrsim 0.2$, making them unreliable. However, at a similar SNR, the sub-band splitting method produces errors $\lesssim 0.2$, which are more accurate and unbiased in-band spectral indices. The broadband spectral indices produce errors $\lesssim 0.2$ even for SNR $\gtrsim 15$, and hence, they are most reliable if there are no higher-order variations in the spectral index. These results may be used to improve the uGMRT observation and data analysis strategies depending on the brightness of the target source.

The Moons And Jupiter Imaging Spectrometer (MAJIS) is the visible and near-infrared imaging spectrometer onboard ESA s Jupiter Icy Moons Explorer (JUICE) mission. Before its integration into the spacecraft, the instrument undergoes an extensive ground calibration to establish its baseline performances. This process prepares the imaging spectrometer for flight operations by characterizing the behavior of the instrument under various operative conditions and uncovering instrumental distortions that may depend on instrumental commands. Two steps of the on-ground calibration campaigns were held at the instrument level to produce the data. Additional in-flight measurements have recently been obtained after launch during the Near-Earth Commissioning Phase. In this article, we present the analyses of these datasets, focusing on the characterization of the spectral performances. First, we describe and analyze the spectral calibration datasets obtained using both monochromatic sources and polychromatic sources coupled with solid and gas samples. Then, we derive the spectral sampling and the spectral response function over the entire field of view. These spectral characteristics are quantified for various operational parameters of MAJIS, such as temperature and spectral binning. The derived on-ground performances are then compared with in-flight measurements obtained after launch and presented in the framework of the MAJIS performance requirements.

Several works have recently applied Jeans modelling to Gaia-based datasets to infer the circular velocity curve for the Milky Way. Such works have consistently found evidence for a continuous decline in the rotation curve beyond $\sim$15kpc possibly indicative of a light dark matter halo. We used Gaia DR3 RVS data, supplemented with Bayesian distances to determine the radial variation of the second moments of the velocity distribution for stars close to the Galactic plane. We have used these profiles to determine the rotation curve using the Jeans equations under the assumption of axisymmetry and explored how they vary with azimuth and above and below the Galactic disk plane. We have applied the same methodology to an N-body simulation of a Milky Way-like galaxy impacted by a satellite akin the Sagittarius dwarf and to the Auriga suite of cosmological simulations. We reveal evidence of disequilibrium and deviations from axisymmetry closer in. We find that the second moment of $V_R$ flattens out at $R \gtrsim 12.5$kpc, and that the second moment of $V_{\phi}$ is different above and below the plane for $R \gtrsim 11$kpc. The simulations indicate that these features are typical of galaxies that have been perturbed by external satellites. They also suggest that the difference between the true circular velocity curve and that inferred from Jeans equations can be as high as 15$\%$, but is likely of order 10$\%$ for the Milky Way. This is of larger amplitude than the systematics inherent to Jeans equations. However, if the density of the tracer population were truncated at large radii, the erroneous conclusion of a steeply declining rotation curve can be reached. We find that steady-state axisymmetric Jeans modelling becomes less robust at large radii, indicating that particular caution is needed when interpreting the rotation curve inferred in those regions.

The accurate characterisation of the stellar magnetism of planetary host stars has been gaining momentum, especially in the context of transmission spectroscopy investigations of exoplanets. Indeed, the magnetic field regulates the amount of energetic radiation and stellar wind impinging on planets, as well as the presence of inhomogeneities on the stellar surface that hinder the precise extraction of the planetary atmospheric absorption signal. We initiated a spectropolarimetric campaign to unveil the magnetic field properties of known exoplanet hosting stars included in the current list of potential Ariel targets. In this work, we focus on HD 63433, a young solar-like star hosting two sub-Neptunes and an Earth-sized planet. These exoplanets orbit within 0.15 au from the host star and have likely experienced different atmospheric evolutionary paths. We analysed optical spectropolarimetric data collected with ESPaDOnS, HARPSpol, and Neo-Narval to compute the magnetic activity indices (log R'_HK , H$\alpha$, and Ca ii infrared triplet), measure the longitudinal magnetic field, and reconstruct the large-scale magnetic topology via Zeeman-Doppler imaging (ZDI). The magnetic field map was then employed to simulate the space environment in which the exoplanets orbit. The reconstructed stellar magnetic field has an average strength of 24 G and it features a complex topology with a dominant toroidal component, in agreement with other stars of a similar spectral type and age. Our simulations of the stellar environment locate 10% of the innermost planetary orbit inside the Alfvén surface and, thus, brief magnetic connections between the planet and the star can occur. The outer planets are outside the Alfvén surface and a bow shock between the stellar wind and the planetary magnetosphere could potentially form.

Making use of the APOGEE DR17 catalogue with high quality data for 143,509 red giant branch stars we explore the strength of different mechanisms that causes a star to radially migrate in the Milky Way stellar disk. At any position in the disk we find stars that are more metal-rich than the local interstellar medium. This is surprising and normally attributed to the migration of these stars after their formation inside their current Galactocentric-radius. Such stars are prime candidates for studying the strength of different migratory processes. We specifically select two types of metal-rich stars: i) super metal-rich stars ([Fe/H] > 0.2) and ii) stars that are more metal-rich than their local environment. For both, we explore the distribution of orbital parameters and ages as evidence of their migration history. We find that most super metal-rich stars have experienced some amount of churning as they have orbits with Rg >= 5 kpc. Furthermore, about half of the super metal-rich stars are on non-circular orbits (ecc > 0.15) and therefore also have experienced blurring. The metallicity of young stars in our sample is generally the same as the metallicity of the interstellar medium, suggesting they have not radially migrated yet. Stars with lower metallicity than the local environment have intermediate to old ages. We further find that super metal-rich stars have approximately the same age distribution at all Galactocentric-radii, which suggests that radial migration is a key mechanism responsible for the chemical compositions of stellar populations in the Milky Way.

In 2007 we were part of a team that discovered the so-called ``Lorimer Burst'', the first example of a new class of objects now known as fast radio bursts (FRBs). These enigmatic events are only a few ms in duration and occur at random locations on the sky at a rate of a few thousand per day. Several thousand FRBs are currently known. While it is now well established that they have a cosmological origin, and about 10\% of all currently known sources have been seen to exhibit multiple bursts, the origins of these enigmatic sources are currently poorly understood. In this article, we review the discovery of FRBs and present some of the highlights from the vast body of work by an international community. Following a brief overview of the scale of the visible Universe in §1, we describe the key moments in radio astronomy (§2) that led up to the discovery of the Lorimer burst (§3). Early efforts to find more FRBs are described in §4 which led to the discovery of the first repeating source (§5). In §6, as we close out on the second decade of FRBs, we outline some of the many open questions in the field and look ahead to the coming years where many surprises are surely in store.

The habitability of exoplanets hosted by M-dwarf stars dramatically depends on their space weather. We present 3D magneto-hydrodynamic simulations to characterise the magneto-plasma environment and thus the habitability of the Earth-like planet Proxima b when it is subject to both calm and extreme (CME-like) space weather conditions. We study the role of the stellar wind and planetary magnetic field, and determine the radio emission arising from the interaction between the stellar wind of Proxima and the magnetosphere of its planet Proxima b. We find that if Prox b has a magnetic field similar to that of the Earth ($B_{\rm p} = B_\oplus \approx 0.32$ G) or larger, the magnetopause standoff distance is large enough to shield the surface from the stellar wind for essentially any planetary tilt but the most extreme values (close to $90^{\circ} $), under a calm space weather. Even if Proxima b is subject to more extreme space weather conditions, the planet is well shielded by an Earth-like magnetosphere ($B_{\rm p} \approx B_\oplus$; $ \approx 23.5^{\circ}$), or if it has tilt smaller than that of the Earth. For calm space weather conditions, the radio emission caused by the day-side reconnection regions can be as high as 7$\times10^{19}$ erg s$^{-1}$ in the super-Alfvénic regime, and is on average almost an order of magnitude larger than the radio emission in the sub-Alfvénic cases, due to the much larger contribution of the bow shock. We also find that the energy dissipation at the bow shock is independent of the angle between the planet's magnetic dipole and the incident stellar wind flow. If Prox b is subject to extreme space weather conditions, the radio emission is more than two orders of magnitude larger than under calm space weather conditions. This result yields expectations for a direct detection--from Earth--in radio of giant planets in close-in orbits.

This is a short account, based on a talk given at the 2024 Moriond Cosmology Conference, of where and why string theory matters in early universe cosmology. It is written for a cosmology audience predisposed to be at best sceptical, and at worst contemptuous, of the notion that either quantum gravity or string theory has any relevance for their discipline. I cover inflation, CMB tensor modes, extended kination or tracker epochs and reheating.

In almost any study involving optical/NIR photometry, understanding the completeness of detection and recovery is an essential part of the work. The recovery fraction is, in general, a function of several variables including magnitude, color, background sky noise, and crowding. We explore how completeness can be modelled, {with the use of artificial-star tests,} in a way that includes all of these parameters \emph{simultaneously} within a neural network (NN) framework. The method is able to manage common issues including asymmetric completeness functions and the bilinear dependence of the detection limit on color index. We test the method with two sample HST (Hubble Space Telescope) datasets: the first involves photometry of the star cluster population around the giant Perseus galaxy NGC 1275, and the second involves the halo-star population in the nearby elliptical galaxy NGC 3377. The NN-based method achieves a classification accuracy of $>$\,94\%, and produces results entirely consistent with more traditional techniques for determining completeness. Additional advantages of the method are that none of the issues arising from binning of the data are present, and that a recovery probability can be assigned to every individual star in the real photometry. Our data, models, and code (called COINTOSS) can be found online on Zenodo at the following link: this https URL.

We investigate the Giant Radio Galaxies' (GRGs)-some of the largest structures powered by supermassive black holes-prevalence within supercluster environments and the influence of such environments on their properties. Utilising two large catalogues of superclusters (401) and GRGs (1446), we establish the existence of 77 GRGs (5.3%) residing in 64 superclusters (16%) within $\rm 0.05 \leq z \leq 0.42$. Among the 77 GRGs found in superclusters, we identify $\sim$70% as residing within galaxy clusters. Within the subset of GRGs not located in superclusters, which constitutes 94.7% of the sample, a mere 21% are associated with galaxy clusters, while the remaining majority are situated in sparser environments. We examine the influence of differing environments -such as cluster versus non-cluster and supercluster versus non-supercluster regions -on the size of GRGs while also exploring the driving factors behind their overall growth. Our findings show that the largest GRGs ($\gtrsim$3 Mpc) grow in underdense environments beyond the confines of dense environments. Moreover, we show that $\sim$ 24% of 1446 GRGs reside in galaxy clusters. We conclude that GRGs preferentially grow in sparser regions of the cosmic web and have a significantly larger median size. Finally, we demonstrate the potential of GRGs as astrophysical probes with specific cases where GRGs, exhibiting polarised emissions and located behind superclusters (acting as natural Faraday screens), were used to estimate magnetic field strengths of the supercluster environment at sub-micro Gauss levels.

Aeronomy, the study of Earth's upper atmosphere and its interaction with the local space environment, has long traced changes in the thermospheres of Earth and other solar system planets to solar variability in the X-ray and extreme ultraviolet (collectively, "XUV") bands. Extending comparative aeronomy to the short-period extrasolar planets may illuminate whether stellar XUV irradiation powers atmospheric outflows that change planetary radii on astronomical timescales. In recent years, near-infrared transit spectroscopy of metastable HeI has been a prolific tracer of high-altitude planetary gas. We present a case study of exoplanet aeronomy using metastable HeI transit observations from Palomar/WIRC and follow-up high-energy data from the Neil Gehrels Swift Observatory that were taken within one month of the WASP-69 system, a K-type main sequence star with a well-studied hot Jupiter companion. Supplemented by archival data, we find that WASP-69's X-ray flux in 2023 was less than 50% of what was recorded in 2016 and that the metastable HeI absorption from WASP-69b was lower in 2023 versus past epochs from 2017-2019. Via atmospheric modeling, we show that this time-variable metastable HeI signal is in the expected direction given the observed change in stellar XUV, possibly stemming from WASP-69's magnetic activity cycle. Our results underscore the ability of multi-epoch, multi-wavelength observations to paint a cohesive picture of the interaction between an exoplanet's atmosphere and its host star.

In this work, we implement Gaussian process regression to reconstruct the expansion history of the universe in a model-agnostic manner, using the Pantheon-Plus SN-Ia compilation in combination with two different BAO measurements (SDSS-IV and DESI DR1). In both the reconstructions, the $\Lambda$CDM model is always included in the 95\% confidence intervals. We find evidence that the DESI LRG data at $z_{\text{eff}} = 0.51$ is not an outlier within our model-independent framework. We study the $\mathcal{O}m$-diagnostics and the evolution of the total equation of state (EoS) of our universe, which hint towards the possibility of a quintessence-like dark energy scenario with a very slowly varying EoS, and a phantom-crossing in higher $z$. The entire exercise is later complemented by considering two more SN-Ia compilations - DES-5YR and Union3 - in combination with DESI BAO. Reconstruction with the DESI BAO + DES-5YR SN data sets predicts that the $\Lambda$CDM model lies outside the 3$\sigma$ confidence levels, whereas with DESI BAO + Union3 data, the $\Lambda$CDM model is always included within 1$\sigma$. We also report constraints on $H_0 r_d$ from our model-agnostic analysis, independent of the pre-recombination physics. Our results point towards an $\approx$ 2$\sigma$ discrepancy between the DESI + Pantheon-Plus and DESI + DES-5YR data sets, which calls for further investigation.

Stars with about 45 to 80% the mass of the Sun, so-called K dwarf stars, have previously been proposed as optimal host stars in the search for habitable extrasolar worlds. These stars are abundant, have stable luminosities over billions of years longer than Sun-like stars, and offer favorable space environmental conditions. So far, the theoretical and experimental focus on exoplanet habitability has been on even less massive, though potentially less hospitable red dwarf stars. Here we present the first experimental data on the responses of photosynthetic organisms to a simulated K dwarf spectrum. We find that garden cress Lepidium sativum under K-dwarf radiation exhibits comparable growth and photosynthetic efficiency as under solar illumination on Earth. The cyanobacterium Chroococcidiopsis sp. CCMEE 029 exhibits significantly higher photosynthetic efficiency and culture growth under K dwarf radiation compared to solar conditions. Our findings of the affirmative responses of these two photosynthetic organisms to K dwarf radiation suggest that exoplanets in the habitable zones around such stars deserve high priority in the search for extrasolar life.

WST - Widefield Spectroscopic Telescope: We summarise the design challenges of instrumentation for a proposed 12m class Telescope that aims to provide a large (>2.5 square degree) field of view and enable simultaneous Multi-object (> 20,000 objects) and Integral Field spectroscopy (inner 3x3 arcminutes field of view), initially at visible wavelengths. For the MOS mode, instrumentation includes the fiber positioning units, fiber runs and the high (R~40,000) and low (R~3,000 - 4,000) resolution spectrographs. For the MUSE like Integral Field Spectrograph, this includes the relay from the Telescope Focal Plane, the multi-stage splitting and slicing and almost 150 identical spectrographs. We highlight the challenge of mass production at a credible cost and the issues of maintenance and sustainable operation.

The intermediate-mass-ratio inspirals (IMRIs) may be surrounded by dark matter (DM) minispikes. The dynamical friction from these DM minispike structures can affect the dynamics and the gravitational wave (GW) emission of the IMRIs. We analyze the effects of general dynamical friction, with a particular contribution from DM particles moving faster than the stellar-mass black hole in an eccentric IMRI. The results show that the dynamical friction caused by these DM particles tends to eccentricify the orbit, and the general dynamical friction is able to increase the eccentricity. We also analyze the effects of general dynamical friction on the GW characteristic strain. The results indicate that the peak value of the characteristic strain occurs at higher frequencies as the power law index of DM minispike $\gamma_\mathrm{sp}$ increases. For the first time, a general analytical relation between the frequency peak value of characteristic strain of GWs and $\gamma_\mathrm{sp}$ is established. Using the analytical relation, the presence of DM and its halo density may be determined potentially from future GW data.

We study the effects of dark matter on the structural properties of neutron stars. In particular we investigate how the presence of a dark matter component influences the mass-radius relation, the value of the maximum mass of a neutron star and others stellar properties. To model ordinary matter we use a state-of-the-art equation of state of $\beta$-stable nuclear matter obtained using the Brueckner-Hartree-Fock quantum many-body approach starting from two-body and three-body nuclear interactions derived from chiral effective field theory. The dark matter component of the star is modeled as a non-self-annihilating system of spin $1/2$ fermions and its equation of state as an ideal relativistic Fermi gas. The equilibrium configurations of these dark matter admixed neutron stars (DANS) are calculated by solving a generalization of the Tolman-Oppenheimer-Volkoff equations to the case where the system consists of two perfect fluids interacting solely through gravity. We find that, depending on the dark matter particle mass $m_\chi$, one can have somehow opposite effects on the stellar properties. In the case $m_\chi = 1\, \mathrm{GeV}$, the stellar gravitational maximum mass $M_{max}$ decreases, whereas in the case $m_\chi = 0.1\, \mathrm{GeV}$, $M_{max}$ increases with respect to the maximum mass of ordinary neutron stars. We also show that the presence of dark matter has indirect sizeable effect on the proton fraction in the ordinary matter fluid and, in the case $m_\chi = 1\, \mathrm{GeV}$, results in a decrease of the threshold gravitational mass $M_{tot}^{durca}$ for having direct URCA processes and fast stellar cooling. Finally we study the stability of dark matter admixed neutron stars with respect to radial perturbations.

Between 2001 and 2023, we obtained high spectral resolution mid-infrared observations of Io using the TEXES instrument at NASA's Infrared Telescope Facility. These observations were centered at 529.8 cm-1 (18.88 {\mu}m) and include several SO2 absorption lines. By modeling the shapes and strengths of these absorption lines, we are able to determine how Io's SO2 atmospheric density varies over the 22-year time period, covering nearly two Jovian years. Previous analysis has shown that the density of Io's atmosphere on the anti-Jovian hemisphere exhibits clear seasonal temporal variability, which can be modeled as the sum of a seasonally-varying frost sublimation component and a constant component, assumed to be volcanic. The new data show that the seasonal pattern repeats during the second Jovian year, confirming the importance of sublimation support. The considerable longitudinal variability in Io's atmospheric density found in previous work is also stable over the second Jovian year with the SO2 column density on the Jupiter-facing hemisphere being 5--8 times lower than the anti-Jovian hemisphere. For the first time, we detect seasonal variability on the Jupiter-facing hemisphere as well. This can also be modeled as a combination of sublimation and a small constant source. The lower atmospheric density on the Jupiter-facing hemisphere can plausibly be explained by the daily Jupiter eclipses, which decrease the surface temperature and therefore reduce the sublimation-driven component of the atmosphere, combined with a lower level of volcanic activity directly emitting SO2 into the atmosphere.

After primordial inflation, the universe may have experienced a prolonged reheating epoch, potentially leading to a phase of matter domination supported by the oscillating inflaton field. During such an epoch, perturbations in the inflaton virialize upon reentering the cosmological horizon, forming inflaton structures. If the primordial overdensities are sufficiently large, these structures collapse to form primordial black holes (PBHs). To occur at a significant rate, this process requires an enhanced primordial power spectrum (PPS) at small scales. The enhancement of the PPS, as well as the formation and tidal interaction of the primordial structures, will in turn source a stochastic gravitational wave background(SGWB) that could be detected by current and/or future gravitational wave detectors. In this paper, we study the SGWB arising from these different sources during slow-reheating, focusing on a PPS that satisfies the requirements necessary for the formation of PBHs with a mass of $M_{\rm PBH}\simeq 10^{21}$ and that constitute the entirety of dark matter in the universe.

Nuclear Star Clusters (NSCs) are amongst the densest stellar systems in the Universe and are found at the centres of many bright spiral and elliptical galaxies, and up to ${\sim}$40% of dwarf galaxies. However, their formation mechanisms, and possible links to globular clusters (GCs), remain debated. This paper uses the EDGE simulations - a collection of zoom-in, cosmological simulations of isolated dwarf galaxies -- to present a new formation mechanism for NSCs. We find that, at a gas spatial and mass resolution of ${\sim}3\,$pc and ${\sim}161$ M$_\odot$, respectively, NSCs naturally emerge in a subset of our EDGE dwarfs with redshift-zero halo masses of $\rm{M}_{\rm{r}200\rm{c}} \sim 5 \times 10^9$ M$_\odot$. These dwarfs are quenched by reionisation, but retain a significant reservoir of gas that is unable to cool and form stars. Sometime after reionisation, the dwarfs then undergo a major (${\sim}$1:1) merger that excites rapid gas cooling, leading to a significant starburst. An NSC forms in this starburst that then quenches star formation thereafter. The result is a nucleated dwarf that has two stellar populations with distinct age: one pre-reionisation and one post-reionisation. Our mechanism is unique for two key reasons. Firstly, the low mass of the host dwarf means that NSCs, formed in this way, can accrete onto galaxies of almost all masses, potentially seeding the formation of NSCs everywhere. Secondly, our model predicts that NSCs should have at least two stellar populations with a large ($\gtrsim$1 billion year) age separation. This yields a predicted colour magnitude diagram for our nucleated dwarfs that has two distinct main sequence turnoffs. Several GCs orbiting the Milky Way, including Omega Centauri and M54, show exactly this behaviour, suggesting that they may, in fact, be accreted NSCs.

We present a computationally efficient galaxy archetype-based redshift estimation and spectral classification method for the Dark Energy Survey Instrument (DESI) survey. The DESI survey currently relies on a redshift fitter and spectral classifier using a linear combination of PCA-derived templates, which is very efficient in processing large volumes of DESI spectra within a short time frame. However, this method occasionally yields unphysical model fits for galaxies and fails to adequately absorb calibration errors that may still be occasionally visible in the reduced spectra. Our proposed approach improves upon this existing method by refitting the spectra with carefully generated physical galaxy archetypes combined with additional terms designed to absorb data reduction defects and provide more physical models to the DESI spectra. We test our method on an extensive dataset derived from the survey validation (SV) and Year 1 (Y1) data of DESI. Our findings indicate that the new method delivers marginally better redshift success for SV tiles while reducing catastrophic redshift failure by $10-30\%$. At the same time, results from millions of targets from the main survey show that our model has relatively higher redshift success and purity rates ($0.5-0.8\%$ higher) for galaxy targets while having similar success for QSOs. These improvements also demonstrate that the main DESI redshift pipeline is generally robust. Additionally, it reduces the false positive redshift estimation by $5-40\%$ for sky fibers. We also discuss the generic nature of our method and how it can be extended to other large spectroscopic surveys, along with possible future improvements.

The assembly and co-evolution of supermassive black holes (SMBH) and their host galaxy stellar population is a key open questions in galaxy evolution. Stellar mass ($M_\star$) and star formation rate (SFR), are inferred by modeling the spectral energy distribution (SED). For galaxies triggering SMBH activity, the active galactic nucleus (AGN) contaminates the light at all wavelengths, hampering the inference of galaxy parameters. Incomplete AGN templates can lead to systematic overestimates of the stellar mass, biasing our understanding of AGN-galaxy co-evolution. This challenge has gained further impetus with the advent of sensitive wide-area surveys with millions of luminous AGN, including by eROSITA, Euclid and LSST. We aim to estimate the accuracy and bias of AGN host galaxy parameters and improve upon existing techniques. This work makes two contributions: 1) a new SED fitting code, GRAHSP, with a flexible, empirically motivated AGN model including a power law continuum emission lines, a FeII forest and a flexible infrared torus. We verify that our model reproduces published X-ray to infrared SEDs of AGN to better than 20\% accuracy. A fully Bayesian fit with nested sampling includes uncertainties in the model and the data, making the inference highly robust. 2) we created a benchmark photometric dataset where pure quasars are merged with non-AGN pure galaxies into a hybrid (Chimera) object but with known galaxy and AGN properties. Comparing the true and retrieved $M_\star$, SFR and AGN luminosities shows that previous codes systematically over-estimate $M_\star$ and SFR by 0.5 dex with a wide scatter of 0.7 dex, at AGN luminosities above 10^44 erg/s. In contrast, GRAHSP shows no bias on $M_\star$ and SFR. GRAHSP also estimates more realistic uncertainties. GRAHSP enables characterization of the environmental conditions conducive to black hole growth. (abridged)

Tornadoes are severe weather phenomena characterized by a violently rotating column of air connecting the ground to a parent storm. Within the United States, hundreds of tornadoes occur every year. Despite this, the dynamics of tornado formation and propagation are not particularly well understood, in part due to the challenge of instrumentation: many existing instruments for measuring atmospheric properties are in-situ detectors, making deployment in or near an active or developing tornado difficult. Here, we combine local atmospheric and cosmic ray air shower simulation to explore the potential for remote measurement of the pressure field within tornado-producing supercell thunderstorms by examining directional variations of the atmospheric muon flux.

We present the discovery of X-ray polarization from the neutron star low-mass X-ray binary and Z-source, GX~340$+$0, using an Imaging X-ray Polarimetry Explorer (IXPE) observation in March 2024. Along with the IXPE observation, we conducted an extensive X-ray and radio monitoring campaign to ascertain the source properties during and around the IXPE observation. The source was within the horizontal branch throughout the multiwavelength campaign. We measured a significant X-ray polarization in 2--8 keV with polarization degree (PD) = $4.02 \pm 0.35$% and polarization angle (PA) = $37.6 \pm 2.5^\circ$. The energy-dependent polarization indicates that in the 2-2.5 keV energy range, the PA is much lower, $\sim9\pm8^\circ$, while other energy bands are consistent with the PA found over 2.5--8 keV. The simultaneous AstroSat-IXPE spectro-polarimetric observations provide some evidence for independent polarization from various spectral components, hinting at a disparity in the PA from the accretion disk and the Comptonized emission, while suggesting an unpolarized emission from the blackbody component. Radio observations in the 0.7--9 GHz frequency range reveal a non-detection of radio emission in 0.7-1.5 GHz and a significant detection in 5.5--9 GHz, suggesting the presence of a spectral break in 1.5-5.5 GHz. Using ATCA observation we place upper limits on the radio polarization at $<$6% on the linear polarization and $<$4% on the circular polarization at 3$\sigma$ level. We discuss the origin of the X-ray polarization and its implications on the geometry of the spectral components.

Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For 180 Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We employed three machine learning algorithms as classifiers: a Random Forest method, developed for this work, and the two classifiers used in our previous study. We studied the statistical behaviour of the s-process elements in the observational sample to investigate if the AGB models systematically under- or overpredict the abundances observed in the Ba stars and show the results in the form of violin plots of the residuals between spectroscopic abundances and model predictions. We find a significant trend in the residuals that implies an underproduction of the elements Nb, Mo, and Ru in the models relative to the observations. This may originate from a process (e.g. the intermediate neutron-capture process, i process) at the metallicity of the Ba stars not yet included in the AGB models. Correlations are found between the residuals of these elements, suggesting a common origin for the deviations. In addition, there is a weak metallicity dependence of their residuals. The s-process temperatures derived with the [Zr/Fe] - [Nb/Fe] thermometer have an unrealistic value for the majority of our stars. The most likely explanation is that at least a fraction of these elements are not produced in a steady-state s process, and instead may be due to processes not included in the AGB models. The mass distribution of the identified models confirms that our sample of Ba stars was polluted by low-mass AGB stars. Most of the matching AGB models require low accreted mass, but a few systems with high accreted mass are needed to explain the observations. (abridged)