Papers for Thursday, Nov 23 2023

We present an extensive study of the effects of neutrino transport in 3-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. We consider a total of 8 BNS configurations, while varying equation of state models, mass ratios and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and post-merger phases, the properties of ejecta and outflow, and the post-merger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties and provides insight into the effects of neutrino trapping in BNS mergers.

We construct an anti-halo void catalogue of $150$ voids with radii $R > 10\,h^{-1}\mathrm{\,Mpc}$ in the Local Super-Volume ($<135\,h^{-1}\mathrm{\,Mpc}$ from the Milky Way), using posterior resimulation of initial conditions inferred by field-level inference with Bayesian Origin Reconstruction from Galaxies (\codefont{BORG}). We describe and make use of a new algorithm for creating a single, unified void catalogue by combining different samples from the posterior. The catalogue is complete out to $135\,h^{-1}\mathrm{\,Mpc}$, with void abundances matching theoretical predictions. Finally, we compute stacked density profiles of those voids which are reliably identified across posterior samples, and show that these are compatible with $\Lambda$CDM expectations once environmental selection (e.g., the estimated $\sim 4\%$ under-density of the Local Super-Volume) is accounted for.

We compare the reduced void probability function (VPF) inside and outside of cosmic voids in the TNG300-1 simulation, both in real and simulated redshift space. The VPF is a special case of the counts-in-cells approach for extracting information of high-order clustering that is crucial for a full understanding of the distribution of galaxies. Previous studies have validated the hierarchical scaling paradigm of galaxy clustering moments, in good agreement with the "negative binomial" model, in redshift surveys, but have also reported that this paradigm is not valid in real space. However, in this work we find that hierarchical scaling can indeed be found in real space inside cosmic voids. This is well fitted by the negative binomial model. We find this result to be robust against changes in void identification, galaxy mass, random dilutions, and redshift. We also obtain that the VPF in real space at high redshift approaches the negative binomial model, and therefore it is similar to the VPF inside voids at the present time. This study points, for the first time, towards evidence of hierarchical scaling of high-order clustering of galaxies in real space inside voids, preserving the pristine structure formation processes of the Universe.

In our previous work, we derive the CR-IGIMF: a new scenario for a variable stellar initial mass function (IMF), which combines numerical results on the role played by cosmic rays in setting the thermal state of star forming gas, with the analytical approach of the integrated galaxy-wide IMF. In this work, we study the implications of this scenario for the properties of local Early-Type galaxies (ETG), as inferred from dynamical, photometric and spectroscopic studies. We implement a library of CR-IGIMF shapes in the framework of the Galaxy Evolution and Assembly (GAEA) model. Our realization includes a derivation of synthetic spectral energy distribution for each model galaxy, allowing a direct derivation of the mass fraction in the mean IMF of low-mass stars (i.e. the dwarf-to-giant ratio - $\rm f_{dg}$) and a comparison with IMF sensitive spectral features. The predictions of the GAEA model implementing the CR-IGIMF confirm our previous findings: it correctly reproduces both the observed excess of z$\sim$0 dynamical mass (mass-to-light ratios) with respect to spectroscopic (photometric) estimates assuming a universal, MW-like, IMF, and the observed increase of [$\alpha$/Fe] ratios with stellar mass in spheroidal galaxies. Moreover, this realization reproduces the increasing trends of $\rm f_{dg}$, and IMF-sensitive line-strengths with velocity dispersion, although the predicted relations are significantly shallower than the observed ones. Our results show that the CR-IGIMF is a promising scenario that reproduces at the same time dynamical, photometric and spectroscopic indications of a varying IMF in local ETGs. The shallow relations found for spectral indices suggest that either a stronger variability as a function of galaxy properties or additional dependences (e.g. as a function of star forming gas metallicity) might be required to match the strength of the observed trends.

The black hole transient GRS~1915+105 entered a new phase of activity in 2018, generally characterised by low X-ray and radio fluxes. This phase has been only interrupted by episodes of strong and variable radio emission, during which high levels of X-ray absorption local to the source were measured. We present 18 epochs of near-infrared spectroscopy (2018--2023) obtained with GTC/EMIR and VLT/X-shooter, spanning both radio-loud and radio-quiet periods. We demonstrate that radio-loud phases are characterised by strong P-Cygni line profiles, indicative of accretion disc winds with velocities of up to $\mathrm{\sim 3000~km~s^{-1}}$. This velocity is consistent with those measured in other black hole transients. It is also comparable to the velocity of the X-ray winds detected during the peak outburst phases in GRS~1915+105, reinforcing the idea that massive, multi-phase outflows are characteristic features of the largest and most powerful black hole accretion discs. Conversely, the evolution of the Br$\gamma$ line profile during the radio-quiet phases follows the expected trend for accretion disc lines in a system that is gradually decreasing its intrinsic luminosity, exhibiting weaker intensities and more pronounced double-peaks.

The emission of hard X-rays associated with white dwarfs (WD) can be generated by the presence of a stellar companion either by the companion's coronal emission or by an accretion disk formed by material stripped from the companion. Recent studies have suggested that a Jupiter-like planet can also be donor of material whose accretion onto the WD can generate hard X-rays. We use the {\sc guacho} code to reproduce the conditions of this WD-planet scenario. With the example of the hard X-ray WD KPD\,0005+5106, we explore different terminal wind velocities and mass-loss rates of a donor planet for a future network of simulations to investigate the luminosity and the spectral and temporal properties of the hard X-ray emission in WD-planet systems. Our simulations show that the material stripped from the planet forms a disk and accretes onto the WD to reach temperatures high enough to generate hard X-rays as usually seen in X-ray binaries with low-mass companions. For high terminal wind velocities, the planet material does not form a disk, but it rather accretes directly onto the WD surface. The simulations reproduce the X-ray luminosity of another X-ray accreting WD (G\,29$-$38), and only for some times reaches the hard X-ray luminosity of KPD\,0005+5106. The X-ray variability is stochastic and does not reproduce the period of KPD\,0005+5106, suggesting that additional physical processes (e.g., hot spots resulting from magnetic channelling of the accreting material) need to be explored.

The most abundant interstellar molecule, molecular Hydrogen (H$_{2}$), is practically invisible in cold molecular clouds. Astronomers typically use carbon monoxide (CO) to trace the bulk distribution and mass of H$_{2}$ in our galaxy and many others. CO observations alone fail to trace a massive component of molecular gas known as "CO-dark" gas. We present an ultra sensitive pilot search for the 18cm hydroxyl (OH) lines in the Andromeda Galaxy (M31) with the 100m Robert C. Byrd Green Bank Telescope. We successfully detected the 1667 and 1665 MHz OH in faint emission. The 1665/1667 MHz line ratio is consistent with the characteristic 5:9 ratio associated with local thermodynamic equilibrium (LTE). To our knowledge, this is the first detection of non-maser 18cm OH emission in another galaxy. We compare our OH and HI observations with archival CO (1-0) observations. Our OH detection position overlaps with the previously discovered Arp Outer Arm in CO. Our best estimates show that the amount of H$_{2}$ traced by OH is 140% higher than the amount traced by CO in this sightline. We show that the amount of dark molecular gas implied by dust data supports this conclusion. We conclude that the 18cm OH lines hold promise as a valuable tool for mapping of the "CO-dark" and "CO-faint" molecular gas phase in nearby galaxies, especially with upcoming multi-beam, phased-array feed receivers on radio telescopes which will allow for drastically improved mapping speeds of faint signals.

Contemporary three-dimensional physics-based simulations of the solar convection zone disagree with observations. They feature differential rotation substantially different from the true rotation inferred by solar helioseismology and exhibit a conveyor belt of convective "Busse" columns not found in observations. To help unravel this so-called "convection conundrum", we use a three-dimensional pseudospectral simulation code to investigate how radially non-uniform viscosity and entropy diffusivity affect differential rotation and convective flow patterns in density-stratified rotating spherical fluid shells. We find that radial non-uniformity in fluid properties enhances polar convection, which, in turn, induces non-negligible lateral entropy gradients that lead to large deviations from differential rotation geostrophy due to thermal wind balance. We report simulations wherein this mechanism maintains differential rotation patterns very similar to the true solar profile outside the tangent cylinder, although discrepancies remain at high latitudes. This is significant because differential rotation plays a key role in sustaining solar-like cyclic dipolar dynamos.

