Papers for Friday, Mar 31 2023

David Kuridze, Lyndsay Fletcher and Hugh Hudson report on the RAS Specialist Discussion Meeting '3D Structure of the Flare Chromosphere'.

The observations of Ic supernovae (Ic/SNe) occurring after the prompt emission of long gamma-ray bursts (GRBs) are addressed within the binary-driven hypernova (BdHN) model. The GRBs originate from a binary composed of a $\sim 10~M_\odot$ carbon-oxygen (CO) star and a companion neutron star (NS). We assume these same progenitors originate the Ic/SN. The binary evolution depends strongly on the binary period, $P_{\rm bin}$. The trigger, given by the CO core collapse, for $P_{\rm bin}$ of up to a few hours leads to an Ic/SN with a fast-spinning NS ($\nu$NS) at its center. For $P_{\rm bin} \sim 4$--$5$ min, BdHN I occur with energies $10^{52}$--$10^{54}$ erg, a contribution by the black hole (BH) created by the NS companion collapse, originates the Mev/GeV radiations. The $\sim$~1 millisecond $\nu$NS originates, by synchrotron radiation, the X-ray afterglow. For $P_{\rm bin} \sim 10$~min, BdHN II occurs with energies of $10^{50}$--$10^{52}$~erg. For $P_{\rm bin} \sim$ hours, BdHN III occurs with energies below $10^{50}$~erg. The $1$--$1000$ ms $\nu$NS, in all BdHNe, originates the X-ray afterglow by synchrotron emission. The SN Ic follows an independent evolution, becoming observable by the nickel decay after the GRB prompt emission. We report $24$ Ic/SNe associated with BdHNe; their optical peak luminosity and their time of occurrence are similar and independent of the associated GRBs. We give four examples of BdHNe and their associated hypernovae. For the first time, we approach new physical processes in BdHNe, identifying seven episodes and their signatures in their spectra.

Uncertainties in stellar population models, both in terms of stellar evolution and stellar spectra, translate into uncertainties in our interpretation of stellar populations in galaxies, since stars are the source of most of the light we receive from them. Observations by JWST are revealing high-redshift galaxies in great detail, which must then be compared to models. One significant source of uncertainty is in the stellar spectra used to generate composite spectra of stellar populations, which are then compared to data. Confidence in theoretical models is important to enable reliable determination of the properties of these galaxies such as their ages and star formation history. Here we present a comparison of spectral synthesis carried out with 6 different stellar spectral libraries using the Binary Population and Spectral Synthesis (BPASS) framework. In photometric colours, the differences between theoretical libraries are relatively small (<0.10 mag), similar to typical observational uncertainties on individual galaxy observations. Differences become more pronounced when detailed spectroscopic properties are examined. Predictions for spectral line indices can vary significantly, with equivalent widths differing by a factor of two in some cases. With these index strengths, some of the libraries yield predictions of ages and metallicities which are unphysical. Many spectral libraries lack wavelength coverage in the ultraviolet, which is of growing importance in the era of JWST observations of distant galaxies, whose flux is dominated by hot, young stars.

We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in F150W and F480M, respectively. At a distance of 10$\,$pc based on its Spitzer parallax, and assuming a random inclination distribution, the physical separation is approximately 1$\,$au. Evolutionary models predict for that an age of 1-5 Gyr, the companion mass is about 4-12.5 Jupiter masses around the 7.5-20 Jupiter mass primary, corresponding to a companion-to-host mass fraction of $q=0.61\pm0.05$. Under the assumption of a Keplerian orbit the period for this extreme binary is in the range of 5-9 years. The system joins a small but growing sample of ultracool dwarf binaries with effective temperatures of a few hundreds of Kelvin. Brown dwarf binaries lie at the nexus of importance for understanding the formation mechanisms of these elusive objects, as they allow us to investigate whether the companions formed as stars or as planets in a disk around the primary.

We report on our study of SN 2022xxf during the first four months of its evolution. The light curves (LCs) display two humps at similar maximum brightness separated by 75d, unprecedented for a broad-lined Type Ic supernova (SN IcBL). SN~2022xxf is the most nearby SN IcBL to date (in NGC~3705, $z = 0.0037$, 20 Mpc). Optical and NIR photometry and spectroscopy are used to identify the energy source powering the LC. Nearly 50 epochs of high S/N-ratio spectroscopy were obtained within 130d, comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution of SN~2022xxf points to typical SN Ic/IcBL, with broad features (up to $\sim14000$ km~s$^{-1}$) and a gradual transition from the photospheric to the nebular phase. However, narrow emission lines (corresponding to $\sim1000-2500$ km~s$^{-1}$) are present from the time of the second rise, suggesting slower-moving circumstellar material (CSM). These lines are subtle, but some are readily noticeable at late times such as in Mg~I $\lambda$5170 and [O~I] $\lambda$5577. Unusually, the near-infrared spectra show narrow line peaks, especially among features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and subsequent rapid LC drop. SN~2022xxf may be related to an emerging number of CSM-interacting SNe Ic, which show slow, peculiar LCs, blue colors, and subtle CSM interaction lines. The progenitor stars of these SNe likely experienced an episode of mass loss shortly prior to explosion consisting of H/He-free material.

Clusters of galaxies are the largest gravitationally-bound systems in the Universe. Their dynamics are dominated by dark matter (DM), which makes them among the best targets for indirect DM searches. We analyze 12 years of data collected by the Fermi Large Area Telescope (Fermi-LAT) in the direction of 49 clusters of galaxies selected for their proximity to the Earth and their high X-ray flux, which makes them the most promising targets. We first create physically motivated models for the DM density around each cluster considering different assumptions for the substructure distribution. Then we perform a combined search for a $\gamma$-ray signal in the {\it Fermi}-LAT data between 500 MeV and 1 TeV. We find a signal of $\gamma$ rays potentially associated with DM that is at a statistical significance of $2.5\sigma-3.0\sigma$ when considering a slope for the subhalo mass distribution $\alpha=1.9$ and minimum mass of $M_{\rm{min}}=10^{-6}$ $M_{\odot}$. The best-fit DM mass and annihilation cross-sections for a $b\bar{b}$ annihilation channel are $m_{\chi}=40-60$ GeV and $\langle \sigma v \rangle = (2-4) \times 10^{-25}$ cm$^3$/s. When we consider $\alpha=2.0$ and $M_{\rm{min}}=10^{-9}$ $M_{\odot}$, the best-fit of the cross section reduces to $\langle \sigma v \rangle = (4-10) \times 10^{-26}$ cm$^3$/s. For both DM substructure models there is a tension between the values of $\langle \sigma v \rangle$ that we find and the upper limits obtained with the non-detection of a $\gamma$-ray flux from Milky Way dwarf spheroidal galaxies. This signal is thus more likely associated with $\gamma$ rays produced in the intracluster region by cosmic rays colliding with gas and photon fields.

We present a detailed study of cool, neutral gas traced by Lya around 4595 z<0.5 galaxies using stacks of background quasar spectra. The galaxies are selected from our MUSEQuBES low-z survey along with data from the literature. These galaxies, with a median stellar mass of log (M*/Msun)= 10.0, are probed by 184 background quasars giving rise to 5054 quasar-galaxy pairs. The median impact parameter is b = 1.5 pMpc (median b/Rvir=10.4) with 204 (419) quasar-galaxy pairs probing b/Rvir < 1 (2). We find excess absorption out to at least ~ 15 Rvir transverse distance and ~ 600 km/s along the line of sight. We show that the median stacked profile for the full sample, dominated by the pairs with b > Rvir, can be explained by a galaxy-absorber two-point correlation function with r0 = 6.7 pMpc and gamma = -1.57. There are strong indications that the inner regions (< Rvir) of the rest equivalent width profile are better explained by a log-linear (or a Gaussian) relation whereas the outer regions are well described by a power-law, consistent with galaxy-absorber large-scale clustering. Using a sub-sample of 339 galaxies (442 quasar-galaxy pairs, median b/Rvir = 1.6) with star formation rate measurements, we find that the Lya absorption is significantly stronger for star-forming galaxies compared to passive galaxies, but only within the virial radius. The Lya absorption at b ~ Rvir for a redshift-controlled sample peaks at M* ~ 10^9 Msun~ (Mhalo ~ 10^11 Msun).

