Papers for Wednesday, May 17 2023

M dwarfs are the most abundant stars in the Universe and are hosts of a rich diversity of planetary companions. In many cases, planets orbiting M dwarfs can be described in remarkable detail. What makes the difference is how deeply we can characterise the host star. This includes to properly model their atmospheres, their abundance of metals, and their activity processes. If they are well described individually, these numerous stars have the potential for providing statistically robust conclusions when combined into larger samples. Carmencita is the input catalogue of nearby M dwarfs for the CARMENES project, which aims to search for potentially habitable Earth-sized planets orbiting them. It contains more than two thousand M dwarfs that are scrutinized by the consortium members from multiple angles. This thesis contributes to the description of each one of these M dwarfs, including astrometry, photometry, activity, kinematics, and multiplicity, but also to the study of the sample as a whole. The empirical observations presented in this study provide an important benchmark for testing and improving theoretical predictions. By taking a careful, individualized approach to the study of M dwarfs, we not only contribute to the study of the Universe's physical processes, but we also pave the way for future discoveries of the potential for life beyond our own planet. Overall, the findings of this study underscore the importance of continued research into the most numerous stars and their planetary systems. We expect that the wealth of data gathered in this thesis will serve as a valuable resource for astronomers and researchers in related fields, and that it will inspire further investigations and new insights into the processes that shape the Universe.

Substantial effort has been devoted to the characterization of transient phenomena from photometric information. Automated approaches to this problem have taken advantage of complete phase-coverage of an event, limiting their use for triggering rapid follow-up of ongoing phenomena. In this work, we introduce a neural network with a single recurrent layer designed explicitly for early photometric classification of supernovae. Our algorithm leverages transfer learning to account for model misspecification, host galaxy photometry to solve the data scarcity problem soon after discovery, and a custom weighted loss to prioritize accurate early classification. We first train our algorithm using state-of-the-art transient and host galaxy simulations, then adapt its weights and validate it on the spectroscopically-confirmed SNe~Ia, SNe~II, and SNe~Ib/c from the Zwicky Transient Facility Bright Transient Survey. On observed data, our method achieves an overall accuracy of $82 \pm 2$% within 3 days of an event's discovery, and an accuracy of $87 \pm 5$% within 30 days of discovery. At both early and late phases, our method achieves comparable or superior results to the leading classification algorithms with a simpler network architecture. These results help pave the way for rapid photometric and spectroscopic follow-up of scientifically-valuable transients discovered in massive synoptic surveys.

In this work, we show that the Cold New Early Dark Energy (Cold NEDE) model in its original form can solve both the Hubble tension and the $S_8$ tension without adding any new ingredients at the fundamental level. So far, it was assumed that the trigger field in the Cold NEDE model is completely subdominant. However, relaxing this assumption and letting the trigger field contribute a mere $0.5\%$ of the total energy density leads to a resolution of the $S_8$ tension while simultaneously improving it as a solution to the $H_0$ tension. Fitting this model to baryonic acoustic oscillations, large-scale-structure, supernovae (including a SH0ES prior), and cosmic microwave background data, we report a preferred NEDE fraction of $f_\mathrm{NEDE}= 0.134^{+0.032}_{-0.025}$ ($68\%$ C.L.), lifting its Gaussian evidence for the first time above $5\sigma$ (up from $4 \sigma$ when the trigger contribution to dark matter is negligible). At the same time, we find the new concordance values $H_0 = 71.71 \pm 0.88 \,\mathrm{km}\, \mathrm{sec}^{-1}\, \mathrm{Mpc}^{-1}$ and $S_8 = 0.793 \pm 0.018$. Excluding large-scale structure data and the SH$_0$ES prior, both Gaussian tensions are reduced below the $2 \sigma$ level.

The full-shape correlations of the Lyman alpha (Ly$\alpha$) forest contain a wealth of cosmological information through the Alcock-Paczy\'{n}ski effect. However, these measurements are challenging to model without robustly testing and verifying the theoretical framework used for analyzing them. Here, we leverage the accuracy and volume of the $N$-body simulation suite \textsc{AbacusSummit} to generate high-resolution Ly$\alpha$ skewers and quasi-stellar object (QSO) catalogs. One of the main goals of our mocks is to aid in the full-shape Ly$\alpha$ analysis planned by the Dark Energy Spectroscopic Instrument (DESI) team. We provide optical depth skewers for six of the fiducial cosmology base-resolution simulations ($L_{\rm box} = 2\,h^{-1}{\rm Gpc}$, $N = 6912^3$) at $z = 2.5$. We adopt a simple recipe based on the Fluctuating Gunn-Peterson Approximation (FGPA) for constructing these skewers from the matter density in an $N$-body simulation and calibrate it against the 1D and 3D Ly$\alpha$ power spectra extracted from the hydrodynamical simulation IllustrisTNG (TNG; $L_{\rm box} = 205\,h^{-1}{\rm Mpc}$, $N = 2500^3$). As an important application, we study the non-linear broadening of the baryon acoustic oscillation (BAO) peak and show the cross-correlation between DESI-like QSOs and our Ly$\alpha$ forest skewers. We find differences on small scales between the Kaiser approximation prediction and our mock measurements of the Ly$\alpha$$\times$QSO cross-correlation, which would be important to account for in upcoming analyses. The \textsc{AbacusSummit} Ly$\alpha$ forest mocks open up the possibility for improved modelling of cross correlations between Ly$\alpha$ and cosmic microwave background (CMB) lensing and Ly$\alpha$ and QSOs, and for forecasts of the 3-point Ly$\alpha$ correlation function. Our catalogues and skewers are publicly available on Globus.

We develop an optimization-based maximum likelihood approach to analyze the cross-correlation of the Cosmic Microwave Background (CMB) and large-scale structure induced by the kinetic Sunyaev-Zeldovich (kSZ) effect. Our main goal is to reconstruct the radial velocity field of the universe. While the existing quadratic estimator (QE) is statistically optimal for current and near-term experiments, the likelihood can extract more signal-to-noise in the future. Our likelihood formulation has further advantages over the QE, such as the possibility of jointly fitting cosmological and astrophysical parameters and the possibility of unifying several different kSZ analyses. We implement an auto-differentiable likelihood pipeline in JAX, which is computationally tractable for a realistic survey size and resolution, and evaluate it on the Agora simulation. We also implement a machine learning-based estimate of the electron density given an observed galaxy distribution, which can increase the signal-to-noise for both the QE and the likelihood method.

Tidal disruptions of stars by stellar-mass black holes are expected to occur frequently in dense star clusters. Building upon previous studies that performed hydrodynamic simulations of these encounters, we explore the formation and long-term evolution of the thick, super-Eddington accretion disks formed. We build a disk model that includes fallback of material from the tidal disruption, accretion onto the black hole, and disk mass losses through winds launched in association with the super-Eddington flow. We demonstrate that bright transients are expected when radiation from the central engine powered by accretion onto the black hole is reprocessed at large radii by the optically-thick disk wind. By combining hydrodynamic simulations of these disruption events with our disk+wind model, we compute light curves of these wind-reprocessed transients for a wide range of stellar masses and encounter penetration depths. We find typical peak bolometric luminosities of roughly $10^{41}-10^{44}\,$erg/s (depending mostly on accretion physics parameters) and temperatures of roughly $10^5-10^6\,$K, suggesting peak emission in the ultraviolet/blue bands. We predict all-sky surveys such as the Vera Rubin Observatory and ULTRASAT will detect up to thousands of these events per year in dense star clusters out to distances of several Gpc.

Neutron stars are born out of core-collapse supernovae, and they are imparted natal kicks at birth as a consequence of asymmetric ejection of matter and possibly neutrinos. Unless the force resulting from the kicks is exerted exactly at their center, it will also cause the neutron star to rotate. In this paper, we discuss the possibility that neutron stars may receive off-center natal kicks at birth, which imprint a natal rotation. In this scenario, the observed pulsar spin and transverse velocity in the Galaxy are expected to correlate. We develop a model of the natal rotation imparted to neutron stars and constrain it by the observed population of pulsars in our Galaxy. At $90\%$ confidence, we find that the location of the off-center kick is $R_{\rm kick}=1.12^{+4.79}_{-0.97}$\,km. Our result is robust when considering pulsars with different observed periods, transverse velocities, and ages. Our constraint can be used as a guide for core-collapse simulations of massive stars.

