Papers for Thursday, Jan 25 2024

Rotational periods derived from autocorrelation (ACF) techniques on stars photometrically similar to the Sun in Kepler data have proven difficult to reliably determine. We investigate various instrumental and astrophysical factors affecting the accuracy of these measurements, including the effects of observational windows and noise, stellar activity and inclination, spectral passbands, and the separate normalization of contiguous segments. We validate that the flux variations due to faculae are very periodic, but starspots are the dominant source of bolometric and visible differential variability in Sun-like stars on rotational timescales. We quantify how much stronger the relative contribution of faculae would have to be to render Sun-like light curves periodic enough to reliably measure with autocorrelation methods. We also quantify how long starspot lifetimes need to be to render pure spot light curves periodic enough. In general, longer observational windows yield more accurate ACF measurements, even when faculae are not present. Due to the enhancement of the relative contribution of faculae, observing stars with intermediate inclinations, during activity minima, and/or through bluer passbands has the effect of strengthening the periodicity of the light curve. We search for other manifestations of faculae in broadband photometry of Sun-like stars and conclude that without absolute flux measurements or restriction to shorter wavelength passbands, differential light curves are uninformative about faculae.

With significantly improved sensitivity, the Einstein Telescope (ET), along with other upcoming gravitational wave detectors, will mark the beginning of precision gravitational wave astronomy. However, the pursuit of surpassing current detector capabilities requires careful consideration of technical constraints inherent in existing designs. The significant improvement of ET lies in the low-frequency range, where it anticipates a one million-fold increase in sensitivity compared to current detectors. Angular control noise is a primary limitation for LIGO detectors in this frequency range, originating from the need to maintain optical alignment. Given the expected improvements in ET's low-frequency range, precise assessment of angular control noise becomes crucial for achieving target sensitivity. To address this, we developed a model of the angular control system of Advanced Virgo, closely matching experimental data and providing a robust foundation for modeling future-generation detectors. Our model, for the first time, enables replication of the measured coupling level between angle and length. Additionally, our findings confirm that Virgo, unlike LIGO, is not constrained by alignment control noise, even under full power operation.

Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $\mu$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$\sigma$ upper limit of 1-6 parts per million, depending on model assumptions).

We present HOMERUN (Highly Optimized Multi-cloud Emission-line Ratios Using photo-ionizatioN), a new approach to modelling emission lines from photoionized gas that can simultaneously reproduce all observed line intensities from a wide range of ionization levels and with high accuracy. Our approach is based on the weighted combination of multiple single-cloud photoionization models and, contrary to previous works, the novelty of our approach consists in using the weights as free parameters of the fit and constraining them with the observed data. One of the main applications of HOMERUN is the accurate determination of gas-phase metallicities and we show that a critical point is to allow for a variation of the N/O and S/O abundance ratios which can significantly improve the quality of the fit and the accuracy of the results. Moreover, our approach provides a major improvement compared to the single-cloud, constant-pressure models commonly used in the literature. By using high-quality literature spectra of H ii regions where 10 to 20 emission lines (including several auroral lines) are detected with high signal-to-noise ratio, we show that all lines are reproduced by the model with an accuracy better than 10%. In particular, the model is able to simultaneously reproduce [O i], [O ii], [O iii], [S ii], and [S iii] emission lines which, to our knowledge, is an unprecedented result. Finally, we show that the gas metallicities estimated with our models for HII regions in the Milky Way are in agreement with the stellar metallicities than the estimates based on the Te-method. Overall, our method provides a new accurate tool to estimate the metallicity and the physical conditions of the ionized gas. It can be applied to many different science cases from HII regions to AGN and wherever there are emission lines from photoionized gas.

The drag force experienced by astronomical objects moving through gaseous media (gas dynamical friction) plays a crucial role in their orbital evolution. Ostriker (1999) derived a formula for gas dynamical friction by linear analysis, and its validity has been confirmed through subsequent numerical simulations. However, the effect of gas accretion onto the objects on the dynamical friction is yet to be understood. In this study, we investigate the Mach number dependence of dynamical friction considering gas accretion through three-dimensional nested-grid simulations. We find that the net frictional force, determined by the sum of the gravitational force exerted by surrounding gas and momentum flux transferred by accreting gas, is independent of the resolution of simulations. Only the gas outside the Bondi-Hoyle-Lyttleton radius contributes to dynamical friction, because the gas inside this radius is eventually absorbed by the central object and returns the momentum obtained through the gravitational interaction with it. In the subsonic case, the front-back asymmetry induced by gas accretion leads to larger dynamical friction than predicted by the linear theory. Conversely, in the slightly supersonic case with the Mach number between 1 and 1.5, the nonlinear effect leads to a modification of the density distribution in a way reducing the dynamical friction compared with the linear theory. At a higher Mach number, the modification becomes insignificant and the dynamical friction can be estimated with the linear theory. We also provide a fitting formula for dynamical friction based on our simulations, which can be used in a variety of applications.

Core-collapse supernovae are candidate sites for the rapid neutron capture process ($r$-process). We explore the effects of enrichment from $r$-process nuclei on the light-curves of hydrogen-rich supernovae (SNe IIP) and assess the detectability of these signatures. We modify the radiation transport code $\texttt{SNEC}$ to include the approximate effects of opacity and radioactive heating from $r$-process elements in the SN ejecta. We present models spanning a range of total $r$-process masses $M_{\rm r}$ and their assumed radial distribution within the ejecta, finding that $M_{\rm r} \gtrsim 10^{-2} M_\odot$ is sufficient to induce appreciable differences in their light-curves as compared to ordinary SNe IIP (without any $r$-process elements). The primary photometric signatures of $r$-process enrichment include a shortening of the plateau phase, coinciding with the hydrogen-recombination photosphere retreating to the $r$-process-enriched layers, and a steeper post-plateau decline associated with a reddening of the SN colors. We compare our $r$-process-enriched models to ordinary IIP models and observational data, showing that yields of $M_{\rm r} \gtrsim 10^{-2} M_\odot$ are potentially detectable across several of the metrics used by transient observers, provided that the $r$-process rich layers are mixed $\gtrsim$ halfway to the ejecta surface. This detectability threshold can roughly be reproduced analytically using a two-zone ("kilonova within a supernova") picture. Assuming that a small fraction of SNe produce a detectable $r$-process yield $M_{\rm r} \gtrsim 10^{-2}M_\odot$, and respecting constraints on the total Galactic production rate, we estimate that $\gtrsim 10^{3}-10^4$ SNe need be observed to find one $r$-enriched event, a feat that may become possible with the Vera Rubin Observatory.

We present the first application of the weighted skew-spectra to analyze non-Gaussian information in galaxy survey data. Using the tree-level galaxy skew-spectra together with the one-loop power spectrum multipoles, we analyze the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering data, and target our search towards the equilateral bispectrum shape of primordial non-Gaussianity. We use the Effective Field Theory model for the galaxy power spectrum and bispectrum, and account for systematic effects, such as the survey geometry. From our likelihood analysis, we find $f_{\rm NL}^{\rm equil}=-34^{+296}_{-334}$ at $68\%$ CL, consistent with previous works, while systematic errors from our treatment of the survey geometry lead to an unreliable estimation of $f_{\rm NL}^{\rm ortho}$. We further constrain the bias and counterterm parameters, while keeping the cosmology fixed to $\textit{Planck }2018$ values. As a check, we also validate our analysis pipeline using the ${\tt Nseries}$ simulation suite.

Studying the nature of spiral arms is essential for understanding the formation of the intricate disc structure of the Milky Way. The European Space Agency's Gaia mission has provided revolutionary observational data that have uncovered detailed kinematical features of stars in the Milky Way. However, so far the nature of spiral arms continues to remain a mystery. Here we present that the stellar kinematics traced by the classical Cepheids around the Perseus and Outer spiral arms in the Milky Way shows strikingly different kinematical properties from each other: the radial and azimuthal velocities of Cepheids with respect to the Galactic centre show positive and negative correlations in the Perseus and Outer arms, respectively. We also found that the dynamic spiral arms commonly seen in an N-body/hydrodynamics simulation of a Milky Way-like galaxy can naturally explain the observed kinematic trends. Furthermore, a comparison with such a simulation suggests that the Perseus arm is being disrupted while the Outer arm is growing. Our findings suggest that two neighbouring spiral arms in distinct evolutionary phases - growing and disrupting phases - coexist in the Milky Way.

Stellar evolution predicts the existence of a mass gap for black hole remnants produced by pair-instability supernova dynamics, whose lower and upper edges are very uncertain. We study the possibility of constraining the location of the upper end of the pair-instability mass gap, which is believed to appear around ${m_\text{min}} \sim130M_\odot$, using gravitational wave observations of compact binary mergers with next-generation ground-based detectors. While high metallicity may not allow for the formation of first-generation black holes on the "far side" beyond the gap, metal-poor environments containing Population III stars could lead to such heavy black hole mergers. We show that, even in the presence of contamination from other merger channels, next-generation detectors will measure the location of the upper end of the mass gap with a relative precision close to $\Delta {m_\text{min}}/{m_\text{min}} \simeq 4\% (N_\text{det}/100 )^{-1/2}$ at 90% C.L., where $N_\text{det} $ is the number of detected mergers with both members beyond the gap. These future observations could reduce current uncertainties in nuclear and astrophysical processes controlling the location of the gap.

