Papers for Tuesday, Nov 14 2023

Studies of the Universe's transition to smoothness in the context of LCDM have all pointed to a transition radius no larger than ~300 Mpc. These are based on a broad array of tracers for the matter power spectrum, including galaxies, clusters, quasars, the Ly-alpha forest and anisotropies in the cosmic microwave background. It is therefore surprising, if not anomalous, to find many structures extending out over scales as large as ~2 Gpc, roughly an order of magnitude greater than expected. Such a disparity suggests that new physics may be contributing to the formation of large-scale structure, warranting a consideration of the alternative FLRW cosmology known as the $R_h=ct$ universe. This model has successfully eliminated many other problems in LCDM. In this paper, we calculate the fractal (or Hausdorff) dimension in this cosmology as a function of distance, showing a transition to smoothness at ~2.2 Gpc, fully accommodating all of the giant structures seen thus far. This outcome adds further observational support for R_h=ct over the standard model.

The common-envelope (CE) phase is a crucial stage in binary star evolution because the orbital separation can shrink drastically while ejecting the envelope of a giant star. Three-dimensional (3D) hydrodynamic simulations of CE evolution are indispensable to learning about the mechanisms that play a role during the CE phase. While these simulations offer great insight, they are computationally expensive. We propose a one-dimensional (1D) model to simulate the CE phase within the stellar evolution code $\texttt{MESA}$ by using a parametric drag force prescription for dynamical drag and adding the released orbital energy as heat into the envelope. We compute CE events of a $0.97\,\mathrm{M}_\odot$ asymptotic giant-branch star and a point mass companion with mass ratios of 0.25, 0.50, and 0.75, and compare them to 3D simulations of the same setup. The 1D CE model contains two free parameters, which we demonstrate are both needed to fit the spiral-in behavior and the fraction of ejected envelope mass of the 1D method to the 3D simulations. For mass ratios of 0.25 and 0.50, we find good-fitting 1D simulations, while for a mass ratio of 0.75, we do not find a satisfactory fit to the 3D simulation as some of the assumptions in the 1D method are no longer valid. In all our simulations, we find that the released recombination energy is important to accelerate the envelope and drive the ejection.

We analyze the physical properties of the gaseous intracluster medium (ICM) at the center of massive galaxy clusters with TNG-Cluster, a new cosmological magnetohydrodynamical simulation. Our sample contains 352 simulated clusters spanning a halo mass range of $10^{14} < {\rm M}_{\rm 500c} / M_\odot < 2 \times 10^{15}$ at $z=0$. We focus on the proposed classification of clusters into cool-core (CC) and non-cool-core (NCC) populations, the $z=0$ distribution of cluster central ICM properties, and the redshift evolution of the CC cluster population. We analyze resolved structure and radial profiles of entropy, temperature, electron number density, and pressure. To distinguish between CC and NCC clusters, we consider several criteria: central cooling time, central entropy, central density, X-ray concentration parameter, and density profile slope. According to TNG-Cluster and with no a-priori cluster selection, the distributions of these properties are unimodal, whereby CCs and NCCs represent the two extremes. Across the entire TNG-Cluster sample at $z=0$ and based on central cooling time, the strong CC fraction is $f_{\rm SCC} = 24\%$, compared to $f_{\rm WCC} = 60\% $ and $f_{\rm NCC} = 16\%$ for weak and non-cool-cores, respectively. However, the fraction of CCs depends strongly on both halo mass and redshift, although the magnitude and even direction of the trends vary with definition. The abundant statistics of simulated high-mass clusters in TNG-Cluster enables us to match observational samples and make a comparison with data. The CC fractions from $z=0$ to $z=2$ are in broad agreement with observations, as are radial profiles of thermodynamical quantities, globally as well as split for CC versus NCC halos. TNG-Cluster can therefore be used as a laboratory to study the evolution and transformations of cluster cores due to mergers, AGN feedback, and other physical processes.

The intracluster medium (ICM) of galaxy clusters encodes the impact of the physical processes that shape these massive halos, including feedback from central supermassive black holes (SMBHs). In this study we examine the gas thermodynamics, kinematics, and the effects of SMBH feedback on the core of Perseus-like galaxy clusters with a new simulation suite: TNG-Cluster. We first make a selection of simulated clusters similar to Perseus based on total mass and inner ICM properties, i.e. cool-core nature. We identify 30 Perseus-like systems among the 352 TNG-Cluster halos at $z=0$. Many exhibit thermodynamical profiles and X-ray morphologies with disturbed features such as ripples, bubbles and shock fronts that are qualitatively similar to X-ray observations of Perseus. To study observable gas motions, we generate XRISM mock X-ray observations and conduct a spectral analysis of the synthetic data. In agreement with existing Hitomi measurements, TNG-Cluster predicts subsonic gas turbulence in the central regions of Perseus-like clusters, with a typical line-of-sight velocity dispersion of 200 km/s. This implies that turbulent pressure contributes $< 10\%$ to the dominant thermal pressure. In TNG-Cluster, such low (inferred) values of ICM velocity dispersion coexist with high-velocity outflows and bulk motions of relatively small amounts of super-virial hot gas, moving up to thousands of km/s. However, detecting these outflows observationally may prove challenging due to their anisotropic nature and projection effects. Driven by SMBH feedback, such outflows are responsible for many morphological disturbances in the X-ray maps of cluster cores. They also increase both the inferred, and intrinsic, ICM velocity dispersion. This effect is somewhat stronger when velocity dispersion is measured from higher-energy lines.

The most massive black holes in our Universe form binaries at the centre of merging galaxies. The recent evidence for a gravitational-wave (GW) background from pulsar timing may constitute the first observation that these supermassive black hole binaries (SMBHBs) merge. Yet, the most massive SMBHBs are out of reach of interferometric detectors and are exceedingly difficult to resolve individually with pulsar timing. These limitations call for unexplored strategies to detect individual SMBHBs in the uncharted frequency band $\lesssim10^{-5}\,\rm Hz$ in order to establish their abundance and decipher the coevolution with their host galaxies. Here we show that SMBHBs imprint detectable long-term modulations on GWs from stellar-mass binaries residing in the same galaxy. We determine that proposed deci-Hz GW interferometers sensitive to numerous stellar-mass binaries will uncover modulations from $\sim\mathcal{O}(10^{-1}$ - $10^4)$ SMBHBs with masses $\sim\mathcal{O}(10^7$ - $10^9)\,\rm M_\odot$ out to redshift $z\sim3.5$. This offers a unique opportunity to map the population of SMBHBs through cosmic time, which might remain inaccessible otherwise.

The 30 Doradus region in the Large Magellanic Cloud (LMC) is the most energetic star-forming region in the Local Group. It is powered by the feedback from the massive stars in R136, the 1-2 Myr old central massive cluster. 30 Doradus has therefore long been regarded as a laboratory for studying star and star cluster formation under conditions reminiscent of the early Universe. We use JWST NIRCam observations to analyse how star formation proceeds in the region. Using selections based on theoretical isochrones on colour-magnitude diagrams, we identify populations of different ages. We select pre-main-sequence (PMS) stars and young stellar objects that show excess emission from warm dust or emission lines. Studying the spatial distribution of the different populations, we find that the youngest PMS stars with ages < 0.5 Myr are located in an elongated structure that stretches towards the north-east from the central cluster. The same structure is found in the sources that show an infrared excess, appears to be overlapping with cold molecular gas, and covers previously investigated sites of ongoing star formation. Pre-main-sequence stars with ages between 1 and 4 Myr and upper main-sequence stars are concentrated in the centre of R136, while older stars are more uniformly distributed across the field and likely belong to the LMC field population. Nonetheless, we find stars with excess emission from on dust or emission lines as far as 100 pc from the centre, indicating extended recent star formation. We interpret the elongated structure formed by the youngest PMS stars to be an indication of the still-ongoing hierarchical assembly of the R136 cluster. Additionally, the lower density of old PMS stars with emission due to ongoing accretion in the central region suggests that feedback from the R136 stars is effective in disrupting the disks of PMS stars.

The most massive galaxy clusters in the Universe host tens to hundreds of massive satellite galaxies, but it is unclear if these satellites are able to retain their own gaseous atmospheres. We analyze the evolution of $\sim90,000$ satellites of stellar mass $\sim10^{9-12.5}\,M_\odot$ around 352 galaxy clusters of mass $M_{\rm vir}\sim10^{14.3-15.4}\,M_\odot$ at $z=0$ from the new TNG-Cluster suite of cosmological magneto-hydrodynamical galaxy cluster simulations. The number of satellites per host increases with host mass, and the mass--richness relation broadly agrees with observations. A halo of mass $M_{\rm 200c}\sim10^{14.5}\,(10^{15})\,M_\odot$ hosts $\sim100\,(300)$ satellites today. Only a minority of satellites retain some gas, and this fraction increases with stellar mass. Lower mass satellites $\sim10^{9-10}\,M_\odot$ are more likely to retain part of their cold interstellar medium, consistent with ram pressure preferentially removing hot extended gas first. At higher stellar masses $\sim10^{10.5-12.5}\,M_\odot$ the fraction of gas-rich satellites increases to unity, and nearly all satellites retain a sizeable portion of their hot, spatially-extended circumgalactic medium (CGM), despite the ejective activity of their supermassive black holes. According to TNG-Cluster, the CGM of these gaseous satellites can be seen in soft X-ray emission (0.5-2.0 keV) that is $\gtrsim10$ times brighter than the local background. This X-ray surface brightness excess around satellites extends to $\sim30-100$ kpc, and is strongest for galaxies with higher stellar masses and larger host-centric distances. Approximately 10 per cent of the soft X-ray emission in cluster outskirts $\sim0.75-1.5R_{\rm 200c}$ originates from satellites. The CGM of member galaxies reflects the dynamics of cluster-satellite interactions and contributes to the observationally-inferred properties of the intracluster medium.

We introduce the new TNG-Cluster project, an addition to the IllustrisTNG suite of cosmological magnetohydrodynamical simulations of galaxy formation. Our objective is to significantly increase the statistical sampling of the most massive and rare objects in the Universe: galaxy clusters with log(M_200c / Msun) > 14.3 - 15.4 at z=0. To do so, we re-simulate 352 cluster regions drawn from a 1 Gpc volume, thirty-six times larger than TNG300, keeping entirely fixed the IllustrisTNG physical model as well as the numerical resolution. This new sample of hundreds of massive galaxy clusters enables studies of the assembly of high-mass ellipticals and their supermassive black holes (SMBHs), brightest cluster galaxies (BCGs), satellite galaxy evolution and environmental processes, jellyfish galaxies, intracluster medium (ICM) properties, cooling and active galactic nuclei (AGN) feedback, mergers and relaxedness, magnetic field amplification, chemical enrichment, and the galaxy-halo connection at the high-mass end, with observables from the optical to radio synchrotron and the Sunyaev-Zeldovich (SZ) effect, to X-ray emission, as well as their cosmological applications. We present an overview of the simulation, the cluster sample, selected comparisons to data, and a first look at the diversity and physical properties of our simulated clusters and their hot ICM.

Galaxy clusters are unique laboratories for studying astrophysical processes and their impact on gas kinematics. Despite their importance, the full complexity of gas motion within and around clusters remains poorly known. This paper is part of a series presenting first results from the new TNG-Cluster simulation, a suite of 352 massive clusters including the full cosmological context, mergers, accretion, baryonic processes, feedback, and magnetic fields. Studying the dynamics and coherence of gas flows, we find that gas motions in cluster cores and intermediate regions are largely balanced between inflows and outflows, exhibiting a Gaussian distribution centered at zero velocity. In the outskirts, even the net velocity distribution becomes asymmetric, featuring a double peak where the second peak reflects cosmic accretion. Across all cluster regions, the resulting net flow distribution reveals complex gas dynamics. These are strongly correlated with halo properties: at a given total cluster mass, unrelaxed, late-forming halos with less massive black holes and lower accretion rates exhibit a more dynamic behavior. Our analysis shows no clear relationship between line-of-sight and radial gas velocities, suggesting that line-of-sight velocity alone is insufficient to distinguish between inflowing and outflowing gas. Additional properties, such as temperature, can help break this degeneracy. A velocity structure function (VSF) analysis indicates more coherent gas motion in the outskirts and more disturbed kinematics towards halo centers. In all cluster regions, the VSF shows a slope close to the theoretical models of Kolmogorov (1/3), except within 50 kpc of the cluster cores, where the slope is significantly steeper. The outcome of TNG-Cluster broadly aligns with observations of the VSF of multiphase gas across different scales in galaxy clusters, ranging from 1 kpc to Megaparsec scales.

Radio relics are diffuse synchrotron sources in the outskirts of merging galaxy clusters energized by the merger shocks. In this paper, we present an overview of the radio relics in massive cluster mergers identified in the new TNG-Cluster simulation. This is a suite of magnetohydrodynamical cosmological zoom-in simulations of 352 massive galaxy clusters with $M_{\rm 500c}= 10^{14.0-15.3}\rm~M_{\odot}$ sampled from a 1 Gpc-size cosmological box. The simulations are performed using the moving-mesh code AREPO with the galaxy formation model and high numerical resolution consistent with the TNG300 run of the IllustrisTNG series. We post-process the shock properties obtained from the on-the-fly shock finder to estimate the diffuse radio emission generated by cosmological shockwaves for a total of $\sim300$ radio relics at redshift $z=0-1$. TNG-Cluster returns a variety of radio relics with diverse morphologies, encompassing textbook examples of double radio relics, single relics, and ``inverted" radio relics that are convex to the cluster center. Moreover, the simulated radio relics reproduce both the abundance and statistical relations of observed relics. We find that extremely large radio relics ($>$ 2 Mpc) are predominantly produced in massive cluster mergers with $M_{\rm 500c}\gtrsim8\times10^{14}~\rm~M_{\odot}$. This underscores the significance of simulating massive mergers to study giant radio relics similar to those found in observations. We release a library of radio relics from the TNG-Cluster simulation, which will serve as a crucial reference for upcoming next-generation surveys.

