Papers for Monday, Apr 29 2024

We present in this paper a compilation of solar atlases from $\lambda$ 3000 {\AA} to $\lambda$ 8800 {\AA} with spectral lines identified by the Moore table and with the corresponding equivalent Lande factors g*. We used two spectra at disk centre ($\mu$ = 1.0), from Delbouille and Kurucz, and two spectra at the limb from Stenflo and Gandorfer, respectively at $\mu$ = 0.145 and $\mu$ = 0.10.

Sigma Orionis is an open cluster in the nearest giant star formation site - Orion. Its youth (3-5 Myr), low reddening, and relative vicinity make it an important benchmark cluster to study stellar and substellar formation and evolution. Young star-forming sites are complex and hierarchical. Precision astrometry from Gaia DR3 enables the exploration of their fine structure. We used the modified convergent point technique to kinematically re-evaluate the members in the Sigma Orionis cluster and its vicinity. We present clear evidence for three kinematically distinct groups in the Sigma Orionis region. The second group, the RV Orionis association, is adjacent to the Sigma Orionis cluster and is composed only of low-mass stars. The third group, the Flame association, whose age is comparable to that of Sigma Orionis, overlaps with the younger NGC 2024 in the Flame Nebula. In total, we have discovered 105 members of this complex not previously found in the literature (82 in Sigma Orionis, 19 in the Flame association, and 4 in the RV Orionis association).

The presence of dense, neutral hydrogen clouds in the hot, diffuse intra-group and intra-cluster medium is an important clue to the physical processes controlling the survival of cold gas and sheds light on cosmological baryon flows in massive halos. Advances in numerical modeling and observational surveys means that theory and observational comparisons are now possible. In this paper, we use the high-resolution TNG50 cosmological simulation to study the HI distribution in seven halos with masses similar to the Fornax galaxy cluster. Adopting observational sensitivities similar to the MeerKAT Fornax Survey (MFS), an ongoing HI survey that will probe to column densities of $10^{18}$ cm$^{-2}$, we find that Fornax-like TNG50 halos have an extended distribution of neutral hydrogen clouds. Within one virial radius, we predict the MFS will observe a total HI covering fraction around $\sim$ 12\% (mean value) for 10 kpc pixels and 6\% for 2 kpc pixels. If we restrict this to gas more than 10 half-mass radii from galaxies, the mean values only decrease mildly, to 10\% (4\%) for 10 (2) kpc pixels (albeit with significant halo-to-halo spread). Although there are large amounts of HI outside of galaxies, the gas seems to be associated with satellites, judging both by the visual inspection of projections and by comparison of the line of sight velocities of galaxies and intracluster HI.

Nested sampling is a promising tool for Bayesian statistical analysis because it simultaneously performs parameter estimation and facilitates model comparison. MultiNest is one of the most popular nested sampling implementations, and has been applied to a wide variety of problems in the physical sciences. However, MultiNest results are frequently unreliable, and accompanying convergence tests are a necessary component of any analysis. Using simple, analytically tractable test problems, I illustrate how MultiNest (1) can produce systematically biased estimates of the Bayesian evidence, which are more significantly biased for problems of higher dimensionality; (2) can derive posterior estimates with errors on the order of $\sim100\%$; (3) is more likely to underestimate the width of a credible interval than to overestimate it - to a minor degree for smooth problems, but much more so when sampling noisy likelihoods. Nevertheless, I show how MultiNest can be used to jump-start Markov chain Monte Carlo sampling or more rigorous nested sampling techniques, potentially accelerating more trustworthy measurements of posterior distributions and Bayesian evidences, and overcoming the challenge of Markov chain Monte Carlo initialization.

The renormalization group equations for large-scale structure (RG-LSS) describe how the bias and stochastic (noise) parameters -- both of matter and biased tracers such as galaxies -- evolve as a function of the cutoff $\Lambda$ of the effective field theory. In previous work, we derived the RG-LSS equations for the bias parameters using the Wilson-Polchinski framework. Here, we extend these results to include stochastic contributions, corresponding to terms in the effective action that are higher order in the current $J$. We show that the RG equations exhibit an interesting, previously unnoticed structure at all orders in $J$, which implies that a single nonlinear bias term immediately generates all stochastic moments through RG evolution. We then derive the nonlinear RG evolution of the (leading-derivative) stochastic parameters for all $n$-point functions, and show that this evolution is controlled by a different, lower scale than the nonlinear scale. This has implications for the optimal choice of the renormalization scale when comparing the theory with data to obtain cosmological constraints.

Accurate modeling of selection effects is a key ingredient to the success of gravitational-wave astronomy. The detection probability plays a crucial role in both statistical population studies, where it enters the hierarchical Bayesian likelihood, and astrophysical modeling, where it is used to convert predictions from population-synthesis codes into observable distributions. We review the most commonly used approximations, extend them, and present some recipes for a straightforward implementation. These include a closed-form expression capturing both multiple detectors and noise realizations written in terms of the so-called Marcum Q-function and a ready-to-use mapping between signal-to-noise ratio thresholds and false-alarm rates from state-of-the-art detection pipelines. The bias introduced by approximating the matched filter signal-to-noise ratio with the optimal signal-to-noise ratio is not symmetric: sources that are nominally below threshold are more likely to be detected than sources above threshold are to be missed. Using both analytical considerations and software injections in detection pipelines, we confirm that including noise realizations when estimating the selection function introduces an average variation of a few %. This effect is most relevant for large catalogs and specific subpopulations of sources at the edge of detectability (e.g. high redshifts).

The spectral energy distribution (SED) of low-luminosity active galactic nuclei (LLAGN) presents unique challenges due to their comparable radiation output to their host galaxies and complex accretion dynamics. We introduce a novel module within the CIGALE framework specifically designed for SED fitting of LLAGN, incorporating both empirical relationships like $L_\mathrm{X}$--$L_\mathrm{12\mu m}$ and physically-based accretion models such as advection-dominated accretion flows (ADAFs) and truncated accretion disks. This allows for more accurate depiction of LLAGN central emissions. Using this module, we analyzed a set of 52 X-ray-detected local galaxies, primarily LINERs and Seyferts, and compared its performance to higher-luminosity AGN from the COSMOS and SDSS datasets. Our results show that the module adeptly estimates bolometric luminosities with high precision, despite significant galaxy contamination. It also introduces a versatile X-ray bolometric correction formula covering a vast range of luminosities. Further, our study explored the $\alpha_\mathrm{ox}$ index, which measures the UV to X-ray emission slope, showing that unlike quasars, LLAGN display either stable or only slightly varying $\alpha_\mathrm{ox}$ values, indicating differing accretion and photon production processes in the low luminosity regime. Additionally, we observed a significant drop of 1.4 dex in specific star formation rates when moving from whole galaxies to a central 9-arcsecond aperture in LLAGN, suggesting potential feedback mechanisms at play. Overall, our findings underscore the importance of a multiwavelength approach in AGN studies, highlighting distinct behaviors of LLAGN compared to quasars, thus enhancing our understanding of LLAGN and providing a framework for future comprehensive AGN population studies.

Stellar halos of galaxies retain crucial clues to their mass assembly history. It is in these galactic components that the remains of cannibalised galactic building blocks are deposited. For the case of the Milky Way, the opportunity to analyse the stellar halo's structure on a star-by-star basis in a multi-faceted approach provides a basis from which to infer its past and assembly history in unrivalled detail. Moreover, the insights that can be gained about the formation of the Galaxy not only help constrain the evolution of our Milky Way, but may also help place constraints on the formation of other disc galaxies in the Universe. This paper includes a summary of work undertaken during a PhD thesis aiming to make progress toward answering the most fundamental question in the field of Galactic archaeology: "How did the Milky Way form?" Through the effort to answer this question, we summarise new insights into aspects of the history of assembly and evolution of our Galaxy and measurements of the structure of various of its Galactic components.

We describe a method for calculating action-angle variables in axisymmetric galactic potentials using Birkhoff normalization, a technique from Hamiltonian perturbation theory. An advantageous feature of this method is that it yields explicit series expressions for both the forward and inverse transformations between the action-angle variables and position-velocity data. It also provides explicit expressions for the Hamiltonian and dynamical frequencies as functions of the action variables. We test this method by examining orbits in a Miyamoto-Nagai model potential and compare it to the popular St\"ackel approximation method. When vertical actions are not too large, the Birkhoff normalization method achieves fractional errors smaller than a part in $10^{3}$ and outperforms the St\"ackel approximation. We also show that the range over which Birkhoff normalization provides accurate results can be extended by constructing Pad\'e approximants from the perturbative series expressions developed with the method. Numerical routines in Python for carrying out the Birkhoff normalization procedure are made available.

An analysis of calibrations of extensive air showers with zenith angles $\theta \le 50^{\circ}$ and energies $E_{\text{SD}} \le 10^{18.5}$ eV was carried out in experiments at the Yakutsk array and Telescope Array. The values of $E_{\text{SD}}$ were determined from particle densities measured with ground-based scintillation detectors at a distance $r = 800$~m from shower axis. Measured densities were compared with the values obtained in simulations preformed with the use of CORSIKA code within the framework of QGSJet-II.04 hadron interaction model for primary protons. For showers with $\theta = 0^{\circ}$ the $E_{\text{SD}}$ estimates of both arrays are very close, and their energy spectra are similar in shape and absolute value. Another Telescope Array calibration, based on measuring the EAS fluorescent radiation with optical detectors, gives underestimated energy $E_{\text{FD}} = 0.787 \times E_{\text{SD}}$.

