Papers for Monday, Jul 10 2023

The controversy "dark matter vs. modified gravity" constitutes a major topic of discussion. It was proposed that dynamical friction could be used to discriminate between the two alternatives. Analytic calculations indicate that, with modified gravity, globular clusters (GCs) of low-mass galaxies experience much stronger dynamical friction than in the equivalent system with Newtonian gravity and dark matter. As a result, in modified gravity the old GCs of low mass galaxies should have already settled in the centers of the galaxies. This is not observed. Here we report on our efforts to verify the analytic results by self-consistent simulations with the MOND-type (modified Newtonian dynamics) gravity. The core stalling mechanism, that was not considered in the analytic calculations, prevents GCs to settle in centers of ultra-diffuse galaxies. For isolated dwarf galaxies, which are gas-rich objects, supernova explosions prevent the GCs from settling.

We analysed the ionised gas kinematics of the dwarf galaxy NGC 4214 using high resolution Fabry-Perot interferometry observations and present a set of narrowband images in the H$\alpha$, [SII] $\lambda$6717 $\r{A}$, [NII] $\lambda$6584 $\r{A}$ and [OIII] $\lambda$5007 $\r{A}$ emission lines. The high-resolution Fabry-Perot observations of the H$\alpha$ emission line, allowed us to derive the velocity field, the velocity dispersion $\sigma$, and the rotation curve of the galaxy. We also present for the first time, three-dimensional kinematic maps of the complexes NGC 4214-I and NGC 4214-II and analysed the kinematics of the ionised gas of two new superbubbles, as well as the supernova remnants previously detected in this galaxy by other authors, in radio, optical and X-ray emission. We computed the expansion velocities of the superbubbles and supernova remnants fitting their velocity profiles and obtained their respective physical parameters. We found that the superbubbles have an expansion velocity of ~50 km s$^{-1}$, dynamical age about $\sim$2 Myr and wind luminosity L$_W$ of ~9X10$^{38}$ erg s$^{-1}$ produced probably by massive stars in OB associations. For supernova remnants, their expansion velocities are between $\sim$48 to $\sim$80 km s$^{-1}$ with ages of about 10$^{4}$ years and kinetic energy of about 10$^{51}$ erg assuming they are in the radiative phase of evolution.

Flattened and kinematically correlated planes of dwarf satellite galaxies have been observed in the Local Volume. The slinging out of satellites during host galaxy mergers has been suggested as a formation mechanism for these peculiar structures. We statistically examined the impact of major mergers on present-time satellite systems for the first time in a full cosmological context using the IllustrisTNG suite of hydrodynamic simulations. Mergers with mass ratios above 1/3 generally have a negligible or adverse impact on the phase-space correlation of observationally motivated satellites. Even high-angular momentum mergers are inefficient at slinging satellites outward due to the extended nature of simulated satellite distributions. Furthermore, any potential merger imprint is partially washed out by post-merger accretion of satellites, while satellites bound to the merging haloes since the merger's beginning are disrupted and stripped of mass - minimising the merger's influence on the present-time distribution of the most massive satellites after 2-5 Gyr. Constraining our sample to satellites bound to their host throughout the full duration of their system's last merger, we recover no particular improvement in their phase-space correlation. Instead, such participant satellites experience a contraction of their radial distribution during and after the merger, resulting in smaller absolute plane heights (but comparable axis ratios). Overall, major mergers do not appear to form correlated planes in a statistical sample. Mergers that efficiently transfer their angular momentum to satellite distributions can marginally enhance their phase-space correlation, but cannot form highly flattened and orbitally coherent configurations as observed in our local Universe.

The star formation rate (SFR) in high redshift galaxies is expected to be time-variable due to competing physical processes. Such stochastic variability might boost the luminosity of galaxies, possibly explaining the over-abundance seen at $z\gtrsim 10$ by JWST. We aim at quantifying the amplitude and timescales of such variability, and identifying the key driving physical processes. We select 245 $z=7.7$ galaxies with stellar mass $5\times 10^{6}\lesssim M_\star/{\rm M}_\odot\lesssim 5\times 10^{10}$ from SERRA, a suite of high-resolution, radiation-hydrodynamic cosmological simulations. After fitting the average SFR trend, $\langle {\rm SFR} \rangle$, we quantify the time-dependent variation, $\delta(t) \equiv \log [\rm SFR/\langle {\rm SFR} \rangle]$ for each system, and perform a periodogram analysis to search for periodicity modulations. We find that $\delta(t)$ is distributed as a zero-mean Gaussian, with standard deviation $\sigma_\delta \simeq 0.24$ (corresponding to a UV magnitude s.d. $\sigma_{\rm UV} \simeq 0.61$) that is independent of $M_\star$. However, the modulation timescale increases with stellar mass: $t_\delta \sim (9, 50, 100)\, \rm Myr$ for $M_\star \sim (0.1, 1, 5)\times 10^9\, {\rm M}_\odot$, respectively. These timescales are imprinted on the SFR by different processes: (i) photoevaporation, (ii) supernova explosions, and (iii) cosmological accretion dominating in low, intermediate, and high mass systems, respectively. The predicted SFR variations cannot account for the required $z\gtrsim 10$ UV luminosity function boost. Other processes, such as radiation-driven outflows clearing the dust, must then be invoked to explain the enhanced luminosity of super-early systems.

Intensity interferometry -- the correlation of spatially separated light intensities -- has historically been an important tool for precision optical astronomical observations. However, due to the extremely narrow field of view, its scope has been limited to studies of the morphology of very bright emission regions, primarily determinations of angular diameters of nearby hot stars. We propose adding an adjustable path extension into the detector optics which creates a primary interference fringe for widely separated sources, allowing maximum source separations parametrically larger than the angular resolution. This Extended-Path Intensity Correlator (EPIC), augmented with advances in single-photon detectors and spectroscopic gratings, would enable ground-based astrometry at microarcsecond-level precision in a field of view as large as several arcseconds. EPIC has the potential to revolutionize astrophysical and cosmological observations requiring high-precision differential astrometry on sources of high surface brightness. We outline how EPIC can be employed to detect the astrometric wobble of Earth-like planets around Sun-like stars at tens to hundreds of parsecs, and expect that EPIC's larger field of view will expand the power of intensity interferometry to a broad range of astronomical applications.

Decaying dark matter (DDM) scenarios have recently re-gained attention due to their potential ability to resolve the well-known clustering (or $S_8$) tension between weak lensing (WL) and cosmic microwave background (CMB) measurements. In this paper, we investigate a well-established model, where the original dark matter (DM) particle decays into a massless and a massive daughter particles. The latter obtains a velocity kick during the decay process resulting in a suppression of the matter power spectrum at scales that are observable with WL shear observations. We perform the first fully nonlinear WL analysis of this two-body decaying dark matter ($\Lambda$DDM) scenario including intrinsic alignment and baryonic feedback processes. We thereby use the cosmic shear band power spectra from the KiDS-1000 data combining it with temperature and polarization data from Planck to constrain the $\Lambda$DDM model. We report new limits on the decay rate and mass splitting parameters that are significantly stronger than previous results, especially for the case of low mass splittings. We also investigate the $S_8$ tension only finding a marginal improvement of 0.3$\sigma$ for $\Lambda$DDM compared to the $\Lambda$CDM case. The improvement is not caused by a shift but a slight bloating of the posterior contours caused by the additional free model parameters. We therefore conclude that the two-body $\Lambda$DDM model does not provide a convincing solution to the $S_8$ tension. Our emulator to model the nonlinear $\Lambda$DDM power spectrum is published as part of the publicly available code DMemu at https://github.com/jbucko/DMemu.

