Papers for Friday, Apr 12 2024

The Multi-channel Photometric Survey Telescope (Mephisto) is a real-time, three-color photometric system designed to capture the color evolution of stars and transients accurately. This telescope system can be crucial in cosmological distance measurements of low-redshift (low-$z$, $z$ $\lesssim 0.1$) Type Ia supernovae (SNe Ia). To optimize the capabilities of this instrument, we perform a comprehensive simulation study before its official operation is scheduled to start. By considering the impact of atmospheric extinction, weather conditions, and the lunar phase at the observing site involving the instrumental features, we simulate the light curves of SNe Ia obtained by the Mephisto. The best strategy in the case of SN Ia cosmology is to take the image at an exposure time of 130 s with a cadence of 3 days. In this condition, Mephisto can obtain hundreds of high-quality SNe Ia to achieve a distance measurement better than $4.5\%$. Given the on-time spectral classification and monitoring of the Lijiang 2.4 m Telescope at the same observatory, Mephisto, in the whole operation, can significantly enrich the well-calibrated sample of supernovae at low-$z$ and improve the calibration accuracy of high-$z$ SNe Ia.

By measuring, modeling and interpreting cosmological datasets, one can place strong constraints on models of the Universe. Central to this effort are summary statistics such as power spectra and bispectra, which condense the high-dimensional dataset into low-dimensional representations. In this work, we introduce a modern set of estimators for computing such statistics from three-dimensional clustering data, and provide a flexible Python implementation; PolyBin3D. Working in a maximum-likelihood formalism, we derive general estimators for the two- and three-point functions, which yield unbiased spectra regardless of the survey mask, weighting scheme, and presence of holes in the window function. These can be directly compared to theory without the need for mask-convolution. Furthermore, we present a numerical scheme for computing the optimal (minimum-variance) estimators for a given survey, which is shown to reduce error-bars on large-scales. Our Python package includes both general "unwindowed'' estimators and idealized equivalents (appropriate for simulations), each of which are efficiently implemented using fast Fourier transforms and Monte Carlo summation tricks, and additionally supports GPU acceleration. These are extensively validated in this work, with Monte Carlo convergence (relevant for masked data) achieved using only a small number of iterations (typically $<10$ for bispectra). This will allow for fast and unified measurement of two- and three-point functions from current and upcoming survey data.

New James Webb Space Telescope (JWST) observations are revealing the first galaxies to be prolific producers of ionizing photons, which we argue gives rise to a tension between different probes of reionization. For hydrogen reionization to proceed there must be enough ionizing photons for all of the hydrogen atoms, including their recombinations. Over the last two decades a consensus has emerged where star-forming galaxies are able to generate these photons, given reasonable values for their number densities, ionizing efficiencies $\xi_{\rm ion}$ (per unit UV luminosity), and escape fractions $f_{\rm esc}$. However, new JWST observations infer high values of $\xi_{\rm ion}$ during reionization and an enhanced abundance of earlier ($z\gtrsim 9$) galaxies, dramatically increasing the number of ionizing photons produced at high $z$. Simultaneously, recent low-$z$ studies predict significant escape fractions for faint reionization-era galaxies. Put together, we show that the galaxies we have directly observed ($M_{\rm UV} < -15$) not only can drive reionization, but would end it too early. That is, our current galaxy observations, taken at face value, imply an excess of ionizing photons and thus a process of reionization in tension with the cosmic microwave background (CMB) and Lyman-$\alpha$ forest. Considering galaxies down to $M_{\rm UV}\approx -11$, below current observational limits, only worsens this tension, requiring on average $f_{\rm esc}\approx 3\%$, far lower than expected for early galaxies from post-reionization studies. We discuss possible avenues to resolve this photon budget crisis, including missing astrophysical or observational selection effects, as well as enhanced recombinations.

Using several variants of the cosmological Simba simulations, we investigate the impact of different feedback prescriptions on the cosmic star formation history. Adopting a global-to-local approach, we link signatures seen in global observables, such as the star formation rate density (SFRD) and the galaxy stellar mass function (GSMF), to feedback effects in individual galaxies. We find a consistent picture: stellar feedback mainly suppresses star formation below halo masses of $M_{\rm H} = 10^{12} \rm \, M_{\odot}$ and before $z = 2$, whereas AGN feedback quenches the more massive systems after $z = 2$. Among Simba's AGN feedback modes, AGN jets are the dominant quenching mechanism and set the shape of the SFRD and the GSMF at late times. AGN-powered winds only suppress the star formation rate in intermediate-mass galaxies ($M_{\rm \star} = 10^{9.5 - 10} \rm \, M_{\odot}$), without affecting the overall stellar mass-assembly significantly. At late times, the AGN X-ray feedback mode mainly quenches residual star formation in massive galaxies. Our analysis reveals that this mode is also necessary to produce the first fully quenched galaxies before $z=2$, where the jets alone are inefficient. These initially highly star-forming galaxies contain relatively large black holes, likely strengthening the X-ray-powered heating and ejection of gas from the dense, central region of galaxies. Such extra heating source quenches the local star formation and produces a more variable accretion rate. More generally, this effect also causes the break down of correlations between the specific star formation rate, the accretion rate and the black hole mass.

The strong tidal force in a supermassive black hole's (SMBH) vicinity, coupled with a higher stellar density at the center of a galaxy, make it an ideal location to study the interaction between stars and black holes. Two stars moving near the SMBH could collide at a very high speed, which can result in a high energy flare. The resulting debris can then accrete onto the SMBH, which could be observed as a separate event. We simulate the light curves resulting from the fallback accretion in the aftermath of a stellar collision near a SMBH. We investigate how it varies with physical parameters of the system. With all other physical parameters of the system held constant, the direction of the relative velocity vector at time of impact plays a large role in determining the overall form of the light curve. One distinctive light curve we notice is characterized by a sustained increase in the luminosity some time after accretion has started. We compare this form to the light curves of some candidate tidal disruption events (TDEs). Stellar collision accretion flares can take on unique appearances that would allow them to be easily distinguished, as well as elucidate underlying physical parameters of the system. There exist several ways to distinguish these events from TDEs, including the much wider range of SMBH masses stellar collisions may exist around.