We present 20,567 Ly$\alpha$ emitters (LAEs) at $z=2.2-7.3$ that are photometrically identified by the SILVERRUSH program in a large survey area up to 25 deg$^2$ with deep images of five broadband filters (grizy) and seven narrowband filters targeting Ly$\alpha$ lines at $z=2.2$, $3.3$, $4.9$, $5.7$, $6.6$, $7.0$, and $7.3$ taken by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey. We select secure $>5\sigma$ sources showing narrowband color excesses via Ly$\alpha$ break screening, taking into account the spatial inhomogeneity of limiting magnitudes. After removing spurious sources by careful masking and visual inspection of coadded and multi-epoch images obtained over the seven years of the surveys, we construct LAE samples consisting of 6,995, 4,642, 726, 6,124, 2,058, 18, and 5 LAEs at $z=2.2$, 3.3, 4.9, 5.7, 6.6, 7.0, and 7.3, respectively, although the $z=7.3$ candidates are tentative. Our LAE catalogs contain 241 spectroscopically confirmed LAEs at the expected redshifts from previous work. We demonstrate that the number counts of our LAEs are consistent with previous studies with similar LAE selection criteria. The LAE catalogs will be made public on our project webpage with detailed descriptions of the content and ancillary information about the masks and limiting magnitudes.

We study the geometry of the transition layer (TL) between the classical Keplerian accretion disk (the TL outer boundary) and the compact object at the TL inner boundary. Our goal is to demonstrate using the hydrodynamical formalism that the TL is created along with a shock due to a discontinuity and to an adjustment of the Keplerian disk motion to a central object. We apply hydrodynamical equations to describe a plasma motion near a central object in the TL. We point out that before matter accretes to a central object the TL cloud is formed between an adjustment radius and the TL inner boundary which is probably a site where the emergent Compton spectrum comes from. Using a generalization of the Randkine-Hugoniot relation and a solution of the azimutal force balance equation we have reproduced the geometric characteristics of TL.

It has been pointed out that the spurious Cosmic Microwave Background (CMB) B-mode polarization signals caused by the absorption of the CMB monopole component due to the Galactic interstellar matter, called the CMB shadow, degrade the accuracy of detecting the CMB B-mode polarization signals imprinted by primordial gravitational waves. We have made a realistic estimation using simulated sky maps of how the CMB shadow affects forthcoming high-precision CMB B-mode experiments for the first time. The Delta-map method, an internal template method taking into account the first-order spatial variation of foregrounds' spectral parameters, is applied as a foreground removal method. We show that if the CMB shadow effects are not taken into account in the foreground removal process, future observations would lead to the false detection of the CMB B-mode polarization signals originating from primordial gravitational waves. We also show that the effect of the CMB shadow can be mitigated by our revised Delta-map method to target the CMB B-mode polarization signals at the level of tensor-to-scalar ratio r=0.001.

We present $\texttt{ExTraSS}$ (EXplosive TRAnsient Spectral Simulator), a newly developed code aimed at generating 3D spectra for supernovae in the nebular phase by using modern multi-dimensional explosion models as input. It is well established that supernovae are asymmetric by nature, and that the morphology is encoded in the line profiles during the nebular phase, months after the explosion. In this work, we use $\texttt{ExTraSS}$ to study one such simulation of a $3.3\,M_\odot$ He-core explosion ($M_\text{ejecta}=1.3\,M_\odot$, $E_\text{kin}=1.05\times10^{51}\,$erg) modelled with the $\texttt{Prometheus-HotB}$ code and evolved to the homologous phase. Our code calculates the energy deposition from the radioactive decay of $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe and uses this to determine the Non-Local-Thermodynamic-Equilibrium temperature, excitation and ionization structure across the nebula. From the physical condition solutions we generate the emissivities to construct spectra depending on viewing angles. Our results show large variations in the line profiles with viewing angles, as diagnosed by the first three moments of the line profiles; shifts, widths, and skewness. We compare line profiles from different elements, and study the morphology of line-of-sight slices that determine the flux at each part of a line profile. We find that excitation conditions can sometimes make the momentum vector of the ejecta emitting in the excited states significantly different from that of the bulk of the ejecta of the respective element, thus giving blueshifted lines for bulk receding material, and vice versa. We compare the 3.3 $M_\odot$ He-core model to observations of the Type Ib supernova SN 2007Y.

Cosmology and astrophysics provide various ways to study the properties of dark matter even if they have negligible non-gravitational interactions with the Standard Model particles and remain hidden. We study a type of hidden dark matter model in which the dark matter is completely decoupled from the Standard Model sector except gravitationally, and consists of a single species with a conserved comoving particle number. This category of hidden dark matter includes models that act as warm dark matter but is more general. In particular, in addition to having an independent temperature from the Standard Model sector, it includes cases in which dark matter is in its own thermal equilibrium or is free-streaming, obeys fermionic or bosonic statistics, and processes a chemical potential that controls the particle occupation number. While the usual parameterization using the free-streaming scale or the particle mass no longer applies, we show that all cases can be well approximated by a set of functions parameterized by only one parameter as long as the chemical potential is nonpositive: the characteristic scale factor at the time of the relativistic-to-nonrelativistic transition. We study the constraints from Big Bang Nucleosynthesis, the cosmic microwave background, the Lyman-$\alpha$ forest, and the smallest halo mass. We show that the most significant phenomenological impact is the suppression of the small-scale matter power spectrum -- a typical feature when the dark matter has a velocity dispersion or pressure at early times. So far, small dark matter halos provide the strongest constraint, limiting the transition scale factor to be no larger than $\sim1.4\times10^{-4}$ times the scale factor at matter-radiation equality.

We study the four-dimensional relationships between magnetic activity, rotation, mass and age for solar-type stars in the age range 5-25Myr. This is the late-pre-main sequence (l-PMS) evolutionary phase when rapid changes in star's interior may lead to the changes in magnetic dynamo mechanisms. We carefully derive rotational periods and spot sizes for 471 members of several l-PMS open clusters using photometric light curves from the Zwicky Transient Facility. Magnetic activity was measured in our previous Chandra-based study, and additional rotational data were obtained from other work. Several results emerge. Mass-dependent evolution of rotation through the l-PMS phase agrees with astrophysical models of stellar angular momentum changes, although the data point to a subpopulation of stars with slower initial rotations than commonly assumed. There is a hint of the onset of unsaturated tachoclinal dependency of X-ray activity on rotation, as reported by Argiroffi et al. (2016), but this result is not confidently confirmed. Both X-ray luminosity and star spot area decrease approximately as t^{-1} for solar mass stars suggesting that spot magnetic fields are roughly constant and l-PMS stars follow the universal solar-scaling law between the X-ray luminosity and surface magnetic flux. Assuming convective dynamos are dominant, theoretical magnetic fluxes fail to reveal the universal law for l-PMS stars that enter late Henyey tracks. Altogether we emerge with a few lines of evidence suggesting that the transition from the turbulent to solar-type dynamo occurs at the later stages of l-PMS evolution as stars approach the Zero-Age Main Sequence.

High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high redshift galaxies. Studying these sub-structures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging data from CANDELS and UVCANDELS, we identify star-forming clumps in an HST/F160W$\leq 25$ AB mag sample of 6767 galaxies at $0.5\leq z\leq 3$ in four fields, GOODS-N, GOODS-S, EGS, and COSMOS. We use a low-pass band filter in Fourier space to reconstruct the background image of a galaxy and detect small-scale features (clumps) on the background-subtracted image. Clumpy galaxies are defined as those having at least one off-center clump that contributes a minimum of 10$\%$ of the galaxy's total rest-frame UV flux. We measure the fraction of clumpy galaxies ($\rm f_{clumpy}$) as a function of stellar mass, redshift, and galaxy environment. Our results indicate that $\rm f_{clumpy}$ increases with redshift, reaching $\sim 65\%$ at $z\sim 1.5$. We also find that $\rm f_{clumpy}$ in low-mass galaxies ($\rm 9.5\leq log(M_*/M_\odot)\leq 10$) is 10$\%$ higher compared to that of their high-mass counterparts ($\rm log(M_*/M_\odot)>10.5$). Moreover, we find no evidence of significant environmental dependence of $\rm f_{clumpy}$ for galaxies at the redshift range of this study. Our results suggest that the fragmentation of gas clouds under violent disk instability remains the primary driving mechanism for clump formation, and incidents common in dense environments, such as mergers, are not the dominant processes.

We discuss the prospect of identifying a white dwarf binary merger by monitoring disappearance of its nearly monochromatic gravitational wave. For a ten-year operation of the laser interferometer space antenna (LISA), the chance probability of observing such an event is roughly estimated to be 20%. By simply using short-term coherent signal integrations, we might determine the merger time with an accuracy of $\sim $3-10 days. Also considering its expected sky localizability $\sim0.1$-$ 0.01 {\rm deg^2}$, LISA might make an interesting contribution to the multi-messenger study on a merger event.

The primary purpose of this work is the provision of accurate, analytic, evolutionary templates for cosmological parameters and fundamental constants in a dynamical cosmology. A flat quintessence cosmology with a dark energy potential that has the mathematical form of the Higgs potential is the specific cosmology and potential addressed in this work. These templates, based on the physics of the cosmology and potential are intended to replace the parameterizations currently used to determine the likelihoods of dynamical cosmologies. Acknowledging that, unlike {\Lambda}CDM, the evolutions are dependent on both the specific cosmology and the dark energy potential the templates are referred to as Specific Cosmology and Potential, SCP, templates. The requirements set for the SCP templates are that they must be accurate, analytic functions of an observable such as the scale factor or redshift. This is achieved through the utilization of a modified beta function formalism that is based on a physically motivated dark energy potential to calculate the beta function. The methodology developed here is designed to be adaptable to other cosmologies and dark energy potentials. The SCP templates are essential tools in determining the relative likelihoods of a range of dynamical cosmologies and potentials. An ultimate purpose is the determination whether dark energy is dynamical or static in a quantitative manner. It is suggested that the SCP templates calculated in this work can serve as fiducial dynamical templates in the same manner as {\Lambda}CDM serves for static dark energy.