The population of binary black hole mergers observed in gravitational waves, together with astrophysical simulations, can help us to understand the properties of the progenitors and the binary formation mechanisms in different astrophysical scenarios. Here we focus on dynamical formation in star clusters. We use the third gravitational-wave transient catalog (GWTC-3) and Rapster, a rapid code to simulate cluster dynamics, to show that it is possible to construct the single-event likelihood of star cluster properties from individual observations. We find that the measured primary mass in a binary black hole merger correlates with the measured star cluster mass, because the mass spectrum of the primary component increases with the mass of the cluster. This trend may be caused by two physical mechanisms: (i) the more efficient production of hierarchical mergers with primary mass above $\sim 40~M_{\odot}$ for cluster masses of $\gtrsim 10^6~M_{\odot}$, and (ii) the suppression of more massive first-generation binaries, which happens because ejected binaries do not merge within the lookback time for cluster masses of $\lesssim 10^5~M_{\odot}$. The formalism presented here can be generalized to infer the population properties of binary progenitors in more realistic scenarios involving multiple formation channels.

Cloud-wind interactions play an important role in long-lived multiphase flows in many astrophysical contexts. When this interaction is primarily mediated by hydrodynamics and radiative cooling, the survival of clouds can be phrased in terms of the comparison between a timescale that dictates the evolution of the cloud-wind interaction, (the dynamical time-scale $\tau_{\rm dyn}$) and the relevant cooling timescale $\tau_{\rm cool}$. Previously proposed survival criteria, which can disagree by large factors about the size of the smallest surviving clouds, differ in both their choice of $\tau_{\rm cool}$ and (to a lesser extent) $\tau_{\rm dyn}$. Here we present a new criterion which agrees with a previously proposed empirical formulae but is based on simple physical principles. The key insight is that clouds can grow if they are able to mix and cool gas from the hot wind faster than it advects by the cloud. Whereas prior criteria associate $\tau_{\rm dyn}$ with the cloud crushing timescale, our new criterion links it to the characteristic cloud-crossing timescale of a hot-phase fluid element, making it more physically consistent with shear-layer studies. We develop this insight into a predictive expression and validate it with hydrodynamic ENZO-E simulations of ${\sim}10^4\, {\rm K}$, pressure-confined clouds in hot supersonic winds, exploring, in particular, high wind/cloud density contrasts, where disagreements are most pronounced. Finally, we illustrate how discrepancies among previous criteria primarily emerged due to different choices of simulation conditions and cooling properties, and discuss how they can be reconciled.

We report the results of an observational campaign using the Effelsberg 100-m telescope of the pulsars J1746$-$2849, J1746$-$2850, J1746$-$2856 and J1745$-$2912 located in the Central Molecular Zone (CMZ) close to the Galactic centre in order to study rotation measure (RM) variations. We report for the first time the RM value of PSR J1746$-$2850 to be $-12234 \pm 181$ rad m$^{-2}$. This pulsar shows significant variations of RM of $300-400$ rad m$^{-2}$ over the course of months to years that suggest a strongly magnetized environment. The structure function analysis of the RM of PSR J1746$-$2850 revealed a steep power-law index of $1.87_{-0.3}^{+0.4}$ comparable to the value expected for isotropic turbulence. This pulsar also showed large dispersion measure (DM) variation of $\sim 50$ pc cm$^{-3}$ in an event lasting a few months where the RM increased by $\sim 200$ rad m$^{-2}$. The large difference in RM between PSR J1746$-$2849 and PSR J1746$-$2850 despite the small angular separation reveals the presence of a magnetic field of at least 70 $\mu$G in the CMZ and can explain the lack of polarization in the radio images of the region. These results contribute to our understanding of the magnetic field in the CMZ and show similarities between the RM behaviours of these pulsars and some fast radio bursts (FRBs).

In this proceeding, we present the 1-dimensional stellar evolution of two rotating population III (Pop III) star models, each having a mass of 25 M$_{\odot}$ at the zero-age main-sequence (ZAMS). The slowly rotating model has an initial angular rotational velocity of 10 per cent of the critical angular rotational velocity. In contrast, the rapidly rotating model has an initial angular rotational velocity of 70 per cent of the critical angular rotational velocity. As an effect of rotationally enhanced mixing, we find that the rapidly rotating model suffers an enormous mass loss due to the deposition of a significant amount of CNO elements toward the surface after the main-sequence phase. We also display the simulated light curves as these models explode into core-collapse supernovae (CCSNe).

Asteroseismology has transformed stellar astrophysics. Red giant asteroseismology is a prime example, with oscillation periods and amplitudes that are readily detectable with time-domain space-based telescopes. These oscillations can be used to infer masses, ages and radii for large numbers of stars, providing unique constraints on stellar populations in our galaxy. The cadence, duration, and spatial resolution of the Roman galactic bulge time-domain survey (GBTDS) are well-suited for asteroseismology and will probe an important population not studied by prior missions. We identify photometric precision as a key requirement for realizing the potential of asteroseismology with Roman. A precision of 1 mmag per 15-min cadence or better for saturated stars will enable detections of the populous red clump star population in the Galactic bulge. If the survey efficiency is better than expected, we argue for repeat observations of the same fields to improve photometric precision, or covering additional fields to expand the stellar population reach if the photometric precision for saturated stars is better than 1 mmag. Asteroseismology is relatively insensitive to the timing of the observations during the mission, and the prime red clump targets can be observed in a single 70 day campaign in any given field. Complementary stellar characterization, particularly astrometry tied to the Gaia system, will also dramatically expand the diagnostic power of asteroseismology. We also highlight synergies to Roman GBTDS exoplanet science using transits and microlensing.

In this work, we present a quintessential interpretation of having a blue-tilted tensor power spectrum for canonical single-field slow-roll inflation to explain the recently observed NANOGracv 15-year signal of Gravitational Waves (GW). We formulate the complete semi-classical description of cosmological perturbation theory in terms of scalar and tensor modes using the Non-Bunch Davies initial condition. We found that the existence of the blue tilt ($n_t$) within the favoured range $1.2<n_t<2.5$ can be explained in terms of a newly derived consistency relation. Further, we compute a new field excursion formula using the Non-Bunch Davies initial condition, that validates the requirement of Effective Field Theory in the sub-Planckian regime, $|\Delta\phi|\ll M_{\rm pl}$ for the allowed value of the tensor-to-scalar ratio, $r<0.06$ from CMB observations. In our study, we refer to this result as Anti Lyth bound as it violates the well-known Lyth bound originally derived for Bunch Davies initial condition. Further, we study the behaviour of the spectral density of GW and the associated abundance with the frequency, which shows that within the frequency domain $10^{-9}{\rm Hz}<f<10^{-7}{\rm Hz}$ the outcome obtained from our analysis is completely consistent with the NANOGracv 15-year signal. Also, we found that the behaviour of GW spectra satisfies the CMB constraints at the low frequency, $f_*\sim 7.7\times 10^{-17}{\rm Hz}$ corresponding to the pivots scale wave number, $k_*\sim 0.05{\rm Mpc}^{-1}$. Finally, the sharp falling behaviour of the GW spectra within the frequency domain $10^{-7}{\rm Hz}<f<1{\rm Hz}$ validates our theory in the comparatively high frequency regime as well.

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are $C^2$ or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than $\Delta \tilde{\Lambda} = 55$ can be captured by the SpEC code for these equations of state.

