Papers for Monday, May 13 2024

The forward-modelling of galaxy surveys has recently gathered interest as one of the primary methods to achieve the precision on the estimate of the redshift distributions required by stage IV surveys. One of the key aspects of forward-modelling is the connection between the physical properties of galaxies and their intrinsic spectral energy distributions (SEDs), achieved through stellar population synthesis (SPS) codes, e.g. FSPS. However, SPS requires many detailed assumptions about the galaxy constituents, for which the model choice or parameters are currently uncertain. In this work, we perform a sensitivity study of the impact that the SED modelling choices variations have on the mean and scatter of the tomographic galaxy redshift distributions. We use the Prospector-$\beta$ model and its SPS parameters to build observed magnitudes of a fiducial sample of galaxies. We then build new samples by varying one SED modelling choice at a time. We model the colour-redshift relation of these galaxy samples using the KiDS-VIKING remapped version (McCullough et al. 2023) of the Masters et al. (2015) SOM. We place galaxies in the SOM cells according to the simulated galaxy colours. We then build color-selected tomographic bins and compare each variant's binned redshift distributions against the estimates obtained for the fiducial model. We find that the SED components related to the IMF, AGN, gas physics, and attenuation law substantially bias the mean and the scatter of the tomographic redshift distributions with respect to those estimated with the fiducial model. For the uncertainty of these choices currently present in the literature, and regardless of any stellar mass function based reweighting strategy applied, the bias in the mean and the scatter of the tomographic redshift distributions is larger than the precision requirements set by Stage IV galaxy surveys, e.g. LSST and Euclid.

We study the impact of black hole nuclear activity on both the global and radial star formation rate (SFR) profiles in X-ray-selected active galactic nuclei (AGN) in the field of miniJPAS, the precursor of the much wider J-PAS project. Our sample includes 32 AGN with $z<0.3$ detected via the \textit{XMM-Newton} and \textit{Chandra} surveys. For comparison, we assembled a control sample of 71 star-forming (SF) galaxies with similar magnitudes, sizes, and redshifts. To derive the global properties of both the AGN and the control SF sample, we used \texttt{CIGALE} to fit the spectral energy distributions derived from the 56 narrowband and 4 broadband filters from miniJPAS. We find that AGN tend to reside in more massive galaxies than their SF counterparts. After matching samples based on stellar mass and comparing their SFRs and specific SFRs (sSFRs), no significant differences appear. This suggests that the presence of AGN does not strongly influence overall star formation. However, when we used miniJPAS as an integral field unit (IFU) to dissect galaxies along their position angle, a different picture emerges. We find that AGN tend to be more centrally concentrated in mass with respect to SF galaxies. Moreover, we find a suppression of the sSFR up to 1R$\mathrm{_e}$ and then an enhancement beyond 1R$\mathrm{_e}$, strongly contrasting with the decreasing radial profile of sSFRs in SF galaxies. This could point to an inside-out quenching of AGN host galaxies. These findings suggest that the reason we do not see differences on a global scale is because star formation is suppressed in the central regions and enhanced in the outer regions of AGN host galaxies. While limited in terms of sample size, this work highlights the potential of the upcoming J-PAS as a wide-field low-resolution IFU for thousands of nearby galaxies and AGN.

The exploration of dark sector interactions via gravitational waves (GWs) from binary inspirals has been a subject of recent interest. We study dark forces using extreme mass ratio inspirals (EMRIs), pointing out two issues of interest. Firstly, the innermost stable circular orbit (ISCO) of the EMRI, which sets the characteristic length scale of the system and hence the dark force range to which it exhibits enhanced sensitivity, probes force mediator masses that complement those studied with supermassive black hole (SMBH) or neutron star binaries. The LISA mission (the proposed $\mu$Ares detector) will probe mediators with masses $m_V \sim 10^{-16}~{\rm eV}$ ($m_V \sim 10^{-18}~{\rm eV}$), corresponding to ISCOs of $10^6 M_\odot$ ($10^8 M_\odot$) central SMBHs. Secondly, while the sensitivity to dark couplings is typically limited by the uncertainty in the binary component masses, independent mass measurements of the central SMBH through reverberation mapping campaigns or the motion of dynamical tracers enable one to break this degeneracy. Our results, therefore, highlight the necessity for coordinated studies, loosely referred to as "multimessenger", between future $\mu{\rm Hz}-{\rm mHz}$ GW observatories and ongoing and forthcoming SMBH mass measurement campaigns, including OzDES-RM, SDSS-RM, and SDSS-V Black Hole Mapper.

Type Ia supernovae (SNe Ia) are standarizable candles whose observed light curves can be used to infer their distances, which can in turn be used in cosmological analyses. As the quantity of observed SNe Ia grows with current and upcoming surveys, increasingly scalable analyses are necessary to take full advantage of these new datasets for precise estimation of cosmological parameters. Bayesian inference methods enable fitting SN Ia light curves with robust uncertainty quantification, but traditional posterior sampling using Markov Chain Monte Carlo (MCMC) is computationally expensive. We present an implementation of variational inference (VI) to accelerate the fitting of SN Ia light curves using the BayeSN hierarchical Bayesian model for time-varying SN Ia spectral energy distributions (SEDs). We demonstrate and evaluate its performance on both simulated light curves and data from the Foundation Supernova Survey with two different forms of surrogate posterior -- a multivariate normal and a custom multivariate zero-lower-truncated normal distribution -- and compare them with the Laplace Approximation and full MCMC analysis. To validate of our variational approximation, we calculate the pareto-smoothed importance sampling (PSIS) diagnostic, and perform variational simulation-based calibration (VSBC). The VI approximation achieves similar results to MCMC but with an order-of-magnitude speedup for the inference of the photometric distance moduli. Overall, we show that VI is a promising method for scalable parameter inference that enables analysis of larger datasets for precision cosmology.

Black holes have been found over a wide range of masses, from stellar remnants with masses of 5-150 solar masses (Msun), to those found at the centers of galaxies with $M>10^5$ Msun. However, only a few debated candidate black holes exist between 150 and $10^5$ Msun. Determining the population of these intermediate-mass black holes is an important step towards understanding supermassive black hole formation in the early universe. Several studies have claimed the detection of a central black hole in $\omega$ Centauri, the Milky Way's most massive globular cluster. However, these studies have been questioned due to the possible mass contribution of stellar mass black holes, their sensitivity to the cluster center, and the lack of fast-moving stars above the escape velocity. Here we report observations of seven fast-moving stars in the central 3 arcseconds (0.08 pc) of $\omega$ Centauri. The velocities of the fast-moving stars are significantly higher than the expected central escape velocity of the star cluster, so their presence can only be explained by being bound to a massive black hole. From the velocities alone, we can infer a firm lower limit of the black hole mass of $\sim$8,200 Msun, making this a compelling candidate for an intermediate-mass black hole in the local universe.

The galaxy system around NGC4490 was recently highlighted to display a flattened, kinematically correlated structure reminiscent of satellite galaxy planes around other hosts. Since known satellite planes are in tension with $\Lambda$CDM expectations from cosmological simulations, we quantitatively assess for the first time the tension posed by the NGC4490 system. We measure the on-sky flattening as the major-to-minor axis ratio b/a of the satellite distribution and their line-of-sight kinematic correlation. Analogs are selected in the IllustrisTNG-50 simulation and their flattening and correlation are similarly measured. We confirm the strong kinematic coherence of all 12 observed objects with available line-of-sight velocities (of 14 in total): the northern ones approach and the southern ones recede relative to the host. The spatial distribution of all 14 objects is substantially flattened with b/a=0.38 (0.26 considering only the 12 objects with velocities). Such extreme arrangements are rare in the simulation at 0.21-0.35%. This would drop further if at least one of the two satellite objects without velocities is confirmed to follow the kinematic trend, and would become zero if both are rejected as non-members. We also identify a likely galaxy pair in the observed system, and find a similar pair in the best-matching simulated analog. Our findings establish NGC4490 as another strong example of a satellite plane in the Local Volume. This emphasizes that planes of satellites are a more general issue faced by $\Lambda$CDM also beyond the Local Group. The tension with typical systems in simulations suggests that the observed one requires a specific formation scenario, potentially connected to the larger-scale galaxy alignment in its vicinity. The presence of galaxy pairs in the observed and a simulated system hints at the importance such groupings may have to understand satellite planes.

We present a study illustrating the effects of the passage of a Large Magellanic Cloud (LMC) mass satellite on the distance and velocity distributions of satellites in $\Lambda+$Cold Dark Matter simulations of Milky Way (MW) sized halos. In agreement with previous studies, we find that during such a passage the velocity distribution develops a high-velocity tail due to the reflex motion of the inner part of the halo, which can bias velocity-based virial halo mass estimates. When the velocity distribution of MW satellites is corrected for effects of the LMC passage, it is consistent with the distributions in halos of masses as low as $M_{\rm 200c}=8\times 10^{11}\, M_\odot$ and as high as $1.5\times 10^{12}\,M_\odot$. We present a new halo mass estimator $M_{\rm 200c}=c\sigma^2_{\rm 3D}r_{\rm med}$, where $c$ is a constant calibrated using satellite systems in the simulated MW-sized halos, $\sigma^2_{\rm 3D}$ is the variance of 3D velocities taken with the sign of the radial velocity of each satellite, and $r_{\rm med}$ is the median halocentric distance of the satellites. We show that the estimator has only $s=8\%$ scatter around the median relation of the estimated and true halo masses and deviates by $<2s$ from the median during the pericentric passage of an LMC-like subhalo. This is because $\sigma^2_{\rm 3D}$ and $r_{\rm med}$ deviate in opposite directions during such passages. We apply the estimator to the MW satellite system and estimate the virial mass of the Milky Way of $M_{\rm 200c}=9.96\pm 1.45\times 10^{11}\, M_\odot$, in good agreement with several recent estimates using other methods.