Thomson scattering of cosmic microwave background (CMB) photons imprints various properties of the baryons around galaxies on the CMB. One such imprint, called patchy screening, is a direct probe of the gas density profile around galaxies. It usefully complements the information from the kinematic and thermal Sunyaev-Zel'dovich effects and does not require individual redshifts. In this paper, we derive new estimators of patchy screening called the "temperature inversion" (TI) and "signed" estimators, analogous to the gradient inversion estimator of CMB lensing. Pedagogically, we clarify the relation between these estimators and the standard patchy screening quadratic estimator (QE). The new estimators trade optimality for robustness to biases caused by the dominant CMB lensing and foreground contaminants, allowing the use of smaller angular scales. We perform a simulated analysis to realistically forecast the expected precision of patchy screening measurements from four CMB experiments, ACT, SPT, Simons Observatory (SO) and CMB-S4, cross-correlated with three galaxy samples from BOSS, unWISE and the simulated Rubin LSST Data Challenge 2 catalog. Our results give further confidence in the first detection of this effect from the ACT$\times$unWISE data in the companion paper and show patchy screening will be a powerful observable for future surveys like SO, CMB-S4 and LSST. Implementations of the patchy screening QE and the TI and signed estimators are publicly available in our LensQuEst and ThumbStack software packages, available at https://github.com/EmmanuelSchaan/LensQuEst and https://github.com/EmmanuelSchaan/ThumbStack , respectively.

There are more than 5000 confirmed and validated planets beyond the solar system to date, more than half of which were discovered by NASA's Kepler mission. The catalog of Kepler's exoplanet candidates has only been extensively analyzed under the assumption of white noise (i.i.d. Gaussian), which breaks down on timescales longer than a day due to correlated noise (point-to-point correlation) from stellar variability and instrumental effects. Statistical validation of candidate transit events becomes increasingly difficult when they are contaminated by this form of correlated noise, especially in the low-signal-to-noise (S/N) regimes occupied by Earth--Sun and Venus--Sun analogs. To diagnose small long-period, low-S/N putative transit signatures with few (roughly 3--9) observed transit-like events (e.g., Earth--Sun analogs), we model Kepler's photometric data as noise, treated as a Gaussian process, with and without the inclusion of a transit model. Nested sampling algorithms from the Python UltraNest package recover model evidences and maximum a posteriori parameter sets, allowing us to disposition transit signatures as either planet candidates or false alarms within a Bayesian framework.

Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$\lambda\lambda$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies.

A fundamental question in galaxy and black-hole evolution remains how galaxies and their supermassive black holes have evolved together over cosmic time. Specifically, it is still unclear how the position of X-ray active galactic nucleus (AGN) host galaxies with respect to the star-forming main sequence (MS) may change with the X-ray luminosity ($L_\mathrm{X}$) of the AGN or the stellar mass ($M_\star$) of the host galaxy. We use data from XMM-SERVS to probe this issue. XMM-SERVS is covered by the largest medium-depth X-ray survey (with superb supporting multiwavelength data) and thus contains the largest sample to date for study. To ensure consistency, we locally derive the MS from a large reference galaxy sample. In our analysis, we demonstrate that the turnover of the galaxy MS does not allow reliable conclusions to be drawn for high-mass AGNs, and we establish a robust safe regime where the results do not depend upon the choice of MS definition. Under this framework, our results indicate that less-massive AGN host-galaxies ($\log M_\star\sim9.5-10.5$ $M_\odot$) generally possess enhanced SFRs compared to their normal-galaxy counterparts while the more-massive AGN host galaxies ($\log M_\star\sim10.5-11.5$ $M_\odot$) lie on or below the star-forming MS. Further, we propose an empirical model for how the placement of an AGN with respect to the MS (SFR$_{norm}$) evolves as a function of both $M_\star$ and $L_\mathrm{X}$.

Weakly magnetized neutron stars in X-ray binaries show complex phenomenology with several spectral components that can be associated with the accretion disk, boundary and/or spreading layer, a corona, and a wind. Spectroscopic information alone is, however, not enough to disentangle these components. Additional information about the nature of the spectral components and in particular the geometry of the emission region can be provided by X-ray polarimetry. One of the objects of the class, a bright, persistent, and rather peculiar galactic Type I X-ray burster was observed with the Imaging X-ray Polarimetry Explorer (IXPE) and the X-ray Multi-Mirror Mission Newton (XMM-Newton). Using the XMM-Newton data we estimated the current state of the source as well as detected strong absorption lines associated with the accretion disk wind. IXPE data showed the source to be significantly polarized in the 2-8 keV energy band with the overall polarization degree (PD) of 1.4% at a polarization angle (PA) of -2 degrees (errors at 68% confidence level). During the two-day long observation, we detected rotation of the PA by about 70 degrees with the corresponding changes in the PD from 2% to non-detectable and then up to 5%. These variations in polarization properties are not accompanied by visible changes in spectroscopic characteristics. The energy-resolved polarimetric analysis showed a significant change in polarization, from being strongly dependent on energy at the beginning of the observation to being almost constant with energy in the later parts of the observation. As a possible interpretation, we suggest the presence of a constant component of polarization, strong wind scattering, or different polarization of the two main spectral components with individually peculiar behavior. The rotation of the PA suggests a 30-degree misalignment of the neutron star spin from the orbital axis.

Knowledge of the formation of neutron stars (NSs) in supernova (SN) explosions is of fundamental importance in wide areas of contemporary astrophysics: X-ray binaries, magnetars, radio pulsars, and, not least, double NS systems which merge and become gravitational wave sources. A recent study by Richardson et al. reported that the NS in the Be-star/X-ray binary SGR 0755-2933 (CPD -29 2176) descended from an ultra-stripped SN. Using the same observational data as Richardson et al., however, we find that the majority of progenitor solutions for SGR 0755-2933 are of normal Type Ib/c SNe, which allows for up to several solar masses of material to be ejected in the SN event. To correctly probe the SN explosion physics and inferring pre-SN conditions in a binary system, a full kinematic analysis based on post-SN data is always needed.

We present results from our deep Far-ultraviolet (FUV) survey using {\em AstroSat}/UVIT of a filamentary structure at $z$ $\sim$ $0.072$. A total of four filaments comprising 58 galaxies were probed in our study. We detect 18 filament galaxies in our FUV observation. All filament galaxies are further classified based on their photometric color, nuclear activity, and morphology. The filaments contain galaxies with mixed stellar population types and structures. We do not detect galaxies in our UVIT survey up to a distance of 0.4~Mpc from the filament axis, implying a dearth of recent star formation in the inner region of filaments. We witness an increase in FUV specific star formation rate (sSFR) of filament galaxies with increasing distance from the filament spine (D$_{\rm fil}$). On the contrary, no such gradient in FUV sSFR with cluster-centric distance is observed in the case of cluster galaxies. The high stellar mass filament galaxies (M$_\star$ $\sim$ 10$^{11}$ M$_\odot$) were more star-forming than cluster galaxies in a fixed mass range. The FUV morphology of some filament galaxies detected in the filament outskirts (D$_{\rm fil}$ $\sim$ 0.9~Mpc) is comparable to or slightly extended than their optical counterpart. The mass assembly of galaxies probed by estimating $(FUV-r)$ color gradients show that more centrally star-forming galaxies reside closer to the filament axis regardless of stellar mass. Our results prove that the likelihood of merger interaction and gas starvation increases on approaching the filament spine. We report a definitive and in-homogeneous impact of filaments on the galaxies residing inside them.

The paradigm of the inflationary universe provides a possible explanation for several observed cosmological properties. In order for such solutions to be successful, the universe must convert the energy stored in the inflaton potential into standard model particles through a process known as reheating. In this paper, we reconsider the reheating process for the case where the inflaton potential respects an approximate (but spontaneously broken) conformal symmetry during the reheating epoch. After reviewing the Effective Field Theory of Reheating, we present solutions for the nonlinear oscillations of the inflaton field, derive the corresponding Hill's equation for the coupled reheating field, and determine the stability diagram for parametric resonance. For this class of models -- the simplest realization being a scalar field with a quartic term -- the expansion of the universe drives the coupled field toward a more unstable part of parameter space, in contrast to the standard case. We also generalize this class of models to include quadratic breaking terms in the potential during the reheating epoch and address the process of stability in that universality class of models.

The detailed design and operation of the Smithsonian Astrophysical Observatory's EBIT are described for the first time, including recent design upgrades that have led to improved system stability and greater user control, increasing the scope of possible experiments. Measurements of emission from highly charged Ar were taken to determine the spatial distribution of the ion cloud and electron beam. An optical setup consisting of two lenses, a narrow band filter, and a CCD camera was used to image visible light, while an X-ray pinhole and CCD camera were used to image X-rays. Measurements were used to estimate an effective electron density of 1.77 x 10$^{10}$ cm$^{-3}$. Additionally, observations of X-ray emission from background EBIT gases were measured with a Silicon Lithium detector. Measurements indicate the presence of Ba and Si, which are both easily removed by dumping the trap every 2 s or less.