Measurements of the mean free path of Lyman-continuum photons in the intergalactic medium during the epoch of reionization can help constrain the nature of the sources as well as sinks of hydrogen-ionizing radiation. A recent approach to this measurement has been to utilize composite spectra of multiple quasars at $z\sim 6$, and infer the mean free path after correcting the spectra for the presence of quasar proximity zones. This has revealed not only a steep drop in the mean free path from $z=5$ to $z=6$, but also potentially a mild tension with reionization simulations. We critically examine such direct measurements of the mean free path for biases due to quasar environment, incomplete reionization, and quasar proximity zones. Using cosmological radiative transfer simulations of reionization combined with one-dimensional radiative transfer calculations of quasar proximity zones, we find that the bias in the mean free path due to overdensities around quasars is minimal at $z\sim 6$. Patchiness of reionization at this redshift also does not affect the measurements significantly. Fitting our model to the data results in a mean free path of $\lambda_{\mathrm{mfp}}=0.90^{+0.66}_{-0.40}$ pMpc at $z=6$, which is consistent with the recent measurements in the literature, indicating robustness with respect to the modelling of quasar proximity zones. We also compare various ways in which the mean free path has been defined in simulations before the end of reionization. Overall, our finding is that recent measurements of the mean free path appear to be robust relative to several sources of potential bias.

The thermal freeze-out mechanism in its classical form is tightly connected to physics beyond the Standard Model around the electroweak scale, which has been the target of enormous experimental efforts. In this work we study a dark matter model in which freeze-out is triggered by a strong first-order phase transition in a dark sector, and show that this phase transition must also happen close to the electroweak scale, i.e. in the temperature range relevant for gravitational wave searches with the LISA mission. Specifically, we consider the spontaneous breaking of a $U(1)^\prime$ gauge symmetry through the vacuum expectation value of a scalar field, which generates the mass of a fermionic dark matter candidate that subsequently annihilates into dark Higgs and gauge bosons. In this set-up the peak frequency of the gravitational wave background is tightly correlated with the dark matter relic abundance, and imposing the observed value for the latter implies that the former must lie in the milli-Hertz range. A peculiar feature of our set-up is that the dark sector is not necessarily in thermal equilibrium with the Standard Model during the phase transition, and hence the temperatures of the two sectors evolve independently. Nevertheless, the requirement that the universe does not enter an extended period of matter domination after the phase transition, which would strongly dilute any gravitational wave signal, places a lower bound on the portal coupling that governs the entropy transfer between the two sectors. As a result, the predictions for the peak frequency of gravitational waves in the LISA band are robust, while the amplitude can change depending on the initial dark sector temperature.

Recent measurements of the ionizing photon mean free path (MFP) based on composite quasar spectra may point to a late end to reionization at $z<6$. These measurements are challenging, however, because they rely on assumptions about the proximity zones of the quasars in the analysis. For example, some of the $z\sim 6$ quasars in the composite might have been close to large-scale regions where reionization was still ongoing ("neutral islands"), and it is unclear how this would affect the measurements. We address the question here with mock MFP measurements from radiative transfer simulations. We find that, even in the presence of neutral islands, the inferred MFP tracks to within $30 \%$ the true attenuation scale of the spatially averaged IGM, which includes opacity from both the ionized medium and the islands. During reionization, this scale is always shorter than the MFP in the ionized medium. The inferred MFP is sensitive at the $< 50\%$ level to assumptions about the quasar environments and lifetimes for realistic models. We demonstrate that future analyses with improved data may require explicitly modeling the effects of neutral islands on the composite spectra, and we outline a method for doing this. Lastly, we quantify the effects of neutral islands on Lyman-series transmission, which has been modeled with optically thin simulations in previous MFP analyses. Neutral islands can suppress transmission at $\lambda_{\rm rest} < 912$ \r{A} significantly, up to a factor of 2 for $z_{\rm qso} = 6$ in a plausible reionization scenario, owing to absorption by many closely spaced lines as quasar light redshifts into resonance. However, the suppression is almost entirely degenerate with the spectrum normalization, thus does not significantly bias the inferred MFP.

Closure invariants in interferometry carry calibration-independent information about the morphology of an observed object. Excepting simple cases, a mapping between closure invariants and morphologies is not well established. We aim to demonstrate that closure invariants can be used to classify the morphology and estimate the morphological parameters using simple Machine Learning models. We consider 6 morphological classes -- point-like, uniform circular disc, crescent, dual disc, crescent with elliptical accretion disc, and crescent with double jet lobes -- described by phenomenological parameters. Using simple logistic regression, multi-layer perceptron (MLP), convolutional neural network, and random forest models on closure invariants obtained from a sparse aperture coverage, we find that all models except logistic regression are able to classify the morphology with an $F_1$ score $\gtrsim 0.8$. The classification accuracy notably improves with greater aperture coverage. We also estimate morphological parameters of uniform circular disc, crescent, and dual disc using simple MLP models, and perform a parametric image reconstruction. The reconstructed images do not retain information about absolute position or intensity scale. The estimated parameters and reconstructed images are found to correspond well with the inputs. However, the prediction accuracy worsens with increasing morphological complexity. This proof-of-concept method opens an independent approach to interferometric imaging under challenging observing conditions such as that faced by the Event Horizon Telescope and Very Long Baseline Interferometry in general, and can complement other methods to robustly constrain an object's morphology.

The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared dust polarization (FIR) polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the FIR total and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism.

The Imaging X-ray Polarimetry Explorer measured with high significance the X-ray polarization of the brightest Z-source Sco X-1, resulting in the nominal 2-8 keV energy band in a polarization degree of 1.0(0.2)% and a polarization angle of 8(6){\deg} at 90% of confidence level. This observation was strictly simultaneous with observations performed by NICER, NuSTAR, and Insight-HXMT, which allowed for a precise characterization of its broad-band spectrum from soft to hard X-rays. The source has been observed mainly in its soft state, with short periods of flaring. We also observed low-frequency quasi-periodic oscillations. From a spectro-polarimetric analysis, we associate a polarization to the accretion disk at <3.2% at 90% of confidence level, compatible with expectations for an electron-scattering dominated optically thick atmosphere at the Sco X-1 inclination of 44{\deg}; for the higher-energy Comptonized component, we obtain a polarization of 1.3(0.4)%, in agreement with expectations for a slab of Thomson optical depth of ~7 and an electron temperature of ~3 keV. A polarization rotation with respect to previous observations by OSO-8 and PolarLight, and also with respect to the radio-jet position angle, is observed. This result may indicate a variation of the polarization with the source state that can be related to relativistic precession or to a change in the corona geometry with the accretion flow.

Superorbital periods that are observed in the brightness of Be/X-ray binaries may be driven by a misaligned and precessing Be star disc. We examine how the precessing disc model explains the superorbital variation of (i) the magnitude of the observed X-ray outbursts and (ii) the observed colour. With hydrodynamical simulations we show that the magnitude of the average accretion rate on to the neutron star, and therefore the X-ray outbursts, can vary by over an order of magnitude over the superorbital period for Be star spin-orbit misalignments $\gtrsim 70^\circ$ as a result of weak tidal truncation. Most Be/X-ray binaries are redder at optical maximum when the disc is viewed closest to face-on since the disc adds a large red component to the emission. However, A0538-66 is redder at optical minimum. This opposite behaviour requires an edge-on disc at optical minimum and a radially narrow disc such that it does not add a large red signature when viewed face-on. For A0538-66, the misalignment of the disc to the binary orbit must be about $70-80^\circ$ and the inclination of the binary orbit to the line of sight must be similarly high, although restricted to $<75^\circ$ by the absence of X-ray eclipses.

The double pulsar system, PSR J0737$-$3039A/B, consists of two neutron stars bound together in a highly relativistic orbit that is viewed nearly edge-on from the Earth. This alignment results in brief radio eclipses of the fast-rotating pulsar A when it passes behind the toroidal magnetosphere of the slow-rotating pulsar B. The morphology of these eclipses is strongly dependent on the geometric orientation and rotation phase of pulsar B, and their time-evolution can be used to constrain the geodetic precession rate of the pulsar. We demonstrate a Bayesian inference framework for modelling eclipse light-curves obtained with MeerKAT between 2019-2023. Using a hierarchical inference approach, we obtained a precession rate of $\Omega_{\rm SO}^{\rm B} = {5.16^{\circ}}^{+0.32^{\circ}}_{-0.34^{\circ}}$ yr$^{-1}$ for pulsar B, consistent with predictions from General Relativity to a relative uncertainty of 6.5%. This updated measurement provides a 6.1% test of relativistic spin-orbit coupling in the strong-field regime. We show that a simultaneous fit to all of our observed eclipses can in principle return a $\sim$1.5% test of spin-orbit coupling. However, systematic effects introduced by the current geometric orientation of pulsar B along with inconsistencies between the observed and predicted eclipse light curves result in difficult to quantify uncertainties. Assuming the validity of General Relativity, we definitively show that the spin-axis of pulsar B is misaligned from the total angular momentum vector by $40.6^{\circ} \pm 0.1^{\circ}$ and that the orbit of the system is inclined by approximately $90.5^{\circ}$ from the direction of our line of sight. Our measured geometry for pulsar B suggests the largely empty emission cone contains an elongated horseshoe shaped beam centered on the magnetic axis, and that it may not be re-detected as a radio pulsar until early-2035.

Particle-in-cell simulations have unveiled that shock-accelerated electrons do not follow a pure power-law distribution, but have an additional low-energy "thermal" part, which owns a considerable portion of the total energy of electrons. Investigating the effects of these thermal electrons on gamma-ray burst (GRB) afterglows may provide valuable insights into the particle acceleration mechanisms. We solve the continuity equation of electrons in the energy space, from which multi-wavelength afterglows are derived by incorporating processes including synchrotron radiation, synchrotron self-absorption, synchrotron self-Compton scattering, and gamma-gamma annihilation. First, there is an underlying positive correlation between temporal and spectral indices due to the cooling of electrons. Moreover, thermal electrons would result in the simultaneous non-monotonic variation in both spectral and temporal indices at multi-wavelength, which could be individually recorded by the 2.5-meter Wide Field Survey Telescope and Vera Rubin Observatory Legacy Survey of Space and Time (LSST). The thermal electrons could also be diagnosed from afterglow spectra by synergy observation in the optical (with LSST) and X-ray bands (with the Microchannel X-ray Telescope on board the Space Variable Objects Monitor). Finally, we use Monte Carlo simulations to obtain the distribution of peak flux ratio ($R_{\rm X}$) between soft and hard X-rays, and of the time delay ($\Delta t$) between peak times of soft X-ray and optical light curves. The thermal electrons significantly raise the upper limits of both $R_{\rm X}$ and $\Delta t$. Thus the distribution of GRB afterglows with thermal electrons is more dispersive in the $R_{\rm X} - \Delta t$ plane.

The present study aims to reinforce the evidence for the ~9 s pulsation in the gamma-ray binary LS 5039, derived with a Suzaku observation in 2007 and that with NuSTAR in 2016 (Yoneda et al 2000). Through a reanalysis of the NuSTAR data incorporating the orbital Doppler correction, the 9.0538 s pulsation was confirmed successfully even in the 3--10 keV range, where it was undetectable previously. This was attained by perceiving an energy-dependent drift in the pulse phase below 10 keV, and correcting the pulse timing of individual photons for that effect. Similarly, an archival 0.7--12 keV data set of LS 5039, taken with the ASCA GIS in 1999 October, was analyzed. The data showed possible periodicity at about 8.882 s, but again the energy-dependent phase drift was noticed below 10 keV. By correcting for this effect, and for the orbital Doppler delays in the LS 5039 system, the 2.8--12 keV periodicity became statistically significant at 8.891+- 0.001 s. The periods measured with ASCA, Suzaku, and NuSTAR approximately follow an average period derivative of dP/dt = 3.0 e-10 s/s. These results provide further evidence for the pulsation in this object, and strengthen the scenario by (Yoneda et al 2000), that the compact object in LS 5039 is a strongly magnetized neutron star.

Super Soft X-ray Sources (SSS) are white dwarf (WD) binaries that radiate almost entirely below $\sim$1~keV. Their X-ray spectra are often complex when viewed with the X-ray grating spectrometers, where numerous emission and absorption features are intermingled and hard to separate. The absorption features are mostly from the WD atmosphere, for which radiative transfer models have been constructed. The emission features are from the corona surrounding the WD atmosphere, in which incident emission from the WD surface is reprocessed. Modeling the corona requires different solvers and assumptions for the radiative transfer, which is yet to be achieved. We chose CAL87, a SSS in the Large Magellanic Cloud, which exhibits emission-dominated spectra from the corona as the WD atmosphere emission is assumed to be completely blocked by the accretion disk. We constructed a radiative transfer model for the corona using the two radiative transfer codes; xstar for a one-dimensional two-stream solver and MONACO for a three-dimensional Monte Carlo solver. We identified their differences and limitations in comparison to the spectra taken with the Reflection Grating Spectrometer onboard the XMM-Newton satellite. We finally obtained a sufficiently good spectral model of CAL87 based on the radiative transfer of the corona plus an additional collisionally ionized plasma. In the coming X-ray microcalorimeter era, it will be required to interpret spectra based on radiative transfer in a wider range of sources than what is presented here.