We propose a simple model that can alleviate the $H_0$ tension while remaining consistent with big bang nucleosynthesis (BBN). It is based on a dark sector described by a standard Lagrangian featuring a $SU(N)$ gauge symmetry with $N\geq3$ and a massive scalar field with a quartic coupling. The scalar acts as dark Higgs leading to spontaneous symmetry breaking $SU(N)\to SU(N\!-\!1)$ via a first-order phase transition \`a la Coleman-Weinberg. This set-up naturally realizes previously proposed scenarios featuring strongly interacting dark radiation (SIDR) with a mass threshold within hot new early dark energy (NEDE). For a wide range of reasonable model parameters, the phase transition occurs between the BBN and recombination epochs and releases a sufficient amount of latent heat such that the model easily respects bounds on extra radiation during BBN while featuring a sufficient SIDR density around recombination for increasing the value of $H_0$ inferred from the cosmic microwave background. Our model can be summarized as a natural mechanism providing two successive increases in the effective number of relativistic degrees of freedom after BBN but before recombination $\Delta N_\mathrm{BBN} \to \Delta N_\mathrm{NEDE} \to \Delta N_\mathrm{IR}$ alleviating the Hubble tension. The first step is related to the phase transition and the second to the dark Higgs becoming non-relativistic. This set-up predicts further signatures, including a stochastic gravitational wave background and features in the matter power spectrum that can be searched for with future pulsar timing and Lyman-$\alpha$ forest measurements.

Most stars form in clusters and groups rather than in isolation. We present $\lesssim 5^{\prime\prime}$ angular resolution ($\sim 2000$ au, or 0.01 pc) Very Large Array NH$_3$ (1,1), (2,2), and (3,3) and 1.3 cm continuum emission observations of the dense gas within the Serpens South protocluster and extended filaments to the north and south. We identify 94 dense cores using a dendrogram analysis of the NH$_3$ (1,1) integrated intensity. Gas temperatures $T_K$ and non-thermal linewidths $\sigma_\mathrm{NT}$ both increase towards the centre of the young stellar cluster, in the dense gas generally and in the cores specifically. We find that most cores (54\%) are super-virial, with gravitationally bound cores located primarily in the filaments. Cores in the protocluster have higher virial parameters by a factor $\sim 1.7$, driven primarily by the increased core $\sigma_\mathrm{NT}$ values. These cores cannot collapse to form stars unless they accrete additional mass or their core internal motions are reduced. The southern filament shows a significant velocity gradient previously interpreted as mass flow toward the cluster. We find more complex kinematics in the northern filament. We find a strong correlation between $\sigma_\mathrm{NT}$ and $T_K$, and argue that the enhanced temperatures and non-thermal motions are due to mechanical heating and interaction between the protocluster-driven outflows and the dense gas. Filament-led accretion may also contribute to the increased $\sigma_\mathrm{NT}$ values. Assuming a constant fraction of core mass ends up in the young stars, future star formation in the Serpens South protocluster will shift to higher masses by a factor $\sim 2$.

Exoplanets in their infancy are ideal targets to probe the formation and evolution history of planetary systems, including the planet migration and atmospheric evolution and dissipation. In this paper, we present spectroscopic observations and analyses of two planetary transits of K2-33b, which is known to be one of the youngest transiting planets (age $\approx 8-11$ Myr) around a pre-main-sequence M-type star. Analysing K2-33's near-infrared spectra obtained by the IRD instrument on Subaru, we investigate the spin-orbit angle and transit-induced excess absorption for K2-33b. We attempt both classical modelling of the Rossiter-McLaughlin (RM) effect and Doppler-shadow analyses for the measurements of the projected stellar obliquity, finding a low angle of $\lambda=-6_{-58}^{+61}$ deg (for RM analysis) and $\lambda=-10_{-24}^{+22}$ deg (for Doppler-shadow analysis). In the modelling of the RM effect, we allow the planet-to-star radius ratio to float freely to take into account the possible smaller radius in the near infrared, but the constraint we obtain ($R_p/R_s=0.037_{-0.017}^{+0.013}$) is inconclusive due to the low radial-velocity precision. Comparison spectra of K2-33 of the 1083 nm triplet of metastable ortho-He I obtained in and out of the 2021 transit reveal excess absorption that could be due to an escaping He-rich atmosphere. Under certain conditions on planet mass and stellar XUV emission, the implied escape rate is sufficient to remove an Earth-mass H/He in $\sim$1 Gyr, transforming this object from a Neptune to a super-Earth.

Multiple analytical and empirical stability criteria have been derived in the literature for two planet systems. But, the dependence of the stability limit on the initial mutual inclination between the inner and outer orbits is not well modeled by previous stability criteria. Here, we derive a semi-analytical stability criteria for two planet systems, at arbitrary inclinations, in which the inner planet is a test particle. Using perturbation theory we calculate the characteristic fractional change in the semi-major axis of the inner binary $\beta=\delta a_1/a_1$ caused by perturbations from the companion. A stability criteria can be derived by setting a threshold on $\beta$ Focusing initially on circular orbits, we derive an analytical expression for $\beta$ for co-planar prograde and retrograde orbits. For non-coplanar configurations, we evaluate a semi-analytical expression. We then generalize to orbits with arbitrary eccentricities and account for the secular effects. Our analytical and semi-analytical results are in excellent agreement with direct N-body simulations. In addition, we show that contours of $\beta\sim0.01$ can serve as criteria for stability. More specifically, we show that (1) retrograde orbits are generally more stable than prograde ones; (2) systems with intermediate mutual inclination are less stable due to vZLK dynamics; and (3) mean-motion resonances (MMRs) can stabilize intermediate inclination secularly unstable regions in phase space, by quenching vZLK secular processes (4) MMRs can destabilize some of the dynamically stable regions. We also point out that these stability criteria can be used to constrain the orbital properties of observed systems and their age. We also point out that these stability criteria can be used to constrain the orbital properties of observed systems (in particular inclination) and their age.

Given the important role turbulence plays in the settling and growth of dust grains in protoplanetary disks, it is crucial that we determine whether these disks are turbulent and to what extent. Protoplanetary disks are weakly ionized near the mid-plane, which has led to a paradigm in which largely laminar magnetic field structures prevail deeper in the disk, with angular momentum being transported via magnetically launched winds. Yet, there has been little exploration on the precise behavior of the gas within the bulk of the disk. We carry out 3D, local shearing box simulations that include all three low-ionization effects (Ohmic diffusion, ambipolar diffusion, and the Hall effect) to probe the nature of magnetically driven gas dynamics 1-30 AU from the central star. We find that gas turbulence can persist with a generous yet physically motivated ionization prescription (order unity Elsasser numbers). The gas velocity fluctuations range from 0.03-0.09 of the sound speed $c_s$ at the disk mid-plane to $\sim c_s$ near the disk surface, and are dependent on the initial magnetic field strength. However, the turbulent velocities do not appear to be strongly dependent on the field polarity, and thus appear to be insensitive to the Hall effect. The mid-plane turbulence has the potential to drive dust grains to collision velocities exceeding their fragmentation limit, and likely reduces the efficacy of particle clumping in the mid-plane, though it remains to be seen if this level of turbulence persists in disks with lower ionization levels.

In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full-Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying $n_s$ and allowing for $w_0w_a$CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within $1\sigma$ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models.

The Dark Energy Spectroscopic Instrument (DESI) will provide unprecedented information about the large-scale structure of our Universe. In this work, we study the robustness of the theoretical modelling of the power spectrum of FOLPS, a novel effective field theory-based package for evaluating the redshift space power spectrum in the presence of massive neutrinos. We perform this validation by fitting the AbacusSummit high-accuracy $N$-body simulations for Luminous Red Galaxies, Emission Line Galaxies and Quasar tracers, calibrated to describe DESI observations. We quantify the potential systematic error budget of FOLPS, finding that the modelling errors are fully sub-dominant for the DESI statistical precision within the studied range of scales. Additionally, we study two complementary approaches to fit and analyse the power spectrum data, one based on direct Full-Modelling fits and the other on the ShapeFit compression variables, both resulting in very good agreement in precision and accuracy. In each of these approaches, we study a set of potential systematic errors induced by several assumptions, such as the choice of template cosmology, the effect of prior choice in the nuisance parameters of the model, or the range of scales used in the analysis. Furthermore, we show how opening up the parameter space beyond the vanilla $\Lambda$CDM model affects the DESI observables. These studies include the addition of massive neutrinos, spatial curvature, and dark energy equation of state. We also examine how relaxing the usual Cosmic Microwave Background and Big Bang Nucleosynthesis priors on the primordial spectral index and the baryonic matter abundance, respectively, impacts the inference on the rest of the parameters of interest. This paper pathways towards performing a robust and reliable analysis of the shape of the power spectrum of DESI galaxy and quasar clustering using FOLPS.

Physical properties of ten millimeter-sized meteoroids from the Southern Delta Aquariids (SDA) shower are derived using optical observations from the Canadian Automated Meteor Observatory between 2020 and 2023. The meteors are found to ablate in two distinct erosion stages, the second stage showing a single, bright leading fragment. Our modelling interprets these observations as evidence for equal masses of compact grains embedded in a porous, low density matrix. The average bulk density of SDA meteors is found to be 1,420 $\pm$ 100 $\rm{kg \ m^{-3}}$, with the compact component having a density of 2,310 $\pm$ 160 $\rm{kg \ m^{-3}}$ and the porous component a density of $700\pm110$ $\rm{kg \ m^{-3}}$. The high bulk density of SDA meteors is comparable to densities found for the Quadrantid and Geminid showers, both of which also have low perihelion distances. This suggests that thermal desorption may play a significant role in the processing of meteoroids.

A major ingredient for kilonova lightcurves is the radioactive heating rate and its dependence on the electron fraction and velocity of the ejecta and, in principle, on the nuclear mass formula. Heating-rate formulae commonly used as the basis for kilonova models are, strictly speaking, incorrect for electron fractions other than $Y_{e} = 0.04$. Here, we introduce new semi-analytical models for kilonovae with better heating rate prescriptions valid for the full parameter space of kilonova velocities and electron fractions. This new prescription produces, on average, dimmer kilonovae at peak that decay more slowly as compared to previously used prescriptions for otherwise identical kilonova physics. We show the dangers of using inappropriate heating rate estimates by simulating realistic observations and inferring the kilonova parameters via a misspecified heating-rate prescription. While providing great fits to the photometry, an incorrect heating-rate prescription fails to recover the input ejecta masses at greater than $5\sigma$. This bias from an incorrect prescription has disastrous consequences for interpreting kilonovae, their use as additional components in gamma-ray burst afterglows, and understanding their role in cosmic chemical evolution or for multi-messenger constraints on the nuclear equation of state. Given the true heating-rate is uncertain, we estimate there is a $\approx 5$ and $\approx 10\%$ systematic uncertainty in the measured ejecta masses and velocities, respectively. For lanthanide-rich ejecta, the systematic uncertainty could be as high as $\approx 50\%$. This systematic uncertainty limits the precision of any measurements that rely on accurate estimates of kilonova ejecta properties.

DESI aims to provide one of the tightest constraints on cosmological parameters by analyzing the clustering of more than thirty million galaxies. However, obtaining such constraints requires special care in validating the analysis methods, and efforts to reduce the computational time required through techniques such as data compression and emulation. In this work, we perform a precision validation of the PyBird power spectrum modelling code with both a traditional, but emulated, Full-Modelling approach and the model-independent Shapefit compression approach. Using cubic simulations, which accurately reproduce the clustering and precision of the DESI survey, we find that the cosmological constraints from Shapefit and Full-Modelling are consistent with each other at the $\sim0.3\sigma$ level. Both Shapefit and Full-Modelling are also consistent with the true $\Lambda$CDM simulation cosmology, even when including the hexadecapole, down to a scale $k_{\mathrm{max}} = 0.20 h \mathrm{Mpc}^{-1}$. For extended models such as the $w$CDM and the $o$CDM models, we find including the hexadecapole can significantly improve the constraints and reduce the systematic errors with the same $k_{\mathrm{max}}$. Furthermore, we also show that the constraints on cosmological parameters with the correlation function evaluated from PyBird down to $s_{\mathrm{min}} = 30 h^{-1} \mathrm{Mpc}$ are unbiased, and consistent with the constraints from the power spectrum.

Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$\alpha$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$\alpha$ flux to estimate the Ly$\alpha$ production rate, we find that $\sim 4\%$ of the Ly$\alpha$ photons are absorbed by water. Combined with direct measurements of water emissions in the $\nu_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$\alpha$ photons in the vicinity of a shock front.

Stars are amongst the most fundamental structures of our Universe. They comprise most of the baryonic and luminous mass of galaxies, synthethise heavy elements, and injec\ t mass, momentum, and energy into the interstellar medium. They are also home to the planets. Since stellar properties are primarily decided by their mass, the so-called \ stellar initial mass function (IMF) is critical to the structuring of our Universe. We review the various physical processes, and theories which have been put forward as well as the numerical simulations which have been carried out to explain the origin of the stellar initial mass function. Key messages from this review are: (1) Gravity and turbulence most likely determine the power-law, high-mass part of the IMF. (2) Depending of the Mach number and the density distribution, several regimes are possible, including $\Gamma _{IMF} \simeq 0$, -0.8, -1 or -1.3 where $d N / d \log M \propto M^{\Gamma_{IMF}}$. These regimes are likely universal, however the transition between these regimes is not. (3) Protostellar jets can play a regulating influence on the IMF by injecting momentum into collapsing clumps and unbinding gas. (4) The peak of the IMF may be a consequence of dust opacity and molecular hydrogen physics at the origin of the first hydrostatic core. This depends weakly on large scale environmental conditions such as radiation, magnetic field, turbulence or metallicity. This likely constitutes one of the reason of the relative universality of the IMF.

We explore the impact of dynamical friction on cosmological scales and show its influence on the evolution of perturbations. In particular, considering smooth and clustering dark energy models, we describe the role played by friction by selecting two main hierarchical models, \textit{i.e.}, the first where the friction term is proportional to the Hubble rate, whereas the second where friction is induced by the dark energy pressure. The second approach generalises the first and translates the idea that pressure is a general relativistic effect, motivating why friction might arise once barotropic dark energy fluids are considered. The corresponding effects of friction are investigated at the level of linear and nonlinear perturbations, using the formalism of the spherical collapse model. Whilst dynamical friction has very small effects and thus it cannot be excluded \emph{a priori}, dissipative pressure friction leads to a substantial slow down in the evolution of perturbations. This can be particularly inferred within the halo mass function, for which we also employ corrections due to dark energy clustering. To this end, in order to discern detectable deviations from the standard cosmological model, we thus highlight where dissipation effects might play a significant role at large scales.

Data-driven models for stellar spectra which depend on stellar labels suffer from label systematics which decrease model performance: the "stellar labels gap". To close the stellar labels gap, we present a stellar label independent model for $\textit{Gaia}$ BP/RP (XP) spectra. We develop a novel implementation of a variational auto-encoder; a $\textit{scatter}$ VAE, which learns to generate an XP spectrum and intrinsic scatter without relying on stellar labels. We demonstrate that our model achieves competitive XP spectra reconstructions in comparison to stellar label dependent models. We find that our model learns stellar properties directly from the data itself. We then apply our model to XP/APOGEE giant stars to study the [$\alpha$/M] information in $\textit{Gaia}$ XP. We provide strong evidence that the XP spectra contain meaningful [$\alpha$/M] information by demonstrating that our model learns the $\alpha$-bimodality $\textit{without relying on stellar label correlations}$, for stars with $T_{\rm eff} <$ 5000 K. We publicly release our trained model, codebase and data. Importantly, our stellar label independent model can be implemented for any/all XP spectra because our model performance scales with training object density, not training label density.

We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70% of the simulated galaxies exhibit an ``inside-out'' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside-out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally- and radially-integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions.

Recent studies of the neutrino-driven wind from proto-neutron stars have indicated that the wind is likely proton-rich for much of its lifetime, and the high flux of neutrinos can induce $\nu$p-process nucleosynthesis allowing for the formation of heavy elements. It has also been shown that gravito-acoustic waves, generated by convection within the proto-neutron star, can significantly alter the dynamics and nucleosynthesis in the wind. Therefore, we present a study of the effects of convection-driven waves on the nucleosynthesis in proton-rich neutrino-driven winds, focusing on the $\nu$p-process. We find that wave effects can strongly impact $\nu$p-process nucleosynthesis even at wave luminosities a factor of $10^{-5}$ smaller than the total neutrino luminosity. The momentum flux of the waves accelerates the wind, reducing the net neutrino heating and the persistent neutron abundance created by p($\bar{\nu}_e,e^+$), which impedes $\nu$p-process nucleosynthesis. However, this effect is generally counteracted by the effects of waves on seed nucleus formation, as the acceleration of the wind and the heating that occurs as these waves shock both favor a more $\alpha$-rich environment with very little heavy seed-nucleus formation. Overall, higher wave luminosities correlate (albeit non-monotonically) with heavier element $\nu$p-processing, up to $A\approx 200$ in some cases. At very high wave luminosities ($\gtrsim 10^{-3}L_\nu$), early shock heating by the waves disrupts $\alpha$ recombination, and drives a suppressed, fast-outflow r-process proceeding up to $A\approx 200$. This occurs despite an assumed neutrino spectrum that predicts a proton-rich wind with equilibrium $Y_e=0.6$.

Understanding the role of magnetic reconnection in the heating and dynamics of the solar atmosphere requires detailed observational data of any observable aspect of the reconnection process, including small-scale features such as plasmoids. We aim to examine the capability of active and upcoming instruments to detect plasmoids generated by reconnection in the corona including low-density regimes. We used the Bifrost code to perform simulations of plasmoid-mediated reconnection in the corona based on a 2D idealised setup: a fan-spine topology with uniform density including thermal conduction. Through forward-modelling of extreme-ultraviolet (EUV) observables, we checked whether our simulated plasmoids could be detected with the instruments of Solar Dynamics Observatory (SDO) and Solar Orbiter (SO), as well as the upcoming Multi-Slit Solar Explorer (MUSE) and Solar-C missions. Short-lived (10-20 s) small-scale (0.2-0.5 Mm) coronal plasmoids are not resolvable with the Atmospheric Imaging Assembly (AIA) onboard SDO, but could be captured with the Extreme Ultraviolet Imager (EUI) of SO. The spatial and temporal high-resolution planned for the MUSE spectrograph (SG) is adequate to obtain full spectral information of these plasmoids. Detection of 0.8 MK plasmoids in the MUSE/SG 171 {\AA} channel should work on full-raster mode in regions with electron densities above 10^9 cm^3 whereas on sit-and-stare mode for lower-density regions. Solar-C could also capture these coronal plasmoids using the EUV High-Throughput Spectroscopic Telescope (EUVST), through rapid changes in Doppler shift and line width in different EUV lines caused by plasmoid motions along the current sheet. With combined spectra of MUSE/SG and Solar-C/EUVST in multiple emission lines, along with high-resolution images from SO/EUI and MUSE/CI, it should be possible to gain new insights about plasmoid formation in the corona.

Some of the Intermediate Mass Black Hole (IMBH) candidates observed at the center of galaxies or in globular clusters and some of the Supermassive Black Holes (SMBHs) seen at the center of many galaxies might be of primordial origin. Indeed, Primordial Black Holes (PBHs) of such mass could have formed when the Universe was $\sim$1-10$^3$ s old, due to the collapse of density fluctuations. In particular, when the Universe was $\sim 1$ s in age, Electron-Positron Annihilation (EPA) took place. We explore the formation of intermediate mass and supermassive PBHs, taking into account the effect of the EPA when the fluctuations have a running-tilt power-law spectrum: when these cross the $10^{-0.5}$-$10^{3.0}$ s Universe horizon they could produce $5\times 10^{3}$ - $5\times 10^{8}M_{\odot}$ PBHs with a density as high as $\sim 10^{10}$/Gpc$^3$. On average, this implies a population of about one thousand PBHs in the Local Group of Galaxies, with the nearest one at about 250 kpc, just under half the distance to the Andromeda galaxy (M31).

Methanol masers at 6.7~GHz are well-known signposts of high-mass star-forming regions. [...] We aim to understand the origin of the ring-like structures outlined by methanol maser emission in a number of sources. This emission could be, a priori, spatially associated with an outflow and/or disc around a high-mass protostar. [...] Using sensitive, three-epoch observations spanning over eight years with the European VLBI Network, we have started the most direct investigations of maser rings using very accurate proper motion measurements with uncertainties below 1\,km~s$^{-1}$. We present full results for the five targets of our sample, G23.207-00.377, G23.389+00.185, G28.817+00.365, G31.047+00.356, and G31.581+00.077, where proper motions show similar characteristics; maser cloudlets do not move inwards towards the centre of the rings but rather outwards. We also include the most circular source, G23.657-00.127, in the discussion as a reference. The magnitude of maser proper motions ranges from a maximum of about 13\,km~s$^{-1}$ to 0.5~km~s$^{-1}$. In two of the five sources with a high number of maser spots (>100), namely G23.207-00.377 and G23.389+00.185, we show that the size of the best elliptical model, fitted to the distribution of persistent masers, increases in time in a manner similar to the case of G23.657-00.127. Moreover, we checked the separations between the pairs of spots from distinct regions, and we were able to assess that G28.817+00.365 and G31.047+00.356 can be interpreted as showing expanding motions. We analysed the profiles of single maser cloudlets and studied their variability. Contrary to single-dish studies, the interferometric data indicate variability of the emission of single-masing cloudlets. In five of the six targets expansion motions prevail. Only in the case of G31.581+00.077 can a scenario of disc-like rotation not be excluded. [...]

Ultra-hot Jupiters are likely doomed by tidal forces to undergo orbital decay and eventual disruption by their stars, but the timescale over which this process unfolds is unknown. We present results from a long-term project to monitor ultra-hot Jupiters transits. We recovered WASP-12 b's orbital decay rate of dP/dt = -29.8 +/- 1.6 ms yr-1, in agreement with prior work. Five other systems initially had promising non-linear transit ephemerides. However, a closer examination of two -- WASP-19 b and CoRoT-2 b, both with prior tentative detections -- revealed several independent errors with the literature timing data; after correction neither planet shows signs of orbital decay. Meanwhile, a potential decreasing period for TrES-1 b, dP/dt = -16 +/- 5 ms yr-1, corresponds to a tidal quality factor Q*' = 160 and likely does not result from orbital decay, if driven by dissipation within the host star. Nominal period increases in two systems, WASP-121 b and WASP-46 b, rest on a small handful of points. Only 1/43 planets (WASP-12 b) in our sample is experiencing detectable orbital decay. For nearly half (20/42) we can rule out dP/dt as high as observed for WASP-12 b. Thus while many ultra-hot Jupiters could still be experiencing rapid decay that we cannot yet detect, a sizeable sub-population of UHJs are decaying at least an order of magnitude more slowly than WASP-12 b. Our reanalysis of Kepler-1658 b with no new data finds that it remains a promising orbital decay candidate. Finally, we recommend that the scientific community take steps to avoid spurious detections through better management of the multi-decade-spanning datasets needed to search for and study planetary orbital decay.

We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$\beta$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of the 14 objects. These results include H$\beta$ time lags, supermassive black hole masses, and velocity-resolved time lags. The velocity-resolved lags allow us to investigate the kinematics of the broad-line region.

The composition and distribution of the gas in a protoplanetary disk plays a key role in shaping the outcome of the planet formation process. Observationally, the recovery of information such as the emission height and brightness temperature from interferometric data is often limited by the imaging processes. To overcome the limitations of image-reconstruction when analyzing gas emission from interferometric observations, we have introduced a parametric model to fit the main observable properties of the gaseous disk component in the visibility plane. This approach is also known as parametric visibility modeling. We applied our parametric visibility modeling to the gas brightness distribution of the molecular line emission from 12CO J=3-2 and 13CO J=3-2 in the disk around MHO 6, a very-low-mass star in the Taurus star-forming Region. To improve the flux fidelity of our parametric models, we combined models with different pixel resolution before the computation of their visibilities, referred to as ``nesting images.'' When we apply our parametric visibility modeling to MHO 6, with independent fits to the emission from its CO isopotologues, the models return the same consistent results for the stellar mass, disk geometry, and central velocity. The surface height and brightness temperature distribution are also recovered. When compared to other disks, MHO 6 surface height is among the most elevated surfaces, consistent with the predictions for disks around very-low-mass stars. This work demonstrates the feasibility of running rapidly iterable parametric visibility models in moderate resolution and sensitivity interferometric observations. More importantly, this methodology opens the analysis of disk's gas morphology to observations where image-based techniques are unable to robustly operate, as in the case of the compact disk around MHO 6.

The changes in colors across a galaxy are intimately connected to the galaxy's formation, growth, quenching history, and dust content. A particularly important epoch in the growth of galaxies is near $z \sim 2$ often referred to as 'cosmic noon', where galaxies on average reach the peak of their star formation. We study a population of 125 cluster galaxies at $z \sim 1.6$ in three Hubble Space Telescope (HST) filters, F475W, F625W, and F160W, roughly corresponding to the rest-frame FUV, NUV, and r band, respectively. By comparing to a control sample of 200 field galaxies at similar redshift, we reveal clear, statistically significant differences in the overall spatially resolved colors and color gradients in galaxies across these two different environments. On average, cluster galaxies have redder UV colors in both the inner and outer regions bounded by $r_{\mathrm{50}}$, as well as an overall wider dispersion of outside-in color gradients. The presence of these observed differences, along with evidence from ancillary data from previous studies, strongly suggests that the environment drives these population-level color differences, by affecting the stellar populations and/or dust content.

The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N$_2$H$^+$ (1$-$0), with a critical density of $2\times10^5$~cm$^{-3}$, and spatial resolution $\sim$0.02~pc. G353.41, at a distance of 2~kpc, is embedded in a larger scale ($\sim$8~pc) filament and has a mass of 2500~M$_{\odot}$ within $1.3\times1.3$~pc$^2$. We extract the N$_2$H$^+$ isolated line component and we decompose it by fitting up to 3 Gaussian velocity components. This allows us to identify velocity structures that are either muddled or impossible to identify in the traditional position-velocity diagram. We identify multiple velocity gradients (VGs) on large and small scales. We find good agreement between the N$_2$H$^+$ and the previously reported DCN core velocities, suggesting that cores are kinematically coupled to the dense gas in which they form. We measure 9 converging V-shaped VGs, located in filaments, that are sometimes associated with cores near their point of convergence. The average timescale associated with the V-shapes are $\sim$67~kyr, or about twice the free-fall time of cores in the same area ($\sim$~33~kyr) but substantially shorter than protostar lifetime estimates ($\sim$~0.5~Myr). We interpret these V-shapes as inflowing gas feeding the regions near cores and we derive their mass accretion rates. This feeding might lead to further filament collapse and formation of new cores. We suggest that the protocluster is collapsing on large scales, but the velocity signature of collapse is slow compared to pure free-fall. Thus these data are consistent with a comparatively slow global protocluster contraction under gravity, and faster core formation within, suggesting the formation of multiple generations of stars over the protocluster lifetime.

Context. Time-series photometry has given astronomers the tools to study time-dependent astrophysical phenomena, from stellar activity to fast radio bursts and exoplanet transits. Transit events in particular are focused primarily on planetary transits, and eclipsing binaries with eclipse geometries that are parameterised with a few variables, while more complex light curves caused by substructure within the transiting object require customized analysis code. Aims. We present Beyond Circular Eclipsers (BeyonCE), which reduces the parameter space encompassed by the transit of circum-secondary disc (CSD) systems with azimuthally symmetric non-uniform optical depth profiles. By rejecting disc geometries that cannot reproduce the measured gradients within their light curves, we can constrain the size and orientation of discs with complex sub-structure. Methods. We map out all the possible geometries of a disc, calculate the gradients for rings crossing the star, then reject those configurations where the measured gradient of the light curve is greater than the theoretical gradient from the given disc orientation. Results. We present the fitting code BeyonCE and demonstrate its effectiveness in considerably reducing the parameter space of discs that contain azimuthally symmetric structure by analyzing the light curves seen towards J1407 and PDS 110 which are attributed to CSD transits.

We present a catalogue of white dwarf candidates constructed from the GALEX and Gaia EDR3 catalogues. The catalogue contains 332,111 candidate binary white dwarf systems and 111,996 candidate single white dwarfs. Where available, the catalogue is augmented with photometry from Pan-STARRS DR1, SDSS DR12 and classifications from StarHorse. We fit photometric data with modeled white dwarf cooling sequences to derive mass, age and effective temperature of the white dwarf as well as mass estimates for the companion. We test our classifications against StarHorse, the Gentile-Fusillo Gaia EDR3 catalogue, and white-dwarf-main-sequence binaries identified in SDSS DR12. This catalogue provides a unique probe of the binarity of white dwarfs as well as the abundance of white-dwarf giant binaries and large mass-ratio stellar binaries which are difficult to probe otherwise.

Accurate age estimation of nearby young moving groups (NYMGs) is important as they serve as crucial testbeds in various fields of astrophysics, including formation and evolution of stars, planets, as well as loose stellar associations. The $\beta$-Pictoris moving group (BPMG), being one of the closest and youngest NYMGs, has been extensively investigated, and its estimated ages have a wide range from $\sim$10 to 25 Myr, depending on the age estimation methods and data used. Unlike other age dating methods, kinematic traceback analysis offers a model-independent age assessment hence the merit in comparing many seemingly discordant age estimates. In this study, we determine the kinematic ages of the BPMG using three methods: probabilistic volume calculation, mean pairwise distance calculation, and covariance matrix analysis. These methods yield consistent results, with estimated ages in the range of 14 to 20 Myr. Implementing corrections to radial velocities due to gravitational redshift and convectional blueshift increases the ages by $\sim2-4$ Myr. Conversely, considering data uncertainties decreases the estimated ages by 1 to 2 Myr. Taken together, our analysis determined the kinematic age of BPMG to be 16.3$^{+3.4}_{-2.1}$ Myr. This age is significantly younger than the commonly accepted age of the BPMG ($\sim$24 Myr) determined primarily from the lithium depletion boundary analysis. This younger kinematic age may point to the discrepancy between the luminosity evolution and lithium depletion models or the presence of unaccounted systematic error in the method. This result underscores the necessity for systematic reevaluations of age-dating methods for nearby, young moving groups.

The properties of relativistic jets, their interaction with the environment, and their emission of radiation can be self-consistently studied by using particle-in-cell (PIC) numerical simulations. Using three-dimensional (3D), relativistic PIC simulations, we present the first self-consistently calculated synthetic spectra of head-on and off-axis emission from electrons accelerated in cylindrical, relativistic plasma jets containing an initial toroidal magnetic field. The jet particles are initially accelerated during the linear stage of growing plasma instabilities, which are the Weibel instability (WI), kinetic Kelvin-Helmholtz instability (kKHI), and mushroom instability (MI). In the nonlinear stage, these instabilities are dissipated and generate turbulent magnetic fields which accelerate particles further. We calculate the synthetic spectra by tracing a large number of jet electrons in the nonlinear stage, near the jet head where the magnetic fields are turbulent. Our results show the basic properties of jitter-like radiation emitted by relativistic electrons when they travel through a magnetized plasma with the plasma waves driven by kinetic instabilities (WI, kKHI, and MI) growing into the nonlinear regime. At low frequencies, the slope of the spectrum is ~ 0.94, which is similar to that of the jitter radiation. The results are relevant to active galactic nuclei/blazars and gamma-ray burst jet emission and set the ground for future studies on synthetic spectra from relativistic jets.

Among the silicon bearing species discovered in the interstellar medium, SiS and SiO stand out as key tracers due to their distinct chemistry and abundances in interstellar and circumstellar environments. Our objective is to enhance the network of Si- and S-bearing chemical reactions for a gas-grain model in molecular clouds, encompassing both low and high metallicities. We have calculated the energies and rate coefficients for 6 neutral atom-diatom reactions involved in the SiCS triatomic system, with a special focus on the C+SiS and S+SiC collisions. We employ the coupled cluster method with single and double substitutions and a perturbative treatment of triple substitutions (CCSD(T)) refined at the explicitly correlated CCSD(T)-F12 level. With these computational results in conjunction with data from the literature, we construct an extended network of neutral-neutral chemical reactions. We performed time-dependent models employing the Nautilus gas-grain code, setting the gas temperature to 10 K and the density to 2x10$^4$ cm$^{-3}$. The temperature dependence for the reactions involving SiS were modelled using $k(T)=\alpha \left( T/300 \right)^{\beta} \exp{(-\gamma/T)}$. The high-metallicity models significantly boost the SiS production. Higher initial abundances of C, S, and Si, roughly $\sim$2, 190, and 210 times higher, respectively, contribute to this. Around 10$^3$ yr, destruction mechanisms become relevant. The proposed production reaction S + SiC $\rightarrow$ C + SiS, mitigates these effects. By expanding the gas reaction network using a high metallicity model, we derived estimates for the abundances of interstellar molecules. The inclusion of neutral-neutral mechanisms, particularly via Si+HS and S+SiC channels, played a pivotal role in determining SiS abundance. These mechanisms carry a significance on a par with the well-known and fast ion-neutral reactions.

Stellar black holes (sBHs) are widely believed to exist in the accretion disks of active galactic nuclei (AGNs). Previous studies often focus on the transient emission produced by embedded sBHs. Here, we explore the possible observational consequences of an AGN accretion disk that contains a population of accreting sBHs. Embedded accreting sBHs change the effective temperature distribution of the AGN accretion disk by heating gas in the outer regions. Two possible observational consequences are presented. First, the spectral energy distribution has a turnover feature at $\sim 4700\ \textrm{\AA}$ when the supermassive black hole (SMBH) mass is $\sim 10^8\ M_{\odot}$, which can help explain the observed shallow spectral shape at wavelengths $>5000\ \textrm{\AA}$ for the Sloan Digital Sky Survey quasar composite spectrum. Second, the half-light radius of a given relatively long wavelength is significantly larger than for an AGN disk without sBHs, which can be tested by microlensing observations. With appropriate sBH distributions, the model can be reconciled with quasar microlensing disk sizes. We propose that the half-light radius-wavelength relation can be utilized to investigate the distributions of embedded sBHs in AGN accretion disks.

We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and proper motion searches as well as presenting hundreds of newly detected stars from our search of Australian SKA Pathfinder observations. The focus of this first version of the catalogue is on objects detected in surveys using SKA precursor instruments; however we will expand this scope in future versions. The 839 objects in the Sydney Radio Star Catalogue are distributed across the whole sky and range from ultracool dwarfs to Wolf-Rayet stars. We find that the radio luminosities of cool dwarfs are lower than the radio luminosities of more evolved sub-giant and giant stars. We use X-ray detections of 530 radio stars by the eROSITA soft X-ray instrument onboard the SRG spacecraft to show that almost all of the radio stars in the catalogue are over-luminous in the radio, indicating that the majority of stars at these radio frequencies are coherent radio emitters. The Sydney Radio Star Catalogue can be found in Vizier or at https://radiostars.org.

In this Letter, we propose a model-independent method to determine the Hubble constant and curvature simultaneously taking advantage of the possibilities of future space-borne gravitational wave (GW) detector DECIGO in combination with the radio quasars as standard rulers. Similarly to the redshift drift in the electromagnetic domain, accelerating expansion of the Universe causes a characteristic phase correction to the gravitational waveform detectable by DECIGO. Hence, one would be able to extract the Hubble parameter $H(z)$. This could be used to recover distance-redshift relation supported by the data not relying on any specific cosmological model. Assuming the FLRW metric, and using intermediate luminosity radio quasars as standard rulers one achieves an interesting opportunity to directly assess $H_0$ and $\Omega_k$ parameters. To test this method we simulated a set of acceleration parameters achievable by future DECIGO. Based on the existing sample of 120 intermediate-luminosity radio-quasars calibrated as standard rulers, we simulated much bigger samples of such standard rulers possible to obtain with VLBI. In the case of $(N=100)$ of radio quasars, which is the size of currently available sample, the precision of cosmological parameters determined would be $\sigma_{H_0}=2.74$ ${\mathrm{~km~s^{-1}~Mpc^{-1}}}$ and $\sigma_{\Omega_k}=0.175$. In the optimistic scenario $(N = 1000)$ achievable by VLBI, the precision of $H_{0}$ would be improved to $1\%$, which is comparable to the result of $\sigma_{H_0} =0.54$ ${\mathrm{~km~s^{-1}~Mpc^{-1}}}$ from \emph{Planck} 2018 TT, TE, EE+lowE+lensing data, and the precision of $\Omega_k$ would be 0.050. Our results demonstrate that such combined analysis, possible in the future, could be helpful to solve the current cosmological issues concerning the Hubble tension and cosmic curvature tension.

I further study the manner by which a pair of opposite jets shape the `keyhole' morphological structure of the core-collapse supernova (CCSN) SN 1997A, now the CCSN remnant (CCSNR) 1987A. By doing so, I strengthen the claim that the jittering-jet explosion mechanism (JJEM) accounts for most, likely all, CCSNe. The `keyhole' structure comprises a northern low-intensity zone closed with a bright rim on its front and an elongated low-intensity nozzle in the south. This rim-nozzle asymmetry is observed in some cooling flow clusters and planetary nebulae that are observed to be shaped by jets. I build a toy model that uses the planar jittering jets pattern, where consecutive pairs of jets tend to jitter in a common plane, implying that the accreted gas onto the newly born neutron star at the late explosion phase flows perpendicular to that plane. This allows for a long-lived jet-launching episode. This long-lasting jet-launching episode launches more mass into the jets that can inflate larger pairs of ears or bubbles, forming the main jets' axis of the CCSNR that is not necessarily related to a possible pre-collapse core rotation. I discuss the relation of the main jets' axis to the neutron star's natal kick velocity.