Potential contribution from gamma-ray sources to the Galactic diffuse gamma rays observed above 100 TeV (sub-PeV energy range) by the Tibet AS{\gamma} experiment is an important key to interpreting recent multi-messenger observations. This paper reveals a surprising fact: none of the 23 Tibet AS{\gamma} diffuse gamma-ray events above 398TeV within the Galactic latitudinal range of |b| < 10 deg. come from the 43 sub-PeV gamma-ray sources reported in the 1LHAASO catalog, which proves that these sources are not the origins of the Tibet AS{\gamma} diffuse gamma-ray events. No positional overlap between the Tibet AS{\gamma} diffuse gamma-ray events and the sub-PeV LHAASO sources currently supports the diffusive nature of the Tibet AS{\gamma} diffuse gamma-ray events, although their potential origin in the gamma-ray sources yet unresolved in the sub-PeV energy range cannot be ruled out.

We report the discovery of two directly imaged, giant planet candidates orbiting the metal-rich DAZ white dwarfs WD 1202-232 and WD 2105-82. JWST's Mid-Infrared Instrument (MIRI) data on these two stars show a nearby resolved source at a projected separation of 11.47 and 34.62 au, respectively. Assuming the planets formed at the same time as their host stars, with total ages of 5.3 and 1.6Gyr, the MIRI photometry is consistent with giant planets with masses about 1-7 Jupiter Masses. The probability of both candidates being false positives due to red background sources is approximately 1 in 3000. If confirmed, these would be the first directly imaged planets that are similar in both age and separation to the giant planets in our own solar system, and they would demonstrate that widely separated giant planets like Jupiter survive stellar evolution. Giant planet perturbers are widely used to explain the tidal disruption of asteroids around metal-polluted white dwarfs. Confirmation of these two planet candidates with future MIRI imaging would provide evidence that directly links giant planets to metal pollution in white dwarf stars.

Near-future facilities observing the high-redshift universe ($2<z<5$) will have an opportunity to take advantage of "multi-tracer" cosmology by observing multiple tracers of the matter density field: Lyman alpha emitters (LAE), Lyman break galaxies (LBG), and CMB lensing $\kappa$. In this work we use Fisher forecasts to investigate the effect of multi-tracers on next-generation facilities. In agreement with previous work, we show that multiple tracers improve constraints primarily from degeneracy breaking, instead of the traditional intuition of sample variance cancellation. Then, we forecast that for both BBN and CMB primary priors, the addition of lensing and LAEs onto a LBG-only sample will gain 25\% or more in many parameters, with the largest gains being factor of $\sim10$ improvement for $f_{\rm EDE}$. We include a preliminary approach towards modelling the impact of radiative transfer (RT) on forecasts involving LAEs by introducing a simplified model at linear theory level. Our results, albeit preliminary, show that the while RT influences LAE-only forecasts strongly, its effect on composite multi-tracer forecasts are limited.

The upcoming Chinese Space Station Telescope (CSST) slitless spectroscopic survey poses a challenge of flux calibration, which requires a large number of flux-standard stars. In this work, we design an uncertainty-aware residual attention network, the UaRA-net, to derive the CSST SEDs with a resolution of R = 200 over the wavelength range of 2500-10000 \AA using LAMOST normalized spectra with a resolution of R = 2000 over the wavelength range of 4000-7000 \AA. With the special structure and training strategy, the proposed model can not only provide accurate predictions of SEDs but also their corresponding errors. The precision of the predicted SEDs depends on effective temperature (Teff), wavelength, and the LAMOST spectral signal-to-noise ratios (SNRs), particularly in the GU band. For stars with Teff = 6000 K, the typical SED precisions in the GU band are 4.2%, 2.1%, and 1.5% at SNR values of 20, 40, and 80, respectively. As Teff increases to 8000 K, the precision increases to 1.2%, 0.6%, and 0.5%, respectively. The precision is higher at redder wavelengths. In the GI band, the typical SED precisions for stars with Teff = 6000 K increase to 0.3%, 0.1%, and 0.1% at SNR values of 20, 40, and 80, respectively. We further verify our model using the empirical spectra of the MILES and find good performance. The proposed method will open up new possibilities for optimal utilization of slitless spectra of the CSST and other surveys.

The EDGE-CALIFA survey provides spatially resolved optical integral field unit (IFU) and CO spectroscopy for 125 galaxies selected from the CALIFA Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different images. By joining these pixel tables to lower dimensional tables that provide radial profiles, integrated spectra, or global properties, it is possible to investigate the dependence of local conditions on large-scale properties. The database is freely accessible and has been utilized in several publications. We illustrate the use of this database and highlight the effects of CO upper limits on the inferred slopes of the local scaling relations between stellar mass, star formation rate (SFR), and H$_2$ surface densities. We find that the correlation between H$_2$ and SFR surface density is the tightest among the three relations.

We systematically explore the influence of the prior of the peak absolute magnitude ($M$) of type Ia supernovae (SNe Ia) on the measurement of the Hubble constant ($H_0$) from SNe Ia observations. We consider five different data-motivated $M$ priors, representing varying levels of dispersion, and assume the spatially-flat $\Lambda$CDM cosmological model. Different $M$ priors lead to relative changes in the mean values of $H_0$ from 2% to 6%. Loose priors on $M$ yield $H_0$ estimates consistent with both the Planck 2018 result and the SH0ES result at the 68% confidence level. We also examine the potential impact of peculiar velocity subtraction on the value of $H_0$, and show that it is insignificant for the SNe Ia observations with redshift $z > 0.01$ used in our analyses.

The Sun is an extraordinary workbench, from which several fundamental astronomical parameters can be measured with high precision. Among these parameters, the solar radius $R_{\odot}$ plays an important role in several aspects, such as in evolutionary models. Despite the efforts in obtaining accurate measurements of $R_{\odot}$, the subject is still debated and measurements are puzzling and/or lacking in many frequency ranges. We aimed to determine the mean, equatorial, and polar radii of the Sun ($R_c$, $R_{eq}$, and $R_{pol}$) in the frequency range 18.1 - 26.1 GHz. We employed single-dish observations from the newly-appointed Medicina "Gavril Grueff" Radio Telescope and the Sardinia Radio Telescope (SRT) throughout 5 years, from 2018 to mid-2023, in the framework of the SunDish project for solar monitoring. Two methods to calculate the radius at radio frequencies are considered and compared. To assess the quality of our radius determinations, we also analysed the possible degrading effects of the antenna beam pattern on our solar maps, using two 2D-models. We carried out a correlation analysis with the evolution of the solar cycle through the calculation of Pearson's correlation coefficient $\rho$. We obtained several values for the solar radius - ranging between 959 and 994 arcsec - and $\rho$, with typical errors of a few arcsec. Our $R_{\odot}$ measurements, consistent with values reported in literature, suggest a weak prolatness of the solar limb ($R_{eq}$ > $R_{pol}$), although $R_{eq}$ and $R_{pol}$ are statistically compatible within 3$\sigma$ errors. The correlation analysis using the solar images from Grueff shows (1) a positive correlation between the solar activity and the temporal variation of $R_c$ (and $R_{eq}$) at all observing frequencies, and (2) a weak anti-correlation between the temporal variation of $R_{pol}$ and the solar activity at 25.8 GHz.

We present ${\sim}0.2$ arcsec ($\sim$80 au) resolution observations of the CO (2-1) and SiO (5-4) lines made with the Atacama large millimeter/submillimeter array toward an extremely young intermediate-mass protostellar source (t$_{\rm dyn}<$1000 years), MMS 1 located in the Orion Molecular Cloud-3 region. We have successfully imaged a very compact CO molecular outflow associated with MMS 1, having deprojected lobe sizes of $\sim$18000 au (red-shifted lobe) and $\sim$35000 au (blue-shifted lobe). We have also detected an extremely compact ($\lesssim$1000 au) and collimated SiO protostellar jet within the CO outflow. The maximum deprojected jet speed is measured to be as high as 93 km s$^{-1}$. The SiO jet wiggles and displays a chain of knots. Our detection of the molecular outflow and jet is the first direct evidence that MMS 1 already hosts a protostar. The position-velocity diagram obtained from the SiO emission shows two distinct structures: (i) bow-shocks associated with the tips of the outflow, and (ii) a collimated jet, showing the jet velocities linearly increasing with the distance from the driving source. Comparisons between the observations and numerical simulations quantitatively share similarities such as multiple-mass ejection events within the jet and Hubble-like flow associated with each mass ejection event. Finally, while there is a weak flux decline seen in the 850 $\mu$m light curve obtained with JCMT/SCUBA 2 toward MMS 1, no dramatic flux change events are detected. This suggests that there has not been a clear burst event within the last 8 years.

