Papers for Monday, Nov 14 2022

Connecting observations of core-collapse supernova explosions to the properties of their massive star progenitors is a long-sought, and challenging, goal of supernova science. Recently, Barker et al. (2022) presented bolometric light curves for a landscape of progenitors from self-consistent neutrino-driven core-collapse supernova (CCSN) simulations. They find a tight relationship between the plateau luminosity of the Type II-P CCSN light curve and the terminal iron core mass of the progenitor. Remarkably, this allows us to constrain progenitor properties with photometry alone. We analyze a large observational sample of Type II-P CCSN light curves and estimate a distribution of iron core masses using the relationship of \citet{barker:2022}. The inferred distribution matches extremely well with the distribution of iron core masses from stellar evolutionary models, and namely, contains high-mass iron cores that suggest contributions from very massive progenitors in the observational data. We use this distribution of iron core masses to infer minimum and maximum mass of progenitors in the observational data. Using Bayesian inference methods to locate optimal initial mass function parameters, we find M$_{\mathrm{min}}=9.8^{+0.37}_{-0.27}$ and M$_{\mathrm{max}}=24.0^{+3.9}_{-1.9}$ solar masses for the observational data.

The large total infrared (TIR) luminosities ($L_{\rm TIR} \gtrsim 10^{12}~L_\odot$) observed in $z \sim 6$ quasars are generally converted into high star formation rates ($SFR \gtrsim 10^2~M_\odot$ yr$^{-1}$) of their host galaxies. However, these estimates rely on the assumption that dust heating is dominated by stellar radiation, neglecting the contribution from the central Active Galactic Nuclei (AGN). We test the validity of this assumption by combining cosmological hydrodynamic simulations with radiative transfer calculations. We find that, when AGN radiation is included in the simulations, the mass (luminosity)-weighted dust temperature in the host galaxies increases from $T\approx 50$ K ($T \approx 70$ K) to $T\approx 80$ K ($T\approx 200$ K), suggesting that AGN effectively heat the bulk of dust in the host galaxy. We compute the AGN-host galaxy $SFR$ from the synthetic spectral energy distribution by using standard $SFR - L_{\rm TIR}$ relations, and compare the results with the "true" values in the simulations. We find that the $SFR$ is overestimated by a factor of $\approx 3$ ($\gtrsim 10$) for AGN bolometric luminosities of $L_{\rm bol} \approx 10^{12}~L_\odot$ ($\gtrsim 10^{13}~ L_\odot$), implying that the star formation rates of $z\sim 6$ quasars can be overestimated by over an order of magnitude.

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions.

Extreme-mass-ratio inspirals (EMRIs) and intermediate-mass-ratio inspirals (IMRIs) are important gravitational-wave (GW) sources for the Laser Interferometer Space Antenna (LISA). So far, their formation and evolution are considered to be independent, but recent theories suggest that stellar-mass black holes (sBHs) and intermediate-mass black hole (IMBHs) can coexist in the accretion disk of an active galactic nucleus (AGN), which indicates that EMRIs and IMRIs may form in the same place. Motivated by the fact that a gas giant migrating in a protoplanetary disk could trap planetesimals close to its orbit, we study in this paper a similar interaction between a gap-opening IMBH in an AGN disk and the sBHs surrounding it. We analyse the torques imposed on the sBHs by the disk as well as by the IMBH, and show that the sBHs can be trapped by the IMBH if they are inside the orbit of the IMBH. Then we implement the torques in our numerical simulations to study the migration of an outer IMBH and an inner sBH, both embedded in an AGN disk. We find that their migration is synchronized until they reach a distance of about ten Schwarzschild radii from the central supermassive black hole, where the pair breaks up due to strong GW radiation. This result indicates that LISA may detect an EMRI and an IMRI within several years from the same AGN. Such a GW source will bring rich information about the formation and evolution of sBHs and IMBHs in AGNs.

In the $\Lambda$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M$_{*}$$\gtrsim$10$\times$$^{7}$M$_{\odot}$, enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that, the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [$\alpha$/Fe]-[Fe/H] compositions of the stellar debris of such mergers $at$ $present$ $day$. These parameters also govern the resulting dynamical state of an accreted galaxy at $z=0$, leading to the expectation that the inner regions of the stellar halo (R$_{\mathrm{GC}}$ $\lesssim$30 kpc) should contain fully phase-mixed debris from both lower and higher mass satellites. In addition, we find that a significant fraction of the lower mass satellites accreted at early times deposit debris in the outer halo (R$_{\mathrm{GC}}$ $>$50 kpc) that are $not$ fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen $in$ $situ$ in higher redshift surveys.