Stochastic temporal wandering of the spin frequency $\nu$ of a rotation-powered pulsar (i.e.~the achromatic component of timing noise unrelated to interstellar propagation) affects the accuracy with which the secular braking torque can be measured. Observational studies confirm that pulsars with anomalous braking indices $\vert n \vert = \vert \nu \ddot{\nu} / \dot{\nu}^2 \vert \gg 1$ exhibit elevated levels of timing noise, where an overdot symbolizes a derivative with respect to time. Here it is shown, through analytic calculations and Monte Carlo simulations involving synthetic data and modern Bayesian timing techniques, that the variance $\langle n^2 \rangle$ of the measured $n$ scales with the square of the timing noise amplitude $\sigma_{\ddot{\nu}}$. The anomalous regime $\langle n^2 \rangle \gg 1$ corresponds to $ \sigma_{\ddot{\nu}}^2 \gg 10^{-60} (\gamma_{\ddot{\nu}}/10^{-6} \, {\rm s^{-1}})^2 (\dot{\nu} / 10^{-14} \, {\rm Hz \, s^{-1}})^4 (\nu / 1 \, {\rm Hz})^{-2} (T_{\rm obs} / 10^8 \, {\rm s}) \, {\rm Hz}^2{\rm s}^{-5 }$, where $\gamma_{\ddot{\nu}}$ is a stellar damping time-scale, and $T_{\rm obs}$ is the total observing time. When the inequality in the above condition is reversed, $n$ is dominated by the secular braking torque, and timing measurements return $n\sim 3$, if the secular braking torque is electromagnetic. The variance $\langle n^2 \rangle$ is greater, when the stochastic process driving spin fluctuations differs from the red noise model (e.g. power-law spectral density) assumed in the timing solution.

A planetary embryo embedded in a gaseous disc can grow by pebble accretion while subjected to a gravitational force from the disc that changes its orbital elements. Usually, that force is considered to arise from the Lindblad and corotation resonances with the embryo. However, more important contributions exist for low-mass planets. Radiative thermal diffusion in the vicinity of embryos yields an additional contribution to the disc's force that damps the eccentricity and inclination much more vigorously than the resonant interaction with the disc, and that in general induces fast inward migration. In addition, the irradiation of the disc by a hot embryo gives rise to an additional contribution that excites eccentricity and inclination, and induces outward migration. Which of the two contributions dominates depends on the embryo's luminosity. We assess the importance of these contributions (termed thermal forces) on the dynamics and growth of a set of pebble-accreting embryos initially of Martian mass, by means of N-body simulations that include analytic expressions for the disc's force. We find very different outcomes for the embryos subjected to thermal forces and those subjected only to resonant forces. Importantly, we find that the median final mass of the embryos subjected to thermal forces is nearly independent of the metallicity, whereas this mass roughly scales with the metallicity when they are subjected only to resonant forces. These results can be explained by the strong damping of eccentricity and inclination at low metallicity, which enhances the embryos' accretion efficiency.

In radio astronomy, signals from radio telescopes are transformed into images of observed celestial objects, or sources. However, these images, called dirty images, contain real sources as well as artifacts due to signal sparsity and other factors. Therefore, radio interferometric image reconstruction is performed on dirty images, aiming to produce clean images in which artifacts are reduced and real sources are recovered. So far, existing methods have limited success on recovering faint sources, preserving detailed structures, and eliminating artifacts. In this paper, we present VIC-DDPM, a Visibility and Image Conditioned Denoising Diffusion Probabilistic Model. Our main idea is to use both the original visibility data in the spectral domain and dirty images in the spatial domain to guide the image generation process with DDPM. This way, we can leverage DDPM to generate fine details and eliminate noise, while utilizing visibility data to separate signals from noise and retaining spatial information in dirty images. We have conducted experiments in comparison with both traditional methods and recent deep learning based approaches. Our results show that our method significantly improves the resulting images by reducing artifacts, preserving fine details, and recovering dim sources. This advancement further facilitates radio astronomical data analysis tasks on celestial phenomena.

Main Sequence Turn-off stars (MSTO) and subgiant stars are good tracers of galactic populations. We present a study of 41,034 MSTO and subgiant stars from the GALAH survey. Using a grid of stellar models that accounts for the variation of O abundances, we determine their ages with a median age uncertainty of $\sim$9.4 per cent. Our analysis reveals that the ages of high-O stars based on O-enhanced models (OEM models) are smaller than those determined with $\alpha$-enhanced models, resulting in a mean fractional age difference of -5.3 per cent at [O/$\alpha$] = 0.2 and -11.0 per cent at [O/$\alpha$] = 0.4. This age difference significantly impacts the age distribution of thick disc and halo stars, leading to a steeper downward trend in the [Fe/H]-age plane from 8 Gyr to 14 Gyr, indicating a shorter formation time-scale and a faster chemical-enhanced history for these populations. We confirm the V-shape of the normalized age-metallicity distribution $p$($\tau$$\mid$[Fe/H]) of thin disc stars, which is presumably a consequence of the second gas infall. Additionally, we find that the halo stars in our sample can be divided into two sequences, a metal-rich sequence (Splash stars) and a metal-poor sequence (accreted stars), with the Splash stars predominantly older than 9 Gyr and the accreted halo stars older than 10 Gyr. Finally, we observe two distinct sequences in the relations between various chemical abundances and age for disc stars, namely a young sequence with ages $<$ $\sim$8 Gyr and an old sequence with ages $>$ $\sim$8 Gyr.

Hot Jupiters have been perceived as loners devoid of planetary companions in close orbital proximity. However, recent discoveries based on space-borne precise photometry have revealed that at least some fraction of giant planets coexists with low-mass planets in compact orbital architectures. We report detecting a 1.446-day transit-like signal in the photometric time series acquired with the Transiting Exoplanet Survey Satellite (TESS) for the WASP-84 system, which is known to contain a hot Jupiter on a circular 8.5-day orbit. The planet was validated based on TESS photometry, and its signal was distilled in radial velocity measurements. The joint analysis of photometric and Doppler data resulted in a multi-planetary model of the system. With a mass of $15\, M_{\oplus}$, radius of $2\, R_{\oplus}$, and orbital distance of 0.024 au, the new planet WASP-84 c was classified as a hot super-Earth with the equilibrium temperature of 1300 K. A growing number of companions to hot Jupiters indicates that a non-negligible part of them must have formed under a quiescent scenario such as disc migration or in-situ formation.

Numerical simulations predict that the spin-down rate of a single rotation-powered neutron star depends on the angle $\alpha$ between its spin and magnetic axes as $P\dot P \propto \mu^2 (k_0 + k_1\sin^2\alpha)$, where $P$ is the star spin period, $\mu$ is its magnetic moment, while $k_0 \sim k_1 \sim 1$. Here we describe a simple observational test for this prediction based on the comparison of spin-down rates of 50 nearly orthogonal (with $\alpha$ close to 90 deg) and 27 nearly aligned (with $\alpha$ close to 0 deg) pulsars. We found, that the apparent pulsar spin-down is consistent with the theory if assumed, that magnetic moments of orthogonal rotators are systematically larger than those of aligned ones for $\sim 0.15..0.2$ dex. Also, as a by-product of the analysis, we provide yet another constraint on the average braking index of radio pulsars as $1 \le n \le 4$ with formal significance not worse than 99\%.

This article reports a sample of 1830 strong [O III] {\lambda}5007 emission-line compact galaxies discovered with the LAMOST spectroscopic survey and the photometric catalog of SDSS. We newly identify 402 spectra of 346 strong [O III]{\lambda}5007 emission-line compact galaxies by finding compact isolated point sources. Combined with the samples in our previous work (Liu et al. 2022), this returns a sample of 1830 unique strong [O III]{\lambda}5007 emission-line compact galaxies with 2033 spectra of z <= 0.53. For the sources with 2{\sigma} [OIII]{\lambda}4363 detections, we calculate the gas-phase metallicity with the direct-Te method, and verify that the strong-line metallicity diagnostics calibrated with the direct-Te method also applies to this sample. The strong [O III]{\lambda}5007 emission-line compact galaxies fall below several Te-calibrated mass-metallicity relations. The N/O measurements of the strong [O iii]{\lambda}5007 emission-line compact galaxies mainly locate at a plateau at low metallicity, indicating the product of primary nucleosynthesis. The Ne3O2 and O32 relation follows a tight linear relation with no redshift evolution. The Ne3O2 anti-correlates with the stellar mass, and at fixed stellar mass the Ne3O2 increase with the redshift. Eight sources with asymmetric [O III]{\lambda}5007 emission-line profiles have been identified, however with no [O III]{\lambda}4363 detection, which proves the rich metal content and complex ionized gas kinematics within the galaxies. Higher-resolution spectroscopy will be necessary to identify the ionized gas components in detail.