Strong gravitational lensing offers a powerful probe of the detailed distribution of matter in lenses, while magnifying and bringing faint background sources into view. Observed strong lensing by massive galaxy clusters, which are often in complex dynamical states, has also been used to map their dark matter substructures on smaller scales. Deep high resolution imaging has revealed the presence of strong lensing events associated with these substructures, namely galaxy-scale sub-halos. However, an inventory of these observed galaxy-galaxy strong lensing (GGSL) events is noted to be discrepant with state-of-the-art {\Lambda}CDM simulations. Cluster sub-halos appear to be over-concentrated compared to their simulated counterparts yielding an order of magnitude higher value of GGSL. In this paper, we explore the possibility of resolving this observed discrepancy by redistributing the mass within observed cluster sub-halos in ways that are consistent within the {\Lambda}CDM paradigm of structure formation. Lensing mass reconstructions from data provide constraints on the mass enclosed within apertures and are agnostic to the detailed mass profile within them. Therefore, as the detailed density profile within cluster sub-halos currently remains unconstrained by data, we are afforded the freedom to redistribute the enclosed mass. We investigate if rearranging the mass to a more centrally concentrated density profile helps alleviate the GGSL discrepancy. We report that refitting cluster sub-halos to the ubiquitous {\Lambda}CDM-motivated Navarro-Frenk-White profile, and further modifying them to include significant baryonic components, does not resolve this tension. A resolution to this persisting GGSL discrepancy may require more careful exploration of alternative dark matter models.

We investigate the dynamics of GN-z11, a luminous galaxy at $z=10.60$, carefully analyzing the public deep integral field spectroscopy (IFS) data taken with JWST NIRSpec IFU. While the observations of the IFS data originally targeted a He II clump near GN-z11, we find that C III]$\lambda\lambda1907,1909$ emission from ionized gas at GN-z11 is bright and spatially extended significantly beyond the point-spread function (PSF) in contrast with the GN-z11's compact UV-continuum morphology. The spatially extended C III] emission of GN-z11 shows a velocity gradient, red- and blue-shifted components in the NE and SW directions, respectively, which cannot be explained by the variation of [C III]$\lambda$1907/C III]$\lambda$1909 line ratios. We perform forward modeling with GalPak$^\mathrm{3D}$, including the effects of PSF smearing and line blending, and find that the best-fit model is a nearly edge-on disk with a rotation velocity of $v_\mathrm{rot}$=$376_{-151}^{+110}$ km s$^{-1}$, a velocity dispersion of $\sigma_v=113_{-48}^{+24}$ km s$^{-1}$, and a ratio of $v_\mathrm{rot}/\sigma_v=3.3_{-1.5}^{+1.7}$, indicative of a rotation-dominated disk at $z=10.6$. Interestingly, the disk rotation velocity is faster than the circular velocity at the virial radius, $v_\mathrm{c}(r_{200})=205_{-34}^{+28}$ km s$^{-1}$, estimated from the stellar mass, suggesting a compact disk produced under weak feedback such predicted in numerical simulations. In fact, the half-light radius of the C III] emitting gas disk is only $294\pm45$ pc, while the one of the stellar UV component is even smaller, $196\pm12$ pc. While higher S/N data are necessary for a conclusion, these observational results would suggest that GN-z11 has a fast-rotating gaseous disk whose center possesses luminous stellar components or AGN providing weak feedback.

LP~349$-$25 is a well studied close stellar binary system comprised of two late M dwarf stars, both stars close to the limit between star and brown dwarf. This system was previously identified as a source of GHz radio emission. We observed LP~349$-$25AB in 11 epochs in 2020$-$2022, detecting both components in this nearby binary system using the Very Long Baseline Array (VLBA). We fit simultaneously the VLBA absolute astrometric positions together with existing relative astrometric observations derived from optical/infrared observations with a set of algorithms that use non-linear least-square, Genetic Algorithm and Markov Chain Montecarlo methods to determine the orbital parameters of the two components. We find the masses of the primary and secondary components to be 0.08188 $\pm$ 0.00061 $M\odot$, and 0.06411 $\pm$ 0.00049 $M_\odot$, respectively, representing one of the most precise mass estimates of any UCDs to date. The primary is an UCD of 85.71$\pm$0.64 M$_{Jup}$, while the secondary has a mass consistent with being a Brown Dwarf of 67.11 $\pm$ 0.51 $M_{Jup}$. This is one of the very few direct detections of a Brown Dwarf with VLBA observations. We also find a distance to the binary system of 14.122 $\pm$ 0.057 pc. Using Stellar Evolutionary Models, we find the model-derived stellar parameters of both stars. In particular, we obtain a model-derived age of 262 Myr for the system, which indicates that LP~349$-$25AB is composed of two pre-main-sequence stars. In addition, we find that the secondary star is significantly less evolved that the primary star.

High-resolution X-ray spectra of $\pi\,$Aqr, a $\gamma\,$Cas-type star, obtained with the Chandra/HETG grating spectrometer, revealed emission lines of H-like ions of Mg, Si, S, and Fe, a strong, hard continuum, and a lack of He-like ions, indicating the presence of very hot thermal plasma. The X-ray light curve showed significant fluctuations, with coherent variability at period of about 3400 seconds in one observation. The hardness ratio was relatively constant except for one observation in which the spectrum was much harder and more absorbed. We interpret the X-ray emission as arising from accretion onto the secondary, which is likely a magnetic white dwarf, an intermediate polar system.

This paper presents an overview of the high-redshift extragalactic science case for the Single Aperture Large Telescope for Universe Studies (SALTUS) far-infrared NASA probe-class mission concept. Enabled by its 14m primary reflector, SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions. SALTUS would be a versatile observatory capable of responding to the scientific needs of the extragalactic community in the 2030s, and a natural follow-on to the near- and mid-IR capabilities of JWST. Key early-universe science goals for SALTUS focus on understanding the role of galactic feedback processes in regulating galaxy growth across cosmic time, and charting the rise of metals and dust from the early universe to the present. This paper summarizes these science cases and the performance metrics most relevant for high-redshift observations.

This paper presents an overview of the Milky Way and nearby galaxies science case for the \textit{Single Aperture Large Telescope for Universe Studies} (SALTUS) far-infrared NASA probe-class mission concept. SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions, thanks to its cold ($<$40~K) 14-m primary mirror. Key Milky Way and nearby galaxies science goals for SALTUS focus on understanding the role of star formation in feedback in the Local Universe. In addition to this science case, SALTUS would open a new window to to of Galactic and extragalactic communities in the 2030s, enable fundamentally new questions to be answered, and be a far-IR analog to the near- and mid-IR capabilities of JWST. This paper summarizes the Milky Way and nearby galaxies science case and plans for notional observing programs in both guaranteed and guest (open) time.

In galaxy formation simulations, it is increasingly common to represent supernovae (SNe) at finite resolution (when the Sedov-Taylor phase is unresolved) via hybrid energy-momentum coupling with some 'terminal momentum' $p_{\rm term}$ (depending weakly on ambient density and metallicity) that represents unresolved work from an energy-conserving phase. Numerical implementations can ensure momentum and energy conservation of such methods, but these require that couplings depend on the surrounding gas velocity field (radial velocity $\langle v_{r} \rangle$). This raises the question of whether $p_{\rm term}$ should also be velocity-dependent, which we explore analytically and in simulations. We show that for simple spherical models, the dependence of $p_{\rm term}$ on $\langle v_{r} \rangle$ introduces negligible corrections beyond those already imposed by energy conservation if $\langle v_{r} \rangle \ge 0$. However, for SNe in some net converging flow ($\langle v_{r} \rangle<0$), naively coupling the total momentum when a blastwave reaches the standard cooling/snowplow phase (or some effective cooling time/velocity/temperature criterion) leads to enormous $p_{\rm term}$ and potentially pathological behaviors. We propose an alternative $\Delta$-Momentum formulation which represents the differential SNe effect and show this leads to opposite behavior of $p_{\rm term}$ in this limit. We also consider a more conservative velocity-independent formulation. Testing in numerical simulations, these directly translate to large effects on predicted star formation histories and stellar masses of massive galaxies, explaining differences between some models and motivating further study in idealized simulations.

There are different classes of pulsating stars in the H-R diagram. While many of those classes are undisputed, some remain a mystery such as the objects historically called "Maia variables." Whereas the presence of such a class was suggested seven decades ago, no pulsational driving mechanism is known that could excite short-period oscillations in these late B to early A-type stars. Alternative hypotheses that would render the reports of variability of those stars erroneous have been proposed, such as incorrect effective temperatures, binarity or rapid rotation, but no certain conclusions have been reached yet. Therefore, the existence of these variables as a homogeneous class of pulsating star is still under discussion. Meanwhile, many new candidates of these variables have been claimed especially by using photometric observations from space telescopes. In this study, we examined 31 objects that are alleged members of this hypothetical group and carried out detailed spectroscopic and photometric analyses to test the proposed hypotheses for their cause of variability. The $T_{eff}$, log$(g)$, $v$ sin $i$, and chemical abundances of the targets were determined and the $TESS$ photometric data were examined. As a result, we found that most of these targets are located inside the $\delta$ Scuti, $\beta$ Cephei, or SPB star instability strips, a few show evidence for binarity, and others exhibit rapid rotation. We give arguments that none of the apparently rapid pulsations in our targets is caused by a star outside any known instability strip. By extrapolation, we argue that most stars proposed as pulsators outside well-established instability domains are misclassified. Hence there is no sufficient evidence justifying the existence of a class of pulsating stars formerly known as the "Maia variables."

We present an expanded and improved deep-learning (DL) methodology for determining centers of star images on HST/WFPC2 exposures. Previously, we demonstrated that our DL model can eliminate the pixel-phase bias otherwise present in these undersampled images; however that analysis was limited to the central portion of each detector. In the current work we introduce the inclusion of global positions to account for the PSF variation across the entire chip and instrumental magnitudes to account for nonlinear effects such as charge transfer efficiency. The DL model is trained using a unique series of WFPC2 observations of globular cluster 47 Tuc, data sets comprising over 600 dithered exposures taken in each of two filters, F555W and F814W. It is found that the PSF variations across each chip correspond to corrections of the order of 100 mpix, while magnitude effects are at a level of about 10 mpix. Importantly, pixel-phase bias is eliminated with the DL model; whereas, with a classic centering algorithm, the amplitude of this bias can be up to 40 mpix. Our improved DL model yields star-image centers with uncertainties of 8-10 mpix across the full field of view of WFPC2.