We study systematics associated with estimating simple stellar population (SSP) parameters -- age, metallicity [M/H], $\alpha$-enhancement [$\alpha$/Fe] and IMF shape -- and associated $M_*/L$ gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs) and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H$_\beta$, $\langle$Fe$\rangle$, TiO$_{\rm 2SDSS}$ and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius $R_e$: (a) ages are younger and [$\alpha$/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about $R_e/2$; (c) this makes $M_*/L$ about $2\times$ larger in the central regions than beyond $R_e/2$. While the models of Conroy et al. (2018) return similar [M/H] and [$\alpha$/Fe] profiles, the age and (hence) $M_*/L$ profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs and S0s: younger ages and less enhanced [$\alpha$/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred $M_*/L$ gradients significantly.

We present new observations of [CII] fine structure line emission from an isolated molecular cloud using the upGREAT instrument onboard SOFIA. These data are analyzed together with archival CO=1-0 and HI 21 cm emission spectra to investigate the role of converging atomic gas flows in the formation of molecular clouds. Bright [CII] emission is detected throughout the mapped area that likely originates from photodissociation regions excited by UV radiation fields produced by newborn stars within the cloud. Upon spatial averaging of the [CII] spectra, we identify weak [CII] emission within velocity intervals where the HI 21 cm line is brightest; these are blue-shifted relative to velocities of the CO and bright [CII] emission by 4 km/s. The brightness temperatures, velocity dispersions, and alignment with HI 21 cm velocities connect this [CII] emission component to the cold, neutral atomic gas of the interstellar medium (CNM). We propose that this CNM feature is an accretion flow onto the far--side of the existing molecular cloud. The mass infall rate is 3.2x10**{-4} Msun/yr. There is no direct evidence of a comparable red--shifted component in the [CII] or HI 21 cm spectral lines that would indicate the presence of a converging flow.

The contemporaneous detection of gravitational waves and gamma rays from the GW170817/GRB 170817A, followed by kilonova emission a day after, confirmed compact binary neutron-star mergers as progenitors of short-duration gamma-ray bursts (GRBs), and cosmic sources of heavy r-process nuclei. However, the nature (and lifespan) of the merger remnant and the energy reservoir powering these bright gamma-ray flashes remains debated, while the first minutes after the merger are unexplored at optical wavelengths. Here, we report the earliest discovery of bright thermal optical emission associated with the short GRB 180618A with extended gamma-ray emission, with ultraviolet and optical multicolour observations starting as soon as 1.4 minutes post-burst. The spectrum is consistent with a fast-fading afterglow and emerging thermal optical emission at 15 minutes post-burst, which fades abruptly and chromatically (flux density $F_{\nu} \propto t^{-\alpha}$, $\alpha=4.6 \pm 0.3$) just 35 minutes after the GRB. Our observations from gamma rays to optical wavelengths are consistent with a hot nebula expanding at relativistic speeds, powered by the plasma winds from a newborn, rapidly-spinning and highly magnetized neutron star (i.e. a millisecond magnetar), whose rotational energy is released at a rate $L_{\rm th} \propto t^{-(2.22\pm 0.14)}$ to reheat the unbound merger-remnant material. These results suggest such neutron stars can survive the collapse to a black hole on timescales much larger than a few hundred milliseconds after the merger, and power the GRB itself through accretion. Bright thermal optical counterparts to binary merger gravitational wave sources may be common in future wide-field fast-cadence sky surveys.

Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGN) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure dominated circum-galactic medium (CGM). But explicit CR-magneto-hydrodynamics (CR-MHD) remains computationally expensive, and it is not clear whether it even makes physical sense in simulations that do not explicitly treat magnetic fields or resolved ISM phase structure. We therefore present an intentionally extremely-simplified 'sub-grid' model for CRs, which attempts to capture the key qualitative behaviors of greatest interest for those interested in simulations or semi-analytic models including some approximate CR effects on galactic (>kpc) scales, while imposing negligible computational overhead. The model is numerically akin to some recently-developed sub-grid models for radiative feedback, and allows for a simple constant parameterization of the CR diffusivity and/or streaming speed; it allows for an arbitrary distribution of sources (proportional to black hole accretion rates or star-particle SNe rates or gas/galaxy star formation rates), and interpolates between the limits where CRs escape the galaxies with negligible losses and those where CRs lose most of their energy catastrophically before escape (relevant in e.g. starburst galaxies). The numerical equations are solved trivially alongside gravity in most codes. We compare this to explicit CR-MHD simulations and discuss where the (many) sub-grid approximations break down, and what drives the major sources of uncertainty.

[Abridged] As is well known, in order to generate magnetic fields of observed amplitudes during inflation, the conformal invariance of the electromagnetic field has to be broken by coupling it either to the inflaton or to the scalar curvature. Couplings to scalar curvature pose certain challenges even in slow roll inflation and it seems desirable to consider couplings to the inflaton. It can be shown that, in slow roll inflation, to generate nearly scale invariant magnetic fields of adequate strengths, the non-conformal coupling to the inflaton has to be chosen specifically depending on the inflationary model at hand. In a recent work, we had found that, when there arise sharp departures from slow roll inflation leading to strong features in the scalar power spectra, there inevitably arise sharp features in the spectra of the electromagnetic fields, unless the non-conformal coupling functions are extremely fine tuned. In particular, we had found that, if there occurs an epoch of ultra slow roll inflation, then the strength of the magnetic field over large scales can be severely suppressed. In this work, we examine whether these challenges can be circumvented in models of inflation involving two fields. We show that the presence of the additional scalar field allows us to construct coupling functions that lead to magnetic fields of required strengths even when there arise intermediate epochs of ultra slow roll inflation. However, we find that the features in the spectra of the magnetic fields that are induced due to the departures from slow roll inflation cannot be completely ironed out. We make use of the code MagCAMB to calculate the effects of the magnetic fields on the anisotropies in the cosmic microwave background and investigate if the spectra with features are broadly consistent with the current constraints.

Current ground-based cosmological surveys, such as the Dark Energy Survey (DES), are predicted to discover thousands of galaxy-scale strong lenses, while future surveys, such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will increase that number by 1-2 orders of magnitude. The large number of strong lenses discoverable in future surveys will make strong lensing a highly competitive and complementary cosmic probe. To leverage the increased statistical power of the lenses that will be discovered through upcoming surveys, automated lens analysis techniques are necessary. We present two Simulation-Based Inference (SBI) approaches for lens parameter estimation of galaxy-galaxy lenses. We demonstrate the successful application of Neural Posterior Estimation (NPE) to automate the inference of a 12-parameter lens mass model for DES-like ground-based imaging data. We compare our NPE constraints to a Bayesian Neural Network (BNN) and find that it outperforms the BNN, producing posterior distributions that are for the most part both more accurate and more precise; in particular, several source-light model parameters are systematically biased in the BNN implementation.

The transport of non-thermal particles across a large-scale magnetic field in the presence of magnetised turbulence has been a long-standing issue in high-energy astrophysics. Of particular interest is the dependence of the parallel and perpendicular mean free paths $\lambda_{\parallel}$ and $\lambda_{\perp}$ on rigidity $\mathcal{R}$. We have revisited this important issue with a view to applications from the transport of Galactic cosmic rays to acceleration at astrophysical shocks. We have run test particle simulations of cosmic ray transport in synthetic, isotropic Kolmogorov turbulence at unprecedentedly low reduced rigidites $r_g/L_c \simeq 10^{-4}$, corresponding to $\mathcal{R} \simeq 10 \, \text{TV}$ for a turbulent magnetic field of $B_{rms} = 4 \, \mu\text{G}$ and correlation length $L_c = 30 \, \text{pc}$. Extracting the (asymptotic) parallel and perpendicular mean free paths $\lambda_{\parallel}$ and $\lambda_{\perp}$, we have found $\lambda_{\parallel} \propto (r_g/L_c)^{1/3}$ as expected for a Kolmogorov turbulence spectrum. In contrast, $\lambda_{\perp}$ has a faster dependence on $r_g/L_c$ for $10^{-2} \lesssim r_g/L_c \lesssim 1$, but for $r_g/L_c \ll 10^{-2}$, also $\lambda_{\perp} \propto (r_g/L_c)^{1/3}$. Our results have important implications for the transport of Galactic cosmic rays.

The transport of high-energy particles in the presence of small-scale, turbulent magnetic fields is a long-standing issue in astrophysics. Analytical theories disagree with numerical simulations at rigidities where the particles' gyroradii are slightly smaller than the correlation length of turbulence. At the same time, extending the numerical simulations to lower rigidities has proven computationally prohibitive. In this letter, we provide a solution to the problem of perpendicular transport in isotropic turbulence at both, high and low rigidities. We also clarify the relation between the perpendicular diffusion of particles and the transport of magnetic field lines. To this end, we have run a large suite of test particle simulations at unprecedentedly low rigidites, making extensive use of graphical processing units (GPUs). We have also developed an analytical model, based on (1) initial particle transport along field lines, (2) the transport of field lines and (3) the eventual decorrelation of particles from field lines. Our numerical results exhibit a non-standard rigidity-dependence for the perpendicular diffusion coefficient at intermediate rigidites. At the lowest rigidities, the standard rigidity-dependence is recovered. The simulated diffusion coefficients are nicely reproduced by our analytical model. We have traced the non-standard rigidity-dependence to a subdiffusive phase in the field line transport. Our study has important implications for the transport of Galactic cosmic rays, acceleration at perpendicular shocks and for high-energy particles in the heliosphere.

Seventy-five billion low-mass stars in our galaxy host at least one small planet in their habitable zone (HZ). The stellar ultraviolet (UV) radiation received by the planets is strong and highly variable, and has consequences for atmospheric loss, composition, and habitability. SPARCS is a NASA-funded mission to characterize the quiescent and flare UV emission from low-mass stars, by observing 10 to 20 low-mass stars, over timescales of days, simultaneously in two UV bands: 153-171 nm and 260-300 nm. SPARCS Sun-synchronous terminator orbit allows for long periods of uninterrupted observations, reaching 10s of days for some targets. The payload consists of a 10 cm-class telescope, a dichroic element, UV detectors and associated electronics, a thermal control system, and an on-board processor. The payload is hosted on a Blue Canyon Technologies 6U CubeSat. SPARCS hosts several technology innovations that have broad applicability to other missions. The payload demonstrates the use of "2D-doped" (i.e., delta- and superlattice-doped) detectors and detector-integrated metal dielectric filters in space. This detector technology provides ~5x larger quantum efficiency than NASA's GALEX detectors. In addition, SPARCS' payload processor provides dynamic exposure control, automatically adjusting the exposure time to avoid flare saturation and to time-resolve the strongest stellar flares. A simple passive cooling system maintains the detector temperature under 238K to minimize dark current. The spacecraft bus provides pointing jitter smaller than 6", minimizing the impact of flat-field errors, dark current, and read-noise. All these elements enable competitive astrophysics science within a CubeSat platform. SPARCS is currently in the final design and fabrication phase (Phase C in the NASA context). It will be launched in 2024, for a primary science mission of one year.

We discuss the influence of extragalactic magnetic fields on the intensity of gamma-ray emission produced in electromagnetic cascades from ultra-high energy cosmic rays propagating in extragalactic space. Both cosmic rays and cascade particles propagate mostly out of galaxies, galactic clusters, and large scale structures as their relative volume is small. Therefore their magnetic fields weakly affect emission produced in cascades. Yet estimates of this influence can be useful searching for dark matter particles, when components of extragalactic gamma-ray background should be known, including cascade gamma-ray emission. To study magnetic field influence on cascade emission we analyze cosmic particle propagation in fields of ~10^(-6) and 10^(-12) G (the former is typical inside galaxies and clusters and the latter is common in voids and outside galaxies and clusters). The calculated spectra of cascade gamma-ray emission are similar in the range of ~10^7-10^9 eV. So analyzing cascade emission in this range it is not necessary to specify models of extragalactic magnetic field.

Coronal mass ejections (CMEs) are the largest type of eruptions on the Sun and the main driver of severe space weather at the Earth. In this study, we implement a force-free spheromak CME description within 3-D magnetohydrodynamic simulations to parametrically evaluate successive interacting CMEs within a representative heliosphere. We explore CME-CME interactions for a range of orientations, launch time variations and CME handedness and quantify their geo-effectiveness via the primary solar wind variables and empirical measures of the disturbance storm time index and subsolar magnetopause standoff distance. We show how the interaction of two moderate CMEs between the Sun and the Earth can translate into extreme conditions at the Earth and how CME-CME interactions at different radial distances can maximise different solar wind variables that induce different geophysical impacts. In particular, we demonstrate how the orientation and handedness of a given CME can have a significant impact on the conservation and loss of magnetic flux, and consequently B$_z$, due to magnetic reconnection with the interplanetary magnetic field. This study thus implicates identification of CME chirality in the solar corona as an early diagnostic for forecasting geomagnetic storms involving multiple CMEs.

The dissipation of tidal inertial waves in planetary and stellar convective regions is one of the key mechanisms that drive the evolution of star-planet/planet-moon systems. In this context, the interaction between tidal inertial waves and turbulent convective flows must be modelled in a realistic and robust way. In the state-of-the-art simulations, the friction applied by convection on tidal waves is modelled most of the time by an effective eddy-viscosity. This approach may be valid when the characteristic length scales of convective eddies are smaller than those of tidal waves. However, it becomes highly questionable in the case where tidal waves interact with potentially stable large-scale vortices, as those observed at the pole of Jupiter and Saturn. They are potentially triggered by convection in rapidly-rotating bodies in which the Coriolis acceleration forms the flow in columnar vortical structures along the direction of the rotation axis. In this paper, we investigate the complex interactions between a tidal inertial wave and a columnar convective vortex. We use a quasi-geostrophic semi-analytical model of a convective columnar vortex. We perform linear stability analysis (LSA) to identify the unstable regime and conduct linear numerical simulations for the interactions between the convective vortex and an incoming tidal inertial wave. We verify that in the unstable regime, an incoming tidal inertial wave triggers the most unstable mode of the vortex leading to turbulent dissipation. For stable vortices, the wave-vortex interaction leads to the momentum mixing while it creates a low-velocity region around the vortex core and a new wave-like perturbation in the form of a progressive wave radiating in the far field. The emission of this secondary wave is the strongest when the wavelength of the incoming wave is close to the characteristic size of the vortex.

We provide a correlation analysis of various signatures associated with traces of the decay of dark matter and galaxy spatial distribution, which play an important role in the context of current and future observations and cosmological constraints. Attention is paid to the constraints that can be obtained for decaying sterile neutrinos when analyzing observations in the context of the Spectr-Roentegn-Gamma (SRG) mission. We examine the correlation spectra of these signatures, which can be obtained both for the eROSITA telescope and for the ART-XC telescope, and investigate the applicability of the multipole approximation.

We search for merger products among the 25 most massive white dwarfs in the Montreal White Dwarf Database 100 pc sample through follow-up spectroscopy and high-cadence photometry. We find an unusually high fraction, 40%, of magnetic white dwarfs among this population. In addition, we identify four outliers in transverse velocity and detect rapid rotation in five objects. Our results show that $56^{+9}_{-10}$\% of the $M\approx1.3~M_{\odot}$ ultramassive white dwarfs form through mergers. This fraction is significantly higher than expected from the default binary population synthesis calculations using the $\alpha$-prescription (with $\alpha \lambda = 2$), and provides further support for efficient orbital shrinkage, such as with low values of the common envelope efficiency.

The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the formation of the BBHs with different formation channels. Of all formation scenarios, the chemically homogeneous evolution (CHE) channel has stood out with distinguishing features, namely, nearly-equal component masses and preferentially high individual spins aligned with the orbital angular momentum. We perform Bayesian inference on the BBH events officially reported in GWTC-3 with astrophysically-predicted priors representing different formation channels of the isolated binary evolution (CEE: common-envelope evolution channel; CHE; SMT: stable mass transfer). Given assumed models, we report strong evidence for GW190517\_055101 being most likely to have formed through the CHE channel. Assuming the BBH events in the subsample are all formed through one of the isolated binary evolution channels, we obtain the lower limits on the local merger rate density of these channels at $11.45 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CEE), $0.18 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CHE), and $0.63 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (SMT) at $90\%$ credible level.