We compare the L-Galaxies semi-analytic model to deep observational data from the UKIDSS Ultra Deep Survey (UDS) across the redshift range $0.5 < z < 3$. We find that the over-abundance of low-mass, passive galaxies at high redshifts in the model can be attributed solely to the properties of `orphan' galaxies, i.e. satellite galaxies where the simulation has lost track of the host dark matter subhalo. We implement a simple model that boosts the star-formation rates in orphan galaxies by matching them to non-orphaned satellite galaxies at a similar evolutionary stage. This straightforward change largely addresses the discrepancy in the low-mass passive fraction across all redshifts. We find that the orphan problem is somewhat alleviated by higher resolution simulations, but the preservation of a larger gas reservoir in orphans is still required to produce a better fit to the observed space density of low-mass passive galaxies. Our findings are also robust to the precise definition of the passive galaxy population. In general, considering the vastly different prescriptions used for orphans in semi-analytic models, we recommend that they are analysed separately from the resolved satellite galaxy population, particularly with JWST observations reigniting interest in the low-mass regime in which they dominate.

Astrometry from {\it Gaia} DR3 has enabled the discovery of a sample of 3000+ binaries containing white dwarfs (WD) and main-sequence (MS) stars in relatively wide orbits, with orbital periods $P_{\rm orb} = (100-1000)$ d. This population was not predicted by binary population synthesis models before {\it Gaia} and -- if the {\it Gaia} orbits are robust -- likely requires very efficient envelope ejection during common envelope evolution (CEE). To assess the reliability of the {\it Gaia} solutions, we measured multi-epoch radial velocities (RVs) of 31 WD+MS binary candidates with $P_{\rm orb} = (40-300)$ d and \texttt{AstroSpectroSB1} orbital solutions. We jointly fit the RVs and astrometry, allowing us to validate the {\it Gaia} solutions and tighten constraints on component masses. We find a high success rate for the {\it Gaia} solutions, with only 2 out of the 31 systems showing significant discrepancies between their {\it Gaia} orbital solutions and our RVs. Joint fitting of RVs and astrometry allows us to directly constrain the secondary-to-primary flux ratio $\mathcal{S}$, and we find $\mathcal{S}\lesssim 0.02$ for most objects, confirming the companions are indeed WDs. We tighten constraints on the binaries' eccentricities, finding a median $e\approx 0.1$. These eccentricities are much lower than those of normal MS+MS binaries at similar periods, but much higher than predicted for binaries formed via stable mass transfer. We present MESA single and binary evolution models to explore how the binaries may have formed. The orbits of most binaries in the sample can be produced through CEE that begins when the WD progenitor is an AGB star, corresponding to initial separations of $2-5$ au. Roughly 50\% of all post-common envelope binaries are predicted to have first interacted on the AGB, ending up in wide orbits like these systems.

We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$\mu$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{M_{\odot}}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{M_{\odot}} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low velocity dispersion ($\bar{\sigma} = 33 \pm 9$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\bar{\sigma} = 11 ^{+8}_{-4}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a Sérsic index, $n$, of $1.3 \pm 0.2$, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang.

The magnetohydrodynamics (MHD) equations plus 'non-ideal' (Ohmic, Hall, ambipolar) resistivities are widely used to model weakly-ionized astrophysical systems. We show that if gradients in the magnetic field become too steep, the implied charge drift speeds become much faster than microphysical signal speeds, invalidating the assumptions used to derive both the resistivities and MHD equations themselves. Generically this situation will excite microscale instabilities that suppress the drift and current. We show this could be relevant at low ionization fractions especially if Hall terms appear significant, external forces induce supersonic motions, or dust grains become a dominant charge carrier. Considering well-established treatments of super-thermal drifts in laboratory, terrestrial, and Solar plasmas as well as conduction and viscosity models, we generalize a simple prescription to rectify these issues, where the resistivities are multiplied by a correction factor that depends only on already-known macroscopic quantities. This is generalized for multi-species and weakly-ionized systems, and leaves the equations unchanged in the drift limits for which they are derived, but restores physical behavior (driving the system back towards slow drift by diffusing away small-scale gradients in the magnetic field) if the limits are violated. This has important consequences: restoring intuitive behaviors such as the system becoming hydrodynamic in the limit of zero ionization; suppressing magnetic structure on scales below a critical length which can comparable to circumstellar disk sizes; limiting the maximum magnetic amplification; and suppressing the effects of the Hall term in particular. This likely implies that the Hall term does not become dynamically important under most conditions of interest in these systems.

Broad Absorption Line Quasars (BALQs) are actively accreting supermassive black holes that have strong outflows characterized by broad absorption lines in their rest-UV spectra. Variability in these absorption lines occurs over months to years depending on the source. WPVS 007, a low-redshift, low-luminosity Narrow-line Seyfert 1 (NLS1) shows strong variability over shorter timescales, providing a unique opportunity to study the driving mechanism behind this variability that may mimic longer scale variability in much more massive quasars. We present the first variability study using {the} spectral synthesis code SimBAL, which provides velocity-resolved changes in physical conditions of the gas using constraints from multiple absorption lines. Overall, we find WPVS 007 to have a highly ionized outflow with a large mass-loss rate and kinetic luminosity. We determine the primary cause of the absorption-line variability in WPVS 007 to be a change in covering fraction of the continuum by the outflow. This study is the first SimBAL analysis where multiple epochs of observation were fit simultaneously, demonstrating the ability of SimBAL to use the time-domain as an additional constraint in spectral models.

Event Horizon Telescope (EHT) images of the horizon-scale emission around the Galactic Center supermassive black hole Sagittarius A* (Sgr A*) favor accretion flow models with a jet component. However, this jet has not been conclusively detected. Using the "best-bet" models of Sgr A* from the EHT collaboration, we assess whether this non-detection is expected for current facilities and explore the prospects of detecting a jet with VLBI at four frequencies: 86, 115, 230, and 345 GHz. We produce synthetic image reconstructions for current and next-generation VLBI arrays at these frequencies that include the effects of interstellar scattering, optical depth, and time variability. We find that no existing VLBI arrays are expected to detect the jet in these best-bet models, consistent with observations to-date. We show that next-generation VLBI arrays at 86 and 115 GHz -- in particular, the EHT after upgrades through the ngEHT program and the ngVLA -- successfully capture the jet in our tests due to improvements in instrument sensitivity and (u,v) coverage at spatial scales critical to jet detection. These results highlight the potential of enhanced VLBI capabilities in the coming decade to reveal the crucial properties of Sgr A* and its interaction with the Galactic Center environment.

Binary black holes (BBHs) born from the evolution of Population III (Pop. III) stars are one of the main high-redshift targets for next-generation ground-based gravitational-wave (GW) detectors. Their predicted initial mass function and lack of metals make them the ideal progenitors of black holes above the upper edge of the pair-instability mass gap, i.e. with a mass higher than $\approx{}134$ (241) M$_\odot$ for stars that become (do not become) chemically homogeneous during their evolution. Here, we investigate the effects of cluster dynamics on the mass function of BBHs born from Pop. III stars, by considering the main uncertainties on Pop. III star mass function, orbital properties of binary systems, star cluster's mass and disruption time. In our dynamical models, at least $\sim$5% and up to 100% BBH mergers in Pop. III star clusters have primary mass $m_1$ above the upper edge of the pair-instability mass gap. In contrast, only $\lesssim {} 3$% isolated BBH mergers have primary mass above the gap, unless their progenitors evolved as chemically homogeneous stars. The lack of systems with primary and/or secondary mass inside the gap defines a zone of avoidance with sharp boundaries in the primary mass - mass ratio plane. Finally, we estimate the merger rate density of BBHs and, in the most optimistic case, we find a maximum of $\mathcal{R}\approx200\,{\rm Gpc^{-3}\,yr^{-1}}$ at $z\sim15$ for BBHs formed via dynamical capture. For comparison, the merger rate density of isolated Pop. III BBHs is $\mathcal{R}\leq{}10\,{\rm Gpc^{-3}\,yr^{-1}}$, for the same model of Pop. III star formation history.

The convergence between astronomy and data sonification represents a significant advancement in the approach and analysis of cosmic information. By surpassing the visual exclusivity in data analysis in astronomy, innovative projects have developed software that goes beyond visual representation, transforming data into auditory and tactile displays. However, it has been evidenced that this novel technique requires specialized training, particularly for audio format data. This work describes the initial development of a platform aimed at providing training for data analysis in astronomy through sonification. The integration of these tools in astronomical education and research opens new horizons, facilitating a more inclusive and multisensory participation in the exploration of space science.

We present the spectral and timing analyses of \textit{AstroSat} observations of the Black Hole X-ray Binary GX 339-4 when the source was close the peak of the outburst in 2024. We find that both the spectral and timing variability of the source is indicative of it in its steep power law (SPL) state during the observations. We used phenomenological and physical models to understand the physics and geometry of accretion during this spectral state of the source. Spectral fits indicate the presence of an accretion disc with a temperature of $kT\sim$0.82 keV and a hot corona with a spectral index of $\sim$2.2 along with a significant contribution from iron line emission from the accretion disc. Strong QPOs were detected at $\sim$4.6 Hz in the Power Density Spectra of the source along with a harmonics feature. Time and phase lag at the QPO frequency are studied and we find a hard lag at the QPO frequency and at the same time a soft lag at the harmonic frequency. We estimate the spin of the black hole and it was found that $a = 0.99 \pm 0.003$. The height of the coronal region is estimated to be about 2.5 $R_{g}$, which is found to be similar to that observed during the previous outbursts of the source. We attempt to discuss the possible physical scenario for the observed spectral and timing features exhibited by the source.

Reverberation mapping is the leading technique used to measure direct black hole masses outside of the local Universe. Additionally, reverberation measurements calibrate secondary mass-scaling relations used to estimate single-epoch virial black hole masses. The Australian Dark Energy Survey (OzDES) conducted one of the first multi-object reverberation mapping surveys, monitoring 735 AGN up to $z\sim4$, over 6 years. The limited temporal coverage of the OzDES data has hindered recovery of individual measurements for some classes of sources, particularly those with shorter reverberation lags or lags that fall within campaign season gaps. To alleviate this limitation, we perform a stacking analysis of the cross-correlation functions of sources with similar intrinsic properties to recover average composite reverberation lags. This analysis leads to the recovery of average lags in each redshift-luminosity bin across our sample. We present the average lags recovered for the H$\beta$, Mg II and C IV samples, as well as multi-line measurements for redshift bins where two lines are accessible. The stacking analysis is consistent with the Radius-Luminosity relations for each line. Our results for the H$\beta$ sample demonstrate that stacking has the potential to improve upon constraints on the $R-L$ relation, which have been derived only from individual source measurements until now.

The scientific study of the Solar System's minor bodies ultimately starts with a search for those bodies. This chapter presents a review of the use of machine learning techniques to find moving objects, both natural and artificial, in astronomical imagery. After a short review of the classical non-machine learning techniques that are historically used, I review the relatively nascent machine learning literature, which can broadly be summarized into three categories: streak detection, detection of moving point sources in image sequences, and detection of moving sources in shift and stack searches. In most cases, convolutional neural networks are utilized, which is the obvious choice given the imagery nature of the inputs. In this chapter I present two example networks: a Residual Network I designed which is in use in various shift and stack searches, and a convolutional neural network that was designed for prediction of source brightnesses and their uncertainties in those same shift-stacks. In discussion of the literature and example networks, I discuss various pitfalls with the use of machine learning techniques, including a discussion on the important issue of overfitting. I discuss various pitfall associated with the use of machine learning techniques, and what I consider best practices to follow in the application of machine learning to a new problem, including methods for the creation of robust training sets, validation, and training to avoid overfitting.

Context: The abundances of the $\alpha$-elements are key for understanding the early chemical enrichment of the Galactic bulge. The elements of interest present lines in different wavelength regions, and some of them show lines only in part of the spectra. In the present work, the CNO trio, the alpha-elements Mg, Si, Ca, and Ti, and odd-Z Na and Al are examined as measured from optical and H-band lines. Aims: The aim of this work is to carry out a detailed comparison of stellar parameters and abundances derived in the optical and near-infrared (H-band). We also inspect the best available lines for a list of bulge stars previously analyzed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) team in the H-band and by our group in the optical. This work is mainly of interest to spectroscopists. Methods: In the present work, we compared the stellar parameters and abundance results derived from APOGEE H-band spectra with optical analyses based on Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope (VLT/UVES) data for eight bulge stars. Results:We point out the most suitable wavelength region for each of the studied elements, and highlight difficulties in the derivation of stellar parameters both in the optical and H-band. The near-infrared will allow observations of a large number of stars in the near future given new instruments soon to be available. The identification of spectral lines in this spectral region and the investigation of their reliability are ongoing efforts worldwide. New instruments will also allow simultaneous observation of H-band and optical.

The time-integrated polarization degree (PD) at prompt optical band of gamma-ray burst (GRB) was predicted to be less than $20\%$, while the time-resolved one can reach as high as $75\%$ in photosphere model. Polarizations in optical band during GRB prompt phase had not been studied under framework of the magnetic reconnection model. Here, a three-segment power laws of the energy spectrum is used to reconstruct the Stokes parameters of the magnetic reconnection model. The multi-wavelength light curves and polarization curves from the optical band to MeV gamma-rays in GRB prompt phase are studied. We found depending mainly on the jet dynamics there is a long lasting high PD phase at all calculated energy bands for the typical parameter sets. The time-resolved PD could be as high as $50\%$, while the time-integrated one is roughly $17\%$) in optical band. It can reach $60\%$ for the time-resolved PD in X-rays and the time-integrated one is around $(30-40)\%$. The polarization angle (PA) evolution is random in both optical and gamma-ray bands for the photosphere model, while it is roughly a constant in the synchrotron models. Therefore, future time-resolved PA observations in the prompt optical or gamma-ray band could distinguish between the photosphere and the synchrotron models.

We study the evolution of primordial magnetic fields until the recombination epoch, which is constrained by the conservation of magnetic helicity density if they are maximally helical and by the Hosking integral if they are non-helical. We combine these constraints with conditions obtained by estimating time scales of energy dissipation processes to describe the evolution of magnetic field strength and magnetic coherence length analytically. The dissipation processes depend on whether magnetic or kinetic energy is dominant, whether the decay dynamics is linear or not, and whether the dominant dissipation term is shear viscosity or drag force. We apply the description to compare constraints on primordial magnetic fields at different epochs in the early universe and argue that magnetogenesis before the electroweak symmetry breaking is not feasible.

Solar active regions (ARs) determine solar polar fields and cause solar cycle variability within the framework of the Babcock-Leighton (BL) dynamo. The contribution of an AR to the polar field is measured by its dipole field, which results from flux emergence and subsequent flux transport over the solar surface. The dipole fields contributed by an AR before and after the flux transport are referred to as the initial and final dipole fields, respectively. For a better understanding and prediction of solar cycles, in this paper, we provide a database including AR's initial and final dipole fields and the corresponding results of their bipolar magnetic region (BMR) approximation from 1996 onwards. We also identify the repeated ARs and provide the optimized transport parameters. Based on our database, we find that although the commonly used BMR approximation performs well for the initial dipole field, it exhibits a significant deviation for the final dipole field. To accurately assess an AR's contribution to the polar field, the final dipole field with its real configuration should be applied. Despite the notable contributions of a few rogue ARs, approximately the top 500 ARs ordered by their final dipole fields are necessary to derive the polar field at the cycle minimum. While flux transport may increase or decrease the dipole field for an individual AR, its collective impact over all ARs in a cycle is a reduction in their total dipole field.

During cluster assembly, a cluster's virialization process leaves behind signatures that can provide information on its dynamical state. However, no clear consensus yet exists on the best way to achieve this. Therefore, we attempt to derive improved recipes for classifying the dynamical state of clusters in observations using cosmological simulations. Cluster halo mass and their subhalos' mass are used to $ 10^{14}M_{\odot} h^{-1}$ and $10^{10}M_{\odot} h^{-1}$ to calculate five independent dynamical state indicators. We experiment with recipes by combining two to four indicators for detecting specific merger stages like recent and ancient mergers. These recipes are made by plotting merging clusters and a control sample of relaxed clusters in multiple indicators parameter space, and then applying a rotation matrix method to derive the best way to separate mergers from the control sample. The success of the recipe is quantified using the success rate and the overlap percentage of the merger and control histograms along the newly rotated $x$-axis. This provides us with recipes using different numbers of combined indicators and for different merger stage. Among the recipes, the stellar mass gap and center offset are the first and second most dominant of the indicators, and using more indicators improves the effectiveness of the recipe. When applied to observations, our results show good agreement with literature values of cluster dynamical state.

Planet-planetesimal interactions cause a planet to migrate, manifesting as a random walk in semi-major axis. In models for Neptune's migration involving a gravitational upheaval, this planetesimal-driven migration is a side-effect of the dynamical friction required to damp Neptune's orbital eccentricitiy. This migration is noisy, potentially causing Trans Neptunian Objects (TNOs) in mean motion resonance to be lost. With Nbody simulations, we validate a previously-derived analytic model for resonance retention and determine unknown coefficients. We identify the impact of random-walk (noisy) migration on resonance retention for resonances up to fourth order lying between 39 au and 75 au. Using a population estimate for the weak 7:3 resonance from the well-characterized Outer Solar System Origins Survey (OSSOS), we rule out two cases: (1) a planetesimal disk distributed between 13.3 and 39.9 au with $\gtrsim$ 30 Earth masses in today's size distribution and $T_{\rm mig} \gtrsim$ 40Myr and (2) a top-heavy size distribution with $\gtrsim$ 2000 Pluto-sized TNOs and $T_{\rm mig} \gtrsim$ 10Myr, where $T_{\rm mig}$ is Neptune's migration timescale. We find that low-eccentricity TNOs in the heavily populated 5:2 resonance are easily lost due to noisy migration. Improved observations of the low-eccentricity region of the 5:2 resonance and of weak mean motion resonances with Rubin Observatory's Legacy Survey of Space and Time (LSST) will provide better population estimates, allowing for comparison with our model's retention fractions and providing strong evidence for or against Neptune's random interactions with planetesimals.

While the appeal of their extraordinary radio luminosity to our curiosity is undiminished, the nature of fast radio bursts (FRBs) has remained unclear. The challenge has been due in part to small sample sizes and limited understanding of telescope selection effects. We here present the first inclusion of the entire set of one-off FRBs from CHIME/FRB Catalog 1 in frbpoppy. Where previous work had to curate this data set, and fit for few model parameters, we have developed full multi-dimensional Markov chain Monte Carlo (MCMC) capabilities for frbpoppy -- the comprehensive, open-science FRB population synthesis code -- that allow us to include all one-off CHIME bursts. Through the combination of these two advances we now find the best description of the real, underlying FRB population, with higher confidence than before. We show that $4\pm3\times10^{3}$ one-off FRBs go off every second between Earth and $z$=1; and we provide a mock catalog based on our best model, for straight-forward inclusion in other studies. We investigate CHIME side-lobe detection fractions, and FRB luminosity characteristics, to show that some bright, local FRBs are still being missed. We find strong evidence that FRB birth rates evolve with the star formation rate of the Universe, even with a hint of a short (0.1$-$1 Gyr) delay time. The preferred contribution of the hosts to the FRB dispersion agrees with a progenitor birth location in the host disk. This population-based evidence solidly aligns with magnetar-like burst sources, and we conclude FRBs are emitted by neutron stars.

The study of the differences detected between the observed and the predicted positions of Uranus taking only the ancient planets into account led to the discovery of planet Neptune in 1846. This event remains one of the best accomplishments ever achieved in the history of Astronomy and Classical Mechanics. In this paper, we study the perturbations in the orbit of Uranus due to Neptune and its effects from a modern numerical point of view of the $N$-body problem. The effects induced by Pluto in the orbit of Neptune, as the historical search for a ninth planet in the Solar System (recently boostered again with the hypothesis of the so-called Planet Nine) back in the days was propelled by some supposed small inconsistencies in the orbit of the ice giants, are also analyzed.

Recent observations have shown that [P/Fe] in the Galactic stars decreases with increasing [Fe/H] for [Fe/H] > -1 whereas it is almost subsolar for [Fe/H]< -2. These [P/Fe] trends with [Fe/H] have not been well reproduced by previous theoretical models incorporating phosphorus (P) enrichment only by core collapse supernoave. We here show, for the first time, that the trends can be naturally explained by our new models incorporating P enrichment by oxygen-neon (ONe) novae which occur at the surface of massive white dwarfs whose masses are larger than 1.25 M_sun with a metallicity-dependence rate. We also show that the observations can be better reproduced by the models by assuming that (i) the total mass of gaseous ejecta per ONe nov (M_ej) is as high as 4 * 10^{-5} M_sun and (ii) the number of such novae per unit mass (N_ONe) is as large as 0.01 at [Fe/H]~-3. The assumed M_ej is consistent with observations, and the high N_ONe is expected from recent theoretical models for ONe nova fractions. We predict that (i) [P/Fe] increases with increasing [Fe/H] for -2 < [Fe/H] < -1 and (ii) [P/Fe] and [Cl/Fe] trends with [Fe/H] are very similar with each other due to very large yields of P and Cl from ONe nova. It is thus worthwhile for future observations to assess the validity of the proposed P enrichment by ONe novae by confirming or ruling out these two predictions.

Large-scale stellar surveys, such as SDSS-V, 4MOST, WEAVE, and PLATO, require accurate atmospheric models and synthetic spectra of stars for accurate analyses of fundamental stellar parameters and chemical abundances. The primary goal of our work is to develop a new approach to solve radiation-hydrodynamics (RHD) and generate model stellar spectra in a self-consistent and highly efficient framework. We build upon the Copenhagen legacy RHD code, the MULTI3D non-local thermodynamic equilibrium (NLTE) code, and the DISPATCH high-performance framework. The new approach allows us to calculate 3D RHD models of stellar atmospheres on timescales of a few thousand CPU hours and to perform subsequent spectrum synthesis in local thermodynamic equilibrium (LTE) or NLTE for the desired physical conditions within the parameter space of FGK-type stars. We compare the 3D RHD solar model with other available models and validate its performance against solar observations, including the centre-to-limb variation of intensities and key solar diagnostic lines of H and Fe. We show that the performance of the new code allows to overcome the main bottleneck in 3D NLTE spectroscopy and enables calculations of multi-dimensional grids of synthetic stellar observables for comparison with modern astronomical observations.

The imaging performance and sensitivity of an X-ray telescope when observing astrophysical sources are primarily governed by the optical design, geometrical uncertainties (figure errors, surface roughness, and mirror alignment inaccuracies), and the reflectivity properties of the X-ray reflecting mirror surface. To thoroughly evaluate the imaging performance of an X-ray telescope with an optical design similar to Wolter-1 optics, which comprises multiple shells with known geometrical uncertainties and mirror reflectivity properties, appropriate computational tools are essential. These tools are used to estimate the angular resolution and effective area for various source energies and locations and, more importantly, to assess the impact of figure errors on the telescope's imaging performance. Additionally, they can also be used to optimize optics geometry by modifying it in reference to the Wolter-1 optics, aiming to minimize the optical aberration associated with the Wolter-1 configuration. In this paper, we introduce DarsakX, a Python-based ray tracing computational tool specifically designed to estimate the imaging performance of a multi-shell X-ray telescope. DarsakX has the capability to simulate the impact of figure errors present in the axial direction of a mirror shell. The geometrical shape of the mirror shells can be defined as a combination of figure error with the base optics, such as Wolter-1 or Conical optics. Additionally, DarsakX allows the exploration of new optical designs involving two reflections similar to Wolter-1 optics but with an improved angular resolution for wide-field telescopes. Developed through an analytical approach, DarsakX ensures computational efficiency, enabling fast processing.

One of the main goals of the NASA's TESS (Transiting Exoplanet Survey Satellite) mission is the discovery of Earth-like planets around nearby M-dwarf stars. Here, we present the discovery and validation of three new short-period Earth-sized planets orbiting nearby M-dwarfs: TOI- 5720b, TOI-6008b and TOI-6086b. We combined TESS data, ground-based multi-color light curves, ground-based optical and near-infrared spectroscopy, and Subaru/IRD RVs data to validate the planetary candidates and constrain the physical parameters of the systems. In addition, we used archival images, high-resolution imaging, and statistical validation techniques to support the planetary validation. TOI-5720b is a planet with a radius of Rp=1.09 Re orbiting a nearby (23 pc) M2.5 host, with an orbital period of P=1.43 days. It has an equilibrium temperature of Teq=708 K and an incident flux of Sp=41.7 Se. TOI-6008b has a period of P=0.86 day, a radius of Rp=1.03 Re, an equilibrium temperature of Teq=707 K and an incident flux of Sp=41.5 Se. The host star (TOI-6008) is a nearby (36 pc) M5 with an effective temperature of Teff=3075 K. Based on the RV measurements collected with Subaru/IRD, we set a 3-sigma upper limit of Mp<4 M_Earth, thus ruling out a star or brown dwarf as the transiting companion. TOI-6086b orbits its nearby (31 pc) M3 host star (Teff=3200 K) every 1.39 days, and has a radius of Rp=1.18 Re, an equilibrium temperature of Teq=634 K and an incident flux of Sp=26.8 Se. Additional high precision radial velocity measurements are needed to derive the planetary masses and bulk densities, and to search for additional planets in the systems. Moreover, short-period earth-sized planets orbiting around nearby M-dwarfs are suitable targets for atmospheric characterization with the James Webb Space Telescope (JWST) through transmission and emission spectroscopy, and phase curve photometry.

Rosetta/OSIRIS took optical measurements of the intensity of scattered light from the coma of 67P/Churyumov-Gerasimenko over a wide range of phase angles. These data have been used to measure the phase angle dependent radiance profile of the dust coma. We want to provide information about the column area densities of the dust coma as seen from Rosetta. This information in combination with the measured OSIRIS phase function can then be used to determine the scattering phase function of the dust particles. We use a simple numerical model to calculate the dust density in the coma. For this we neglect all forces but solar gravitation and radiation pressure. As this cannot describe particles close to the surface of the comet, we assume starting conditions at a sufficient distance. We evaluate the column area density as observed from Rosetta/OSIRIS and compare the results for different spacecraft positions, dust sizes and surface activity distributions. We find the phase angle dependence of the column area density to be largely independent of particle size and spacecraft positions. The determining factor is the activity distribution across the surface, especially the activity on the night side. For models with no night side activity, we find the column area density at high phase angles to be roughly two orders of magnitude larger than at low phase angles. The radiance profile measured from inside a cometary coma results from the combined effects of a phase angle dependent column area density and the scattering phase function. The radiance profile is therefore strongly dependent on the surface activity distribution, and - unless the dust emission is isotropic - any attempt to infer particle properties (as expressed through the scattering phase function) from such data must take into account and de-bias for this spatial variation of the dust column area density.

We present direct N-body simulations of black-hole-only clusters with up to $2 \cdot 10^4$ compact objects, zero natal spin and no primordial binaries as predicted by various primordial black hole (PBH) Dark Matter models. The clusters' evolution is computed using ${\tt NBODY6\!+\!+GPU}$, including the effects of the tidal field of the galaxy, the kicks of black hole mergers and orbit-averaged energy loss by gravitational radiation of binaries. We investigate clusters with four initial mass distributions, three of which attempt to model a generic PBH scenario using a lognormal mass distribution and a fourth one that can be directly linked to a monochromatic PBH scenario when accretion is considered. More specifically, we dive into the clusters' internal dynamics, describing their expansion and evaporation, along with the resultant binary black hole mergers. We also compare several simulations with and without black hole merger kicks and find modelling implications for the probability of hierarchical mergers.

Large-scale cosmological simulations suggest that feedback from active galactic nuclei (AGN) plays a crucial role in galaxy evolution. In this study, we directly test this hypothesis utilising SDSS spectra of a sample of 48 low redshift (z<0.14) type 2 quasars (QSO2s). We characterised the kinematics of the warm ionised gas by performing a non-parametric analysis of the [OIII]$\lambda 5007$ emission line, as well as constrain the properties of the young stellar populations (YSP) of their host galaxies through spectral synthesis modelling. These analyses revealed that 85% of the QSO2s display gas velocity dispersions larger than that of the stellar component of their host galaxies, indicating the presence of AGN-driven outflows. Comparing the gas kinematics to the AGN properties, we found a positive correlation between gas velocity dispersion and 1.4 GHz radio luminosity but not with AGN bolometric luminosity or Eddington ratio, suggesting that, either the radio luminosity is the key factor in driving outflows or that the outflows themselves are shocking the ISM and producing synchrotron emission. We found that 98% of the sample host YSPs to varying degrees, with star formation rates (SFR) $0 \le SFR \le 92 \mbox{ M}_{\odot} \mbox{yr}^{-1}$, averaged over 100 Myr. We compared the gas kinematics and outflow properties to the SFRs to establish possible correlations which may suggest that the presence of the outflowing gas is impacting SF and find that none exists, leading to the conclusion that, on the scales probed by the SDSS fibre, AGN-driven outflows do not impact SF on the timescales probed in this study. However, we found a positive correlation between the light-weighted stellar ages of the QSO2s and the black hole mass, which may indicate that successive AGN episodes lead to the suppression of SF over the course of galaxy evolution.

Quasi-periodic eruptions (QPEs) are intense repeating soft X-ray bursts with recurrence times about a few hours to a few weeks from galactic nuclei. Though the debates on the origin of QPEs have not completely settled down, more and more analyses favor the interpretation that QPEs are the result of collisions between a stellar mass object (a stellar mass black hole or a main sequence star) and an accretion disk around a supermassive black hole (SMBH) in galactic nuclei. If this interpretation is correct, QPEs will be invaluable in probing the orbits of stellar mass objects in the vicinity of SMBHs, and further inferring the formation of extreme mass ratio inspirals (EMRIs), one of the major targets of spaceborne gravitational wave missions. In this work, we extended the EMRI orbital analysis in Paper I arXiv:2401.11190 to all the known QPE sources with more than $6$ flares observed. Among all the analyzed 5 QPE sources, two distinct EMRI populations are identified: 4 EMRIs are of low orbital eccentricity (consistent with 0) which should be born in the wet EMRI formation channel, and 1 mildly eccentric EMRI (with $e= 0.25^{+0.18}_{-0.20}$ at 2-$\sigma$ confidence level) is consistent with the predictions of both the dry loss-cone formation channel and the Hills mechanism.

Y Cam is classified as one of the oscillating Eclipsing Algol (oEA) systems, which feature a $\delta$ Scuti-type pulsating component alongside mass transfer phenomena. oEA systems are invaluable for probing the evolutionary processes and internal structures of binary components offering insights through their binary variations and oscillating. In this study, we conducted a comprehensive investigation of Y Cam utilizing high-quality photometric TESS data, and high-resolution ELODIE spectra. Through our analysis, we examined the radial velocity variation, performed binary modeling, and calculated the effective temperature values of binary components. The fundamental stellar parameters, such as mass and radius, were determined with an accuracy of $\sim$ $2-6$ %. Furthermore, we examined the orbital period variation to assess the amount of mass transfer using the available minima times of the system and three new minima times obtained from TESS light curves. Analyzing the pulsation structure of the system with the TESS data revealed the dominant pulsation period and amplitude of the pulsating component to be 0.066 d and 4.65 mmag, respectively. Notably, we observed frequency modulations with the orbital period's frequency, along with variations in the amplitude of the highest amplitude frequency across different orbital phases. Remarkably, the amplitude reaches its peak at phases 0.5 and 1. These findings indicate a candidate of a tidally tilted pulsator. Consequently, we investigated the evolutionary status of the binary components using MESA binary evolution models, determining the age of the system to be 3.28 $\pm$ 0.09 Gyr.

Various stellar objects experience a velocity kick at some point in their evolution. These include neutron stars and black holes at their birth or binary systems when one of the two components goes supernova. For most of these objects, the magnitude of the kick and its impact on the object dynamics remains a topic of debate. We investigate how kicks alter the velocity distribution of objects born in the Milky Way disc, both immediately after the kick and at later times, and whether these kicks are encoded in the observed population of Galactic neutron stars. We simulate the Galactic trajectories of point masses on circular orbits in the disc after being perturbed by an isotropic kick, with a Maxwellian distribution of magnitudes with $\sigma=265$ km/s. Then, we simulate the motion of these point masses for $200$ Myr. These trajectories are then evaluated, either for the Milky Way population as a whole or for those passing within two kiloparsecs of the Sun, to get the time evolution of the velocities. During the first $20$ Myr, the bulk velocity of kicked objects becomes temporarily aligned to the cylindrical radius, implying an anisotropy in the velocity orientations. Beyond this age, the velocity distribution shifts toward lower values and settles to a median of $\sim200$ km/s. Around the Sun, the distribution also loses its upper tail, primarily due to unbound objects escaping the Galaxy. We compare this to the velocities of Galactic pulsars and find that pulsars show a similar evolution with characteristic age. The shift of the velocity distribution is due to bound objects spending most of their orbits at larger radii after the kick. They are, therefore, decelerated by the Galactic potential. We find the same deceleration to be predicted for nearby objects and the total population and conclude it is also observed in Galactic pulsars.

In stars that lie on the main sequence in the Hertzsprung Russel diagram, like our sun, hydrogen is fused to helium in a number of nuclear reaction chains and series, such as the proton-proton chain and the carbon-nitrogen-oxygen cycles. Precisely determined thermonuclear rates of these reactions lie at the foundation of the standard solar model. This review, the third decadal evaluation of the nuclear physics of hydrogen-burning stars, is motivated by the great advances made in recent years by solar neutrino observatories, putting experimental knowledge of the proton-proton chain neutrino fluxes in the few-percent precision range. The basis of the review is a one-week community meeting held in July 2022 in Berkeley, California, and many subsequent digital meetings and exchanges. Each of the relevant reactions of solar and quiescent stellar hydrogen burning is reviewed here, from both theoretical and experimental perspectives. Recommendations for the state of the art of the astrophysical S-factor and its uncertainty are formulated for each of them. Several other topics of paramount importance for the solar model are reviewed, as well: recent and future neutrino experiments, electron screening, radiative opacities, and current and upcoming experimental facilities. In addition to reaction-specific recommendations, also general recommendations are formed.

The Mid-InfraRed Instrument (MIRI) on board the James Webb Space Telescope (JWST) probes the chemistry and dust mineralogy of the inner regions of protoplanetary disks. The observed spectra are unprecedented in their detail, complicating interpretations which are mainly based on manual continuum subtraction and 0D slab models. We investigate the physical conditions under which the gas emits in protoplanetary disks. Based on MIRI spectra, we apply a full Bayesian analysis that provides the posterior distributions of dust and molecular properties. For doing so, we introduce the Dust Continuum Kit with Line emission from Gas (DuCKLinG), a model describing the molecular line emission and the dust continuum simultaneously without large computational cost. The dust model is based on work by Juhasz et al. (2009, 2010). The molecular emission is based on LTE slab models, but with radial gradients in column densities and temperatures. The model is compared to observations using Bayesian analysis. We benchmark this model to a complex thermo-chemical ProDiMo model and fit the MIRI spectrum of GWLup. We find that the retrieved molecular conditions from DuCKLinG fall within the true values from ProDiMo. The column densities retrieved by Grant et al. (2023) fall within the retrieved ranges in this study for all examined molecules (CO2, H2O, HCN, and C2H2). Similar overlap is found for the temperatures with only the temperature range of HCN not including the previously found value. This discrepancy may be due to the simultaneous fitting of all molecules compared to the step-by-step fitting of the previous study. There is statistically significant evidence for radial temperature and column density gradients for H2O and CO2 compared to the constant temperature and column density assumed in the 0D slab models. Additionally, HCN and C2H2 emit from a small region with near constant conditions.

Studies of surface brightness (SB) fluctuations in the intracluster medium (ICM) present an indirect estimate of turbulent pressure support and associated Mach numbers. While high resolution X-ray spectroscopy offer means to directly constrain line of sight gas motions, including those due to turbulence, such observations are relatively expensive and will be limited to nearby, bright clusters. In this respect, SB fluctuations are the most economical means to constrain turbulent motions at large cluster radii across a range of redshifts and masses. To forecast what current and future X-ray and SZ facilities may achieve in SB fluctuation studies, I review and synthesize matters of accuracy and precision with respect to calculating power spectra of SB fluctuations, from which turbulent properties are derived. Balance concerns of power spectrum accuracy and precision across a range of spatial scales, I propose the use of three annuli with: [0,0.4] $R_{500}$, [0.4,1] $R_{500}$, and [1,1.5] $R_{500}$. Adopting these three regions, I calculate the uncertainty in the hydrostatic mass bias, $\sigma_{b_{\mathcal{M}}}$, can be achieved for various instruments in several scenarios. I find that $\textit{Lynx}$ and AtLAST are competitive in their constraints at $R_{500}$, while AtLAST should perform better when constraining $\sigma_{b_{\mathcal{M}}}$ at $R_{200}$.

Aims. We present further development of the rolling root mean square (rRMS) algorithm. These improvements consist of an increase in computational speed and an estimation of the uncertainty on the recovered diagnostics. This improvement is named the cross root mean square (xRMS) algorithm. Methods. We used the quantile method to recover the statistics of the line profiles in order to study the evolution of the prominence observed by IRIS on 1 October 2019. We then introduced the improvements to rRMS. These improvements greatly increased the computational speed, and this increase in speed allowed us to use a large model grid. Thus, we utilised a grid of 23 940 models to recover the thermodynamic diagnostics. We used the 'good' (but not 'best') fitting models to recover an estimate of the uncertainty on the recovered diagnostics. Results. The maximum line-of-sight (LOS) velocities were found to be 70 km/s. The line widths were mostly 0.4 Å with the asymmetries of most pixels around zero. The central temperature of the prominence was found to range from 10 kK to 20 kK, with uncertainties of approximately +/-5 to +/-15 kK. The central pressure was around 0.2 dyn/cm2, with uncertainties of +/-0.2 to +/-0.3 dyn/cm2. The ionisation degree ranged from 1 to 1000, with uncertainties mostly in the range +/-10 to +/-100. The electron density was mostly 10^10 /cm3, with uncertainties of mostly +/-10^9. Conclusions. The new xRMS algorithm finds an estimation of the errors of the recovered thermodynamic properties. To our knowledge, this is the first attempt at systematically determining the errors from forward modelling. The large range of errors found may hint at the degeneracies present when using a single ion and/or species from forward modelling. In the future, co-aligned observations of more than one ion and/or species should be used to attempt to constrain this problem.

A series of Multichannel Quantum Defect Theory-based computations have been performed, in order to produce the cross sections of rotational transitions (excitations $N_{i}^{+}-2 \rightarrow$ $N_{i}^{+}$, de-excitations $N_{i}^{+}$ $\rightarrow$ $N_{i}^{+}-2$, with $N_{i}^{+}=2$ to $10$) and of their competitive process, the dissociative recombination, induced by collisions of HD$^+$ ions with electrons in the energy range $10^{-5}$ to 0.3 eV. Maxwell anisotropic rate coefficients, obtained from these cross sections in the conditions of the Heidelberg Test Storage Ring (TSR) experiments ($k_{B}T_{t}=2.8$ meV and $k_{B}T_{l}=45$ $\mu$eV), have been reported for those processes in the same electronic energy range. Maxwell isotropic rate coefficients have been as well presented for electronic temperatures up to a few hundreds of Kelvins. Very good overall agreement is found between our results for rotational transitions and the former theoretical computations as well as with experiment. Furthermore, owing to the full rotational computations performed, the accuracy of the resulting dissociative recombination cross sections is considerably improved.

Context. While Class II stars have already accreted most of their mass, the continued inflow of fresh material via Bondi-Hoyle accretion acts as an additional mass reservoir for their circumstellar disks. This may explain the observed accretion rates of pre-main- sequence (PMS) stars, as well as observational inconsistencies in the mass and angular momentum balance of their disks. Aims. Using a new simulation that reproduces the stellar IMF, we want to quantify the role of Bondi-Hoyle accretion in the formation of Class II disks, as well address the prospect of its observational detection with the James Webb Space Telescope (JWST). Methods. We study the mass and angular momentum of the accreting gas using passively-advected tracer particles in the simulation, and we carry out radiative transfer calculations of near-infrared scattering to generate synthetic JWST observations of Bondi-Hoyle trails of PMS stars. Results. Gas accreting on Class II PMS stars approximately 1 Myr after their formation has enough mass and angular momentum to strongly affect the evolution of the preexisting disks. The accreted angular momentum is large enough to also explain the observed size of Class II disks. The orientation of the angular momentum vector can differ significantly from that of the previously accreted gas, which may result in a significant disk warping or misalignment. We also predict that JWST observations of Class II stars will be able to detect Bondi-Hoyle trails with a 90% success rate with only 2 min exposure time, if stars with accretion rates \dot{M} > 5e-10 Msol/yr and luminosity of L > 0.5 Lsol are selected.

We present a new method to distinguish between different states (e.g., high and low, quiescent and flaring) in astronomical sources with count data. The method models the underlying physical process as latent variables following a continuous-space Markov chain that determines the expected Poisson counts in observed light curves in multiple passbands. For the underlying state process, we consider several autoregressive processes, yielding continuous-space hidden Markov models of varying complexity. Under these models, we can infer the state that the object is in at any given time. The state predictions from these models are then dichotomized with the help of a finite-mixture model to produce state classifications. We apply these techniques to X-ray data from the active dMe flare star EV Lac, splitting the data into quiescent and flaring states. We find that a first-order vector autoregressive process efficiently separates flaring from quiescence: flaring occurs over 30-40% of the observation durations, a well-defined persistent quiescent state can be identified, and the flaring state is characterized by higher temperatures and emission measures.

ExoplANETS-A is an EU Horizon-2020 project with the primary objective of establishing new knowledge on exoplanet atmospheres. Intimately related to this topic is the study of the host-stars radiative properties in order to understand the environment in which exoplanets lie. The aim of this work is to exploit archived data from space-based observatories and other public sources to produce uniform sets of stellar data that can establish new insight on the influence of the host star on the planetary atmosphere. We have compiled X-ray and UV luminosities, which affect the formation and the atmospheric properties of the planets, and stellar parameters, which impact the retrieval process of the planetary-atmosphere's properties and its errors. Our sample is formed of all transiting-exoplanet systems observed by HST or Spitzer. It includes 205 exoplanets and their 114 host-stars. We have built a catalogue with information extracted from public, online archives augmented by quantities derived by the Exoplanets-A work. With this catalogue we have implemented an online database which also includes X-ray and OHP spectra and TESS light curves. In addition, we have developed a tool, exoVOSA, which is able to fit the spectral energy distribution of exoplanets. We give an example of using the database to study the effects of the host-star high-energy emission on the exoplanet atmosphere. The sample has a planet radius valley which is located at 1.8 Earth radii, in agreement with previous studies. Multiplanet systems in our sample were used to test the photoevaporation model and we find that out of 14 systems, only one significant case poses a contradiction to it (K2-3). In summary, the exoplanet and stellar resources compiled and generated by ExoplANETS-A form a sound basis for current JWST observations and for future work in the era of Ariel.

We present a new magnetic atmosphere model code for obtaining synthetic spectral fluxes of hydrogen-rich magnetic white dwarfs. To date, observed spectra have been analyzed with models that neglet the magnetic field effects on the atomic populations. In this work, we incorporate to state-of-art theory into the evaluation of numerical densities of atoms, free electrons and ions in local thermodynamical equilibrium under the action of a magnetic field. The energy distribution of atoms is rigorously evaluated for arbitrary field strength. This energy pattern includes from tightly bound states to metastable or truly bound, highly excited states embedded in the continuum, i.e., over the first Landau level. Finite nuclear mass effects and the coupling between the internal atomic structure and the motion of the atom across the magnetic field are also considered. Synthetic fluxes are generated with integrations of numerical solutions of polarized radiative transfer over the visible stellar disk using a spherical $t$-design method. The atmosphere code is tested with observations from the Sloan Digital Sky Survey for a group of known magnetic white dwarfs. Physical stellar parameters are obtained from least-squares fits to the observed energy distribution and compared with results of previous works. We show that the use of zerofield ionization equilibrium in spectral analysis can lead to underestimated effective temperatures for high magnetic white dwarfs.

Galactic cosmic ray (GCR) particles have a significant impact on the particle-induced background of X-ray observatories, and their flux exhibits substantial temporal variability, potentially influencing background levels. In this study, we present one-day binned high-energy reject rates derived from the Chandra-ACIS and XMM-Newton EPIC-pn instruments, serving as proxies for GCR particle flux. We systematically analyze the ACIS and EPIC-pn reject rates and compare them with the AMS proton flux. Our analysis initially reveals robust correlations between the AMS proton flux and the ACIS/EPIC-pn reject rates when binned over 27-day intervals. However, a closer examination reveals substantial fluctuations within each 27-day bin, indicating shorter-term variability. Upon daily binning, we observe finer. temporal structures in the datasets, demonstrating the presence of recurrent variations with periods of $\sim$ 25 days and 23 days in ACIS and EPIC-pn reject rates, respectively, spanning the years 2014 to 2018. Notably, during the 2016--2017 period, we additionally detect periodicities of $\sim$13.5 days and 9 days in the ACIS and EPIC-pn reject rates, respectively. Intriguingly, we observe a time lag of $\sim$ 6 days between the AMS proton flux and the ACIS/EPIC-pn reject rates during the second half of 2016. This time lag is not visible before 2016 and aftern2017. The underlying physical mechanisms responsible for this time lag remain a subject of ongoing investigation.

Planets in multi-planet systems are expected to migrate inward as near-resonant chains, thus allowing them to undergo gravitational planet-planet interactions and possibly maintain a non-zero obliquity. The TRAPPIST-1 system is in such a near-resonant configuration, making it plausible that TRAPPIST-1e has a non-zero obliquity. In this work, we use the ExoCAM GCM to study the possible climates of TRAPPIST-1e at varying obliquities and atmospheric compositions. We vary obliquity from 0$^\circ$ to 90$^\circ$ and the partial pressure of carbon dioxide from 0.0004 bars (modern Earth-like) to 1 bar. We find that models with a higher obliquity are hotter overall and have a smaller day-night temperature contrast than the lower obliquity models, which is consistent with previous studies. Most significantly, the super-rotating high-altitude jet becomes sub-rotating at high obliquity, thus impacting cloud and surface temperature patterns. As the amount of carbon dioxide increases, the climate of TRAPPIST-1e becomes hotter, cloudier, and less variable. From modeled thermal phase curves, we find that the impact of obliquity could potentially have observable consequences due to the effect of cloud coverage on the outgoing longwave radiation.

We introduce a new diagnostic for the null test of dynamical dark energy. This diagnostic is useful, especially when we include anisotropic baryon acoustic oscillation (BAO) data in an analysis, to quantify the deviations from the standard $\Lambda$CDM model. This null test is independent of any late-time cosmological model or parametrization. With this, we study the evidence for dynamical dark energy in light of Dark Energy Spectroscopic Instrument (DESI) 2024 data combined with cosmic microwave background (CMB) observations of the Planck 2018 mission and local $H_0$ measurements. We find low (around 0.6$\sigma$) to moderate (around 1$\sigma$) evidence for dynamical dark energy which is not that significant. Although we get individual deviations at around 2$\sigma$ at the effective redshift 0.51 of the DESI 2024 data, the average deviations combining all redshift points are almost independent of the inclusion or exclusion of the data at effective redshift 0.51 (the differences are only around 0.3$\sigma$). We get almost similar results for other non-DESI BAO data. The evidences are almost independent of the early-time physics and these evidences have very low dependence on $H_0$ values. The evidence of dynamical dark energy, obtained in this analysis, is not consistent with the DESI key paper results with the $w_0w_a$CDM model, but interestingly, these are almost consistent with the DESI's results with $z_pw_p$CDM model.

We measure the gas-phase abundances of the elements He, N, O, Ne, S, Ar, and Fe in the Lyman-continuum emitting region of the Sunburst Arc, a highly magnified galaxy at redshift $z=2.37$. We detect the temperature-sensitive auroral lines [SII]$\lambda\lambda4069,4076$, [OII]$\lambda\lambda7320,7330$, [SIII]$\lambda6312$, [OIII]$\lambda4363$, and [NeIII]$\lambda3343$ in a stacked spectrum of 5 multiple images of the Lyman-continuum emitter (LCE), from which we directly measure the electron temperature in the low, intermediate, and high ionization zones. We also detect the density-sensitive doublets of [OII]$\lambda\lambda3727,3729$, [SII]$\lambda\lambda6717,6731, and [ArIV]$\lambda\lambda4713,4741$, which constrain the density in both the low- and high-ionization gas. With these temperature and density measurements, we measure gas-phase abundances with similar rigor as studies of local galaxies. We measure a gas-phase metallicity for the LCE of $12+\log(\textrm{O}/\textrm{H}) = 7.97 \pm 0.05$, and find an enhanced nitrogen abundance $\log(\textrm{N}/\textrm{O}) = -0.65^{+0.16}_{-0.25}$. This nitrogen abundance is consistent with enrichment from a population of Wolf-Rayet stars, additional signatures of which are reported in a companion paper. Abundances of sulfur, argon, neon, and iron are consistent with local low-metallicity HII regions and low-redshift galaxies. This study represents the most complete chemical abundance analysis of a galaxy at Cosmic Noon to date, which enables direct comparisons between local HII regions and those in the distant universe.

In this work, we present independent determinations of cosmological parameters and new constraints on $f(T)$ cosmologies, employing two new catalogs related to HII galaxy Hubble and CMB distance priors, along with the local standard measurements, SNIa, $H(z)$ measurements, growth rate data (RSD), and BAO baselines. We found that the marginalised best-fit C.L. $H_0$ and $\sigma_8$ parameters within these cosmologies can relax the current cosmological tensions using HIIG data. This produces a larger range of admissible values for the current Hubble constant, and when all baselines are considered, the uncertainty bands for $H_0$ and the matter density parameter reduce significantly.

We discuss a class of theories that predict a fermionic dark matter candidate from gauge anomaly cancellation. As an explicit example, we study the predictions in theories where the global symmetry associated with baryon number is promoted to a local gauge symmetry. In this context the symmetry-breaking scale has to be below the multi-TeV scale in order to be in agreement with the cosmological constraints on the dark matter relic density. The new physical "Cucuyo" Higgs boson in the theory has very interesting properties, decaying mainly into two photons in the low mass region, and mainly into dark matter in the intermediate mass region. We study the most important signatures at the Large Hadron Collider, evaluating the experimental bounds. We discuss the correlation between the dark matter relic density, direct detection and collider constraints. We find that these theories are still viable, and are susceptible to being probed in current, and future high-luminosity, running.

We characterize in a novel manner the physical properties of the low temperature Fermi gas in the degenerate domain as a function of temperature and chemical potential. For the first time we obtain low temperature $T$ results in the domain where several fermions are found within a de Broglie spatial cell. In this regime, the usual high degeneracy Sommerfeld expansion fails. The other known semi-classical Boltzmann domain applies when fewer than one particle is found in the de Broglie cell. The relative errors of the three approximate methods (Boltzmann limit, Sommerfeld expansion, and the new domain of several particles in the de Broglie cell) are quantified. In order to extend the understanding of the Sommerfeld expansion, we developed a novel characterization of the Fermi distribution allowing the separation of the finite and zero temperature phenomena.

Quadratic quasi-normal modes, generated at second order in black hole perturbation theory, are a promising target for testing gravity in the nonlinear regime with next-generation gravitational wave detectors. While their frequencies have long been known, their amplitudes remain poorly studied. We introduce regular variables and compute amplitudes for Schwarzschild black holes with the Leaver algorithm. We find a nonlinear ratio $\mathcal{R}\simeq0.154e^{-0.068i}$ for the most excited $\ell=4$ mode, matching results from Numerical Relativity. We also predict new low-frequency $\ell=2$ quadratic modes.

We consider primordial nucleosynthesis in the presence of hypothetical quasi-stable doubly charged particles. Existence of $X^{--}$ with macroscopic lifetimes will lead to the formation of its bound states with $^4$He and other light elements, significantly facilitating the subsequent formation of lithium nuclei. From observational constraints on maximum allowable amount of lithium, that we update in this work, we derive strong constraints on the abundance and lifetime of $X^{--}$. In a likely cosmological freeze-out scenario with temperatures initially exceeding the mass of $X^{--}$, the BBN constrains the lifetime of these particles to be less than about 100 seconds. For parametrically long lifetimes, lithium abundance data constrain $X^{--}$ abundance to be less than $10^{-9}$ relative to protons, regardless of whether these particles decay or remain stable. Stable particles could saturate the dark matter density only if their mass is comparable to or in excess of $10^{10}$ GeV, and most of $X^{--}$ will be found in bound states with beryllium nuclei, so that chemically they would appear as abnormally heavy helium isotopes.

We investigate the influence of the reheating temperature of the visible sector on the freeze-in dark matter (DM) benchmark model for direct detection experiments, where DM production is mediated by an ultralight dark photon. Here we consider a new regime for this benchmark: we take the initial temperature of the thermal Standard Model (SM) bath to be below the DM mass. The production rate from the SM bath is drastically reduced due to Boltzmann suppression, necessitating a significant increase in the portal coupling between DM and the SM to match the observed relic DM abundance. This enhancement in coupling strength increases the predicted DM-electron scattering cross section, making freeze-in DM more accessible to current direct detection experiments.

Solar filaments are well-known tracers of polarity inversion lines that separate two opposite magnetic polarities on the solar photosphere. Because observations of filaments began long before the systematic observations of solar magnetic fields, historical filament catalogs can facilitate the reconstruction of magnetic polarity maps at times when direct magnetic observations were not yet available. In practice, this reconstruction is often ambiguous and typically performed manually. We propose an automatic approach based on a machine-learning model that generates a variety of magnetic polarity maps consistent with filament observations. To evaluate the model and discuss the results we use the catalog of solar filaments and polarity maps compiled by McIntosh. We realize that the process of manual compilation of polarity maps includes not only information on filaments, but also a large amount of prior information, which is difficult to formalize. In order to compensate for the lack of prior knowledge for the machine-learning model, we provide it with polarity information at several reference points. We demonstrate that this process, which can be considered as the user-guided reconstruction or super-resolution, leads to polarity maps that are reasonably close to hand-drawn ones, and additionally allows for uncertainty estimation.

In this paper I show the equivalence, under appropriate assumptions, of two alternative methods to account for the presence of selection biases (also called selection effects) in population studies: one is to include the selection effects in the likelihood directly; the other follows the procedure of first inferring the observed distribution and then removing selection effects a posteriori. Moreover, I investigate a potential bias allegedly induced by the latter approach: I show that this procedure, if applied under the appropriate assumptions, does not produce the aforementioned bias.

Ultrahigh-energy neutrinos (UHE$\nu$s) can be used as a valuable probe of superheavy dark matter above $\sim 10^9$ GeV, the latter being difficult to probe with collider and direct detection experiments due to the feebly interacting nature. Searching for radio emissions originating from the interaction of UHE$\nu$s with the lunar regolith enables us to explore energies beyond $10^{12}$ GeV, which astrophysical accelerators cannot achieve. Taking into account the interaction of UHE$\nu$s with the cosmic neutrino background and resulting standard neutrino cascades to calculate the neutrino flux on Earth, we investigate sensitivities of such lunar radio observations to very heavy dark matter. We show that the proposed ultra-long wavelength lunar radio telescope, as well as the existing low-frequency array, can provide the most stringent constraints on decaying or annihilating superheavy dark matter with masses at $\gtrsim 10^{12}$ GeV. The limits are complementary to or even stronger than those from other UHE$\nu$ detectors, such as the IceCube-Gen2 radio array and GRAND.

The latest molecular data - potential energy curves and Rydberg$/$valence interactions - characterizing the super-excited electronic states of CO are reviewed, in order to provide inputs for the study of their fragmentation dynamics. Starting from this input, the main paths and mechanisms for CO$^+$ dissociative recombination are analyzed; its cross sections are computed using a method based on Multichannel Quantum Defect Theory. Convoluted cross sections, giving both isotropic and anisotropic Maxwellian rate-coefficients, are compared with merged-beam and storage-ring experimental results. The calculated cross sections underestimate the measured ones by a factor of $2$, but display a very similar resonant shape. These facts confirm the quality of our approach for the dynamics, and call for more accurate and more extensive molecular structure calculations.

We propose new classes of inflation models based on the modular symmetry, where the modulus field $\tau$ serves as the inflaton. We establish a connection between modular inflation and modular stabilization, wherein the modulus field rolls towards a fixed point along the boundary of the fundamental domain. We find the modular symmetry strongly constrain the possible shape of the potential and identify some parameter space where the inflation predictions agree with cosmic microwave background observations. The tensor-to-scalar ratio is predicted to be smaller than $10^{-6}$ in our models, while the running of spectral index is of the the order of $10^{-4}$.

We examine and compare the gravitational lensing, in the strong field limit, for the spherically symmetric string-inspired Euler-Heisenberg black holes, characterized by additional parameters $Q^2$ and $\alpha-\beta$, representing magnetic charge and coupling constant, respectively. Our analysis reveals a reduction in the photon sphere radius $x_{ps}$, critical impact parameter $u_{ps}$ and angular position $\theta_\infty$ with increasing magnitude of $Q^2$ and $\alpha-\beta$. Consequently, the value of these quantities is consistently lower than that of its GR equivalents. Further, the ratio $r_{mag}$ of the flux of the first image to all others decreases with $Q^2$ and $\alpha-\beta$. Unlike Schwarzschild black holes, string-inspired Euler-Heisenberg black holes have a smaller deflection angle $\alpha_D$, which decreases even more as $Q^2$ increases. Moreover, the time delay for Sgr A* and M87* can reach up to $~11.302$ and $~17085.1$ minutes, respectively, at $Q^2=0.1$ and $\alpha-\beta=-1$, deviating from Schwarzschild black holes by $~0.194$ and $~293.6$ minutes which are not very significant. For Sgr A* and M87*, we determine $\theta_\infty$ to range within $(23.81, 26.28)~\mu as$ and $(17.89, 19.78)~\mu as$ respectively, with angular separations $s$ ranging from $(3.33-5.67)~nas$ for Sgr A* and $(2.51-4.26)~nas$ for M87*. EHT bounds on the $\theta_{sh}$ of Sgr A* and M87* within the $1\sigma$ region, bound the parameters $Q^2$ and $\alpha-\beta$ as: for Sgr A* $0.29278 \le Q^2 \le 0.60778$ and for M87* $0 < Q^2 \le 0.08473$, but in both the cases we found no bound on the parameter $\alpha-\beta$. We also estimate the parameters $\alpha-\beta$ and $Q^2$ associated with string-inspired Euler-Heisenberg black holes using the EHT observation results of Sgr A* and M87*.

Primordial neutrino-antineutrino asymmetries can be constrained through big-bang nucleosynthesis (BBN) relic abundances and cosmic microwave background (CMB) anisotropies, both observables being sensitive to neutrino properties. The latter constraint, which is due to gravitational effects from all neutrino flavors, is very minute since it is at least quadratic in the asymmetries. On the contrary, the constraints from primordial abundances presently dominate, although these abundances are almost only sensitive to the electron flavor asymmetry. It is generally assumed that neutrino asymmetries are sufficiently averaged by flavor oscillations prior to BBN, which allows to constrain a common primordial neutrino asymmetry at the epoch of BBN. This simplified approach suffers two caveats that we deal with in this article, combining a neutrino evolution code and BBN calculation throughout the MeV era. First, flavor "equilibration" is not true in general, therefore an accurate dynamical evolution of asymmetries is needed to connect experimental observables to the primordial asymmetries. Second, the approximate averaging of asymmetries through flavor oscillations is associated to a reheating of the primordial plasma. It is therefore crucial to correctly describe the interplay between flavor equilibration and neutrino decoupling, as an energy redistribution prior to decoupling does not significantly alter the final effective number of neutrino species' value. Overall, we find that the space of allowed initial asymmetries is generically unbound when using currently available primordial abundances and CMB measurements. We forecast constraints using future CMB experiment capabilities, which should reverse this experimental misfortune.

Polarization wiggling is a direct gravitational effect exerted by a gravitational field on any electromagnetic radiation traversing it. This effect is investigated in linear gravity for spacetimes with flat or conformally flat backgrounds. First, we show how the polarization wiggle rates can be computed in the conformal frame and transformed to the physical frame. Then it is shown that polarization wiggling is not sensitive to scalar perturbations to the metric, while vector and tensor perturbations do induce polarization wiggling. This poses two natural questions: Can polarized electromagnetic radiation be used to measure vectorial and tensorial components of gravitational fields directly? And if so, how? Next, polarization wiggling is studied for an arbitrary vector perturbation to the spacetime metric. In a stationary spacetime, the polarization wiggle rate is proportional to the difference in frame dragging rates between radiation emission and measurement events. We show how this can be used to measure the angular momentum of a gravitational source if the emitter orbits the gravitational source on a known orbit. Next, the effect of a gravitational tensor mode with arbitrary polarization is studied. Finally, the polarization wiggling effect induced by a gravitational tensor mode is analyzed. Two cases are demonstrated: A spacetime with a flat Minkowski background and an expanding cosmology with a conformally flat background. In both cases, the polarization wiggling frequency equals the frequency of the gravitational tensor mode, while the other state parameters are encoded in the polarization wiggling amplitude and phase of the polarized radiation. We show that measurements of polarization wiggling frequency, amplitude and phase of polarized radiation from different directions enables all state parameters of a gravitational tensor mode to be determined.

In this work, we study interaction between dark energy and dark matter, where dark energy is described by a massive vector field, and dark matter is modelled as a fluid. We present new interaction term, which affects only perturbations and can give interesting phenomenology. Then we present a general Lagrangian for the interacting vector dark energy with dark matter. For the dark energy, we choose Proca theory with $G_{3}$ term to study its phenomenological consequence. For this model, we explore both background and perturbation dynamics. We also present the no-ghost condition for tensor modes and scalar modes. Subsequently, we also study the evolution of the overdensities of both baryon and cold dark matter in the high$-k$ limit. We show that the effective gravitational coupling is modified for cold dark matter and baryon. We also choose a simple concrete model and numerically show a suppression in the growth of cold dark matter overdensity.

We study thermal and non-thermal resonant leptogenesis in a general setting where a heavy scalar $\phi$ decays to right-handed neutrinos (RHNs) whose further out-of-equilibrium decay generates the required lepton asymmetry. Domination of the energy budget of the Universe by the $\phi$ or the RHNs alters the evolution history of the primordial gravitational waves (PGW), of inflationary origin, which re-enter the horizon after inflation, modifying the spectral shape. The decays of $\phi$ and RHNs release entropy into the early Universe while nearly degenerate RHNs facilitate low and intermediate scale leptogenesis. We show that depending on the coupling $y_R$ of $\phi$ to radiation species, RHNs can achieve thermal abundance before decaying, which gives rise to thermal leptogenesis. A characteristic damping of the GW spectrum resulting in two knee-like features or one knee-like feature would provide evidence for low-scale thermal and non-thermal leptogenesis respectively. We explore the parameter space for the lightest right-handed neutrino mass $M_1\in[10^2,10^{14}]$ GeV and washout parameter $K$ that depends on the light-heavy neutrino Yukawa couplings $\lambda$, in the weak ($K < 1$) and strong ($K > 1$) washout regimes. The resulting novel features compatible with observed baryon asymmetry are detectable by future experiments like LISA and ET. By estimating signal-to-noise ratio (SNR) for upcoming GW experiments, we investigate the effect of the scalar mass $M_\phi$ and reheating temperature $T_\phi$, which depends on the $\phi-N$ Yukawa couplings $y_N$.

General-relativistic equilibria of differentially rotating stars are expected in a number of astrophysical scenarios, from core-collapse supernovae to the remnant of binary neutron-star mergers. The latter, in particular, have been the subject of extensive studies where they were modeled with a variety of laws of differential rotation with varying degree of realism. Starting from accurate and fully general-relativistic simulations of binary neutron-star mergers with various equations of state and mass ratios, we establish the time when the merger remnant has reached a quasi-stationary equilibrium and extract in this way realistic profiles of differential rotation. This allows us to explore how well traditional laws reproduce such differential-rotation properties and to derive new laws of differential rotation that better match the numerical data in the low-density Keplerian regions of the remnant. In this way, we have obtained a novel and somewhat surprising result: the dynamical stability line to quasi-radial oscillations computed from the turning-point criterion can have a slope that is not necessarily negative with respect to the central rest-mass density, as previously found with traditional differential-rotation laws. Indeed, for stellar models reproducing well the properties of the merger remnants, the slope is actually positive, thus reflecting remnants with angular momentum at large distances from the rotation axis, and hence with cores having higher central rest-mass densities and slower rotation rates.