Purpose: The Hard X-ray Modulation Telescope is China's first X-ray astronomy satellite launched on June 15th, 2017, dubbed Insight-HXMT. Active and passive thermal control measures are employed to keep devices at suitable temperatures. In this paper, we analyzed the on-orbit thermal monitoring data of the first 5 years and investigated the effect of thermal deformation on the point spread function (PSF) of the telescopes. Methods: We examined the data of the on-orbit temperatures measured using 157 thermistors placed on the collimators, detectors and their support structures and compared the results with the thermal control requirements. The thermal deformation was evaluated by the relative orientation of the two star sensors installed on the main support structure. its effect was estimated with evolution of the PSF obtained with calibration scanning observations of the Crab nebula. Conclusion: The on-orbit temperatures met the thermal control requirements thus far, and the effect of thermal deformation on the PSF was negligible after the on-orbit pointing calibration.

Giant Low Surface Brightness Galaxies (GLSBGs) are fundamentally distinct from normal galaxies (LSBGs) in star formation and evolution. In this work, we collected 27 local GLSBGs. They have high stellar masses (M*>10^10 Msolar) and low SFRs. With the specific SFRs lower than the characteristic value of the local star-forming (SF) galaxies of M*=10^10 Msolar(sSFR < 0.1 Gyr^-1), GLSBGs deviate from the SF main sequence (MS) defined for local SFGs respectively by E07 and S16 at the high M* regime. They are HI-rich systems with HI gas mass fractions (fHI) higher than the S16 MS galaxies, but have little molecular gas (H2), implying a low efficiency of HI-to-H2 transition due to the low HI surface densities that are far lower than the minimum of 6 - 8 Msolar pc^-2 required for shielding the formed H2 from photodissociation. For GLSBGs, the inner, bulge-dominated part with lower SFRs and higher M* is the main force pulling the entire GLSBG off from the MS, while the outer, disk-dominated part with relatively higher SFRs and lower M* reduces the deviations from the MS. For some cases, the outer, disk-dominated parts even tend to follow the MS. In the aspect of NUV - r versus g - r colors, the outer, disk-dominated parts are blue and behave similarly to the normal star-forming galaxies while the inner, bulge-dominated parts are in statistics red, indicating an inside-out star formation mechanism for the GLSBGs. They show few signs of external interactions in morphology, excluding the recent major merger scenario.

We present the results obtained by analysing new Astrosat/UVIT far ultraviolet (FUV) image of the collisional-ring galaxy Cartwheel. The FUV emission is principally associated with the star-forming outer ring, with no UV detection from the nucleus and inner ring. A few sources are detected in the region between the inner and the outer rings, all of which lie along the spokes. The FUV fluxes from the detected sources are combined with aperture-matched multi-band photometric data from archival images to explore the post-collision star formation history of the Cartwheel. The data were corrected for extinction using Av derived from the Balmer decrement ratios and commonly used extinction curves. We find that the ring regions contain stellar populations of wide range of ages, with the bulk of the FUV emission coming from non-ionizing stars, formed over the last 20 to 150 Myr, that are ~25 times more massive than the ionizing populations. On the other hand, regions belonging to the spokes have negligible current star formation, with the age of the dominant older population systematically increasing as its distance from the outer ring increases. The presence of populations of a wide range of ages in the ring suggests that the stars formed in the wave in the past were dragged along it to the current position of the ring. We derive an average steady star formation rate, SFR=5 Msun/yr, over the past 150 Myr, with an increase to ~18 Msun/yr in the recent 10 Myr.

There have been significant developments in the period estimation tools and methods for analysing high energy pulsars in the past few decades. However, these tools lack well-standardised methods for calculating uncertainties in period estimation and other recovered parameters for Poisson--dominated data. Error estimation is important for assigning confidence intervals to the models we study, but due to their high computational cost, errors in the pulsar periods were largely ignored in the past. Furthermore, existing literature has often employed semi-analytical techniques that lack rigorous mathematical foundations or exhibit a predominant emphasis on the analysis of white noise and time series data. We present results from our numerical and analytical study of the error distribution of the recovered parameters of high energy pulsar data using the $Z_n^2$ method. We comprehensively formalise the measure of error for the generic pulsar period with much higher reliability than some common methods. Our error estimation method becomes more reliable and robust when observing pulsars for few kilo-seconds, especially for typical pulsars with periods ranging from a few milliseconds to a few seconds. We have verified our results with observations of the \emph{Crab} pulsar, as well as a large set of simulated pulsars. Our codes are publicly available for use.

The lateral density data obtained for different secondaries of an extensive air shower (EAS) from an array of detectors are usually described by some suitable lateral density functions (LDFs). Analyzing non-vertical simulated EASs generated with the CORSIKA code, it is found that the lateral and polar density distributions of electrons and muons are asymmetric in the ground plane. It means that typical expressions for symmetric lateral density functions (SLDFs) (\emph{e.g.} the Nishimura-Kamata-Greisen function) are inadequate to reconstruct the lateral and polar dependencies of such asymmetric electron or muon densities accurately. In order to provide a more consistent LDF for non-vertical shower reconstruction in the ground plane, the paper considers the issue of the modification of the SLDF analytically. The asymmetry arising from additional attenuation and correction of the positional coordinates (radial and polar) of cascade particles causes a gap length between the center of concentric equidensity ellipses and the EAS core. A toy function is introduced as a basic LDF to describe the asymmetric lateral and polar density distributions of electrons or muons of EASs, thereby predicting the gap length parameter. Consequently, the desired LDF describing the asymmetric density distributions of electrons and muons of EASs has emerged. We compare results from detailed simulations with the predictions of the analytical parametrization. The LDF derived in this work is found to be well-suited to reconstruct EASs in the ground plane directly.

Increased eccentricity of a black hole binary leads to reduced merger times. With n-body simulations and analytic approximations including the effects of general relativity (GR), we show that even a low mass companion orbiting a black hole binary can cause significant eccentricity oscillations of the binary as a result of the Kozai-Lidov mechanism. A companion with a mass as low as about 1% of the binary mass can drive the binary eccentricity up to >~ 0.8, while a mass of a few percent can drive eccentricities greater than 0.98. For low mass companions, this mechanism requires the companion to be on an orbit that is closer to retrograde than to prograde to the binary orbit and this may occur through capture of the third body. The effects of GR limit the radial range for the companion for which this mechanism works for the closest binaries. The merger timescale may be reduced by several orders of magnitude for a captured companion mass of only a few percent of the binary mass.

We delved into the assembly pathways and environments of compact groups (CGs) of galaxies using mock catalogues generated from semi-analytical models (SAMs) on the Millennium simulation. We investigate the ability of SAMs to replicate the observed CG environments and whether CGs with different assembly histories tend to inhabit specific cosmic environments. We also analyse whether the environment or the assembly history is more important in tailoring CG properties. We find that about half of the CGs in SAMs are non-embedded systems, 40% are inhabiting loose groups or nodes of filaments, while the rest distribute evenly in filaments and voids, in agreement with observations. We observe that early-assembled CGs preferentially inhabit large galaxy systems (~ 60%), while around 30% remain non-embedded. Conversely, lately-formed CGs exhibit the opposite trend. We also obtain that lately-formed CGs have lower velocity dispersions and smaller crossing times than early-formed CGs, but mainly because they are preferentially non-embedded. Those lately-formed CGs that inhabit large systems do not show the same features. Therefore, the environment plays a strong role in these properties for lately-formed CGs. Early-formed CGs are more evolved, displaying larger velocity dispersions, shorter crossing times, and more dominant first-ranked galaxies, regardless of the environment. Finally, the difference in brightness between the two brightest members of CGs is dependent only on the assembly history and not on the environment. CGs residing in diverse environments have undergone varied assembly processes, making them suitable for studying their evolution and the interplay of nature and nurture on their traits.

We present a start-to-end simulation aimed at studying the long-term fate of high mass X-ray binaries and whether a Thorne-\.Zytkow object (T\.ZO) might ultimately be produced. We analyze results from a 3D hydrodynamical simulation that models the eventual fate of LMC X-4, a compact high mass X-ray binary system, after the primary fills its Roche lobe and engulfs the neutron star companion. We discuss the outcome of this engulfment within the standard paradigm of T\.ZO formation. The post-merger angular momentum content of the stellar core is a key ingredient, as even a small amount of rotation can break spherical symmetry and produce a centrifugally supported accretion disk. Our findings suggest the inspiraling neutron star, upon merging with the core, can accrete efficiently via a disk at high rates ($\approx 10^{-2}M_\odot/{\rm s}$), subsequently collapsing into a black hole and triggering a bright transient with a luminosity and duration typical of an ultra-long gamma-ray burst. We propose that the canonical framework for T\.ZO formation via common envelope needs to be revised, as the significant post-merger accretion feedback will unavoidably unbind the vast majority of the surrounding envelope.

Unveiling the nature of progenitors is crucial for understanding the origin and the mechanism of core-collapse and thermonuclear supernovae (SNe). While several methods have been developed to derive stellar properties so far, many questions remain poorly understood. In this paper we demonstrate an observational approach to constrain progenitors of supernova remnants (SNRs) using abundances of carbon (C), nitrogen (N), and oxygen (O) in shock-heated circumstellar material (CSM). Our calculations with stellar evolution codes indicate that a total amount of these CNO elements will provide a more sensitive determination of the progenitor masses than the conventional method based on ejecta abundances. If the CNO lines (particularly those of C and N) are detected and measured their abundance ratios accurately, they can provide relatively robust constraint on the progenitor mass (and in some cases the rotation velocity) of SNRs. Since our method requires a better energy resolution and larger effective area in the soft X-ray band ($<1$~keV), XRISM launched on September 7, 2023 and next-generation microcalorimeter missions such as Athena, Lynx, LEM, and HUBS will bring a new insight into link between the progenitors and their remnants.

The Galclaim software is blue designed to identify association between astrophysical transient sources and host galaxy by computing the probability of chance alignment. It is distributed as an open source Python software. It is already used to identify, confirm or reject host galaxy candidates of GRBs and to validate or invalidate transient candidates in astrophysical observations. Such tools are also very useful to characterise archived transient candidates in large sky survey telescopes.

In this paper, we investigate the astrophysical processes of stellar-mass black holes (sMBHs) embedded in advection-dominated accretion flows (ADAFs) of supermassive black holes (SMBHs) in low-luminosity active galactic nuclei (AGNs). The sMBH is undergoing Bondi accretion at a rate lower than the SMBH. Outflows from the sMBH-ADAF dynamically interact with their surroundings and form a cavity inside the SMBH-ADAF, thereby quenching the accretion onto the SMBH. Rejuvenation of the Bondi accretion is rapidly done by turbulence. These processes give rise to quasi-periodic episodes of sMBH activities and create flickerings from relativistic jets developed by the Blandford-Znajek mechanism if the sMBH is maximally rotating. Accumulating successive sMBH-outflows trigger viscous instability of the SMBH-ADAF, leading to a flare following a series of flickerings. Recently, the similarity of near-infrared flare's orbits has been found by GRAVITY/VLTI astrometric observations of Sgr A$\!^{*}$: their loci during the last 4-years consist of a ring in agreement with the well-determined SMBH mass. We apply the present model to Sgr A$\!^{*}$, which shows quasi-periodic flickerings. A SMBHH of $\sim 40 M_{\odot}$ is preferred orbiting around the central SMBH of Sgr A$\!^{*}$ from fitting radio to X-ray continuum. Such an extreme mass ratio inspiraling (EMRI) provides an excellent laboratory for LISA, Taiji and Tianqin detection of mHz gravitational waves with strains of $\sim 10^{-17}$, as well as their polarization.

Strong iron lines are a common feature of the optical spectra of active galactic nuclei (AGNs) and quasars from $z\sim 6-7$ to the local Universe, and [Fe/Mg] ratios do not show cosmic evolution. During active episodes, accretion disks surrounding supermassive black holes (SMBHs) inevitably form stars in the self-gravitating part and these stars accrete with high accretion rates. In this paper, we investigate the population evolution of accretion-modified stars (AMSs) to produce irons and magnesium in AGNs. The AMSs as a new type of stars are allowed to have any metallicity but without significant loss from stellar winds since the winds are choked by the dense medium of the disks and return to the core stars. Mass functions of the AMS population show a pile-up or cutoff pile-up shape in top-heavy or top-dominant forms if the stellar winds are strong, consistent with the narrow range of supernovae (SN) explosions driven by the known pair-instability. This provides an efficient way to produce metals. Meanwhile, SN explosions support an inflated disk as a dusty torus. Furthermore, the evolving top-heavy initial mass functions (IMFs) lead to bright luminosity in infrared bands in dusty regions. This contributes a new component in infrared bands which is independent of the emissions from the central part of accretion disks, appearing as a long-term trending of the NIR continuum compared to optical variations. Moreover, the model can be further tested through reverberation mapping of emission lines, including LIGO/LISA detections of gravitational waves and signatures from spatially resolved observations of GRAVITY+/VLTI.