We follow up emission line galaxies identified through the near-infrared slitless HST/WFC3 WISP survey with VLT/FORS2 optical spectroscopy. Over 4 WISP fields, we targetted 85 of 138 line emission objects at $0.4<z<2$ identified in WFC3 spectroscopy. Half the galaxies are fainter than $H_{AB}=24$mag, and would not have been included in many well-known surveys based on broad-band magnitude selection. We confirm 95% of the initial WFC3 grism redshifts in the 38 cases where we detect lines in FORS2 spectroscopy. However, for targets which exhibited a single emission line in WFC3, up to 65% at $z<1.28$ did not have expected emission lines detected in FORS2 and hence may be spurious (although this false-detection rate improves to 33% using the latest public WISP emission line catalogue). From the Balmer decrement the extinction of the WISP galaxies is consistent with $A($H$\alpha)=1$mag. From SED fits to multi-band photometry including Spitzer $3.6\mu$m, we find a median stellar mass of $\log_{10}(M/M_{\odot})=8.94$. Our emission-line-selected galaxies tend to lie above the star-forming main sequence (i.e. higher specific star formation rates). Using [OIII], [OII] and H$\beta$ lines to derive gas-phase metallicities, we find typically sub-solar metallicities, decreasing with redshift. Our WISP galaxies lie below the $z=0$ mass-metallicity relation, and galaxies with higher star formation rates tend to have lower metallicity. Finally, we find a strong evolution with redshift of the H$\alpha$ rest-frame equivalent width, $EW_0($H$\alpha)\propto (1+z)^{1.85\pm0.24}$ with higher $EW_0$ galaxies having larger [OIII]/H$\beta$ and O32 ratios on average, suggesting lower metallicity or higher ionisation parameter in these extreme emission line galaxies.