Although red clump stars function as reliable standard candles, their surface characteristics (i.e. T$_{\text{eff}}$, log~$g$, and [Fe/H]) overlap with those of red giant branch stars, which are not standard candles. Recent results have revealed that spectral features containing carbon (e.g. CN molecular bands) carry information correlating with the ''gold-standard'' asteroseismic classifiers that distinguish red clump from red giant branch stars. However, the underlying astrophysical processes driving the correlation between these spectroscopic and asteroseismic quantities in red giants remain inadequately explored. This study aims to enhance our understanding of this ''spectro-seismic'' effect, by refining the list of key spectral features predicting red giant evolutionary state. In addition, we conduct further investigation into those key spectral features to probe the astrophysical processes driving this connection. We employ the data-driven The Cannon algorithm to analyse high-resolution ($R\sim80,000$) Veloce spectra from the Anglo-Australian Telescope for 301 red giant stars (where asteroseismic classifications from the TESS mission are known for 123 of the stars). The results highlight molecular spectroscopic features, particularly those containing carbon (e.g. CN), as the primary indicators of the evolutionary states of red giant stars. Furthermore, by investigating CN isotopic pairs (that is, $^{12}$C$^{14}$N and $^{13}$C$^{14}$N) we find statistically significant differences in the reduced equivalent widths of such lines, suggesting that physical processes that change the surface abundances and isotopic ratios in red giant stars, such as deep mixing, are the driving forces of the ''spectro-seismic'' connection of red giants.

We search for MeV gamma-ray emission between 0.75-30 MeV from five galaxy clusters, viz. Coma, VIRGO, SPT-CL J2012-5649, Bullet, and El Gordo, using archival data from the COMPTEL telescope. For this purpose we use three search templates: point source, radial disk and radial Gaussian. We do not detect any signals from Coma, SPT-CL J2012-5649, Bullet and El Gordo clusters with the 95\% c.l. photon energy flux limit $\sim 10^{-10} ~\rm{erg /cm^2/s} $. For VIRGO, we detect a non-zero signal between 0.92 to 1.63 MeV having marginal significance of about $2.5\sigma$, with the observed energy flux $ \sim 10^{-9}~\rm{ergs/cm^2/s}$. However, we do not confirm the previously reported evidence in literature for a gamma-ray line from Coma and VIRGO clusters between 5-7 MeV.

A quantitative analytical model based on the semblance method between the modulation factor with solar phenomena is proposed. Different Local Interstellar Spectra (LIS) have been computed to introduce into a transport equation solution which in turn have been introduced into the atmospheric yield function (Caballero-Lopez & Moraal 2012), that allows to compute a Cosmic Rays (CR) Modulation Factor. The results were as expected. &here are correlation between modulation factor and sunspots, and anticorrelation between modulation factor and mean magnetic field. A transport equation's solution is necessary to compute atmospheric yield function, in this case the used transport equation's solutions were convection-diffusion and force field. Both solutions offer similar models, yet the force field solution shows a higher correlation value in the semblance than the convection-diffusion solution. Several LIS were also computed because they are introduced into the transport equation solutions. The used LIS were Lagner, Potgieter and Webber LIS in 2003, Burguer and Potgieter LIS in 2000, Garcia-Munoz, Mason and Simpson LIS in 1975 and Ghelfi, Barao, Derome and Maurin LIS in 2017. Those LIS were used because they have a model for different nuclear composition: Helium and Hydrogen. The LIS with more changes when is introduced into the semblance is Garcia-Munoz, Mason and Simpson in 1975.

Compact groups of dwarf galaxies (CGDs) have been observed at low redshifts ($z<0.1$) and are direct evidence of hierarchical assembly at low masses. To understand the formation of CGDs and the galaxy assembly in the low-mass regime, we search for analogues of compact (radius $\leq 100$ kpc) groups of dwarfs ($7 \leq \log[M_{\ast}/{\rm M}_\odot] \leq 9.5$) in the IllustrisTNG highest-resolution simulation. Our analysis shows that TNG50-1 can successfully produce CGDs at $z=0$ with realistic total and stellar masses. We also find that the CGD number density decreases towards the present, especially at $z \lesssim 0.26$, reaching $n \approx 10^{-3.5}$ $\rm cMpc^{-3}$ at $z = 0$. This prediction can be tested observationally with upcoming surveys targeting the faint end of the galaxy population and is essential to constrain galaxy evolution models in the dwarf regime. The majority of simulated groups at $z \sim 0$ formed recently ($\lesssim 1.5 \ \rm Gyr$), and CGDs identified at $z \leq 0.5$ commonly take more than 1 Gyr to merge completely, giving origin to low- to intermediate-mass ($8 \leq \log[M_{\ast}/{\rm M}_\odot] \leq 10$) normally star-forming galaxies at $z=0$. We find that halos hosting CGDs at $z = 0$ formed later when compared to halos of similar mass, having lower stellar masses and higher total gas fractions. The simulations suggest that CGDs observed at $z \sim 0$ arise from a late hierarchical assembly in the last $\sim 3$ Gyr, producing rapid growth in total mass relative to stellar mass and creating dwarf groups with median halo masses of $\sim 10^{11.3}$ $\rm M_\odot$ and B-band mass-to-light ratios mostly in the range $10 \lesssim M/L \lesssim 100$, in agreement with previous theoretical and observational studies.

We have undertaken a deep investigation of a well defined sample of 136 PNe located in a 10x10 degree central region of the Galactic Bulge observed with the ESO VLT and supplemented by archival HST imagery. These studies have provided precise morphologies, major axes position angles and the most robust sample of consistently derived chemical abundances available to date. Using these data we have statistically confirmed, at 5-sigma, the precise PNe population that provides the PNe alignment of major axes previously suggested in the Galactic Bulge, revealed a partial solution to the sulfur anomaly and uncovered interesting morphological, abundance and kinematic features. We summarise the most significant findings here with detailed results appearing in a series of related publications.

We report the direct detection of Lyman Continuum (LyC) emission from 9 galaxies and 1 Active Galactic Nuclei (AGN) at $z$ $\sim$ 1.1-1.6 in the GOODS-North field using deep observations from the Ultraviolet Imaging Telescope (UVIT) onboard AstroSat. The absolute escape fraction of the sources estimated from the far-ultraviolet (FUV) and H$\alpha$ line luminosities using Monte Carlo (MC) analysis of two Inter-Galactic Medium (IGM) models span a range $\sim$ 10 - 55 $\%$. The restframe UV wavelength of the sources falls in the extreme-ultraviolet (EUV) regime $\sim$ 550-700 \AA, the shortest LyC wavelength range probed so far. This redshift range remains devoid of direct detections of LyC emission due to the instrumental limitations of previously available facilities. With UVIT having a very low detector noise, each of these sources are detected with an individual signal-to-noise ratio (SNR) $>$ 3 while for the stack of six sources, we achieve an SNR $\sim$ 7.4. The LyC emission is seen to be offset from the optical centroids and extended beyond the UVIT PSF of 1.$^{\prime\prime}6$ in most of the sources. This sample fills an important niche between GALEX and Cosmic Origins Spectrograph (COS) at low-$z$, and HST WFC3 at high-$z$ and is crucial in understanding the evolution of LyC leakers.

Stellar coronal sources have been observed in the past not only for their astrophysical interest in the field of binary system evolution and interaction, but also for their invaluable roles as benchmarks for plasma spectral models and as calibration sources for high resolution spectroscopic X-ray instruments. These include the gratings on-board Chandra and XMM-Newton, as well as the new generation of high resolution capable-detectors recently flown on-board XRISM and planned for the future also on-board the Athena and the LEM missions. In our previous paper exploiting Chandra/HETG observations of the prototypical coronal source Capella, it has been shown that the centroid energies of the many X-ray emission lines detected in the spectrum of this object change as a function of time due to the Doppler modulation within the binary. This is an effect that needs to be corrected while performing calibrations of high resolution X-ray instruments. In this paper, we extend our previous work on Capella to other known stellar coronal sources which have been observed with the Chandra/HETG (11 objects in total). We measure in several objects clear trends in the velocity shifts along the orbit of the primary star, meaning that in these sources one of the two star components is largely dominating the high energy emission. In a number of systems the trend in the velocity shift is not obvious. This can be ascribed to the fact that both stellar components contribute significantly to the X-ray emission.

The circumplanetary environments in our Solar System host a stunning array of moon and ring systems. Study of these environs has yielded valuable insights into planetary system formation and evolution, and there is every reason to believe that we will have much to learn from the moons and rings that are likely to exist in exoplanetary systems as well. This has motivated a small but growing number of researchers to investigate questions related to the formation, stability, long-term viability, composition, and observability of these exomoons and exorings. Still, due to a range of significant observational challenges, we remain at a relatively early stage of this work. As a result, we continue to face a number of difficult, unanswered questions, but this also means there are myriad opportunities for fundamental contributions to the field. In this review we will examine a variety of important issues for the astronomical community to consider, with an aim of providing a comprehensive understanding of ongoing efforts to identify and characterize exomoons and exorings, while also increasing interest and engagement. We begin with an overview of what we expect from systems hosting moons and/or rings in terms of their architectures, habitability, and observational signatures. We then highlight the contributions from a variety of works that have been aimed at detecting and characterizing them. We conclude by examining the outlook for finding these objects and discussing a number of ongoing challenges that we will want to overcome in the years ahead.

Convective storms (CS) on Saturn appear as bright clouds, irregular in shape and rapidly changing over time, accompanied by lightning phenomena. They cover a wide range of size scales, from 500 km or less and up to sizes that encircle the planet. In this chapter we focus on medium or synoptic-scale events (typical sizes 2,000 - 7,000 km) and planetary-scale major storms known as Great White Spots (GWS). In general, the storms emerge in westward jets and in the peak of the strong eastward equatorial jet. They have been observed preferentially in the northern hemisphere and seem to occur in the same latitude in intervals of 30 and 60 years coinciding with the times of maximum insolation within the seasonal cycle. Moist convection models locate the origin of the storms in the water clouds (at about 10 bars of pressure) and to some extent in ammonia clouds, and their thermodynamics explain the possible periodicity of the phenomena. Once triggered, the storms grow and expand, and their interaction with the zonal wind system gives rise to anticyclonic ovals and to disturbances that propagate zonally covering an entire latitude band in the case of the GWS.