Pulsar wind nebulae are a possible final stage of the circumstellar evolution of massive stars, where a fast rotating, magnetised neutron star produces a powerful wind that interacts with the supernova ejecta. The shape of these so called plerionic supernova remnants is influenced by the distribution of circumstellar matter at the time of the explosion, itself impacted by the magnetic field of the ambient medium responsible for the expansion of the circumstellar bubble of the progenitor star. To understand the effects of magnetization on the circumstellar medium and resulting pulsar nebulae, we conduct 2D magnetohydrodynamical simulations. Our models explore the impact of the interstellar medium magnetic field on the morphology of a supernova remnant and pulsar wind nebula that develop in the circumstellar medium of massive star progenitor in the warm phase of the Milky Ways interstellar medium. Our simulations reveal that the jet like structures formed on both sides perpendicularly to the equatorial plane of the pulsar, creating complex radio synthetic synchrotron emissions. This morphology is characterized by a rectangular like remnant, which is typical of the circumstellar medium of massive stars in a magnetized medium, along with the appearance of a spinning top structure within the projected rectangle. We suggest that this mechanism may be partially responsible for the complex morphologies observed in pulsar wind nebulae that do not conform to the typical torus, jet or bow shock, tail shapes observed in most cases.

Methane (CH4) is a primarily biogenic greenhouse gas. As such, it represents an essential biosignature to search for life on exoplanets. Atmospheric CH4 abundance on Earth-like inhabited exoplanets is likely controlled by marine biogenic production and atmospheric photochemical consumption. Such interactions have been previously examined for the case of the early Earth where primitive marine ecosystems supplied CH4 to the atmosphere, showing that the atmospheric CH4 response to biogenic CH4 flux variations is nonlinear, a critical property when assessing CH4 reliability as a biosignature. However, the contributions of atmospheric photochemistry, metabolic reactions, or solar irradiance to this nonlinear response are not well understood. Using an atmospheric photochemical model and a marine microbial ecosystem model, we show that production of hydroxyl radicals from water vapor photodissociation is a critical factor controlling the atmospheric CH4 abundance. Consequently, atmospheric CH4 partial pressure (pCH4) on inhabited Earth-like exoplanets orbiting Sun-like stars (F-, G-, and K-type stars) would be controlled primarily by stellar irradiance. Specifically, irradiance at wavelengths of approximately 200-210 nm is a major controlling factor for atmospheric pCH4 when the carbon dioxide partial pressure is sufficiently high to absorb most stellar irradiance at 170-200 nm. Finally, we also demonstrated that inhabited exoplanets orbiting near the outer edge of K-type stars' habitable zones are better suited for atmospheric pCH4 buildup. Such properties will valuably support future detection of life signatures.

We study the problem of reconstruction of high-energy cosmic rays mass composition from the experimental data of extensive air showers. We develop several machine learning methods for the reconstruction of energy spectra of separate primary nuclei at energies 1-100 PeV, using the public data and Monte-Carlo simulations of the KASCADE experiment from the KCDC platform. We estimate the uncertainties of our methods, including the unfolding procedure, and show that the overall accuracy exceeds that of the method used in the original studies of the KASCADE experiment.

Planets orbiting young stars are thought to experience atmospheric evaporation as a result of the host stars' high magnetic activity. We study the evaporation history and expected future of the three known transiting exoplanets in the young multiplanet system K2-198. Based on spectroscopic and photometric measurements, we estimate an age of the K-dwarf host star between 200 and 500 Myr, and calculate the high-energy environment of these planets using eROSITA X-ray measurements. We find that the innermost planet K2-198c has likely lost its primordial envelope within the first few tens of Myr regardless of the age at which the star drops out of the saturated X-ray regime. For the two outer planets, a range of initial envelope mass fractions is possible, depending on the not-yet-measured planetary mass and the stars' spin-down history. Regarding the future of the system, we find that the outermost planet K2-198b is stable against photoevaporation for a wide range of planetary masses, while the middle planet K2-198d is only able to retain an atmosphere for a mass range between ~7 and 18 Earth-masses. Lower-mass planets are too susceptible to mass loss, and a very thin present-day envelope for higher-mass planets is easily lost with the estimated mass-loss rates. Our results support the idea that all three planets started out above the radius valley in the (sub-)Neptune regime and were then transformed into their current states by atmospheric evaporation, but also stress the importance of measuring planetary masses for (young) multiplanet systems before conducting more detailed photoevaporation simulations.

Grids of zero-age horizontal branch (ZAHB) models are presented, along with a suitable interpolation code, for -2.5 <= [Fe/H] <= -0.5, in steps of 0.2 dex, assuming Y = 0.25 and 0.29, [O/Fe] = +0.4 and +0.6, and [m/Fe] = 0.4 for all of the other alpha elements. The HB populations of 37 globular clusters (GCs) are fitted to these ZAHBs to derive their apparent distance moduli, (m-M)_V. With few exceptions, the best estimates of their reddenings from dust maps are adopted. The distance moduli are constrained using the prediction that (M_F606W-M_F814W)_0 colours of metal-poor, main-sequence stars at M_F606W >~ 5.0 have very little sensitivity to [Fe/H]. Intrinsic (M_F336W-M_F606W)_0 colours of blue HB stars, which provide valuable connections between GCs with exclusively blue HBs and other clusters of similar metallicity that also have red HB components, limit the uncertainties of relative (m-M)_V values to within +/- 0.03-0.04 mag. The ZAHB-based distances agree quite well with the distances derived by Baumgardt & Vasiliev (2021, MNRAS, 505, 5957). Their implications for GC ages are briefly discussed. Stellar rotation and mass loss appear to be more important than helium abundance variations in explaining the colour-magnitude diagrams of second-parameter GCs (those with anomalously very blue HBs for their metallicities).

The warm inflationary scenario is investigated in the context of affine gravity formalism. A general framework is provided for studying different single-field potentials. Using the sphaleron mechanism we explain the continuous dissipation of the inflaton field into radiation, leading to the $\Gamma=\Gamma_0 T^3$ dissipation coefficient. The treatment is performed in the weak and strong dissipation limits. We consider the quartic potential as a case study to provide a detailed study. Moreover, in this study, we discuss various constraints on inflationary models in general. We compare the theoretical results of the quartic potential model within warm inflation with the observational constraints from Planck $2018$ and BICEP/Keck 2018, as presented by the tensor-to-scalar ratio, spectral index and the perturbation spectrum.

We fit a dynamical model to Kepler systems that contain four or more transiting planets using the analytic method AnalyticLC, and obtain physical and orbital parameters for 101 planets in 23 systems, of which 95 are of mass significance better than 3 sigma, and 46 are without previously reported mass constraints nor upper limits. In addition, we compile a list of 71 KOIs that display significant Impact Parameter Variations (TbVs), complementing our previously published work on two- and three-transiting planet systems. Together, these works include the detection of significant TbV signals of 130 planets, which is, to our knowledge, the largest catalog of this type to date. The results indicate that the typical detectable TbV rate in the Kepler population is of order 10^{-2} yr^{-1}, and that rapid TbV rates (>~0.05 yr^{-1}) are observed only in systems that contain a transiting planet of an orbital period less than ~20 days. The observed TbV rates are only weakly correlated with orbital period within Kepler's <~100 days-period planets. If this extends to longer periods, it implies a limit on the utility of the transit technique for long-period planets. The TbVs we find may not be detectable in direct impact parameter measurements but rather are inferred from the full dynamics of the system, encoded in all types of transit variations. Finally, we find evidence that the mutual inclinations distribution is qualitatively consistent with the previously suggested AMD (angular momentum deficit) model using an independent approach.

Estimation of planetary orbital and physical parameters from light-curve data relies heavily on the accurate interpretation of Transit Timing Variations (TTV) measurements. In this letter, we review the process of TTV measurement and compare two fitting paradigms - one that relies on making transit-by-transit timing estimates and then fitting a TTV model to the observed timings and one that relies on fitting a global flux model to the entire light-curve data simultaneously. The latter method is achieved either by solving for the underlying planetary motion (often referred to as "photodynamics"), or by using an approximate or empirical shape of the TTV signal. We show that across a large range of the transit SNR regime, the probability distribution function (PDF) of the mid-transit time significantly deviates from a Gaussian, even if the flux errors do distribute normally. Treating the timing uncertainties as if they are distributed normally leads, in such a case, to a wrong interpretation of the TTV measurements. We illustrate these points using numerical experiments and conclude that a fitting process that relies on a global flux fitting rather than the derived TTVs, should be preferred.

We propose a measure, the joint differential entropy of eigencolours, for determining the spatial complexity of exoplanets using only spatially unresolved light curve data. The measure can be used to search for habitable planets, based on the premise of a potential association between life and exoplanet complexity. We present an analysis using disk-integrated light curves from Earth, developed in previous studies, as a proxy for exoplanet data. We show that this quantity is distinct from previous measures of exoplanet complexity due to its sensitivity to spatial information that is masked by features with large mutual information between wavelengths, such as cloud cover. The measure has a natural upper limit and appears to avoid a strong bias toward specific planetary features. This makes it a candidate for being used as a generalisable measure of exoplanet habitability, since it is agnostic regarding the form that life could take.

Asteroseismology is our only means of measuring stellar rotation in their interiors, rather than at their surfaces. Some techniques for measurements of this kind -- "rotational inversions" -- require the shapes of linear response kernels computed from reference stellar models to be representative of those in the stars they are intended to match. This is not the case in evolved stars exhibiting gravitoacoustic mixed modes: we show that the action of the asteroseismic surface term -- systematic errors in the modelling of near-surface layers -- changes the shapes of their inversion kernels. Corrections for the surface term are not ordinarily considered necessary for rotational inversions. We show how this may have caused previous estimates of red-giant envelope rotation rates from mixed-mode asteroseismic inversions to have been unintentionally contaminated by core rotation as a result, with errors comparable to the entire reported estimates. We derive a mitigation procedure for this hitherto unaccounted systematic error, and demonstrate its viability and effectiveness. We recommend this mitigation be applied when revising existing rotational inversions. Finally, we discuss both the prospects for applying such mitigation to the harder problem of inversions for stellar structure (rather than rotation), as well as the broader implications of this systematic error with regards to the longstanding problem of internal angular momentum transport.

Planetary systems around solar analogs inform us about how planets form and evolve in Solar System-like environments. We report the detection and characterization of two planetary systems around the solar analogs TOI-1736 and TOI-2141 using TESS photometry data and spectroscopic data obtained with the SOPHIE instrument on the 1.93 m telescope at the Observatoire de Haute-Provence (OHP). We performed a detailed spectroscopic analysis of these systems to obtain the precise radial velocities (RV) and physical properties of their host stars. TOI-1736 and TOI-2141 each host a transiting sub-Neptune with radii of $2.44\pm0.18$ R$_{\oplus}$ and $3.05\pm0.23$ R$_{\oplus}$, orbital periods of $7.073088(7)$ d and $18.26157(6)$ d, and masses of $12.8\pm1.8$ M$_{\oplus}$ and $24\pm4$ M$_{\oplus}$, respectively. TOI-1736 shows long-term RV variations that are consistent with a two-planet solution plus a linear trend of $-0.177$ ms$^{-1}$d$^{-1}$. We measured an RV semi-amplitude of $201.1\pm0.7$ ms$^{-1}$ for the outer companion, TOI-1736 c, implying a projected mass of $m_{c}\sin{i}=8.09\pm0.20$ M$_{\rm Jup}$. From the GAIA DR3 astrometric excess noise, we constrained the mass of TOI-1736 c at $8.7^{+1.5}_{-0.6}$ M$_{\rm Jup}$. This planet is in an orbit of $570.2\pm0.6$ d with an eccentricity of $0.362\pm0.003$ and a semi-major axis of $1.381\pm0.017$ au, where it receives a flux of $0.71\pm0.08$ times the bolometric flux incident on Earth, making it an interesting case of a supergiant planet that has settled into an eccentric orbit in the habitable zone of a solar analog. Our analysis of the mass-radius relation for the transiting sub-Neptunes shows that both TOI-1736 b and TOI-2141 b likely have an Earth-like dense rocky core and a water-rich envelope.

Tidal disruption events\,(TDEs) provide a valuable probe in studying the dynamics of stars in the nuclear environments of galaxies. Recent observations show that TDEs are strongly overrepresented in post-starburst or "green valley" galaxies, although the underlying physical mechanism remains unclear. Considering the possible interaction between stars and active galactic nucleus\,(AGN) disk, the TDE rates can be greatly changed compared to those in quiescent galactic nuclei. In this work, we revisit TDE rates by incorporating an evolving AGN disk within the framework of the "loss cone" theory. We numerically evolve the Fokker-Planck equations by considering the star-disk interactions, in-situ star formation in the unstable region of the outer AGN disk and the evolution of the accretion process for supermassive black holes\,(SMBHs). We find that the TDE rates are enhanced by about two orders of magnitude shortly after the AGN transitions into a non-active stage. During this phase, the accumulated stars are rapidly scattered into the loss cone due to the disappearance of the inner standard thin disk. Our results provide an explanation for the overrepresentation of TDEs in post-starburst galaxies.

Binary stars are believed to be key determinants in understanding globular cluster evolution. In this paper, we present the Multi-band photometric analyses of five variables in the nearest galactic globular cluster M4, from the observations of CASE, M4 Core Project with HST for four variables (V48, V49, V51, and V55) and the data collected from T40 and C18 Telescopes of Wise Observatory for one variable (NV4). The light curves of five binaries are analyzed using the Wilson-Devinney method (WD) and their fundamental parameters have been derived. A period variation study was carried out using times of minima obtained from the literature for four binaries and the nature of the variation observed is discussed. The evolutionary state of systems is evaluated using M-R diagram, correlating with that of a few well-studied close binaries in other globular clusters. Based on the data obtained from the Gaia DR3 database, a three-dimensional Vector-Point Diagrams (VPD) was built to evaluate the cluster membership of the variables, and two out of them (V49 and NV4) were found to be not cluster members.