In the second paper of this series, we developed a new distance determination method using the median $J$ magnitude of carbon-rich asymptotic giant branch stars (CS) as standard candles and the Magellanic Clouds as the fundamental calibrators. The $J$-band CS luminosity function was modeled using a modified Lorentzian distribution whose parameters were used to determined whether the LMC or SMC was the most suitable calibrator. In this third paper of the series, we expand our sample of galaxies and introduce a more robust method to determine the parameters of the Lorentzian model. The new fitting method uses an un-binned maximum likelihood estimator to determine the parameters of the Lorentzian model resulting in parameter errors that are significantly smaller compared to the second paper. We test our method in NGC 6822, IC 1613, NGC 3109 and WLM. We also estimate the distances to the same sample of galaxies via the tip of the red giant branch (TRGB) detection method. Our results from the CS measurements agree well with those obtained from the TRGB.

Yellow hypergiants (YHGs) are often presumed to represent a transitional post-red supergiant (RSG) phase for stars $\sim$30-40 \msun. Here we present visual-wavelength echelle spectra of six YHG candidates in the Galactic cluster Westerlund 1, and we compare them to known YHGs, IRC +10420 and Hen3-1979. We find that the six YHG candidates do not exhibit any metallic emission lines, nor do they show strong H$\alpha$ emission, and as such do not meet the criteria necessary to be classified as YHGs. In conjunction with their moderate luminosities of \logl = 4.7-5.4 estimated from optical/infrared photometry, we suggest instead that they are normal yellow supergiants (YSGs) with more modest initial masses around 15-20 \msun. This adds additional support to the hypothesis that Wd1 is a multi-age cluster with an older age than previously assumed, and is not a $\sim$5 Myr old cluster caught at a very specific transitional point when single-star evolution might yield Wolf-Rayet stars, luminous blue variables (LBVs), RSGs, and YHGs in the same cluster. Nevertheless, the population of YSGs in Wd1 is very unusual, with YSGs outnumbering RSGs, but with both spanning a large luminosity range. Here, we discuss evolutionary scenarios that might have led to the high fraction of YSGs. The number of YSGs and their significant luminosity spread cannot be explained by simple population synthesis models with single or binary stars. Even with multiple ages or a large age spread, the high YSG/RSG ratio remains problematic. We suggest instead that the objects may experience a prolonged YSG phase due to evolution in triple systems.

The enigmatic open clusters serve as a constant reminder of the mysteries of the universe, helping to confront astronomical theories. Unknown to many, these clusters often possess tails with inappropriate labels, serving as the tell-tale signs of their historical journey. But unlike typical tails, these extensions can either precede or follow the body, yet they consistently unfold a cosmic mystery to be solved. I present a succinct survey of this subject matter, detailing the intrepid efforts of astronomers who have dared to challenge our knowledge about these creatures, and offer a novel proposal for their nomenclature, while not disregarding the philosophical ramifications.

The gamma-ray burst (GRB) 221009A, with its extreme brightness, has provided the opportunity to explore GRB prompt and afterglow emission behavior on short time scales with high statistics. In conjunction with detection up to very high-energy gamma-rays, studies of this event shed light on the emission processes at work in the initial phases of GRBs emission. Using INTEGRAL/IBIS's soft gamma-ray detector, PICsIT (200-2600 keV), we studied the temporal and spectral evolution during the prompt phase and the early afterglow period. We found a "flux-tracking" behavior with the source spectrum "softer" when brighter. However the relationship between the spectral index and the flux changes during the burst. The PICsIT light curve shows afterglow emission begins to dominate at ~ T0 + 630s and decays with a slope of 1.6 +/- 0.2, consistent with the slopes reported at soft X-rays.

Core-collapse supernovae are predicted to produce gravitational waves (GWs) that may be detectable by Advanced LIGO/Virgo. These GW signals carry information from the heart of these catacylsmic events, where matter reaches nuclear densities. Recent studies have shown that it may be possible to infer properties of the proto-neutron star (PNS) via gravitational waves generated by hydrodynamic perturbations of the PNS. However, we lack a comprehensive understanding of how these relationships may change with the properties of core-collapse supernovae. In this work, we build a self-consistent suite of over 1000 exploding core-collapse supernovae from a grid of progenitor masses and metallicities combined with six different nuclear equations of state. Performing a linear perturbation analysis on each model, we compute the resonant gravitational-wave frequencies of the PNS, and we motivate a time-agnostic method for identifying characteristic frequencies of the dominant gravitational-wave emission. From this, we identify two characteristic frequencies, of the early- and late-time signal, that measure the surface gravity of the cold remnant neutron star, and simultaneously constrain the hot nuclear equation of state. However, we find that the details of the core-collapse supernova model, such as the treatment of gravity or the neutrino transport, and whether it explodes, noticeably change the magnitude and evolution of the PNS eigenfrequencies.

We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiting $K$'-band magnitude of 15.8, which is at least 4 magnitudes deeper than previous spectroscopic searches for binaries at larger radii within the central nuclear star cluster. From this primary dataset, over 600 new radial velocities are extracted and reported, increasing by a factor of 3 the number of such measurements. We find no significant periodic signals in our sample of 28 stars, of which 16 are massive, young (main-sequence B) stars and 12 are low-mass, old (M and K giant) stars. Using Monte Carlo simulations, we derive upper limits on the intrinsic binary star fraction for the young star population at 47% (at 95% confidence) located $\sim$20 mpc from the black hole. The young star binary fraction is significantly lower than that observed in the field (70%). This result is consistent with a scenario in which the central supermassive black hole drives nearby stellar binaries to merge or be disrupted and may have important implications for the production of gravitational waves and hypervelocity stars.

We present the first study of galaxy evolution in $\ddot{\mu}$ based cosmologies. We find that recent JWST observations of massive galaxies at extremely high redshifts are consistent with such a cosmology. However, the low redshift Universe is entirely divergent from the $\ddot{\mu}$ cosmic star formation rate density. We thus propose that our Universe was at one point dominated by a Primordial Bovine Herd (PBH) which later decayed producing dark energy. Note that we do not detail the mechanisms by which this decay process takes place. Despite its vanishingly small probability for existence, a $\ddot{\mu}$ based cosmological model marries the disparate findings in the high and low redshift Universe.

What's in a name, a poet once asked. To which we present this work, where we investigate the importance of a paper title in ensuring its best outcome. We queried astronomy papers using NASA ADS and ranked 6000 of them in terms of cheekiness level. We investigate the correlation between citation counts and (i) the presence of a colon, and (ii) cheekiness ranking. We conclude that colon matters in the anatomy of a paper title. So does trying to be cheeky, but we find that too much cheekiness can lead to cringefests. Striking the right balance is therefore crucial. May we recommend aiming for a level 4 cheekiness on a scale of 1-5.

Ueta & Otsuka (2021) proposed a method, named as the "Proper Plasma Analysis Practice", to analyze spectroscopic data of ionized nebulae. The method is based on a coherent and simultaneous determination of the reddening correction and physical conditions in the nebulae. The same authors (Ueta & Otsuka 2022, UO22) reanalyzed the results of Galera-Rosillo et al. (2022, GR22) on nine of the brightest planetary nebulae in M31. They claim that, if standard values of the physical conditions are used to compute the extinction instead of their proposed method, extinction correction is underestimated by more than 50% and hence, ionic and elemental abundance determinations, especially the N/O ratio, are incorrect. Several tests were performed to assess the accuracy of the results of GR22, when determining: i) the extinction coefficient, ii) the electron temperature and density, and iii) the ionic abundances. In the latter case, N+ /H+ ionic abundance was recalculated using both H_alpha and H_beta as the reference H I emissivity. The analysis shows that the errors introduced by adopting standard values of the plasma conditions by GR22 are small, within their quoted uncertainties. On the other hand, the interstellar extinction in UO22 is found to be overestimated for five of the nine nebulae considered. This propagates into their analysis of the properties of the nebulae and their progenitors. The python notebook used to generate all the results presented in this paper are of public access on a Github repository. The results from GR22 are proven valid and the conclusions of the paper hold firmly. Although the PPAP is, in principle, a recommended practice, we insist that it is equally important to critically assess which H I lines are to be included in the determination of the interstellar extinction coefficient, and to assert that physical results are obtained for the undereddened line ratios.