We use the latest parallaxes measurements from Gaia DR3 to obtain a geometric calibration of the tip of the red giant branch (TRGB) in Cousins $I$ magnitudes as a standard candle for cosmology. We utilise the following surveys: SkyMapper DR3, APASS DR9, ATLAS Refcat2, and Gaia DR3 synthetic photometry to obtain multiple zero-point calibrations of the TRGB magnitude, $M_{I}^{TRGB}$. Our sample contains Milky Way halo stars at high galactic latitudes ($|b| > 36$) where the impact of metallicity, dust, and crowding are minimised. The magnitude of the TRGB is identified using Sobel edge detection, but this approach introduced a systematic offset. To address this issue, we utilised simulations with PARSEC isochrones and showed how to calibrate and remove this bias. Applying our method within the colour range where the slope of the TRGB is relatively flat for metal-poor halo stars (1.55 $<$ $(BP-RP)$ $<$ 2.25), we find a weighted average $M_{I}^{TRGB} = -4.042 \pm 0.041$ (stat) $\pm0.031$ (sys) mag. A geometric calibration of the Milky Way TRGB has the benefit of being independent of other distance indicators and will help probe systematics in the local distance ladder, leading to improved measurements of the Hubble constant.

Classical T Tauri stars are low-mass pre-main sequence stars with an active circumstellar environment. In this work we present the identification and study of 260 Classical T Tauri stars using LAMOST Data Release 6, among which 104 stars are newly identified. We distinguish Classical T Tauri stars from Giants and main-sequence dwarfs based on the log g values and the presence of H (alpha) emission line and infrared excess that arises from the circumstellar accretion disk. We estimated the mass and age of 210 stars using the Gaia color-magnitude diagram. The age is from 0.1 to 20 Myr, where 90% of the stars have age below 10 Myr and the mass ranges between 0.11 to 1.9 M(solar). From the measured H(alpha) equivalent widths, we homogeneously estimated the mass accretion rates for 172 stars, with most values ranging from 10^-7 to 10^-10 M(solar) yr^-1. The mass accretion rates are found to follow a power law distribution with the mass of the star, having a relation of the form Macc proportional to M(star)^1.43 +/- 0.26, in agreement with previous studies.

The European Southern Observatory (ESO) implemented a new paradigm called Distributed Peer Review (DPR) as part of its proposal evaluation process in Period 110. Under DPR, Principal Investigators who submit proposals agree to review a certain number of proposals submitted by their peers and accept that their own proposal(s) are reviewed by their peers who have also submitted proposals in the same cycle. This article presents a brief overview of the DPR process at ESO, and its outcomes based on data from periods 110 and 111.

In this paper, we found 2939 S-type stars from LAMOST Data Release 10 using two machine-learning methods, and 2306 of them were reported for the first time. The main purpose of this work is to study how to divide S-type stars into intrinsic and extrinsic stars with photometric data and LAMOST spectra. Using infrared photometric data, we adopted two methods to distinguish S-type stars, i.e., XGBoost algorithm and color-color diagrams. We trained XGBoost model with 15 input features consisting of colors and absolute magnitudes of Two Micron All Sky Survey (2MASS), AllWISE, AKARI, and IRAS, and found that the model trained by input features with 2MASS, AKARI, and IRAS data has the highest accuracy of 95.52%. Furthermore, using this XGBoost model, we found four color-color diagrams with six infrared color criteria to divide S-type stars, which has an accuracy of about 90%. Applying the two methods to the 2939 S-type stars, 381 (XGBoost)/336 (color-color diagrams) intrinsic and 495 (XGBoost)/82 (color-color diagrams) extrinsic stars were classified, respectively. Using these photometrically classified intrinsic and extrinsic stars, we retrained XGBoost model with their blue and red medium-resolution spectra, and the 2939 stars were divided into 855 intrinsic and 2056 extrinsic stars from spectra with an accuracy of 94.82%. In addition, we also found four spectral regions of Zr I (6451.6A), Ne II (6539.6A), H{\alpha} (6564.5A), and Fe I (6609.1A) and C I (6611.4A) are the most important features, which can reach an accuracy of 92.1% when using them to classify S-type stars.

Solar flares are energetic events in the solar atmosphere that are often linked with solar radio bursts (SRBs). SRBs are observed at metric to decametric wavelengths and are classified into five spectral classes (Type I--V) based on their signature in dynamic spectra. The automatic detection and classification of SRBs is a challenge due to their heterogeneous form. Near-realtime detection and classification of SRBs has become a necessity in recent years due to large data rates generated by advanced radio telescopes such as the LOw Frequency ARray (LOFAR). In this study, we implement congruent deep learning models to automatically detect and classify Type III SRBs. We generated simulated Type III SRBs, which were comparable to Type IIIs seen in real observations, using a deep learning method known as Generative Adversarial Network (GAN). This simulated data was combined with observations from LOFAR to produce a training set that was used to train an object detection model known as YOLOv2 (You Only Look Once). Using this congruent deep learning model system, we can accurately detect Type III SRBs at a mean Average Precision (mAP) value of 77.71%.

Perturbation theory is a powerful tool for studying large-scale structure formation in the universe and calculating observables such as the power spectrum or bispectrum. However, beyond linear order, typically this is done by assuming a simplification in the time-dependence of gravitational-coupling kernels between the matter and velocity fluctuations. Though the true dependencies are known for Lambda cold dark matter cosmologies, they are ignored due to the computational costs associated with considering them in full and, instead, are replaced by simpler dependencies valid for an Einstein--de-Sitter cosmology. Here we develop, implement and demonstrate the effectiveness of a new numerical method for finding the full dynamical evolution of these kernels to all perturbative orders based upon spectral methods using Chebyshev polynomials. This method is found to be orders of magnitude more efficient than direct numerical solvers while still producing highly accurate and reliable results. A code implementation of the Chebyshev spectral method is then presented and characterised. The code has been made publicly available alongside this paper. We expect our method to be of use for interpretation of upcoming galaxy clustering measurements.

We examine the Galactic warp in a sample of all classical Cepheids with Gaia DR3 radial velocity. In each radial bin, we determine (1) the inclined plane normal to the mean orbital angular momentum of the stars and (2) that best fitting their positions. We find no warping inside $R\approx 11$ kpc; for larger $R$ the disc is increasingly inclined, reaching $i\sim 3^{\circ}$ at $R \ge 14$ kpc. With larger $R$ the azimuth of the warp's ascending node shifts from $\varphi_{\mathrm{lon}}\approx-15^\circ$ at 11 kpc by about $14^{\circ}$/kpc in the direction of Galactic rotation, implying a leading spiral of nodes, the general behaviour of warped galaxies. From the method of fitting planes to the positions we also obtain $\dot{\varphi}_{\mathrm{lon}}$ and find prograde precession of $\dot{\varphi}_{\mathrm{lon}} \sim 12$ km/s/kpc at 12 kpc decreasing to $\sim 6$ km/s/kpc at 14 kpc and beyond. This would unwind the leading spiral of nodes in $\sim 100$ Myr, suggesting that our instantaneous measurements of $\dot{\varphi}_{\mathrm{lon}}$ reflect transient behaviour. This is consistent with existing simulations, which show oscillations in $\dot{\varphi}_{\mathrm{lon}}$ overlaying a long-term retrograde differential precession which generates the leading spiral of nodes.

Triaxial dynamical models of massive galaxies observed in the ATLAS3D project can provide new insights into the complex evolutionary processes that shape galaxies. The ATLAS3D survey is ideal as the sample comprises a good mix of fast and slow rotators with vastly different mass assembly histories. We present a detailed dynamical study with our triaxial modelling code DYNAMITE, which models galaxies as a superposition of their stellar orbits. The models allow us to constrain the intrinsic shape of the stellar component, the distributions of the visible and invisible matter and the orbit distribution in these nearby early-type galaxies and to relate it with different evolutionary scenarios. Triaxial modelling is essential for these galaxies to understand their complex kinematical features.