Tidal orbital decay plays a vital role in the evolution of hot Jupiter systems. As of now, this was only observationally confirmed for the WASP-12 system. There are a few other candidates, including WASP-4 b, but no conclusive result could be obtained for these systems as of yet. In this study, we present an analysis of new TESS data of WASP-4 b together with archival data, taking the light-time effect (LTE), induced by the second planetary companion, into account as well. We make use of three different Markov-Chain-Monte-Carlo models; a circular orbit with a constant orbital period, a circular orbit with a decaying orbit, and an elliptical orbit with apsidal precession. This analysis is repeated for four cases. The first case features no LTE correction, with the remaining three cases featuring three different timing correction approaches. Comparison of these models yields no conclusive answer to the cause of WASP-4\,b's apparent transit timing variations. A broad range of values of the orbital decay and apsidal precession parameters are possible, depending on the LTE correction. This work highlights the importance of continued photometric and spectroscopic monitoring of hot Jupiters.

The cosmic string contributes to our understanding and revelation of the fundamental structure and evolutionary patterns of the universe, unifying our knowledge of the cosmos and unveiling new physical laws and phenomena. Therefore, we anticipate the detection of Stochastic Gravitational Wave Background (SGWB) signals generated by cosmic strings in space-based detectors. We have analyzed the detection capabilities of individual space-based detectors, Lisa and Taiji, as well as the joint space-based detector network, Lisa-Taiji, for SGWB signals produced by cosmic strings, taking into account other astronomical noise sources. The results indicate that the Lisa-Taiji network exhibits superior capabilities in detecting SGWB signals generated by cosmic strings and can provide strong evidence. The Lisa-Taiji network can achieve an uncertainty estimation of $\Delta G\mu/G\mu<0.5$ for cosmic string tension $G\mu\sim4\times10^{-17}$, and can provide evidence for the presence of SGWB signals generated by cosmic strings at $G\mu\sim10^{-17}$, and strong evidence at $G\mu\sim10^{-16}$. Even in the presence of only SGWB signals, it can achieve a relative uncertainty of $\Delta G\mu/G\mu<0.5$ for cosmic string tension $G\mu<10^{-18}$, and provide strong evidence at $G\mu\sim10^{-17}$.

Far-infrared (FIR) observations from the \textit{Herschel Space Observatory} are used to estimate the IR properties of UV-selected galaxies. We stack the PACS (100, 160 $\mu \mathrm{m}$) and SPIRE (250, 350 and 500$\mu \mathrm{m}$) maps of the Chandra deep field south (CDFS) on a source list of galaxies selected in the rest-frame ultraviolet (UV) in a redshift range of $0.6-1.2$. This source list is created using observations from the XMM-OM telescope survey in the CDFS using the UVW1 (2910 {\AA}) filter. The stacked data are binned according to the UV luminosity function (LF) of these sources, and the average photometry of the UV-selected galaxies is estimated. By fitting modified black bodies and IR model templates to the stacked photometry, average dust temperatures and total IR luminosity are determined. The luminosity-weighted average temperatures do not show significant evolution between the redshift bins centred at 0.7 and 1.0. The infrared excess (IRX), unobscured, and obscured SFR values are obtained from the UV and IR luminosities. Dust attenuation is constant for UV luminosities above \num{9e10} $\mathrm{L_\odot}$, but increases as UV luminosity decreases below this threshold. It remains constant as a function of IR luminosities at fixed redshift across the luminosity range of our sources. In comparison to local luminous infrared galaxies (LIRGs) with similar SFRs, the higher redshift star-forming galaxies in the sample show a lesser degree of dust attenuation. Finally, the inferred dust attenuation is used to correct the unobscured star formation rate density (SFRD) in the redshift range of 0.6 to 1.2. The dust-corrected SFRDs are found to be consistent with measurements from IR-selected samples at the same redshifts.

Observations of magnetic helicity transportation through the solar photosphere reflect the interaction of turbulent plasma movements and magnetic fields in the solar dynamo process. In this chapter, we have reviewed the research process of magnetic helicity inferred from the observed solar magnetic fields in the photosphere and also the solar morphological configurations with solar cycles. After introducing some achievements in the study of magnetic helicity, some key points would like to be summarized. The magnetic (current) helicity in the solar surface layer presents a statistical distribution similar to that of the sunspot butterfly diagram, but its maximum value is delayed from the extreme value of the sunspot butterfly diagram and corresponds in the phase with the statistical eruption of solar flares. During the spatial transport of magnetic (current) helicity from the interior of the sun into the interplanetary space at the time-space scale of the solar cycle, it shows the statistical distribution and the fluctuation with the hemispheric sign rule. These show that the current helicity and magnetic helicity transport calculation methods are complementary to each other. We also notice that the study of the inherent relationship between magnetic helicity and the solar cycle still depends on the observed accuracy of the solar magnetic field.

Only recently, complex models that include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infalls, and stellar disk in a cosmological context became available to study the dynamics of disk galaxies such as the Milky Way (MW). We use a MW model from a high-resolution hydrodynamical cosmological simulation named GARROTXA to establish the relationship between the vertical disturbances seen in its galactic disk and multiple perturbations, from the dark matter halo, satellites and gas. We calculate the bending modes in the galactic disk in the last 6 Gyr of evolution. To quantify the impact of dark matter and gas we compute the vertical acceleration exerted by these components onto the disk and compare them with the bending behavior with Fourier analysis. We find complex bending patterns at different radii and times, such as an inner retrograde mode with high frequency, as well as an outer slower retrograde mode excited at different times. The amplitudes of these bending modes are highest during the early stages of the thin disk formation and reach up to 8.5 km s-1 in the late disk evolution. We find that the infall of satellite galaxies leads to a tilt of the disk, and produces anisotropic gas accretion with subsequent star formation events, and supernovae, creating significant vertical accelerations onto the disk plane. The misalignment between the disk and the inner stellar/dark matter triaxial structure, formed during the ancient assembly of the galaxy, creates a strong vertical acceleration on the stars. We conclude that several agents trigger the bending of the stellar disk and its phase spirals in this simulation, including satellite galaxies, dark sub-halos, misaligned gaseous structures, and the inner dark matter profile, which coexist and influence each other, making it challenging to establish direct causality.

In this contribution, we briefly present the equation-of-state modelling for application to neutron stars and discuss current constraints coming from nuclear physics theory and experiments. To assess the impact of model uncertainties, we employ a nucleonic meta-modelling approach and perform a Bayesian analysis to generate posterior distributions for the equation of state with filters accounting for both our present low-density nuclear physics knowledge and high-density neutron-star physics constraints. The global structure of neutron stars thus predicted is discussed in connection with recent astrophysical observations.

Regularized imaging spectroscopy was introduced for the construction of electron flux images at different energies from count visibilities recorded by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). In this work we seek to extend this approach to data from the Spectrometer/Telescope for Imaging X-rays (STIX) on-board the Solar Orbiter mission. Our aims are to demonstrate the feasibility of regularized imaging spectroscopy as a method for analysis of STIX data, and also to show how such analysis can lead to insights into the physical processes affecting the nonthermal electrons responsible for the hard X-ray emission observed by STIX. STIX records imaging data in an intrinsically different manner from RHESSI. Rather than sweeping the angular frequency plane in a set of concentric circles (one circle per detector), STIX uses $30$ collimators, each corresponding to a specific angular frequency. In this paper we derive an appropriate modification of the previous computational approach for the analysis of the visibilities observed by STIX. This approach also allows for the observed count data to be placed into non-uniformly-spaced energy bins. We show that the regularized imaging spectroscopy approach is not only feasible for analysis of the visibilities observed by STIX, but also more reliable. Application of the regularized imaging spectroscopy technique to several well-observed flares reveals details of the variation of the electron flux spectrum throughout the flare sources. We conclude that the visibility-based regularized imaging spectroscopy approach is well-suited to analysis of STIX data. We also use STIX electron flux spectral images to track, for the first time, the behavior of the accelerated electrons during their path from the acceleration site in the solar corona toward the chromosphere

A large number of white dwarf and donor masses in cataclysmic variables have been found via modelling the primary eclipse, a method that relies on untested assumptions. Recent measurements of the mass of the white dwarf in the cataclysmic variable GY Cnc, obtained via modelling its ultraviolet spectrum, conflict with the mass obtained via modelling the eclipse light curve. Here we measure the radial velocity of the absorption lines from the donor star in GY Cnc to be $K_{\rm abs} = 280 \pm 2$ kms$^{-1}$, in excellent agreement with the prediction based on the masses derived from modelling the eclipse light curve. It is possible that the white dwarf mass derived from the ultraviolet spectrum of GY Cnc is affected by the difficulty of disentangling the white dwarf spectrum from the accretion disc spectrum.

We report on a simultaneous observational campaign with both Swift/XRT and NuSTAR targeting the symbiotic X-ray binary IGR J16194-2810. The main goal of the campaign was to investigate the possible presence of cyclotron scattering absorption features in the broad-band spectrum of the source, and help advance our understanding of the process of neutron star formation via the accretion-induced collapse of a white dwarf. The 1-30 keV spectrum of the source, as measured during our campaign, did not reveal the presence of any statistically significant absorption feature. The spectrum could be well described using a model comprising a thermal black-body hot component, most likely emerging from the surface of the accreting neutron star, and a power-law with no measurable cut-off energy (and affected by a modest absorption column density). Compared to previous analyses in the literature, we could rule out the presence of a colder thermal component emerging from an accretion disk, compatible with the idea that IGR J16194-2810 is a wind-fed binary (as most of the symbiotic X-ray binaries). Our results were strengthened by exploiting the archival XRT and INTEGRAL data, extending the validity of the spectral model used up to 0.3-40 keV and demonstrating that IGR J16194-2810 is unlikely to undergo significant spectral variability over time in the X-ray domain.

The Galactic Center Excess (GCE) $\gamma$-ray emission detected with the Large Area Telescope onboard the {\it Fermi Gamma-ray Space Telescope} has been considered as a possible sign for dark matter (DM) annihilation, but other possibilities such as the millisecond pulsar (MSP) origin have also been suggested. As a spectral fitting method, constructed based on properties of $\gamma$-ray MSPs, has been developed, we apply this method to the study of the GCE emission for the purpose of probing the MSP origin for the GCE. A number of $\sim$1660 MSPs can provide a fit to the spectrum of the GCE emission upto $\sim$10\,GeV, but the higher energy part of the spectrum requires additional emission components. We further carry out a stacking analysis of 30--500\,GeV data for relatively nearby $\gamma$-ray MSPs, and the resulting flux upper limits are still lower than those of the GCE emission. We consider the single DM annihilation channel $\tau^{+}\tau^{-}$ or channel $b\bar{b}$,or the combination of the two for comparison, and find they generally can provide better fits than MSPs. Combination of MSPs plus a DM channel are also tested, and MSPs plus the DM channel $b\bar{b}$ can always provide better fits. Comparing this combination case to the pure DM channel $b\bar{b}$, the MSP contribution is found to be marginally needed.

We present a photometric variability survey of young planetary-mass objects using the New Technology Telescope in the Js and Ks bands. Surface gravity plays an important role in the atmospheric structure of brown dwarfs, as young low gravity L dwarfs have a higher variability rate than field L dwarfs. In this study, we extend variability studies to young T-type planetary-mass objects and investigate the effects of surface gravity on the variability of L and T dwarfs across a large sample. We conduct continuous monitoring for 18 objects with spectral types from L5 to T8 and detect four new variables and two variable candidates. Combining with previous variability surveys of field and young L and T objects, we find that young objects tend to be more variable than field objects within peak-to-peak variability amplitude ranges of 0.5-10 per cent and period ranges of 1.5-20 hr. For the first time, we constrain the variability rate of young T dwarfs to be 56 per cent compared to 25 per cent for field T dwarfs. Both field and young samples have higher variability rates at the L/T transition than outside the L/T transition. The differences in the variability rates between field and young samples are about 1 sigma and therefore larger sample sizes are needed to confirm and refine the results. Besides the L/T transition, young L dwarfs with strong variability tend to assemble in a narrow spectral type range of L6-L7.5. This work supports the critical role of surface gravity on the atmospheric structure from L to T spectral types.

Models of neutron stars are considered in the case of a uniform density distribution. An algebraic equation, valid for any equation of state, is obtained. This equation allows one to find the approximate mass of a star of a given density without resorting to the integration of differential equations. The solutions presented in the paper for various equations of state, including more realistic ones, differ from the exact solutions obtained by the numerical integration of differential equations by at most ~20%

The emergence of supermassive black holes (SMBHs) in the early universe remains a topic of profound interest and debate. In this paper, we investigate the formation and growth of the first SMBHs within the framework of Modified Gravity (MOG), where gravity exhibits increased strength. We explore how MOG, as an alternative to the standard model, may offer novel insights into the emergence of SMBHs and potentially reconcile the discrepancies observed in the accretion and growth processes. We examine the dynamics of gas and matter in this modified gravitational framework, shedding light on the unique interplay between gravity and the formation of SMBHs.