Kinetic exospheric models revealed that the solar wind is accelerated by an ambipolar electric field up to supersonic velocities. The presence of suprathermal Strahl electrons at the exobase can further increase the velocity to higher values, leading to profiles comparable to the observations in the fast and slow wind at all radial distances. Such suprathermal electrons are observed at large distances and recently at low distances as well. Those suprathermal electrons were introduced into the kinetic exospheric model using Kappa distributions. Here, the importance of the exobase's altitude is also underlined for its ability to maintain the electric potential to a higher level for slower winds, conversely to what is induced through the effect of a lower kappa index only. In fact, the exobase is located at lower altitude in the coronal holes where the density is smaller than in the other regions of the corona, allowing the wind originating from the holes to be accelerated from lower distances to higher velocities. The new observations of Parker Solar Probe (PSP) and Solar Orbiter (SolO) from launch to mid-2023 are here used to determine the characteristics of the plasma in the corona so that the model fits best to the averaged observed profiles for the slow and fast winds. The observations at low radial distances show suprathermal electrons already well present in the Strahl in the antisunward direction and a deficit in the sunward direction, confirming the exospheric feature of almost no incoming particles.

Based on 2-minute cadence TESS data, we investigate pulsations of TIC 65138566 and TIC 139729335 and discover them to be two new HADS stars with equally spaced g modes. We recognize the radial fundamental mode $f_1$ = 18.3334 c/d and the first overtone $f_3$ = 23.6429 c/d for TIC 65138566, and identify the highest peak $f_1$ = 19.0955 c/d as the radial fundamental mode for TIC 139729335. For g modes, both stars display a regular period spacing of 2413 s. Through detailed seismological analysis, we deduce that these period spacing patterns correspond to modes with $\ell$ = 1. Moreover, our analysis reveals that with the increase in masses and metallicities, the star should display a higher degree of evolution to match a specific period spacing $\Pi_0$. Conversely, the star should have a lower extent of evolution to match the radial fundamental mode. These two contradictory behaviors allow us to precisely obtain stellar physical parameters. TIC 65138566 and TIC 139729335 are determined to be two main sequence stars that have almost the same range of masses and metallicities, with $M$ = 1.36 $\pm$ 0.06 $M_{\odot}$ and $Z$ = 0.005 $\pm$ 0.002. The hydrogen abundance in the core of TIC 65138566 is estimated to be about 0.28, while TIC 139729335 has a slightly higher value of around 0.31. Finally, we suggest that the HADS $\delta$ Scuti-$\gamma$ Doradus star TIC 308396022 is a main sequence star with $M$ = 1.54 $\pm$ 0.08 $M_{\odot}$, $Z$ = 0.007 $\pm$ 0.001, and $X_{\rm c}$ = 0.18 $\pm$ 0.02.

It is well-established that the magnetic Penrose process (MPP) could be highly efficient (efficiency can even exceed $100\%$) for extracting the energy from a Kerr black hole, if it is immersed in a mG order magnetic field. Considering the exact solution of the magnetized Kerr spacetime, here we derive the exact expression of efficiency ($\eta_{\rm MPP}$) for MPP, which is valid for both the Kerr black hole (BH) as well as Kerr superspinar (SS), and also from the weak magnetic field to an ultra-strong magnetic field $(B)$ which can even distort the original Kerr geometry. We show that although the value of $\eta_{\rm MPP}$ increases upto a certain value of ultra-strong magnetic field ($B_p$), it decreases to zero for $B > B_p$, in case of the Kerr BHs. One intriguing feature that emerges is, $\eta_{\rm MPP}$ acquires the maximum value for the Kerr parameter $a_* \approx 0.786$ (unlike $a_*=1$ for the ordinary PP), decreases for the range $0.786 < a_* \leq 1$, and reaches to $20.7\%$ for $a_*=1$ with a few limitations. This indicates that the BH starts to expel the effect of magnetic field for $a_* > 0.786$, and is fully expelled from the extremal Kerr BH due to the gravitational Meissner effect. As a special case of MPP, we also study the ordinary Penrose process (PP) for the magnetized Kerr spacetime. We show that the efficiency of PP decreases with increasing the magnetic field for the Kerr BH. In fact, the MPP for Kerr BHs, Kerr SSs and the ordinary PP for Kerr SSs can be superefficient for the astrophysical applications to powering engines in the high-energy sources like active galactic nuclei and quasars, in the weak magnetic fields. It is almost impossible to extract the energy from a BH (SS) through MPP (PP) in the ultra-strong magnetic fields.

In Very Long Baseline Interferometry (VLBI), signals from multiple antennas combine to create a sparsely sampled virtual aperture, its effective diameter determined by the largest antenna separation. The inherent sparsity makes VLBI imaging an ill-posed inverse problem, prompting the use of algorithms like the Multiobjective Evolutionary Algorithm by Decomposition (MOEA/D), as proposed in the first paper of this series. This study focuses on extending MOEA/D to polarimetric and time dynamic reconstructions, particularly relevant for the VLBI community and the Event Horizon Telescope Collaboration (EHTC). MOEA/D's success in providing a unique, fast, and largely unsupervised representation of image structure serves as the basis for exploring these extensions. The extension involves incorporating penalty terms specific to total intensity imaging, time-variable, and polarimetric variants within MOEA/D's multiobjective, evolutionary framework. The Pareto front, representing non-dominated solutions, is computed, revealing clusters of proximities. Testing MOEA/D with synthetic datasets representative of EHTC's main targets demonstrates successful recovery of polarimetric and time-dynamic signatures despite sparsity and realistic data corruptions. MOEA/D's extension proves effective in the anticipated EHTC setting, offering an alternative and independent claim to existing methods. It not only explores the problem globally but also eliminates the need for parameter surveys, distinguishing it from Regularized Maximum Likelihood (RML) methods. MOEA/D emerges as a novel and useful tool for robustly characterizing polarimetric and dynamic signatures in VLBI datasets with minimal user-based choices. Future work aims to address the last remaining limitation of MOEA/D, specifically regarding the number of pixels and numerical performance, to establish it within the VLBI data reduction pipeline.

The following results of dissertation are submitted for defense: 1. Precise measurements of coordinates and proper motion of the pulsar PSR 0329+54 using the VLBI method. 2. Establishing of the reason for the discrepancy between the coordinates of pulsars measured by VLBI and timing methods, which comes down to the influence of low-frequency noise at the time of pulse arrivals (TOAs) from the pulsar. 3. A special method for processing timing observations that allows you to correct TOA - coordinates of pulsars. 4. Theoretical dependencies of the behavior of dispersion of pulsar parameter depending on the observation interval and the type of correlated noise, on the basis of which it became possible to propose a new time scale BPT based on the orbital motion of pulsar in binary system stable over long periods of time (more than 10 years). 5. A theory that explains spontaneous changes in the rotational frequency of pulsars through their interaction with the passing gravitating mass and the theoretical power spectrum of low-frequency fluctuations of the pulsar rotational phase caused by gravitational disturbances from the passage of bodies near the pulsar.

GK Dra is a detached eclipsing binary system containing two early-F stars, one evolved, in an orbit with a period of 9.974 d and a small eccentricity. Its eclipsing nature was discovered using Hipparcos data, and pulsations were found in follow-up ground-based data. Extensive observations have been obtained using the Transiting Exoplanet Survey Satellite (TESS), and we use these and published spectroscopy to perform a detailed reanalysis of the system. We determine masses of $1.421 \pm 0.012$ and $1.775 \pm 0.028$ Msun, and radii of $1.634 \pm 0.011$ and $2.859 \pm 0.028$ Rsun. The secondary component is more massive, larger, and slightly cooler than its companion; the eclipses are total. The properties of the system can be matched by theoretical predictions for an age of 1.4 Gyr and a slightly sub-solar metallicity. We measure 15 significant pulsation frequencies in the TESS light curve, of which three are in the frequency domain of $\gamma$ Doradus pulsations and the remaining 12 are $\delta$ Scuti pulsations; the system is thus a hybrid pulsator. The strongest pulsation can be definitively assigned to the secondary star as it has been detected in radial velocities of this object. TESS will observe GK Dra again for ten consecutive sectors in the near future.

CW Eri is a detached eclipsing binary system of two F-type stars with an orbital period of 2.728 d. Light curves from two sectors of observations with the Transiting Exoplanet Survey Satellite (TESS) and previously published radial velocity data are analysed to determine the system's physical properties to high precision. We find the masses of the two stars to be $1.568 \pm 0.016$ Msun and $1.314 \pm 0.010$ Msun, the radii to be $2.105 \pm 0.007$ Rsun and $1.481 \pm 0.005$ Rsun, and the system's orbit to have an eccentricity of $0.0131 \pm 0.0007$. The quality of the TESS photometry allows the definition of a new high-precision orbital ephemeris, however no evidence of pulsation is found. We derive a distance to the system of $191.7\pm 3.8$ pc, a value consistent with the Gaia DR3 parallax which yields a distance of $187.9^{+0.6}_{-0.9}$ pc. The measured parameters of both stellar components are found to be in agreement with theoretical predictions for a solar chemical composition and an age of 1.7 Gyr.