IceCube has recently reported the detection of $\sim 1-10 \,{\rm TeV}$ neutrinos from the nearby active galaxy, NGC 1068. The lack of TeV-scale emission from this source suggests that these neutrinos are generated in the dense corona that surrounds NGC 1068's supermassive black hole. In this paper, we present a physical model for this source, including the processes of pair production, pion production, synchrotron, and inverse Compton scattering. We have also performed a new analysis of Fermi-LAT data from the direction of NGC 1068, finding that the gamma-ray emission from this source is very soft but bright at energies below $\sim 1 \, {\rm GeV}$. Our model can predict a gamma-ray spectrum that is consistent with Fermi-LAT observations, but only if the magnetic field within the corona of this active galactic nucleus (AGN) is quite high, namely $B\gtrsim 6 \, {\rm kG}$. To explain the observed neutrino emission, this source must accelerate protons with a total power that is comparable to its intrinsic X-ray luminosity. In this context, we consider two additional nearby active galaxies, NGC 4151 and NGC 3079, which have been identified as promising targets for IceCube.

We use CEERS JWST/NIRCam imaging to measure rest-frame near-IR light profiles of $>$500 $M_\star>10^{10}~M_\odot$ galaxies in the redshift range $0.5<z<2.3$. We compare the resulting rest-frame 1.5-2$\mu$m half-light radii ($R_{\rm{NIR}}$) with stellar half-mass radii (\rmass) derived with multi-color light profiles from CANDELS HST imaging. In general agreement with previous work, we find that $R_{\rm{NIR}}$ and \rmass~are up to 40\%~smaller than the rest-frame optical half-light radius $R_{\rm{opt}}$. The agreement between $R_{\rm{NIR}}$ and \rmass~is excellent, with negligible systematic offset ($<$0.03 dex) up to $z=2$ for quiescent galaxies and up to $z=1.5$ for star-forming galaxies. We also deproject the profiles to estimate \rmassd, the radius of a sphere containing 50\% of the stellar mass. We present the $R-M_\star$ distribution of galaxies at $0.5<z<1.5$, comparing $R_{\rm{opt}}$, \rmass~and \rmassd. The slope is significantly flatter for \rmass~and \rmassd~ compared to $R_{\rm{opt}}$, mostly due to downward shifts in size for massive star-forming galaxies, while \rmass~and \rmassd~do not show markedly different trends. Finally, we show rapid size evolution ($R\propto (1+z)^{-1.7\pm0.1}$) for massive ($M_\star>10^{11}~M_\odot$) quiescent galaxies between $z=0.5$ and $z=2.3$, again comparing $R_{\rm{opt}}$, \rmass~and \rmassd. We conclude that the main tenets of the size evolution narrative established over the past 20 years, based on rest-frame optical light profile analysis, still hold in the era of JWST/NIRCam observations in the rest-frame near-IR.

We present a statistical analysis of the He ii 4686 emission line in the spectra of the black hole and Wolf-Rayet (WR) star of the high-mass X-ray binary IC10 X-1. This line is visibly skewed, and the third moment (skewness) varies with the binary's orbital phase. We describe a new method of extracting such weak/faint features lying barely above a noisy continuum. Using the moments of these features, we have been able to decompose these skewed lines into two symmetric Gaussian profiles as a function of the orbital phase. The astrophysical implications of this decomposition are significant due to the complex nature of wind-accretion stream interactions in such binary systems. Previous studies have already shown a 0.25 phase lag in the radial velocity curve of the star and the X-ray eclipse, which indicates that the He ii emitters might be in the stellar wind, hence not tracing the star's orbital motion. Results from this work further suggest the existence of two separate emitting regions, one in the stellar wind in the shadow of the WR star, and another in the accretion stream that impacts the black hole's outer accretion disk; and the observed skewed He ii lines can be reproduced by superposition of the two corresponding time-dependent Gaussian emission profiles.

The Superpressure Balloon-borne Imaging Telescope (SuperBIT) is a diffraction-limited, wide-field, 0.5 m, near-infrared to near-ultraviolet observatory designed to exploit the stratosphere's space-like conditions. SuperBIT's 2023 science flight will deliver deep, blue imaging of galaxy clusters for gravitational lensing analysis. In preparation, we have developed a weak lensing measurement pipeline with modern algorithms for PSF characterization, shape measurement, and shear calibration. We validate our pipeline and forecast SuperBIT survey properties with simulated galaxy cluster observations in SuperBIT's near-UV and blue bandpasses. We predict imaging depth, galaxy number (source) density, and redshift distribution for observations in SuperBIT's three bluest filters; the effect of lensing sample selections is also considered. We find that in three hours of on-sky integration, SuperBIT can attain a depth of b = 26 mag and a total source density exceeding 40 galaxies per square arcminute. Even with the application of lensing-analysis catalog selections, we find b-band source densities between 25 and 30 galaxies per square arcminute with a median redshift of z = 1.1. Our analysis confirms SuperBIT's capability for weak gravitational lensing measurements in the blue.

Collisionless simulations of structure formation with significant local primordial non-Gaussianities at Mpc scales have shown that a non-Gaussian tail favouring underdensities, with a negative $f_{\rm NL}$ parameter, can significantly change the merging history of galaxy-sized dark matter halos, which then typically assemble later than in vanilla $\Lambda$CDM. Moreover, such a small-scale negative $f_{\rm NL}$ could have interesting consequences for the cosmological $S_8$ tension. Here, we complement our previous work on collisionless simulations with new hydrodynamical simulations of galaxy formation in boxes of 30 Mpc/$h$, using the {\sc RAMSES} code. In particular, we show that all feedback prescriptions being otherwise identical, simulations with a negative $f_{\rm NL} \sim -1000$ on small scales, hence forming galaxies a bit later than in vanilla $\Lambda$CDM, allow to form simulated galaxies with more disky kinematics than in the vanilla case. Therefore, such small-scale primordial non-Gaussianities could potentially help alleviate, simultaneously, tensions in cosmology and galaxy formation. These hydrodynamical simulations on small scales will need to be complemented with larger box simulations with scale-dependent non-Gaussianities, to statistically confirm these trends and explore their observational consequences in further detail.

Mass loss is one of the key parameters that determine stellar evolution. Despite the progress we have achieved over the last decades we still cannot match the observational derived values with theoretical predictions. Even worse, there are certain phases, such as the B[e] supergiants (B[e]SGs) and the Luminous Blue Variables (LBVs), where significant mass is lost through episodic or outburst activity. This leads to various structures around them that permit dust formation, making these objects bright IR sources. The ASSESS project aims to determine the role of episodic mass in the evolution of massive stars, by examining large numbers of cool and hot objects (such as B[e]SGs/LBVs). For this, we initiated a large observing campaign to obtain spectroscopic data for $\sim$1000 IR selected sources in 27 nearby galaxies. Within this project we successfully identified 7 B[e] supergiants (one candidate) and 4 Luminous Blue Variables of which 6 and 2, respectively, are new discoveries. We used spectroscopic, photometric, and light curve information to better constrain the nature of the reported objects. We particularly note the presence of B[e]SGs at metallicity environments as low as 0.14 Z$_{\odot}$.

Star formation, together with the associated chemical and energy feedback, is one of the most important processes in galaxy evolution. The star formation activity in galaxies defines and affects many of their fundamental properties, such as stellar mass, morphology and chemical enrichment levels. Simple models for star formation in cosmological hydrodynamical simulations have shown to be successful in reproducing the star formation rate (SFR) levels and shapes of different types of galaxies. However, with the advent of high-resolution simulations and more detailed observations, more sophisticated star formation models are needed; in particular, to better understand the relation between star formation and the amount of gas in the atomic and molecular phases. In this work, we apply a novel star formation model, recently developed to work in the context of hydrodynamical simulations, to the study of the SFR in Milky Way-mass galaxies. The new implementation describes the formation of molecular hydrogen from atomic material, considering also possible dependencies with the chemical abundance of the gas. This allows to implement various star formation models, where the SFR of a gas cloud is determined by the atomic and/or molecular gas phases, and to compare their predictions to recent observational results.