We investigate the coalescence of massive black hole ($M_{\rm BH}\gtrsim 10^{6}~\rm M_{\odot}$) binaries (MBHBs) at $6<z<10$ by adopting a suite of cosmological hydrodynamical simulations of galaxy formation, zoomed-in on biased ($ >3 \sigma$) overdense regions ($M_h\sim 10^{12}~\rm M_{\odot}$ dark matter halos at $z = 6$) of the Universe. We first analyse the impact of different resolutions and AGN feedback prescriptions on the merger rate, assuming instantaneous mergers. Then, we compute the halo bias correction factor due to the overdense simulated region. Our simulations predict merger rates that range between 3 - 15 $\rm yr^{-1}$ at $z\sim 6$, depending on the run considered, and after correcting for a bias factor of $\sim 20-30$. For our fiducial model, we further consider the effect of delay in the MBHB coalescence due to dynamical friction. We find that 83 per cent of MBHBs will merge within the Hubble time, and 21 per cent within 1 Gyr, namely the age of the Universe at $z > 6$. We finally compute the expected properties of the gravitational wave (GW) signals and find the fraction of LISA detectable events with high signal-to-noise ratio (SNR $>$ 5) to range between 66-69 per cent. However, identifying the electro-magnetic counterpart of these events remains challenging due to the poor LISA sky localization that, for the loudest signals ($\mathcal M_c\sim 10^6~\rm M_{\odot}$ at $z=6$), is around 10 $\rm deg^2$.

X-ray polarimetry is a powerful tool to probe the geometry of the hot X-ray corona in active galactic nuclei (AGN). Here, we present our results on the characterisation of the X-ray polarization of the radio-quiet Seyfert-type AGN IC 4329A at a redshift of $z$ = 0.016. This is based on observations carried out by the {\it Imaging X-ray Polarimeter (IXPE)}. {\it IXPE} observed IC 4329A on January 5, 2023, for a total observing time of 458 ks. From the model-independent analysis, we found a polarization degree ($\Pi_{X}$) of 3.7$\pm$1.5$\%$ and a polarization position angle ($\Psi_{X}$) of 61$^{\circ}$$\pm$12$^{\circ}$ in the 2$-$8 keV energy range (at 1$\sigma$ confidence level). This is also in agreement with the values of $\Pi_{X}$ and $\Psi_{X}$ of 4.7$\pm$2.2$\%$ and 71$^{\circ}$ $\pm$14$^{\circ}$ respectively obtained from spectro-polarimetric analysis of the I, Q and U Stokes spectra in the 2$-$8 keV energy band (at the 90$\%$ confidence). The value of $\Pi_X$ in the 2-8 keV band obtained from the model-independent analysis is lower than the minimum detectable polarization (MDP) value of 4.5$\%$. However, $\Pi_X$ obtained from spectro-polarimetric analysis in the 2-8 keV band is larger than the MDP value. In the 3-5 keV band, we found $\Pi_X$ of 6.5 $\pm$ 1.8, which is larger than the MDP value of 5.5$\%$. The observed moderate value of $\Pi_{X}$ obtained from the analysis of the {\it IXPE} data in the 3$-$5 keV band argues against a spherical lamp$-$post geometry for the X-ray corona in IC 4329A; however, considering simulations, the observed polarization measurements tend to favour a conical shape geometry for the corona. This is the first time measurement of X-ray polarization in IC 4329A. Measurements of the X-ray polarization in many such radio-quiet AGN will help in constraining the geometry of the X-ray corona in AGN.

In the expanding Universe, the dimensional quantities like the wavelength and the temperature of photons are cosmologically redshifted by the relative difference between the observed and emitted ones. Therefore, it can be physically meaningful to talk about the constancy or variability of any dimensional constant (not only of dimensionless one) when the Universe is expanding. It has been known that one can measure the temporal variation of the fine structure constant $\alpha$ in the emission and absorption lines of quasar spectra when the speed of light varies for cosmic time, even though this statement is model dependent. Current observations based on the alkali doublet method and on the many-multiplet one show superficially contradictory results. The former finds no statistically significant evidence for a time dependence of $\alpha$, while the latter does. The so-called meVSL model can reconcile these results naturally without any contradiction.

The superposition of many astrophysical gravitational waves (GW) signals below typical detection thresholds baths detectors in a stochastic gravitational wave background (SGWB). In this work we present a Fourier space approach to compute the frequency-domain distribution of stochastic gravitational wave backgrounds produced by discrete sources. The expressions for the moment generating function and the distribution of observed (discrete) Fourier modes are provided. The results are then applied to the SGWB originating from the mergers of compact stellar remnants (black holes and neutron stars) in the Universe, which are found to exhibit a $-4$ power-law tail. This tail is verified in the signal-to-noise ratio distribution of GWTC events. Furthermore, the extent to which the subtraction of bright (loud) mergers gaussianizes the resulting confusion noise of unresolved sources is illustrated. The power-law asymptotic tail for the SGWB, and an exponentially decaying tail in the case of the confusion background, are also derived analytically. Our results generalize to any background of gravitational waves emanating from discrete sources.

The majority of all stars are members of a binary system. The evolution of such binary stars and their subsequent production of pairs of compact objects in tight orbits, such as double neutron stars and double black holes, play a central role in modern astrophysics, Binary evolution leads to the formation of different types of violent cosmic events such as novae, supernova explosions, gamma-ray bursts, mass transfer and accretion processes in X-ray binaries, and the formation of exotic radio millisecond pulsars. In some cases, the binary systems terminate as spectacular collisions between neutron stars and/or black holes. These collisions lead to powerful emission of gravitational waves, as detected by LIGO since 2015. The coming decade is expected to reveal a large number of discoveries of binary compact systems, as well as their progenitors and merger remnants, from major instruments such as the radio Square-Kilometre Array; the gravitational wave observatories LIGO-Virgo-KAGRA-IndIGO and LISA; the astrometric space observatory Gaia; the James Webb Space Telescope; and the X-ray space observatories eXTP, STROBE-X, and Athena. In this light, it is important to have a modern textbook on the physics of binary stars evolution, from ordinary stars to X-ray binaries and gravitational wave sources. The scope of this book is that the reader (student or educated expert) will learn the physics of binary interactions, from stellar birth to compact objects, and relate this knowledge to the latest observations. The reader will learn about stellar structure and evolution, and detailed binary interactions covering a broad range of phenomena, including mass transfer and orbital evolution, formation and accretion onto compact objects (white dwarfs, neutron stars and black holes), and their observational properties. Exercises are provided throughout the book.

Observational studies have identified several sub-structures in different regions of the Magellanic Clouds, the nearest pair of interacting dwarf satellites of the Milky Way. By studying the metallicity of the sources in these sub-structures, we aim to shed light on the possible origin of these sub-structures. Spectroscopic metallicities exist only for a few thousand sources, mostly giant stars located in specific regions of the galaxies. These metallicities come from different instruments at various spectral resolutions, and systematic uncertainties hamper comparisons and draw firm conclusions about their origin. The third data release of \textit{Gaia} has provided us with $\sim$ 0.17 million XP spectra of the different stellar populations in the SMC alone as faint as $\sim$ 18 mags in the G band, which are spread across $\sim$ 10$^\circ$ from the SMC centre. We aim to determine the metallicities of these sources based on synthetic Str\"{o}mgren photometry derived from XP spectra and produce a high-resolution metallicity map of the SMC. Our metallicity gradient estimate of the SMC turns out to be --0.062 $\pm$ 0.009 dex/deg. This is comparable with the previous estimates, which also validate our method of metallicity estimation. We aim to apply this method to other stellar populations and to the LMC to create a high-resolution metallicity map of the Magellanic Clouds.

We present new period-luminosity and period-luminosity-metallicity relations for Galactic RR Lyrae stars based on a sample of 28 pulsators located at distances up to $1.5$ kpc from the Sun. Near-infrared photometry was obtained at the Cerro Armazones Observatory and parallaxes were taken from the Gaia Early Data Release 3. Relations were determined for the 2MASS $JHK_s$ bands and the $W_{JK}$ Wesenheit index. We compare our results with other calibrations available in the literature and obtain very good agreement with the photometry of RR Lyraes from the Large Magellanic Cloud anchored using the distance to the Cloud, which based on detached eclipsing binaries. We find that the dependence of absolute magnitudes on metallicity of $0.070\pm 0.042$ mag/dex ($J-$ band) to $0.087 \pm 0.031$ mag/dex ($W_{JK}$ index) for the population of fundamental pulsators (RRab) that is in agreement with previously published phenomenological works. We perform a refined determination of distance to the LMC based on our new calibration and photometry from Szewczyk et al. (2008). We study the dependence of the fitted parameters of fiducial relations and the LMC distance on the systematic parallax offset.

Fast radio bursts (FRBs) are a recently discovered class of GHz-band, ms-duration, Jy-level-flux astrophysical transients, which origin is still a mystery. Exploring their gamma-ray counterpart is crucial for constraining their origin and emission mechanism. Thanks to more than 13 years of gamma-ray data collected by the Fermi-Large Area Telescope, and to more than 1000 FRB events, one of the largest sample created as of today, we perform the largest and deepest search for gamma-ray emission from FRB sources to date. In addition to the study of individual FRB events on different time-scales (from few seconds up to several years), we performed, for the first time, a stacking analysis on the full sample of FRB events as well as a search for triplet photons in coincidence with the radio event. We do not detect significant emission, reporting the most stringent constraints, on short time scales, for the FRB-like emission from SGR 1935+2154 with $E<10^{41}$ erg, corresponding to a factor $<10^7$ with respect to the emitted radio energy. For the stacked signal of steady emission from all repeaters, the obtained upper limit (UL) on the FRBs luminosity ($L<1.6\times10^{43}$ erg s$^{-1}$) is more than two orders of magnitudes lower than those derived from the individual sources. Finally, no individual or triplet photons have been significantly associated with FRB events. We derived the LAT ms energy sensitivity to be $E<10^{47}$ (D$_L$/150 Mpc)$^2$ erg, ruling out a gamma-ray-to-radio energy ratio greater than $10^9$ on ms timescales. The results reported here represent the most stringent UL reported so far on the high-energy emission from FRBs on short and long time scales, as well as on cumulative emission and individual photon searches. While the origin of FRBs is still unclear, our work provides important constraints for FRB modeling, which might shed light on their emission mechanism.