A new generative technique is presented in this paper that uses Deep Learning to reconstruct stellar spectra based on a set of stellar parameters. Two different Neural Networks were trained allowing the generation of new spectra. First, an autoencoder is trained on a set of BAFGK synthetic data calculated using ATLAS9 model atmospheres and SYNSPEC radiative transfer code. These spectra are calculated in the wavelength range of Gaia RVS between 8 400 and 8 800 {\AA}. Second, we trained a Fully Dense Neural Network to relate the stellar parameters to the Latent Space of the autoencoder. Finally, we linked the Fully Dense Neural Network to the decoder part of the autoencoder and we built a model that uses as input any combination of $T_{eff}$, $\log g$, $v_e \sin i$, [M/H], and $\xi_t$ and output a normalized spectrum. The generated spectra are shown to represent all the line profiles and flux values as the ones calculated using the classical radiative transfer code. The accuracy of our technique is tested using a stellar parameter determination procedure and the results show that the generated spectra have the same characteristics as the synthetic ones.

We present a signal-foreground separation algorithm for filtering observational data to extract spectral distortions of the cosmic microwave background (CMB). Our linear method, called the least response method (LRM), is based on the idea of simultaneously minimizing the response to all possible foregrounds with poorly defined spectral shapes and random noise while maintaining a constant response to the signal of interest. This idea was introduced in detail in our previous paper. Here, we have expanded our analysis by taking into consideration all the main foregrounds. We draw a detailed comparison between our approach and the moment internal linear combination method, which is a modification of the internal linear combination technique previously used for CMB anisotropy maps. We demonstrate advantages of LRM and evaluate the prospects for measuring various types of spectral distortions. Besides, we show that LRM suggests the possibility of its improvements if we use an iterative approach with sequential separation and partial subtraction of foreground components from the observed signal. In addition, we estimate the optimal temperature that the telescope's optical system should have in order to detect the chemical type $\mu$ distortions. We present a design of an instrument where, according to our estimates, the optimal contrast between its thermal emission and the CMB allows us to measure such distortions.

We present candl, an automatically differentiable python likelihood for analysing Cosmic Microwave Background (CMB) power spectrum measurements. candl is powered by JAX, which makes it fast and easy to calculate derivatives of the likelihood. This facilitates, for example, robust Fisher matrices without finite-difference methods. We show the benefits of candl through a series of example calculations, covering forecasting, robustness tests, and gradient-based Markov chain Monte Carlo sampling. These also include optimising the band power bin width to minimise parameter errors of a realistic mock data set. Moreover, we calculate the correlation of parameter constraints from correlated and partially overlapping subsets of the SPT-3G 2018 $TT/TE/EE$ data release. In a traditional analysis framework, these tasks are slow and require careful fine-tuning to obtain stable results. As such, a fully differentiable pipeline allows for a higher level of scrutiny; we argue that this is the paradigm shift required to leverage incoming data from ground-based experiments, which will significantly improve the cosmological parameter constraints from the Planck mission. candl comes with the latest primary and lensing power spectrum data from the South Pole Telescope and Atacama Cosmology Telescope collaborations and will be used as part of the upcoming SPT-3G $TT/TE/EE$ and $\phi\phi$ data releases. Along with the core code, we release a series of auxiliary tools, which simplify common analysis tasks and interface the likelihood with other cosmological software. candl is pip-installable and publicly available on GitHub: https://github.com/Lbalkenhol/candl .

Context. Abundances of s- and r-process elements in Sun-like stars constrain nucleosynthesis in extreme astrophysical events, such as compact binary mergers and explosions of highly magnetised rapidly rotating massive stars. Aims. We measure solar abundances of yttrium (Y) and europium (Eu) using 3D non-local thermal equilibrium (NLTE) models. We use the model to determine the abundance of Y, and also explore the model's ability to reproduce the solar centre-to-limb variation of its lines. In addition, we determine the Eu abundance using solar disc-centre and integrated flux spectra. Methods. We developed an NLTE model of Eu and updated our model of Y with collisional data from detailed quantum-mechanical calculations. We used the IAG spatially resolved high-resolution solar spectra to derive the solar abundances of Y across the solar disc and of Eu for integrated flux and at disc centre using a set of carefully selected lines and a 3D radiation-hydrodynamics model of the solar atmosphere. Results. We find 3D NLTE solar abundances of A(Y)$_{\textrm{3D NLTE}}$=$2.30 \pm 0.03_{\textrm{stat}} \pm 0.07_{\textrm{syst}}$ dex based on observations at all angles and A(Eu)$_{\textrm{3D NLTE}}$=$0.57 \pm 0.01_{\textrm{stat}} \pm 0.06_{\textrm{syst}}$ dex based on the integrated flux and disc-centre intensity. 3D NLTE modelling offers the most consistent abundances across the solar disc, and resolves the problem of severe systematic bias in Y and Eu abundances inherent to 1D LTE, 1D NLTE, and 3D LTE modelling.

The SAAO Cape Town campus was declared a National Heritage Site in December 2018, just short of its 200th anniversary, but is now in a run-down condition. As the former Royal Observatory, it is the oldest scientific institution in South Africa and probably in all Africa. It has a fascinating and well-documented history and surely deserves better. For many years maintenance has been neglected and many of the old telescopes and buildings are in a poor state. They are beginning to show signs of serious decay. Some examples are given.

With the recent radioastronomical detection of cis-trans-carbonic acid (H$_2$CO$_3$) in a molecular cloud toward the galactic center, the more stable but currently unobserved cis-cis conformer is shown here to have strong IR features. While the higher-energy cis-trans-carbonic acid was detected at millimeter and centimeter wavelengths, owing to its larger dipole moment, the vibrational structure of cis-cis-carbonic acid is more amenable to its observation at micron wavelengths. Even so, both conformers have relatively large IR intensities, and some of these fall in regions not dominated by polycyclic aromatic hydrocarbons. Water features may inhibit observation near the 2.75 $\mu$m hydride stretches, but other vibrational fundamentals and even overtones in the 5.5 $\mu$m to 6.0 $\mu$m range may be discernible with JWST data. This work has employed high-level, accurately benchmarked quantum chemical anharmonic procedures to compute exceptionally accurate rotational spectroscopic data compared to experiment. Such performance implies that the IR absorption and even cascade emission spectral features computed in this work should be accurate and will provide the needed reference for observation of either carbonic acid conformer in various astronomical environments.

The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at high redshifts using a large sample of AGNs selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we use more than 600 soft X-ray selected high-z sources in the Chandra Deep fields, the Chandra COSMOS Legacy survey and the XMM-XXL northern field. We derive the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. For modeling the parametric form of the XLF and the absorption function, we use a Bayesian methodology allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity dependent density evolution model cannot be securely ruled out. A large fraction ($60\%)$ of our sources are absorbed by column densities of $\rm N_H \geq 10^{23} cm^{-2} $, while $17$\% of the sources are CTK. Our results favor a scenario where both the ISM of the host and the AGN torus contribute to the obscuration. The derived BHAD is in agreement with the simulations, if one takes into account that the X-ray AGN are hosted by massive galaxies, while it differs from the one derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties.

We present a flare star catalog from four years of non-targeted millimeter-wave survey data from the South Pole Telescope (SPT). The data were taken with the SPT-3G camera and cover a 1500-square-degree region of the sky from $20^{h}40^{m}0^{s}$ to $3^{h}20^{m}0^{s}$ in right ascension and $-42^{\circ}$ to $-70^{\circ}$ in declination. This region was observed on a nearly daily cadence from 2019-2022 and chosen to avoid the plane of the galaxy. A short-duration transient search of this survey yields 111 flaring events from 66 stars, increasing the number of both flaring events and detected flare stars by an order of magnitude from the previous SPT-3G data release. We provide cross-matching to Gaia DR3, as well as matches to X-ray point sources found in the second ROSAT all-sky survey. We have detected flaring stars across the main sequence, from early-type A stars to M dwarfs, as well as a large population of evolved stars. These stars are mostly nearby, spanning 10 to 1000 parsecs in distance. Most of the flare spectral indices are constant or gently rising as a function of frequency at 95/150/220 GHz. The timescale of these events can range from minutes to hours, and the peak $\nu L_{\nu}$ luminosities range from $10^{27}$ to $10^{31}$ erg s$^{-1}$ in the SPT-3G frequency bands.