Strong lensing in galaxy clusters probes properties of dense cores of dark matter halos in mass, studies the distant universe at flux levels and spatial resolutions otherwise unavailable, and constrains cosmological models independently. The next-generation large scale sky imaging surveys are expected to discover thousands of cluster-scale strong lenses, which would lead to unprecedented opportunities for applying cluster-scale strong lenses to solve astrophysical and cosmological problems. However, the large dataset challenges astronomers to identify and extract strong lensing signals, particularly strongly lensed arcs, because of their complexity and variety. Hence, we propose a framework to detect cluster-scale strongly lensed arcs, which contains a transformer-based detection algorithm and an image simulation algorithm. We embed prior information of strongly lensed arcs at cluster-scale into the training data through simulation and then train the detection algorithm with simulated images. We use the trained transformer to detect strongly lensed arcs from simulated and real data. Results show that our approach could achieve 99.63 % accuracy rate, 90.32 % recall rate, 85.37 % precision rate and 0.23 % false positive rate in detection of strongly lensed arcs from simulated images and could detect almost all strongly lensed arcs in real observation images. Besides, with an interpretation method, we have shown that our method could identify important information embedded in simulated data. Next step, to test the reliability and usability of our approach, we will apply it to available observations (e.g., DESI Legacy Imaging Surveys) and simulated data of upcoming large-scale sky surveys, such as the Euclid and the CSST.