Fast Radio Bursts (FRBs) are a powerful and mysterious new class of transient that are luminous enough to be detected at cosmological distances. By associating FRBs to host galaxies, we can measure intrinsic and environmental properties that test FRB origin models, in addition to using them as precise probes of distant cosmic gas. The 110-antenna Deep Synoptic Array (DSA-110) is a radio interferometer built to maximize the rate at which it can simultaneously detect and localize FRBs. Here, we present the first sample of FRBs and host galaxies discovered by the DSA-110. This sample of 11 FRBs is the largest uniform sample of localized FRBs to date and is selected based on association to host galaxies identified in optical imaging by Pan-STARRS1 and follow-up spectroscopy at the Palomar and Keck observatories. These FRBs have not been observed to repeat and their radio properties (dispersion, temporal scattering, energy) are similar to that of the known non-repeating FRB population. Most host galaxies have ongoing star formation, as has been identified before for FRB hosts. In contrast to prior work, a large fraction (four of eleven) of the new sample are more massive than 10$^{11}$\ M$_{\odot}$ and most had elevated star formation rates more than 100 Myr in their past. The distribution of star-formation history across this host-galaxy sample shows that the delay-time distribution is wide, spanning from $\sim100$\,Myr to $\sim10$\,Gyr. This requires the existence of one or more progenitor formation channels associated with old stellar populations, such as the binary evolution of compact objects.

This paper reports estimated stellar parameters of 1,153 Kepler red giant branch stars determined with asteroseismic modeling. We use radial-mode oscillation frequencies, gravity-mode period spacings, Gaia luminosities, and spectroscopic data to characterize these stars. Compared with previous studies, we find that the two additional observed constraints, i.e., the gravity-mode period spacing and luminosity, significantly improve the precision of fundamental stellar parameters. The typical uncertainties are 2.9% for the mass, 11% for the age, 1.0% for the radius, 0.0039 dex for the surface gravity, and 0.5\% for the helium core mass, making this the best-characterized large sample of red-giant stars available to date. With better characterizations for these red giants, we recalibrate the seismic scaling relations and study the surface term on the red-giant branch. We confirm that the surface term depends on the surface gravity and effective temperature, but there is no significant correlation with metallicity.

We investigate the influence of magnetic field amplification on the core-collapse supernovae in highly magnetized progenitors through three-dimensional simulations. By considering rotating models, we observe a strong correlation between the exponential growth of the magnetic field in the gain region and the initiation of shock revival, with a faster onset compared to the non-rotating model. We highlight that the mean magnetic field experiences exponential amplification as a result of $\alpha$-effect in the dynamo process, which works efficiently with the increasing kinetic helicity of the turbulence within the gain region. Our findings indicate that the significant amplification of the mean magnetic fields leads to the development of locally intense turbulent magnetic fields, particularly in the vicinity of the poles, thereby promoting the revival of the shock by neutrino heating.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) has discovered more than 650 new pulsars, which account for 20% of our known Galactic pulsar population. In this paper, we estimate the prospect of a pulsar survey with a radio telescope array to be planned -- the FAST Array (FASTA), consists of six "FAST-type" telescopes. Such a sensitive radio telescope array would be a powerful instrument in probing the pulsar population deep into our Galaxy as well as in nearby galaxies. We simulate the FASTA pulsar discovery prospects with different Galactic pulsar population models and instrumental parameter combinations. We find that FASTA could detect tens of thousands of canonical pulsars and well-over thousands of millisecond pulsars. We also estimate the potential yield if the FASTA is used to search for pulsars from the nearby spiral galaxy M31, and find that it would probably discover around a hundred new radio pulsars.

Radio galaxy morphological classification is one of the critical steps when producing source catalogues for large-scale radio continuum surveys. While many recent studies attempted to classify source radio morphology from survey image data using deep learning algorithms (i.e., Convolutional Neural Networks), they concentrated on model robustness most time. It is unclear whether a model similarly makes predictions as radio astronomers did. In this work, we used Local Interpretable Model-agnostic Explanation (LIME), an state-of-the-art eXplainable Artificial Intelligence (XAI) technique to explain model prediction behaviour and thus examine the hypothesis in a proof-of-concept manner. In what follows, we describe how \textbf{LIME} generally works and early results about how it helped explain predictions of a radio galaxy classification model using this technique.

Blazars emit across all electromagnetic wavelengths. While the so-called one-zone model has described well both quiescent and flaring states, it cannot explain the radio emission and fails in more complex data sets, such as AP Librae. In order to self-consistently describe the entire electromagnetic spectrum emitted by the jet, extended radiation models are necessary. Notably, kinetic descriptions of extended jets can provide the temporal and spatial evolution of the particle species and the full electromagnetic output. Here, we present the initial results of a newly developed hadro-leptonic extended-jet code: ExHaLe-jet. As protons take much longer than electrons to lose their energy, they can transport energy over much larger distances than electrons and are therefore essential for the energy transport in the jet. Furthermore, protons induce injection of additional pairs through pion and Bethe-Heitler pair production, which can explain a dominant leptonic radiation signal while still producing neutrinos. In this talk, we discuss the differences between leptonic and hadronic dominated SED solutions, the SED shapes, evolution along the jet flow, and jet powers. We also highlight the important role of external photon fields, such as the accretion disk and the BLR.

The Small Magellanic Cloud (SMC) is the host of a rich system of globular clusters (GCs) that span a wide age range. The chemical composition of the SMC clusters is still poorly understood, despite their significance to chemical evolution studies. Here, we provide the first detailed chemical study of evolved giants in three distinct clusters, NGC 121 (10.5 Gyr), NGC 339 (6 Gyr), and NGC 419 (1.4 Gyr). Results are based on high-resolution spectra obtained with FLAMES at the Very Large Telescope. The chemical fingerprints of these clusters closely resemble those of SMC field stars, supporting the SMC's specific history of chemical enrichment relative to the Milky Way. The approximately solar-scaled [alpha/Fe] observed in all three clusters, independent of their [Fe/H], demonstrate the SMC's low star formation efficiency. Compared to their Milky Way counterparts, elements primarily produced by massive stars are severely underrepresented. Particularly, the young cluster NGC 419's extremely low [Zn/Fe] shows that hypernovae have contributed relatively little during the past two Gyr. The three GCs have high [Eu/Fe] values regardless of their age. This suggests that the production of the r-process elements in the SMC was extremely efficient up to 1.5 Gyr ago, with an enrichment timescale comparable to that from Type Ia supernovae. When the properties of the oldest SMC object NGC 121 are compared to those of in-situ Milky Way clusters and accreted clusters linked to the Gaia-Enceladus merger event, it is shown that the SMC had already attained the same metallicity as Gaia-Enceladus but with lower [Fe/H] ratios at the age of NGC 121. This suggests that the chemical enrichment histories of the early SMC and Gaia-Enceladus differed, and that the SMC probably had a lower mass in its early ages than Gaia-Enceladus.