The variations of elemental abundance and their ratios along the Galactocentric radius result from the chemical evolution of the Milky Way disks. The $\rm ^{12}C/^{13}C$ ratio in particular is often used as a proxy to determine other isotopic ratios, such as $\rm ^{16}O/^{18}O$ and $\rm ^{14}N/^{15}N$. Measurements of $\rm ^{12}CN$ and $\rm ^{13}CN$ (or $\rm C^{15}N$) -- with their optical depths corrected via their hyper-fine structure lines -- have traditionally been exploited to constrain the Galactocentric gradients of the CNO isotopic ratios. Such methods typically make several simplifying assumptions (e.g. a filling factor of unity, the Rayleigh-Jeans approximation, and the neglect of the cosmic microwave background) while adopting a single average gas phase. However, these simplifications introduce significant biases to the measured $\rm ^{12}C/^{13}C$ and $\rm ^{14}N/^{15}N$. We demonstrate that exploiting the optically thin satellite lines of $\rm ^{12}CN$ constitutes a more reliable new method to derive $\rm ^{12}C/^{13}C$ and $\rm ^{14}N/^{15}N$ from CN isotopologues. We apply this satellite-line method to new IRAM 30-m observations of $\rm ^{12}CN$, $\rm ^{13}CN$, and $\rm C^{15}N$ $N=1\to0$ towards 15 metal-poor molecular clouds in the Galactic outer disk ($R_{\rm gc} > $ 12 kpc), supplemented by data from the literature. After updating their Galactocentric distances, we find that $\rm ^{12}C/^{13}C$ and $\rm ^{14}N/^{15}N$ gradients are in good agreement with those derived using independent optically thin molecular tracers, even in regions with the lowest metallicities. We therefore recommend using optically thin tracers for Galactic and extragalactic CNO isotopic measurements, which avoids the biases associated with the traditional method.