It has been suggested that a trail of diffuse galaxies, including two dark matter deficient galaxies (DMDGs), in the vicinity of NGC1052 formed because of a high-speed collision between two gas-rich dwarf galaxies, one bound to NGC1052 and the other one on an unbound orbit. The collision compresses the gas reservoirs of the colliding galaxies, which in turn triggers a burst of star formation. In contrast, the dark matter and pre-existing stars in the progenitor galaxies pass through it. Since the high pressures in the compressed gas are conducive to the formation of massive globular clusters (GCs), this scenario can explain the formation of DMDGs with large populations of massive GCs, consistent with the observations of NGC1052-DF2 (DF2) and NGC1052-DF4. A potential difficulty with this `mini bullet cluster' scenario is that the observed spatial distributions of GCs in DMDGs are extended. GCs experience dynamical friction causing their orbits to decay with time. Consequently, their distribution at formation should have been even more extended than that observed at present. Using a semi-analytic model, we show that the observed positions and velocities of the GCs in DF2 imply that they must have formed at a radial distance of 5-10kpc from the center of DF2. However, as we demonstrate, the scenario is difficult to reconcile with the fact that the strong tidal forces from NGC1052 strip the extendedly distributed GCs from DF2, requiring 33-59 massive GCs to form at the collision to explain observations.

We present the largest and most homogeneous collection of near-infrared (NIR) spectra of Type Ia Supernovae (SNe Ia): 339 spectra of 98 individual SNe obtained as part of the Carnegie Supernova Project-II. These spectra, obtained with the FIRE spectrograph on the 6.5 m Magellan Baade telescope, have a spectral range of 0.8-2.5 $\mu$m. Using this sample, we explore the NIR spectral diversity of SNe Ia and construct a template of spectral time series as a function of the light-curve-shape parameter, color stretch sBV. Principal Component Analysis is applied to characterize the diversity of the spectral features and reduce data dimensionality to a smaller subspace. Gaussian process regression is then used to model the subspace dependence on phase and light-curve shape and the associated uncertainty. Our template is able to predict spectral variations that are correlated with sBV , such as the hallmark NIR features: Mg II at early times and the H-band break after peak. Using this template reduces the systematic uncertainties in K-corrections by $\sim$90% compared to those from the Hsiao template (Hsiao 2009). These uncertainties are on the level of 4 $\times$ 10$^{-4}$ mag on average. We have also explored a neural network approach using a conditional variational autoencoder that produces promising results for characterizing supernova spectra, though requires a larger data set to assemble comparable quality. This template can serve as the baseline spectral energy distribution for light-curve fitters and can identify peculiar spectral features that might point to compelling physics. The results presented here will substantially improve future SN Ia cosmological experiments, for both nearby and distant samples.

Hierarchical clustering is a common algorithm in data analysis. It is unique among many clustering algorithms in that it draws dendrograms based on the distance of data under a certain metric, and group them. It is widely used in all areas of astronomical research, covering various scales from asteroids and molecular clouds, to galaxies and galaxy cluster. This paper systematically reviews the history and current status of the development of hierarchical clustering methods in various branches of astronomy. These applications can be grouped into two broad categories, one revealing the intrinsic hierarchical structure of celestial systems and the other classifying large samples of celestial objects automatically. By reviewing these applications, we can clarify the conditions and limitations of the hierarchical clustering algorithm, and make more reasonable and reliable astronomical discoveries.

The TYPHOON program is producing an atlas of spectroscopic data cubes of 44 large-angular-sized galaxies with complete spatial coverage from 3650-9000 A. This survey provides an unparalleled opportunity to study variations in the interstellar medium (ISM) properties within individual HII regions across the entire star-forming disks of nearby galaxies. This can provide key insights into the spatial distribution and resolved properties of the ISM to understand how efficiently metals are mixed and redistributed across spirals and dwarf galaxies. In this Proceeding, we present early science results from six nearby spiral galaxies as part of the TYPHOON program from Grasha et al. (2022). We use HIIPhot to identify the HII regions within the galaxy based on the surface brightness of the H-alpha emission line and measure variations of the HII region oxygen abundance. In this initial work, we find that while the spiral pattern plays a role in organizing the ISM, it alone does not establish the relatively uniform azimuthal variations we observe across all the galaxies. Differences in the metal abundances are more likely driven by the strong correlations with the local physical conditions. We find a strong and positive correlation between the ionization parameter and the local abundances as measured by the relative metallicity offset $\Delta$(O/H), indicating a tight relationship between local physical conditions and their localized enrichment of the ISM. These variations can be explained by a combination of localized, star formation-driven self-enrichment and large-scale mixing-driven dilution due to the passing of spiral density waves.

We explain the overall equilibrium-temperature-dependent trend in the exoplanet mass-radius diagram, using the escape mechanisms of hydrogen and relevant volatiles, and the chemical equilibrium calculation of molecular hydrogen (H$_2$) break-up into atomic hydrogen (H). We identify two Cosmic Hydrogen and Ice Loss Lines (CHILLs) in the mass-radius diagram. Gas disks are well known to disperse in ten million years. However, gas-rich planets may lose some or almost all gas on a much longer timescale. We thus hypothesize that most planets that are born out of a hydrogen-gas-dominated nebular disk begin by possessing a primordial H$_2$-envelope. This envelope is gradually lost due to escape processes caused by host-stellar radiation.

The Population III (Pop III) stars could die as supernovae that chemically enrich the early universe and trigger the primeval galaxy formation. Here, we model the effect of Pop III supernova remnants (SNRs) on the formation of primeval galaxies using high-resolution radiation-hydrodynamical simulations with ENZO. We find that SNRs from a top-heavy Pop~III initial mass function (IMF) produce more metals, which leads to more efficient cooling of gas, thus an earlier episode of Pop~II star formation in the primeval galaxies. These Pop~II stars have metallicities of $\sim 10^{-3}-10^{-2}$ Zsun, and their mass function follows a power law in which indices depend on types of Pop III SNRs. Furthermore, the morphology of our primeval galaxies looks irregular similar to the local dwarf galaxies. These primeval galaxies have bolometric luminosities up to $6.0 \times 10^5 $ Lsun that is possibly detected by the James Webb Space Telescope ({JWST}) if they form at the redshift of $z \sim 10$.

We present a spatially resolved spectroscopic study for the metal poor HII galaxy J084220+115000 using MEGARA Integral Field Unit observations at the Gran Telescopio Canarias. We estimated the gas metallicity using the direct method for oxygen, nitrogen and helium and found a mean value of 12+$\log$(O/H)=$8.03\pm$0.06, and integrated electron density and temperature of $\sim161$ cm$^{-3}$ and $\sim15400$ K, respectively. The metallicity distribution shows a large range of $\Delta$(O/H) = 0.72 dex between the minimum and maximum (7.69$\pm$0.06 and 8.42$\pm$0.05) values, unusual in a dwarf star-forming galaxy. We derived an integrated $\log$(N/O) ratio of $-1.51\pm0.05$ and found that both N/O and O/H correspond to a primary production of metals. Spatially resolved maps indicate that the gas appears to be photoionized by massive stars according to the diagnostic line ratios. Between the possible mechanisms to explain the starburst activity and the large variation of oxygen abundance in this galaxy, our data support a possible scenario where we are witnessing an ongoing interaction triggering multiple star-forming regions localized in two dominant clumps.

This is the first in a series of papers that present sets of different results for 136 compact, known planetary nebulae within a 10 $\times$ 10 degree region of the Galactic bulge. We use a large, previously unpublished sample of our own extensive ESO 8 m VLT deep imaging and spectroscopic data. This is combined with archival deep $\textit{HST}$ imaging where available to provide a detailed morphological classification and study. The influence of angular resolution and sensitivity when assigning a morphology is discussed. A large fraction (68%) of the sample are shown to be bipolars and the implications for this in the context of planetary nebulae progenitors are explored. Four new planetary nebula central stars are also identified which are not in $\textit{Gaia}$. This is based on both VLT and deep archival Pan-STARRS broad-band imagery. Some 11 putative central stars previously reported, based on $\textit{Gaia}$ positions, are also not the true central star. In other cases the $\textit{Gaia}$ central stars reported in the literature are actually based on the overall centroid position of a very compact planetary nebula rather than the actual central star within it. $\textit{Gaia}$ parallax distances and kinematic ages for PNe in this sample are provided where possible based on fresh angular size measures from the new VLT imagery and $\textit{Gaia}$ distances and literature expansion velocities when available. All these results are discussed within the context of the overall characteristics of the Galactic bulge and its planetary nebula population.

We report X-ray observations of the High Mass X-ray Binary (HMXB) pulsar 4U 1907+09. Spectral and Timing analysis of the source has been performed using NuSTAR observation. Timing analysis of the photon events revealed the coherent X-ray pulsation of the source with a pulse period of $442.92\;\pm\;0.03$ s. It is observed that the source is spinning down at a rate of $0.1971(4) s yr^{-1}$. The pulse profile is characterized by a decaying amplitude of the secondary peak and relative growth in the amplitude of the primary peak with an increase in energy. The broad-band spectral coverage of NuSTAR has been used to observe multiple absorption features in the X-ray continuum of the source. We confirm the presence of two prominent cyclotron absorption features at $\sim 17$ keV and $\sim 38$ keV respectively. In addition, we have detected an absorption-line feature at $\sim 8$ keV, with an equivalent width of $\sim 1.3$ keV. The variation of the spectral parameters with pulse phase has been observed using phase-resolved spectroscopy and the relevant variabilities of the parameters have been discussed with the underlying physical implications. The continuum evolution and variations in spectral parameters have also been studied by time-resolved spectroscopy.

Fast radio bursts (FRBs) are millisecond-timescale bursts of coherent radio emission that are luminous enough to be detectable at cosmological distances. In this review I describe the discovery of FRBs, subsequent advances in our understanding of them, and future prospects. Thousands of potentially observable FRBs reach Earth every day; they probably originate from highly magnetic and/or rapidly rotating neutron stars in the distant Universe. Some FRBs repeat, with this sub-class often occurring in highly magnetic environments. Two repeaters exhibit cyclic activity windows, consistent with an orbital period. One nearby FRB was from a Galactic magnetar during an X-ray outburst. The host galaxies of some FRBs have been located, providing information about the host environments and the total baryonic content of the Universe.

Context. Some of of the most prominent sources for energetic particles in our Solar System are huge eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs), which usually drive shocks that accelerate charged particles up to relativistic energies. In particular, energetic electron beams can generate radio bursts through the plasma emission mechanism. The main types of bursts associated with CME shocks are type II and herringbone bursts. However, it is currently unknown where early accelerated electrons that produce metric type II bursts and herringbones propagate and when they escape the solar atmosphere. Aims. Here, we investigate the acceleration location, escape, and propagation directions of electron beams during the early evolution of a strongly expanding CME-driven shock wave associated with herrinbgone bursts. Methods. We used ground-based radio observations from the Nan\c{c}ay Radioheliograph combined with space-based extreme-ultraviolet and white-light observations from the Solar Dynamics Observatory and and the Solar Terrestrial Relations Observatory. We produced a three-dimensional (3D) representation of the electron acceleration locations which, combined with results from magneto-hydrodynamic (MHD) models of the solar corona, was used to investigate the origin of the herringbone bursts observed. Results. Multiple herringbone bursts are found close to the CME flank in plane-of-sky images. Some of these herringbone bursts have unusual inverted J shapes and opposite drifting herringbones also show opposite senses of circular polarisation. By using a 3D approach combined with the radio properties of the observed bursts, we find evidence that the first radio emission in the CME eruption most likely originates from electrons that initially propagate in regions of low Alfv\'en speeds and along closed magnetic field lines forming a coronal streamer.

One of the big challenges in exoplanet science is to determine the atmospheric makeup of extrasolar planets, and to find biosignatures that hint at the existence of biochemical processes on another world. The biomarkers we are trying to detect are gases in the exoplanet atmosphere like oxygen or methane, which have deep absorption features in the visible and near-infrared spectrum. Here we establish the ultimate quantum limit for determining the presence or absence of a spectral absorption line, for a dim source in the presence of a much brighter stellar source. We characterise the associated error exponent in both the frameworks of symmetric and asymmetric hypothesis testing. We found that a structured measurement based on spatial demultiplexing allows us to decouple the light coming from the planet and achieve the ultimate quantum limits. If the planet has intensity $\epsilon \ll 1$ relative to the star, we show that this approach significantly outperforms direct spectroscopy yielding an improvement of the error exponent by a factor $1/\epsilon$. We find the optimal measurement, which is a combination of interferometric techniques and spectrum analysis.

We investigate the formation and eruption of hot channels associated with the M6.5 class flare (SOL2015-06-22T18:23) occurring in NOAA AR 12371 on 2015 June 22. Two flare precursors are observed before the flare main phase. Observations in 94 {\AA} and 131 {\AA} by SDO/AIA have revealed the early morphology of the first hot channel as a group of hot loops, which is termed as seed hot channel. A few seed hot channels are formed above the polarity inversion line (PIL) and the formation is associated with footpoint brightenings' parallel motion along the PIL, which proceeds into the early stage of the flare main phase. During this process, seed hot channels build up and rise slowly, being accelerated at the peak of the second precursor. They merge in the process of acceleration forming a larger hot channel, which then forms an "inverted {\gamma}" shape kinking structure. Before the flare peak, the second kinking hot channel with negative crossing appears near the first kinking hot channel that has erupted. The eruption of these two hot channels produce two peaks on the main flare's GOES light curve. The footpoint brightenings' propagation along the PIL indicate that the first kinking hot channel may be formed due to zipper reconnection. The occurrence of merging between seed hot channels observed by AIA is supported by the extrapolated nonlinear force-free field models. The observed writhing motion of the first kinking hot channel may be driven by the Lorentz force.

We present measurements of the specific angular momentum content (j) of galaxies drawn from the Simba cosmological hydrodynamic simulations. For the stellar, HI and baryonic matter components we demonstrate the existence of extremely tight relations between j and the mass contained within the radius at which the HI mass surface density decreases to 1 Msun/pc^2. These relations are broadly consistent with a variety of empirical measurements. We confirm the observational result that the scatter in the stellar j-M relation is driven largely by HI content, and measure the dependence of its scatter on the deviations of galaxies from other important scaling relations. For a given stellar mass, HI-rich/poor galaxies have more/less-than-average stellar specific angular momentum. A similar, yet weaker, correlation exists for HI mass fraction. Overall, our results demonstrate the utility of the Simba simulations as a platform for understanding and contextualising the data and results from forthcoming large galaxy surveys.

In recent years, Galactic archaeology has become a particularly vibrant field of astronomy, with its main focus set on the oldest stars of our Galaxy. In most cases, these stars have been identified as the most metal-poor. However, the struggle to find these ancient fossils has produced an important bias in the observations - in particular, the intermediate metal-poor stars (-2.5<[Fe/H]< -1.5) have been frequently overlooked. The missing information has consequences for the precise study of the chemical enrichment of our Galaxy, in particular for what concerns neutron-capture elements and it will be only partially covered by future multi-object spectroscopic surveys such as WEAVE and 4MOST. Measuring at Intermediate Metallicity Neutron Capture Elements (MINCE) is gathering the first high-quality spectra (high S/N ratio and high resolution) for several hundreds of bright and metal-poor stars, mainly located in our Galactic halo. We compiled our selection mainly on the basis of Gaia data and determined the stellar atmospheres of our sample and the chemical abundances of each star. In this paper, we present the first sample of 59 spectra of 46 stars. We measured the radial velocities and computed the Galactic orbits for all stars. We found that 8 stars belong to the thin disc, 15 to disrupted satellites, and the remaining cannot be associated to the mentioned structures, and we call them halo stars. For 33 of these stars, we provide abundances for the elements up to zinc. We also show the chemical evolution results for eleven chemical elements, based on recent models. Our observational strategy of using multiple telescopes and spectrographs to acquire high S/N and high-resolution spectra has proven to be very efficient since the present sample was acquired over only about one year of observations. Finally, our target selection strategy proved satisfactory for our purposes.

Following from our recent work, we present results of a detailed analysis of a representative sample of nearby galaxies. The photometric parameters of the morphological components are obtained from bulge-disk decompositions, using GALFIT software. The previously obtained method and library of numerical corrections for dust, decomposition and projection effects, are used to correct the measured (observed) parameters to intrinsic values. Observed and intrinsic galaxy dust and star-formation related scaling relations are presented, to emphasize the scale of the biases introduced by these effects. To understand the extent to which star-formation is distributed in the young stellar disks of galaxies, star-formation connected relations which rely on measurements of scale-lengths and fluxes / luminosities of H$\alpha$ images, are shown. The mean dust opacity, dust-to-stellar mass and dust-to-gas ratios of the sample, together with the main characteristics of the intrinsic relations are found to be consistent with values found in the literature.

Context: Supernovae (SNe) inject large amounts of energy and chemically enriched materials into their surrounding interstellar medium and, in some instances, into molecular clouds (MCs). The interaction of a supernova remnant (SNR) with a MC plays a crucial role in the evolution of the cloud's physical and chemical properties. Despite their importance, only a handful of studies have been made addressing the molecular richness in MCs impacted by SNRs. (Sub)millimter wavelength observations of SNR-MC can be used to build a census of their molecular richness, which in turn can motivate various chemical and physical models aimed at explaining the chemical evolution of the clouds. Aims: We carried out multi-molecule/multi-transition observations toward the region F abutting the SNR W28. We used the detected lines to constrain the physical conditions of this region. Methods: We used the APEX Telescope to observe molecular lines in the frequency from $213\rm{-}374\, \textrm{GHz}$. We used non-LTE RADEX modeling to interpret the observational data. Results: We detected emission from multiple molecular species, namely CH$_{3}$OH, H$_{2}$CO, SO, SiO, CN, CCH, NO, CS, HCO$^+$, HCN, HNC, N$_2$H$^+$, CO, and from isotopologues of some of them. We report the first detection of thermally excited (non-maser) CH$_{3}$OH emission toward a SNR. Employing non-LTE RADEX modeling of multiple H$_{2}$CO and CH$_{3}$OH lines, we constrain the kinetic temperature from 60 to 100$\,$K and the gas density from $9\times 10^{5}$ to $5\times 10^{6}\,\rm cm^{-3}$. We obtained an ortho-para ratio $\sim$2 for H$_{2}$CO, which indicates that formaldehyde is most likely formed on dust grain surfaces. Conclusions: Our results show that molecules as complex as H$_{2}$CO and CH$_{3}$OH can be detected in SNR-MC interactions. This could motivate chemical modelling to explore their formation pathways.

Many stars of different spectral types with planets in the habitable zone are known to emit flares. Until now, studies that address the long-term impact of stellar flares and associated Coronal Mass Ejections (CMEs) assumed that the planet's interior remains unaffected by interplanetary CMEs, only considering the effect of plasma/UV interactions on the atmosphere of planets. Here, we show that the magnetic flux carried by flare-associated CMEs results in planetary interior heating by Ohmic dissipation and leads to a variety of interior--exterior interactions. We construct a physical model to study this effect and apply it to the TRAPPIST-1 star whose flaring activity has been constrained by Kepler observations. Our model is posed in a stochastic manner to account for uncertainty and variability in input parameters. Particularly for the innermost planets, our results suggest that the heat dissipated in the silicate mantle is both of sufficient magnitude and longevity to drive geological processes and hence facilitate volcanism and outgassing of the TRAPPIST-1 planets. Furthermore, our model predicts that Joule heating can further be enhanced for planets with an intrinsic magnetic field compared to those without. The associated volcanism and outgassing may continuously replenish the atmosphere and thereby mitigate the erosion of the atmosphere caused by the direct impact of flares and CMEs. To maintain consistency of atmospheric and geophysical models, the impact of stellar flares and CMEs on atmospheres of close-in exoplanetary systems needs to be studied in conjunction with the effect on planetary interiors.

We present an investigation on escape fractions of UV photons from a unique sample of lensed low-mass emission line selected galaxies at z < 3.5 found in the SHARDS Hubble Frontier Fields medium-band survey. We have used this deep imaging survey to locate 42 relatively low-mass galaxies, down to $log(M_{*}/M_{\odot}) = 7$, between redshifts 2.4 < z < 3.5 which are candidate line emitters. Using deep multi-band Hubble UVIS imaging we investigate the flux of escaping ionizing photons from these systems, obtaining 1$\sigma$ upper limits of $f^{rel}_{esc}$ ~7% for individual galaxies, and < 2% for stacked data. We measure potential escaping Lyman-continuum flux for two low-mass line emitters with values at $f^{\rm rel}_{\rm esc} = 0.032^{+0.081}_{-0.009}$ and $f^{\rm rel}_{\rm esc} = 0.021^{+0.101}_{-0.006}$, both detected at the ~3.2$\sigma$ level. A detailed analysis of possible contamination reveals a < 0.1% probability that these detections result from line-of-sight contamination. The relatively low Lyman-continuum escape fraction limit, and the low fraction of systems detected, is an indication that low-mass line emitting galaxies may not be as important a source of reionization as hoped if these are analogs of reionization sources. We also investigate the structures of our galaxy sample, finding no evidence for a correlation of escape fraction with asymmetric structure.

Capture into mean motion resonance (MMR) is an important dynamical mechanism as it shapes the final architecture of a planetary system. We simulate systems of two or three planets undergoing migration with varied initial parameters such as planetary mass and disk surface density and analyse the resulting resonant chains. In contrast to previous studies, our results show that the disk properties have the dominant impact on capture into mean motion resonance, while the total planetary mass barely affects the final system configuration as long as the planet does not open a gap in the disk. We confirm that the adiabatic resonant capture is the correct framework to understand the conditions leading to MMR formation, since its predictions are qualitatively similar to the numerical results. However, we find that the eccentricity damping can facilitate the capture in a given resonance. We find that under typical disk conditions, planets tend to be captured into 2:1 or 3:2 MMRs, which agrees well with the observed exoplanet MMRs. Our results predict two categories of systems: those that have uniform chains of wide resonances (2:1 or 3:2 MMRs) and those that have a more compact inner pair than the outer pair such as 4:3:2 chains. Both categories of resonant chains are present in observed exoplanet systems. On the other hand, chains with a wider inner pair than the outer one are very rare and emerge from stochastic capture. Our work here can be used to link current configuration of exoplanetary systems to the formation conditions within protoplanetary disks.

We present the first look at the spectral and timing analysis of Narrow-Line Seyfert 1 Galaxies (NLS1s) with the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spectrum-Roentgen-Gamma (SRG) mission. The sample of about 1,200 NLS1s was obtained via cross-match of the first eROSITA All-Sky Survey (eRASS1) catalogue with the catalogue of spectroscopically selected NLS1s from the Sloan Digital Sky Survey (SDSS) DR12 by Rakshit et al. [ApJS, 229, 39 (2017)]. The X-ray spectral analysis is based on a simple power-law fit. The photon index distribution has a mean value of about 2.81+-0.03, as expected from previous X-ray studies of NLS1s. Interestingly, it is positively skewed, and about 10 percent of the sources are located in the super-soft tail of photon indices larger than 4. These sources are of further interest as their source counts run into the X-ray background at values at around 1 keV. We argue that ionised outflows have been detected by eROSITA and may account for some of the extreme spectral steepness, which is supported by correlations found between the photon index and optical outflows parameters. We analysed the intrinsic X-ray variability of the eRASS1 to eRASS3 light curves of the sample but do not find significant variability neither during the individual survey scans nor between them.

Classical T Tauri stars (CTTSs) are young stellar objects that accrete materials from their accretion disc influenced by their strong magnetic field. The magnetic pressure truncates the disc at a few stellar radii and forces the material to leave the disc plane and fall onto the stellar surface by following the magnetic field lines. However, this global scheme may be disturbed by the presence of a companion interacting gravitationally with the accreting component. This work is aiming to study the accretion and the magnetic field of the tight eccentric binary DQ Tau, composed of two equal-mass ($\sim$ 0.6 \msun ) CTTSs interacting at different orbital phases. We investigated the variability of the system using a high-resolution spectroscopic and spectropolarimetric monitoring performed with ESPaDOnS at the CFHT. We provide the first ever magnetic field analysis of this system, the Zeeman-Doppler imaging revealed a stronger magnetic field for the secondary than the primary (1.2 kG and 0.5 kG, respectively), but the small-scale fields analysed through Zeeman intensification yielded similar strengths (about 2.5 kG). The magnetic field topology and strengths are compatible with the accretion processes on CTTSs. Both components of this system are accreting, with a change of the main accretor during the orbital motion. In addition, the system displays a strong enhancement of the mass accretion rate at periastron and apastron. We also discovered, for the first time in this system, the apsidal motion of the orbital ellipse.

We present a homogeneously-selected sample of 15779 candidate binary systems with main sequence primary stars and orbital periods shorter than 5 days. The targets were selected from TESS full-frame image lightcurves on the basis of their tidally-induced ellipsoidal modulation. Spectroscopic follow-up suggests a sample purity of $83 \pm 13$ per cent. Injection-recovery tests allow us to estimate our overall completeness as $28 \pm 3$ per cent with $P_\mathrm{orb} < 3$ days and to quantify our selection effects. We estimate that $39 \pm 4$ per cent of our sample are contact binary systems, and we disentangle the period distributions of the contact and detached binaries. We derive the orbital period distribution of the main sequence binary population at short orbital periods, finding a distribution continuous with the log-normal distribution previously found for solar-type stars at longer periods, but with a significant steepening at $P_\mathrm{orb} \lesssim 3$ days, and a pile-up of contact binaries at $P_\mathrm{orb} \approx 0.4$ days. Companions in the period range 1--5 days are found to be an order of magnitude more frequent around stars hotter than $\approx 6250 K$ (the Kraft break) when compared to cooler stars, suggesting that magnetic braking plays an important role in shaping the temperature and period distributions. We detect resolved tertiary companions to $9.0 \pm 0.2$ per cent of our binaries with a median separation of 3200 AU. The frequency of tertiary companions rises to $29 \pm 5$ per cent among the systems with the shortest ellipsoidal periods. This large binary sample with quantified selection effects will be a powerful resource for future studies of detached and contact binary systems with $P_\mathrm{orb} < 5$ days.

Based on almost 20 years of data collected by the high-resolution spectrometer SPI on board the International Gamma-Ray Astrophysics Laboratory (INTEGRAL) we present constraints on a decaying dark matter particle manifesting itself via a narrow line-like spectral feature. Our ON-OFF type analysis of the Milky Way observations allowed us to constrain the lifetime to be $\gtrsim 10^{20}-10^{21}$ yrs for DM particles with masses $40\,\text{keV}\,<\,M_{\text{DM}}\,<\,14\,\text{MeV}$. Within this mass range our analysis also reveals 32 line-like features detected at $\geq 3\sigma$ significance, 29 of which coincide with known instrumental and astrophysical lines. In particular, we report on the detection of the electron-positron annihilation (511 keV) and $^{26}$Al (1809 keV) lines with spatial profiles consistent with previous results in the literature. For the particular case of the sterile neutrino DM we report the limits on the mixing angle as a function of sterile neutrino mass. We discuss the dominant impact of systematic uncertainties connected to the strongly time-variable INTEGRAL/SPI instrumental background as well as the ones connected to the uncertainties of MW DM density profile measurements on the derived results.

Various WIMP dark matter searches using the full data set of XMASS-I, a single-phase liquid xenon detector, are reported in this paper. Stable XMASS-I data taking accumulated a total live time of 1590.9 days between November 20, 2013 and February 1, 2019 with an analysis threshold of ${\rm 1.0\,keV_{ee}}$. In the latter half of data taking a lower analysis threshold of ${\rm 0.5\,keV_{ee}}$ was also available through a new low threshold trigger. Searching for a WIMP signal in the detector's 97~kg fiducial volume yielded a limit on the WIMP-nucleon scattering cross section of ${\rm 1.4\times 10^{-44}\, cm^{2}}$ for a ${\rm 60\,GeV/c^{2}}$ WIMP at the 90$\%$ confidence level. We also searched for WIMP induced annual modulation signatures in the detector's whole target volume, containing 832~kg of liquid xenon. For nuclear recoils of a ${\rm 8\,GeV/c^{2}}$ WIMP this analysis yielded a 90\% CL cross section limit of ${\rm 2.3\times 10^{-42}\, cm^{2}}$. At a WIMP mass of ${\rm 0.5\, GeV/c^{2}}$ the Migdal effect and Bremsstrahlung signatures were evaluated and lead to 90\% CL cross section limits of ${\rm 1.4\times 10^{-35}\, cm^{2}}$ and ${\rm 1.1\times 10^{-33}\, cm^{2}}$ respectively.

TOI-1695 is a V-mag=13 M-dwarf star from the northern hemisphere at 45$\,$pc from the Sun, around which a 3.134-day periodic transit signal from a super-Earth candidate was identified in TESS photometry. With a transit depth of 1.3$\,$mmag, the radius of candidate TOI-1695.01 was estimated by the TESS pipeline to be 1.82$\,$R$_\oplus$ with an equilibrium temperature of $\sim$620$\,$K. We successfully detect a reflex motion of the star and establish it is due to a planetary companion at an orbital period consistent with the photometric transit period thanks to a year-long radial-velocity monitoring of TOI-1695 by the SPIRou infrared spectropolarimeter. We use and compare different methods to reduce and analyse those data. We report a 5.5-$\sigma$ detection of the planetary signal, giving a mass of $5.5$$\pm$$1.0\,$M$_\oplus$ and a radius of 2.03$\pm$0.18$\,$R$_\oplus$. We derive a mean equilibrium planet temperature of 590$\pm$90$\,$K. The mean density of this small planet of 3.6$\pm$1.1$\,$g$\,$cm$^{-3}$ is similar (1.7-$\sigma$ lower) than that of the Earth. It leads to a non-negligible fraction of volatiles in its atmosphere with $f_{H,He}$=0.28$^{+0.46}_{-0.23}$% or $f_\text{water}$=23$\pm$12%. TOI-1695$\,$b is a new sub-Neptune planet at the border of the M-dwarf radius valley that can help test formation scenarios for super-Earth/sub-Neptune-like planets.

An exomoon will produce transit timing variations (TTVs) upon the parent planet and their undersampled nature causes half of such TTVs to manifest within a frequency range of 2 to 4 cycles, irrespective of exomoon demographics. Here, we search through published Kepler TTV data for such signals, applying a battery of significance and robustness checks, plus independent light curve analyses for candidate signals. Using the original transit times, we identify 11 (ostensibly) single-planets with a robust, significant and fast ($P_{TTV}<4$ cycles) TTV signal. However, of these, only 5 are recovered in an independent analysis of the original photometry, underscoring the importance of such checks. The surviving signals are subjected to an additional trifecta of statistical tests to ensure signal significance, predictive capability and consistency with an exomoon. KOI-3678.01, previously validated as Kepler-1513b, is the only case that passes every test, exhibiting a highly significant (>20 sigma) TTV signal with a periodicity, amplitude and shape consistent with that caused by an exomoon. Our analysis finds that this planet is $8.2_{-0.5}^{+0.7}$ $R_{\oplus}$ orbiting at $0.53_{-0.03}^{+0.04}$ AU around a late G-type dwarf. After forecasting the planetary mass, we expect it to be capable of maintaining at least a 0.3 $M_{\oplus}$ exomoon for 5 Gyr, and the TTV signal corresponds to a moon mass as low as 0.75 Lunar masses. We thus encourage follow-up observations and dynamical analysis of this unique signal, but caution skepticism until such data can be obtained.

The outburst activity in the MAXI~J1348--630 has sparked a great deal of controversy, whether the source contains a black hole (BH). Here, we present the results of our analysis of the outburst of MAXI J348--630 using Swift/XRT data. We find that energy spectra in all spectral states can be modeled using a combination of Comptonization and Gaussian iron-line components. In addition, we show that the X-ray photon index, Gamma, is correlated with the mass accretion rate, Mdot. We find that Gamma increases monotonically with Mdot from the LHS to the HSS, and then saturated at Gamma~3. This index behavior is similar to that exhibited by a number of other BH candidates. This result represents observational evidence of the presence of a BH in MAXI~J1348--630. We also show that the value of Gamma is correlated with the QPO frequency, \nu_{L}. Based on this correlation, we applied a scaling method to estimate a BH mass of 14.8 +/- 0.9 M_{\odot}, using the BH binary XTE~J1550--564 as a reference source. The recent discovery of a giant dust scattering ring around MAXI~J1348--630 by SRG/eROSITA has refined distance estimates to this X-ray source. With this distance, we were able to estimate the disk inclination i = (65+/- 7)^0 using the scaling technique for the correlation between Gamma and normalization proportional to Mdot. An initial decrease in kT_s occurred simultaneously with an increase in the illumination fraction, f. At the start of the outburst, the Compton cloud (or "corona") is very extended and, thus, the seed photons injected to the corona from the relatively far-away disk region, where kT_s is about 0.2--0.4 keV. While Mdot increases (or luminosity increases), the corona contracts, thus increasing the seed photon temperature, kT_s.

The search for alternative cosmological models is largely motivated by the growing discordance between the predictions of LCDM and the ever improving observations, such as the disparity in the value of H_0 measured at low and high redshifts. One model, in particular, known as the R_h=ct universe, has been highly successful in mitigating or removing all of the inconsistencies. In this picture, however, the anisotropies in the cosmic microwave background (CMB) would have emerged at a redshift z ~ 16, rather than via fluctuations in the recombination zone at z~1080. We demonstrate here that a CMB created in the early Universe, followed by scattering through a Pop III generated dust screen, cannot yet be ruled out by the current data. Indeed, the Planck measurements provide a hint of a ~2-4% frequency dependence in the CMB power spectrum, which would be naturally explained as a variation in the optical depth through the dust, but not a Thomson scattering-dominated recombination environment. Upcoming measurements should be able to easily distinguish between these two scenarios, e.g., via the detection of recombination lines at z~1080, which would completely eliminate the dust reprocessing idea.

A spectroscopic investigation of the lithium resonance doublet in $\xi$ Boo A and $\xi$ Boo B in terms of both abundance and isotopic ratio is presented. We obtained new $R$=130\,000 spectra with a signal-to-noise ratio (S/N) per pixel of up to 3200 using the 11.8m LBT and PEPSI. From fits with synthetic line profiles based on 1D-LTE MARCS model atmospheres and 3D-NLTE corrections, we determine the abundances of both isotopes. For $\xi$ Boo A, we find A(Li) = 2.40$\pm$0.03\,dex and $^6$Li/$^7$Li < 1.5$\pm$1.0\,\%\ in 1D-LTE, which increases to $\approx$2.45 for the 3D-NLTE case. For $\xi$ Boo B we obtain A(Li) = 0.37$\pm$0.09\,dex in 1D-LTE with an unspecified $^6$Li/$^7$Li level. Therefore, no $^6$Li is seen on any of the two stars. We consider a spot model for the Li fit for $\xi$ Boo B and find A(Li) = 0.45$\pm$0.09\,dex. The $^7$Li abundance is 23 times higher for $\xi$ Boo A than the Sun's, but three times lower than the Sun's for $\xi$ Boo B while both fit the trend of single stars in the similar-aged M35 open cluster. Effective temperatures are redetermined from the TiO band head strength. We note that the best-fit global metallicities are --0.13$\pm$0.01\,dex for $\xi$ Boo A but +0.13$\pm$0.02\,dex for $\xi$ Boo B. Lithium abundance for the K5V benchmark star 61\,Cyg\,A was obtained to A(Li)$\approx$0.53\,dex when including a spot model but to $\approx$0.15\,dex without a spot model.

The Tip of the Red Giant Branch (TRGB) is an apparent discontinuity in the color-magnitude diagram (CMD) along the giant branch due to the end of the red giant evolutionary phase and is used to measure distances in the local universe. In practice, the tip is often fuzzy and its localization via edge detection response (EDR) relies on several methods applied on a case-by-case basis. It is hard to evaluate how individual choices affect a distance estimation using only a single host field while also avoiding confirmation bias. To devise a standardized approach, we compare unsupervised, algorithmic analyses of the TRGB in multiple halo fields per galaxy, up to 11 fields for a single host and 50 fields across 10 galaxies, using high signal-to-noise stellar photometry obtained by the GHOSTS survey with the Hubble Space Telescope. We first optimize methods for the lowest field-to-field dispersion including spatial filtering to remove star forming regions, smoothing and weighting of the luminosity function, selection of the RGB by color, and tip selection based on the number of likely RGB stars and the ratio of stars above versus below the tip ($R$). We find $R$, which we call the tip `contrast', to be the most important indicator of the quality of EDR measurements; we find that field-to-field EDR repeatability varies from 0.3 mag to $\leq$ 0.05 mag for $R=4$ to 7, respectively, though less than half the fields reach the higher quality. Further, we find that $R$, which varies with the age/metallicity of the stellar population based on models, correlates with the magnitude of the tip (and after accounting for low internal extinction), i.e., a tip-contrast relation with slope of $-0.023\pm0.0046$ mag/ratio, a $\sim 5\sigma$ result that improves standardization of the TRGB. We discuss the value of consistent TRGB standardization across rungs for robust distance ladder measurements.