Emission lines indicative of active accretion have been seen on a handful of low-mass companions (M < 30 MJup) to stars. Line variability is ubiquitous on stellar accretors but has never been characterized in detail on low-mass companions and can give insights on the accretion mechanism at play. We investigate the emission line variability of two low-mass companions (M<30 MJup) to stars to understand their accretion mechanisms. Using J-band observations, we analyze the short to long-term variability of the HI Paschen {\beta} emission line (1.282 {\mu}m) for GQ Lup b and GSC 06214-00210 b. Archival spectroscopic observations are also examined to extend the time span. We compare their line profiles and intensities to more massive accretors and magnetospheric accretion and shock models. Both objects have HI Paschen {\beta} flux variability that is moderate at short timescales (< 50 %) and increases at longer timescales (~1000 % on decade timescales), resembling classical T Tauri stars. GQ Lup b's line profiles are compatible with magnetospheric accretion. GSC 06214-00210 b's profiles are reproduced by both magnetospheric accretion and shock models, except for the brightest epoch for which the shock model is highly favored. Both companions have C/O values broadly consistent with solar values. While magnetospheric accretion is favored for GQ Lup b, higher resolution (R > 10000) observations are required to disentangle the two (non-exclusive) line formation mechanisms. The similarity in variability behavior may support similar accretion mechanisms between these low-mass companions and classical T Tauri stars. The significant variability observed at months and longer timescales could explain the low yield of H{\alpha} imaging campaigns.

This study presents the largest mid-infrared spectral catalog of Herbig Ae/Be stars to date, containing the Spitzer Infrared Spectrograph spectra of 126 stars. Based on the catalog analysis, two prominent infrared vibrational modes of C\textsubscript{60} bands at 17.4 $\mu m$ and 18.9 $\mu m$ are detected in the spectra of nine sources, while 7.0 $\mu m$ feature is identified in the spectra of HD 319896. The spectral index analysis and the comparison of the known sources with C\textsubscript{60} features indicated that there exist two different types of emission classes among the sample of stars. The infrared spectra of six Herbig Ae/Be stars in this study resemble that of reflection nebulae, and their association with previously known reflection nebulae is confirmed. In the case of three Herbig Ae/Be stars we report the tentative evidence of C\textsubscript{60} emission features originating from the circumstellar disk or nearby diffused emission region. The detection fraction of C\textsubscript{60} in the total HAeBe star sample is $\sim$ 7\%, whereas the detection fraction is 30\% for HAeBe stars associated with nebulosity. In the catalog, C\textsubscript{60} is exclusively present in the circumstellar regions of B type Herbig Ae/Be stars, with no evidence of its presence detected in stars with later spectral types. The present study has increased the number of young stellar objects and reflection nebulae detected with C\textsubscript{60} multifold, which can help in understanding the excitation and formation pathway of the species.

In this dataset we provide a comprehensive collection of magnetograms (images quantifying the strength of the magnetic field) from the National Aeronautics and Space Administration's (NASA's) Solar Dynamics Observatory (SDO). The dataset incorporates data from three sources and provides SDO Helioseismic and Magnetic Imager (HMI) magnetograms of solar active regions (regions of large magnetic flux, generally the source of eruptive events) as well as labels of corresponding flaring activity. This dataset will be useful for image analysis or solar physics research related to magnetic structure, its evolution over time, and its relation to solar flares. The dataset will be of interest to those researchers investigating automated solar flare prediction methods, including supervised and unsupervised machine learning (classical and deep), binary and multi-class classification, and regression. This dataset is a minimally processed, user configurable dataset of consistently sized images of solar active regions that can serve as a benchmark dataset for solar flare prediction research.

NASA's Kepler Space Telescope has been instrumental in the task of finding the presence of exoplanets in our galaxy. This search has been supported by computational data analysis to identify exoplanets from the signals received by the Kepler telescope. In this paper, we consider building upon some existing work on exoplanet identification using residual networks for the data of the Kepler space telescope and its extended mission K2. This paper aims to explore how deep learning algorithms can help in classifying the presence of exoplanets with less amount of data in one case and a more extensive variety of data in another. In addition to the standard CNN-based method, we propose a Siamese architecture that is particularly useful in addressing classification in a low-data scenario. The CNN and ResNet algorithms achieved an average accuracy of 68% for three classes and 86% for two-class classification. However, for both the three and two classes, the Siamese algorithm achieved 99% accuracy.

The origin of the gamma-ray emission from M87 is currently a matter of debate. This work aims to localize the VHE (100 GeV-100 TeV) gamma-ray emission from M87 and probe a potential extended hadronic emission component in the inner Virgo Cluster. The search for a steady and extended gamma-ray signal around M87 can constrain the cosmic-ray energy density and the pressure exerted by the cosmic rays onto the intra-cluster medium, and allow us to investigate the role of the cosmic rays in the active galactic nucleus feedback as a heating mechanism in the Virgo Cluster. H.E.S.S. telescopes are sensitive to VHE gamma rays and have been utilized to observe M87 since 2004. We utilized a Bayesian block analysis to identify M87 emission states with H.E.S.S. observations from 2004 until 2021, dividing them into low, intermediate, and high states. Because of the causality argument, an extended ($\gtrsim$kpc) signal is allowed only in steady emission states. Hence, we fitted the morphology of the 120h low state data and found no significant gamma-ray extension. Therefore, we derived for the low state an upper limit of 58"(corresponding to $\approx$4.6kpc) in the extension of a single-component morphological model described by a rotationally symmetric 2D Gaussian model at 99.7% confidence level. Our results exclude the radio lobes ($\approx$30 kpc) as the principal component of the VHE gamma-ray emission from the low state of M87. The gamma-ray emission is compatible with a single emission region at the radio core of M87. These results, with the help of two multiple-component models, constrain the maximum cosmic-ray to thermal pressure ratio $X_{{CR,max.}}$$\lesssim$$0.32$ and the total energy in cosmic-ray protons (CRp) to $U_{CR}$$\lesssim$5$\times10^{58}$ erg in the inner 20kpc of the Virgo Cluster for an assumed CRp power-law distribution in momentum with spectral index $\alpha_{p}$=2.1.

Gravity-mode asteroseismology has significantly improved our understanding of mixing in intermediate mass stars. However, theoretical pulsation periods of stellar models remain in tension with observations, and it is often unclear how the models of these stars should be further improved. Inversions provide a path forward by directly probing the internal structure of these stars from their pulsation periods, quantifying which parts of the model are in need of improvement. This method has been used for solar-like pulsators, but has not yet been applied to main-sequence gravity-mode pulsators. Our aim is to determine whether structure inversions for gravity-mode pulsators are feasible. We focus on the case of slowly rotating SPB stars. We computed and analyzed dipole mode kernels for three variables pairs: $(\rho,c), (N^2,c)$, and $(N^2,\rho)$. We assessed the potential of these kernels by predicting the oscillation frequencies of a model after perturbing its structure. We then tested two inversion methods, RLS and SOLA, using a model grid computed with MESA and GYRE. We find that changing the stellar structure affects the oscillation frequencies in a nonlinear way. The oscillation modes for which this nonlinear dependency is the strongest are in resonance with the near-core peak in the buoyancy frequency. The near core region of the star can be probed with SOLA, while RLS requires fine tuning to obtain accurate results. Both RLS and SOLA are strongly affected by the nonlinear dependencies on the structure differences, as these methods are based on a first-order approximation. These inversion methods need to be modified for meaningful applications of inversions to SPB stars. Our results show that inversions of gravity-mode pulsators are possible in principle, but that the typical inversion methods developed for solar-like oscillators are not applicable. [abridged]