Polarimetric features during the prompt phase of Gamma-ray Bursts (GRBs) have been essential for elucidating the debated emission mechanisms and gaining insight into the inner structure of GRBs. However, the potential impact of photon-Axion-Like-Particle (ALP) mixing in extragalactic magnetic fields, leading to significant modifications to the initial polarization state, has been overlooked in discussions concerning prompt phase constraints. In this work, we first examine the statistical characteristics of linear polarization degree ($\Pi_{L}$) in GRBs, by utilizing data from polarimetric missions focusing on sub-MeV emissions. Our analysis, conducted with a restricted sample of GRBs spanning various redshifts, reveals a diverse distribution of $\Pi_{L}$, which currently shows no correlation with the GRBs' spectral parameters or properties of candidate host galaxies. We then explore alternations to the initial $\Pi_{L}$ due to photon-ALP mixing within a domain-like structure of the intergalactic magnetic field (${\bf B}_{\rm IGM} $). With the existence of ALPs with $m_{a}$$~$$\lesssim$$~$$10^{-14}$$~$eV and $g_{a\gamma}~$$\simeq$$~0.5\times10^{-11}$, the mixing leads to a decrease in the polarization degree of initially fully linearly polarized photons, while it induces a certain degree of polarization to initially unpolarized photons. To ensure that the effect of mixing is small enough to be negligible, the mixing term $\Delta_{a\gamma} \equiv 1/2\ g_{a\gamma} {\bf B}_{\rm IGM}$ should be less than $1.5\times 10^{-4}$ Mpc$^{-1}$. Currently, the number of GRBs with both sub-MeV polarization measurement and redshift confirmation remains very limited. Certification of redshift for GRBs with low $\Pi_{L}$ would further constrain the parameter space of ALPs or provide an independent means to determine the upper limit on ${\bf B}_{\rm IGM}$.