A link between exposures of water (H${}_{2}$O) ice with traces of an ammoniated compound (e.g., a salt) and the probable effusion of a water-rich cryolava onto the surface of Pluto has been established in previous investigations (Dalle Ore et al. 2019). Here we present the results from the application of a machine learning technique and a radiative transfer model to a water-ice-rich exposure in Kiladze area and surroundings on Pluto. We demonstrate the presence of an ammoniated material suggestive of an undetermined but relatively recent emplacement event. Kiladze lies in a region of Pluto's surface that is structurally distinct from that of the areas where similar evidence points to cryovolcanic activity at some undetermined time in the planet's history. Although the Kiladze depression superficially resembles an impact crater, a close inspection of higher-resolution images indicates that the feature lacks the typical morphology of a crater. Here we suggest that a cryolava water carrying an ammoniated component may have come onto the surface at the Kiladze area via one or more volcanic collapses, as in a resurgent volcanic caldera complex. Large regions east of Kiladze also exhibit the presence of H${}_{2}$O ice and have graben-like structures suggestive of cryovolcanic activity, but with existing data are not amenable to the detailed search that might reveal an ammoniated component.

We explore the early evolution of flame ignition and spreading on the surface of a neutron star in three-dimensions, in the context of X-ray bursts. We look at the nucleosynthesis and morphology of the burning front and compare to two-dimensional axisymmetric simulations to gauge how important a full three-dimensional treatment of the flame is for the early dynamics. Finally, we discuss the progress toward full-star resolved flame simulations.

We measure and analyze galaxy clustering and the dependence on luminosity, color, age, stellar mass and specific star formation rate using Baryon Oscillation Spectroscopic Survey (BOSS) galaxies at $0.48<z<0.62$. We fit the monopole and quadrupole moments of the two-point correlation function (2PCF) and its projection on scales of $0.1$ -- $60.2h^{-1}$Mpc, after having split the catalog in a variety of ways. We find that the clustering dependence is consistent with previous well-established results showing the broad trends expected: For example, that brighter, redder, older, more massive and quenched galaxies are more strongly clustered. We also investigate the dependence on additional parameters previously derived from stellar population synthesis model fits to the spectra. We find that galaxy clustering depends on look-back formation time at a low level, while it has little dependence on metallicity. To understand the physics behind these trends, we fit the clustering with a simulation-based emulator to simultaneously model cosmology and galaxy bias using a Halo Occupation Distribution framework. After marginalizing parameters determining the background cosmology, galaxy bias, and a scaling parameter to decouple halo velocity field, we find that the growth rate of large scale structure as determined by the redshift-space distortions is consistent with previous analysis using the full sample and is independent of the galaxy selection. This demonstrates that cosmological inference using small scale clustering measurements is robust to changes in the catalog selection.

Planetesimals -- asteroids and comets -- are the building blocks of planets in protoplanetary discs and the source of dust, ice and gas in debris discs. Along with planets they comprise the left-over material after star formation that constitutes a planetary system. Planets influence the dynamics of planetesimals, sculpting the orbits of debris belts to produce asymmetries or gaps. We can constrain the architecture of planetary systems, and infer the presence of unseen planetary companions, by high spatial resolution imaging of debris discs. HD~16743 is a relatively young F-type star that hosts a bright edge-on debris disc. Based on far-infrared \textit{Herschel} observations its disc was thought to be stirred by a planetary companion. Here we present the first spatially resolved observations at near-infrared and millimetre wavelengths with \textit{HST} and ALMA, revealing the disc to be highly inclined at $87\fdg3~^{+1\fdg9}_{-2\fdg5}$ with a radial extent of 157.7$^{+2.6}_{-1.5}$~au and a FWHM of 79.4$^{+8.1}_{-7.8}$~au ($\Delta R/R = 0.5$). The vertical scale height of the disc is $0.13~\pm~0.02$, significantly greater than typically assumed unstirred value of 0.05, and could be indicative of stirring of the dust-producing planetesimals within the disc by bodies at least a few times the mass of Pluto up to 18.3~$M_{\oplus}$ in the single object limit.

Fast radio bursts (FRBs) are astronomical transients with millisecond timescales. Although most of the FRBs are not observed to repeat, a few of them are detected to repeat more than hundreds of times. There exist a large variety of physical properties among these bursts, suggesting heterogeneous mechanisms of FRBs. In this paper, we conduct a categorisation on the extremely frequently repeating FRB 20201124A with the assistance of machine learning, as such techniques have the potential to use subtle differences and correlations that humans are unaware of to better classify bursts. The research is carried out by applying the unsupervised Uniform Manifold Approximation and Projection (UMAP) model on the FRB 20201124A data provided by Five-hundred-meter Aperture Spherical radio Telescope (FAST). The algorithm eventually categorises the bursts into three clusters. In addition to the two categories in previous work based on waiting time, a new way for categorisation has been found. The three clusters are either high energy, high frequency, or low frequency, reflecting the distribution of FRB energy and frequency. Importantly, a similar machine learning result is found in another frequently repeating FRB20121102A, implying a common mechanism among this kind of FRB. This work is one of the first steps towards the systematical categorisation of the extremely frequently repeating FRBs.

We analysed state-of-the-art observations of the solar atmosphere to investigate the dependence of the \ca brightness of several solar features on spectral bandwidth and spatial resolution of the data. Specifically, we study data obtained at the Swedish Solar Telescope with the CRiSP and CHROMIS instruments. The analyzed data, which are characterized by spectral bandwidth of 0.12 \AA \ and spatial resolution of 0.078\arcsec, were acquired close to disc center by targeting a quiet Sun area and an active region. We convolved the original observations with Gaussian kernels to degrade their spectral bandwidth and spatial resolution to the instrumental characteristics of the most prominent series of \ca observations available to date. We then studied the effect of data degradation on the observed regions and on parameters derived from \ca line measurements that are largely employed as diagnostics of the solar and stellar chromospheres. We find that the effect of degrading the spectral resolution of \ca observations and line profiles depends on both the employed bandwidth and observed solar region. Besides, we found that the spatial degradation impacts the data characterized by a broad bandwidth to a larger extent compared to those acquired with a narrow band. However, the appearance of the observed solar regions is only slightly affected by the spatial resolution of data with bandwidths up to 1 \AA \ and in the range [3,10] \AA. Finally, we derived relationships that can be used to intercalibrate results from observations taken with different instruments in diverse regions of the solar atmosphere.

Debris discs around main sequence stars have been extensively characterised from infrared to millimetre wavelengths through imaging, spectroscopic, and total intensity (scattered light and/or thermal emission) measurements. Polarimetric observations have only been used sparingly to interpret the composition, structure, and size of dust grains in these discs. Here we present new multi-wavelength aperture polarisation observations with parts-per-million sensitivity of a sample of twelve bright debris discs, spanning a broad range of host star spectral types, and disc properties. These measurements were mostly taken with the HIgh Precision Polarimetric Instrument on the Anglo-Australian Telescope. We combine these polarisation observations with the known disc architectures and geometries of the discs to interpret the measurements. We detect significant polarisation attributable to circumstellar dust from HD 377 and HD 39060, and find tentative evidence for HD 188228 and HD 202628.

Tests of gravity are important to the development of our understanding of gravitation and spacetime. Binary pulsars provide a superb playground for testing gravity theories. In this chapter we pedagogically review the basics behind pulsar observations and pulsar timing. We illustrate various recent strong-field tests of the general relativity (GR) from the Hulse-Taylor pulsar PSR~B1913+16, the double pulsar PSR~J0737$-$3039, and the triple pulsar PSR~J0337+1715. We also overview the inner structure of neutron stars (NSs) that may influence some gravity tests, and have used the scalar-tensor gravity and massive gravity theories as examples to demonstrate the usefulness of pulsar timing in constraining specific modified gravity theories. Outlooks to new radio telescopes for pulsar timing and synergies with other strong-field gravity tests are also presented.

Different approaches have been adopted to study short- and long-term stellar magnetic activity, and although the mechanisms by which low-mass stars generate large-scale magnetic fields are not well understood, it is known that stellar rotation plays a key role. There are stars that show a cyclical behaviour in their activity which can be explained by solar dynamo or $\alpha\Omega$ dynamo models. However, when studying late-type dwarf stars, it is necessary to implement other indicators to analyse their magnetic activity. In the present work, we perform a comparative study between the best-known activity indicators so far defined from the Ca II and H$\alpha$ lines to analyse M-dwarf stars. We studied a sample of 29 M stars with different chromospheric activity levels and spectral classes ranging from dM0 to dM6. To do so, we employed 1796 wide range spectra from different instruments with a median time span of observations of 21 yr. In addition, we complemented our data with photometric observations from the TESS space mission for better stellar characterisation and short-term analysis. We obtained a good and significant correlation ($rho = 0.91$) between the indexes defined from the two lines for the whole set of stars in the sample. However, we found that there is a deviation for faster rotators (with $P_{rot} < 4$ days) and higher flare activity (at least one flare per day). There is an overall positive correlation between Ca II and H$\alpha$ emission in dM stars, except during flare events. In particular, we found that low-energy high-frequency flares could be responsible for the deviation in the linear trend in fast-rotator M dwarfs. This implies that the rotation period could be a fundamental parameter to study the stellar activity and that the rotation could drive the magnetic dynamo in low-mass active stars.