In this work, we search for observational evidence of higher-order secular perturbations in three eclipsing binaries. These are slightly eccentric binaries, and they form the inner pairs of tight, compact, hierarchical triple star systems. We analyze simultaneously the high precision satellite ($Kepler$ and $TESS$) light curves, eclipse timing variations, combined spectral energy distributions (through catalog passband magnitudes) and, where available, radial velocities of KICs 9714358, 5771589 and TIC 219885468. Besides the determination of robust astrophysical and dynamical properties of the three systems, we find evidence that the observed unusual eclipse timing variations of KIC 9714358 are a direct consequence of the octupole-order secular eccentricity perturbations forced by an unusual, resonant behaviour between the lines of the apsides of the inner and outer orbital ellipses. We also show that, despite its evident cyclic eclipse depth variations, KIC~5771589 is an almost perfectly coplanar system (to within $0.3^\circ$), and we explain the rapid eclipse depth variations with the grazing nature of the eclipses. Finally, we find that the inner pair of TIC~219885468 consists of two twin stars and, hence, in this triple there are no octupole order three-body perturbations. Moreover, we show that this triple is also coplanar on the same level as the former one, but due to its deep eclipses, it does not exhibit eclipse depth variations. We intend to follow this work up with further analyses and a quantitative comparison of the theoretical and the observed perturbations.

The \ion{H}{I}-rich ultra-diffuse galaxies (HUDGs) offer a unique case for studies of star formation laws (SFLs) as they host low star formation efficiency (SFE) and low-metallicity environments where gas is predominantly atomic. We collect a sample of six HUDGs in the field and investigate their location in the extended Schmidt law($\Sigma_{\text {SFR }} \propto \left(\Sigma_{\text{star}}^{0.5} \Sigma_{\text{gas}}\right)^{1.09}$). They are consistent with this relationship well (with deviations of only 1.1 sigma). Furthermore, we find that HUDGs follow the tight correlation between the hydrostatic pressure in the galaxy mid-plane and the quantity on the x-axis ($\rm log(\Sigma_{star}^{0.5}\Sigma_{gas})$) of the extended Schmidt law. This result indicates that these HUDGs can be self-regulated systems that reach the dynamical and thermal equilibrium. In this framework, the stellar gravity compresses the disk vertically and counteracts the gas pressure in the galaxy mid-plane to regulate the star formation as suggested by some theoretical models.

In this work, we study the silicate dust content in the disk of one of the youngest eruptive stars, V900 Mon, at the highest angular resolution probing down to the inner 10 au of said disk, and study the historical evolution of the system traced in part by a newly discovered emission clump. We performed high-angular resolution mid-infrared interferometric observations of V900 Mon with MATISSE/VLTI with a spatial coverage ranging from 38-m to 130-m baselines, and compared them to archival MIDI/VLTI data. We also mined and re-analyzed archival optical and infrared photometry of the star to study its long-term evolution since its eruption in the 1990s. We complemented our findings with integral field spectroscopy data from MUSE/VLT. The MATISSE/VLTI data suggest a radial variation of the silicate feature in the dusty disk, whereby at large spatial scales ($\geq10$ au) the protostellar disk's emission is dominated by large-sized ($\geq1\,\mu m$) silicate grains, while at smaller spatial scales and closer to the star ($\leq5$ au), silicate emission is absent suggesting self-shielding. We propose that the self-shielding may be the result of small dust grains at the base of the collimated CO outflow previously detected by ALMA. A newly discovered knot in the MUSE/VLT data, located at a projected distance approximately 27,000 au from the star, is co-aligned with the molecular gas outflow at a P.A. of $250^o$ ($\pm5^o$) consistent with the position angle and inclination of the disk. The knot is seen in emission in H$\alpha$, [N II], and the [S II] doublet and its kinematic age is about 5150 years. This ejected material could originate from a previous eruption.

We present Spright, a Python package that implements a fast and lightweight mass-density-radius relation for small planets. The relation represents the joint planetary radius and bulk density probability distribution as a mean posterior predictive distribution of an analytical three-component mixture model. The analytical model, in turn, represents the probability for the planetary bulk density as three generalised Student's t-distributions with radius-dependent weights and means based on theoretical composition models. The approach is based on Bayesian inference and aims to overcome the rigidity of simple parametric mass-radius relations and the danger of overfitting of non-parametric mass-radius relations. The package includes a set of pre-trained and ready-to-use relations based on two M dwarf catalogues, one FGK star catalogue, and two theoretical composition models for water-rich planets. The inference of new models is easy and fast, and the package includes a command line tool that allows for coding-free use of the relation, including the creation of publication-quality plots. Additionally, we study whether the current mass and radius observations of small exoplanets support the presence of a population of water-rich planets positioned between rocky planets and sub-Neptunes. The study is based on Bayesian model comparison and shows somewhat strong support against the existence of a water-world population around M dwarfs. However, the results of the study depend on the chosen theoretical water-world density model. A more conclusive result requires a larger sample of precisely characterised planets and community consensus on a realistic water world interior structure and atmospheric composition model.

Dimorphos' oblate shape challenges formation models. Landslides on Didymos, induced by YORP effect, probably created a debris ring from which Dimorphos would have formed. Nonetheless, ring-derived satellites typically form with a prolate lemon shape. In light of the newest Dimorphos shape model, we revisit our previous work, Madeira et al. (2023a), and conducted new 1D simulations with material deposition in an extended region from 1.0 to 1.5 Didymos radii. An instantaneous landslide leads to a fast formation of a prolate Dimorphos directly from the ring. Now, if Didymos progressively deposits mass, Dimorphos grows through low-velocity impacts with large impactor-to-target mass ratio (pyramidal regime growth). Even during rapid, day-scale depositions, Dimorphos experiences one of these impacts, while for slower depositions, the satellite formation is primarily via pyramidal impacts. This process might reshape the satellite into an oblate shape (Leleu et al., 2018) or even in a contact-binary shape, a scenario worthy of investigation that should be studied in the future with more suitable tools. Our conclusions can be applied to Dinkinesh's satellite, recently discovered by NASA's Lucy mission.

Mergers of neutron stars (NSs) and black holes (BHs) are nowadays observed routinely thanks to gravitational-wave (GW) astronomy. In the isolated binary-evolution channel, a common-envelope (CE) phase of a red supergiant (RSG) and a compact object is crucial to sufficiently shrink the orbit and thereby enable a merger via GW emission. Here, we use the outcomes of two three-dimensional (3D) magneto-hydrodynamic CE simulations of an initially 10.0 solar-mass RSG with a 5.0 solar-mass BH and a 1.4 solar-mass NS, respectively, to explore the further evolution and final fate of the post-CE binaries. Notably, the 3D simulations reveal that the post-CE binaries are likely surrounded by circumbinary disks (CBDs), which contain substantial mass and angular momentum to influence the subsequent evolution. The binary systems in MESA modelling undergo another phase of mass transfer (MT) and we find that most donor stars do not explode in ultra-stripped supernovae (SNe), but rather in Type Ib/c SNe. The final orbits of our models with the BH companion are too wide, and NS kicks are actually required to sufficiently perturb the orbit and thus facilitate a merger via GW emission. Moreover, by exploring the influence of CBDs, we find that mass accretion from the disk widens the binary orbit, while CBD-binary resonant interactions can shrink the separation and increase the eccentricity depending on the disk mass and lifetime. Efficient resonant contractions may even enable NS/BH to merge with the remnant He stars before a second SN explosion, which may be observed as gamma-ray burst-like transients, luminous fast blue optical transients and Thorne-\.Zytkow objects. For the surviving post-CE binaries, the CBD-binary interactions may significantly increase the GW-induced double compact merger fraction. We conclude that accounting for CBD may be crucial to better understand observed GW mergers.

ANAIS is a direct dark matter detection experiment whose goal is to confirm or refute in a model independent way the positive annual modulation signal claimed by DAMA/LIBRA. Consisting of 112.5 kg of NaI(Tl) scintillators, ANAIS-112 is taking data at the Canfranc Underground Laboratory in Spain since August, 2017. Results corresponding to the analysis of three years of data are compatible with the absence of modulation and incompatible with DAMA/LIBRA. However, testing this signal relies on the knowledge of the scintillation quenching factors (QF), which measure the relative efficiency for the conversion into light of the nuclear recoil energy with respect to the same energy deposited by electrons. Previous measurements of the QF in NaI(Tl) show a large dispersion. Consequently, in order to better understand the response of the ANAIS-112 detectors to nuclear recoils, a specific neutron calibration program has been developed. This program combines two different approaches: on the one hand, QF measurements were carried out in a monoenergetic neutron beam; on the other hand, the study presented here aims at the evaluation of the QF by exposing directly the ANAIS-112 crystals to neutrons from low activity $^{252}$Cf sources, placed outside the lead shielding. Comparison between these onsite neutron measurements and detailed GEANT4 simulations will be presented, confirming that this approach allows testing different QF models.

PKS 1830$-$211 is a $\gamma$-ray emitting, high-redshift (z $= 2.507 \pm 0.002$), lensed flat-spectrum radio quasar. During the period mid-February to mid-April 2019, this source underwent a series of strong $\gamma$-ray flares that were detected by both AGILE-GRID and Fermi-LAT, reaching a maximum $\gamma$-ray flux of $F_{\rm E>100 MeV}\approx 2.3\times10^{-5}$ ph cm$^{-2}$ s$^{-1}$. Here we report on a coordinated campaign from both on-ground (Medicina, OVRO, REM, SRT) and orbiting facilities (AGILE, Fermi, INTEGRAL, NuSTAR, Swift, Chandra), with the aim of investigating the multi-wavelength properties of PKS 1830$-$211 through nearly simultaneous observations presented here for the first time. We find a possible break in the radio spectra in different epochs above 15 GHz, and a clear maximum of the 15 GHz data approximately 110 days after the $\gamma$-ray main activity periods. The spectral energy distribution shows a very pronounced Compton dominance (> 200) which challenges the canonical one-component emission model. Therefore we propose that the cooled electrons of the first component are re-accelerated to a second component by, e.g., kink or tearing instability during the $\gamma$-ray flaring periods. We also note that PKS 1830$-$211 could be a promising candidate for future observations with both Compton satellites (e.g., e-ASTROGAM) and Cherenkov arrays (CTAO) which will help, thanks to their improved sensitivity, in extending the data availability in energy bands currently uncovered.

We present a study on the coupling of haze and clouds in the atmosphere of WASP-39b. We developed a cloud microphysics model simulating the formation of Na2S and MgSiO3 condensates over photochemical hazes in gas giant atmospheres. We apply this model to WASP-39b, recently observed with the JWST to study how these heterogeneous components may affect the transit spectrum. We simulate both morning and evening terminators independently and average their transit spectra. While MgSiO3 formation has negligible impact on the spectrum, Na2S condensates produce gray opacities in the water band, in agreement with HST and JWST observations. Moreover, the formation of Na2S on the morning side depletes the atmosphere of its sodium content, decreasing the strength of the Na line. Combining morning and evening profiles results in a good fit of the Na observations. These nominal results assume a small Na2S/haze contact angle (5.7{\deg}). Using a larger value (61{\deg}) reduces the cloud density and opacity, but the effect on the Na profile and spectral line remains identical. In addition, the presence of haze in the upper atmosphere reproduces the UV-visible slope observed in the HST and VLT data and contributes to the opacity between the water bands at wavelengths below 2 microns. The averaged spectra are rather insensitive to the variation of eddy diffusion and haze mass flux tested in this study, though the UV-visible slope, probing the haze layer above the clouds, is affected. Finally, our disequilibrium chemistry model, including photochemistry, reproduces the SO2 and CO2 absorption features observed.

We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the Illustris-TNG suite. We create an HOD mock with Illustris-like correlations using our method, constraining the associated IA model parameters, with the $\chi^2_{\rm dof}$ between our model's correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation--position and orientation--orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to $0.1~h^{-1}{\rm Mpc}$. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.

Young solar-type stars frequently produce superflares, serving as a unique window into the young Sun-Earth environments. Large solar flares are closely linked to coronal mass ejections (CMEs) associated with filament/prominence eruptions, but its observational evidence for stellar superflares remains scarce. Here, we present a 12-day multi-wavelength campaign observation of young solar-type star EK Draconis (G1.5V, 50-120 Myr age) utilizing TESS, NICER, and Seimei telescope. The star has previously exhibited blueshifted H$\alpha$ absorptions as evidence for a filament eruption associated with a superflare. Our simultaneous optical and X-ray observations identified three superflares of $1.5\times10^{33}$ -- $1.2\times10^{34}$ erg. We report the first discovery of two prominence eruptions on a solar-type star, observed as blueshifted H$\alpha$ emissions at speed of 690 and 430 km s$^{-1}$ and masses of $1.1\times10^{19}$ and $3.2\times10^{17}$ g, respectively. The faster, massive event shows a candidate of post-flare X-ray dimming with the amplitude of up to $\sim$10 \%. Several observational aspects consistently point to the occurrence of a fast CME associated with this event. The comparative analysis of the estimated length scales of flare loops, prominences, possible dimming region, and starspots provides the overall picture of the eruptive phenomena. Furthermore, the energy partition of the observed superflares in the optical and X-ray bands is consistent with flares from the Sun, M-dwarfs, and close binaries, yielding the unified empirical relations. These discoveries provide profound implications of impact of these eruptive events on the early Venus, Earth and Mars and young exoplanets.