GG Tau is one of the most studied multiple young stellar systems: GG Tau A is a hierarchical triple surrounded by a massive disc and its companion, GG Tau B, is also a binary. Despite numerous observational attempts, an understanding of the geometry of the GG Tau A system is still elusive. We provide new astrometric measures of the system and we run a set of hydrodynamical simulations with two representative orbits to test how they impact a disc composed of dust and gas. We test the dynamical evolution of the two scenarios on short and long timescales. We obtain synthetic flux emission from our simulations and we compare them with 1300 $\mu$m ALMA dust continuum emission and 1.67 $\mu$m SPHERE dust scattering images to infer the most likely orbital arrangement. We extend the analysis of the binary orbital parameters using six new epochs from archival data, showing that the current measurements alone are not capable of breaking the degeneracy between families of coplanar and misaligned orbits. We found that the time-scale for the onset of the disc eccentricity growth, $\tau_{ecc}$, is a fundamental time-scale for the morphology of the system. Results from numerical simulations show that the best match between is obtained with the misaligned configuration ($\Delta\theta= 30^\circ$) on timescales shorter than $\tau_{ecc}$. The results exhibit an almost circular cavity and dust ring. However, for both scenarios, the cavity size and its eccentricity quickly grow for timescales longer than $\tau_{ecc}$ and the models do not reproduce the observed morphology anymore. This implies that either the age of the system is shorter than $\tau_{ecc}$ or that the disc eccentricity growth is not triggered or dissipated. This finding raises questions on the future evolution of the GG Tau A system and, more in general, on the time evolution of eccentric binaries and their circumbinary discs.

We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided between direct and indirect impact. For the most massive models ($60$ and $85$ $M_\odot$) with $\dot M\gtrsim2\times10^{-7}$ $M_\odot$ yr$^{-1}$, the impact is direct because lower mass loss make stars remove less envelope and therefore remain more massive and less chemically enriched at their surface at the end of their MS phase. For the less massive models ($20$ and $25$ $M_\odot$) with $\dot M\lesssim2\times10^{-8}$ $M_\odot$ yr$^{-1}$, the impact is indirect because lower mass loss make the stars keep high rotational velocities for a longer period of time, then extending the H-core burning lifetime and reaching the end of the MS with higher surface enrichment. Given that the conditions at the H-depletion change, the stars will lose more mass during their He-core burning stages anyways. For $M_\text{zams}=20$ to $40$ $M_\odot$, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to BSG and finally RSG, with larger mass fractions of helium compared to old evolution models. New models also set down to $M_\text{zams}=85\,M_\odot$ the minimal initial mass required for a single star to become WR at metallicity $Z=0.002$. New values for $\dot M$ need to be complemented with upgrades in additional features such as convective core overshooting and distribution of rotational velocities, besides more detailed observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low metallicity environments.