As part of an intensive effort to observe the predicted 2022 Tau Herculids outburst, we recorded almost 800 individual meteor streaks on May 30th and 31st, 2022, using a high-sensitivity Sony~$\alpha$7 camera. The video recordings were obtained under perfect conditions at the McDonald Observatory, Texas, USA. The meteor sample is dominated by the predicted Tau Herculids shower, however, we also noted significant activity of sporadic meteors and other possible weak showers. We found that the time of the maximum activity matched very well the predictions, while we note the large fraction of faint meteors that were not detectable visually. We determined the radiant, and the time evolution of the activities and currently we are working on the determination of the brightness statistics.

Here we present a continuation of the work that was based on video observations of the Tau Her 2022 outburst from the McDonald Observatory, Texas, US. On the night of the maximum in 2022 we detected 626 individual Tau Her meteors, for which we estimated photovisual magnitudes and analysed their distribution on the sky to determine the radiant position in an innovative way. The derived population index is $5.56 \pm 1.83$, while the radiant position for the mid-time of the observations is: $\mathrm{RA} = 209.71 ^\circ \pm 1.01 ^\circ$, $\mathrm{Dec} = 27.73 ^\circ \pm 0.07 ^\circ$. Both measurements are in good agreement with other results in the literature, although our population index seems to be higher than most of the published values. We speculate this might be due to the higher sensitivity of our equipment to fainter meteors.

We present a study of spectral line width measurements from the Extreme Ultraviolet Imaging Spectrometer (EIS) on {\it Hinode}. We used spectral line profiles of Fe {\sc xvi} 262.984 {\AA}, Fe {\sc xiv} 264.787 {\AA}, Fe {\sc xiv} 270.519 {\AA}, Fe {\sc xiv} 274.203 {\AA}, and Fe {\sc xv} 284.160 {\AA}, and studied 11 active regions. Previous studies of spectral line widths have shown that in hot loops in the cores of active regions, the observed non-thermal velocities are smaller than predicted from models of reconnection jets in the corona or shock heating associated with Alfv\'{e}n waves. The observed line widths are also inconsistent with models of chromospheric evaporation due to coronal nanoflares. We show that recent advances in higher resolution Alfv\'{e}n wave turbulence modeling enables us to obtain non-thermal velocities similar to those measured in active regions. The observed non-thermal velocities for the 11 active regions in our study are in the range of 17$-$30 $\rm km ~ s^{-1}$, consistent with the spectral line non-thermal widths predicted from our model of 16 interacting flux tubes, which are in the range of ~15$-$37 $\rm km ~ s^{-1}$.

Interest in the study of magnetic fields and the properties of interstellar dust, explored through increasingly capable far-IR/submillimeter polarimetry, along with maturing detector technology, have set the stage for a transformative leap in polarization mapping capability using a cryogenic space telescope. We describe the approach pursued by the proposed Probe far-Infrared Mission for Astrophysics (PRIMA) to make ultra-deep maps of intensity and polarization in four bands in the 91-232 micron range. A simple, polarimetry-optimized PRIMA Polarimetric Imager (PPI) is designed for this purpose, consisting of arrays of single-polarization Kinetic Inductance Detectors oriented with three angles which allow measurement of Stokes I, Q, and U in single scans. In this study, we develop an end-to-end observation simulator to perform a realistic test of the approach for the case of mapping a nearby galaxy. The observations take advantage of a beam-steering mirror to perform efficient, two-dimensional, crossing scans. Map making is based on 'destriping' approaches demonstrated for Herschel/SPIRE and Planck. Taking worst-case assumptions for detector sensitivity including 1/f noise, we find excellent recovery of simulated input astrophysical maps, with I, Q, and U detected at near fundamental limits. We describe how PPI performs detector relative calibration and mitigates the key systematic effects to accomplish PRIMA polarization science goals.

Large-scale coronal plasma evolutions can be adequately described by magnetohydrodynamics (MHD) equations. However, full multi-dimensional MHD simulations require substantial computational resources. Given the low plasma $\beta$ in the solar corona, in many coronal studies, it suffices to approximate the magnetic field to remain topologically fixed and effectively conduct one-dimensional (1D) hydrodynamic (HD) simulations instead. This approach is often employed in studies of coronal loops and their liability to form condensations related to thermal instability. While 1D HD simulations along given and fixed field line shapes are convenient and fast, they are difficult to directly compare with multi-dimensional phenomena. Therefore, it is more convenient to solve volume-filling, multi-dimensional versions of the MHD equations where we freeze the magnetic field, transforming it into frozen-field HD (ffHD) equations for simulation. We have incorporated this ffHD module into our open-source MPI-AMRVAC code and tested it using a two-dimensional (2D) evaporation--condensation model to study prominence formation due to radiative losses. The 2D ffHD results are compared with those from actual 2D MHD and pseudo-2D HD simulations, analyzing the differences and their causes. Pseudo-2D studies account for the known {flux tube} expansion effects. Overall, the performance of 2D ffHD is close to that of 2D MHD and pseudo-2D HD. The 2D tests conducted in this paper will be extended in follow-up studies to 3D simulations based on analytical or observational approaches.

The discovery of a star formed out of pair-instability supernova ejecta would have massive implications for the Population III star initial mass function and the existence of stars over 100 Msun, but none have yet been found. Recently, the star LAMOST J1010+2358 was claimed to be a star that formed out of gas enriched by a pair-instability supernova. We present a non-LTE abundance analysis of a new high-resolution Keck/HIRES spectrum of J1010+2358. We determined the carbon and aluminum abundances needed to definitively distinguish between enrichment by a pair-instability and core-collapse supernova. Our new analysis demonstrates that J1010+2358 does not have the unique abundance pattern of a a pair-instability supernova, but was instead enriched by the ejecta of a low mass core-collapse supernova. Thus, there are still no known stars displaying unambiguous signatures of pair-instability supernovae.

Modelling the interaction between ionizing photons emitted from massive stars and their environment is essential to further our understanding of galactic ecosystems. We present a hybrid Radiation-Hydrodynamics (RHD) scheme that couples an SPH code to a grid-based Monte Carlo Radiative Transfer code. The coupling is achieved by using the particle positions as generating sites for a Voronoi grid, and applying a precise mapping of particle-interpolated densities onto the grid cells that ensures mass conservation. The mapping, however, can be computationally infeasible for large numbers of particles. We introduce our tree-based algorithm for optimizing coupled RHD codes. Astrophysical SPH codes typically utilize tree-building procedures to sort particles into hierarchical groups (referred to as nodes) for evaluating self-gravity. Our algorithm adaptively walks the gravity tree and transforms the extracted nodes into pseudo-SPH particles, which we use for the grid construction and mapping. This method allows for the temporary reduction of fluid resolution in regions that are less affected by the radiation. A neighbour-finding scheme is implemented to aid our smoothing length solver for nodes. We show that the use of pseudo-particles produces equally accurate results that agree with benchmarks, and achieves a speed-up that scales with the reduction in the final number of particle-cell pairs being mapped.

Within the mass range of $10^{16}-5\times 10^{18}$ g, primordial black holes (PBHs) persist as plausible candidates for dark matter. Our study involves a reassessment of the constraints on PBHs through a comparative analysis of the cosmic X-ray background (CXB) and the emissions arising from their Hawking evaporation. We identify previously overlooked radiation processes across the relevant energy bands, potentially refining the bounds on PBHs. These processes encompass the direct emission from Hawking radiation, in-flight annihilation, the final state of radiation, and positronium annihilation. Thorough consideration is given to all these processes and their respective emission fractions, followed by a precise calculation of the $\mathcal{D}$ factor for observations directed towards the high-latitude Galactic contribution. Furthermore, we integrate the flux originating from extragalactic sources, both of which contribute to the measured isotropic flux. Through a comparative analysis of data derived from previous CXB observations utilizing an Active Galactic Nuclei (AGN) double power-law model, we establish the most stringent constraints for PBHs, thereby excluding the possibility of PBHs constituting the entire dark matter mass within the range of $2.5\times 10^{17}-3 \times 10^{17} $g.

In a strong gravitational lens, perturbations by low-mass dark matter halos can be detected by differences between the measured image fluxes relative to the expectation from a smooth model for the mass distribution which contains only the gravitational effects of the main deflector. The abundance of these low-mass structures can be used to constrain the properties of dark matter. Traditionally only the lensed quasar positions have been to predict the smooth-model flux ratios. We demonstrate that significant additional information can be gained by using the lensed quasar host galaxy which appears as an extended arc and constrains the smooth-model over a much larger angular area. We simulate Hubble Space Telescope-quality mock observations based on the lensing system WGD2038-4008 and we compare the model-predicted flux ratio precision and accuracy for two cases; one of which the inference is based only on the lensed quasar image positions, and the other based on the extended arcs as well as lensed quasar image positions. For our mock lens systems we include both elliptical, and higher order m=3 and m=4 multipole terms in the smooth-mass distributions with amplitudes based on the optically measured shapes of massive elliptical galaxies. We find that the extended arcs improve the precision of the model-predicted flux ratios by a factor of 6-8, depending on the strength of the multipole terms. Furthermore, with the extended arcs, we are also able to accurately recover the m=3, 4 mass multipole strengths and angles $a_3/a$, $a_4/a$, $\phi_3-\phi_0$, and $\phi_4-\phi_0$ to a precision of 0.002, 0.002, $3^\circ$ and $3^\circ$, respectively. This work implies that lensed arcs can constrain deviations from ellipticity in strong lens systems, and potentially lead to more robust constraints on substructure properties from flux ratios.

The theoretical debris supply rate from a tidal disruption of stars can exceed about one hundred times of the Eddington accretion rate for a $10^{6-7}M_{\odot}$ supermassive black hole (SMBH). It is believed that a strong wind will be launched from the disk surface due to the radiation pressure in the case of super-Eddington accretion, which may be one of the mechanisms for formation of the envelope as observed in tidal disruption events (TDEs). In this work, we explore the evolution of the envelope that formed from the optical thick winds by solving the global solution of the slim-disk model. Our model can roughly reproduce the typical temperature, luminosity and size of the photosphere for TDEs. Based on \texttt{CLOUDY} modeling, we find that, if only considering the radiation-driven disk wind, the emission line luminosities are normally much lower than the typical observational results due to the limited atmosphere mass outside the envelope. We propose that the ejection of the outflow from the self-collision of the stellar debris during the circularization may provide enough matter outside the disk-wind photosphere. Our calculated spectra can roughly reproduce the main properties of several typical emission lines (e.g., $\rm H\alpha$, $\rm H\beta$ and \ion{He}{2}), which was applied well to a TDE candidate AT2018dyb.