Stellar clusters (SC) are fundamental building blocks of galaxies and are among the most studied astronomical objects in the Cosmos. The recent association of diffuse $\gamma$-ray emission detected by different experiments with a dozen young SCs suggests the presence of some process able to accelerate particles at least up to hundreds of TeV. In this Ph.D. thesis, we investigate the capability of young massive stellar clusters (YMSC) to produce cosmic rays under the assumption that particles are accelerated at the cluster wind termination shock. The study is divided into three parts. First, we focus on the specific case of Cygnus OB2. We model the observed $\gamma$-ray emission (in a pure hadronic scenario) assuming different models for particle propagation in the neighborhood of the cluster. We found that particles accelerated by Cygnus OB2 can account for both the $\gamma$-ray spectrum and the radial morphology at very high energy. In the second part, we compute the diffuse $\gamma$-ray emission expected by the unresolved population of Galactic YMSC. For this purpose, we build a synthetic population of YMSC based on the properties of local SC. Under the assumption of a pure hadronic emission, we found that YMSCs can significantly contribute to the observed diffuse $\gamma$-ray emission at a few TeV. The final part of the work is dedicated to understand the impact of CRs produced by SCs on the ionization rate of molecular clouds close to those SCs. We found that the ionization rate can significantly differ from the expected value in clouds located in the unperturbed interstellar medium. We show that the measured value of ionization rate, paired with $\gamma$-ray observations, can be used to constrain particle diffusion in the vicinity of the stellar cluster.

The diffusion of gas through porous material is important to understand the physical processes underlying cometary activity. We study the diffusion of a rarefied gas (Knudsen regime) through a packed bed of monodisperse spheres via experiments and numerical modelling, providing an absolute value of the diffusion coefficient and compare it to published analytical models. The experiments are designed to be directly comparable to numerical simulations, by using precision steel beads, simple geometries, and a trade-off of the sample size between small boundary effects and efficient computation. For direct comparison, the diffusion coefficient is determined in Direct Simulation Monte Carlo (DSMC) simulations, yielding a good match with experiments. This model is further-on used on a microscopic scale, which cannot be studied in experiments, to determine the mean path of gas molecules and its distribution, and compare it against an analytical model. Scaling with sample properties (particle size, porosity) and gas properties (molecular mass, temperature) is consistent with analytical models. As predicted by these, results are very sensitive on sample porosity and we find that a tortuosity $q(\varepsilon)$ depending linearly on the porosity $\varepsilon$ can well reconcile the analytical model with experiments and simulations. Mean paths of molecules are close to those described in the literature, but their distribution deviates from the expectation for small path lengths. The provided diffusion coefficients and scaling laws are directly applicable to thermophysical models of idealised cometary material.

We present a study of the co-evolution of a population of primordial star-forming minihalos at Cosmic Dawn. In this study, we highlight the influence of individual Population III stars on the ability of nearby minihalos to form sufficient molecular hydrogen to undergo star formation of their own. In the absence of radiation, we find the minimum halo mass required to bring about collapse and star formation to be 10^5 Msun, which then increases to 10^6 Msun after two stars have formed. We find an inverse relationship between the mass of a halo and the time required for it to recover its molecular gas after being disrupted by radiation from a nearby star. We also take advantage of the extremely high resolution to investigate the effects of major and minor mergers on the gas content of star-forming minihalos. Contrary to previous claims of fallback of supernova ejecta, we find that minihalos evacuated after hosting Pop III stars primarily recover gas through mergers with undisturbed halos. We identify an intriguing type of major merger between recently evacuated halos and gas-rich ones, finding that these "dry" mergers accelerate star formation instead of suppressing it like their low redshift counterparts. We attribute this to the gas-poor nature of one of the merging halos resulting in no significant rise in temperature or turbulence and instead inducing a rapid increase in central density and hydrostatic pressure. This constitutes a novel formation pathway for Pop III stars and establishes major mergers as potentially the primary source of gas, thus redefining the role of major mergers at this epoch.