We present a study of CH$_3$CH$_2$CCH, CH$_3$CH$_2$CN, CH$_2$CHCCH, and CH$_2$CHCN in TMC-1 using the QUIJOTE$^1$ line survey. We confirm the presence of CH$_3$CH$_2$CCH in TMC-1, which was previously reported as tentative by our group. From a detailed study of the ethynyl and cyanide derivatives of CH$_2$CH$_2$ and CH$_3$CH$_3$ in TMC-1, we found that the CH$_2$CHCCH/CH$_2$CHCN and CH$_3$CH$_2$CCH/CH$_3$CH$_2$CN abundance ratios are 1.5$\pm$0.1 and 4.8$\pm$0.5, respectively. The derived CH$_2$CHCCH/CH$_3$CH$_2$CCH abundance ratio is 15.3$\pm$0.8, and that of CH$_2$CHCN over CH$_3$CH$_2$CN is 48$\pm$5. All the single substituted isotopologs of vinyl cyanide have been detected, and we found that the first and second carbon substitutions in CH$_2$CHCN provide a $^{12}$C/$^{13}$C ratio in line with that found for other three-carbon bearing species such as HCCNC and HNCCC. However, the third $^{13}$C isotopolog, CH$_2$CH$^{13}$CN, presents an increase in its abundance similar to that found for HCCCN. Finally, we observed eight $b$-type transitions of CH$_2$CHCN, and we find that their intensity cannot be fitted adopting the dipole moment $\mu_b$ derived previously. These transitions involve the same rotational levels as those of the $a$-type transitions. From their intensity, we obtain $\mu_b$=0.80$\pm$0.03\,D, which is found to be in between earlier values derived in the laboratory using intensity measurements or the Stark effect. Our chemical model indicates that the abundances of CH$_3$CH$_2$CCH, CH$_3$CH$_2$CN, CH$_2$CHCCH, and CH$_2$CHCN observed in TMC-1 can be explained in terms of gas-phase reactions.

We report the discovery of 20 new dwarf galaxies in the Local Volume identified as optical counterparts to the Five-hundred-meter Aperture Spherical radio Telescope (FAST) All Sky HI Survey (FASHI) sources. The galaxies have a median stellar mass of $7.8\times 10^6~M_{\odot}$ and a median HI mass of $1.0\times 10^7~M_{\odot}$. Most of them are field galaxies while three are probable members of M 101 and M 106 groups. We also found seven FASHI radio sources to be probable dark HI clouds in nearby groups. Together with four other known HI clouds in the local groups, their mean-square radial velocity difference of 49 km s$^{-1}$ with respect to the host galaxies yields an average total mass of these groups to be ($2.7\pm1.0)\times10^{11}~M_{\odot}$ on the projected scale of 90 kpc.

Hub-filament systems (HFSs) being the potential sites of formation of star clusters and high mass stars, provide a test bed for the current theories that attempt to explain star formation globally. It is thus important to study a large number of HFSs using both intensity and velocity information to constrain these objects better observationally. We present here a study of the hub-filament system associated with G6.55-0.1 using newly obtained observations of radio continuum and $J$=2--1 transition of CO, $^{13}$CO, and C$^{18}$O. The radio continuum maps show multiple peaks that coincide with far-infrared dust continuum peaks indicating the presence of more than one young massive stars in the hub of the HFS. We used the velocity information from the C$^{18}$O(2--1) map to (a) show that the source G6.55-0.1 is not physically associated with the SNR W28 and (b) disentangle and identify the velocity components genuinely associated with G6.55-0.1. Among the velocity-coherent structures identified, the two filaments at 13.8 and 17.3 km s$^{-1}$ contribute a total mass accretion rate of $\sim$3000 M$_{\odot}$ Myr$^{-1}$ to the hub. Both the filaments also show V-shaped structure, characteristic of gravitational collapse, in their velocity profile at the location of the hub. Estimated mass per unit length of the segments of the filaments are smaller than the critical line masses derived from virial equilibrium considerations. This suggests that while the filaments are not gravitationally collapsing as a whole, the spectra from the hub indicate that the inner parts are dynamically decoupled and collapsing to form stars.