We utilise the unprecedented depth and resolution of recent early-release science (ERS) JWST observations to define the near-infrared counterparts of sub-millimetre selected galaxies (SMGs). We identify 45 SCUBA-2 SMG positions within The Cosmic Evolution Early Release Science Survey (CEERS) JWST/NIRCam fields. Through an analysis of multi-wavelength $p$-values, NIRCam colours and predicted SCUBA-2 fluxes, we define 43 JWST/NIRCam counterparts to the SCUBA-2 SMGs, finding a 63 per cent agreement with those identified in prior $HST$ studies. Using EaZy-py we fit the available HST and JWST observations to quantify the photometric redshifts of the NIRCam-SMGs, establishing a broad range of redshift from $z$$\approx$0.2$-$5.4 with a median of $z$$\approx$2.29, in agreement with other studies of SMGs. We analyse their rest-frame optical and near-infrared morphological properties (e.g. effective radius (R$_{\rm e}$), S\'ersic index ($n$), CAS, Gini and M$_{20}$), finding, on average, late-type disc-like morphologies with large scatter into the intermediate and merger regions of the non-parametric parameter space. For the non-merging galaxies, we find a median rest-frame optical size and S\'ersic index (and $1\sigma$ scatter) of R$_{\rm e}$=3.10$\pm$1.67kpc and $n$=0.96$\pm$0.66. Whilst in the rest-frame near-infrared we establish more compact, higher S\'ersic index morphologies (R$_{\rm e}$=1.64$\pm$0.97, $n$=1.85$\pm$0.63). We further establish that both the rest-frame optical and near-infrared effective radii correlate negatively (at a 2$\sigma$ level) with redshift whilst the S\'ersic index remains constant with cosmic time. Our results are consistent with the picture of inside-out galaxy evolution, with more centrally concentrated older stellar populations, and more extended, younger star-forming regions whose stellar emission is heavily attenuated in the central regions.

Exomoons have so far eluded ongoing searches. Several studies have exploited transit and transit timing variations and high-resolution spectroscopy to identify potential exomoon candidates. One method of detecting and confirming these exomoons is to search for signals of planet-moon interactions. In this work, we present the first radio observations of the exomoon candidate system WASP 69b. Based on the detection of alkali metals in the transmission spectra of WASP-69b, it was deduced that the system might be hosting an exomoon. WASP 69b is also one of the exoplanet systems that will be observed as part of JWST cycle-1 GTO. This makes the system an excellent target to observe and follow up. We observed the system for 32 hrs at 150 MHz and 218 MHz using the upgraded Giant Metrewave Radio Telescope (uGMRT). Though we do not detect radio emission from the systems, we place strong $3\sigma$ upper limits of 3.3 mJy at 150 MHz and 0.9 mJy at 218 MHz. We then use these upper limits to estimate the maximum mass loss from the exomoon candidate.

We are carrying out the GPPS survey by using the FAST, the most sensitive systematic pulsar survey in the Galactic plane. In addition to about 500 pulsars already discovered through normal periodical search, we report here the discovery of 76 new transient radio sources with sporadic strong pulses, detected by using the newly developed module for a sensitive single pulse search. Their small DM values suggest that they all are the Galactic RRATs. More radio pulses have been detected from 26 transient radio sources but no periods can be found due to a limited small number of pulses from all FAST observations. The following-up observations show that 16 transient sources are newly identified as being the prototypes of RRATs with a period already determined from more detected sporadic pulses, 10 sources are extremely nulling pulsars, and 24 sources are weak pulsars with sparse strong pulses. On the other hand, 48 previously known RRATs have been detected by the FAST. Except for 1 RRAT with four pulses detected in a session of five minute observation and 4 RRATs with only one pulse detected in a session, sensitive FAST observations reveal that 43 RRATs are just generally weak pulsars with sporadic strong pulses or simply very nulling pulsars, so that the previously known RRATs always have an extreme emission state together with a normal hardly detectable weak emission state. This is echoed by the two normal pulsars J1938+2213 and J1946+1449 with occasional brightening pulses. Though strong pulses of RRATs are very outstanding in the energy distribution, their polarization angle variations follow the polarization angle curve of the averaged normal pulse profile, suggesting that the predominant sparse pulses of RRATs are emitted in the same region with the same geometry as normal weak pulsars.

The work describes a system for converting VLBI observation data using the algorithms of coherent dedispersion and compensation of two-bit signal sampling. Coherent dedispersion is important for processing pulsar observations to obtain the best temporal resolution, while correction for signal sampling makes it possible to get rid of a number of parasitic effects that interfere with the analysis of the diffraction pattern of pulsars. A pipeline has been established that uses the developed converter and the ASC Software Correlator, which will allow reprocessing all archived data of Radioastron pulsar observations and to conduct a search for giant pulses, which requires the best temporal resolution.

Magnetic field is ubiquitous in interstellar media and channels turbulent flow from kilo-parsec to sub-parsec scales in both diffuse ISM and molecular clouds. The determination of magnetohydrodynamic (MHD) turbulence and 3D magnetic field properties in ISM is notoriously difficult. In this study, we establish a statistical recipe `Y-parameter' based on the recent development of turbulence statistical theory, namely the turbulence anisotropy analysis, to reconstruct 3D magnetic fields in MHD simulations and understand the mode-decomposition of MHD turbulence. In our analysis, we used 25 MHD turbulence datacubes simulated using ZEUS-MP and Athena++ codes. We found that the anisotropy of Stokes parameters can act as a diagnostic for retrieving the magnetic field inclination in ISM and identifying dominating mode. It is supported by the value space separation of the Y-parameter for decomposed Alfvenic and compressible MHD cubes and which is decreasing and increasing with mean field inclination angle, $\theta_{\lambda}$, respectively. In the total cube analysis, Y$\sim 1.5$ (with Y$>1.5$ for A-mode and Y$<1.5$ for C-mode) provides a statistical demarcation to obtain the dominant fraction of MHD turbulence modes. Furthermore, we have found that the (i) if Y$\gtrsim 2.5$, $10^\circ<\theta_{\lambda}< 30^\circ$ and A mode or $5^\circ<\theta_{\lambda}< 10^\circ$ and C-mode (ii) if Y$ \lesssim 1.0 $, $\theta_{\lambda} \lesssim 5^\circ$ and C-mode or $\theta_{\lambda} \gtrsim 60^\circ$ and A-mode (iii) $40^\circ \lesssim \theta_{\lambda} \lesssim 60^\circ$ with A-mode or $\theta_{\lambda} \gtrsim 70^\circ$ with C-mode, if Y-parameter is in the intermediate range. As a consequence, in the future with the availability of vast radio polarisation surveys, this technique can play a leading role in detecting 3D magnetic field in ISM and characterizing the nature of the interstellar turbulence.

The ARCADE radio excess and EDGES measurement remain puzzling. A link between the two has been previously considered, however, in this work we highlight an important related effect that was not analyzed in detail before. By performing cosmological thermalization calculations with soft photon injection using {\tt CosmoTherm}, we show that for the 21 cm signal generation the interplay between enhanced radio spectral distortions and the associated heating can hide a significant radio excess before the reionzation era. We illustrate this effect for a simple power-law soft photon source in decaying particle scenarios. Even if simplistic, the uncovered link between CMB spectral distortions and 21 cm cosmology should apply to a much broader range of scenarios. This could significantly affect the constraints derived from existing and future 21 cm observations on the evolution of the ambient radio background. In particular, scenarios that would be ruled out by existing data without heating could become viable solutions once the heating is accounted for in the modelling. Our calculations furthermore highlight the importance of global 21 cm observations reaching into the dark ages, where various scenarios can potentially be distinguished.

We study the effect of plasma composition on the dynamics and morphology of the relativistic astrophysical jets. Our work is based on a relativistic total variation diminishing (TVD) simulation code. We use a relativistic equation of state in the simulation code which accounts for the thermodynamics of a multispecies plasma which is a mixture of electrons, positrons, and protons. To study the effect of plasma composition we consider various jet models. These models are characterized by the same injection parameters, same jet kinetic luminosity, and the same Mach numbers. The evolution of these models shows that the plasma composition affects the jet head propagation speed, the structure of the jet head, and the morphology despite fixing the initial parameters. We conclude that the electron-proton jets are the slowest and show more pronounced turbulent structures in comparison to other plasma compositions. The area and locations of the hot-spots also depend on the composition of jet plasma. Our results also show that boosting mechanisms are also an important aspect of multi-dimensional simulations which are also influenced by the change in composition.

X Persei is a persistent low-luminosity X-ray pulsar of period of $\sim$835 s in a Be binary system. The field strength at the neutron star surface is not known precisely, but indirect signs indicate a magnetic field above $10^{13}$ G, which makes the object one of the most magnetized known X-ray pulsars. Here we present the results of observations X Persei performed with the Imaging X-ray Polarimetry Explorer (IXPE). The X-ray polarization signal was found to be strongly dependent on the spin phase of the pulsar. The energy-averaged polarization degree in 3-8 keV band varied from several to $\sim$20 per cent over the pulse with a positive correlation with the pulsed X-ray flux. The polarization angle shows significant variation and makes two complete revolutions during the pulse period resulting in nearly nil pulse-phase averaged polarization. Applying the rotating vector model to the IXPE data we obtain the estimates for the rotation axis inclination and its position angle on the sky as well as for the magnetic obliquity. The derived inclination is close to the orbital inclination reported earlier for X Persei. The polarimetric data imply a large angle between the rotation and magnetic dipole axes, which is similar to the result reported recently for the X-ray pulsar GRO J1008$-$57. After eliminating the effect of polarization angle rotation over the pulsar phase using the best-fitting rotating vector model, the strong dependence of the polarization degree with energy was discovered with its value increasing from 0% at $\sim$2 keV to 30% at 8 keV.

We investigate the exponential $f(Q)$ symmetric teleparallel gravitation, namely $f(Q)=Q+\alpha Q_0(1-e^{-\beta\sqrt{Q/Q_0}})$ using \texttt{ME-GADGET} code to probe the structure formation with box sizes $L_{\mathrm{box}}=10/100$ Mpc$/h$ and middle resolution $N_p^{1/3}=512$. To reproduce viable cosmology within the aforementioned modified gravity theory, we first perform Markov Chain Monte Carlo (MCMC) sampling on OHD/BAO/Pantheon datasets and constrain a parameter space. Furthermore, we also derive theoretical values for deceleration parameter $q(z)$, statefinder pair $\{r,s\}$ and effective gravitational constant $G_{\mathrm{eff}}$, perform $Om(z)$ diagnostics. While carrying out N-body+SPH simulations, we derive CDM+baryons over density/temperature/mean molecular weight fields, matter power spectrum (both 2/3D, with/without redshift space distortions), bispectrum, two-point correlation function and halo mass function. Results for small and big simulation box sizes are therefore properly compared, halo mass function is related to the Seth-Tormen theoretical prediction and matter power spectrum to the standard \texttt{CAMB} output.

Primordial Black Holes (PBHs) may form in the early Universe, from the gravitational collapse of large density perturbations, generated by large quantum fluctuations during inflation. Since PBHs form from rare over-densities, their abundance is sensitive to the tail of the primordial probability distribution function (PDF) of the perturbations. It is therefore important to calculate the full PDF of the perturbations, which can be done non-perturbatively using the 'stochastic inflation' framework. In single field inflation models generating large enough perturbations to produce an interesting abundance of PBHs requires violation of slow roll. It is therefore necessary to extend the stochastic inflation formalism beyond slow roll. A crucial ingredient for this are the stochastic noise matrix elements of the inflaton potential. We carry out analytical and numerical calculations of these matrix elements for a potential with a feature which violates slow roll and produces large, potentially PBH generating, perturbations. We find that the transition to an ultra slow-roll phase results in the momentum induced noise terms becoming larger than the field noise whilst each of them falls exponentially for a few e-folds. The noise terms then start rising with their original order restored, before approaching constant values which depend on the nature of the slow roll parameters in the post transition epoch. This will significantly impact the quantum diffusion of the coarse-grained inflaton field, and hence the PDF of the perturbations and the PBH mass fraction.

Production of gravitationally coupled light moduli fields must be suppressed in the early universe, so that its decay products do not alter Big Bang Nucleosynthesis (BBN) predictions for light elements. On the other hand, the moduli quanta can be copiously produced non-thermally during preheating after the end of inflation. In this work, we study the production of moduli in the $\alpha$-attractor inflationary model through parametric resonances. For our case, where the inflationary potential at its minimum is quartic, the inflaton field self-resonates, and subsequently induces large production of moduli particles. We find that this production is suppressed for small values of $\alpha$. Combining semi-analytical estimation and numerical lattice simulations, we infer the parametric dependence on $\alpha$ and learn that $\alpha$ needs to be $\lesssim 10^{-8}\,m_{\rm Pl}^2$ to be consistent with BBN. This in turn predicts an upper bound on the energy scale of inflation and on the reheating temperature. Further, it implies an extremely small tensor-to-scalar ratio that quantifies the amplitude of primordial gravitational waves over large scales.

Context. Warm corona is a possible explanation for Soft X-ray Excess in Active Galactic Nuclei (AGN). This paper contains self consistent modeling of both: accretion disk with optically thick corona, where the gas is heated by magneto-rotational instability dynamo (MRI), and cooled by radiation which undergoes free-free absorption and Compton scattering. Aims. We determine the parameters of warm corona in AGN using disk-corona structure model that takes into account magnetic and radiation pressure. We aim to show the role of thermal instability (TI) as a constraint for warm, optically thick X-ray corona in AGN. Methods. With the use of relaxation code, the vertical solution of the disk driven by MRI together with radiative transfer in hydrostatic and radiative equilibrium is calculated, which allows us to point out how TI affects the corona for wide range of global parameters. Results. We show that magnetic heating is strong enough to heat upper layers of the accretion disk atmosphere, which form the warm corona covering the disk. Magnetic pressure does not remove TI caused by radiative processes operating in X-ray emitting plasma. TI disappears only in case of accretion rates higher than 0.2 of Eddington, and high magnetic field parameter $\alpha_{\rm B}$ > 0.1. Conclusions. TI plays the major role in the formation of the warm corona above magnetically driven accretion disk in AGN. The warm, Compton cooled corona, responsible for soft X-ray excess, resulted from our model has typical temperature in the range of 0.01 - 2 keV and optical depth even up to 50, which agrees with recent observations.

I employ an optimization-based inference methodology together with an Ising model, in an intentionally ineffectual manner, to get away with murdering an obstreperous scientific collaborator. The antics of this collaborator, hereafter "Conan O'Brien," were impeding the publication of an important manuscript. With my tenure date looming, I found myself desperate. Luckily, I study inference, a computational means to find a solution to a physical problem, based on available measurements (say, a dead body) and a dynamical model assumed to give rise to those measurements (a murderer). If the measurements are insufficient and/or the model is incomplete, one obtains multiple "degenerate" solutions to the problem. Degenerate solutions are all equally valid given the information available, and thus render meaningless the notion of one "correct" solution. Typically in scientific research, degeneracy is undesirable. Here I describe the opposite situation: a quest to create degenerate solutions in which to cloak myself. Or even better: to render measurements incompatible with a solution in which I am the murderer. Moreover, I show how one may sabotage an inference procedure to commit an untraceable crime. I sit here now, typing victoriously, a free woman. Because you won't believe me anyway. And even if you do, you'll never prove a thing.

In this work, we adopt a cosmological model-independent approach for the first time to test the question of whether the mass density power-law index($\gamma$) of the strong gravitational lensing system(SGLS) evolves with redshift, and the JLA SNe Ia sample and the quasar sample from Risaliti \& Lusso (2019) are used to provide the luminosity distances to be calibrated. Our work is based on the flat universe assumption and the cosmic distance duality relation. A reliable data-matching method is used to pair SGLS-SNe and SGLS-quasar. By using the maximum likelihood method to constrain the luminosity distance and $\gamma$ index, we obtain the likelihood function values for the evolved and non-evolved cases, and then use the Akaike weights and the BIC selection weights to compare the advantages and disadvantages of these two cases. We find that the $\gamma$ index is slightly more likely to be a non-evolutionary model for $\gamma=2$ in the case of the currently used samples with low redshift ($z_l<\sim$0.66). With Akaike weights, the relative probability is 66.3\% versus 33.7\% and 69.9\% versus 30.1\% for the SGLS+SNe Ia sample and SGLS+quasar sample, respectively, and with BIC selection weights, the relative probability is 87.4\% versus 12.6\% and 52.0\% versus 48.0\% for the two samples. In the evolving case for the relatively low redshift lens (SGLS+SNe Ia), with redshift 0.0625 to 0.659, $\gamma= 2.058^{+0.041}_{-0.040}-0.136^{+0.163}_{-0.165}z$. At high redshift (SGLS+quasar ), with redshift 0.0625 to 1.004, $\gamma= 2.051^{+0.076}_{-0.077}-0.171^{+0.214}_{-0.196}z$. Although not the more likely model, this evolved $\gamma$ case also fits the data well, with a negative and mild evolution for both low and high redshift samples.

White-light superflares from ultra cool stars are thought to be resulted from magnetic reconnection, but the magnetic dynamics in a fully convective star is not clear yet. In this paper, we report a stellar superflare detected with the Ground Wide Angle Camera (GWAC), along with rapid follow-ups with the F60A, Xinglong 2.16m and LCOGT telescopes. The effective temperature of the counterpart is estimated to be $2200\pm50$K by the BT-Settl model, corresponding to a spectral type of L0. The $R-$band light curve can be modeled as a sum of three exponential decay components, where the impulsive component contributes a fraction of 23\% of the total energy, while the gradual and the shallower decay phases emit 42\% and 35\% of the total energy, respectively. The strong and variable Balmer narrow emission lines indicate the large amplitude flare is resulted from magnetic activity. The bolometric energy released is about $6.4\times10^{33}$ ergs, equivalent to an energy release in a duration of 143.7 hours at its quiescent level. The amplitude of $\Delta R=-8.6 $mag ( or $\Delta V=-11.2$ mag), placing it one of the highest amplitudes of any ultra cool star recorded with excellent temporal resolution. We argue that a stellar flare with such rapidly decaying and huge amplitude at distances greater than 1 kpc may be false positive in searching for counterparts of catastrophic events such as gravitational wave events or gamma-ray bursts, which are valuable in time-domain astronomy and should be given more attention.

Low-mass stars and sub-stellar objects are essential in tracing the initial mass function (IMF). We study the nearby young $\sigma$ Orionis cluster (d$\sim$408 pc; age$\sim$1.8 Myr) using deep NIR photometric data in J, W and H-bands from WIRCam on the Canada-France-Hawaii Telescope. We use the water absorption feature to photometrically select the brown dwarfs and confirm their nature spectroscopically with the IRTF-SpeX. Additionally we select candidate low-mass stars for spectroscopy and analyze their membership and that of literature sources using astrometry from Gaia DR3. We obtain the near-IR spectra for 28 very low-mass stars and brown dwarfs and estimate their spectral type between M3-M8.5 (mass ranging between 0.3-0.01 M$_{\odot}$). Apart from these, we also identify 5 new planetary mass candidates which require further spectroscopic confirmation of youth. We compile the comprehensive catalog of 170 spectroscopically confirmed members in the central region of the cluster, for a wide mass range of $\sim$19-0.004 M$_{\odot}$. We estimate the star/BD ratio to be $\sim$4, within the range reported for other nearby star forming regions. With the updated catalog of members we trace the IMF down to 4 M$_\mathrm{Jup}$ and we find that a two-segment power-law fits the sub-stellar IMF better than the log-normal distribution.

A significant fraction of debris discs consist of a bright ring beyond which extends a wide halo. Such a halo should be made of small grains produced in the ring of parent bodies (PB) and pushed on high-e orbits by radiation pressure. It has been shown that, under several simplifying assumptions, the surface brightness (SB) of this halo should radially decrease as $r^{-3.5}$ in scattered light. We aim to revisit the halo phenomenon and focus on two so far unexplored issues: 1) How the unavoidable presence of small unbound grains, non-isotropic scattering phase functions (SPF) and finite instrument resolution affect scattered light SB profiles, and 2) How the halo phenomenon manifests itself at longer wavelengths. We find that unbound grains account for a significant fraction of the halo's luminosity in scattered light, and can significantly flatten the SB radial profile. Realistic size-dependent SPFs also have an effect, resulting here again in shallower SB profiles. For edge-on discs, non-resolving the vertical profile can also flatten the projected SB. We show that roughly half of the observationally-derived halo profiles found in the literature are compatible with our new results, and that roughly half of the remaining systems are probably shaped by additional processes. We also propose that, in future observational studies, the characteristics of PB belt and halos should be fitted separately. In thermal emission, wide halos should remain detectable up to the far-IR and, with the exception of the $\sim 8-15\mu$m domain, the halo accounts for more than half of the system's total flux up to $\lambda\sim80-90\mu$m. The halo's contribution strongly decreases in the sub-mm to mm but still represents a few percents of the system's luminosity at $\lambda\sim 1$mm. For unresolved systems, the presence of a halo can also affect the determination of the disc's radius from its SED.

We study the non-linear dynamics of axion inflation, capturing for the first time the inhomogeneity and full dynamical range, till the end of inflation. Accounting for inhomogeneous effects during backreaction leads to a number of new relevant results, compared to spatially homogeneous studies: {\it i)} the number of extra efoldings beyond slow roll inflation increases very rapidly with the coupling, {\it ii)} oscillations of the inflaton velocity are attenuated, {\it iii)} the tachyonic gauge field helicity spectrum is smoothed out (i.e.~the spectral oscillatory features disappear), broadened, and shifted to smaller scales, and {\it iv)} the non-tachyonic helicity is excited, reducing the chiral asymmetry, now scale dependent. Our results are expected to impact strongly on the phenomenology and observability of axion inflation, including gravitational wave generation and primordial black hole production.

Stream-stream collision may be an important pre-peak energy dissipation mechanism in tidal disruption events (TDEs). We perform local three-dimensional radiation hydrodynamic simulations in a wedge geometry including the gravity to study stream self-crossing, with emphasis on resolving the collision and following the subsequent outflow. We find that the collision can contribute to pre-peak optical emissions by converting $\gtrsim5\%$ of stream kinetic energy to radiation, yielding prompt emission of $\sim10^{42-44}\rm erg~s^{-1}$. The radiative efficiency is sensitive to stream mass fallback rates, and strongly depends on the downstream gas optical depth. Even for a sub-Eddington ($10\%$) mass fallback rate, the strong radiation pressure produced in the collision can form a local super-Eddington region near the collision site, where a fast, aspherical outflow is launched. Higher mass fallback rate usually leads to more optically-thick outflow and lower net radiative efficiency. For $\dot{M}\gtrsim0.1\dot{M}_{\rm Edd}$, the estimated photosphere size of the outflow can expand by one to two orders of magnitudes reaching $\sim10^{14}\rm cm$. The average gas temperature at this photospheric surface is a few $\times10^{4}$K, roughly consistent with inferred pre-peak photosphere properties for some optical TDEs. We find that the dynamics is sensitive to collision angle and collision radius, but the radiative efficiency or outflow properties show more complex dependency than is often assumed in ballistic models.

The dust properties of the line-of-sight materials in neutron star low-mass X-ray binaries (LMXBs) can be probed by X-ray observations and laboratory experiments. We use a Markov chain Monte Carlo (MCMC) method to conduct a spectral analysis of Chandra ACIS-S/HETG archival data of a sample of LMXBs, including GX 5-1 and GX 13+1. Our MCMC-based analysis puts constraints on the Si K-edge dust properties of the outflowing disk winds in this sample. Further X-ray observations of other LMXBs will help us better understand the grain features of dense outflows and accretion flows in neutron star binary systems.

Stellar metallicity is a critical factor to characterize the stellar coronae because it directly affects the radiative energy loss from the atmosphere. By extending theoretical relations for solar coronal loops introduced by \cite{Rosner1978}, we analytically derive scaling relations for stellar coronal loops with various metallicities. In order to validate the derived relations, we also perform magnetohydrodyamic simulations for the heating of coronal loops with different metallicities by changing radiative loss functions according to the adopted elemental abundances. The simulation results nicely explain the generalized analytical scaling relations and show a strong dependence of the thermodynamical and radiative properties of the loops on metallicity. Higher density and temperature are obtained in lower-metallicity coronae because of the inefficient radiative cooling, provided that the surface condition is unchanged. Thus, it is estimated that the X-ray radiation from metal-poor coronae is higher because of their denser coronal gas. The generalized scaling laws can also be used as a tool to study the condition of high-energy radiation around magnetically active stars and their impact on planetary environments.

NGC 7582 (z = 0.005264; D = 22.5 Mpc) is a highly variable, changing-look AGN. In this work, we explore the X-ray properties of this source using XMM-Newton and NuSTAR archival observations in the 3-40 keV range, from 2001 to 2016. NGC 7582 exhibits a long-term variability between observations but also a short-term variability in two observations that has not been studied before. To study the variability, we perform a time-resolved spectral analysis using a phenomenological model and a physically-motivated model (uxclumpy). The spectral fitting is achieved using a nested sampling Monte Carlo method. uxclumpy enables testing various geometries of the absorber that may fit AGN spectra. We find that the best model is composed of a fully covering clumpy absorber. From this geometry, we estimate the velocity, size and distance of the clumps. The column density of the absorber in the line of sight varies from Compton-thin to Compton-thick between observations. Variability over the timescale of a few tens of kilo-seconds is also observed within two observations. The obscuring clouds are consistent with being located at a distance not larger than 0.6 pc, moving with a transverse velocity exceeding $\sim 700$ km s$^{-1}$. We could put only a lower limit on the size of the obscuring cloud being larger than $10^{13}$ cm. Given the sparsity of the observations, and the limited exposure time per observation available, we cannot determine the exact structure of the obscuring clouds. The results are broadly consistent with comet-like obscuring clouds or spherical clouds with a non-uniform density profile.

The growth of galaxies in the early Universe is driven by accretion of circum- and inter-galactic gas. Simulations predict that steady streams of cold gas penetrate the dark matter halos of galaxies, providing the raw material necessary to sustain star formation. We report a filamentary stream of gas that extends for 100 kiloparsecs and connects to the massive radio galaxy 4C 41.17. The stream is detected using sub-millimeter observations of the [CI] line of atomic carbon, a tracer of neutral atomic or molecular hydrogen gas. The galaxy contains a central gas reservoir that is fueling a vigorous starburst. Our results show that the raw material for star formation can be present in cosmic streams outside galaxies.

While observations of molecular gas at cosmic noon and beyond have focused on the gas within galaxies (i.e., the interstellar medium; ISM), it is also crucial to study the molecular gas reservoirs surrounding each galaxy (i.e., in the circumgalactic medium; CGM). Recent observations of galaxies and quasars hosts at high redshift (z>2) have revealed evidence for cold gaseous halos of scale r_CGM~10kpc, with one discovery of a molecular halo with r_CGM~200kpc and a molecular gas mass one order of magnitude larger than the ISM of the central galaxy. As a follow-up, we present deep ACA and ALMA observations of CO(3-2) from this source and two other quasar host galaxies at z~2.2. While we find evidence for CO emission on scales of r~10kpc, we do not find evidence for molecular gas on scales larger than r>20 kpc. Therefore, our deep data do not confirm the existence of massive molecular halos on scales of ~100 kpc for these X-ray selected quasars. As an interesting by-product of our deep observations, we obtain the tentative detection of a negative continuum signal on scales larger than r>200kpc, which might be tracing the Sunyaev-Zeldovich effect associated with the halo heated by the active galactic nucleus (AGN). If confirmed with deeper data, this could be direct evidence of the preventive AGN feedback process expected by cosmological simulations.

Measurements of neutron star mass and radius or tidal deformability deliver unique insight into the equation of state (EOS) of cold dense matter. EOS inference is very often done using generalized parametric or non-parametric models which deliver no information on composition. In this paper we consider a microscopic nuclear EOS model based on a field theoretical approach. We show that current measurements from NICER and gravitational wave observations constrain primarily the symmetric nuclear matter EOS. We then explore what could be delivered by measurements of mass and radius at the level anticipated for future large-area X-ray timing telescopes. These should be able to place very strong limits on the symmetric nuclear matter EOS, in addition to constraining the nuclear symmetry energy that determines the proton fraction inside the neutron star.

Ultralight dark matter (ULDM) is proposed as a theoretical candidate of dark matter particles with masses of approximately $10^{-22}$ eV. The interactions between ULDM particles and standard model particles would cause variations in pulse arrival times of the pulsars, which makes the pulsar timing array (PTA) can be used to indirectly detect ULDM. In this work, we use the gamma-ray PTA composed of 29 millisecond pulsars observed by the Fermi Large Area Telescope (Fermi-LAT) to test four ultralight dark matter effects, including gravitational effects for generalized ULDM with different Spin-0/1, the fifth-force coupling effect of Dark Photon and the modified gravitational effect of the Spin-2 ULDM. The gamma-ray pulsar timing is not affected by the ionized interstellar medium and suffers relatively simple noises, unlike that of the radio band. No significant signals of ULDM are found based on the Fermi-LAT PTA for all four kinds of ULDM models. Constraints on ULDM parameters are set with the 95% confidence level, which provides a complementary check of the non-detection of ULDM for radio PTAs and direct detection experiments.

We have found evidence of bulk velocities following active galactic nucleus (AGN) bubbles in the Virgo cluster and galaxy motions in the Centaurus cluster. In order to increase the sample and improve our understanding of the intracluster medium (ICM), we present the results of a detailed mapping of the Ophiuchus cluster with {\it XMM-Newton} to measure bulk flows through very accurate Fe~K measurements. To measure the gas velocities we use a novel EPIC-pn energy scale calibration, which uses the Cu K$\alpha$ instrumental line as reference for the line emission. We created 2D spectral maps for the velocity, metallicity, temperature, density, entropy and pressure with a spatial resolution of 0.25$'$ ($\sim 26$~kpc). The ICM velocities in the central regions where AGN feedback is most important are similar to the velocity of the brightest cluster galaxy (BCG). We have found a large interface region where the velocity changes abruptly from blueshifted to redshifted gas which follows a sharp surface brightness discontinuity. We also found that the metallicities and temperatures do not change as we move outwards from the giant radio fossil previously identified in radio observations of the cluster. Finally, we have found a contribution from the kinetic component of $<25\%$to the total energy budget for large distances.

We study the formation of the Stokes profiles of the He I multiplet at 10830 A when relaxing two of the approximations that are often considered in the modeling of this multiplet, namely the lack of self-consistent radiation transfer and the assumption of equal illumination of the individual multiplet components. This He I multiplet is among the most important ones for the diagnostic of the outer solar atmosphere from spectropolarimetric observations, especially in prominences, filaments, and spicules. However, the goodness of these approximations is yet to be assessed, especially in situations where the optical thickness is of the order or larger than one, and radiation transfer has a significant impact in the local anisotropy and the ensuing spectral line polarization. This issue becomes particularly relevant in the ongoing development of new inversion tools which take into account multi-dimensional radiation transfer effects. To relax these approximations we generalize the multi-term equations for the atomic statistical equilibrium to allow for differential illumination of the multiplet components and implement them in a one-dimensional radiative transfer code. We find that, even for this simple geometry and relatively small optical thickness, both radiation transfer and differential illumination effects have a significant impact on the emerging polarization profiles. This should be taken into account in order to avoid potentially significant errors in the inference of the magnetic field vector.

We study gravitational wave microlensing by primordial black holes (PBHs), accounting for the effect of a particle dark matter minihalo surrounding them. Such minihaloes are expected when PBHs make up only a fraction of all dark matter. We find that the LIGO-Virgo detections imply a $1\sigma$ bound on the abundance of PBHs heavier than $50 M_{\odot}$. The next generation observatories can potentially probe PBHs as light as $0.01 M_\odot$ and down to $2\times10^{-4}$ fraction of all dark matter. We also show that these detectors can distinguish between dressed and naked PBHs, providing a novel way to study the distribution of particle dark matter around black holes and potentially shed light on the origins of black holes.

Recent numerical studies have revealed the physically intriguing fact that charged black holes whose charge-to-mass ratios are larger than the critical value $(Q/M)_{\text{crit}}=\sqrt{2(9+\sqrt{6})}/5$ can support hairy matter configurations which are made of scalar fields with a non-minimal negative coupling to the Gauss-Bonnet invariant of the curved spacetime. Using {\it analytical} techniques, we explore the physical and mathematical properties of the composed charged-black-hole-nonminimally-coupled-linearized-massless-scalar-field configurations in the near-critical $Q/M\gtrsim (Q/M)_{\text{crit}}$ regime. In particular, we derive an analytical resonance formula that describes the charge-dependence of the dimensionless coupling parameter $\bar\eta_{\text{crit}}=\bar\eta_{\text{crit}}(Q/M)$ of the composed Einstein-Maxwell-nonminimally-coupled-scalar-field system along the {\it existence-line} of the theory, a critical border that separates bald Reissner-Nordstr\"om black holes from hairy charged-black-hole-scalar-field configurations. In addition, it is explicitly shown that the large-coupling $-\bar\eta_{\text{crit}}(Q/M)\gg1$ analytical results derived in the present paper for the composed Einstein-Maxwell-scalar theory agree remarkably well with direct numerical computations of the corresponding black-hole-field resonance spectrum.

This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided.

For much of February 2023, the world was in panic as repeated balloon-like unidentified flying objects (UFOs) were reported over numerous countries by governments that often responded with military action. As a result, most of these craft either escaped or were destroyed, making any further observation of them nearly impossible. These were not the first time balloon-like objects have loomed over Earth, nor are they likely to be the last. This has prompted us to push for a better understanding of UFOs. First we demonstrate that the distribution of balloon incidents and other UFO reports are consistent with being drawn from the same geographic distribution, and further that both of these distributions are consistent with the areas of the Earth that feature the jet stream. Second we show that there are more UFO sightings during meteor showers, as we would expect if meteor showers, already a known source of extraterrestrial material, are being used to provide some manner of distraction to help alien craft enter the Earth's atmosphere without drawing undue attention. These links between alleged balloon incidents, UFO reports, and meteor showers establish a transport pipeline for alien craft from interplanetary and possibly interstellar space to the Earth's surface.

The range of the U bosonic coupling constants in neutron star matter is a very interesting but still unsolved problem which has multifaceted influences in nuclear physics, particle physics, astrophysics and cosmology. The combination of the theoretical numerical simulation and the recent observations provides a very good opportunity to solve this problem. In the present work, the range of the U bosonic coupling constants is inferred based on the three relations of the mass-radius, mass-frequency and mass-tidal deformability in neutron star containing hyperons using the GM1, TM1 and NL3 parameter sets under the two flavor symmetries of the SU(6) and SU(3) in the framework of the relativistic mean field theory. Combined with observations from PSRs J1614-2230, J0348+0432, J2215-5135, J0952-0607, J0740+6620, J0030-0451, J1748-2446ad, XTE J1739-285, GW170817 and GW190814 events, our numerical results show that the U bosonic coupling constants may tend to be within the range from 0 to 20 GeV$^{-2}$ in neutron star containing hyperons. Moreover, the numerical results of the three relations obtained by the SU(3) symmetry are better in accordance with observation data than those obtained by the SU(6) symmetry. The results will help us to improve the strict constraints of the equation of state for neutron stars containing hyperons.

Understanding the equation of state of dense QCD matter remains a major challenge in both nuclear physics and astrophysics. Neutron star observations from electromagnetic and gravitational wave spectra provide critical insights into the behavior of dense neutron-rich matter. The next generation of telescopes and gravitational wave observatories will offer even more detailed observations of neutron stars. Utilizing deep learning techniques to map neutron star mass and radius observations to the equation of state allows for its accurate and reliable determination. This work demonstrates the feasibility of using deep learning to extract the equation of state directly from neutron star observational data, and to also obtain related nuclear matter properties such as the slope, curvature, and skewness of the nuclear symmetry energy at saturation density. Most importantly, we show that this deep learning approach is able to reconstruct \textit{realistic} equations of state, and deduce \textit{realistic} nuclear matter properties. This highlights the potential of artificial neural networks in providing a reliable and efficient means to extract crucial information about the equation of state and related properties of dense neutron-rich matter in the era of multi-messenger astrophysics.

We demonstrate that it is possible to test models of gravity, such as Palatini $f(R)$ and Eddington-inspired Born-Infeld models, using seismic data from Earth. By incorporating additional limitations on Earth's moment of inertia and mass given from observational data, the models' parameters can be restricted to a $2\sigma$ level of accuracy. Our novel tool provides that $\beta\lesssim 10^9 \text{m}^2$ for Palatini and $\epsilon\lesssim 4\cdot 10^9 \text{m}^2$ for Eddington-inspired Born-Infeld gravity. We also discuss further enhancements to the proposed method, aimed at imposing even more stringent constraints on modified gravity proposals.

In this paper, we would like to examine whether a novel Starobinsky-Bel-Robinson gravity model admits exponential inflationary solutions with or without spatial anisotropies. As a result, we are able to derive exact de Sitter as well as Bianchi type I inflationary solutions to this Starobinsky-Bel-Robinson model. However, stability analysis using the dynamical system approach indicates that both of the obtained inflationary solutions turn out to be unstable. Finally, we point out that a stable de Sitter inflationary solution can be obtained by flipping the coefficient's sign of $R^2$ term in the Starobinsky-Bel-Robinson gravity.

In this work we consider the observational properties of compact boson stars with self-interactions orbited by isotropically emitting (hot-spot) sources and optically thin accretion disks. We consider two families of boson stars supported by quartic and sixth-order self-interaction potentials, and choose three samples of each of them in growing compactness; only those with large enough compactness are capable to hold light-rings, namely, null bound orbits. For the hot-spots, using inclination angles $\theta=\{20^\circ, 50^\circ, 80^\circ \}$ we find a secondary track plunge-through image of photons crossing the interior of the boson star, which can be further decomposed into additional images if the star is compact enough. For accretion disks we find that the latter class of stars actually shows a sequence of additional secondary images in agreement with the hot-spot analysis, a feature absent in typical black hole space-times. Furthermore, we also find a shadow-like central brightness depression for some of these stars in both axial observations and at the inclination angles above. We discuss our findings in relation to the capability of boson stars to effectively act as black hole mimickers in their optical appearances as well as potential observational discriminators.

Gravitational-wave signals from compact binary coalescences are most efficiently identified through matched filter searches, which match the data against a pre-generated bank of gravitational-wave templates. Although different techniques for performing the matched filter, as well as generating the template bank, exist, currently all modelled gravitational-wave searches use templates that restrict the component spins to be aligned (or anti-aligned) with the orbital angular momentum. This means that current searches are less sensitive to gravitational-wave signals generated from binaries with generic spins (precessing), suggesting that, potentially, a significant fraction of signals may remain undetected. In this work we introduce a matched filter search that is sensitive to signals generated from precessing binaries and can realistically be used during a gravitational-wave observing run. We take advantage of the fact that a gravitational-wave signal from a precessing binary can be decomposed into a power series of five harmonics, to show that a generic-spin template bank, which is only $\sim 3\times$ larger than existing aligned-spin banks, is needed to increase our sensitive volume by $\sim 100\%$ for neutron star black hole binaries with total mass larger than $17.5\, M_{\odot}$ and in-plane spins $>0.67$. In fact, our generic spin search performs as well as existing aligned-spin searches for neutron star black hole signals with insignificant in-plane spins, but improves sensitivity by $\sim60\%$ on average across the full generic spin parameter space. We anticipate that this improved technique will identify significantly more gravitational-wave signals, and, ultimately, help shed light on the unknown spin distribution of binaries in the universe.

The physical interpretation of the exact solutions of the Einstein field equations is, in general, a challenging task, part of the difficulties lying in the significance of the coordinate system. We discuss the extension of the International Astronomical Union (IAU) reference system to the exact theory. It is seen that such an extension, retaining some of its crucial properties, can be achieved in a special class of spacetimes, admitting non-shearing congruences of observers which, at infinity, have zero vorticity and acceleration. As applications, we consider the FLRW, Kerr and NUT spacetimes, the van Stockum rotating dust cylinder, spinning cosmic strings and, finally, we debunk the so-called Balasin-Grumiller (BG) model, and the claims that the galaxies' rotation curves can be explained through gravitomagnetic effects without the need for Dark Matter. The BG spacetime is shown to be completely inappropriate as a galactic model: its dust is actually static with respect to the asymptotic inertial frame, its gravitomagnetic effects arise from unphysical singularities along the axis (a pair of NUT rods, combined with a spinning cosmic string), and the rotation curves obtained are merely down to an invalid choice of reference frame -- the congruence of zero angular momentum observers, which are being dragged by the singularities.

Galactic binaries, and notably double white dwarfs systems, will be a prominent source for the future LISA and Einstein Telescope detectors. Contrarily to the black holes observed by the current LIGO-Virgo-KAGRA network, such objects bear intense magnetic fields, that are naturally expected to leave some imprints on the gravitational wave emission. The purpose of this work is thus to study those imprints within the post-Newtonian (PN) framework, particularly adapted to double white dwarfs systems. To this end, we construct an effective action that takes into account the whole electromagnetic structure of a star, and then specify it to dipolar order. With this action at hand, we compute the acceleration and Noetherian quantities for generic electric and magnetic dipoles, at a relative 2PN order. Finally, focusing on physically relevant systems, we show that the magnetic effects on the orbital frequency, energy and angular momentum is significant, confirming previous works conclusions.