When accreting X-ray pulsars (XRPs) undergo bright X-ray outbursts, their luminosity-dependent spectral and timing features can be analysed in detail. The XRP GRO J1750-27 recently underwent one of such episodes, during which it was observed with $NuSTAR$ and monitored with $NICER$. Such a data set is rarely available, as it samples the outburst over more than a month at a luminosity that is always exceeding ${\sim}5\times10^{37}\,$erg/s. This value is larger than the typical critical luminosity value, where a radiative shock is formed above the neutron star's surface. Our data analysis of the joint spectra returns a highly ($N_H\sim(5-8)\times10^{22}\,$cm$^{-2}$) absorbed spectrum showing a K$\alpha$ iron line, a soft blackbody component likely originating from the inner edge of the accretion disk, and confirms the discovery of one of the deepest cyclotron lines, at a centroid energy of ${\sim}44\,$keV corresponding to a magnetic field strength of $4.7\times10^{12}\,$G. This value is independently supported by the best-fit physical model for spectral formation in accreting XRPs which, in agreement with recent findings, favours a distance of $14$ kpc and also reflects a bulk-Comptonization dominated accretion flow. Contrary to theoretical expectations and observational evidence from other similar sources, the pulse profiles as observed by $NICER$ through the outburst raise, peak and decay remain remarkably steady. The $NICER$ spectrum, including the iron K$\alpha$ line best-fit parameters, also remain almost unchanged at all probed outburst stages, similar to the pulsed fraction behaviour. We argue that all these phenomena are linked and interpret them as resulting from a saturation effect of the accretion column's emission, which occurs in the high-luminosity regime.