Experimental 21 cm cosmology aims to detect the formation of the first stars during the cosmic dawn and the subsequent epoch of reionization by utilizing the 21 cm hydrogen line transition. While several experiments have published results that begin to constrain the shape of this signal, a definitive detection has yet to be achieved. In this paper, we investigate the influence of uncertain antenna-sky interactions on the possibility of detecting the signal. This paper aims to define the level of accuracy to which a simulated antenna beam pattern is required to agree with the actual observing beam pattern of the antenna to allow for a confident detection of the global 21 cm signal. By utilising singular value decomposition, we construct a set of antenna power patterns that incorporate minor, physically motivated variations. We take the absolute mean averaged difference between the original beam and the perturbed beam averaged over frequency ($\Delta D$) to quantifying this difference, identifying the correlation between $\Delta D$ and antenna temperature. To analyse the impact of $\Delta D$ on making a confident detection, we utilize the REACH Bayesian analysis pipeline and compare the Bayesian evidence $\log \mathcal{Z}$ and root-mean-square error for antenna beams of different $\Delta D$ values. Our calculations suggest that achieving an agreement between the original and perturbed antenna power pattern with $\Delta D$ better than -35 dB is necessary for confident detection of the global 21 cm signal. Furthermore, we discuss potential methods to achieve the required high level of accuracy within a global 21~cm experiment.

We present a search for transient radio sources on time-scales of seconds to hours at 144 MHz using the LOFAR Two-metre Sky Survey (LoTSS). This search is conducted by examining short time-scale images derived from the LoTSS data. To allow imaging of LoTSS on short time-scales, a novel imaging and filtering strategy is introduced. This includes sky model source subtraction, no cleaning or primary beam correction, a simple source finder, fast filtering schemes and source catalogue matching. This new strategy is first tested by injecting simulated transients, with a range of flux densities and durations, into the data. We find the limiting sensitivity to be 113 and 6 mJy for 8 second and 1 hour transients respectively. The new imaging and filtering strategies are applied to 58 fields of the LoTSS survey, corresponding to LoTSS-DR1 (2% of the survey). One transient source is identified in the 8 second and 2 minute snapshot images. The source shows one minute duration flare in the 8 hour observation. Our method puts the most sensitive constraints on/estimates of the transient surface density at low frequencies at time-scales of seconds to hours; $<4.0\cdot 10^{-4} \; \text{deg}^{-2}$ at 1 hour at a sensitivity of 6.3 mJy; $5.7\cdot 10^{-7} \; \text{deg}^{-2}$ at 2 minutes at a sensitivity of 30 mJy; and $3.6\cdot 10^{-8} \; \text{deg}^{-2}$ at 8 seconds at a sensitivity of 113 mJy. In the future, we plan to apply the strategies presented in this paper to all LoTSS data.

We have witnessed different values of the Hubble constant being found in the literature in the past years. Albeit, early measurements often result in an $H_0$ much smaller than those from late-time ones, producing a statistically significant discrepancy, and giving rise to the so-called Hubble tension. The trouble with the Hubble constant is often treated as a cosmological problem. However, the Hubble constant can be a laboratory to probe cosmology and particle physics models. In our work, we will investigate if the possibility of explaining the $H_0$ trouble using non-thermal dark matter production aided by phantom-like cosmology is consistent with the Cosmic Background Radiation (CMB) and Baryon Acoustic Oscillation (BAO) data. We performed a full Monte Carlo simulation using CMB and BAO datasets keeping the cosmological parameters $\Omega_b h^2$, $\Omega_c h^2$, $100\theta$, $\tau_{opt}$, and $w$ as priors and concluded that a non-thermal dark matter production aided by phantom-like cosmology yields at most $H_0=70.5$ km s$^{-1}$Mpc$^{-1}$ which is consistent with some late-time measurements. However, if $H_0> 72$ km s$^{-1}$ Mpc$^{-1}$ as many late-time observations indicate, an alternative solution to the Hubble trouble is needed. Lastly, we limited the fraction of relativistic dark matter at the matter-radiation equality to be at most 1\%.

The arrival of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), Euclid-Wide and Roman wide area sensitive surveys will herald a new era in strong lens science in which the number of strong lenses known is expected to rise from $\mathcal{O}(10^3)$ to $\mathcal{O}(10^5)$. However, current lens-finding methods still require time-consuming follow-up visual inspection by strong-lens experts to remove false positives which is only set to increase with these surveys. In this work we demonstrate a range of methods to produce calibrated probabilities to help determine the veracity of any given lens candidate. To do this we use the classifications from citizen science and multiple neural networks for galaxies selected from the Hyper Suprime-Cam (HSC) survey. Our methodology is not restricted to particular classifier types and could be applied to any strong lens classifier which produces quantitative scores. Using these calibrated probabilities, we generate an ensemble classifier, combining citizen science and neural network lens finders. We find such an ensemble can provide improved classification over the individual classifiers. We find a false positive rate of $10^{-3}$ can be achieved with a completeness of $46\%$, compared to $34\%$ for the best individual classifier. Given the large number of galaxy-galaxy strong lenses anticipated in LSST, such improvement would still produce significant numbers of false positives, in which case using calibrated probabilities will be essential for population analysis of large populations of lenses.

The very high energy gamma-ray source HESS J1809-193 has been detected by the LHAASO and HAWC observatory beyond 100 TeV energy. It is an interesting candidate for exploring the underlying mechanisms of gamma-ray production due to the presence of supernova remnants, pulsar and molecular clouds close to it. We have considered the injection of the energetic cosmic rays from a past explosion, whose reminiscent may be SNR G011.0-00.0, which is located within the extended gamma-ray source HESS J1809-193. We explain the multi-wavelength data from the region of HESS J1809-193 with synchrotron, inverse Compton, bremsstrahlung emission of cosmic ray electrons and secondary gamma-ray production in interactions of cosmic ray protons with the cold protons in the local molecular clouds within a time-dependent framework including the diffusion loss of cosmic rays. The observational data has been modelled with the secondary photons produced by the time-evolved cosmic ray spectrum, assuming the age of the explosion is 4500 years.

Using deep JWST imaging from JADES, JEMS and SMILES, we characterize optically-faint and extremely red galaxies at $z>3$ that were previously missing from galaxy census estimates. The data indicate the existence of abundant, dusty and post-starburst-like galaxies down to $10^8$M$_\odot$, below the sensitivity limit of Spitzer and ALMA. Modeling the NIRCam and HST photometry of these red sources can result in extreme, high values for both stellar mass and star formation rate (SFR); however, including 7 MIRI filters out to 21$\mu$m results in decreased mass (median 0.6 dex for log$_{10}$M$^*$/M$_{\odot}>$10), and SFR (median 10$\times$ for SFR$>$100 M$_{\odot}$/yr). At $z>6$, our sample includes a high fraction of little red dots (LRDs; NIRCam-selected dust-reddened AGN candidates). We significantly measure older stellar populations in the LRDs out to rest-frame 3$\mu$m (the stellar bump) and rule out a dominant contribution from hot dust emission, a signature of AGN contamination to stellar population measurements. This allows us to measure their contribution to the cosmic census at $z>3$, below the typical detection limits of ALMA ($L_{\rm IR}<10^{12}L_\odot$). We find that these sources, which are overwhelmingly missed by HST and ALMA, could effectively double the obscured fraction of the star formation rate density at $4<z<6$ compared to some estimates, showing that prior to JWST, the obscured contribution from fainter sources could be underestimated. Finally, we identify five sources with evidence for Balmer breaks and high stellar masses at $5.5<z<7.7$. While spectroscopy is required to determine their nature, we discuss possible measurement systematics to explore with future data.

Aims. We aim to characterize the gas properties in the cluster outskirts ($R_{500}<r<R_{200}$) and in the detected inter-cluster filaments ($>R_{200}$) of the A3391/95 system and to compare them to predictions. Methods. We performed X-ray image and spectral analyses using the eROSITA PV data to assess the gas morphology and properties in the outskirts and the filaments in the directions of the previously detected Northern and Southern Filament of the A3391/95 system. We took particular care of the foreground. Results. In the filament-facing outskirts of A3391 and the Northern Clump, we find higher temperatures than typical cluster outskirts profiles, with a significance of $1.6-2.8\sigma$, suggesting heating due to their connections with the filaments. We confirm SB excess in the profiles of the Northern, Eastern, and Southern Filaments. From spectral analysis, we detect hot gas of ~1 keV for the Northern and Southern Filaments. The filament metallicities are below 10% solar metallicity and the $n_e$ range between 2.6 and $6.3\times10^{-5}~\mathrm{cm^{-3}}$. The characteristic properties of the Little Southern Clump (~1.5$R_{200}$ from A3395S in the Southern Filament) suggest that it is a small galaxy group. Excluding the LSC from the analysis of the Southern Filament decreases the gas density by 30%. This shows the importance of taking into account any clumps to avoid overestimation of the gas measurement in the outskirts and filament regions. Conclusions. The $n_e$ of the filaments are consistent with the WHIM properties as predicted by cosmological simulations, but the temperatures are higher, close to the upper WHIM temperature limit. As both filaments are short and located in a denser environment, stronger gravitational heating may be responsible for this temperature enhancement. The metallicities are low, but still within the expected range from the simulations.

We present a catalog of 689 galaxy cluster candidates detected at significance $\xi>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with $\sim94\%$ purity. The sample has an approximately redshift-independent mass threshold at redshift $z>0.25$ and spans $1.5 \times 10^{14} < M_{500c} < 9.1 \times 10^{14}$ $M_\odot/h_{70}$ \ and $0.03<z\lesssim1.6$ in mass and redshift, respectively; 21\% of the confirmed clusters are at $z>1$. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered $\xi$ by a median value of 0.032, or $\sim0.8\%$ of the $\xi=4$ threshold value, and $\sim7\%$ of candidates have a predicted contamination greater than $\Delta \xi = 1$. With the exception of a small number of systems $(<1\%)$, an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample will be a key component in upcoming astrophysical and cosmological analyses of clusters. The SPTpol millimeter-wave maps and associated data products used to produce this sample are available at https://pole.uchicago.edu/public/Data/Releases.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu.

We report CCD photometric observations of the globular cluster NGC 5897, in the Johnson system filters B, V , R, and I. With the values for these magnitudes we obtain various colour indices and produce several colour-magnitude diagrams. We present eight colour-magnitude diagrams: V vs B-V , B vs B-V , V vs V-I, I vs V-I, R vs R-I, I vs R-I, V vs V-R, and R vs V-R. In all of these diagrams we can clearly see the Giant Branch, the Horizontal Branch and the beginning of the Main Sequence. To the left of the Main Sequence turn-off point we detect a somewhat large number of Blue Straggler stars. We determine the mean value of the visual magnitude of the HB as $16.60 \pm 0.46$. This value is fainter than the value found by other authors.

Wind Roche-Lobe Overflow (WRLOF) is a mass-transfer mechanism for stellar binaries wherein the wind acceleration zone of the donor star exceeds its Roche lobe radius, allowing material to be transferred to the accretor. WRLOF may explain characteristics observed in blue lurkers and blue stragglers, such as their fast rotation rate. While WRLOF has been implemented in rapid population synthesis codes, it has yet to be explored thoroughly in detailed binary models such as MESA, and over a wide range of initial binary configurations. We incorporate WRLOF accretion in MESA using the POSYDON infrastructure and investigate wide low-mass binaries at solar metallicity, and perform a parameter study over initial orbital period and star mass. In most of the models where we consider angular momentum transfer during accretion, the accretor is spun up to the critical rotation rate and develops a boosted wind. Balanced by boosted wind loss, the accretor only gains $\sim 2\%$ of its total mass due to wind accretion, but can maintain a near-critical rotation rate during WRLOF. Notably, the mass transfer efficiency is significantly smaller than in previous studies in which the rotation of the accretor star is ignored. We compare our results to observational data of blue lurkers in M67 and find that the WRLOF mechanism can qualitatively explain their rapid rotation speed, their location on the HR diagram and their orbital periods.

Exoplanet characterization missions planned for the future will soon enable searches for life beyond our solar system. Critical to the search will be the development of life detection strategies that can search for biosignatures while maintaining observational efficiency. In this work, we adopted a newly developed biosignature decision tree strategy for remote characterization of Earth-like exoplanets. The decision tree offers a step-by-step roadmap for detecting exoplanet biosignatures and excluding false positives based on Earth's biosphere and its evolution over time. We followed the pathways for characterizing a modern Earth-like planet and an Archean Earth-like planet and evaluated the observational trades associated with coronagraph bandpass combinations of designs consistent with The Habitable Worlds Observatory (HWO) precursor studies. With retrieval analyses of each bandpass (or combination), we demonstrate the utility of the decision tree and evaluated the uncertainty on a suite of biosignature chemical species and habitability indicators (i.e., the gas abundances of H$_2$O, O$_2$, O$_3$, CH$_4$, and CO$_2$). Notably for modern Earth, less than an order of magnitude spread in the 1-$\sigma$ uncertainties were achieved for the abundances of H$_2$O and O$_2$, planetary surface pressure, and atmospheric temperature with three strategically placed bandpasses (two in the visible and one in the near-infrared). For the Archean, CH$_4$ and H$_2$O were detectable in the visible with a single bandpass.

We consider the effects of backreaction on axion-SU(2) dynamics during inflation. We use the linear evolution equations for the gauge field modes and compute their backreaction on the background quantities numerically using the Hartree approximation. We find a new dynamical attractor solution for the axion field and the vacuum expectation value of the gauge field, where the latter has an opposite sign with respect to the chromo-natural inflation solution. Our findings are of particular interest to the phenomenology of axion-SU(2) inflation, redefining parts of the viable parameter space. In addition, the backreaction effects lead to characteristic oscillatory features in the primordial gravitational wave background that are potentially detectable with upcoming gravitational wave detectors.