Neutron stars (NSs) and white dwarfs (WDs) are characterized by different geometric and physical properties, but their observed properties are often similar, making them difficult to distinguish. Therefore, it is desirable to search for their spectral features that could be easily identified from observations. We present spectral and timing signatures of NSs and WDs hosted in accreting X-ray binaries that can be easily identified from X-ray observations. We perform spectral and timing analysis of 4U~1636--53 and SS~Cygni, as typical representatives of such NS and WD binaries, based on their X-ray observations by RXTE, ASCA, Suzaku and BeppoSAX uising {\it Comptonization} spectral model. As a result, we formulate a criterion that makes it easy to distinguish NS from WD in such binaries: NS X-rays exhibits clear quasi-stable behavior with the index $\Gamma\to2$ and is characterized by quasi periodic oscillations (QPOs) at $\nu_{QPO} >0.5$~Hz, although WD X-rays is stable with $\Gamma \to1.85$ and is accompanied by QPOs at $\nu_{QPO}<0.05$~Hz during source outbursts. In addition, we revealed that in 4U~1636--53 the mHz QPOs anti-correlate with the plasma temperature, $T_e$ of Compton cloud (or the corona around a NS. This allowed us to associate mHz-QPOs with the corona dynamics during outburst cycle. The above index effect, now well established for 4U~1636--53 and SS~Cygni using extensive observations, has previously been found in other low-mass X-ray NS and WD binaries and agrees well with the criterion for distinguishing NSs and WDs presented here.

We measure the signal of secondary halo bias as a function of a variety of intrinsic and environmental halo properties, and characterize its statistical significance as a function of cosmological redshift. Using fixed and paired $N$-body simulations of dark-matter halos -- the \texttt{UNIT} simulation -- with masses above $10^{11}M_{\odot}h^{-1}$ identified over a wide range of cosmological redshifts ($0<z<5$), we explore the behavior of the scaling relations among different halo properties. We include novel environmental properties based on the halo distribution as well as the underlying dark-matter field. We implement an object-by-object estimator of large-scale effective bias and test its validity against standard approaches. With a bias assigned to each tracer, we perform a statistical analysis aiming at characterizing the distribution of bias and the signal of secondary halo bias. We show how the halo scaling relations linking direct probes of the halo potential well do not depend on the environment. On the contrary, links between halo mass and the so called set of secondary halo properties are sensitive to the cosmological environment. We show that the signal of secondary bias derives statistically from secondary correlations beyond the standard link to halo mass. We show that the secondary bias arise through non-local and/or environmental properties related either to the halo distribution or to the properties of the underlying dark-matter field. Properties such as the tidal field and the local Mach number generates the signals of secondary bias with the highest significance. We propose applications of the assignment of individual bias for the generation of mock catalogs containing the signal of secondary bias, as well as a series of cosmological analyses aiming at mining large galaxy data sets.

Direct imaging spectroscopy with future space-based telescopes will constrain terrestrial planet atmospheric composition and potentially detect biosignature gases. One promising indication of life is abundant atmospheric O2. However, various non-biological processes could also lead to O2 accumulation in the atmospheres of potentially habitable planets around Sun-like stars. In particular, the absence of non-condensible background gases such as N2 could result in appreciable H escape and abiotic O2 buildup, so identifying background atmosphere composition is crucial for contextualizing any O2 detections. Here, we perform retrievals on simulated directly imaged terrestrial planets using rfast, a new exoplanet atmospheric retrieval suite with direct imaging analysis capabilities. By simulating Earth-analog retrievals for varied atmospheric compositions, cloud properties, and surface pressures, we determine what wavelength range, spectral resolution, and signal-to-noise ratio (S/N) are necessary to constrain background gases' identity and abundance. We find N2 backgrounds can be uniquely identified with S/N$\sim$20 observations, provided that wavelength coverage extends beyond $\sim$1.6 $\mu$m to rule out CO-dominated atmospheres. Additionally, there is a low probability of O2-dominated atmospheres due to an O2-N2 degeneracy that is only totally ruled out at S/N$\sim$40. If wavelength coverage is limited to 0.2-1.1 $\mu$m, then although all other cosmochemically plausible backgrounds can be readily excluded, N2 and CO backgrounds cannot be distinguished. Overall, our simulated retrievals and associated integration time calculations suggest that near-infrared coverage to at least 1.6 $\mu$m and apertures approaching 8m are needed to confidently rule out O2 biosignature false positives within feasible integration times

New observations are presented of four evolved objects that display long, multi-year variations in their light curves. These are interpreted as good evidence of their binary nature, with the modulation caused by the barycenter motion of the evolved star resulting in a periodic obscuration by a circumbinary disk. Although protoplanetary nebulae (PPNe) commonly possess bipolar nebulae, which are thought to be shaped by a binary companion, there are very few PPNe in which a binary companion has been found. Three of the objects in this study appear to be PPNe, IRAS 07253-2001, 08005-2356, and 17542-0603, with long periods of 5.2, 6.9, and 8.2 yrs, respectively. The binary nature of IRAS 08005-2356 has recently been confirmed by a radial velocity study. Two samples, one of PPNe and the other of post-AGB star candidates, are investigated for further evidence on how common is a long-period light curve variation. Both samples suggest such light variations are not common. The fourth object, IRAS 20056+1834 (QY Sge), is an obscured RV Tau variable of the RVb subclass, with a long period of 3.9 yrs and pulsation periods of 102.9 and 51.5 days. The period of this object is seen to vary by 2%. Evidence is presented for a recent mass ejection in IRAS 17542-0603.

The cosmological 21-cm signal from neutral hydrogen, which is considered as a promising tool, is being used to observe and study the cosmic dawn (CD) and epoch of reionization (EoR). A significant part of this thesis focuses on the semi-analytical modeling of the global HI 21-cm signal from CD considering several physical processes. Further, it investigates the nature of galaxies that dominate during CD and EoR using current available observations. In our work, we study the redshift evolution of the primordial magnetic field (PMF) during the dark ages and cosmic dawn, and prospects of constraining it in light of EDGES 21-cm signal in the `colder IGM' background. We find that the IGM heating rate due to the PMF enhances compared to the standard scenario. However, PMF is an unlikely candidate for explaining the rise of the EDGES absorption signal at lower redshift. We further consider, in detail, the heating of the IGM owing to cosmic ray protons generated by the supernovae from both early Pop~III and Pop~II stars. We show that the EDGES signal can be well fitted by the cosmic ray heating along with the Lyman-$\alpha$ coupling and the dark matter-baryon interaction. We, further, explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a MCMC sampler, we find that high contribution from low-mass halos along with high photon escape fractions are required to simultaneously reproduce the existing constraints. With the extreme effort in building more advanced and sophisticated telescopes, the future 21-cm signal detection would be able to provide better constraints on the amplitude of PMF and the efficiencies on cosmic ray protons, and consequently on early star formation rates.

Age is one of the most fundamental parameters of stars, yet it is one of the hardest to determine as it requires modelling various aspects of stellar formation and evolution. When we compare the ages derived from isochronal and dynamical traceback methods for six young stellar associations, we find a systematic discrepancy. Specifically, dynamical traceback ages are consistently younger by an average of $\langle\Delta_{\rm Age}\rangle = 5.5 \pm 1.1$ Myr. We rule out measurement errors as the cause of the age mismatch and propose that $\Delta_{\rm Age}$ indicates the time a young star remains bound to its parental cloud before moving away from its siblings. In this framework, the dynamical traceback "clock" starts when a stellar cluster or association begins to expand after expelling most of the gas, while the isochronal "clock" starts earlier when most stars form. The difference between these two age-dating techniques is a powerful tool to constraint evolutionary models, as isochronal ages cannot be younger than dynamical traceback ages. Measuring the $\Delta_{\rm Age}$ accurately and understanding its variations across different environments will provide further information on the impact of local conditions and stellar feedback on the formation and dispersal of stellar clusters.

The Galactic global magnetic field is thought to play a vital role in shaping Galactic structures such as spiral arms and giant molecular clouds. However, our knowledge of magnetic field structures in the Galactic plane at different distances is limited, as measurements used to map the magnetic field are the integrated effect along the line of sight. In this study, we present the first-ever tomographic imaging of magnetic field structures in a Galactic spiral arm. Using optical stellar polarimetry over a $17' \times 10'$ field of view, we probe the Sagittarius spiral arm. Combining these data with stellar distances from the $Gaia$ mission, we can isolate the contributions of five individual clouds along the line of sight by analyzing the polarimetry data as a function of distance. The observed clouds include a foreground cloud ($d < 200$ pc) and four clouds in the Sagittarius arm at 1.23 kpc, 1.47 kpc, 1.63 kpc, and 2.23 kpc. The column densities of these clouds range from 0.5 to $2.8 \times 10^{21}~\mathrm{cm}^{-2}$. The magnetic fields associated with each cloud show smooth spatial distributions within their observed regions on scales smaller than 10 pc and display distinct orientations. The position angles projected on the plane-of-sky, measured from the Galactic north to east, for the clouds in increasing order of distance are $135^\circ$, $46^\circ$, $58^\circ$, $150^\circ$, and $40^\circ$, with uncertainties of a few degrees. Notably, these position angles deviate significantly from the direction parallel to the Galactic plane.

We present a measurement of the Hubble Constant $H_0$ using the gravitational wave event GW190412, an asymmetric binary black hole merger detected by LIGO/Virgo, as a dark standard siren. This event does not have an electromagnetic counterpart, so we use the statistical standard siren method and marginalize over potential host galaxies from the Dark Energy Spectroscopic Instrument (DESI) survey. GW190412 is well-localized to 12 deg$^2$ (90% credible interval), so it is promising for a dark siren analysis. The dark siren value for $H_0=85.4_{-33.9}^{+29.1}$ km/s/Mpc, with a posterior shape that is consistent with redshift overdensities. When combined with the bright standard siren measurement from GW170817 we recover $H_0=77.96_{-5.03}^{+23.0}$ km/s/Mpc, consistent with both early and late-time Universe measurements of $H_0$. This work represents the first standard siren analysis performed with DESI data, and includes the most complete spectroscopic sample used in a dark siren analysis to date.

The bispectrum is an important statistics helpful for measuring the primordial non-Gaussianity parameter $f_{\mathrm{NL}}$ to less than order unity in error, which would allow us to distinguish between single and multi-field inflation models. The Spectro-Photometer for the History of the Universe, Epoch of Reionization and Ices Explorer (SPHEREx) mission is particularly well-suited for making this measurement with its $\sim$100-band all-sky observations in the near-infrared. Consequently, the SPHEREx data will contain galaxies with spectroscopic-like redshift measurements as well as those with much larger errors. In this paper, we evaluate the impact of photometric redshift errors on $f_{\mathrm{NL}}$ constraints in the context of an updated multi-tracer forecast for SPHEREx, finding that the azimuthal averages of the first three even bispectrum multipoles are no longer sufficient for capturing most of the information (as opposed to the case of spectroscopic surveys shown in the literature). The final SPHEREx result with all five galaxy samples and six redshift bins is however not severely impacted because the total result is dominated by the samples with the best redshift errors, while the worse samples serve to reduce cosmic variance. Our fiducial result of $\sigma_{f_{\mathrm{NL}}} = 0.7$ from bispectrum alone is increased by $18\%$ and $3\%$ when using $l_{\mathrm{max}}=0$ and 2 respectively. We also explore the impact on parameter constraints when varying the fiducial redshift errors, as well as using subsets of multi-tracer combinations or triangles with different squeezing factors. Note that the fiducial result here is not the final SPHEREx capability, which is still on target for being $\sigma_{f_{\mathrm{NL}}} = 0.5$ once the power spectrum will be included.

Outflows in the Active Galactic Nucleus (AGN) are considered to play a key role in the host galaxy evolution through transfer of a large amount of energy. A Narrow Line Region (NLR) in the AGN is composed of ionized gas extending from pc-scales to kpc-scales. It has been suggested that shocks are required for ionization of the NLR gas. If AGN outflows generate such shocks, they will sweep through the NLR and the outflow energy will be transferred into a galaxy-scale region. In order to study contribution of the AGN outflow to the NLR-scale shock, we measure the [\ion{Fe}{2}]$\lambda12570$/[\ion{P}{2}]$\lambda11886$ line ratio, which is a good tracer of shocks, using near-infrared spectroscopic observations with WINERED (Warm INfrared Echelle spectrograph to Realize Extreme Dispersion and sensitivity) mounted on the New Technology Telescope. Among 13 Seyfert galaxies we observed, the [\ion{Fe}{2}] and [\ion{P}{2}] lines were detected in 12 and 6 targets, respectively. The [\ion{Fe}{2}]/[\ion{P}{2}] ratios in 4 targets were found to be higher than 10, which implies the existence of shocks. We also found that the shock is likely to exist where an ionized outflow, i.e., a blue wing in [\ion{S}{3}]$\lambda9533$, is present. Our result implies that the ionized outflow present over a NLR-scale region sweeps through the interstellar medium and generates a shock.

The current studies of microlensing planets are limited by small number statistics. Follow-up observations of high-magnification microlensing events can efficiently form a statistical planetary sample. Since 2020, the Korea Microlensing Telescope Network (KMTNet) and the Las Cumbres Observatory (LCO) global network have been conducting a follow-up program for high-magnification KMTNet events. Here, we report the detection and analysis of a microlensing planetary event, KMT-2023-BLG-1431, for which the subtle (0.05 magnitude) and short-lived (5 hours) planetary signature was characterized by the follow-up from KMTNet and LCO. A binary-lens single-source (2L1S) analysis reveals a planet/host mass ratio of $q = (0.72 \pm 0.07) \times 10^{-4}$, and the single-lens binary-source (1L2S) model is excluded by $\Delta\chi^2 = 80$. A Bayesian analysis using a Galactic model yields estimates of the host star mass of $M_{\rm host} = 0.57^{+0.33}_{-0.29}~M_\odot$, the planetary mass of $M_{\rm planet} = 13.5_{-6.8}^{+8.1}~M_{\oplus}$, and the lens distance of $D_{\rm L} = 6.9_{-1.7}^{+0.8}$ kpc. The projected planet-host separation of $a_\perp = 2.3_{-0.5}^{+0.5}$ au or $a_\perp = 3.2_{-0.8}^{+0.7}$, subject to the close/wide degeneracy. We also find that without the follow-up data, the survey-only data cannot break the degeneracy of central/resonant caustics and the degeneracy of 2L1S/1L2S models, showing the importance of follow-up observations for current microlensing surveys.

Space safety and astronomy are at odds. The problem posed by space debris and derelict satellites in the low Earth orbit is an existential threat to all space operations. These dangerous objects in space are more easily tracked with ground-based LiDAR if they are highly reflective, especially in the near-infrared (NIR) range. At the same time, reflective objects in orbit are the bane of ground-based astronomers, causing light pollution and marring images with bright streaks. How can this tension be resolved? The hypothesis tested is that a near-infrared-transparent (NIRT) coating which is opaque in the visible light range and transparent in the NIR range is a promising candidate for use in satellite construction. This experiment tests whether typical spacecraft surfaces such as anodized aluminum or multi-layer insulation (MLI) with a NIRT coating applied will absorb visible light and reflect NIR. The findings confirm the efficacy of the NIRT coating for this purpose, reducing visible light reflection by 47% (+/-3%) and increasing reflection in the NIR by 7% (+/-2%). This promising novel NIRT coating may help provide a path forward to resolve the tension between astronomy and the space industry.

Fast radio bursts (FRBs) can present narrow-band spectra with a variety of polarization properties. We study spectral properties from perspectives of intrinsic radiation mechanisms by invoking bunching mechanism and perturbations on charged bunches moving at curved trajectories. The narrow-band emission likely reflects some quasi-periodic structure on the bulk of bunches, which may be due to quasi-periodically sparking in a "gap" or quasi-monochromatic Langmuir waves. A spike would appears at the spectrum if the perturbations can induce a monochromatic oscillation of bunches, however, it is hardly to create a narrow-band spectrum because a large amplitude scattering of Lorentz factor can let the spike wider. Monochromatic Alfv\'en waves as the perturbations do not generate the narrow-band spectrum from 100 MHz to 10 GHz. Both the bunching mechanism and perturbations scenarios share the same polarization properties with a normally distributed bulk of bunches. We investigate absorption effects including Landau damping and curvature self-absorption in magnetosphere, which are significant at low frequencies. Subluminous O-mode photons can not escape from the magnetosphere, leading to a height-dependent lower frequency cut-off, so that the spectrum is narrow-band when the cut-off is close to the characteristic frequency of curvature radiation. The spectral index is 5/3 at low frequency bands due to the curvature self-absorption but not as steep as the observations.

It is widely accepted that X-ray emission in luminous active galactic nuclei (AGNs) originates from hot corona. To prevent the corona from over-cooling by strong X-ray emission, steady heating to the corona is essential, for which the most promising mechanisms is the magnetic reconnection. Detailed studies of the coupled disc and corona, in the frame of magnetic field transferring accretion-released energy from the disc to the corona, reveal that the thermal electrons can only produce X-ray spectrum with $\Gamma_{\rm 2-10\,keV}>2.1$, which is an inevitable consequence of the radiative coupling of the thermal corona and disc. In the present work, we develop the magnetic-reconnection-heated corona model by taking into account the potential non-thermal electrons accelerated in the magnetic reconnection process, in addition to the thermal electrons. We show that the features of the structure and spectrum of the coupled disc and corona can be affected by the fraction of magnetic energy allocated to thermal electrons. Furthermore, we investigate the effects of the power-law index and energy range of non-thermal electrons and the magnetic field on the spectrum. It is found that the X-ray spectrum from the Comptonization of the hybrid electrons can be flatter than that from thermal electrons only, in agreement with observations. By comparing with the observed hard X-ray data, we suggest that a large fraction ($>40\%$) of the magnetic energy be allocated to the non-thermal electrons in the luminous and flat X-ray spectrum AGNs.

Cosmic shear and cosmic magnification reflect the same gravitational lensing field. Each of these two probes are affected by different systematics. We study the auto- and cross-correlations of the cosmic shear from the China Space Survey Telescope (CSST) and cosmic magnification of supernovae from Large Synoptic Survey Telescope (LSST). We want to answer, to what extent, by adding the magnification data we can remove the systematic bias in cosmic shear measurement. We generate the mock shear/magnification maps based on the correlation between of different tomographic bins. After obtaining the corrected power spectra, we adopt the Markov Chain Monte Carlo (MCMC) technique to fit the theoretical models, and investigate the constraints on the cosmological and nuisance parameters. We find that the with only cosmic shear data, there are $1\sigma$ bias in $\sigma_8$ and intrinsic alignment model parameters. By adding the magnification data, we are able to remove these biases perfectly.

We have been performing daily observations of bright microquasars at 1.2-20~GHz with the Northern sector of RATAN-600 radio telescope for more than ten years. During the 2019--2021 observations we recorded bright flares, which we call giant flares because fluxes reach record levels -- above 20~Jy -- during these events. In this paper we report the results of intraday variations of the Cygnus~X-3} microquasar in multi-azimuth observations made with the "North sector with a flat-sheet reflector" during giant flares of Cygnus X-3. These were the first such observations made simultaneously at several frequencies on a short time scale (10 minutes). Observational data consists of 31 measurement made within +/-2.7 hours of the culmination of the object. We are the first to discover the evolution of the spectrum of the flare emission of Cygnus~X-3 on a time scale comparable to the orbital period of the binary. The measurement data allowed us to determine the temporal and spectral parameters of radio emission, which are typical for synchrotron flare emission in relativistic jets. Evolution of the radio emission of X-ray binaries on short time scales is a key to understanding the formation of jet outbursts in the process of mass accretion of the matter of the donor star onto the relativistic object.

The stellar mass Tully-Fisher relation (STFR) and its scatter encode valuable information about the processes shaping galaxy evolution across cosmic time. However, we are still missing a proper quantification of the STFR slope and scatter dependence on the baryonic tracer used to quantify rotational velocity, on the velocity measurement radius and on galaxy integrated properties. We present a catalogue of stellar and ionised gas (traced by H$\alpha$ emission) kinematic measurements for a sample of galaxies drawn from the MaNGA Galaxy Survey, providing an ideal tool for galaxy formation model calibration and for comparison with high-redshift studies. We compute the STFRs for stellar and gas rotation at 1, 1.3 and 2 effective radii ($R_e$). The relations for both baryonic components become shallower at 2$R_e$ compared to 1$R_e$ and 1.3$R_e$. We report a steeper STFR for the stars in the inner parts ($\leq 1.3 R_e$) compared to the gas. At 2$R_e$, the relations for the two components are consistent. When accounting for covariances with integrated v/$\sigma$, scatter in the stellar and gas STFRs shows no strong correlation with: optical morphology, star formation rate surface density, tidal interaction strength or gas accretion signatures. Our results suggest that the STFR scatter is driven by an increase in stellar/gas dispersional support, from either external (mergers) or internal (feedback) processes. No correlation between STFR scatter and environment is found. Nearby Universe galaxies have their stars and gas in statistically different states of dynamical equilibrium in the inner parts ($\leq 1.3 R_e$), while at 2$R_{e}$ the two components are dynamically coupled.

Radiative transfer (RT) simulations are a powerful tool that enables the calculation of synthetic images of a wide range of astrophysical objects. These simulations are often based on the Monte Carlo (MC) method, as it provides the needed versatility that allows the consideration of the diverse and often complex conditions found in those objects. However, this method faces fundamental problems in the regime of high optical depths which may result in noisy images and underestimated flux values. In this study, we propose an advanced MCRT method, i.e., an enforced minimum scattering order that is aimed at providing a minimum quality of determined flux estimates. For that purpose, we extended our investigations of the scattering order problem and derived an analytic expression for the minimum number of interactions that depends on the albedo and optical depth of the system, which needs to be considered to achieve a certain coverage of the scattering order distribution. The method is based on the utilization of this estimated minimum scattering order and enforces the consideration of a sufficient number of interactions during a simulation. Moreover, we identified two notably distinct cases that shape the kind of complexity that arises in MCRT simulations: the albedo-dominated and the optical depth-dominated case. Based on that, we analyzed implications regarding the best usage of a stretching method as a means to alleviate the scattering order problem. We find that its most suitable application requires taking into account the albedo and the optical depth. Then, we argue that the derived minimum scattering order can be used to assess the performance of a stretching method with regard to the scattering orders its usage promotes. Finally, we stress the need for developing advanced pathfinding techniques to fully solve the problem of MCRT simulations in the regime of high optical depths.

The combination of high-contrast imaging and medium to high spectral resolution spectroscopy offers new possibilities for the detection and characterization of exoplanets. The molecular mapping technique uses the difference between the planetary and stellar spectra. While traditional post-processing techniques are quickly limited by speckle noise at short angular separation, it efficiently suppresses speckles. Its performance depends on multiple parameters such as the star magnitude, the adaptive optics residual halo, the companion spectrum, the telluric absorption, as well as the telescope and instrument properties. Exploring this parameter space through end-to-end simulations to predict potential science cases and to optimize future instrument designs is very time-consuming, making it difficult to draw conclusions. We propose to define an efficient methodology for such an analysis. Explicit expressions of the estimates of signal and noise are derived, and they are validated through comparisons with end-to-end simulations. They provide an understanding of the instrumental dependencies, and help to discuss optimal instrumental choices with regard to the targets of interest. They are applied in the case of ELT/HARMONI, as a tool to predict the contrast performance in various observational cases. We confirm the potential of molecular mapping for high-contrast detections, especially for cool planets at short separations. We provide guidelines based on quantified estimates for design trade-offs of future instruments. We discuss the planet detection performances of HARMONI observing modes. While they nicely cover the appropriate requirements for high detection capability of warm exoplanets, a transmission extended down to J band would be beneficial. A contrast of a few 1E-7 at 50mas should be within reach on bright targets in photon noise regime with molecular mapping.

Magnetic fields have been measured recently in the core of red giant stars thanks to their effects on stellar oscillation frequencies. The search for magnetic signatures in pulsating stars, such as $\gamma$ Doradus or Slowly Pulsation B stars, requires to adapt the formalism developed for the slowly rotating red giants to rapidly rotating stars. We perform a theoretical analysis of the effects of an arbitrary magnetic field on high radial order gravity and Rossby modes in a rapidly rotating star. The magnetic effects are treated as a perturbation. For high radial order modes, the contribution of the radial component of the magnetic field is likely to dominate over the azimuthal and latitudinal components. The rotation is taken into account through the traditional approximation of rotation. General expressions of the frequency shift induced by an arbitrary radial magnetic field are derived. Approximate analytical forms are obtained in the high-order high-spin-parameter limits for the modes most frequently observed in $\gamma$ Dor stars. We propose simple methods to detect seismic magnetic signatures and measure possible magnetic fields in such stars. These methods offer new possibilities to look for internal magnetic fields in future observations, such as the ones of the PLATO mission, or to revisit existing Kepler or TESS data.

The investigation of emission line regions within active galaxies (AGNs) has a rich and extensive history, now extending to the use of AGNs and quasars as "standardizable" cosmological indicators, shedding light on the evolution of our universe. As we enter the era of advanced observatories, such as the successful launch of JWST and the forthcoming Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), the landscape of AGN exploration across cosmic epochs is poised for exciting advancements. In this work, we delve into recent developments in AGN variability research, anticipating the substantial influx of data facilitated by LSST. The article highlights recent strides made by the AGN Polish Consortium in their contributions to LSST. The piece emphasizes the role of quasars in cosmology, dissecting the intricacies of their calibration as standard candles. The primary focus centers on the relationship between the broad-line region size and luminosity, showcasing recent breakthroughs that enhance our comprehension of this correlation. These breakthroughs encompass a range of perspectives, including spectroscopic analyses, photoionization modeling, and collaborative investigations with other cosmological tools. The study further touches on select studies, underlining how the synergy of theoretical insights and advancements in observational capabilities has yielded deeper insights into these captivating cosmic entities.

We study the long-term spin period evolution of the Be/X-ray binary pulsar GX 304-1 and discover a new torque reversal after nearly three years of spinning up. The averaged spin-up rate before the torque reversal is $\sim 1.3 \times 10^{-13}$ Hz s$^{-1}$ which changes to an averaged spin-down rate of $\sim -5 \times 10^{-14}$ Hz s$^{-1}$. The pulsar is detected at low luminosities (about ${2}\times 10^{35}$ erg s$^{-1}$) near periastron passages during the torque reversal suggesting that accretion is not quenched during this period. The long-term optical observations of the companion star suggest that the activity of the companion star may have decreased compared to the period when X-ray outbursts were earlier detected from the pulsar. The spin-up rates estimated during regular bright outbursts of the pulsar are observed to decrease systematically as the pulsar enters a low activity state after the outbursts and undergoes torque reversal. We explore plausible mechanisms to explain the torque reversal and long-term spin-down in this pulsar.

The Hubble parameter $H_0$, is not a univocally-defined quantity: it relates redshifts to distances in the near Universe, but is also a key parameter of the $\Lambda$CDM standard cosmological model. As such, $H_0$ affects several physical processes at different cosmic epochs, and multiple observables. We have counted more than a dozen $H_0$'s which are expected to agree if a) there are no significant systematics in the data and their interpretation and b) the adopted cosmological model is correct. With few exceptions (proverbially confirming the rule) these determinations do not agree at high statistical significance; their values cluster around two camps: the low (68 km/s/Mpc) and high (73 km/s/Mpc) camp. It appears to be a matter of anchors: the shape of the Universe expansion history agrees with the model, it is the normalizations that disagree. Beyond systematics in the data/analysis, if the model is incorrect there are only two viable ways to "fix" it: by changing the early time ($z\gtrsim 1100$) physics and thus the early time normalization, or by a global modification, possibly touching the model's fundamental assumptions (e.g., homogeneity, isotropy, gravity). None of these three options has the consensus of the community. The research community has been actively looking for deviations from $\Lambda$CDM for two decades; the one we might have found makes us wish we could put the genie back in the bottle.

Dark matter may be composed of ultra-light bosons whose de Broglie wavelength in galaxies is of order 1 kpc. The standard model for this fuzzy dark matter (FDM) is a complex scalar field that obeys the Schr\"odinger-Poisson equations. The wavelike nature of FDM leads to fluctuations in the gravitational field that can pump energy into the stellar components of a galaxy. Heuristic arguments and theoretical analyses suggest that these fluctuations can be modelled by replacing FDM with a system of quasiparticles (QPs). We test this hypothesis by comparing self-consistent simulations of a Schr\"odinger field with those using a system of QPs in one spatial dimension. Simulations of pure FDM systems allow us to derive a phenomenological relation between the number of QPs that is required to model FDM with a given de Broglie wavelength. We also simulate systems of FDM and stars and find that the FDM pumps energy into the stars whether it is described by QPs or a Schr\"odinger field with the FDM adiabatically contracting and the stellar system adiabatically expanding. However, we find that QPs overestimate dynamical heating.

Accurate particle simulations are essential for the next generation of experiments in astroparticle physics. The Monte Carlo simulation library PROPOSAL is a flexible tool to efficiently propagate high-energy leptons and photons through large volumes of media, for example in the context of underground observatories. It is written as a C++ library, including a Python interface. In this paper, the most recent updates of PROPOSAL are described, including the addition of electron, positron, and photon propagation, for which new interaction types have been implemented. This allows the usage of PROPOSAL to simulate electromagnetic particle cascades, for example in the context of air shower simulations. The precision of the propagation has been improved by including rare interaction processes, new photonuclear parametrizations, deflections in stochastic interactions, and the possibility of propagating in inhomogeneous density distributions. Additional technical improvements regarding the interpolation routine and the propagation algorithm are described.

The astrochemistry of the important biogenic element phosphorus (P) is still poorly understood, but observational evidence indicates that P-bearing molecules are likely associated with shocks. We study P-bearing molecules, as well as some shock tracers, towards one of the chemically richest hot molecular core, G31.41+0.31, in the framework of the project "G31.41+0.31 Unbiased ALMA sPectral Observational Survey" (GUAPOS), observed with the Atacama Large Millimeter Array (ALMA). We have observed the molecules PN, PO, SO, SO2, SiO, and SiS, through their rotational lines in the spectral range 84.05-115.91 GHz, covered by the GUAPOS project. PN is clearly detected while PO is tentatively detected. The PN emission arises from two regions southwest of the hot core peak, "1" and "2", and is undetected or tentatively detected towards the hot core peak. the PN and SiO lines are very similar both in spatial emission morphology and spectral shape. Region "1" is in part overlapping with the hot core and it is warmer than region "2", which is well separated from the hot core and located along the outflows identified in previous studies. The column density ratio SiO/PN remains constant in regions "1" and "2", while SO/PN, SiS/PN, and SO2/PN decrease by about an order of magnitude from region "1" to region "2", indicating that SiO and PN have a common origin even in regions with different physical conditions. Our study firmly confirms previous observational evidence that PN emission is tightly associated with SiO and it is likely a product of shock-chemistry, as the lack of a clear detection of PN towards the hot-core allows to rule out relevant formation pathways in hot gas. We propose the PN emitting region "2" as a new astrophysical laboratory for shock-chemistry studies

We perform nonlocal thermodynamic equilibrium (NLTE) model atmosphere analyses of the three hottest hypervelocity stars (space velocities between $\approx$ 1500-2800 km s$^{-1}$) known to date, which were recently discovered spectroscopically and identified as runaways from Type Ia supernovae. The hottest of the three (J0546$+$0836, effective temperature $T_\mathrm{eff}$ = 95,000 $\pm$ 15,000 K, surface gravity log g = $5.5 \pm 0.5$) has an oxygen-dominated atmosphere with a significant amount of carbon (C = $0.10 \pm 0.05$, O = $0.90 \pm 0.05$, mass fractions). Its mixed absorption+emission line spectrum exhibits photospheric absorption lines from O V and O VI as well as O III and O IV emission lines that are formed in a radiation-driven wind with a mass-loss rate of the order of $10^{-8}$ $M_\odot$ yr$^{-1}$. Spectroscopically, J0546$+$0836 is a [WC]-PG1159 transition-type pre-white dwarf. The second object (J0927$-$6335) is a PG1159-type white dwarf with a pure absorption-line spectrum dominated by C III/C IV and O III/O IV. We find $T_\mathrm{eff}$ = 60,000 $\pm$ 5000 K, log g = $7.0 \pm 0.5$, and a carbon- and oxygen-dominated atmosphere with C = $0.47 \pm 0.25$, O = $0.48 \pm 0.25$, and possibly a minute amount of helium (He = $0.05 \pm 0.05$). Comparison with post-AGB evolutionary tracks suggests a mass of $M\approx0.5$ $M_\odot$ for both objects, if such tracks can safely be applied to these stars. We find the third object (J1332$-$3541) to be a relatively massive ($M=0.89 M_\odot$) hydrogen-rich (DAO) white dwarf with $T_\mathrm{eff}$ = 65,657 $\pm$ 2390 K, log g = $8.38 \pm 0.08$, and abundances H = $0.65 \pm 0.04$ and He = $0.35 \pm 0.04$. We discuss our results in the context of the "dynamically driven double-degenerate double-detonation" (D$^6$) scenario proposed for the origin of these stars.

Radio-loud active galactic nuclei (RLAGN) play an important role in the evolution of galaxies through the effects on their environment. The two major morphological classes are core-bright (FRI) and edge-bright (FRII) sources. With the LOw-Frequency ARray (LOFAR) we compare the FRI and FRII evolution down to lower flux densities and with larger samples than before with the aim to examine the cosmic space density evolution for FRIs and FRIIs by analyzing their space density evolution between L_150~10^24.5 W/Hz and L_150~10^28.5 W/Hz and up to z=2.5. We construct radio luminosity functions (RLFs) from FRI and FRII catalogues based on recent data from LOFAR at 150MHz to study the space densities as a function of radio luminosity and redshift. To partly correct for selection biases and completeness, we simulate how sources appear at a range of redshifts. We report a space density enhancement from low to high redshift for FRI and FRII sources brighter than L_150~10^27 W/Hz. This is possibly related to the higher gas availability in the earlier denser universe. The constant FRI/FRII space density ratio evolution as a function of radio luminosity and redshift in our results suggests that the jet-disruption of FRIs might be primarily caused by events occurring on scales within the host galaxy, rather than being driven by changes in the overall large-scale environment. Remaining selection biases in our results also highlight the need to resolve more sources at angular scales below 40 arcsec and therefore strengthens the motivation for the further development and automation of the calibration and imaging pipeline of LOFAR data to produce images at sub-arcsecond resolution.

One of the goals of Space Weather studies is to achieve a better understanding of impulsive phenomena, such as Coronal Mass Ejections (CMEs), in order to improve our ability to forecast them and mitigate the risk to our technologically driven society. The essential part of achieving this goal is to assess the performance of forecasting models. To this end, the quality and availability of suitable data are of paramount importance. In this work, we have merged already publicly available data of CMEs from both in-situ and remote instrumentation in order to build a database of CME properties. To evaluate the accuracy of such a database and confirm the relationship between in-situ and remote observations, we have employed the drag-based model (DBM) due to its simplicity and inexpensive cost of computational resources. In this study, we have also explored the parameter space for the drag parameter and solar wind speed using a Monte Carlo approach to evaluate how well the DBM determines the propagation of CMEs for the events in the dataset. The dataset of geoeffective CMEs constructed as a result of this work provides validation of the initial hypothesis about DBM, and solar wind speed and also yields further insight into CME features like arrival time, arrival speed, lift-off time, etc. Using a data-driven approach, this procedure allows us to present a homogeneous, reliable, and robust dataset for the investigation of CME propagation. On the other hand, possible CME events are identified where DBM approximation is not valid due to model limitations and higher uncertainties in the input parameters, those events require more thorough investigation.

Context: Coronal mass ejections (CMEs) are rapid eruptions of magnetized plasma that occur on the Sun, which are known as the main drivers of adverse space weather. Accurately tracking their evolution in the heliosphere in numerical models is of utmost importance for space weather forecasting. Aims: The main objective of this paper is to implement the Regularized Biot-Savart Laws (RBSL) method in a new global corona model COCONUT. This approach has the capability to construct the magnetic flux rope with an axis of arbitrary shape. Methods: We present the implementation process of the RBSL flux rope model in COCONUT, which is superposed onto a realistic solar wind reconstructed from the observed magnetogram around the minimum of solar activity. Based on this, we simulate the propagation of an S-shaped flux rope from the solar surface to a distance of 25 solar radii. Results: Our simulation successfully reproduces the birth process of a CME originating from a sigmoid in a self-consistent way. The model effectively captures various physical processes and retrieves the prominent features of the CMEs in observations. In addition, the simulation results indicate that the magnetic topology of the CME flux rope at around 20 solar radii deviates from a coherent structure, and manifests as a mix of open and closed field lines with diverse footpoints. Conclusions: This work demonstrates the potential of the RBSL flux rope model in reproducing CME events that are more consistent with observations. Moreover, our findings strongly suggest that magnetic reconnection during the CME propagation plays a critical role in destroying the coherent characteristic of a CME flux rope.

METIS, the Mid-Infrared ELT Imager and Spectrograph, is one of the four first-generation ELT instruments scheduled to see first light in 2028. Its two main science modules are supported by an adaptive optics system featuring a pyramid sensor with 90x90 subapertures working in the H and K bands. During the PDR and FDR phases, extensive simulations were carried out to support the sensing, reconstruction, and control concept of METIS single-conjugate adaptive optics (SCAO) system. We present details on the implementation of the COMPASS-based environment used for the simulations, the metrics used for analyzing our performance expectations, an overview of the main results, and some details on special cases like non-common path aberrations (NCPA) and water vapor seeing, as well as the low-wind effect.

The properties of sub-second time variability of the X-ray emission of the black-hole binary GRS 1915+105 are very complex and strictly connected to its patterns of variability observed on long time scales. A key aspect for determining the geometry of the accretion flow is the study of time lags between emission at different energies, as they are associated to key time scales of the system. In particular, it is important to examine the lags associated to the strong low-frequency Quasi-periodic Oscillations (QPOs), as the QPOs provide unambiguous special frequencies to sample the variability. We have analyzed data from an observation with the AstroSat satellite, in which the frequency of the low-frequency QPO varies smoothly between 2.5 and 6.6 Hz on a time scale of ~10 hours. The derived phase lags show the same properties and evolution of those observed on time scales of a few hundred days, indicating that changes in the system geometry can take place on times below one day. We fit selected energy spectra of the source and rms and phase-lag spectra of the QPO with a time-variable Comptonization model, as done previously to RossiXTE data of the same source, and find that indeed the derived parameters match those obtained for variations on much longer time scales.

We have monitored the Didymos-Dimorphos binary system in imaging polarimetric mode before and after the impact from the Double Asteroid Redirection Test (DART) mission. A previous spectropolarimetric study showed that the impact caused a dramatic drop in polarisation. Our longer-term monitoring shows that the polarisation of the post-impact system remains lower than the pre-impact system even months after the impact, suggesting that some fresh ejecta material remains in the system at the time of our observations, either in orbit or settled on the surface. The slope of the post-impact polarimetric curve is shallower than that of the pre-impact system, implying an increase in albedo of the system. This suggests that the ejected material is composed of smaller and possibly brighter particles than those present on the pre-impact surface of the asteroid. Our polarimetric maps show that the dust cloud ejected immediately after the impact polarises light in a spatially uniform manner (and at a lower level than pre-impact). Later maps exhibit a gradient in polarisation between the photocentre (which probes the asteroid surface) and the surrounding cloud and tail. The polarisation occasionally shows some small-scale variations, the source of which is not yet clear. The polarimetric phase curve of Didymos-Dimorphos resembles that of the S-type asteroid class.

Ram pressure stripping (RPS) is one of the most invoked mechanisms to explain the observed differences between cluster and field galaxies. In the local Universe, its effect on the galaxy star forming properties has been largely tackled and the general consensus is that this process first compresses the gas available in the galaxy disks, boosting the star formation for a limited amount of time, and then removes the remaining gas leading to quenching. Much less is known on the effect and preponderance of RPS at higher redshift, due to the lack of statistical samples. Exploiting VLT/MUSE observations of galaxies at 0.2<z<0.55 and the catalog of ram pressure stripped galaxies by Moretti et al., we compare the global star formation rate-mass (SFR-M*) relation of 29 cluster galaxies undergoing RPS to that of 26 field and cluster undisturbed galaxies that constitute our control sample. Stripping galaxies occupy the upper envelope of the control sample SFR-M* relation, showing a systematic enhancement of the SFR at any given mass. The boost is >3sigma when considering the SFR occurring in both the tail and disk of galaxies. The enhancement is retrieved also on local scales: considering spatially resolved data, ram pressure stripped galaxies overall have large {\Sigma}SFR values, especially for Sigma_*>10^7.5M_sun kpc~2. RPS seems to leave the same imprint on the SFR-M* and Sigma_SFR-Sigma_* relations both in the Local Universe and at z~0.35.

This study combines cosmic chronometer (CC) Hubble parameter data with growth rate (f) observations to constrain the $\Omega_{\rm m0}$ parameter. Utilizing a consistency relation independent of cosmological models, we employ Gaussian process regression to reconstruct Hubble parameter and growth rate values. The resulting $\Omega_{\rm m0}h^2$ constraint is $\Omega_{\rm m0}h^2=0.139\pm0.017$. Incorporating $H_0$ measurements, we find $\Omega_{\rm m0}$ values from CC data ($0.308\pm0.057$), tip of the Red Giant Branch ($0.285\pm0.038$), and SHOES measurements ($0.259\pm0.033$). Interestingly, a higher mean $H_0$ correlates with a lower mean $\Omega_{\rm m0}$. In summary, our cosmological model-independent approach offers valuable constraints on $\Omega_{\rm m0}$, affirming the consistency of FLRW and Newtonian perturbation theory.

Transiting planets with orbital periods longer than 40 d are extremely rare among the 5000+ planets discovered so far. The lack of discoveries of this population poses a challenge to research into planetary demographics, formation, and evolution. Here, we present the detection and characterization of HD88986b, a potentially transiting sub-Neptune, possessing the longest orbital period among known transiting small planets (< 4 R$_{\oplus}$) with a precise mass measurement ($\sigma_M/M$ > 25%). Additionally, we identified the presence of a massive companion in a wider orbit around HD88986. Our analysis reveals that HD88986b, based on two potential single transits on sector 21 and sector 48 which are both consistent with the predicted transit time from the RV model, is potentially transiting. The joint analysis of RV and photometric data show that HD88986b has a radius of 2.49$\pm$0.18 R$_{\oplus}$, a mass of 17.2$^{+4.0}_{-3.8}$ M$_{\oplus}$, and it orbits every 146.05$^{+0.43}_{-0.40}$ d around a subgiant HD88986 which is one of the closest and brightest exoplanet host stars (G2V type, R=1.543 $\pm$0.065 R$_{\odot}$, V=$6.47\pm 0.01$ mag, distance=33.37$\pm$0.04 pc). The nature of the outer, massive companion is still to be confirmed; a joint analysis of RVs, Hipparcos, and Gaia astrometric data shows that with a 3$\sigma$ confidence interval, its semi-major axis is between 16.7 and 38.8 au and its mass is between 68 and 284 M$_{Jup}$. HD88986b's wide orbit suggests the planet did not undergo significant mass loss due to extreme-ultraviolet radiation from its host star. Therefore, it probably maintained its primordial composition, allowing us to probe its formation scenario. Furthermore, the cold nature of HD88986b (460$\pm$8 K), thanks to its long orbital period, will open up exciting opportunities for future studies of cold atmosphere composition characterization.

Exploiting the Persistent Homology technique and an associated complementary representation which enables us to construct a synergistic pipeline for different topological features quantified by Betti curves in reducing the degeneracy between cosmological parameters, we investigate the footprint of summed massive neutrinos ($M_{\nu}$) in different density fields simulated by the publicly available Quijote suite. Evolution of topological features in the context of super-level filtration on three-dimensional density fields, reveals remarkable indicators for constraining the $M_{\nu}$ and $\sigma_8$. The abundance of 2-holes is more sensitive to the presence of $M_{\nu}$, also the persistence of topological features plays a crucial role in cosmological inference and reducing the degeneracy associated with $M_{\nu}$ simulation rather than their birth thresholds when either the total matter density ($m$) field or those part including only cold dark matter+baryons ($cb$) is utilized. Incorporating the Betti-1 and Betti-2 for $cb$ part of $M^+_{\nu}$ simulation marginalized over the thresholds implies $5\%$ variation compared to the massless neutrinos simulation. The constraint on $M_{\nu}$ from $\beta_k$ and its joint analysis with birth threshold and persistency of topological features for total mass density field smoothed by $R=5$ Mpc h$^{-1}$ at zero redshift reach to $0.0172$ eV and $0.0152$ eV, at $1\sigma$ confidence interval, respectively.

Context: High-resolution spectrographs fed by adaptive optics (AO) provide a unique opportunity to characterize directly imaged exoplanets. Observations with such instruments allow us to probe the atmospheric composition, spin rotation, and radial velocity of the planet, thereby helping to reveal information on its formation and migration history. The recent upgrade of the Cryogenic High-Resolution Infrared Echelle Spectrograph (CRIRES+) at the VLT makes it a highly suitable instrument for characterizing directly imaged exoplanets. Aims: In this work, we report on observations of $\beta$ Pictoris b with CRIRES+ and use them to constrain the planets atmospheric properties and update the estimation of its spin rotation. Methods: The data were reduced using the open-source \textit{pycrires} package. We subsequently forward-modeled the stellar, planetary, and systematic contribution to the data to detect molecules in the planet's atmosphere. We also used atmospheric retrievals to provide new constraints on its atmosphere. Results: We confidently detected water and carbon monoxide in the atmosphere of $\beta$ Pictoris b and retrieved a slightly sub-solar carbon-to-oxygen ratio, which is in agreement with previous results. The interpretation is hampered by our limited knowledge of the C/O ratio of the host star. We also obtained a much improved constraint on its spin rotation of $19.9 \pm 1.0$ km/s, which gives a rotation period of $8.7 \pm 0.8$ hours, assuming no obliquity. We find that there is a degeneracy between the metallicity and clouds, but this has minimal impact on the retrieved C/O, $v\sin{i}$, and radial velocity. Our results show that CRIRES+ is performing well and stands as a highly useful instrument for characterizing directly imaged planets.

We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$\Lambda$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calculations, including the omission of screening terms and the use of perturbative kernels based on the Einstein-de Sitter universe; we explore different infrared-resummation schemes, tracer bias models and a linear treatment of massive neutrinos; we employ two models for redshift space distortions, the Taruya-Nishimishi-Saito prescription and the Effective Field Theory of Large-Scale Structure. This work further provides a preliminary validation of the codes being considered by Euclid for the spectroscopic clustering probe in beyond-$\Lambda$CDM scenarios. We calculate and compare the $\chi^2$ statistic to assess the different modelling choices. This is done by fitting the spectroscopic clustering predictions to measurements from numerical simulations and perturbation theory-based mock data. We compare the behaviour of this statistic in the beyond-$\Lambda$CDM cases, as a function of the maximum scale included in the fit, to the baseline $\Lambda$CDM case. We find that the Einstein-de Sitter approximation without screening is surprisingly accurate for all cases when comparing to the halo clustering monopole and quadrupole obtained from simulations. Our results suggest that the inclusion of multiple redshift bins, higher-order multipoles, higher-order clustering statistics (such as the bispectrum) and photometric probes such as weak lensing, will be essential to extract information on massive neutrinos, modified gravity and dark energy.

PSR J1012$-$4235 is a 3.1ms pulsar in a wide binary (37.9 days) with a white dwarf companion. We detect, for the first time, a strong relativistic Shapiro delay signature in PSR J1012$-$4235. Our detection is the result of a timing analysis of data spanning 13 years and collected with the Green Bank, Parkes, and MeerKAT Radio Telescopes and the Fermi $\gamma$-ray space telescope. We measured the orthometric parameters for Shapiro delay and obtained a 22$\sigma$ detection of the $h_{\rm 3}$ parameter of 1.222(54) $\mu$s and a 200$\sigma$ detection of $\varsigma$ of 0.9646(49). With the assumption of general relativity, these measurements constrain the pulsar mass ($M_{\rm p}=1.44^{+0.13}_{-0.12}$M$_{\odot}$), the mass of the white dwarf companion ($M_{\rm c} = 0.270^{+0.016}_{-0.015}$M$_{\odot}$ ), and the orbital inclination ($i=88.06^{+0.28}_{-0.25} \deg$). Including the early $\gamma$-ray data in our timing analysis facilitated a precise measurement of the proper motion of the system of 6.58(5) mas yr$^{-1}$. We also show that the system has unusually small kinematic corrections to the measurement of the orbital period derivative, and therefore has the potential to yield stringent constraints on the variation of the gravitational constant in the future.

We explore fitting gamma-ray burst spectra with three physically-motivated models, and thus revisit the viability of synchrotron radiation as the primary source of GRB prompt emission. We pick a sample of 100 bright GRBs observed by the Fermi Gamma-ray Burst Monitor (GBM), based on their energy flux values. In addition to the standard empirical spectral models used in previous GBM spectroscopy catalogs, we also consider three physically-motivated models; (a) a Thermal Synchrotron model, (b) a Band model with a High-energy Cutoff, and (c) a Smoothly Broken Power Law (SBPL) model with a Multiplicative Broken Power Law (MBPL). We then adopt the Bayesian information criterion (BIC) to compare the fits obtained and choose the best model. We find that 42% of the GRBs from the fluence spectra and 23% of GRBs from the peak-flux spectra have one of the three physically-motivated models as their preferred one. From the peak-flux spectral fits, we find that the low-energy index distributions from the empirical model fits for long GRBs peak around the synchrotron value of -2/3, while the two low-energy indices from the SBPL+MBPL fits of long GRBs peak close to the -2/3 and -3/2 values expected for a synchrotron spectrum below and above the cooling frequency.

We present a photometric characterization of 208 ultra-cool dwarfs (UCDs) with spectral types between M4 and L4, from 20-second and 2-minute cadence TESS light curves. We determine rotation periods for 87 objects (42 percent) and identify 778 flare events in 103 UCDs (49.5 percent). For 777 flaring events (corresponding to 102 objects), we derive bolometric energies between 2.1e30 and 1.1e34 erg , with 56 superflare events. No transiting planets or eclipsing binaries were identified. We find that the fraction of UCDs with rotation and flaring activity is, at least, 20 percent higher in M4-M6 spectral types than in later UCDs (M7-L4). For spectral types between M4 and L0, we measure the slope of the flare bolometric energy-duration correlation to be gamma = 0.497 +/- 0.058, which agrees with that found in previous studies for solar-type and M dwarfs. Moreover, we determine the slope of the flare frequency distribution to be alpha = -1.75 +/- 0.04 for M4-M5 dwarfs, alpha = -1.69 +/- 0.04 and alpha = -1.72 +/- 0.1 for M6-M7 and M8-L0 dwarfs, respectively, which are consistent with previous works that exclusively analysed UCDs. These results support the idea that independently of the physical mechanisms that produce magnetic activity, the characteristics of the rotational modulation and flares are similar for both fully-convective UCDs and partially-convective solar-type and early-M stars. Based on the measured UCD flare distributions, we find that UV radiation emitted from flares does not have the potential to start prebiotic chemistry.

We propose a method to determine the mass and spin of primordial black holes based on measuring the energy and emission rate at the dipolar and quadrupolar peaks in the primary photon Hawking spectrum, applicable for dimensionless spin parameters $\tilde{a}\gtrsim 0.6$. In particular, we show that the ratio between the energies of the two peaks is only a function of the black hole spin, while the ratio between their emission rates depends also on the line-of-sight inclination. The black hole mass and distance from the Earth may then be inferred from the absolute values of the peak energies and emission rates. This method is relevant for primordial black holes born with large spin parameters that are presently still in the early stages of their evaporation process.

The early post-merger phase of a binary neutron-star coalescence is shaped by characteristic rotational velocities as well as violent density oscillations and offers the possibility to constrain the properties of neutron star matter by observing the gravitational wave emission. One possibility to do so is the investigation of gravitational wave damping through the bulk viscosity which originates from violations of weak chemical equilibrium. Motivated by these prospects, we present a comprehensive report about the implementation of the self-consistent and second-order formulation of the equations of relativistic hydrodynamics for dissipative fluids proposed by M\"uller, Israel and Stewart. Furthermore, we report on the results of two test problems, namely the viscous damping of linear density oscillations of isolated nonrotating neutron stars and the viscous migration test, both of which confirm our implementation and can be used for future code tests. Finally, we present fully general-relativistic simulations of viscous binary neutron-star mergers. We explore the structural and thermal properties of binary neutron-star mergers with a constant bulk viscosity prescription and investigate the impact of bulk viscosity on dynamical mass ejection. We find that inverse Reynolds numbers of order $\sim 1\%$ can be achieved for the highest employed viscosity thereby suppressing the dynamically ejected mass by a factor of $\sim 5$ compared to the inviscid case.

The speed of sound squared (SSS) $s^2$ in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of $s^2$ has been predicted for massive NSs using various NS Equation of State (EOS) models. However, the nature, cause, location and size of the peak in $s^2$ profiles are still very EOS model dependent. In this work, we investigate systematically $s^2$ profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman--Oppenheimer--Volkoff (TOV) equations and analyzing $s^2$ profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no $s^2$ peak to large ones for massive NSs possessing a peak in their $s^2$ profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. These features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs.

Accurate noise estimation from gravitational wave (GW) data is critical for Bayesian inference. However, recent studies on ringdown signal, such as those by Isi et al. [1], Cotesta et al. [2], and Isi and Farr [3], have encountered disagreement in noise estimation, leading to inconsistent results. The key discrepancy between these studies lies in the usage of different noise estimation methods, augmented by the usage of different sampling rates. We achieved consistent results across various sampling rates by correctly managing noise estimation, shown in the case study of the GW150914 ringdown signal. By conducting a time-domain Bayesian inference analysis on GW data, starting from the peak of the signal, we discovered that the first overtone mode is weakly supported by the amplitude distribution, with a confidence level of $1.6{\sigma}$, and is slightly disfavored by the log-Bayes factor. Overall, in our time-domain analysis we conclude there is no strong evidence for overtones in GW150914.

Neutron stars (NSs) serve as laboratories for probing strongly interacting matter at the most extreme densities. Their inner cores are expected to be dense enough to host deconfined quark matter. Utilizing state-of-the-art theoretical and multi-messenger constraints, we statistically determine the bulk properties of dense NS matter. We show that the speed of sound can be expressed in terms of the slope and curvature of the energy per particle. We demonstrate that the restoration of conformal symmetry requires changing the sign of the curvature of the bulk energy per particle as a function of energy density. Furthermore, we find that such a sign change is closely related to the peak in the speed of sound. We argue that the curvature of the energy per particle may serve as an approximate order parameter that signifies the onset of strongly coupled conformal matter in the NS core.