The effect of metallicity on the theoretical and empirical period-luminosity (PL) relations of Cepheid variables is not well understood and remains a highly debated issue. Here, we examine empirical colour-magnitude diagrams (CMDs) of Classical and Type-II Cepheids in the Magellanic Clouds and compare those with the theoretically predicted instability strip (IS) edges. We explore the effects of incorporating turbulent flux, turbulent pressure, and radiative cooling into the convection theory on the predicted IS at various metallicities using MESA-RSP. We find that the edges become redder with the increasing complexity of convection physics incorporated in the fiducial convection sets, and are similarly shifted to the red with increasing metallicity. The inclusion of turbulent flux and pressure improves the agreement of the red edge of the IS, while their exclusion leads to better agreement with observations of the blue edge. About 90% of observed stars are found to fall within the predicted bluest and reddest edges across the considered variations of turbulent convection parameters. Furthermore, we identify and discuss discrepancies between theoretical and observed CMDs in the low effective temperature and high luminosity regions for stars with periods greater than ~ 20 days. These findings highlight the potential for calibrating the turbulent convection parameters in stellar pulsation models or the prediction of a new class of rare, long-period, 'red Cepheids', thereby improving our understanding of Cepheids and their role in cosmological studies.

A hypothetical photon mass, $m_{\gamma}$, can produce a frequency-dependent vacuum dispersion of light, which leads to an additional time delay between photons with different frequencies when they propagate through a fixed distance. The dispersion measure--redshift measurements of fast radio bursts (FRBs) have been widely used to constrain the rest mass of the photon. However, all current studies analyzed the effect of the frequency-dependent dispersion for massive photons in the standard $\Lambda$CDM cosmological context. In order to alleviate the circularity problem induced by the presumption of a specific cosmological model based on the fundamental postulate of the masslessness of photons, here we employ a new model-independent smoothing technique, Artificial Neural Network (ANN), to reconstruct the Hubble parameter $H(z)$ function from 34 cosmic-chronometer measurements. By combining observations of 32 well-localized FRBs and the $H(z)$ function reconstructed by ANN, we obtain an upper limit of $m_{\gamma} \le 3.5 \times 10^{-51}\;\rm{kg}$, or equivalently $m_{\gamma} \le 2.0 \times 10^{-15}\;\rm{eV/c^2}$ ($m_{\gamma} \le 6.5 \times 10^{-51}\;\rm{kg}$, or equivalently $m_{\gamma} \le 3.6 \times 10^{-15}\;\rm{eV/c^2}$) at the $1\sigma$ ($2\sigma$) confidence level. This is the first cosmology-independent photon mass limit derived from extragalactic sources.

We report the first detection of evolving Low-Frequency Quasi-periodic Oscillation (LFQPO) frequencies in hard X-rays upto $100$ keV with $AstroSat/LAXPC$ during `unusual' outburst phase of Swift J1727.8-1613 in hard-intermediate state (HIMS). The observed LFQPO in $20 - 100$ keV has a centroid $\nu_{_{\rm QPO}}=1.43$ Hz, a coherence factor $Q= 7.14$ and an amplitude ${\rm rms_{_{\rm QPO}}} = 10.95\%$ with significance $\sigma = 5.46$. Type-C QPOs ($1.09-2.6$ Hz) are found to evolve monotonically during HIMS of the outburst with clear detection in hard X-rays ($80 - 100$ keV), where ${\rm rms_{_{\rm QPO}}}$ decreases ($\sim 12-3\%$) with energy. Further, $\nu_{_{\rm QPO}}$ is seen to correlate (anti-correlate) with low (high) energy flux in $2-20$ keV ($15-50$ keV). Wide-band ($0.7 - 40$ keV) energy spectrum of $NICER/XTI$ and $AstroSat/LAXPC$ is satisfactorily described by the `dominant' thermal Comptonization contribution ($\sim 88\%$) in presence of a `weak' signature of disk emissions ($kT_{\rm in} \sim 0.36$ keV) indicating the harder spectral distribution. Considering source mass $M_{\rm BH}=10M_\odot$ and distance $1.5 < {\rm d~(kpc)} < 5$, the unabsorbed bolometric luminosity is estimated as $\sim 0.03-0.92\%L_{\rm Edd}$. Finally, we discuss the implications of our findings in the context of accretion dynamics around black hole X-ray binaries.

Magnetic fields in the stellar interiors are key candidates to explain observed core rotation rates inside solar-like stars along their evolution. Recently, asteroseismic estimates of radial magnetic field amplitudes near the hydrogen-burning shell (H-shell) inside about 24 red-giants (RGs) have been obtained by measuring frequency splittings from their power spectra. Using general Lorentz-stress (magnetic) kernels, we investigated the potential for detectability of near-surface magnetism in a 1.3 $M_{\odot}$ star of super-solar metallicity as it evolves from a mid sub-giant to a late sub-giant into an RG. Based on these sensitivity kernels, we decompose an RG into three zones - deep core, H-shell, and near-surface. The sub-giants instead required decomposition into an inner core, an outer core, and a near-surface layer. Additionally, we find that for a low-frequency g-dominated dipolar mode in the presence of a typical stable magnetic field, ~25% of the frequency shift comes from the H-shell and the remaining from deeper layers. The ratio of the subsurface tangential field to the radial field in H-burning shell decides if subsurface fields may be potentially detectable. For p-dominated dipole modes close to $\nu_\rm{max}$, this ratio is around two orders of magnitude smaller in subgiant phases than the corresponding RG. Further, with the availability of magnetic kernels, we propose lower limits of field strengths in crucial layers in our stellar model during its evolutionary phases. The theoretical prescription outlined here provides the first formal way to devise inverse problems for stellar magnetism and can be seamlessly employed for slow rotators.