Previous studies have found that red giants in close binary systems undergoing spin-orbit resonance exhibit an enhanced level of magnetic activity from measurements of the indices of photometric variability $S_{ph}$ and chromospheric emission $S_{CaII}$. Here, we complement the previous works by measuring two other indicators of chromospheric activity: the near-ultraviolet (NUV) excess $\Delta m_{NUV}$ using GALEX data, as well as the H$\alpha$ chromospheric index $S_{H\alpha}$ using spectroscopic data from LAMOST. We consider a sample of 4465 single and binary red giants observed by Kepler, and measure $\Delta m_{NUV}$ and $S_{H\alpha}$ for 842 and 3362 targets, respectively. We investigate the correlations between $\Delta m_{NUV}$, $S_{H\alpha}$, $S_{ph}$ and $S_{CaII}$, which probe magnetic activity at different heights from the photosphere to the upper chromosphere. We also find that red giants exhibiting low-amplitude oscillations tend to exhibit larger $\Delta m_{NUV}$ values, but no larger $S_{H\alpha}$ values. Importantly, we show that red giants in a close-binary configuration with spin-orbit resonance or tidal locking display significantly larger $\Delta m_{NUV}$ values and $S_{H\alpha}$ values than single red giants or red giants in binary systems with no special tidal configuration. This result reinforces previous claims that tidal locking leads to larger magnetic fields. We provide criteria to classify the active red giants (single or binary), based on their rotation period and activity indices. Since $\sim$ 90 millions stars have UV photometric observations from GALEX, and $\sim$ 10 million stars have LAMOST spectra, where the signal-to-noise ratio is higher in the vicinity of the H$\alpha$ line than the CaII H & K lines, the NUV excess and the H$\alpha$ index represent highly valuable activity indicators that could help identifying tidally-interacting stellar systems.

Hybrid neutron stars, the compact objects consisting hadronic matter and strange quark matter, can be considered as the probes for the scalar tensor gravity. In this work, we explore the scalarization of hybrid neutron stars in the scalar tensor gravity. For the hadronic phase, we apply a piecewise polytropic equation of state constrained by the observational data of GW170817 and the data of six low-mass X-ray binaries with thermonuclear burst or the symmetry energy of the nuclear interaction. In addition, to describe the strange quark matter inside the hybrid neutron star, different MIT bag models are employed. We study the effects of the value of bag constant, the mass of s quark, the perturbative quantum chromodynamics correction parameter, and the density jump at the surface of quark-hadronic phase transition on the scalarization of hybrid neutron stars. Our results confirm that the scalarization is more sensitive to the value of bag constant, the mass of s quark, and the density jump compared to the perturbative quantum chromodynamics correction parameter.

The Space-based multi-band astronomical Variable Objects Monitor (SVOM) is a Chinese-French mission dedicated to the study of the transient sky. It is scheduled to start operations in 2024. ECLAIRs is a coded-mask telescope with a large field of view. It is designed to detect and localize gamma-ray bursts (GRBs) in the energy range from 4 keV up to 120 keV. In 2021, the ECLAIRs telescope underwent various calibration campaigns in vacuum test-chambers to evaluate its performance. Between 4 and 8 keV, the counting response of the detection plane shows inhomogeneities between pixels from different production batches. The efficiency inhomogeneity is caused by low-efficiency pixels (LEPs) from one of the two batches, together with high-threshold pixels (HTPs) whose threshold was raised to avoid cross-talk effects. In addition, some unexpected noise was found in the detection plane regions close to the heat pipes. We study the impact of these inhomogeneities and of the heat-pipe noise at low energies on the ECLAIRs onboard triggers. We propose different strategies in order to mitigate these impacts and to improve the onboard trigger performance. We analyzed the data from the calibration campaigns and performed simulations with the ground model of the ECLAIRs trigger software in order to design and evaluate the different strategies. Most of the impact of HTPs can be corrected for by excluding HTPs from the trigger processing. To correct for the impact of LEPs, an efficiency correction in the shadowgram seems to be a good solution. An effective solution for the heat-pipe noise is selecting the noisy pixels and ignoring their data in the 4--8 keV band during the data analysis.

Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectrum of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm$^{\rm -1}$ range originate from CH(D) in-plane and out-of-plane (OOP) wags and bends, as well as bending motions including the C$\equiv$C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm$^{\rm -1}$ for d$_1$-phenylacetylene, a one-to-one match is seen when comparing computations and experiment with all features assigned to combination bands and overtones. The C$\equiv$C stretch observed in phenylacetylene is not observed in d$_1$-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these observations it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general.

Optical phase curves of hot Jupiters can reveal global scattering properties. We implement a Bayesian inference framework for optical phase curves with flux contributions from: reflected light from a potentially inhomogeneous atmosphere, thermal emission, ellipsoidal variations, Doppler beaming, and stellar rotation via a Gaussian process in the time domain. We probe for atmospheric homogeneity and time-variability using the reflected light inferences for highly precise Kepler light curves of five hot Jupiters. We also investigate the scattering properties which constrain the most likely condensates in the inhomogeneous atmospheres. Cross validation prefers inhomogeneous albedo distributions for Kepler-7 b and Kepler-41 b, and a weak preference for inhomogeneity for KOI-13 b. None of the five planets exhibit significant variations in geometric albedo on one-year timescales, in agreement with theoretical expectations. We show that analytic reflected light phase curves with isotropic multiple scattering are in excellent agreement with full Rayleigh multiple scattering calculations, allowing for accelerated and analytic inference. In a case study of Kepler-41 b, we identify perovskite, forsterite, and enstatite as possible scattering species consistent with the reflected light phase curves, with condensate particle radii in the range 0.01-0.1 micron.

We explore the joint detection prospects of short gamma-ray bursts (sGRBs) and their gravitational wave (GW) counterparts by the current and upcoming high-energy GRB and GW facilities from binary neutron star (BNS) mergers. We consider two GW detector networks: (1) A four-detector network comprising LIGO Hanford, Livingston, Virgo, and Kagra, (IGWN4) and (2) a future five-detector network including the same four detectors and LIGO India (IGWN5). For the sGRB detection, we consider existing satellites Fermi and Swift and the proposed all-sky satellite Daksha. Most of the events for the joint detection will be off-axis, hence, we consider a broad range of sGRB jet models predicting the off-axis emission. Also, to test the effect of the assumed sGRB luminosity function, we consider two different functions for one of the emission models. We find that for the different jet models, the joint sGRB and GW detection rates for Fermi and Swift with IGWN4 (IGWN5) lie within 0.07-0.62$\mathrm{\ yr^{-1}}$ (0.8-4.0$\mathrm{\ yr^{-1}}$) and 0.02-0.14$\mathrm{\ yr^{-1}}$ (0.15-1.0$\mathrm{\ yr^{-1}}$), respectively, when the BNS merger rate is taken to be 320$\mathrm{\ Gpc^{-3}~yr^{-1}}$. With Daksha, the rates increase to 0.2-1.3$\mathrm{\ yr^{-1}}$ (1.3-8.3$\mathrm{\ yr^{-1}}$), which is 2-9 times higher than the existing satellites. We show that such a mission with higher sensitivity will be ideal for detecting a higher number of fainter events observed off-axis or at a larger distance. Thus, Daksha will boost the joint detections of sGRB and GW, especially for the off-axis events. Finally, we find that our detection rates with optimal SNRs are conservative, and noise in GW detectors can increase the rates further.

Non-parametric morphology statistics have been used for decades to classify galaxies into morphological types and identify mergers in an automated way. In this work, we assess how reliably we can identify galaxy post-mergers with non-parametric morphology statistics. Low-redshift (z<0.2), recent (t_post-merger < 200 Myr), and isolated (r > 100 kpc) post-merger galaxies are drawn from the IllustrisTNG100-1 cosmological simulation. Synthetic r-band images of the mergers are generated with SKIRT9 and degraded to various image qualities, adding observational effects such as sky noise and atmospheric blurring. We find that even in perfect quality imaging, the individual non-parametric morphology statistics fail to recover more than 55% of the post-mergers, and that this number decreases precipitously with worsening image qualities. The realistic distributions of galaxy properties in IllustrisTNG allow us to show that merger samples assembled using individual morphology statistics are biased towards low mass, high gas fraction, and high mass ratio. However, combining all of the morphology statistics together using either a linear discriminant analysis or random forest algorithm increases the completeness and purity of the identified merger samples and mitigates bias with various galaxy properties. For example, we show that in imaging similar to that of the 10-year depth of the Legacy Survey of Space and Time (LSST), a random forest can identify 89% of mergers with a false positive rate of 17%. Finally, we conduct a detailed study of the effect of viewing angle on merger observability and find that there may be an upper limit to merger recovery due to the orientation of merger features with respect to the observer.

Gravitational-wave astronomy has developed enormously over the last decade with the first detections across different frequency bands, but has yet to access $0.1-10$ $\mathrm{Hz}$ gravitational waves. Gravitational waves in this band are emitted by some of the most enigmatic sources, including intermediate-mass binary black hole mergers, early inspiralling compact binaries, and possibly cosmic inflation. To tap this exciting band, we propose the construction of a detector based on pulsar timing principles, the Artificial Precision Timing Array (APTA). We envision APTA as a solar system array of artificial "pulsars"$-$precision-clock-carrying satellites that emit pulsing electromagnetic signals towards Earth or other centrum. In this fundamental study, we estimate the clock precision needed for APTA to successfully detect gravitational waves. Our results suggest that a clock relative uncertainty of $10^{-17}$, which is currently attainable, would be sufficient for APTA to surpass LISA's sensitivity in the decihertz band and observe $10^3-10^4$ $\mathrm{M}_\odot$ black hole mergers. Future atomic clock technology realistically expected in the next decade would enable the detection of an increasingly diverse set of astrophysical sources, including stellar-mass compact binaries that merge in the LIGO-Virgo-KAGRA band, extreme-mass-ratio inspirals, and Type Ia supernovae. This work opens up a new area of research into designing and constructing artificial gravitational-wave detectors relying on the successful principles of pulsar timing.