We study Monge-Amp\`ere gravity (MAG) as an effective theory of cosmological structure formation through optimal transport theory. MAG is based on the Monge-Amp\`ere equation, a nonlinear version of the Poisson equation, that relates the Hessian determinant of the potential to the density field. We explain how MAG emerges from a conditioned system of independent and indistinguishable Brownian particles, through the large deviation principle, in the continuum limit. To numerically explore this highly non-linear theory, we develop a novel N-body simulation method based on semi-discrete optimal transport. Our results obtained from the very first N-body simulation of Monge-Amp\`ere gravity with over 100 millions particles show that on large scales, Monge-Amp\`ere gravity is similar to the Newtonian gravity but favours the formation of anisotropic structures such as filaments. At small scales, MAG has a weaker clustering and is screened in high-density regions. Although here we study the Monge-Amp\`ere gravity as an effective rather than a fundamental theory, our novel highly-performant optimal transport algorithm can be used to run high-resolution simulations of a large class of modified theories of gravity, such as Galileons, in which the equations of motion are second-order and of Monge-Amp\`ere type.

A substantial fraction of quasars display broad absorption lines (BALs) in their rest-frame ultraviolet spectra. While the origin of BALs is thought to be related to the accretion disc wind, it remains unclear whether the observed ratio of BAL to non-BAL quasars is due to orientation. We conducted observations of 48 BAL quasars and the same number of non-BAL quasars at 322 MHz using the Giant Metrewave Radio Telescope. Combined with previous flux measurements ranging from MHz to GHz frequencies, we compared continuum radio spectra between the two quasar groups. These data offer insights into low-frequency radio properties that have been difficult to investigate with previous observations only at GHz frequencies. Our results present that $73\pm13$ per cent of the BAL quasars exhibit steep or peaked spectra, a higher proportion than $44 \pm 14$ per cent observed in the non-BAL quasars. In contrast, there are no discernible differences between the two quasar groups in the radio luminosity, peak frequency, and spectral index distributions of sources with steep or peaked spectra and sources with flat or inverted spectra. Generally, as the jet axis and line of sight become closer to parallel, quasars exhibit flat or inverted spectra rather than steep or peaked ones. Therefore, these results suggest that BAL quasars are more frequently observed farther from the jet axis than non-BAL quasars. However, given that a certain proportion of BAL quasars exhibit flat or inverted spectra, more than the simple orientation scenario is required to elucidate the radio properties of BAL quasars.

The Gaia satellite has provided the community with three releases containing astrometrical and photometric data as well as by products, such as stellar parameters and variability indicators. By selecting in the Gaia database, one can select stars with the requested characteristics, such as high speed. At present any selection is based on available Gaia releases including a subset of the observations. This, for some stars, can show some limitations, for example there is still not a sufficient number of observations to detect binarity. We investigated a star selected in Gaia EDR3 for its high speed that appears unbound to the Galaxy. We requested high-quality spectra to derive more information on the star. {From the spectroscopic investigation we confirm the low metallicity content of the star, and we derive a detailed chemical composition. The star is poor in carbon and very rich in oxygen: [(C+N+O)/Fe]=+0.65. From the two spectra observed we conclude that the star is in a binary system and from the investigation of the ionisation balance we derive that the star is closer than implied by the Gaia DR3 parallax, and thus has a a lower intrinsic luminosity. The star is probably still unbound, but there is the possibility that it is bound to the Galaxy. Its low carbon abundance suggests that the star was formed in a dwarf galaxy.

The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called ``torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common. Contrary to theories suggesting deep origins for these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the Sun's outer 5-10 percent, the ``Near-Surface Shear Layer". Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the Magneto-Rotational Instability renowned in accretion-disk theory and observed in laboratory experiments. Together, these two facts prompt the general question: Is it possible that the solar dynamo is a near-surface instability? Here, we report strong affirmative evidence in stark contrast to traditional paradigms focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes. State-of-the-art numerical simulations corroborate these estimates and, strikingly, reproduce hemispherical magnetic current helicity laws. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, impacting Earth's electromagnetic infrastructure.

The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a temperature of 2510 $\pm$ 50 K and a metallicity of $0.13\substack{+0.12 \\ -0.11}$. Recent work suggests this system is associated with the Hyades, placing it an older age than previous estimates. Both metallicities agree with current $[Fe/H]$ Hyades measurements (0.11 -- 0.21). Using stellar evolutionary models, we obtain significantly higher masses for the objects, of 0.30 $\pm$ 0.15 $M_\odot$ and 0.08 $\pm$ 0.01 $M_\odot$ (84 $\pm$ 11 $M_{Jup}$) respectively. Using the RVs and a new astrometry point from Keck/NIRC2, we find that the system is likely an edge-on, moderately eccentric ($0.41\substack{+0.27 \\ -0.08}$) configuration. We also estimate the C/O ratio of both objects using custom grid models, obtaining 0.42 $\pm$ 0.10 (primary) and 0.55 $\pm$ 0.10 (companion). From these results, we confirm that this system most likely went through a binary star formation process in the Hyades. The significant changes in this system's parameters since its discovery highlight the importance of high resolution spectroscopy for both orbital and atmospheric characterization of directly imaged companions.

The star formation rate (SFR) is a crucial astrophysical tracer for understanding the formation and evolution of galaxies, determining the interaction between interstellar medium properties and star formation. The mainstream approach to study the stellar content in galaxies relies on pure stellar population synthesis models. However, these methods fail to account for the contamination of SFR caused by nebular gas radiation. Recent studies have indicated that neglecting nebular radiation contamination appears to be non-negligible in galaxies with intense star-forming activities and at relatively high redshifts, potentially leading to overestimation of stellar masses. However, there is currently limited targeted research, particularly regarding galaxies at higher redshifts (z < 0.5). In this investigation, we employ the BPT diagram to select a sample of 2575 star-formation galaxies (SFG) from the SDSS-DR18 dataset, all within specified signal-to-noise ratio (S/N) ranges. Using the spectroscopic fitting tool FADO, which is capable of excluding nebular radiation contributions in spectral fitting, we conduct a tentative investigation of the SFR of star-forming galaxies in SDSS- DR18 with redshifts z < 0.5. Our results show that 45% of the samples show H{\alpha}{\lambda}6563 obtained from FADO fitting to be smaller than that derived solely from the pure stellar population synthesis model qsofitmore, particularly pronounced between redshifts 0.2 and 0.4. We find that the contribution of nebulae is significant and exhibits an evolutionary trend with redshift. We anticipate that by combining optical and near-infrared spectral data, the influence of nebulae may become more prominent in star-forming galaxies at higher redshifts.

We introduce galmoss, a python-based, torch-powered tool for two-dimensional fitting of galaxy profiles. By seamlessly enabling GPU parallelization, galmoss meets the high computational demands of large-scale galaxy surveys, placing galaxy profile fitting in the LSST-era. It incorporates widely used profiles such as the S\'ersic, Exponential disk, Ferrer, King, Gaussian, and Moffat profiles, and allows for the easy integration of more complex models. Tested on 8,289 galaxies from the Sloan Digital Sky Survey (SDSS) g-band with a single NVIDIA A100 GPU, galmoss completed classical S\'ersic profile fitting in about 10 minutes. Benchmark tests show that galmoss achieves computational speeds that are 6 $\times$ faster than those of default implementations.

For segmented telescopes, achieving fine co-focus adjustment is essential for realizing co-phase adjustment and maintenance, which involves adjusting the millimeter-scale piston between segments to fall within the capture range of the co-phase detection system. CGST proposes using a SHWFS for piston detection during the co-focus adjustment stage. However, the residual piston after adjustment exceeds the capture range of the broadband PSF phasing algorithm$(\pm 30 \mu m) $, and the multi-wavelength PSF algorithm requires even higher precision in co-focus adjustment. To improve the co-focus adjustment accuracy of CGST, a fine co-focus adjustment based on cross-calibration is proposed. This method utilizes a high-precision detector to calibrate and fit the measurements from the SHWFS, thereby reducing the impact of atmospheric turbulence and systematic errors on piston measurement accuracy during co-focus adjustment. Simulation results using CGST demonstrate that the proposed method significantly enhances adjustment accuracy compared to the SHWFS detection method. Additionally, the residual piston after fine co-focus adjustment using this method falls within the capture range of the multi-wavelength PSF algorithm. To verify the feasibility of this method, experiments were conducted on an 800mm ring segmented mirror system, successfully achieving fine co-focus adjustment where the remaining piston of all segments fell within $\pm 15 \mu m$.

Synchrotron observation serves as a fundamental tool for studying magnetic fields in various astrophysical settings, yet its ability to unveil three-dimensional (3D) magnetic fields-including plane-of-the-sky orientation, inclination angle relative to the line of sight, and magnetization-remains largely underexplored. Inspired by the latest insights into anisotropic magnetohydrodynamic (MHD) turbulence, we found that synchrotron emission's intensity structures inherently reflect this anisotropy, carrying detailed information about 3D magnetic fields. Capitalizing on this foundation, we integrate a machine learning approach-Convolutional Neural Network (CNN)-to extract this latent information, thereby facilitating the exploration of 3D magnetic fields. The model is trained on synthetic synchrotron emission maps, derived from 3D MHD turbulence simulations encompassing a range of sub-Alfv\'enic to super-Alfv\'enic conditions. We show that the CNN model is physically interpretable and the CNN is capable of reconstructing 3D magnetic field topology and assessing magnetization. In addition, we test our methodology against noise and resolution effects. We show that this CNN-based approach maintains a high degree of robustness in tracing 3D magnetic fields, even when the low spatial frequencies of the synchrotron image are absent. This renders the method particularly suitable for application to interferometric data lacking single-dish measurements.

A compact source, G0.02467-0.0727, was detected in ALMA \threemm observations in continuum and very broad line emission. The continuum emission has a spectral index $\alpha\approx3.3$, suggesting that the emission is from dust. The line emission is detected in several transitions of CS, SO, and SO$_2$ and exhibits a line width FWHM $\approx160$ \kms. The line profile appears Gaussian. The emission is weakly spatially resolved, coming from an area on the sky $\lesssim1"$ in diameter ($\lesssim10^4$ AU at the distance of the Galactic Center; GC). The centroid velocity is $v_{LSR}\approx40$-$50$ \kms, which is consistent with a location in the Galactic Center. With multiple SO lines detected, and assuming local thermodynamic equilibrium (LTE) conditions, $T_\mathrm{LTE} = 13$ K, which is colder than seen in typical GC clouds, though we cannot rule out low-density, subthermally excited, warmer gas. Despite the high velocity dispersion, no emission is observed from SiO, suggesting that there are no strong ($\gtrsim10~\mathrm{km~s}^{-1}$) shocks in the molecular gas. There are no detections at other wavelengths, including X-ray, infrared, and radio. We consider several explanations for the Millimeter Ultra-Broad Line Object (MUBLO), including protostellar outflow, explosive outflow, collapsing cloud, evolved star, stellar merger, high-velocity compact cloud, intermediate mass black hole, and background galaxy. Most of these conceptual models are either inconsistent with the data or do not fully explain it. The MUBLO is, at present, an observationally unique object.

In this paper, we explore the possibility of detecting M-dwarf flares using data from the Zwicky Transient Facility data releases (ZTF DRs). We employ two different approaches: the traditional method of parametric fit search and a machine learning algorithm originally developed for anomaly detection. We analyzed over 35 million ZTF light curves and visually scrutinized 1168 candidates suggested by the algorithms to filter out artifacts, occultations of a star by an asteroid, and known variable objects of other types. Our final sample comprises 134 flares with amplitude ranging from 0.2 to 4.6 magnitudes, including repeated flares and complex flares with multiple components. Using Pan-STARRS DR2 colors, we also assigned a corresponding spectral subclass to each object in the sample. For 13 flares with well-sampled light curves, we estimated the bolometric energy. Our results show that the ZTF's cadence strategy is suitable for identifying M-dwarf flares and other fast transients, allowing for the extraction of significant astrophysical information from their light curves.

Core-collapse supernova remnants are the nebular leftover of defunct massive stars which have died during a supernova explosion, mostly while undergoing the red supergiant phase of their evolution. The morphology and emission properties of those remnants are a function of the distribution of circumstellar material at the moment of the supernova, the intrisic properties of the explosion, as well as those of the ambient medium. By means of 2.5 dimensional numerical magnetohydrodynamics simulations, we model the long term evolution of supernova remnants generated by runaway rotating massive stars moving into a magnetised interstellar medium. Radiative transfer calculations reveal that the projected non-thermal emission of the supernova remnants decreases with time, i.e. older remnants are fainter than younger ones. Older (80 kyr) supernova remnants whose progenitors were moving with space velocity corresponding to a Mach number M = 1 (v_star = 20 km/s ) in the Galactic plane of the ISM (nISM = 1/cm3 ) are brighter in synchrotron than when moving with a Mach number M = 2 (v_star = 40 km/s ). We show that runaway red supergiant progenitors first induce an asymmetric non thermal 1.4 GHz barrel like synchrotron supernova remnants (at the age of about 8 kyr), before further evolving to adopt a Cygnus loop like shape (at about 80 kyr). It is conjectured that a significative fraction of supernova remnants are currently in this bilateral-to-Cygnus-loop evolutionary sequence, and that this should be taken into account in the data interpretation of the forthcoming Cherenkov Telescope Array (CTA) observatory.

We apply current analytical knowledge on the characteristic mass and linear evolution of miniclusters down to redshift $z=0$ to the hypothetical minicluster distribution of the Milky Way. Using the mass-radius relation and a core-halo relation for stable soliton solutions composed of axion-like particles (ALPs), we connect the galactic minicluster mass distribution to that of their ALP star cores. We consider different temperature evolutions of the ALP field with masses in the range $10^{-12}\,\mathrm{eV} \leq m_a \leq 10^{-3}\,$eV and infer the abundance and properties of QCD axion- and ALP stars in our galaxy. We re-evaluate detection prospects for collisions of neutron stars with both ALP stars and miniclusters as well as relativistic ALP bursts, so-called Bosenovae. Our analysis shows that the collision rates between miniclusters and neutron stars can become as large as $\sim 10^5\,$yr$^{-1}$ galaxy$^{-1}$, but that the fraction of encounters that can lead to resonance between ALP mass and magnetosphere plasma frequency is generally well below $\sim 1\,$yr$^{-1}$ galaxy$^{-1}$, depending on the ALP model. We confirm previous results that merger rates of ALP stars are extremely small $< 10^{-12}\,$yr$^{-1}$ galaxy$^{-1}$, while their host miniclusters can merge much more frequently, up to $\sim 10^3\,$yr$^{-1}$ galaxy$^{-1}$ for the QCD axion. We find that Bosenovae and parametric resonance are much more likely to lead to observable signatures than neutron star encounters. We also suggest that a combination of accretion and parametric resonance can lead to observable radio lines for a wide range of ALP masses $m_a$ and photon-couplings $g_{a\gamma\gamma}$.

The Lagoon Nebula (M8) is host to multiple regions with recent and ongoing massive star formation. With M8-Main and M8 East, two prominent regions of massive star formation have been studied in detail over the past years, while large parts of the nebula have received little attention. These largely unexplored regions comprise a large sample of molecular clumps that are affected by the presence of massive O- and B-type stars. We establish an inventory of species observed towards 37 known molecular clumps in M8 by conducting an unbiased line survey for each clump. For this, we used APEX and the IRAM 30m telescope for pointed on-off observations on the clumps. These observations cover bandwidths of 53GHz and 40GHz in frequency ranges from 210GHz to 280GHz and from 70GHz to 117GHz, respectively. Temperatures are derived from rotational transitions of CH3CN, CH3C2H and para-H2CO. Additional archival data from the Spitzer, Herschel, MSX, APEX, WISE, JCMT and AKARI telescopes are used to derive physical parameters of the dust emission by fitting spectral energy distributions to the observed flux densities. Across the observed M8 region, we identify 346 transitions from 70 different molecular species, including isotopologues. We detect tracers of photo-dissociation regions across all the clumps and 38% of these clumps show signs of star formation. We find that PDR tracers are most abundant in clumps with relatively lower H2 column densities. When comparing M8 clumps to ATLASGAL sources at similar distances, we find them to be slightly less massive and have compatible luminosities and radii. This possibly indicates a fragmentation of the gas caused by the O- and B-type stars. In contrast, dust temperatures of the clumps in M8 are found to be increased by approximately 5K (25%) indicating substantial external heating of the clumps by radiation of the present massive stars.

Recent measurements using catalogues of quasars and radio galaxies have shown that the dipole anisotropy in the large-scale distribution of matter is about twice as large as is expected in the standard $\Lambda$CDM model, indeed in any cosmology based on the Friedman-Lema\^itre-Robertson-Walker (FLRW) metric. This expectation is based on the kinematic interpretation of the dipole anisotropy of the cosmic microwave background, i.e. as arising due to our local peculiar velocity. The effect of aberration and Doppler boosting on the projected number counts on the sky of cosmologically distant objects in a flux-limited catalogue can then be calculated and confronted with observations. This fundamental consistency test of FLRW models proposed by Ellis & Baldwin (1984) was revisited by Dalang & Bonvin (2022) [arXiv:2111.03616] who argued that redshift evolution of the sources can significantly affect the expected matter dipole. In this note we demonstrate that the Ellis & Baldwin test is in fact robust to such effects, hence the > 5$\sigma$ dipole anomaly uncovered by Secrest et al. (2021, 2022) [arXiv:2009.14826, arXiv:2206.05624] remains an outstanding challenge to the $\Lambda$CDM model.

Finite-temperature one-loop renormalization of the Standard Model, coupled with dynamic metrics, is conducted in this study. The entire analysis is coherently carried out by using the refined background field method, applied in the spirit of the Coleman-Weinberg technique. The general form of the propagator, introduced in our previous work to facilitate Feynman diagram computation in a general curved background, proves useful in the presence of time-dependent temperature. Its utilization allows for the analysis of a FLRW background to essentially reduce to that of a constant finite-T flat spacetime. Notably, this method also provides an effective tool for addressing one of the longstanding problems in QCD, the Linde problem. The implications of our findings for cosmology, particularly the cosmological constant problem and Hubble tension, are discussed.

We statistically investigate high-frequency whistler waves (with frequencies higher than $\sim 10$ % of the local elect ron cyclotron frequency) at Earth's bow shock using Magnetospheric Multi-Scale (MMS) spacecraft observations. We focus specifically on the wave power within the shock transition layer, where we expect electron acceleration via stochastic sh ock drift acceleration (SSDA) to occur associated with efficient pitch-angle scattering by whistler waves. We find that the wave power is positively correlated with both the Alfv\'en Mach number in the normal incidence frame $M_{\rm A}$ and in the de Hoffmann-Teller frame $M_{\rm A}/\cos \theta_{Bn}$. The empirical relation with $M_{\rm A}/\cos \theta_{Bn}$ is compared with the theory of SSDA that predicts a threshold wave power proportional to $(M_{\rm A}/\cos \theta_{Bn})^{-2}$. The result suggests that the wave power exceeds the theoretical threshold for $M_{\rm A} / \cos \theta_{Bn} \gtrsim 30-60$, beyond which efficient electron acceleration is expected. This aligns very well with previous statistical analysis of electron acceleration at Earth's bow shock (M. Oka, G eophys.~Res.~Lett., 33, 5, 2006). Therefore, we consider that this study provides further support for SSDA as the mechanism of electron acceleration at Earth's bow shock. At higher-Mach-number astrophysical shocks, SSDA will be able to inject electrons into the diffusive shock acceleration process for subsequent acceleration to cosmic-ray energies.

In this work, we embark on the thermodynamic investigation concerning a family of primary charged black holes within the context of shift and parity symmetric Beyond Horndeski gravity. Employing the Euclidean approach, we derive the functional expression for the free energy and derive the first thermodynamic law, offering a methodology to address the challenge of extracting the thermal quantities in shift-symmetric scalar tensor theories characterized by linear time dependence in the scalar field. Following the formal analysis, we provide some illustrative examples focusing on the thermal evaporation of these fascinating objects.

The GQuEST (Gravity from the Quantum Entanglement of Space-Time) experiment uses tabletop-scale Michelson laser interferometers to probe for fluctuations in space-time. We present an interferometer design featuring a novel photon counting readout method that provides unprecedented sensitivity, as it is not subject to the interferometric standard quantum limit. We evaluate the potential of this design to measure space-time fluctuations motivated by recent `geontropic' quantum gravity models. The accelerated accrual of statistical power offered by the photon counting readout enables GQuEST to detect the predicted quantum gravity phenomena within measurement times at least 100 times shorter than equivalent conventional interferometers. The GQuEST design thus enables a fast and sensitive search for signatures of quantum gravity in a laboratory-scale experiment.