The process leading to the formation of the terrestrial planet remains elusive. In a previous publication, we have shown that, if the first generation of planetesimals forms in a ring at about 1 AU and the gas disk's density peaks at the ring location, planetary embryos of a few martian masses can grow and remain in the ring. In this work, we extend our simulations beyond the gas-disk stage, covering 200 Myr and accounting for the phase of giant planet instability, assumed to happen at different times. About half of the simulations form a pair of Venus and Earth analogues and, independently, about 10% form a Mars analogue. We find that the timing of the giant planet instability affects statistically the terrestrial system's excitation state and the timing of the last giant impacts. Hence a late instability (about 60 to 100 Myr after the Solar system's birth) is more consistent with a late Moon-formation time, as suggested by radioactive chronometers. However, the late veneer mass (LVM: mass accreted after the last giant impact) of Earth-sized planets suffering a giant impact after 80 My is usually an order of magnitude lower than the value inferred from geochemistry. In addition, the final angular momentum deficit (AMD) of the terrestrial planets tends to be too high. We tested the effect on the final AMD of the generation of debris during collisions and found that it is too small to change these conclusions. We argue that the best-case scenario is that the Moon-forming event occurred between 50 and 80 My, possibly just following the giant planet instability.

High precision spectrographs might exhibit temporal variations of their reference velocity or nightly zero point (NZP). One way to monitor the NZP is to measure bright stars, which are assumed to have an intrinsic radial velocity variation much smaller than the instrument's precision. While this method is effective in most cases, it does not fully propagate the uncertainty arising from NZP variations. We present a new method to correct for NZP variations in radial-velocity time series. This method uses Gaussian Processes based on ancillary information to model these systematic effects. It enables us to propagate the uncertainties of this correction into the overall error budget. Another advantage of this approach is that it relies on ancillary data collected simultaneously with the spectra rather than solely on dedicated observations of constant stars. We applied this method to the SOPHIE spectrograph at the Haute-Provence Observatory using a few instrument's housekeeping data, such as the internal pressure and temperature variations. Our results demonstrate that this method effectively models the red noise of constant stars, even with a limited amount of housekeeping data, while preserving the signals of exoplanets. Using both simulations with mock planets and real data, we found that this method improves the false-alarm probability of detections by several orders of magnitude. By simulating numerous planetary signals, we were able to detect up to 10 percent more planets with small amplitude radial velocity signals. We used this new correction to reanalysed the planetary system around HD158259 and improved the detection of the outermost planets. We also suggest decreasing the observing cadence of the constant stars to optimise telescope time for scientific targets.

While Tycho's supernova remnant is one of the most studied type Ia Galactic supernova remnants, a global view of the physical properties of its ejecta is lacking, to understand its mysteries. In particular, the spatial distribution of the Si-rich ejecta line-of-sight velocity presents a large-scale unexplained asymmetry, with the north dominantly blueshifted and the south redshifted. To investigate the origin of this line-of-sight velocity asymmetry in the ejecta, we carry out a detailed X-ray spatially-resolved spectral analysis of the entire shocked ejecta in Tycho's SNR to determine the physical properties of its various components. This study is based on the archival deep X-ray observations from the Chandra space telescope. The spatially-resolved spectral analysis in 211 regions over the entire SNR is based on a tesselation method applied to the line-of-sight velocity map. A Bayesian tool is used to conduct the fitting, using a nested sampling algorithm. It allows us to obtain a complete view of the statistical landscape. We provide maps of the physical parameters of the various components across the SNR ejecta. The Doppler shift map confirms spectrally the large-scale north-south asymmetry in the line-of-sight velocity. We reveal different spatial distributions of temperature and ionization time for IMEs and for iron-rich ejecta, but none of these maps shows structure associated to the large-scale north-south asymmetry in the line-of-sight velocity distribution. The abundance maps show spatial variations, depending on the element, perhaps due to an origin in different layers during the explosion. We compare these abundances with some nucleosynthesis models. In addition, we observe for the first time an emission line at 0.654 keV possibly related to oxygen. Its spatial distribution differs from the other elements, so that this line may arise in the ambient medium.

Classical Cepheids (DCEPs) are crucial for calibrating the extragalactic distance ladder, ultimately enabling the determination of the Hubble constant through the PL and PW relations they exhibit. Hence it's vital to understand how the PL and PW relations depend on metallicity. This is the purpose of the C-MetaLL survey within which this work is situated. DCEPs are also very important tracers of the young populations placed along the Galactic disc. We aim to enlarge the sample of DCEPs with accurate abundances from high-resolution spectroscopy. Our goal is to extend the range of measured metallicities towards the metal-poor regime to better cover the parameter space. We observed objects in a wide range of Galactocentric radii, allowing us to study in detail the abundance gradients present in the Galactic disc. We present the results of the analysis of 331 spectra obtained for 180 individual DCEPs with a variety of high-resolution spectrographs. We derived accurate atmospheric parameters, radial velocities, and abundances for up to 29 different species. The iron abundances range between 0.5 and -1 dex with a rather homogeneous distribution in metallicity. The sample presented in this paper was complemented with that already published in the context of the C-MetaLL survey, resulting in a total of 292 pulsators whose spectra have been analysed in a homogeneous way. These data were used to study the abundance gradients of the Galactic disc in a range of Galactocentric radii spanning the range 5-20 kpc. For most of the elements we found a clear negative gradient, with a slope of -0.071\pm0.003 dex kpc^-1 for [Fe/H] case. Through a qualitative fit with the Galactic spiral arms we shown how our farthest targets (R_GC>10 kpc) trace both the Outer and OSC arms. The homogeneity of the sample will be of pivotal importance for the study of the metallicity dependance of the DCEP PL relations.

In this paper, a kinetic energy recovery system for large telescopes is presented, with the Atacama Large Aperture Submm Telescope (AtLAST) as a possible target application. The system consists of supercapacitors integrated in the DC-link of motor inverters through a bidirectional DC-DC converter. The optimal system design, based on the energy flow analysis within the telescope's power electronics, is introduced. The proposed system is simulated as part of the telescope's drives, providing not only a significant reduction in energy consumption of the telescope due to motion, but also remarkably reducing (or shaving) grid power peaks. We find that the system presented here can contribute to making both current and future observatories more sustainable.

We report direct observations of a fast magnetosonic forward--reverse shock pair observed by Solar Orbiter on March 8, 2022 at the short heliocentric distance of 0.5 au. The structure, sharing some features with fully-formed stream interaction regions (SIRs), is due to the interaction between two successive coronal mass ejections (CMEs), never previously observed to give rise to a forward--reverse shock pair. The scenario is supported by remote observations from the STEREO-A coronographs, where two candidate eruptions compatible with the in-situ signatures have been found. In the interaction region, we find enhanced energetic particle activity, strong non-radial flow deflections and evidence of magnetic reconnection. At 1~au, well radially-aligned \textit{Wind} observations reveal a complex event, with characteristic observational signatures of both SIR and CME--CME interaction, thus demonstrating the importance of investigating the complex dynamics governing solar eruptive phenomena.

Promptly after high-resolution experiments harbinger the field of precision cosmology low- and high-redshift observations abruptly gave rise to a tension in the measurement of the present-day expansion rate of the Universe ($H_0$) and the clustering of matter ($S_8$). The statistically significant discrepancies between the locally measured values of $H_0$ and $S_8$ and the ones inferred from observations of the cosmic microwave background assuming the canonical $\Lambda$ cold dark matter (CDM) cosmological model have become a new cornerstone of theoretical physics. $\Lambda_s$CDM is one of the many beyond Standard Model setups that have been proposed to simultaneously resolve the cosmological tensions. This setup relies on an empirical conjecture, which postulates that $\Lambda$ switched sign (from negative to positive) at a critical redshift $z_c \sim 2$. We reexamine a stringy model that can describe the transition in the vacuum energy hypothesized in $\Lambda_s$CDM. The model makes use of the Casimir forces driven by fields inhabiting the incredible bulk of the dark dimension scenario. Unlike the $\Lambda_s$CDM setup the model deviates from $\Lambda$CDM in the early universe due to the existence of relativistic neutrino-like species. Using the Boltzmann solver CLASS in combination with MontePython we confront predictions of the stringy model to experimental data (from the Planck mission, Pantheon+ supernova type Ia, BAO, and KiDS-1000). We show that the string-inspired model provides a satisfactory fit to the data and can resolve the cosmological tensions.

We explore halo assembly bias on cluster scales using large samples of superclusters. Leveraging the largest-ever X-ray galaxy cluster and supercluster samples obtained from the first SRG/eROSITA all-sky survey, we construct two subsamples of galaxy clusters which consist of supercluster members (SC) and isolated clusters (ISO) respectively. After correcting the selection effects on redshift, mass, and survey depth, we compute the excess in the concentration of the intracluster gas of isolated clusters with respect to supercluster members, defined as $\delta c_{\rm gas} \equiv c_{\rm gas,ISO}/c_{\rm gas,SC}-1$, to investigate the environmental effect on the concentration of clusters, an inference of halo assembly bias on cluster scales. We find that the average gas mass concentration of isolated clusters is a few percent higher than that of supercluster members, with a significance of $2\sigma$. The result on $\delta c_{\rm gas}$ varies with the overdensity ratio $f$ in supercluster identification, cluster mass proxies, and mass ranges, but remains positive in all the measurements. We measure slightly larger $\delta c_{\rm gas}$ when adopting a higher $f$ in supercluster identification. $\delta c_{\rm gas}$ is also larger for low-mass clusters, and almost negligible for high-mass clusters. We perform weak lensing analyses to compare the total mass concentration of the two classes and find a similar trend as obtained from gas mass concentration. Our results are consistent with the prediction of HAB on cluster scales, where halos located in denser environments are less concentrated, and this trend is stronger for low-mass halos than for massive halos. These phenomena can be interpreted by the fact that clusters in denser environments such as superclusters have experienced more mergers than isolated clusters in their assembling history.

Weakly interacting massive particles (WIMPs) are well-motivated candidates for dark matter. One signature of galactic WIMPs is the annual modulation expected in a detector's interaction rate, which arises from Earth's rotation around the Sun. Over two decades, the DAMA/LIBRA experiment has observed such modulation with 250 kg of NaI(Tl) scintillators, in accordance with WIMP expectations but inconsistent with the negative results of other experiments. The signal, however, depends on the target, so in order to either validate or refute the DAMA signal it is necessary to replicate the experiment using the same material. This is the goal of the ANAIS-112 experiment, currently underway since August 2017 with 112.5 kg of NaI(Tl). In this work, we present a reanalysis of three years of data employing an improved analysis chain to enhance the experimental sensitivity. The results presented here are consistent with the absence of modulation and inconsistent with DAMA's observation at nearly 3$\sigma$ C.L., with the potential to reach a 5$\sigma$ level within 8 years of data collection. Additionally, we explore the impact of different scintillation quenching factors in the comparison between ANAIS-112 and DAMA/LIBRA.

The Archives of Photographic PLates for Astronomical USE (APPLAUSE) project is aimed at digitising astronomical photographic plates from three major German plate collections, making them accessible through integration into the International Virtual Observatory (IVO). Photographic plates and related materials (logbooks, envelopes, etc.) were scanned with commercial flatbed scanners. Astrometric and photometric calibrations were carried out with the developed PyPlate software, using Gaia EDR3 data as a reference. The APPLAUSE data publication complies with IVO standards. The latest data release contains images and metadata from 27 plate collections from the partner institutes in Hamburg, Bamberg, and Potsdam, along with digitised archives provided by Tautenburg, Tartu, and Vatican observatories. Altogether, over two billion calibrated measurements extracted from about 70,000 direct photographic plates can readily be used to create long-term light curves. For instance, we constructed the historic light curve of the enigmatic dipping star KIC 8462852. We found no evidence of previously assumed variations on timescales of decades in our light curve. Potential uses of APPLAUSE images for transient sources can be appreciated by following the development of the nova shell of GK Per (1901) over time and the change in brightness of two extragalactic supernovae. The database holds about 10,000 spectral plates. We made use of objective prism plates to follow the temporal changes of Nova DN Gem through 1912 and 1913, highlighting an outburst in early 1913.

The star formation rate (SFR), the number of stars formed per unit of time, is a fundamental quantity in the evolution of the Universe. While turbulence is believed to play a crucial role in setting the SFR, the exact mechanism remains unclear. Turbulence promotes star formation by compressing the gas, but also slows it down by stabilizing the gas against gravity. Most widely-used analytical models rely on questionable assumptions, including: $i)$ integrating over the density PDF, a one-point statistical description that ignores spatial correlation, $ii)$ selecting self-gravitating gas based on a density threshold that often ignores turbulent dispersion, $iii)$ assuming the freefall time as the timescale for estimating SFR without considering the need to rejuvenate the density PDF, $iv)$ assuming the density PDF to be lognormal. Improving upon the only existing model that incorporates the spatial correlation of the density field, we present a new analytical model. We calculate the time needed to rejuvenate density fluctuations of a given density and spatial scale, revealing that it is generally much longer than the freefall time, rendering the latter inappropriate for use. We make specific predictions regarding the role of the Mach number, $ M $, and the driving scale of turbulence divided by the mean Jeans length. At low to moderate Mach numbers, turbulence does not reduce and may even slightly promote star formation by broadening the PDF. However, at higher Mach numbers, most density fluctuations are stabilized by turbulent dispersion, leading to a steep drop in the SFR as the Mach number increases.

Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Within the gravo-turbulent framework, we assume that turbulence shapes the ISM, creating density fluctuations that, if gravitationally unstable, will collapse and form stars. The goal of this work is to quantify how different regimes of turbulence, characterized by the strength and compressibility of the driving, shape the density field. We are interested in the outcome in terms of SFR and how it compares with existing analytical models for the SFR. We run a series of hydrodynamical simulations of turbulent gas. The simulations are first conducted without gravity, so that the density and velocity are shaped by the turbulence driving. Gravity is then switched on, and the SFR is measured and compared with analytical models. The physics included in these simulations is very close to the one assumed in the classical gravo-turbulent SFR analytical models, which makes the comparison straightforward. We found that the existing analytical models convincingly agree with simulations at low Mach number, but we measure a much lower SFR in the simulation with a high Mach number. We develop, in a companion paper, an updated physically-motivated SFR model that reproduces well the inefficient high Mach regime of the simulations. Our work demonstrates that accurate estimations of the turbulent-driven replenishment time of dense structures and the dense gas spatial distribution are necessary to correctly predict the SFR in the high Mach regime. The inefficient high-Mach regime is a possible explanation for the low SFR found in dense and turbulent environments such as the centers of our Milky Way and other galaxies.