Cosmic birefringence is the in-vacuo rotation of the linear polarization plane experienced by photons of the Cosmic Microwave Background (CMB) radiation when theoretically well-motivated parity-violating extensions of Maxwell electromagnetism are considered. If the angle, parametrizing such a rotation is dependent on the photon's direction, then this phenomenon is called Anisotropic Cosmic Birefringence (ACB). In this paper, we perform for the first time a tomographic treatment of the ACB, by considering photons emitted both at the recombination and reionization epoch. This allows one to extract additional and complementary information about the physical source of cosmic birefringence with respect to the isotropic case. We focus here on the case of an axion-like field $\chi$, whose coupling with the electromagnetic sector induces such a phenomenon, by using an analytical and numerical approach (which involves a modification of the CLASS code). We find that the anisotropic component of cosmic birefringence exhibits a peculiar behavior: an increase of the axion mass implies an enhancement of the anisotropic amplitude, allowing to probe a wider range of masses with respect to the purely isotropic case. Moreover, we show that at large angular scales, the interplay between the reionization and recombination contributions to ACB is sensitive to the axion mass, so that at sufficiently low multipoles, for sufficiently light masses, the reionization contribution overtakes the recombination one, making the tomographic approach to cosmic birefringence a promising tool for investigating the properties of this axion-like field.