The Mercator projection is sometimes confused with another mapping technique, specifically the central cylindrical projection, which projects the Earth's surface onto a cylinder tangent to the equator, as if a light source is at the Earth's center. Accidentally, this misconception is rather close to a truth. The only operation that the map needs is a free bending in a uniform gravitational field if the map's material is dense and soft enough to produce a catenary profile. The north and south edges of the map should be parallel and placed in the same plane at the appropriate distance. In this case, the bent map been projected onto this plane gives the Mercator projection. This property is rather curious, since it allows to make such a sophisticated one-to-one mapping as the Mercator projection using simple tools available in the workroom.

We investigated a bulk viscous fluid universe with cosmological constant {\Lambda} by assuming that the bulk viscosity to be proportional to the Hubble parameter. We found that for an expanding universe, the (relative) matter density will be always greater than a non-zero constant, and tends to this non-zero constant in the future. We show that the bulk viscosity model has a significantly better fitting to the combined SNeIa + CMB + BAO + H(z) data than the {\Lambda}CDM model. Generally, the evolution or values of some cosmological parameters predicted by the bulk viscosity model do not deviate significantly from which are obtained from the {\Lambda}CDM model since the bulk viscosity coefficient obtained from the astronomical observational data is so small. We also made a statefinder analysis of the bulk viscosity model and found that the evolution of the {r, s} parameters behaves in such a way that 0 < s < 1, 0.945 < r <1, indicating the bulk viscosity model is different from the {\Lambda}CDM model.

The Penrose process for the decay of electrically charged particles in a Reissner-Nordstr\"om-anti-de Sitter black hole spacetime is studied. To extract large quantities of energy one needs to mount a recursive Penrose process where particles are confined and can bounce back to suffer ever again a decaying process in the black hole electric ergoregion. In an asymptotically anti-de Sitter (AdS) spacetime, two situations of confinement are possible. One situation uses a reflecting mirror at some radius, which obliges the energetic outgoing particles to return to the decaying point. The other situation uses the natural AdS property that sends back at some intrinsic returning radius those outgoing energetic particles. In addition, besides the conservation laws the decaying process must obey, one has to set conditions at the decaying point for the particles debris. These conditions restrain the possible scenarios, but there are still a great number of available scenarios for the decays. Within these, we choose two scenarios, scenario 1 and scenario 2, that pertain to the masses and electric charges of the final particles. Thus, in the mirror situation we find that scenario 1 leads to a black hole energy factory, and scenario 2 ends in a black hole bomb. In the no mirror situation, i.e., pure Reissner-Nordstr\"om-AdS, scenario 1 leads again to a black hole energy factory, but scenario 2 yields no bomb. This happens because the volume in which the particles are confined increases to infinity along the chain of decays, leading to a zero value of the extracted energy per unit volume and the bomb is demined. The whole treatment performed here involves no backreaction on the black hole mass and electric charge, nevertheless we speculate that the end state of the recursive process is a Reissner-Nordstr\"om-AdS black hole with very short hair, i.e., with one particle at rest at some definite radius.

The BOREXINO experiment has been collecting solar neutrino data since 2011, providing the opportunity to study the variation of the event rate over a decade. We find that at 96 \% C.L., the rate of low energy events shows a time modulation favoring a correlation with a flux from Jupiter. We present a new physics model based on dark matter of mass 1-4 GeV captured by Jupiter that can account for such modulation. We discuss how this scenario can be tested.

We calculate, via variational techniques, single- and two-photon Rydberg microwave transitions, as well as scalar and tensor polarizabilities of sodium atom using the parametric one-electron valence potential, including the spin-orbit coupling. The trial function is expanded in a basis set of optimized Slater-type orbitals, resulting in highly accurate and converged eigen-energies up to $n=60$. We focus our studies on the microwave band 90-150 GHz, due to its relevance to laser excitation in the Earth's upper-atmospheric sodium layer for wavelength-dependent radiometry and polarimetry, as precise microwave polarimetry in this band is an important source of systematic uncertainty in searches for signatures of primordial gravitational waves within the anisotropic polarization pattern of photons from the cosmic microwave background. We present the most efficient transition coefficients in this range, as well as the scalar and tensor polarizabilities compared with available experimental and theoretical data.

We present the early release of the atlas of continuous gravitational waves covering frequencies from 20 Hz to 1500 Hz and spindowns from -5e-10 to 5e-10 Hz/s. Compared to the previous atlas release we have greatly expanded the parameter space, and we now also provide polarization-specific data - both for signal-to-noise ratios and for the upper limits. Continuous wave searches are computationally difficult and take a long time to complete. The atlas enables new searches to be performed using modest computing power. To allow new searches to start sooner, we are releasing this data early, before our followup stages have completed.

Understanding the late-time acceleration of the universe and its subtleties is one of the biggest mysteries in cosmology. A lot of different approaches have been put forward to deal with this, ranging from the conventional cosmological constant to various models of dark energy and beyond. Recently one very interesting approach to explaining the late time acceleration has been put forward, where the the expansion of the universe is driven by mergers with other "baby" universes and has been shown to be quite viable as well from the point of view of recent observational data. So in this work we examine the possibility of various rip scenarios and other future cosmological singularities in such "multiversal" scenario, probing such singularities for the first time in a multi universe scenario. We examine two models of such a baby universe merging cosmology, and show that remarkably no rip scenario or future cosmological singularity is possible in such models.

The structure of the Standard Model (SM) of particle physics points toward grand unified theories (GUTs) where strong and electroweak interactions are unified in a non-Abelian GUT group. The spontaneous breaking of the GUT symmetry to the SM symmetry, together with cosmic inflation, generically leads to metastable topological defects, the most prominent example being cosmic strings. The gravitational-wave background (GWB) produced by a cosmic string network is one of the candidates for an explanation of the GWB recently observed by pulsar timing array (PTA) experiments. We review some properties of the predicted GWB with emphasis on potential implications for GUT model building. The most striking prediction is a GWB in the LIGO-Virgo-KAGRA band that could be discovered in the near future.

Gravitational radiation-reaction phenomena occurring in the dynamics of inspiralling compact binary systems are investigated at the first post-Newtonian order beyond the quadrupole approximation in the context of Einstein-Cartan theory, where quantum spin effects are modeled via the Weyssenhoff fluid. We exploit balance equations for the energy and angular momentum to determine the binary orbital decay until the two bodies collide. Our framework deals with both quasi-elliptic and quasi-circular trajectories, which are then smoothly connected. Key observables like the laws of variation of the orbital phase and frequency characterizing the quasi-circular motion are derived analytically. We conclude our analysis with an estimation of the spin contributions at the merger, which are examined both in the time domain and the Fourier frequency space through the stationary wave approximation.

Neutron stars contain neutron-rich matter with around 5% protons at nuclear saturation density. In this Letter, we consider equilibrium between bulk phases of matter based on asymmetric nuclear matter calculations using chiral effective field theory interactions rather than, as has been done in the past, by interpolation between the properties of symmetric nuclear matter and pure neutron matter. Neutron drip (coexistence of nuclear matter with pure neutrons) is well established, but from earlier work it is unclear whether proton drip (equilibrium between two phases, both of which contain protons and neutrons) is possible. We find that proton drip is a robust prediction of any physically reasonable equation of state, but that it occurs over a limited region of densities and proton fractions. An analytical model based on expanding the energy in powers of the proton density, rather than the neutron excess, is able to account for these features of the phase diagram.

A precise modelling of the dynamics of bubbles nucleated during first-order phase transitions in the early Universe is pivotal for a quantitative determination of various cosmic relics, including the stochastic background of gravitational waves. The equation of motion of the bubble front is affected by the out-of-equilibrium distributions of particle species in the plasma which, in turn, are described by the corresponding Boltzmann equations. In this work we provide a solution to these equations by thoroughly incorporating the non-linearities arising from the population factors. Moreover, our methodology relies on a spectral decomposition that leverages the rotational properties of the collision integral within the Boltzmann equations. This novel approach allows for an efficient and robust computation of both the bubble speed and profile. We also refine our analysis by including the contributions from the electroweak gauge bosons. We find that their impact is dominated by the infrared modes and proves to be non-negligible, contrary to the naive expectations.

We develop a compressible liquid-drop model to describe the crust of neutron stars for which the role of the nuclear clusters, the neutron gas, and the electrons are clearly identified. The novelty relies on the contribution of the neutron gas, which is qualitatively adjusted to reproduce 'ab initio' predictions in dilute neutron matter. We relate the properties of dilute neutron matter to the ones of neutron stars crust and we compare the mean field approximation to an improved approach which better describes dilute neutron matter. The latter is quite sensitive to the unitary limit, a universal feature of Fermi systems having a large value of the scattering length and a small interaction range. While the impact of the accurate description of dilute neutron matter is important in uniform matter (up to 30\% corrections with respect to a mean field calculations), we find a reduction of this impact in the context of the crust of neutron stars due to the additional matter components (nuclear clusters and electrons). In agreement with our previous works, dilute neutron matter is however a necessary ingredient for accurate predictions of the properties of the crust of neutron stars.