The distribution of baryons provides a significant way to understand the formation of galaxy clusters by revealing the details of its internal structure and changes over time. In this paper, we present theoretical studies on the scaled profiles of physical properties associated with the baryonic components, including gas density, temperature, metallicity, pressure and entropy as well as stellar mass, metallicity and satellite galaxy number density in galaxy clusters from $z=4$ to $z=0$ by tracking their progenitors. These mass-complete simulated galaxy clusters are coming from THE THREE HUNDRED with two runs: GIZMO-SIMBA and Gadget-X. Through comparisons between the two simulations, and with observed profiles which are generally available at low redshift, we find that (1) the agreements between the two runs and observations are mostly at outer radii $r \gtrsim 0.3r_{500}$, in line with the self-similarity assumption. While Gadget-X shows better agreements with the observed gas profiles in the central regions compared to GIZMO-SIMBA; (2) the evolution trends are generally consistent between the two simulations with slightly better consistency at outer radii. In detail, the gas density profile shows less discrepancy than the temperature and entropy profiles at high redshift. The differences in the cluster centre and gas properties imply different behaviours of the AGN models between Gadget-X and GIZMO-SIMBA, with the latter, maybe too strong for this cluster simulation. The high-redshift difference may be caused by the star formation and feedback models or hydrodynamics treatment, which requires observation constraints and understanding.

We investigate the formation of large peaks in the inflationary curvature power spectrum from double-well potentials. In such scenarios, the initial CMB spectrum is created at large field values. Subsequently, the inflaton will cross one of the minima and will decelerate rapidly as it reaches the local maximum at the origin, either falling back or crossing it. During this final phase, a significant peak in the curvature power spectrum can be generated. Our analysis reveals that this class of models produces more pronounced peaks than the more commonly studied quasi-inflection point scenarios with less tuning for the model parameters. Finally, we construct an explicit theoretically motivated inflationary scenario that is consistent with the latest CMB observations and capable of generating sufficiently large curvature perturbations for primordial black holes.

Several analyses have shown that LCDM-based models cannot jointly describe the clustering (GC) and galaxy-galaxy lensing (GGL) of galaxies in the SDSS-III BOSS survey, which is commonly known as the 'lensing-is-low problem'. In this work, we show that an extension of Subhalo Abundance Matching, dubbed SHAMe, successfully solves this problem. First, we show that this model accurately reproduces the GC and GGL of a mock galaxy sample in the TNG300 hydrodynamic simulation with analogous properties to BOSS galaxies. Then, we switch our attention to observed BOSS galaxies at z=0.31-0.43, and we attempt to reproduce their GC and GGL by evaluating SHAMe on two different simulations: one adopting best-fitting cosmological parameters from Planck and the other from weak gravitational lensing surveys (Low S8), where the amplitude of matter fluctuations is lower for the latter. We find excellent agreement between SHAMe predictions and observations for both cosmologies, indicating that the lensing-is-low problem originates from approximations in previous theoretical descriptions of the data. The main difference between SHAMe results in these cosmologies is the level of galaxy assembly bias, which is approximately 20 and 10% for Planck and Low S8, respectively. These results highlight the dangers of employing oversimplified models to analyse current large-scale structure datasets, and the need for realistic yet flexible descriptions of the galaxy-halo connection.

Photochemical hazes are expected to form in hot Jupiter atmospheres and may explain the strong scattering slopes and muted spectral features observed in the transmission spectra of many hot Jupiters. Absorption and scattering by photochemical hazes have the potential to drastically alter temperature structure and atmospheric circulation of these planets but have previously been neglected in general circulation models (GCMs). We present GCM simulations of hot Jupiter HD 189733b that include photochemical hazes as a radiatively active tracer fully coupled to atmospheric dynamics. The influence of haze radiative feedback strongly depends on the assumed haze optical properties. For soot hazes, two distinct thermal inversions form, separated by a local temperature minimum around 10$^{-5}$ bar caused by upwelling on the dayside mixing air with low haze abundance upwards. The equatorial jet broadens and slows down. The horizontal distribution of hazes remains relatively similar to simulations with radiatively passive tracers. For Titan-type hazes, the equatorial jet accelerates and extends to much lower pressures, resulting in a dramatically different 3D distribution of hazes compared to radiatively passive or soot hazes. Further experimental and observational studies to constrain the optical properties of photochemical hazes will therefore be crucial for understanding the role of hazes in exoplanet atmospheres. In the dayside emission spectrum, for both types of hazes the amplitude of near-infrared features is reduced, while the emitted flux at longer wavelengths ($>$4 $\mu$m) increases. Haze radiative feedback leads to increased phase curve amplitudes in many infrared wavelength regions, mostly due to stronger dayside emission.

We calculate the strange star properties in the framework of the Field Correlator Method. We find that for the values of the gluon condensate $G_2=0.006\;{\rm GeV}^4$ and $G_2=0.0068\;{\rm GeV}^4$, which give a critical temperature $T_c\sim170\;{\rm MeV}$ at $\mu_c=0$, the sequences of strange stars are compatible with some of the semi-empirical mass-radius relations and data obtained from astrophysical observations.

We investigate general aspects of the mass defects of strange stars in the context of the Field Correlator Method, without magnetic field. The main parameters of the model that enter the corresponding nonperturbative equation of state of the quark gluon plasma are the gluon condensate $G_2$ and the large distance static $Q{\bar Q}$ potential $V_1$. We calculate mass defects of stellar configurations in the central density range $11<\log\rho_c<18$. In general, the mass defects are strongly dependent on the model parameters. For a large range of values of $G_2$ and $V_1$, we obtain anomalous mass defects with magnitudes around $10^{53}\,$erg\,, of the same order of the observed energies of gamma-ray bursts and neutrino emissions in SN1987A, and of the theoretically predicted energies of the quark-novae explosions.

We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that requires the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss and how the activation functions of the neural model are dynamically accounted for. The performance of our GA-PINN model is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with our GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative $l^{2}$ error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ``ground truth'', and with the error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by our GA-PINN model is comparable to that obtained with a shock-capturing scheme and significantly higher than that achieved by a baseline PINN algorithm. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved ones, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN method discussed could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions.

Some analyses of the third gravitational wave catalogue released by the LIGO-Virgo-KAGRA collaboration (LVK) suggest an excess of black holes around $15-20 M_{\odot}$. In order to investigate this feature, we introduce two flexible population models, a semi-parametric one and a non-parametric one. Both make use of reversible jump Markov chain Monte-Carlo to optimise their complexity. We also illustrate how the latter can be used to efficiently perform model selection. Our parametric model broadly agrees with the fiducial analysis of the LVK, but finds a peak of events at slightly larger masses. Our non-parametric model shows this same displacement. Moreover, it also suggests the existence of an excess of black holes around $20 M_{\odot}$. We assess the robustness of this prediction by performing mock injections and running hierarchical analyses on those. We find that such a feature might be due to statistical fluctuations, given the small number of events observed so far, with a $5\%$ probability. We estimate that with a few hundreds of observations, as expected for O4, our non-parametric model will, be able to robustly determine the presence of this excess. It will then allow for an efficient agnostic inference of the properties of black holes.

In this essay it is proved that there are black holes that are dangerously cold. In particular, by analyzing the emission spectra of highly charged black holes we reveal the fact that near-extremal black holes whose Bekenstein-Hawking temperatures lie in the regime $T_{\text{BH}}\lesssim m^6_e/e^3$ may turn into horizonless naked singularities, thus violating the cosmic censorship principle, if they emit a photon with the characteristic thermal energy $\omega=O(T_{\text{BH}})$ [here $\{m_e,e\}$ are respectively the proper mass and the electric charge of the electron, the lightest charged particle]. We therefore raise here the conjecture that, in the yet unknown quantum theory of gravity, the temperatures of well behaved black-hole spacetimes are fundamentally bounded from below by the relation $T_{\text{BH}}\gtrsim m^6_e/e^3$.

The study of spinning bodies moving in curved spacetime has relevance to binary black hole systems with large mass ratios, as well as being of formal interest. At zeroth order in a binary's mass ratio, the smaller body moves on a geodesic of the larger body's spacetime. Post-geodesic corrections describing forces driving the small body's worldline away from geodesics must be incorporated to model the system accurately. An important post-geodesic effect is the gravitational self-force, which describes the small body's interaction with its own spacetime curvature. This effect includes the backreaction due to gravitational-wave emission that leads to the inspiral of the small body into the black hole. When a spinning body orbits a black hole, its spin couples to spacetime curvature. This introduces another post-geodesic correction known as the spin-curvature force. An osculating geodesic integrator that includes both the backreaction due to gravitational waves and spin-curvature forces can be used to generate a spinning-body inspiral. In this paper, we use an osculating geodesic integrator to combine the leading backreaction of gravitational waves with the spin-curvature force. Our analysis only includes the leading orbit-averaged dissipative backreaction, and examines the spin-curvature force to leading order in the small body's spin. This is sufficient to build generic inspirals of spinning bodies, and serves as a foundation for further work examining how to include secondary spin in large-mass-ratio waveform models.

Direct detection experiments are still one of the most promising ways to unravel the nature of dark matter. To fully understand how well these experiments constrain the dark matter interactions with the Standard Model particles, all the uncertainties affecting the calculations must be known. It is especially critical now because direct detection experiments recently moved from placing limits only on the two elementary spin independent and spin dependent operators to the complete set of possible operators coupling dark matter and nuclei in non-relativistic theory. In our work, we estimate the effect of nuclear configuration-interaction uncertainties on the exclusion bounds for one of the existing xenon-based experiments for all fifteen operators. We find that for operator number 13 the $\pm1\sigma$ uncertainty on the coupling between the dark matter and nucleon can reach more than 50% for dark matter masses between 10 and 1000 GeV. In addition, we discuss how quantum computers can help to reduce this uncertainty.