Einstein's field equations with a background source term that induces perturbations and the applications of this new formalism to a compact dense matter relativistic star are presented. We introduce a new response kernel in the field equations between the metric and fluid perturbations. A source term which drives or induces the sub-hydro mesoscopic scales perturbations in the astrophysical system, is of importance here. Deterministic as well as stochastic perturbations for the radial case are worked out as solutions of field equations. We also touch upon polar perturbations that are deterministic with oscillatory parts. The stochastic perturbation are of significance in terms of two point or point separated correlations which form the building blocks for studying equilibrium and non-equilibrium statistical mechanics for the system. Our main aim is to build a theory for intermediate scale physics, for dense exotic matter and investigate structure of the compact astrophysical objects. Specifically, turbulence which connects various scales in superfluid matter in the dense stars is an area gaining importance. The work presented here is the starting point for new theoretical frame work that touches upon various yet unexplored scales where mechanical and dynamical effects interior to the matter of the star are of significance. Thus there is scope for extending studies in asteroseismology to mesoscopic effects at the new intermediate scales in the cold dense matter fluid . This is expected to enable us to probe astrophysical features at refined scales through theoretical formulations as well as for observational consequences.

The momentum distribution of particles accelerated at strong non-relativistic shocks may be influenced by the spatial distribution of the flow speed around the shock. This phenomenon becomes evident in the cosmic-ray modified shock, where the particle spectrum itself determines the flow velocity profile upstream. However, what if the flow speed is not uniform downstream as well? Hydrodynamics indicates that its spatial variation over the length scales involved in the acceleration of particles in supernova remnants (SNRs) could be noticeable.} {In the present paper, we address this issue, initially following Bell's approach to particle acceleration and then by solving the kinetic equation. We obtained an analytical solution for the momentum distribution of particles accelerated at the cosmic-ray modified shock with spatially variable flow speed downstream.} {We parameterized the downstream speed profile to illustrate its effect on two model cases, the test particle and non-linear acceleration at the shock.The resulting particle spectrum is generally softer in Sedov SNRs because the flow speed distribution reduces the overall shock compression accessible to particles with higher momenta. On the other hand, the flow structure in young SNRs could lead to harder spectra. The diffusive properties of particles play a crucial role as they determine the distance from the shock, and, as a consequence, the flow speed that particles encounter downstream. We discuss the effect of the plasma velocity gradient to be (partially) responsible for the evolution of the radio index and for the high-energy break visible in gamma rays from some SNRs. We expect that the effect from the gradient of the flow velocity downstream could be prominent in regions of SNRs with higher diffusion coefficient and lower magnetic field, i.e. where acceleration of particles is not very efficient.

We use Faraday rotation measurements from the latest catalog LoTSS DR2 from LOFAR to probe intergalactic magnetic fields. To identify the extragalactic component of the observed rotation measure (RM) we use two different techniques: residual rotation measure (RRM) and close radio pairs. For the RRM approach, we conclude that, despite smaller measurement errors in the LOFAR data, robust and conservative treatment of the systematic uncertainties in the Galactic contribution to RM results in the constraint on a homogeneous volume-filling magnetic field at the level 2.4 nG, slightly weaker than previous constraints from NVSS data, and does not allow to probe the presence of over-magnetized bubbles predicted by the AGN feedback model of the IllustrisTNG. Analyzing close radio pairs we found that in only 0.5% of our mock realizations of observed data, the expected contribution from the over-magnetized bubbles does not exceed LoTSS DR2 data.

We apply a methodology to build a sample of extreme emission line galaxies (EELGs) using integral field spectroscopy data. In this work we follow the spectroscopic criteria corresponding for EELG selection and use the MUSE Hubble Ultra-Deep field survey, which includes the deepest spectroscopic survey ever performed. Objects in the primary (extended) sample were detected requiring a rest-frame equivalent width EWo $\geqslant$ 300A (200A $\leq$ EWo $\leq$ 300A) in any of the emission lines of [OII]$\lambda\lambda$3726,29, [OIII]$\lambda\lambda$5007,4959, or H$\alpha$. A detailed closer inspection of the spectra of the candidates selected has been performed on a one by one basis, in order to confirm their classification. For this sample, the line fluxes, physical properties and chemical abundances of the EELGs have been derived as well as their spatially resolved structure and kinematics. Four (five) of the galaxies in the primary (extended) sample, $\sim$57$\%$ ($\sim$83$\%$), were spatially resolved. Three (none) of them present a clear pattern compatible with rotation. We have shown how our entire EELGs sample share the same loci defined by high-redshift galaxies (z $\approx$ 6-8) for the mass-metallicity relation, illustrating their role as local analogs.

It is well established that black holes reside in the central regions of virtually all types of known galaxies. Recent observational and numerical studies however challenge this picture, suggesting that intermediate-mass black holes in dwarf galaxies may be found on orbits far from the center. In particular, constant-density cores minimize orbital energy losses due to dynamical friction, and allow black holes to settle on stable off-center orbits. Using controlled simulations, we study the dynamics of off-center black holes in dwarf galaxies with such cores. We propose a new scenario to describe off-center mergers of massive black holes, starting with a Jacobi capture. We focus on initially circular and co-planar black hole orbits and explore a large parameter space of black hole masses and orbital parameters. We find that Jacobi captures are a complex and chaotic phenomenon that occurs in about 13% of cases in this simplified setup, and we quantify how the likelihood of capture depends on the simulation parameters. We note that this percentage is likely an upper limit of the general case. Nevertheless, we show that Jacobi captures in cored dwarf galaxies can facilitate the formation of off-center black hole binaries, and that this process is sufficiently common to have a substantial effect on black hole growth. While our setup only allows for temporary captures, we expect dissipative forces from baryons and post-Newtonian corrections to maintain the captures over time and to lead to the formation of stable binary systems. This motivates future studies of the effectiveness of such dissipative forces, within stripped nuclei or globular clusters, in forming stable bound systems.

We test the potential of Bayesian synthesis of upcoming multi-instrument data to extract orbital parameters and individual light curves of close binary supermassive black holes (CB-SMBH) with subparsec separations. Next generation (ng) interferometers, will make possible the observation of astrometric wobbles in CB-SMBH. Combining them with periodic variable time-domain data from surveys like the Vera C. Rubin Legacy Survey of Space and Time (LSST), allows for a more information on CB-SMBH candidates compared to standalone observational methods. Our method reliably determines binary parameters and component fluxes from binary total flux across long-term, intermediate and short-term binary dynamics and observational configurations, assuming ten annual observations, even in short period "q-accrete" objects. Expected CB-SMBH astrometric wobbles constructed from binary dynamical parameters, might serve in refining observational strategies for CB-SMBH. Combination of inferred mass ratio, light curves of binary components, and observed photocenter wobbles can be a proxy for the activity states of CB-SMBH components.

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to SNe Ia. The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of $^{56}\text{Ni}$ in SNe Ia is produced for different temperatures, densities, and proton to neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of $\alpha$-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and $^{44}\text{Ti}$ and $^{56}\text{Ni}$ abundances by each model. The consequences of $\mathtt{y}_e \neq 0.5$ in the initial composition are also explored where $\mathtt{y}_e$ is the electron fraction. The simulated results show that $^{56}\text{Ni}$ production decreases by decreasing $\mathtt{y}_e$ as expected, and that the $^{43}\text{Sc}$ is a key isotope in proton and neutron channels toward $^{56}\text{Ni}$ production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the $^{56}\text{Ni}$ abundance in SNe Ia for $\mathtt{y}_e \lesssim 0.5$ initial conditions.

Small rocky planets, as well as larger planets that suffered extensive volatile loss, tend to be drier and have thinner atmospheres as compared to Earth. Such planets probably outnumber worlds better endowed with volatiles, being the most common habitable planets. For the subgroup of fast rotators following eccentric orbits, atmospheres suffer radiative forcing and their heat capacity provides a method for gauging atmospheric thickness and surface conditions. We further explore the model presented in a previous paper and apply it to real and hypothetical exoplanets in the habitable zone of various classes of stars, simulating atmospheric and orbital characteristics. For planetary eccentricities e ~0.3, the forcing-induced hypothetical temperature variation would reach ~80 K for airless planets and ~10 K for planets with substantial atmospheres. For Kepler-186 f and Kepler-442 b, assuming e ~0.1, temperature variations can reach ~24 K. We also consider habitable exomoons in circular orbits around gas giants within the habitable zone, which suffer radiative forcing due to their epicyclic motion. We study several combinations of parameters for the characterization of planets (mass, eccentricity and semi-major axis) and exomoons (mass, orbital radius, albedo and atmospheric characteristics) for different stellar types. For e ~0.3, exomoon temperature varies up to ~90 K, while for ~0.6 variations can reach ~200 K. Such exomoons may plausibly retain their volatiles by continued volcanic activity fueled by tidal dissipation. Although currently undetectable, such effects might be within reach of future ELT-class telescopes and space missions with mid-infrared and coronagraphic capabilities.

The development of the latest generation of Imaging Atmospheric Cherenkov Telescopes (IACTs) over recent decades has led to the discovery of new extreme astrophysical phenomena in the very-high-energy (VHE, E > 100 GeV) gamma-ray regime. Time-domain and multi-messenger astronomy are inevitably connected to the physics of transient VHE emitters, which show unexpected (and mostly unpredictable) flaring or exploding episodes at different timescales. These transients often share the physical processes responsible for the production of the gamma-ray emission, through cosmic-ray acceleration, magnetic reconnection, jet production and/or outflows, and shocks interactions. In this review, we present an up-to-date overview of the VHE transients field, spanning from novae to supernovae, neutrino counterparts or fast radio bursts, among others, and we outline the expectations for future facilities.

The identification of extragalactic fast optical transients (eFOTs) as potential multi-messenger sources is one of the main challenges in time-domain astronomy. However, recent developments have allowed for probes of rapidly-evolving transients. With the increasing number of alert streams from optical time-domain surveys, the next paradigm is building technologies to rapidly identify the most interesting transients for follow-up. One effort to make this possible is the fitting of objects to a variety of eFOT lightcurve models such as kilonovae and $\gamma$-ray burst (GRB) afterglows. In this work, we describe a new framework designed to efficiently fit transients to light curve models and flag them for further follow-up. We describe the pipeline's workflow and a handful of performance metrics, including the nominal sampling time for each model. We highlight as examples ZTF20abwysqy, the shortest long gamma ray burst discovered to date, and ZTF21abotose, a core-collapse supernova initially identified as a potential kilonova candidate.

Global magnetic fields of early-type stars are commonly characterised by the mean longitudinal magnetic field $\langle B_{\rm z} \rangle$ and the mean field modulus $\langle B \rangle$, derived from the circular polarisation and intensity spectra, respectively. Observational studies often report a root mean square (rms) of $\langle B_{\rm z} \rangle$ and an average value of $\langle B \rangle$. In this work, I used numerical simulations to establish statistical relationships between these cumulative magnetic observables and the polar strength, $B_{\rm d}$, of a dipolar magnetic field. I show that in the limit of many measurements randomly distributed in rotational phase, $\langle B_{\rm z} \rangle_{\rm rms}$=$0.179^{+0.031}_{-0.043}$$B_{\rm d}$ and $\langle B \rangle_{\rm avg}$=$0.691^{+0.020}_{-0.023}$$B_{\rm d}$. The same values can be recovered with only three measurements, provided that the observations are distributed uniformly in the rotational phase. These conversion factors are suitable for ensemble analyses of large stellar samples, where each target is covered by a small number of magnetic measurements.

Mercury is notoriously difficult to form in solar system simulations, due to its small mass and iron-rich composition. Smooth particle hydrodynamics simulations of collisions have found that a Mercury-like body could be formed by one or multiple giant impacts, but due to the chaotic nature of collisions it is difficult to create a scenario where such impacts will take place. Recent work has found more success forming Mercury analogues by adding additional embryos near Mercury's orbit. In this work, we aim to form Mercury by simulating the formation of the solar system in the presence of the giant planets Jupiter and Saturn. We test out the effect of an inner disk of embryos added on to the commonly-used narrow annulus of initial material. We form Mercury analogues with core-mass fractions (CMF) $> 0.4$ in $\sim 10\%$ of our simulations, and twice that number of Mercury analogues form during the formation process, but are unstable and do not last to the end of the simulations. Mercury analogues form at similar rates for both disks with and without an inner component, and most of our Mercury analogues have lower CMF than that of Mercury, $\sim 0.7$, due to significant accretion of debris material. We suggest that a more in-depth understanding of the fraction of debris mass that is lost to collisional grinding is necessary to understand Mercury's formation, or some additional mechanism is required to stop this debris from accreting.

Context. Recent studies of the optical depth comparing [12CII] and [13CII] line profiles in Galactic star-forming regions revealed strong self-absorption in [12CII] by low excitation foreground material, implying a large column density of C+ corresponding to an equivalent AV of a few, up to about 10 mag. Aims. As the nature and origin of such a large column of cold C+ foreground gas are difficult to explain, it is essential to constrain the physical conditions of this material. Methods. We conducted high-resolution observations of [OI] 63 um and [OI] 145 um lines in M17 SW and Mon R2. The [OI] 145 um transition traces warm PDR-material, while the [OI] 63 um line traces foreground material as manifested by absorption dips. Results. Comparison of both [OI] line profiles with [CII] isotopic lines confirms warm PDR-origin background emission and a significant column of cold foreground material causing self-absorption visible in [12CII] and [OI] 63 um profiles. In M17 SW, the C+ and O column densities are comparable for both layers. Mon R2 exhibits larger O columns compared to C+, indicating additional material where the carbon is neutral or in molecular form. Small-scale spatial variation of the foreground absorption profiles and the large column density (around 1E18 cm-2 ) of the foreground material suggest emission from high-density regions associated with the cloud complex, not a uniform diffuse foreground cloud. Conclusions. The analysis confirms that the previously detected intense [CII] foreground absorption is attributable to a large column of low excitation dense atomic material, where carbon is ionized, and oxygen is in neutral atomic form.

At the moment, there are two neutron star X-ray binaries with massive red supergiants as donors. De et al. (2023) proposed that the system SWIFT J0850.8-4219 contains a neutron star at the propeller stage. We study this possibility by applying various models of propeller spin-down. We demonstrate that the duration of the propeller stage is very sensitive to the regime of rotational losses. Only in the case of a relatively slow propeller model proposed by Davies and Pringle (1981), the duration of the propeller is long enough to provide a significant probability to observe the system at this stage. Future determination of the system parameters (orbital and spin periods, magnetic field of the compact object, etc.) will allow putting strong constraints on the propeller behavior.

Improved polarization measurements at frequencies below 70 GHz with degree-level angular resolution are crucial for advancing our understanding of the Galactic synchrotron radiation and the potential polarized anomalous microwave emission and ultimately benefiting the detection of primordial $B$ modes. In this study, we present sensitivity-improved 40 GHz polarization maps obtained by combining the CLASS 40 GHz and WMAP $Q$-band data through a weighted average in the harmonic domain. The decision to include WMAP $Q$-band data stems from similarities in the bandpasses. Leveraging the accurate large-scale measurements from WMAP $Q$ band and the high-sensitivity information from CLASS 40 GHz band at intermediate scales, the noise level at $\ell\in[30, 100]$ is reduced by a factor of $2-3$ in the map space. A pixel domain analysis of the polarized synchrotron spectral index ($\beta_s$) using WMAP $K$ band and the combined maps (mean and 16/84th percentile across the $\beta_s$ map: $-3.08_{-0.20}^{+0.20}$) reveals a stronger preference for spatial variation (PTE for a uniform $\beta_s$ hypothesis smaller than 0.001) than the results obtained using WMAP $K$ and $Ka$ bands ($-3.08_{-0.14}^{+0.14}$). The cross-power spectra of the combined maps follow the same trend as other low-frequency data, and validation through simulations indicates negligible bias introduced by the combination method (sub-percent level in the power spectra). The products of this work are publicly available on $\mathtt{LAMBDA}$.

Gaia BH3 is the first observed dormant black hole (BH) with a mass of $\approx{30}$ M$_\odot$ and represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3 focusing on the old and metal-poor stellar population of the Milky Way halo. We use our open-source population synthesis code SEVN to evolve $3.2 \times 10^8$ binaries exploring 16 sets of parameters. We find that systems like Gaia BH3 form preferentially from binaries initially composed of a massive star ($40-60$ M$_\odot$) and a low mass companion ($<1$ M$_\odot$) in a wide ($P>10^3$ days) and eccentric orbit ($e>0.6$). Such progenitors do not undergo any Roche-lobe overflow episode during their entire evolution, so that the final orbital properties of the BH-star system are essentially determined at the core collapse of the primary star. Low natal kicks ($\approx$ 10~km/s) significantly favour the formation of Gaia BH3-like systems, but high velocity kicks up to $\approx 220$ km/s are also allowed. We estimate the formation efficiency for Gaia BH3-like systems in old ($t>10$ Gyr) and metal-poor ($Z<0.01$) populations to be $4 \times 10^{-8}$ M$_\odot^{-1}$, representing $\approx 3\%$ of the whole simulated BH-star population. We expect up to $\sim 3000$ BH-star systems in the Galactic halo formed through isolated evolution, of which $\sim 100$ are compatible with Gaia BH3 alike. Gaia BH3-like systems represent a common product of isolated binary evolution at low metallicity ($Z<0.01$), but given the steep density profile of the Galactic halo we do not expect more than one in the halo at the observed distance of Gaia BH3. Considering the estimated formation efficiency for both the isolated and dynamical formation channel, we conclude that they are almost equally likely to explain the origin of Gaia BH3.

In our quest to understand the generation of cosmological perturbations, we face two serious obstacles: we do not have direct information about the environment experienced by primordial perturbations during inflation, and our observables are practically limited to correlators of massless fields, heavier fields and derivatives decaying exponentially in the number of e-foldings. The flexible and general framework of open systems has been developed precisely to face similar challenges. Building on previous work, we develop a Schwinger-Keldysh path integral description for an open effective field theory of inflation, describing the possibly dissipative and non-unitary evolution of the Goldstone boson of time translations interacting with an unspecified environment, under the key assumption of locality in space and time. Working in the decoupling limit, we study the linear and interacting theory in de Sitter and derive predictions for the power spectrum and bispectrum that depend on a finite number of effective couplings organised in a derivative expansion. The smoking gun of interactions with the environment is an enhanced but finite bispectrum close to the folded kinematical limit. We demonstrate the generality of our approach by matching our open effective theory to an explicit model. Our construction provides a standard model to simultaneously study phenomenological predictions as well as quantum information aspects of the inflationary dynamics.

The spontaneous breaking of a $U(1)$ symmetry via an intermediate discrete symmetry may yield a hybrid topological defect of "domain walls bounded by cosmic strings". The decay of this defect network leads to a unique gravitational wave signal spanning many orders in observable frequencies, that can be distinguished from signals generated by other sources. We investigate the production of gravitational waves from this mechanism in the context of the type-I two-Higgs-doublet model extended by a $U(1)_R$ symmetry, that simultaneously accommodates the seesaw mechanism, anomaly cancellation, and eliminates flavour-changing neutral currents. The gravitational wave spectrum produced by the string-bounded-wall network can be detected for $U(1)_R$ breaking scale from $10^{12}$ to $10^{15}$ GeV in forthcoming interferometers including LISA and Einstein Telescope, with a distinctive $f^{3}$ slope and inflexion in the frequency range between microhertz and hertz.

After cosmic inflation, coherent oscillations of the inflaton field about a monomial potential $V(\phi)\sim \phi^k$ result in an expansion phase characterized by a stiff equation-of-state $w\simeq(k-2)/(k+2)$. Sourced by the oscillating inflaton condensate, parametric (self)resonant effects can induce the exponential growth of inhomogeneities eventually backreacting and leading to the fragmentation of the condensate. In this work, we investigate realizations of inflation giving rise to such dynamics, assuming an inflaton weakly coupled to its decay products. As a result, the transition to a radiation-dominated universe, i.e. reheating, occurs after fragmentation. We estimate the consequences on the production of gravitational waves by computing the contribution induced by the stiff equation-of-state era in addition to the signal generated by the fragmentation process for $k=4,6,8,10$. Our results are independent of the reheating temperature provided that reheating is achieved posterior to fragmentation. Our work shows that the dynamics of such weakly-coupled inflaton scenario can actually result in characteristic gravitational wave spectra with frequencies from Hz to GHz, in the reach of future gravitational wave observatories, in addition to the complementarity between upcoming detectors in discriminating (post)inflation scenarios. We advocate the need of developing high-frequency gravitational wave detectors to gain insight into the dynamics of inflation and reheating.

Leptogenesis typically requires the introduction of heavy particles whose out-of-equilibrium decays are essential for generating a matter-antimatter asymmetry, according to one of Sakharov's conditions. We demonstrate that in Dirac leptogenesis, scatterings between the light degrees of freedom - Standard Model particles plus Dirac neutrinos - are sufficient to generate the asymmetry. Due to its vanishing source term in the Boltzmann equations, the asymmetry of right-handed neutrinos solely arises through wash-in processes. Sakharov's conditions are satisfied because the right-handed neutrino partners are out of equilibrium. Consequently, heavy degrees of freedom never needed to be produced in the early universe, allowing for a reheating temperature well below their mass scale. Considering a minimal leptoquark model, we discuss the viable parameter space along with the observational signature of an increased number of effective neutrinos in the early universe.

The expansion and upgrade of the global network of ground-based gravitational wave detectors promises to improve our capacity to infer the sky-localization of transient sources, enabling more effective multi-messenger follow-ups. At the same time, the increase in the signal-to-noise ratio of detected events allows for more precise estimates of the source parameters. This study aims to assess the performance of advanced-era networks of ground-based detectors, focusing on the Hanford, Livingston, Virgo, and KAGRA instruments. We use full Bayesian parameter estimation procedures to predict the scientific potential of a network. Assuming a fixed LIGO configuration, we find that the addition of the Virgo detector is beneficial to the sky localization starting from a binary neutron star horizon distance of 20 Mpc and improves significantly from 40 Mpc onwards for both a single and double LIGO detector network, reducing the inferred mean sky-area by up to 95%. Similarly, the KAGRA detector tightens the constraints, starting from a sensitivity range of 10 Mpc. Looking at highly-spinning binary black holes, we find significant improvements with increasing sensitivity in constraining the intrinsic source parameters when adding Virgo to the two LIGO detectors. Finally, we also examine the impact of the low-frequency cut-off data on the signal-to-noise ratio. We find that existing 20 Hz thresholds are sufficient and propose a metric to monitor this to study detector performance. Our findings quantify how future enhancements in detector sensitivity and network configurations will improve the localization of gravitational wave sources and allow for more precise identification of their intrinsic properties.

A Bayesian method is used in this extensive work to generate a large set of minimally constrained equations of state (EOSs) for matters in neutron stars (NS). These EOSs are analyzed for their correlations with key NS properties, such as the tidal deformability, radius, and maximum mass, within the mass range of $1.2-2M_\odot$. The observed connections between the pressure of $\beta$-equilibrated matter and the properties of neutron stars at different densities offer significant insights into the behavior of NS matter in a nearly model-independent manner. The study also examines the influence of various factors on the correlation of symmetry energy parameters, such as slope and curvature parameters at saturation density ($\rho_0=0.16 ~\text{fm}^{-3}$) with the tidal deformability and radius of neutron stars. This study investigates the robustness of the observed correlations by considering the distributions and interdependence of symmetry energy parameters. Furthermore, the utilization of Principal Component Analysis (PCA) is employed to unveil the complicated relationship between various nuclear matter parameters and properties of neutron stars. This analysis highlights the importance of employing multivariate analysis techniques in order to comprehend the variety in tidal deformability and radius observed across distinct masses of NS. This comprehensive study aims to establish a connection between the parameters of nuclear matter and the properties of neutron stars, providing significant insights into the behavior of NS matter across different circumstances.

We investigate a Schwarzschild metric exhibiting a signature change across the event horizon, which gives rise to what we term a Lorentzian-Euclidean black hole. The resulting geometry is regularized by employing the Hadamard partie finie technique, which allows us to prove that the metric represents a solution of vacuum Einstein equations. In this framework, we introduce the concept of atemporality as the dynamical mechanism responsible for the transition from a regime with a real-valued time variable to a new one featuring an imaginary time. We show that this mechanism prevents the occurrence of the singularity and, by means of the regularized Kretschmann invariant, we discuss in which terms atemporality can be considered as the characteristic feature of this black hole.

We analyse the search behaviour of genetic programming for symbolic regression in practically relevant but limited settings, allowing exhaustive enumeration of all solutions. This enables us to quantify the success probability of finding the best possible expressions, and to compare the search efficiency of genetic programming to random search in the space of semantically unique expressions. This analysis is made possible by improved algorithms for equality saturation, which we use to improve the Exhaustive Symbolic Regression algorithm; this produces the set of semantically unique expression structures, orders of magnitude smaller than the full symbolic regression search space. We compare the efficiency of random search in the set of unique expressions and genetic programming. For our experiments we use two real-world datasets where symbolic regression has been used to produce well-fitting univariate expressions: the Nikuradse dataset of flow in rough pipes and the Radial Acceleration Relation of galaxy dynamics. The results show that genetic programming in such limited settings explores only a small fraction of all unique expressions, and evaluates expressions repeatedly that are congruent to already visited expressions.

We analyze the effective four-dimensional dynamics of the extra-dimensional moduli fields in curved braneworlds having nested warping, with particular emphasis on the doubly warped model which is interesting in the light of current collider constraints on the mass of the Kaluza-Klein graviton. The presence of a non-zero brane cosmological constant ($\Omega$) naturally induces an effective moduli potential in the four-dimensional action, which shows distinct features in dS ($\Omega>0$) and AdS ($\Omega<0$) branches. For the observationally interesting case of dS 4-branes, a metastable minimum in the potential arises along the first modulus, with no minima along the higher moduli. The underlying nested geometry also leads to interesting separable forms of the non-canonical kinetic terms in the Einstein frame, where the brane curvature directly impacts the kinetic properties of only the first modulus. We subsequently explore the ability of curved multiply warped geometries to drive inflation with an in-built exit mechanism, by assuming predominant slow roll along each modular direction on a case-by-case basis. We find slow roll on top of the metastable plateau along the first modular direction to be the most viable scenario, with the higher-dimensional moduli parametrically tuning the height of the potential without significant impact on the inflationary observables. On the other hand, while slow roll along the higher moduli can successfully inflate the background and eventually lead to an exit, consistency with observations seemingly requires unphysical hierarchies among the extra-dimensional radii, thus disfavouring such scenarios.

The origin of neutrino masses remains unknown. Both the vacuum mass and the dark mass generated by the neutrino interactions with DM particles or fields can fit the current oscillation data. The dark mass squared is proportional to the DM number density and therefore varies on the galactic scale with much larger values around the Galactic Center. This affects the group velocity and the arrival time delay of core-collapse supernovae neutrinos. This time delay, especially for the $\nu_e$ neutronization peak with a sharp time structure, can be used to distinguish the vacuum and dark neutrino masses. For illustration, we explore the potential of DUNE which is sensitive to $\nu_e$. Our simulations show that DUNE can distinguish the two neutrino mass origins at more than $5\sigma\,$C.L., depending on the observed local value of neutrino mass, the neutrino mass ordering, the DM density profile, and the SN location.

Every signal propagating through the universe is diffracted by the gravitational fields of intervening objects, aka gravitational lenses. Diffraction is most efficient when caused by compact lenses, which invariably produce additional images of a source. The signals associated with additional images are generically faint, but their collective effect may be detectable with coherent sources, such as gravitational waves (GWs), where both amplitude and phase are measured. Here, I describe lens stochastic diffraction (LSD): Poisson-distributed fluctuations after GW events caused by compact lenses. The amplitude and temporal distribution of these signals encode crucial information about the mass and abundance of compact lenses. Through the collective stochastic signal, LSD offers an order-of-magnitude improvement over single lens analysis for objects with mass $\gtrsim 10^3 M_\odot$. This gain can improve limits on compact dark-matter halos and allows next-generation instruments to detect supermassive black holes, given the abundance inferred from quasar luminosity studies.

We consider the nuclear absorption of dark matter as an alternative to the typical indirect detection search channels of dark matter decay or annihilation. In this scenario, an atomic nucleus transitions to an excited state by absorbing a pseudoscalar dark matter particle and promptly emits a photon as it transitions back to its ground state. The nuclear excitation of carbon and oxygen in the Galactic Center would produce a discrete photon spectrum in the $\mathcal{O}(10)$ MeV range that could be detected by gamma-ray telescopes. Using the \texttt{BIGSTICK} large-scale shell-model code, we calculate the excitation energies of carbon and oxygen. We constrain the dark matter-nucleus coupling for current COMPTEL data, and provide projections for future experiments AMEGO-X, e-ASTROGAM, and GRAMS for dark matter masses from $\sim$ 10 to 30 MeV. We find the excitation process to be very sensitive to the dark matter mass and find that the future experiments considered would improve constraints on the dark matter-nucleus coupling within an order of magnitude.