Within the dynamically cold low inclination portion of the Classical Kuiper Belt, there exists a population of weakly bound binary systems with a number of unusual properties; most notable of which is their extremely wide orbital separations; beyond 7% of their Hill radii. The stability and evolution of these Ultra-Wide Trans-Neptunian Binaries (TNBs) have, in the past, been studied extensively under the assumption that the primary evolving mechanisms are interactions between the binary components and impacting Trans-Neptunian Objects (TNOs). Here, we instead study their evolution as driven by the gravitational perturbations of close passing but non-impacting TNOs. By simulating these passages, we show that the aggregate effects of encounters over billions of years have a significant effect on Kuiper Belt binary evolution. Such processes can lead to tight binaries widening significantly over time, approaching and sometimes surpassing the separation of the widest known TNBs. We also find that the eccentricity and inclination distributions of observed Ultra-Wide TNBs can be sampled from such widened binaries. While we are unable to produce enough wide binaries to explain their abundance, the orbital properties of ones we do produce are consistent with known wide binaries.

The linear matter power spectrum $P(k,z)$ connects theory with large scale structure observations in cosmology. Its scale dependence is entirely encoded in the matter transfer function $T(k)$, which can be computed numerically by Boltzmann solvers, and can also be computed semi-analytically by using fitting functions such as the well-known Bardeen-Bond-Kaiser-Szalay (BBKS) and Eisenstein-Hu (EH) formulae. However, both the BBKS and EH formulae have some significant drawbacks. On the one hand, although BBKS is a simple expression, it is only accurate up to $10\%$, which is well above the $1\%$ precision goal of forthcoming surveys. On the other hand, while EH is as accurate as required by upcoming experiments, it is a rather long and complicated expression. Here, we use the Genetic Algorithms (GAs), a particular machine learning technique, to derive simple and accurate fitting formulae for the transfer function $T(k)$. When the effects of massive neutrinos are also considered, our expression slightly improves over the EH formula, while being notably shorter in comparison.