The Fisher information matrix is a quantity of fundamental importance for information geometry and asymptotic statistics. In practice, it is widely used to quickly estimate the expected information available in a data set and guide experimental design choices. In many modern applications, it is intractable to analytically compute the Fisher information and Monte Carlo methods are used instead. The standard Monte Carlo method produces estimates of the Fisher information that can be biased when the Monte-Carlo noise is non-negligible. Most problematic is noise in the derivatives as this leads to an overestimation of the available constraining power, given by the inverse Fisher information. In this work we find another simple estimate that is oppositely biased and produces an underestimate of the constraining power. This estimator can either be used to give approximate bounds on the parameter constraints or can be combined with the standard estimator to give improved, approximately unbiased estimates. Both the alternative and the combined estimators are asymptotically unbiased so can be also used as a convergence check of the standard approach. We discuss potential limitations of these estimators and provide methods to assess their reliability. These methods accelerate the convergence of Fisher forecasts, as unbiased estimates can be achieved with fewer Monte Carlo samples, and so can be used to reduce the simulated data set size by several orders of magnitude.

We offer a survey of the matter-antimatter evolution within the primordial Universe. While the origin of the tiny matter-antimatter asymmetry has remained one of the big questions in modern cosmology, antimatter itself has played a large role for much of the Universe's early history. In our study of the evolution of the Universe we adopt the position of the standard model $\Lambda$-CDM Universe implementing the known baryonic asymmetry. We present the composition of the Universe across its temperature history while emphasizing the epochs where antimatter content is essential to our understanding. Special topics we address include the heavy quarks in quark-gluon plasma (QGP), the creation of matter from QGP, the free-streaming of the neutrinos, the vanishing of the muons, the magnetism in the electron-positron cosmos, and a better understanding of the environment of the Big Bang Nucleosynthesis (BBN) producing the light elements. We suggest but do not explore further that the methods used in exploring the early Universe may also provide new insights in the study of exotic stellar cores, magnetars, as well as gamma-ray burst (GRB) events. We describe future investigations required in pushing known physics to its extremes in the unique laboratory of the matter-antimatter early Universe.

The Effective Field Theory (EFT) of perturbations on an arbitrary background geometry with a timelike scalar profile has been recently constructed in the context of scalar-tensor theories. Unlike General Relativity, the regular Hayward metric is realized as an exact background metric in the Effective Field Theory with timelike scalar profile without resorting to special matter field, such as nonlinear electrodynamics. The fundamental quasinormal mode for axial graviational perturbations of this black hole has been considered recently with the help of various methods. Here we make a further step in this direction and find that, unlike the fundamental mode, a few first overtones deviate from their Schwarzschild limit at a much higher rate. This outburst of overtones occurs because the overtones are extremely sensitive to the least change of the near-horizon geometry. The analytical formula for quasinormal modes is obtained in the eikonal regime. In addition, we calculated grey-body factors and showed that regular Hayward black hole with a scalar hair has smaller grey-body factor than the Schwarzschild one. Integration of the wave-like equation in time-domain shows that the power-law tails following the ring-down phase at late times are indistinguishable from the Schwarzschild ones.

The merits of a perturbation theory based on a mean to osculating transformation that is pure periodic in the fast angle are investigated. The exact separation of the purely short-period effects of the perturbed Keplerian dynamics from the long-period mean frequencies is achieved by a non-canonical transformation, which, therefore, cannot be computed by Hamiltonian methods. For this case, the evolution of the mean elements strictly adheres to the average behavior of the osculating orbit. However, due to the inescapable truncation of perturbation solutions, the fact that this theory confines the long-period variations of the semimajor axis into the mean variation equations, how tiny they may be, can have adverse effects in the accuracy of long-term semi-analytic propagations based on it

We study the Hamiltonian structure of unimodular-like theories, where the cosmological constant (or other supposed constants of nature) are demoted from fixed parameters to classical constants of motion. No new local degrees of freedom are present as a result of a $U(1)$ gauge invariance of the theory. Hamiltonian analysis of the action reveals that the only possible gauge fixing that can be enforced is setting the spatial components of the four-volume time vector ${\cal T}^{i}\approx0$. As a consequence of this, the gauge-fixed unimodular path integral is equivalent to the minisuperspace unimodular path integral. However, should we break the $U(1)$ gauge invariance, two things happen: a massless propagating degree of freedom appears, and the (gauge-invariant) zero-mode receives modified dynamics. The implications are investigated, with the phenomenology depending crucially on the target ``constant''.

(1) Purpose: Conditions of formation of compound nuclear system needed for synthesis of heavy nuclei in pycnonuclear reactions in compact stars are studied on a quantum mechanical basis. (2) Methods: Method of multiple internal reflections is generalized for pycnoreactions in compact stars with new calculations of quasibound spectra and spectra of zero-point vibrations. (3) Results: Peculiarities of the method are analyzed for reaction with isotopes of Carbon. The developed method takes into account continuity and conservation of quantum flux (describing pycnonuclear reaction) inside the full spacial region of reaction including nuclear region. This gives appearance of new states (called as quasibound states), in which compound nuclear systems of Magnesium are formed with the largest probability. These states have not been studied yet in synthesis of elements in stars. Energy spectra of zero-point vibrations and spectra of quasibound states are estimated with high precision for reactions with isotopes of Carbon. At the first time influence of plasma screening on quasibound states and states of zero-point vibrations in pycnonuclear reactions has been studied. (4) Conclusion: The probability of formation of compound nuclear system in quasibound states in pycnonuclear reaction is essentially larger than the probability of formation of this system in states of zero-point vibrations studied by Zel'dovich and followers. So, synthesis of Magnesium from isotopes of Carbon is more probable through the quasibound states than through the states of zero-point vibrations in compact stars. Energy spectra of zero-point vibrations are changed essentially after taking plasma screening into account. Analysis shows that from all studied isotopes of Magnesium only \isotope[24]{Mg} is stable after synthesis at energy of relative motion of 4.881~MeV of incident nuclei \isotope[12]{C}.

Mercury is the closest planet to the Sun, possesses a weak intrinsic magnetic field and has only a very tenuous atmosphere (exosphere). These three conditions result in a direct coupling between the plasma emitted from the Sun (namely the solar wind) and Mercury's surface. The planet's magnetic field leads to a non-trivial pattern of plasma precipitation onto the surface, that is expected to contribute to the alteration of the regolith over geological time scales. The goal of this work is to study the solar wind plasma precipitation onto the surface of Mercury from a geographical perspective, as opposed to the local-time-of-day approach of previous precipitation modeling studies. We employ solar wind precipitation maps for protons and electrons from two fully-kinetic numerical simulations of Mercury's plasma environment. These maps are then integrated over two full Mercury orbits (176 Earth days). We found that the plasma precipitation pattern at the surface is most strongly affected by the upstream solar wind conditions, particularly by the interplanetary magnetic field direction, and less by Mercury's 3:2 spin-orbit resonance. We also found that Mercury's magnetic field is able to shield the surface from roughly 90% of the incoming solar wind flux. At the surface, protons have a broad energy distribution from below 500 eV to more than 1.5 keV; while electrons are mostly found in the range 0.1-4 keV. These results will help to better constrain space weathering and exosphere source processes at Mercury, as well as to interpret observations by the ongoing ESA/JAXA BepiColombo mission.

We present simulation results of a gradual solar energetic particle (SEP) event detected on 2021 October 9 by multiple spacecraft, including BepiColombo (Bepi) and near-Earth spacecraft such as the Advanced Composition Explorer (ACE). A peculiarity of this event is that the presence of a high speed stream (HSS) affected the low-energy ion component ($\lesssim 5$ MeV) of the gradual SEP event at both Bepi and ACE, despite the HSS having only a modest solar wind speed increase. Using the EUHFORIA (European Heliospheric FORecasting Information Asset) magnetohydrodynamic model, we replicate the solar wind during the event and the coronal mass ejection (CME) that generated it. We then combine these results with the energetic particle transport model PARADISE (PArticle Radiation Asset Directed at Interplanetary Space Exploration). We find that the structure of the CME-driven shock was affected by the non-uniform solar wind, especially near the HSS, resulting in a shock wavefront with strong variations in its properties such as its compression ratio and obliquity. By scaling the emission of energetic particles from the shock to the solar wind compression at the shock, an excellent match between the PARADISE simulation and in-situ measurements of $\lesssim 5$ MeV ions is obtained. Our modelling shows that the intricate intensity variations observed at both ACE and Bepi were influenced by the non-uniform emission of energetic particles from the deformed shock wave and demonstrates the influence of even modest background solar wind structures on the development of SEP events.