Topological properties of the spectrum of shallow-water waves on a rotating spherical body are established. Particular attention is paid to its spectral flow, i.e. the modes whose frequencies transit between the Rossby and inertia-gravity wavebands as the zonal wave number is varied. Organising the modes according to the number of zeros of their meridional velocity, we conclude that the net number of modes transiting between the shallow-water wavebands on the sphere is null, in contrast with the Matsuno spectrum. This difference can be explained by a miscount of zeros under the $\beta$-plane approximation. We corroborate this result with the analysis of Delplace et al (2017) by showing that the curved metric discloses a pair of degeneracy points in the Weyl symbol of the wave operator, non-existent under the $\beta$-plane approximation, each of them bearing a Chern number $-1$.

We revisit supernova (SN) bounds on a hidden sector consisting of millicharged particles $\chi$ and a massless dark photon. Unless the self-coupling is fine-tuned to be small, rather than exiting the SN core as a gas, the particles form a relativistic fluid and subsequent dark QED fireball, streaming out against the drag due to the interaction with matter. Novel bounds due to excessive energy deposition in the mantle of low-energy SNe can be obtained. The cooling bounds from SN~1987A are unexpectedly not affected in the free-streaming regime. The inclusion of $\chi\bar{\chi}\rightarrow \rm SM$ substantially modifies the constraints in the trapping regime, which can only be treated hydrodynamically. Our results can be adapted to generic sub-GeV self-interacting dark sectors.

From 1998 to 2017, neutron monitors located at an altitude of 4300 m on the Tibetan plateau detected 127 long-duration bursts of high-energy radiation in association with thunderclouds. These bursts typically lasted for 10 to 40 minutes, and 89\% of them occurred between 10:00 and 24:00 local time. They were also found to be more likely to occur at night, especially during 18:00$-$06:00 local time period. The observed diurnal and seasonal variations in burst frequency were consistent with the frequencies of lightning and precipitation on the Tibetan plateau. Based on 19 years of data, the present study suggests that an annual variation in burst frequency has a periodicity of $\sim$16 years and a lag of $\sim$3 years relative to solar activity.

Einstein-aether theory provides a model to test the validity of local Lorentz invariance in gravitational interactions. The speed of gravitational waves as measured from the binary neutron star event GW170817 sets stringent limits on Einstein-aether theory, but only on a combination of the theory's free parameters. For this reason, a significant part of the theory's parameter space remains unconstrained by observations. Motivated by this, we explore the propagation of gravitational waves in Einstein-aether theory over an inhomogeneous background (i.e., gravitational wave lensing) as a potential mechanism to break the degeneracies between the theory's free parameters, and hence enable new constraints on the theory to be obtained. By bringing the field equations into the form of the so-called kinetic matrix and applying a formalism known as the propagation eigenstate framework, we find that the speed of gravitational waves is modified by inhomogeneities in the aether field. However, the modification is common to both gravitational polarizations and vanishes in the limit in which gravitational waves propagate with luminal speed. This lens-dependent gravitational wave speed contrasts with the lens-induced birefringence observed in other theories beyond general relativity, like Horndeski's theory. While the potential to improve tests based on gravitational-wave speed is limited, our formalism sets the basis to fully describe signal propagation over inhomogeneous spacetimes in Einstein-aether theory and other extensions of general relativity.

We examine the Universe's transparency to gamma rays within a Lorentz Invariance Violation (LIV) framework, focusing on photon subluminal quadratic corrections driven by a high-energy scale. Based on an explicit calculation, we provide a new expression for the cross section that overcomes the limitations of previous approaches and refines existing constraints for the LIV scale, while we introduce a new approximation that may be useful in LIV scenarios beyond effective field theory. These improvements appear essential for setting constraints on LIV effects with future observations at ultra-high energies, where previous approximations may fall short.

At high energies, the dynamics of a plasma with charged fermions can be described in terms of chiral magnetohydrodynamics. Using direct numerical simulations, we demonstrate that chiral magnetic waves (CMWs) can produce a chiral asymmetry $\mu_5 = \mu_\mathrm{L} - \mu_\mathrm{R}$ from a spatially fluctuating (inhomogeneous) chemical potential $\mu = \mu_\mathrm{L} + \mu_\mathrm{R}$, where $\mu_\mathrm{L}$ and $\mu_\mathrm{R}$ are the chemical potentials of left- and right-handed electrically charged fermions, respectively. If the frequency of the CMW is less than or comparable to the characteristic growth rate of the chiral dynamo instability, the magnetic field can be amplified on small spatial scales. The growth rate of this small-scale chiral dynamo instability is determined by the spatial maximum value of $\mu_5$ fluctuations. Therefore, the magnetic field amplification occurs during periods when $\mu_5$ reaches temporal maxima during the CMW. If the small-scale chiral dynamo instability leads to a magnetic field strength that exceeds a critical value, which depends on the resistivity and the initial value of $\mu$, magnetically-dominated turbulence is produced. Turbulence gives rise to a large-scale dynamo instability, which we find to be caused by the magnetic alpha effect. Our results have consequences for the dynamics of certain high-energy plasmas, such as the early Universe or proto-neutron stars.

We investigate how vector and isovector interactions can be determined within the density regime of neutron stars, while fulfilling nuclear and astrophysics constrains. We make use of the Chiral Mean Field (CMF) model, a SU(3) nonlinear realization of the sigma model within the mean-field approximation, for the first time within a Bayesian analysis framework. We show that neutron-matter $\chi$EFT constraints at low density are only satisfied if the vector-isovector mixed interaction term is included, e.g., a $\omega^2\rho^2$ term. We also show the behavior of the model with respect to the conformal limit. We demonstrate that the CMF model is able to predict a value for the parameter $d_c$ related to the trace anomaly and its derivative takes values below 0.2 above four times saturation density within a hadronic model that does not include a phase transition to deconfined matter. We compare these effects with results from other (non-chiral) Relativistic Mean Field models to assess how different approaches to incorporating the same physical constraints affect predictions of neutron-star properties and dense matter equations of state. We also include data from the gravitation wave event GW230529 detected by the LIGO-Virgo-Kagra collaboration and the most recent radius measurement of PSR J0437-4715 from the NASA NICER mission. Our analysis reveals that this new NICER measurement leads to an average reduction of approximately $\sim 0.15$ km radius in the posterior of the neutron-star mass-radius relationship.