We present a novel implementation for active galactic nucleus (AGN) feedback through ultra-fast winds in the code gizmo. Our feedback recipe accounts for the angular dependence of radiative feedback upon black hole spin. We self-consistently evolve in time i) the gas accretion process from resolved scales to an unresolved AGN disc, ii) the evolution of the spin of the massive black hole (MBH), iii) the injection of AGN-driven winds into the resolved scales, and iv) the spin-induced anisotropy of the overall feedback process. We test our implementation by following the propagation of the wind-driven outflow into a homogeneous medium, and we compare the results against simple analytical models. Then, we consider an isolated galaxy setup and there we study the impact of the AGN feedback on the evolution of the MBH and the of the host galaxy. We find that: i) AGN feedback limits the gas inflow that powers the MBH, with a consequent weak impact on the host galaxy characterized by a star formation (SF) suppression of about a factor of two in the nuclear region; ii) the impact of AGN feedback on the host galaxy and on MBH growth is primarily determined by the AGN luminosity, rather than by its angular pattern set by the MBH spin; iii) the imprint of the angular pattern of the AGN radiation emission manifest in a more clear way at high accretion rates. At such high rates the more isotropic angular patterns, proper to higher spin values, sweep away gas in the nuclear region more easily, hence causing a slower MBH mass and spin growths and a higher quenching of the SF. We argue that the influence of spin-dependent anisotropy of AGN feedback on MBH and galaxy evolution is likely to be relevant in those scenarios characterized by high and prolonged MBH accretion episodes and by high AGN wind-galaxy coupling. Such conditions are more frequently met in galaxy mergers and/or high redshift galaxies.

We present the $\nu\phi$MTH, a Mirror Twin Higgs (MTH) model realizing asymmetric reheating, baryogenesis and twin-baryogenesis through the out-of-equilibrium decay of a right-handed neutrino without any hard $\mathbb{Z}_2$ breaking. The MTH is the simplest Neutral Naturalness solution to the little hierarchy problem and predicts the existence of a twin dark sector related to the Standard Model (SM) by a $\mathbb{Z}_2$ symmetry that is only softly broken by a higher twin Higgs vacuum expectation value. The asymmetric reheating cools the twin sector compared to the visible one, thus evading cosmological bounds on $\Delta N_{\mathrm{eff}}$. The addition of (twin-)colored scalars allows for the generation of the visible baryon asymmetry and, by the virtue of the $\mathbb{Z}_2$ symmetry, also results in the generation of a twin baryon asymmetry. We identify a unique scenario with top-philic couplings for the new scalars that can satisfy all cosmological, proton decay and LHC constraints; yield the observed SM baryon asymmetry; and generate a wide range of possible twin baryon DM fractions, from negligible to unity. The viable regime of the theory contains several hints as to the possible structure of the Twin Higgs UV completion. Our results motivate the search for the rich cosmological and astrophysical signatures of twin baryons, and atomic dark matter more generally, at cosmological, galactic and stellar scales.

In this dissertation, the nature of Dark Energy (DE) is examined from both theoretical and phenomenological perspectives. The possibility of DE being a dynamic quantity in quantum field theory (QFT) in curved spacetime is studied. The primary aim is to go beyond the usual approach that relies on ad hoc fields and instead treat DE as a quantum vacuum under appropriate QFT renormalization. Specifically, the dynamic behavior of DE could arise from quantum vacuum fluctuations in the Universe, evolving alongside the background expansion. Thus, the evolution of the vacuum energy density can be expressed in terms of the Hubble function and its derivatives, $\rho_{\rm vac} =\rho_{\rm vac}(H)$. This approach yields a significant revelation: the equation of state of the quantum vacuum, derived from first principles, deviates from its traditional constant value of $w_{\rm vac}=-1$. Additionally, a new inflationary mechanism emerges in this context, rooted in the quantum effects in curved spacetime. Moreover, the thesis displays a phenomenological exploration of two related models that go beyond the $\Lambda$CDM model: the Brans-Dicke model with a cosmological constant and the Running Vacuum Model, which is related to the QFT calculations. These models have been tested under different datasets and scenarios to determine the constraints on their free parameters. The results of the fits are presented and discussed in relation to cosmological tensions concerning $H_0$ and $\sigma_8$. The conclusions drawn from this thesis indicate promising signals of the dynamic behavior of quantum vacuum, potentially impacting the cosmological constant problem and the cosmological tensions.

Strong gravitational lensing can produce copies of gravitational-wave signals from the same source with the same waveform morphologies but different amplitudes and arrival times. Some of these strongly-lensed gravitational-wave signals can be demagnified and become sub-threshold. We present TESLA-X, an enhanced approach to the original GstLAL-based TargetEd Subthreshold Lensing seArch (TESLA) method, for improving the detection efficiency of these potential sub-threshold lensed signals. TESLA-X utilizes lensed injections to generate a targeted population model and a targeted template bank. We compare the performance of a full template bank search, TESLA, and TESLA-X methods via a simulation campaign, and demonstrate the performance of TESLA-X in recovering lensed injections, particularly targeting a mock event. Our results show that the TESLA-X method achieves a maximum of $\sim 20\%$ higher search sensitivity compared to the TESLA method within the sub-threshold regime, presenting a step towards detecting the first lensed gravitational wave. TESLA-X will be employed for the LIGO-Virgo-KAGRA's collaboration-wide analysis to search for lensing signatures in the fourth observing run.

We perform a dynamical system analysis and a Bayesian model selection for a new set of interacting scenarios in the framework of modified holographic Ricci dark energy models (MHR-IDE). The dynamical analysis shows a modified radiation epoch and a late-time attractor corresponding to dark energy. We use a combination of background data such as type Ia supernovae, cosmic chronometers, cosmic microwave background, and baryon acoustic oscillations measurements. We find evidence against all the MHR-IDE scenarios studied with respect to $\Lambda$CDM, when the full joint analysis is considered.

In this work, we investigate whether the compact object at the center of the Milky Way is a naked singularity described by the $\textit{q}$-metric spacetime. Our fitting of the astrometric and spectroscopic data for the S2 star implies that similarly to the Schwarzschild black hole, the $\textit{q}$-metric naked singularity offers a satisfactory fit to the observed measurements. Additionally, it is shown that the shadow produced by the naked singularity is consistent with the shadow observed by the Event Horizon Telescope collaboration for Sgr-A*. It is worth mentioning that the spatial distribution of the S-stars favors the notion that the compact object at the center of our Galaxy can be described by an almost static spacetime. Based on these findings, the $\textit{q}$-metric naked singularity turns up as a compelling candidate for further investigation.

We study the generalizations of the original Alcubierre warp drive metric to the case of curved spacetime background. We find that the presence of a horizon is essential when one moves from spherical coordinates to Cartesian coordinates in order to avoid additional singularities. For the specific case of Schwarzschild black hole, the horizon would be effectively absent for the observers inside the warp bubble, implying that warp drives may provide a safe route to cross horizons. Moreover, we discover that the black hole's gravitational field can decrease the amount of negative energy required to sustain a warp drive, which may be instrumental for creating microscopic warp drives in lab experiments. A BEC model is also introduced to propose possible test in the Analogue Gravity framework.

Small-scale interplanetary magnetic flux ropes (SMFRs) are similar to ICMEs in magnetic structure, but are smaller and do not exhibit ICME plasma signatures. We present a computationally efficient and GPU-powered version of the single-spacecraft automated SMFR detection algorithm based on the Grad-Shafranov (GS) technique. Our algorithm is capable of processing higher resolution data, eliminates selection bias caused by a fixed $\avg{B}$ threshold, has improved detection criteria demonstrated to have better results on an MHD simulation, and recovers full 2.5D cross sections using GS reconstruction. We used it to detect 512,152 SMFRs from 27 years (1996 to 2022) of 3-second cadence \emph{Wind} measurements. Our novel findings are: (1) the radial density of SMFRs at 1 au (${\sim}1$ per $\si{10^6\kilo\meter}$) and filling factor (${\sim}$35\%) are independent of solar activity, distance to the heliospheric current sheet (HCS), and solar wind plasma type, although the minority of SMFRs with diameters greater than ${\sim}$0.01 au have a strong solar activity dependence; (2) SMFR diameters follow a log-normal distribution that peaks below the resolved range ($\gtrsim 10^4$ km), although the filling factor is dominated by SMFRs between $10^5$ to $10^6$ km; (3) most SMFRs at 1 au have strong field-aligned flows like those from PSP measurements; (4) in terms of diameter $d$, SMFR poloidal flux $\propto d^{1.2}$, axial flux $\propto d^{2.0}$, average twist number $\propto d^{-0.8}$, current density $\propto d^{-0.8}$, and helicity $\propto d^{3.2}$. Implications for the origin of SMFRs and switchbacks are briefly discussed. The new algorithm and SMFR dataset are made freely available.

The radiation mechanism of fast radio bursts (FRBs) has been extensively studied but still remains elusive. In the search for dark matter candidates, the QCD axion and axionlike particles (ALPs) have emerged as prominent possibilities. These elusive particles can aggregate into dense structures called axion stars through Bose-Einstein condensation (BEC). Such axion stars could constitute a significant portion of the mysterious dark matter in the universe. When these axion stars grow beyond a critical mass, usually through processes like accretion or merging, they undergo a self-driven collapse. Traditionally, for spherically symmetric axion clumps, the interaction between axions and photons does not lead to parametric resonance, especially when the QCD axion-photon coupling is at standard levels. Nevertheless, our study indicates that even QCD axion stars with typical coupling values can trigger stimulated decay during their collapse, rather than producing relativistic axions through self-interactions. This process results in short radio bursts, with durations of around 0.1 seconds, and can be potentially observed using radio telescopes like FAST or SKA. Furthermore, we find that collapsing axion stars for ALPs with specific parameters may emit radio bursts lasting just milliseconds with a peak luminosity of $1.60\times10^{42}\rm{erg/s}$, matching the characteristics of the observed non-repeating FRBs.

In this work, we examine the propagation of gravitational waves in cosmological and astrophysical spacetimes in the context of Einstein--Gauss-Bonnet gravity, in view of the GW170817 event. The perspective we approach the problem is to obtain a theory which can produce a gravitational wave speed that is equal to that of light in the vacuum, or at least the speed can be compatible with the constraints imposed by the GW170817 event. As we show, in the context of Einstein--Gauss-Bonnet gravity, the propagation speed of gravity waves in cosmological spacetimes can be compatible with the GW170817 event, and we reconstruct some viable models. However, the propagation of gravity waves in spherically symmetric spacetimes violates the GW170817 constraints, thus it is impossible for the gravitational wave that propagates in a spherically symmetric spacetime to have a propagating speed which is equal to that of light in the vacuum. The same conclusion applies to the Einstein--Gauss-Bonnet theory with two scalars. We discuss the possible implications of our results on spherically symmetric spacetimes.

The strongly-coupled system like the quark-hadron transition (if it is of first order) is becoming an active play-yard for the physics of cosmological first-order phase transitions. However, the traditional field theoretic approach to strongly-coupled first-order phase transitions is of great challenge, driving recent efforts from holographic dual theories with explicit numerical simulations. These holographic numerical simulations have revealed an intriguing linear correlation between the phase pressure difference (pressure difference away from the wall) to the non-relativistic terminal velocity of an expanding planar wall, which has been reproduced analytically alongside both cylindrical and spherical walls from perfect-fluid hydrodynamics in our previous study but only for a bag equation of state. We have also found in our previous study a universal quadratic correlation between the wall pressure difference (pressure difference near the bubble wall) to the non-relativistic terminal wall velocity regardless of wall geometries. In this paper, we will generalize these analytic relations between the phase/wall pressure difference and terminal wall velocity into a more realistic equation of state beyond the simple bag model, providing the most general predictions so far for future tests from holographic numerical simulations of strongly-coupled first-order phase transitions

We study chirality production in the pseudoscalar inflation model of magnetogenesis taking into account the Schwinger effect and particle collisions in plasma in the relaxation time approximation. We consider the Schwinger production of one Dirac fermion species by an Abelian gauge field in two cases: (i) the fermion carries only the weak charge with respect to the U(1) group and (ii) it is also charged with respect to another strongly coupled gauge group. While the gradient-expansion formalism is employed for the description of the evolution of gauge field, plasma is described by hydrodynamical approach which allows us to determine the number, energy density, and chirality of produced fermions. It is found that while chirality production is very efficient for both, weakly and strongly interacting fermions, the resulting gauge field is typically stronger in the case of strongly interacting fermions due to suppression of the Schwinger conductivity by particle collisions.

For the analysis of gravitational-wave signals, fast and accurate gravitational-waveform models are required. These enable us to obtain information on the system properties from compact binary mergers. In this article, we introduce the NRTidalv3 model, which contains a closed-form expression that describes tidal effects, focusing on the description of binary neutron star systems. The model improves upon previous versions by employing a larger set of numerical-relativity data for its calibration, by including high-mass ratio systems covering also a wider range of equations of state. It also takes into account dynamical tidal effects and the known post-Newtonian mass-ratio dependence of individual calibration parameters. We implemented the model in the publicly available LALSuite software library by augmenting different binary black hole waveform models (IMRPhenomD, IMRPhenomX, and SEOBNRv5_ROM). We test the validity of NRTidalv3 by comparing it with numerical-relativity waveforms, as well as other tidal models. Finally, we perform parameter estimation for GW170817 and GW190425 with the new tidal approximant and find overall consistent results with respect to previous studies.