The solar wind (SW) is a supersonic outflow of plasma from the solar corona, with the latitudinal speed and density profiles varying with the solar activity. The SW protons charge exchange with the inflowing interstellar neutral atoms and create energetic neutral atoms (ENAs), which bring information on the physical state of the plasma within the boundary region of the heliosphere. The speed of the ENAs depends on their energies, and consequently observations at different energies provide information on different epochs backwards in time. Therefore, understanding the history of the evolution of the SW is important to understand this information. In this paper, we extend the work by \citet{porowski_etal:22a}, who provided the WawHelioIon 3DSW model of the time evolution of latitudinal profiles of the SW speed and density based on results of analysis of interplanetary scintillations (IPS). Based on results of Principal Component Analysis, we seek for correlation between selected solar proxies and the structure of the SW obtained from IPS and show that it is possible to reproduce the evolution of the SW structure during the past three solar cycles based on the proxies. With this, we extend the history of the evolution of the SW structure back to 1976, i.e., to the epoch when observations of the key proxies -- the inclination of the SW current sheet and the solar polar magnetic fields -- became available. We point out the potential of the use of the proxies for forecasting the structure of the SW into the future.

flowMC is a Python library for accelerated Markov Chain Monte Carlo (MCMC) leveraging deep generative modeling. It is built on top of the machine learning libraries JAX and Flax. At its core, flowMC uses a local sampler and a learnable global sampler in tandem to efficiently sample posterior distributions. While multiple chains of the local sampler generate samples over the region of interest in the target parameter space, the package uses these samples to train a normalizing flow model, then uses it to propose global jumps across the parameter space. The flowMC sampler can handle non-trivial geometry, such as multimodal distributions and distributions with local correlations. The key features of flowMC are summarized in the following list: * Since flowMC is built on top of JAX, it supports gradient-based samplers through automatic differentiation such as MALA and Hamiltonian Monte Carlo (HMC). * flowMC uses state-of-the-art normalizing flow models such as Rational-Quadratic Splines to power its global sampler. These models are very efficient in capturing important features within a relatively short training time. * Use of accelerators such as GPUs and TPUs are natively supported. The code also supports the use of multiple accelerators with SIMD parallelism. * By default, Just-in-time (JIT) compilations are used to further speed up the sampling process. * We provide a simple black box interface for the users who want to use flowMC by its default parameters, yet provide at the same time an extensive guide explaining trade-offs while tuning the sampler parameters. The tight integration of all the above features makes flowMC a highly performant yet simple- to-use package for statistical inference.

Many astrophysical and terrestrial scenarios involving magnetic fields can be approached in axial geometry. Although the smoothed particle hydrodynamics (SPH) technique has been successfully extended to magneto-hydrodynamics (MHD), a well-verified, axisymmetric MHD scheme based on such technique does not exist yet. In this work we fill that gap in the scientific literature and propose and check a novel axisymmetric MHD hydrodynamic code, that can be applied to physical problems which display the adequate geometry. We show that the hydrodynamic code built following these axisymmetric hypothesis is able to produce similar results than standard 3D-SPMHD codes with equivalent resolution but with much lesser computational load.

An equation of pair correlation function has been derived from the first two members of BBGKY hierarchy in a weakly coupled inhomogeneous self gravitating system in quasi thermal equilibrium. This work may be useful to study the thermodynamic properties of the central region of a star cluster which is older than a few or more central relaxation time.

We examine indirect detection of dark matter that annihilates into dark glueballs, which in turn decay into the Standard Model via a range of portals. This arises if the dark matter candidate couples to a confining gauge force without light flavours, representative of many possible complex dark sectors. Such Hidden Valley scenarios are being increasingly considered due to non-detection of minimal models as well as theoretical motivations such as the Twin Higgs solution to the little hierarchy problem. Study of dark glueballs in indirect detection has previously been hampered by the difficulty of modeling their production in dark showers. We use the recent GlueShower code to produce the first constraints on dark matter annihilating via dark glueballs into the Standard Model across photon, antiproton, and positron channels. We also fit the Galactic Centre Excess and use this observation, combined with other astrophysical constraints, to show how multi-channel observations can constrain UV and IR details of the theory, namely the exact decay portal and hadronization behaviour respectively. This provides unique complementary discovery and diagnostic potential to Hidden Valley searches at colliders. It is interesting to note that thermal WIMPs annihilating to $\mathcal{O}(10~\mathrm{GeV})$ dark glueballs and then the Standard Model via the Twin-Higgs-like decay portal can account for the Galactic Centre Excess while respecting other constraints.

The Caral civilization developed on the north-central coast of Peru and had an occupation period between 2870 and 1970 years BC. The first studies carried out in the field of archaeoastronomy showed evidence of possible astronomical orientations in some buildings of its capital city, the Ciudad Sagrada de Caral. However, methodological issues cast doubt on these conclusions. A recent study carried out a more general statistical analysis, which covered a total of 55 architectural structures distributed in ten urban settlements that were part of this civilization, thus managing to identify topographic and astronomical orientation patterns. Based on this evidence, we propose to carry out a new study focused on the capital city, with the objective of analyzing the orientation pattern of the city, placing emphasis on the analysis of the most important religious and administrative structures in order to determine their functionality and their possible links with relevant astronomical objects. The study will include field work to measure the various structures and the subsequent statistical analysis of the data, using declination histograms, density functions and probability tests.

We revisit the constraint results of different dynamical dark energy models including the Chevallier-Polarski-Linder (CPL) model with $w(z)=w_{0}+w_{1}\frac{z}{1+z}$ and the other two models with the logarithm parametrization of $w(z)=w_{0}+w_{1}\left(\frac{\ln (2+z)}{1+z}-\ln 2\right)$ and the oscillating parametrization of $w(z)=w_{0}+w_{1}\left(\frac{\sin(1+z)}{1+z}-\sin(1)\right)$. The advantage over the CPL model is that the latter two parametrizations for dark energy can explore the whole evolution history of the universe properly. Using the current latest mainstream observations including the cosmic microwave background and the baryon acoustic oscillation as well as the type Ia supernovae, we perform the $\chi^2$ statistic analysis to global fit these models, finding that the logarithm parametrization and the oscillating parametrization are slightly preferred against the CPL scenario. We constrain the total neutrino mass in these dynamical dark energy models. We find that, compared with those in the CPL model, much looser constraints on $\sum m_{\nu}$ are obtained in the logarithm model and the oscillating model. Consideration of the possible mass ordering of neutrinos reveals that the most stringent constraint on $\sum m_{\nu}$ appears in the degenerate hierarchy case. In addition, we confirm that the normal hierarchy case is slightly favored over the inverted one.

We provide a comprehensive study of observable spectra from dark matter pair-annihilation or decay into sterile (right-handed) neutrinos. This occurs, for instance, in neutrino portal dark matter models, where a sterile neutrino acts as the portal between dark matter and the Standard Model sector. The subsequent decays of right-handed neutrinos produce detectable Standard Model particles, notably photons, positrons, and neutrinos. We study the phenomenology of models where the right-handed neutrino masses are below the GeV scale, as well as models where they are at, or significantly heavier than, the TeV scale. In both instances, and for different reasons, the standard tools, including Monte Carlo simulations, are both inadequate and inaccurate. We present the complete framework to compute the relevant branching ratios for right-handed neutrino decays and the spectra of secondary photons, positrons, and neutrinos for a broad range of dark matter and right-handed neutrino masses. We discuss the general features of such signals, and compare the spectra to standard signals from dark matter annihilation/decay into bottom quarks. Additionally, we provide open source code1 that can be used to compute such spectra. The code is available at https://github.com/LoganAMorrison/blackthorn.

Discrepancies between observations at early and late cosmic epochs, and the vacuum energy problem associated with the interpretation of cosmological constant, are questioning the $\Lambda$CDM model. Motivated by these conceptual and observational facts, extensions of Einstein's gravity are recently intensively considered in view of curing unsolved issues suffered by General Relativity at ultraviolet and infrared scales. Here, we provide a short overview of some aspects of $f(R)$ gravity, focusing, in particular, on cosmological applications. Specifically, Noether symmetries are adopted as a criterion to select viable models and investigate the corresponding dynamics. We thus find solutions to the cosmological field equations, analyzing the behaviour of selected models from the matter-dominated to the present epoch. Moreover, constraints coming from energy conditions and the so-called swampland criteria are also considered. In particular, we qualitatively discuss the possibility of $f(R)$ gravity to account for fixing cosmic tensions.

First order phase transitions are well-motivated and extensively studied sources of gravitational waves (GWs) from the early Universe. The vacuum energy released during such transitions is assumed to be transferred primarily either to the expanding walls of bubbles of true vacuum, whose collisions source GWs, or to the surrounding plasma, producing sound waves and turbulence, which act as GW sources. In this Letter, we study an alternative possibility that has so far not been considered: the released energy gets transferred primarily to feebly interacting particles that do not admit a fluid description but simply free-stream individually. We develop the formalism to study the production of GWs from such configurations, and demonstrate that such GW signals have qualitatively distinct characteristics compared to conventional sources and are potentially observable with near-future GW detectors.