Context: Tidal and rotational deformation of fluid giant extra-solar planets may impact their transit light curves, making the $k_2$ Love number observable in the upcoming years. Studying the sensitivity of $k_2$ to mass concentration at depth is thus expected to provide new constraints on the internal structure of gaseous extra-solar planets. Aims: We investigate the link between the mean polar moment of inertia $N$ of a fluid, stably layered extra-solar planet and its $k_2$ Love number, aiming at obtaining analytical relationships valid, at least, for some particular ranges of the model parameters. We also seek a general, approximate relationship useful to constrain $N$ once observations of $k_2$ will become available. Methods: For two-layer fluid extra-solar planets, we explore the relationship between $N$ and $k_2$ by analytical methods, for particular values of the model parameters. We also explore approximate relationships valid over all the possible range of two-layer models. More complex planetary structures are investigated by the semi-analytical propagator technique. Results: A unique relationship between $N$ and $k_2$ cannot be established. However, our numerical experiments show that a `rule of thumb' can be inferred, valid for complex, randomly layered stable planetary structures. The rule robustly defines the upper limit to the values of $N$ for a given $k_2$, and agrees with analytical results for a polytrope of index one and with a realistic non-rotating model of the tidal equilibrium of Jupiter.

Metal-poor stars are key for studying the formation and evolution of the Galaxy. Evidence of the early mergers that built up the Galaxy remains in the distributions of abundances, kinematics, and orbital parameters of its stars. Several substructures resulting from these mergers have been tentatively identified in the literature. We conduct a global analysis of the chemodynamic properties of metal-poor stars. Our aim is to identify signs of accreted and in situ stars in different regions of the parameter space and to investigate their differences and similarities. We selected a sample of about 6600 metal-poor stars with [Fe/H] $\leq$ -0.8 from DR3 of the GALAH survey. We used unsupervised machine learning to separate stars in a parameter space made of two normalised orbital actions, plus [Fe/H] and [Mg/Fe], without additional a priori cuts on stellar properties. We divided the halo stars in four main groups. All groups exhibit a significant fraction of in situ contamination (ISC). Accreted stars of these groups have very similar chemical properties, except for those of the group of stars with very retrograde orbits. This points to at most two main sources of accreted stars in the current sample, the major one related to Gaia-Enceladus (GE) and the other possibly related to Thamnos and/or Sequoia. Stars of GE are r-process enriched at low metallicities, but a contribution of the s-process appears with increasing metallicity. A flat trend of [Eu/Mg] as a function of [Fe/H] suggests that only core collapse supernovae contributed to r-process elements in GE. To better characterise accreted stars in the low metallicity regime, high precision abundances and guidance from chemical evolution models are needed. It is possible that ISC in samples of accreted stars has been underestimated. This can have important consequences for attempts to estimate the properties of the original systems.

The origin of current angular momentum (AM) of the black hole (BH) in X-ray binary (XRB) is still unclear, which is related with the birth and/or the growth of the BH. Here we collect the spin parameters $a_{*}$ measured in BH XRBs and find an apparent bimodal distribution centered at $\sim$ 0.17 and 0.83. We find a positive relation between the spin parameter and the orbital period/orbital separation through combining distinct XRB categories, including neutron star (NS) low-mass X-ray binaries (LMXBs), Roche-lobe overflow (RLOF) BH XRBs and wind-fed BH XRBs. It seems that the AM of the compact star and the binary orbit correlates by combining the different XRB systems. These positive relations imply that accretion process is a common mechanism for spinning up the compact star in these diverse XRB systems. We infer that the low and high spin BH XRBs may experience different evolution and accretion history, which corresponds to the bimodal distribution of the BH spin parameters. The low spin BHs ($a_{*}<0.3$) are similar to the NS LMXBs, the compact star of which is spun-up by the low-level accretion, and the high spin BHs ($a_{*}>0.5$) had experienced a short hypercritical accretion ($\gg \dot{M}_\mathrm{Edd}$) period, during which, the BH spin dramatically increased.

Neutrino interactions with stellar material are widely believed to be fundamental to the explosion of massive stars. However, this important process has remained difficult to confirm observationally. We propose a new method to verify it using X-ray observations of the supernova remnant Cassiopeia A. The elemental composition in its Fe-rich ejecta that could have been produced at the innermost region of the supernova, where neutrinos are expected to interact, allows us to examine the presence of neutrino interactions. Here we demonstrate that the amount of Mn produced without neutrino nucleosynthesis processes (i.e., the $\nu$- and $\nu$p-process) is too small to explain the Mn/Fe mass ratio we measure (0.14--0.67\%). This result supports the operation of significant neutrino interactions in the Cassiopeia A supernova. If the observed Mn/Fe mass ratio purely reflects the production at the innermost region of the supernova, this would be the first robust confirmation of neutrino-matter interactions in an individual supernova. We further show that the Mn/Fe mass ratio has the potential to constrain supernova neutrino parameters (i.e., total neutrino luminosity, neutrino temperature). Future spatially-resolved, high-resolution X-ray spectroscopy will allow us to investigate the details of neutrino-supernova astrophysics through its signatures in elemental composition not only in Cassiopeia A but also in other remnants.

We present the detailed fundamental stellar parameters of the close visual binary system; HD39438 for the first time. We used Al-Wardat's method for analyzing binary and multiple stellar systems (BMSSs). The method implements Kurucz's plane parallel model atmospheres to construct synthetic spectral energy distributions for both components of the system. It then combines the results of the spectroscopic analysis with the photometric analysis and then compares them with the observed ones to construct the best synthetic spectral energy distributions for the combined system. The analysis gives the precise fundamental parameters of the individual components of the system. Based on the positions of the components of HD39438 on the H-R diagram, and evolutionary and isochrones tracks, we found that the system belongs to the main sequence stars with masses of 1.24 and 0.98 solar masses for the components A and B, respectively, and age of 1.995 Gyr for both components. The main result of HD39438 is new dynamical parallax, which is estimated to be 16.689+- 0.03 mas.

The matter power spectrum $P(k)$ is one of the main quantities connecting observational and theoretical cosmology. Although for a fixed redshift this can be numerically computed very efficiently by Boltzmann solvers, an analytical description is always desirable. However, accurate fitting functions for $P(k)$ are only available for the concordance model. Taking into account that forthcoming surveys will further constrain the parameter space of cosmological models, it is also of interest to have analytical formulations for $P(k)$ when alternative models are considered. Here, we use the genetic algorithms, a machine learning technique, to find a parametric function for $P(k)$ considering several possible effects imprinted by modifications of gravity. Our expression for the $P(k)$ of modified gravity shows a mean accuracy of around 1-2% when compared with numerical data obtained via modified versions of the Boltzmann solver CLASS, and thus it represents a competitive formulation given the target accuracy of forthcoming surveys.

We present a comprehensive physical and chemical study of the fragmentation and star formation activity towards the massive clump AGAL G338.9188+0.5494 harbouring the extended green object EGO 338.92+0.55(b). The presence of an EGO embedded in a massive clump, suggests, at clump scale, that high-mass star formation is occurring. The main goal of this work is to find evidence of such high-mass star formation, but at core scale. Using mm observations of continuum and lines obtained from the ALMA database at Bands 6 and 7, we study the substructure of the massive clump. The angular resolution of the data is about 0.5'', which allow us to resolve structures of about 0.01pc ($\sim$ 2000 au) at the distance of 4.4 kpc. The continuum emission at 340 GHz reveals that the molecular clump is fragmented in five cores, labeled from C1 to C5. The $^{12}$CO J=3--2 emission shows the presence of molecular outflows related to three of them. The analysis of the CH$_3$CN and CH$_3$CCH emissions suggests temperatures of about 340 and 72~K, respectively, for C1, showing that the methyl cyanide would trace a gas layer closer to the protostar than the methyl acetylene. The obtained mass of core C1 ranges from 3 to 10 M$_{\odot}$. We found that the discovered molecular outflow arising from core C1 should be the main responsible for the 4.5 $\mu$m extended emission. The average mass and energy of such a molecular outflow is about 0.5 M$_{\odot}$~and $10^{46}$~erg, respectively, which suggest that 10 M$_{\odot}$ is the most likely mass value for core C1. Additionally we found that the region is chemically very rich with several complex molecular species. Particularly, from the analysis of the CN emission we found strong evidence that such a radical is indirectly tracing the molecular outflows, more precisely the border of the cavity walls carved out by such outflows.

We investigate the high energy emission activities of two bright gamma-ray pulsars, PSR~J2021+4026 and Vela. For PSR~J2021+4026, the state changes in the gamma-ray flux and spin-down rate have been observed. We report that the long-term evolution of the gamma-ray flux and timing behavior of PSR~J2021+4026 suggests a new gamma-ray flux recovery at around MJD~58910 and a flux decrease around MJD~59500. During this epoch, the staying time, the gamma-ray flux difference and spin-down rate are smaller than previous epochs in the same state. The waiting timescale of the quasi-periodic state changes is similar to the waiting timescale of the glitch events of the Vela pulsar. For the Vela pulsar, the quench of the radio pulse was in a timescale of $\sim0.2$~s after the 2016 glitch, and the glitch may disturb the structure of the magnetosphere. Nevertheless, we did not find any evidence for a long-term change in the gamma-ray emission properties using years of $Fermi$-LAT data, and therefore, no long-term magnetosphere structural change. We also conduct searching for photons above 100~GeV using 15-year $Fermi$-LAT data, and found none. Our results provide additional information for the relation between the state change of the gamma-ray emission and the glitch event.

{The variability study of 6.7\,GHz methanol masers has become a useful way to improve our understanding of the physical conditions in high-mass star-forming regions.} {Based on the single-dish monitoring using the Irbene telescopes, we selected three sources with close sky positions.} {We imaged them using the European Very Long Baseline Interferometer Network and searched available data on VLBI archives to follow detailed changes in their structures and single maser spot variability.} {All three targets show a few groups of maser cloudlets of a typical size of 3.5\,mas and the majority of them show linear or arched structures with velocity gradients of order 0.22\kms\,mas$^{-1}$. The cloudlets and overall source morphologies are remarkably stable on time scales of 7-15\,yr supporting a scenario of variability due to changes in the maser pumping rate.}

Recent research has emphasized the benefits of accurately reconstructing the initial Lagrangian positions of biased tracers from their positions at a later time, to gain cosmological information. A weighted semi-discrete optimal transport algorithm can achieve the required accuracy, provided the late-time positions are known, with minimal information about the background cosmology. The algorithm's performance relies on knowing the masses of the biased tracers, and depends on how one models the distribution of the remaining mass that is not associated with these tracers. We demonstrate that simple models of the remaining mass result in accurate retrieval of the initial Lagrangian positions, which we quantify using pair statistics and the void probability function. This is true even if the input positions are affected by redshift-space distortions. The most sophisticated models assume that the masses of the tracers, and the amount and clustering of the missing mass are known; we show that the method is robust to realistic errors in the masses of the tracers and remains so as the model for the missing mass becomes increasingly crude.

Spherical flows are a classic problem in astrophysics which are typically studied from a global perspective. However, much like with accretion discs, there are likely many instabilities and small scale phenomena which would be easier to study from a local perspective. For this purpose, we develop a local model for a spherically contracting/expanding gas cloud, in the spirit of the shearing box, $\beta$-plane and expanding box models which have had extensive use in studies of accretion discs, planets and stellar winds respectively. The local model consists of a, spatially homogeneous, periodic box with a time varying aspect ratio, along with a scale factor (analogous to that in FRW/Newtonian cosmology) relating the box coordinates to the physical coordinates of the global problem. We derive a number of symmetries and conservation laws exhibited by the local model. Some of these reflect symmetries of the periodic box, modified by the time dependant geometry, while others are local analogues for symmetries of the global problem. The energy, density and vorticity in the box also generically increase(/decrease) as a consequence of the collapse(/expansion). We derive a number of nonlinear solutions, including a local analogue of uniform density zonal flows, which grow as a consequence of angular momentum conservation. Our model is closely related to the accelerated expanding box model of Tenerani \& Velli and is an extension of the isotropic model considered by Robertson \& Goldreich.

SW Sex stars are an informal sub-class of eclipsing nova-like cataclysmic variables. We report 934 new eclipse times measured over the past 17 years for HS 0728+6738 (V482 Cam), SW Sex, DW UMa, HS 0129+2933 (TT Tri), V1315 Aql, PX And, HS 0455+8315, HS 0220+0603, BP Lyn, BH Lyn, LX Ser, UU Aqr, V1776 Cyg, RW Tri, 1RXS J064434.5+334451, AC Cnc, V363 Aur, and BT Mon. When combined with published eclipse times going back in some cases many decades, we show that these binary systems exhibit a range of behaviors, including increasing, decreasing, and possibly oscillating orbital periods. Nevertheless, the duration of these observations is still not long enough to be able to make reliable quantitative statements about their long term behaviors. In addition to these long term trends, we also observed rapid and unusual decreases in the orbital periods of SW Sex and RW Tri during 2017 and 2018, respectively.

In this work, the pulsation analysis is performed on 83 high-amplitude $\delta$ Scuti stars, which have been observed by the Transiting Exoplanet Survey Satellite (TESS). The results show that 49 of these HADS show single-mode pulsation, 27 of them show radial double-modes pulsation (in which 22 of them pulsate with the fundamental and first overtone modes and 5 of them pulsate with the first and second overtone modes), and 7 of them show radial triple-modes pulsation (3 of which are newly confirmed triple-mode HADS). The histogram of the fundamental periods and the ratios between the fundamental and first overtone periods show bimodal structures, which might be caused by the stellar evolution in this specific phase. Most of the radial triple-mode HADS have a fundamental amplitude of 41-54 mmag, and 50% of them have similar amplitudes of the fundamental and first overtone pulsation modes. All these hints require further confirmation not only in observations with more HADS samples, but also in theoretical models with suitable treatments of stellar evolution and pulsation.

We investigate the peculiar nature of strange stars through an analysis of different quark models, i.e. vBag model and CFL model equation of states at different parameter sets, and focus on understanding the equation of state governing the intriguing central compact object (CCO) within the supernova remnant HESS J1731-347, with a mass and radius of $M = 0.77^{+0.20}_{-0.17} M_{\odot}$ and $R = 10.4^{+0.86}_{-0.78}$ km, respectively. Additionally, we compare the radial oscillations of two models to determine the frequency of the HESS J1731-347 compact object at its maximum mass. The frequencies of radial oscillations are computed for each of the four EoSs considered. In total, the 10 lowest radial frequencies for each of those EoSs have been computed. By delving into these aspects, we aim at deepening our understanding of strange stars and their connection to the observed HESS J1731-347 mass-radius relationship.

We provide the first comprehensive study of hyperons in neutron star mergers and quantify their specific impact. We discuss the thermal behavior of hyperonic equations of state~(EoSs) as a distinguishing feature from purely nucleonic models in the remnants of binary mergers. Finite temperature enhances the production of hyperons, which leads to a reduced pressure as highly degenerate nucleons are depopulated. This results in a characteristic increase of the dominant postmerger gravitational-wave frequency by up to $\sim150$~Hz compared to purely nucleonic EoS models. By our comparative approach we can directly link this effect to the occurrence of hyperons. Although this feature is generally weak, it is in principle measurable if the EoS and stellar parameters of cold neutron stars are sufficiently well determined. Considering that the mass-radius relations of purely nucleonic and hyperonic EoSs may be indistinguishable and the overall challenge to infer the presence of hyperons in neutron stars, these findings are important as a new route to answer the outstanding question about hyperonic degrees of freedom in high-density matter.

We report on the spectroscopic confirmation of 68 new bright ($G=13.5-17.2$ mag) and blue (pre-)white dwarfs (WDs). This finding has allowed us to almost double the number of the hottest ($T_{\mathrm{eff}} \geq 60$kK) known WDs brighter than $G=16$ mag. We increased the number of known ultra-high excitation (UHE) WDs by 20%, found one unambiguous close binary system consisting of one DA WD with an irradiated low-mass companion, one DAO, and one DOA WD that are likely in their transformation phase of becoming pure DA WDs, one rare, naked O(H) star, two DA and two DAO WDs with $T_{\mathrm{eff}}$ possibly in excess of 100kK, three new DOZ WDs, and three of our targets are central stars of (possible) planetary nebulae. Using non-local thermodynamic equilibrium models, we derived the atmospheric parameters of these stars and by fitting their spectral energy distribution we derived their radii, luminosities, and gravity masses. In addition, we derived their masses in the Kiel and Hertzsprung-Russell diagram (HRD). We find that Kiel, HRD, and gravity mass agree only in half of the cases. This is not unexpected and we attribute this to the neglect of metal opacities, possibly stratified atmospheres, as well as possible uncertainties of the parallax zero point determination. Furthermore, we carried out a search for photometric variability in our targets using archival data, finding that 26% of our targets are variable. This includes 15 new variable stars, with only one of them being clearly an irradiation effect system. Strikingly, the majority of the variable stars exhibit non-sinusoidal light-curve shapes, which are unlikely explained in terms of close binary systems. We propose that a significant fraction of all (not just UHE) WDs develop spots when entering the WD cooling phase. We suggest that this could be related to the on-set of weak magnetic fields and possibly diffusion.

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($\tau_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe.

We do a detailed analysis of a well-theoretically motivated interacting dark energy scalar field model with a time-varying interaction term. Using current cosmological datasets from CMB, BAO, Type Ia Supernova, $H(z)$ measurements from cosmic chronometers, angular diameter measurements from Megamasers, growth measurements, and local SH0ES measurements, we found that dark energy component may act differently than a cosmological constant at early times. The observational data also does not disfavor a small interaction between dark energy and dark matter at late times. When using all these datasets in combination, our value of $H_0$ agrees well with SH0ES results but in 3.5$\sigma$ tension with Planck results. We also did AIC and BIC analysis, and we found that the cosmological data prefer coupled quintessence model over $\Lambda$CDM, although the chi-square per number of degrees of freedom test prefers the latter.

We examine the dependence of polycyclic aromatic hydrocarbon (PAH) band intensity ratios as a function of the average number of carbon atoms and assess their effectiveness as tracers for PAH size, utilising the data, models, and tools provided by the NASA Ames PAH Infrared Spectroscopic Database. To achieve this, we used spectra from mixtures of PAHs of different ionisation fractions, following a size distribution. Our work, congruent with earlier findings, shows that band ratios that include the 3.3 ${\mu}$m PAH band provide the best PAH size tracers for small-to-intermediate sized PAHs. In addition, we find that band ratios that include the sum of the 15-20 ${\mu}$m PAH features (I$_{\Sigma_{15-20}}$) and the 6.2 or 7.7 ${\mu}$m bands also serve as good tracers for PAH size in the case of small-to-intermediate sized PAHs, for objects under a similar PAH size distribution as with the presented models. For different PAH size distributions, the application of a scaling factor to the I$_{6.2}$/I$_{\Sigma_{15-20}}$ ratio can provide estimates for the size of the small-to-intermediate PAH population within sources. Employment of the I$_{6.2}$/I$_{\Sigma_{15-20}}$ and I$_{7.7}$/I$_{\Sigma_{15-20}}$ ratios can be of particular interest for JWST observations limited only to $\sim$ 5-28 ${\mu}$m MIRI(-MRS) coverage.

We examine the backreaction effect of the enhanced small-scale scalar perturbations from the sound speed resonance (SSR) mechanism for primordial black hole formation in Dirac-Born-Infeld (DBI) inflation. Within the perturbative regime, the backreaction effect of perturbations on the background dynamics can be described by an effective action after integrating out the perturbation sector. Starting with the effective field theory of a specific DBI inflation model that incorporates SSR, we obtain the one-loop effective action by integrating out the scalar perturbations at the quadratic level. Using the effective Friedmann equations derived from this one-loop effective action, we solve the Hubble parameter with backreaction and the effective perturbation dynamics on this background as well. Our numerical findings reveal that, for a viable parameter space, the backreaction effect results in a relative correction to the Hubble parameter of approximately $10^{-6}$, whereas the relative correction to the slow-roll parameter can vary between $-0.4$ and $0.1$, before gradually converging to $10^{-6}$. Furthermore, our results show that the backreaction effect on SSR sound speed causes a slight reduction in the resonant peak of curvature power spectrum, which can be negligible and does not spoil the SSR mechanism.

A considerable effort has been dedicated recently to the construction of generic equations of state (EOSs) for matter in neutron stars. The advantage of these approaches is that they can provide model-independent information on the interior structure and global properties of neutron stars. Making use of more than $10^6$ generic EOSs, we asses the validity of quasi-universal relations of neutron star properties for a broad range of rotation rates, from slow-rotation up to the mass-shedding limit. In this way, we are able to determine with unprecedented accuracy the quasi-universal maximum-mass ratio between rotating and nonrotating stars and reveal the existence of a new relation for the surface oblateness, i.e., the ratio between the polar and equatorial proper radii. We discuss the impact that our findings have on the imminent detection of new binary neutron-star mergers and how they can be used to set new and more stringent limits on the maximum mass of nonrotating neutron stars, as well as to improve the modelling of the X-ray emission from the surface of rotating stars.

To unify the standard model of particle physics and general relativity, we may require a quantum description of gravity, which will change our notion of spacetime at very high energies. In this dissertation we explore possible traces of new physics beyond special relativity, using the propagation of high energy astroparticles. For this purpose, the two ways of going beyond Lorentz invariance are presented, a breaking of the Lorentz invariance (Lorentz invariance violation or LIV or its deformation (doubly special relativity or DSR), emphasizing their conceptual and phenomenological differences. For the study of LIV, the work focuses on the prediction of modifications in the expected neutrino flux on Earth, both from astrophysical and cosmogenic origin (from the interaction of cosmic rays with the background radiation during their propagation). For the study of DSR we focus instead on the search for anomalies in the time of flight of massless particles (time delays) and on the study of the expected flux of gamma rays on Earth. The results obtained show the possibility of using astroparticle observations as a window to quantum gravity phenomenology, at energies attainable at present and/or in the very near future.

New Early Dark Energy introduces a new phase of dark energy that decays in a fast-triggered phase transition around matter-radiation equality. The presence of a trigger mechanism sets it apart from other early dark energy models. Here, we will argue that New Early Dark Energy offers a simple and natural framework to extend $\Lambda$CDM while also providing a pathway to resolving the $H_0$ tension alongside its smaller cousin, the $S_8$ tension. At the microscopic level, we discuss the possibility that the trigger is either given by an ultralight scalar field or a dark sector temperature. In both cases, it prompts the transition of an $\mathrm{eV}$-mass scalar field from its false to its true minimum. Furthermore, we argue that the same phase transition could give rise to a dynamic process for generating neutrino masses.

Supernova neutrino boosted dark matter (SN$\nu$ BDM) and its afterglow effect have been shown to be a promising signature for beyond Standard Model (bSM) physics. The time-evolution feature of SN$\nu$ BDM allows for possibly direct inference of DM mass $m_\chi$, and results in significant background suppression with improving sensitivity. This paper extends the earlier study and provides a general framework for computing the SN$\nu$ BDM fluxes for a supernova that occurs at any location in our galaxy. A bSM $U(1)_{L_\mu-L_\tau}$ model with its gauge boson coupling to both DM and the second and third generation of leptons is considered, which allows for both DM-$\nu$ and DM-$e$ interactions. Detailed analysis of the temporal profile, angular distribution, and energy spectrum of the SN$\nu$ BDM are performed. Unique signatures in SN$\nu$ BDM allowing extraction of $m_\chi$ and detail features that contain information of the underlying interaction type are discussed. Expected sensitivities on the above new physics model from Super-Kamiokande, Hyper-Kamiokande, and DUNE detections of BDM events induced by the next galactic SN are derived and compared with the existing bounds.

We investigate neutrino non-radiative two-body decay in vacuum, in relation to SN1987A. In a full $3\nu$ decay framework, we perform a detailed likelihood analysis of the 24 neutrino events from SN1987A observed by Kamiokande-II, IMB, and Baksan. We consider both normal and inverted neutrino mass orderings, and the possibility of strongly hierarchical and quasi-degenerate neutrino mass patterns. The results of the likelihood analysis show that the sensitivity is too low to derive bounds in the case of normal mass ordering. On the contrary, in the case of inverted mass ordering we obtain the bound $\tau/m \ge 2.4 \times 10^{5}$ s/eV ($1.2 \times 10^{5}$) s/eV at 68 $\%$ (90 $\%$) CL on the lifetime-to-mass ratio of the mass eigenstates $\nu_2$ and $\nu_1$.

In this paper, we study a possible early universe source for the recent observation of a stochastic gravitational wave background at the NANOGrav pulsar timing array. The source is a tachyonic instability in a dark gauge field induced by an axion-like particle (ALP), a known source for gravitational waves. We find that relative to the previous analysis with the NANOGrav 12.5-year data set, the current 15-year data set favors parameter space with a relatively larger axion mass and decay constant. This favored parameter space is heavily constrained by $\Delta N_{\rm eff}$ and overproduction of ALP dark matter. While there are potential mechanisms for avoiding the second problem, evading the $\Delta N_{\rm eff}$ constraint remains highly challenging. In particular, we find that the gravitational wave magnitude is significantly suppressed with respect to the gauge boson dark radiation, which implies that successfully explaining the NANOGrav observation requires a large additional dark radiation, violating the cosmological constraints.

We implement an eccentric search for compact binary mergers based on particle swarm optimization. The orbital eccentricity is an invaluable input for understanding the formation scenarios of the binary mergers and can play a pivotal role in finding their electromagnetic counterparts. Current modelled searches rely on pre-computed template banks that are computationally expensive and resistant towards expanding the search parameter space dimensionality. On the other hand, particle swarm optimization offers a straightforward algorithm that dynamically selects template points while exploring an arbitrary dimensional parameter space. Through extensive evaluation using simulated signals from spin-aligned eccentric binary mergers, we discovered that the search exhibits a remarkable autonomy in capturing the effects of both eccentricity and spin. We describe our search pipeline and revisit some of the merger candidates from the gravitational wave transient catalogs.