Papers for Monday, May 15 2023

Identifying the Milky Way's very high energy hadronic cosmic-ray accelerators -- the PeVatrons -- is a critical problem. While gamma-ray observations reveal promising candidate sources, neutrino detection is needed for certainty, and this has not yet been successful. Why not? There are several possibilities, as we delineated in a recent paper [T. Sudoh and J. F. Beacom, Phys. Rev. D 107, 043002 (2023)]. Here we further explore the possibility that the challenges arise because PeVatrons have a large angular extent, either due to cosmic-ray propagation effects or due to clusters of sources. We show that while extended neutrino sources could be missed in the commonly used muon-track channel, they could be discovered in the all-flavor shower channel, which has a lower atmospheric-neutrino background flux per solid angle. Intrinsically, showers are quite directional and would appear so in water-based detectors like the future KM3NeT, even though they are presently badly smeared by light scattering in ice-based detectors like IceCube. Our results motivate new shower-based searches as part of the comprehensive approach to identifying the Milky Way's hadronic PeVatrons.

The physical conditions giving rise to high escape fractions of ionising radiation (LyC $f_{\rm{esc}}$) in star-forming galaxies - most likely protagonists of cosmic reionisation - are not yet fully understood. Using the properties of the Lyman-$\alpha$ line profile associated with LyC escape, we select potential LyC leakers and non-leakers from a compiled sample of 1422 MUSE-Wide and MUSE HUDF Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<6.7. We perform spectral stacking to obtain high signal-to-noise detections of rest-frame UV absorption and emission lines, and find that the stacks with LyC-leaker candidates show (i) strong nebular OIII]1666, [SiIII]1883 and [CIII]1907+CIII]1909 emission, suggesting high ionisation parameters due to an elevated production rate of ionising photons coming from young, metal-poor stars; (ii) high equivalent widths of HeII1640 (~1-3 A), possibly indicating a hard ionising spectrum alongside with a high ionising photon production efficiency; (iii) SiII*1533 emission, revealing the presence of neutral hydrogen off the line of sight, thus implying a highly anisotropic interstellar medium (ISM); (iv) high CIV1548,1550 to [CIII]1907+CIII]1909 ratios (CIV/CIII] > 0.75), partly associated with the increased ISM transparency. In contrast, the stacks with non-leakers show weaker nebular emission lines, low HeII1640 equivalent widths (<1 A), and low CIV/CIII] (<0.25), suggesting a low ionisation state of the ISM and a high neutral hydrogen content. Consequently, our results substantiate that the CIV/CIII] ratio can be used as an indirect tracer of $f_{\rm{esc}}$, providing a promising tool for identification of ionising sources among star-forming galaxies in the epoch of reionisation.

Recent observations of radio relics - diffuse radio emission in galaxy clusters - have revealed that these sources are not smooth but consist of structures in the form of threads and filaments. We investigate the origin of these filamentary structures and the role of projection effects. To this end, we have developed a tool that extracts the filamentary structures from background emission. Moreover, it is capable of studying both two-dimensional and three-dimensional objects. We apply our structure extractor to, both, observations and cosmological simulations of radio relics. Using Minkowski functionals, we determine the shape of the identified structures. In our 2D analysis, we find that the brightest structures in the observed and simulated maps are filaments. Our analysis of the 3D simulation data shows that radio relics do not consist of sheets but only of filaments and ribbons. Furthermore, we did not find any measurable projection effects that could hide any sheet-like structures in projection. We find that, both, the magnetic field and the shock front consist of filaments and ribbons that cause filamentary radio emission.

We present the state-of-the-art single-zone nuclear reaction network WinNet that is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis where only a few nuclei are considered to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. Once the manuscript has been accepted for publication, WinNet will be publicly available and open source.

The early-science observations made by the James Webb Space Telescope (JWST) have revealed an excess of ultra-massive galaxy candidates that appear to challenge the standard cosmological model ($\Lambda$CDM). Here, we argue that any modifications to $\Lambda$CDM that can produce such ultra-massive galaxies in the early Universe would also affect the UV galaxy luminosity function (UV LF) inferred from the Hubble Space Telescope (HST). The UV LF covers the same redshifts ($z\approx 7-10$) and host-halo masses $(M_\mathrm{h}\approx 10^{10}-10^{12}\, M_\odot$) as the JWST candidates, but tracks star-formation rate rather than stellar mass. We consider beyond-$\Lambda$CDM power-spectrum enhancements and show that any departure large enough to reproduce the abundance of ultra-massive JWST candidates is in conflict with the HST data. Our analysis, therefore, severely disfavors a cosmological explanation for the JWST abundance problem. Looking ahead, we determine the maximum allowable stellar-mass function and provide projections for the high-$z$ UV LF given our constraints on cosmology from current HST data.

Stars in the Galactic disc, including the Solar system, have deviated from their birth orbits and have experienced radial mixing and vertical heating. By performing hydrodynamical simulations of a galactic disc, we investigate how much tracer particles, which are initially located in the disc to mimic newborn stars and the thin and thick disc stars, are displaced from initial near-circular orbits by gravitational interactions with giant molecular clouds (GMCs). To exclude the influence of other perturbers that can change the stellar orbits, such as spiral arms and the bar, we use an axisymmetric form for the entire galactic potential. First, we investigate the time evolution of the radial and vertical velocity dispersion $\sigma_R$ and $\sigma_z$ by comparing them with a power law relation of $\sigma \propto t^{\beta}$. Although the exponents $\beta$ decrease with time, they keep large values of 0.3 $\sim$ 0.6 for 1 Gyr, indicating fast and efficient disc heating. Next, we find that the efficient stellar scattering by GMCs also causes a change in angular momentum for each star and, therefore, radial migration. This effect is more pronounced in newborn stars than old disc stars; nearly 30 per cent of stars initially located on the galactic mid-plane move more than 1 kpc in the radial direction for 1 Gyr. The dynamical heating and radial migration drastically occur in the first several hundred Myr. As the amplitude of the vertical oscillation increases, the time spent in the galactic plane, where most GMCs are distributed, decreases, and the rate of an increase in the heating and migration slows down.

Aims. Deep millimeter surveys are necessary to probe the dust-obscured galaxies at high redshift. We conducted a large observing program at 1.2 and 2 mm with the NIKA2 camera installed on the IRAM 30-meter telescope. This NIKA2 Cosmological Legacy Survey (N2CLS) covers two emblematic fields: GOODS-N and COSMOS. We introduce the N2CLS survey and present new 1.2 and 2 mm number count measurements based on the tiered N2CLS observations from October 2017 to May 2021. Methods. We develop an end-to-end simulation that combines an input sky model with the instrument noise and data reduction pipeline artifacts. This simulation is used to compute the sample purity, flux boosting, pipeline transfer function, completeness, and effective area of the survey. We used the 117 deg$^2$ SIDES simulations as the sky model, which include the galaxy clustering. Our formalism allows us to correct the source number counts to obtain galaxy number counts, the difference between the two being due to resolution effects caused by the blending of several galaxies inside the large beam of single-dish instruments. Results. The N2CLS-May2021 survey reaches an average 1-$\sigma$ noise level of 0.17 and 0.048 mJy on GOODS-N over 159 arcmin$^2$, and 0.46 and 0.14 mJy on COSMOS over 1010 arcmin$^2$, at 1.2 and 2 mm, respectively. For a purity threshold of 80%, we detect 120 and 67 sources in GOODS-N and 195 and 76 sources in COSMOS, at 1.2 and 2 mm, respectively. Our measurement connects the bright single-dish to the deep interferometric number counts. After correcting for resolution effects, our results reconcile the single-dish and interferometric number counts and are further accurately compared with model predictions.

Massive neutrinos modify the expansion history of the universe and suppress the structure formation below their free streaming scale. Cosmic microwave background (CMB) observations at small angular scales can be used to constrain the total mass $\Sigma m_\nu$ of the three neutrino flavors. However, at these scales, the CMB-measured $\Sigma m_\nu$ is degenerate with $\tau$, the optical depth to reionization, which quantifies the damping of CMB anisotropies due to the scattering of CMB photons with free electrons along the line of sight. Here we revisit the idea to use 21-cm power spectrum observations to provide direct estimates for $\tau$. A joint analysis of CMB and 21-cm data can alleviate the $\tau-\Sigma m_\nu$ degeneracy, making it possible to measure $\Sigma m_\nu$ with unprecedented precision. Forecasting for the upcoming Hydrogen Epoch of Reionization Array (HERA), we find that a $\lesssim\mathcal{O}(10\%)$ measurement of $\tau$ is achievable, which would enable a $\gtrsim 5\sigma$ measurement of $\Sigma m_\nu=60\,[{\rm meV}]$, for any astrophysics model that we considered. Precise estimates of $\tau$ also help reduce uncertainties in other cosmological parameters, such as $A_s$, the amplitude of the primordial scalar fluctuations power spectrum.

Notwithstanding the tremendous growth of the exoplanetary field in the last decade, limited attention has been paid to the planets around binary stars. Circumbinary planets (CBPs) have been discovered primarily around Main Sequence (MS) stars. No exoplanet has been found orbiting double white dwarf (DWD) binaries yet. We modelled the long-term evolution of CBPs, throughout the life stages of their hosts, from MS to white dwarf (WD). Our goal is to provide the community with both theoretical constraints on CBPs evolution beyond the MS and the occurrence rates of planet survival. We further developed the publicly available Triple Evolution Simulation (TRES) code, implementing a variety of physical processes affecting substellar bodies. We then used this code to simulate the evolution, up to one Hubble time, of two synthetic populations of circumbinary giant planets. Each population has been generated using different priors for the planetary orbital parameters. In our simulated populations we identified several evolutionary categories, such as survived, merged, and destabilised systems. Our primary focus is those systems where the planet survived the entire system evolution and orbits a DWD binary, which we call "Magrathea" planets. We found that a significant fraction of simulated CBPs survive and become Magratheas. In the absence of multi-planet migration mechanisms, this category of planets is characterised by long orbital periods. Magrathea planets are a natural outcome of triple systems evolution, and they could be relatively common in the Galaxy. They can survive the death of their binary hosts if they orbit far enough to avoid engulfment and instabilities. Our results can ultimately be a reference to orient future observations of this uncharted class of planets and to compare different theoretical models.

Non-Gaussian distributions in cosmology are commonly evaluated with Monte Carlo Markov-chain methods, as the Fisher-matrix formalism is restricted to the Gaussian case. The Metropolis-Hastings algorithm will provide samples from the posterior distribution after a burn-in period, and the corresponding convergence is usually quantified with the Gelman-Rubin criterion. In this paper, we investigate the convergence of the Metropolis-Hastings algorithm by drawing analogies to statistical Hamiltonian systems in thermal equilibrium for which a canonical partition sum exists. Specifically, we quantify virialisation, equipartition and thermalisation of Hamiltonian Monte Carlo Markov-chains for a toy-model and for the likelihood evaluation for a simple dark energy model constructed from supernova data. We follow the convergence of these criteria to the values expected in thermal equilibrium, in comparison to the Gelman-Rubin criterion. We find that there is a much larger class of physically motivated convergence criteria with clearly defined target values indicating convergence. As a numerical tool, we employ physics-informed neural networks for speeding up the sampling process.

Within the convection zone of a rotating star, the presence of the Coriolis force stabilizes long-wavelength convective modes. These modes, which would have been unstable if the star lacked rotation, are called overstable convective modes or thermal Rossby waves. We demonstrate that the Sun's rotation rate is sufficiently rapid that the lower half of its convection zone could possess overstable modes. Further, we present an analytic solution for atmospheric waves that reside within a polytropic stratification. We explore in detail the properties of the overstable and unstable wave modes that exist when the polytrope is weakly unstable to convective overturning. Finally, we discuss how the thermal Rossby waves that reside within the convection zone of a star might couple with the prograde branch of the $g$ modes that are trapped within the star's radiative zone. We suggest that such coupling might enhance the photospheric visibility of a subset of the Sun's $g$ modes.

We present a spectroscopic analysis of a sample of 48 M dwarf stars ($0.2 M_{\odot}< M < 0.6 M_{\odot}$) from the Hyades open cluster using high-resolution H-band spectra from the SDSS/APOGEE survey. Our methodology adopts spectrum synthesis with LTE MARCS model atmospheres, along with the APOGEE DR17 line list, to determine effective temperatures, surface gravities, metallicities, and projected rotational velocities. The median metallicity obtained for the Hyades M dwarfs is [M/H]= 0.09$\pm$0.03 dex, indicating a small internal uncertainty and good agreement with optical results for Hyades red-giants. Overall, the median radii are larger than predicted by stellar models by 1.6$\pm$2.3\% and 2.4$\pm$2.3\%, relative to a MIST and DARTMOUTH isochrone, respectively. We emphasize, however, that these isochrones are different and the fractional radius inflation for the fully- and partially-convective regimes have distinct behaviors depending on the isochrone. Using a MIST isochrone there is no evidence of radius inflation for the fully convective stars, while for the partially convective M-dwarfs the radii are inflated by 2.7$\pm$2.1\%, which is in agreement with predictions from models that include magnetic fields. For the partially-convective stars, rapid-rotators present on average higher inflation levels than slow-rotators. The comparison with SPOTS isochrone models indicates that the derived M dwarf radii can be explained by accounting for stellar spots in the photosphere of the stars, with 76\% of the studied M dwarfs having up to 20\% spot coverage, and the most inflated stars with $\sim$20 -- 40\% spot coverage.

The recent observation of a low-mass $z=5.2$ and an intermediate-mass $z=7.3$ (JADES-GS-z7-01-QU) quenched galaxy with JWST / NIRSpec is the first evidence of halted star formation above $z\sim 5$. Here we show how bursty star formation at high redshift gives rise to temporarily quenched, or miniquenched galaxies in the mass range $M_{\star} = 10^7-10^9 \ M_{\odot}$ using three models of galaxy formation: the periodic box simulation IllustrisTNG, the zoom-in simulation VELA and an empirical halo model. The main causes for mini-quenching are stellar feedback, lack of gas accretion onto galaxies and galaxy-galaxy interactions. The abundance of mini-quenching events agrees across the three models: the population first appears below $z\sim 8$, after which the fraction of miniquenched galaxies increases with cosmic time, from $\sim 0.5$% at $z=7$ to $\sim 1-2$% at $z=4$, corresponding to comoving number densities of $8.0\times 10^{-6}$ Mpc$^{-3}$ and $5.4\times 10^{-4}$ Mpc$^{-3}$, respectively. The star formation rate duty cycle ($f_{\mathrm{duty}}\sim 99.56^{+0.4}_{-4.5}$% at $z=7$) inferred for VELA galaxies is consistent therewith. Star formation histories (SFHs) in VELA suggest that mini-quenching at $z=4-8$ is short-lived with a duration of $\sim 20-40$ Myr, which is close to the free-fall timescale of the inner halo. However, mock spectral energy distributions of miniquenched galaxies in IllustrisTNG and VELA do not match JADES-GS-z7-01-QU photometry, unless their SFHs are artificially altered to be more bursty on timescales of $\sim 40$ Myr. Studying miniquenched galaxies might aid in calibrating the sub-grid models governing galaxy formation, as these may not generate sufficient burstiness at high redshift to explain the SFH inferred for JADES-GS-z7-01-QU.

We report on a timing and spectral analysis of a 50-ks NuSTAR observation of IGR J16320-4751 (= AX J1631.9-4752); a high-mass X-ray binary hosting a slowly-rotating neutron star. In this observation from 2015, the spin period was 1,308.8+/-0.4 s giving a period derivative dP/dt ~ 2E-8 s s-1 when compared with the period measured in 2004. In addition, the pulsed fraction decreased as a function of energy, as opposed to the constant trend that was seen previously. This suggests a change in the accretion geometry of the system during the intervening 11 years. The phase-averaged spectra were fit with the typical model for accreting pulsars: a power law with an exponential cutoff. This left positive residuals at 6.4 keV attributable to the known iron K-alpha line, as well as negative residuals around 14 keV from a candidate cyclotron line detected at a significance of 5-sigma. We found no significant differences in the spectral parameters across the spin period, other than the expected changes in flux and component normalizations. A flare lasting around 5 ks was captured during the first half of the observation where the X-ray emission hardened and the local column density decreased. Finally, the binary orbital period was refined to 8.9912+/-0.0078 d thanks to Swift/BAT monitoring data from 2005-2022.

Using the data of the SRG/eROSITA all-sky survey, we stacked a sample of ~40 galaxy cluster images in the 0.3--2.3 keV band, covering the radial range up to $10\times R_{\rm 500c}$. The excess emission on top of the galactic and extragalactic X-ray backgrounds and foregrounds is detected up to $\sim 3\times R_{\rm 500c}$. At these distances, the surface brightness of the stacked image drops below $\sim 1$% of the background. The density profile reconstructed from the X-ray surface brightness profile agrees well (within $\sim30$%) with the mean gas profile found in numerical simulations, which predict the local gas overdensity of $\sim$ 20--30 at $3\times R_{\rm 500c}$ and the gas fraction close to the universal value of $\frac{\Omega_b}{\Omega_m}\approx 0.15$ in the standard $\Lambda$CDM model. Taking at face value, this agreement suggests that up to $\sim 3\times R_{\rm 500c}$ the X-ray signal is not strongly boosted by the gas clumpiness, although a scenario with a moderately inhomogeneous gas cannot be excluded. A comparison of the derived gas density profile with the electron pressure profile based on the SZ measurements suggests that by $r\sim 3\times R_{\rm 500c}$ the gas temperature drops by a factor of $\sim$ 4--5 below the characteristic temperature of a typical cluster in the sample within $R_{\rm 500c}$, while the entropy keeps growing up to this distance. Better constraints on the gas properties just beyond $3\times R_{\rm 500c}$ should be possible with a sample larger than used for this pilot study.

The IceCube project transformed a cubic kilometer of transparent, natural Antarctic ice into a Cherenkov detector. It discovered neutrinos of TeV-PeV energy originating beyond our Galaxy with an energy flux that exceeds the one of high-energy gamma rays of extragalactic origin. Unlike at any other wavelength of light, extragalactic neutrinos outshine the nearby sources in our own Milky way. Updated measurements of the diffuse cosmic neutrino flux indicate that the high-energy gamma rays produced by the neutral pions that accompany cosmic neutrinos lose energy in the sources and are likely to be observed at MeV energy, or below. After the reanalysis of 10 years of archival data with an improved data selection and enhanced data analysis methods, the active galaxy NGC 1068 emerged as the hottest spot in the neutrino sky map. It is also the most significant source in a search at the positions of 110 preselected high-energy gamma-ray sources. Additionally, we find evidence for neutrino emission from the active galaxies PKS 1424+240 and TXS 0506+056. TXS 0506+056 had already been identified as a neutrino source in a multimessenger campaign triggered by a neutrino of 290 TeV energy and, by the independent observation of a neutrino burst in 2014 from this source in archival IceCube data. The observations point to active galaxies as the sources of cosmic neutrinos, and cosmic rays, with the gamma-ray-obscured dense cores near the supermassive black holes at their center as the sites where neutrinos originate, typically within $10\sim100$ Schwarzschild radii.

We present results of long-term photometric monitoring of SS433 which proves a secular evolutionary increase of the orbital period of SS433 at a rate of $(1.14\pm 0.25)\times 10^{-7}$ s~s$^{-1}$. Using a physical model of non-conservative mass transfer in SS433 through a supercritical accretion disc around the compact companion, we reliably confirm that the binary mass ratio in SS433, $q=M_x/M_v$ is $\gtrsim 0.8$. For an optical star mass $M_v\sim 10 M_\odot$ the compact object in SS433 is a black hole with mass $M_{BH}\gtrsim 8 M_\odot$. We discuss evolutionary implications of the found orbital period increase in SS433 -- a secular change in the orbital separation and a size of the Roche lobe of the optical star. We show that for the mass-loss rate $dM_v/dt\sim 10^{-4}-3\times 10^{-5} M_\odot$ per year and an optical star mass $M_v \sim 10-15 M_\odot$ the found orbital period increase implies the corresponding orbital separation increase while the Roche lobe size can shrink or expand around a mean constant value depending on the optical star mass-loss rate which may be modulated with the precessional period.

We studied the H$_2$ column density probability distribution function (N-PDF) based on molecular emission lines using the Nobeyama 45-m Cygnus X CO survey data. Using the DENDROGRAM and SCIMES algorithms, we identified 124 molecular clouds in the $^{13}$CO data. From these identified molecular clouds, an N-PDF was constructed for 11 molecular clouds with an extent of more than 0.4 deg$^2$. From the fitting of the N-PDF, we found that the N-PDF could be well-fitted with one or two log-normal distributions. These fitting results provided an alternative density structure for molecular clouds from a conventional picture. We investigated the column density, dense molecular cloud cores, and radio continuum source distributions in each cloud and found that the N-PDF shape was less correlated with the star-forming activity over a whole cloud. Furthermore, we found that the log-normal N-PDF parameters obtained from the fitting showed two impressive features. First, the log-normal distribution at the low-density part had the same mean column density ($\sim$ 10$^{21.5}$ cm$^{-2}$) for almost all the molecular clouds. Second, the width of the log-normal distribution tended to decrease with an increasing mean density of the structures. These correlations suggest that the shape of the N-PDF reflects the relationship between the density and turbulent structure of the whole molecular cloud but is less affected by star-forming activities.

We present the results of spectroscopic observations of two eclipsing WR+OB-type systems - CQ Cep and CX Cep, performed in 2020-2023 with a low-resolution slit spectrograph TDS ($\lambda\lambda= 3660-7410$\AA, $R = 1300-2500$) on 2.5-m telescope of the SAI MSU Caucasian Mountain Observatory. For CQ Cep, the radial velocity curves of a WN6 star are constructed, the problem of visibility of spectroscopic traces of an OB star is discussed and the components' mass ratio $q\sim 0.6$ is estimated. For CX Cep, the radial velocity curves are constructed for both the WN5 and O5 components enabling their masses and circular orbit elements to be refined. The comparison of the radial velocity curves of these systems obtained in different epochs allowed us to derive the orbital period change rate $\dot{P}$ by the spectroscopic method, which is found to be in good agreement with estimates obtained by comparing the moments of primary eclipse minima: $\dot{P} = -0.0151\pm0.0013$ s yr$^{-1}$ for CQ Cep and $\dot{P} = 0.054\pm0.009$ s yr$^{-1}$ for CX Cep. The prospects of applicability of the spectroscopic dynamical method for studying the orbital evolution of Galactic WR+OB binaries and related objects are considered. We also discuss the effect of finite sizes of stars with stellar wind mass loss in close binary systems on their orbital evolution.

While the abundance of elemental deuterium is relatively low (D/H ~ a few 1E-5), orders of magnitude higher D/H abundance ratios have been found for many interstellar molecules, enhanced by deuterium fractionation. In cold molecular clouds (T < 20K) deuterium fractionation is driven by the H2D+ ion, whereas at higher temperatures (T > 20-30K) gas-phase deuteration is controlled by reactions with CH2D+ and C2HD+. While the role of H2D+ in driving cold interstellar deuterium chemistry is well understood, thanks to observational constraints from direct measurements of H2D+, deuteration stemming from CH2D+ is far less understood, caused by the absence of direct observational constraints of its key ions. Therefore, making use of chemical surrogates is imperative for exploring deuterium chemistry at intermediate temperatures. Formed at an early stage of ion-molecule chemistry, directly from the dissociative recombination of CH3+ (CH2D+), CH (CD) is an ideal tracer for investigating deuterium substitution initiated by reactions with CH2D+. This paper reports the first detection of CD in the interstellar medium, carried out using the APEX 12m telescope toward the widely studied low-mass protostellar system IRAS 16293-2422. Gas-phase chemical models reproducing the observed CD/CH abundance ratio of 0.016 suggests that it reflects `warm deuterium chemistry' (which ensues in moderately warm conditions of the interstellar medium) and illustrates the potential use of the CD/CH ratio in constraining the gas temperatures of the envelope gas clouds it probes.

We aim at characterizing the low-mass (sub)stellar population of the central portion (2.4 pc$^2$) of the $\sim$2 Myr old cluster NGC 2244 using near infrared spectroscopy. By studying this cluster, characterized by a low stellar density and numerous OB stars, we aim at exploring the effect that OB stars may have on the production of BDs. We obtain near infrared HK spectroscopy of 85 faint candidate members of NGC 2244. We derive the spectral type and extinction by comparison with spectral templates. We evaluate cluster membership using three gravity-sensitive spectral indices based on the shape of the $H$-band. Furthermore, we evaluate the infrared excess from Spitzer of all the candidate members of the cluster. Finally, we estimate the mass of all the candidate members of the cluster and derive the initial mass function, star-to-BD number ratio and disk fraction. The initial mass function is well represented by a power law ($dN/dM \propto M^{-\alpha}$) below 0.4 $M_\odot$, with a slope $\alpha$ = 0.7-1.1 depending on the fitted mass range. We calculate a star-to-BD number ratio of 2.2-2.8. We find the low-mass population of NGC 2244 to be consistent with nearby star-forming regions, although it is at the high-end of BD production. We find BDs in NGC 2244 to be on average closer to OB stars than to low-mass stars, which could potentially be the first evidence of OB stars affecting the formation of BDs. We find a disk fraction of all the members with spectral type later than K0 of 39$\pm$9% which is lower than typical values found in nearby star-forming regions of similar ages.

The Sun and other solar-type stars have magnetic fields that permeate their interior and surface, extends through the interplanetary medium, and is the main driver of stellar activity. Stellar magnetic activity affects physical processes and conditions of the interplanetary medium and orbiting planets. Coronal mass ejections (CMEs) are the most impacting of these phenomena in near-Earth space weather, and consist of plasma clouds, with magnetic field, ejected from the solar corona. Precisely predicting the trajectory of CMEs is crucial in determining whether a CME will hit a planet and impact its magnetosphere and atmosphere. Despite the rapid developments in the search for stellar CMEs, their detection is still very incipient. In this work we aim to better understand the propagation of CMEs by analysing the influence of initial parameters on CME trajectories, such as position, velocities, and stellar magnetic field's configuration. We reconstruct magnetograms for Kepler-63 (KIC 11554435) and Kepler-411 (KIC 11551692) from spot transit mapping, and use a CME deflection model, ForeCAT, to simulate trajectories of hypothetical CMEs launched into the interplanetary medium from Kepler-63 and Kepler-411. We apply the same methodology to the Sun, for comparison. Our results show that in general, deflections and rotations of CMEs decrease with their radial velocity, and increase with ejection latitude. Moreover, magnetic fields stronger than the Sun's, such as Kepler-63's, tend to cause greater CME deflections.

The zeta ring is the innermost component of the Uranian ring system. It is of scientific interest because its morphology changed significantly between the Voyager 2 encounter in 1986 and subsequent Earth-based observations around 2007. It is also of practical interest because some Uranus mission concepts have the spacecraft pass through the inner flank of this ring. Recent re-examinations of the Voyager 2 images have revealed additional information about this ring that provide a more complete picture of the ring's radial brightness profile and phase function. These data reveal that this ring's brightness varies with phase angle in a manner similar to other tenuous rings, consistent with it being composed primarily of sub-millimeter-sized particles. The total cross section of particles within this ring can also be estimated from these data, but translating that number into the actual risk to a spacecraft flying through this region depends on a number of model-dependent parameters. Fortunately, comparisons with Saturn's G and D rings allows the zeta-ring's particle number density to be compared with regions previously encountered by the Voyager and Cassini spacecraft. Finally, these data indicate that the observed changes in the zeta-ring's structure between 1986 and 2007 are primarily due to a substantial increase in the amount of dust at distances between 38,000 km and 40,000 km from Uranus' center.

GW170817 is the unique gravitational-wave (GW) event that is associated to the electromagnetic (EM) counterpart GRB 170817A. NGC 4993 is identified as the host galaxy of GW170817/GRB 170817A. In this paper, we particularly focus on the spatially resolved properties of NGC 4993. We present the photometric results from the comprehensive data analysis of the high spatial-resolution images in the different optical bands. The morphological analysis reveals that NGC 4993 is a typical early-type galaxy without significant remnants of major galaxy merger. The spatially resolved stellar population properties of NGC 4993 suggest that the galaxy center has passive evolution with the outskirt formed by gas accretion. We derive the merging rate of the compact object per galaxy by a co-evolution scenario of supermassive black hole and its host galaxy. If the galaxy formation is at redshift 1.0, the merging rate per galaxy is $3.2\times 10^{-4}$ to $7.7\times 10^{-5}$ within the merging decay time from 1.0 to 5.0 Gyr. The results provide the vital information for the ongoing GW EM counterpart detections. The HST data analysis presented in this paper can be also applied for the Chinese Space Station Telescope (CSST) research in the future.

This manuscript develops a simultaneous navigation and gravity estimation strategy around a small body. The scheme combines dynamical model compensation with a mascon gravity fit. Dynamical compensation adds the unmodeled acceleration to the filter state. Consequently, the navigation filter is able to generate an on-orbit position-unmodeled acceleration dataset. The available measurements correspond to the landmarks-based navigation technique. Accordingly, an on-board camera is able to provide landmark pixels. The aforementioned position-unmodeled acceleration dataset serves to train a mascon gravity model on-board while in flight. The training algorithm finds the optimal mass values and locations using Adam gradient descent. By a careful choice of the mascon variables and constraints projection, the masses are ensured to be positive and within the small body shape. The numerical results provide a comprehensive analysis on the global gravity accuracy for different estimation scenarios.

Thorne-$\.Z$ytkow objects (T$\.Z$O) are potential end products of the merger of a neutron star with a non-degenerate star. In this work, we have computed the first grid of evolutionary models of T$\.Z$Os with the MESA stellar evolution code. With these models, we predict several observational properties of T$\.Z$Os, including their surface temperatures and luminosities, pulsation periods, and nucleosynthetic products. We expand the range of possible T$\.Z$O solutions to cover $3.45 \lesssim \log \left(T/K\right) \lesssim 3.65$ and $4.85 \lesssim \log \left(L/L_{\odot}\right) \lesssim 5.5$. Due to the much higher densities our T$\.Z$Os reach compared to previous models, if T$\.Z$Os form we expect them to be stable over a larger mass range than previously predicted, without exhibiting a gap in their mass distribution. Using the GYRE stellar pulsation code we show that T$\.Z$Os should have fundamental pulsation periods of 1000--2000 days, and period ratios of $\approx$0.2--0.3. Models computed with a large 399 isotope fully-coupled nuclear network show a nucleosynthetic signal that is different to previously predicted. We propose a new nucleosynthetic signal to determine a star's status as a T$\.Z$O: the isotopologues $^{44}\rm{Ti} \rm{O}_2$ and $^{44}\rm{Ti} \rm{O}$, which will have a shift in their spectral features as compared to stable titanium-containing molecules. We find that in the local Universe (~SMC metallicities and above) T$\.Z$Os show little heavy metal enrichment, potentially explaining the difficulty in finding T$\.Z$Os to-date.

Context. Filamentary structures in the interstellar medium are closely related to star formation. Dense gas mass fraction (DGMF) or clump formation efficiency in large-scale filaments possibly determine their hosting star formation activities. Aims. We aim to automatically identify large-scale filaments, characterize them, investigate their association with Galactic structures, and study their DGMFs. Methods. We use a modified minimum spanning tree (MST) algorithm to chain parsec-scale 13CO clumps previously extracted from the SEDIGISM (Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium) survey. The MST connects nodes in a graph such that the sum of edge lengths is minimum. Modified MST also ensures velocity coherence between nodes, so the identified filaments are coherent in position-position-velocity (PPV) space. Results. We generate a catalog of 88 large-scale ($>10pc$) filaments in the inner Galactic plane (with $-60^\circ < l < 18^\circ and $|b| < 0.5^\circ$). These SEDIGISM filaments are larger and less dense than MST filaments previously identified from the BGPS and ATLASGAL surveys. We find that eight of the filaments run along spiral arms and can be regarded as "bones" of the Milky Way. We also find three bones associated with the Local Spur in PPV space. By compiling 168 large-scale filaments with available DGMF across the Galaxy, an order of magnitude more than previously investigated, we find that DGMFs do not correlate with Galactic location, but bones have higher DGMFs than other filaments.

The Gaia-Multi-Peak (GMP) technique can identify large numbers of dual or lensed active galactic nuclei (AGN) candidates at sub-arcsec separation, allowing us to study both multiple super-massive black holes (SMBH) in the same galaxy and rare, compact lensed systems. The observed samples can be used to test the predictions of the models of SMBH merging once 1) the selection function of the GMP technique is known, and 2) each system has been classified as dual AGN, lensed AGN, or AGN/star alignment. Here we show that the GMP selection is very efficient for separations above 0.15" when the secondary (fainter) object has magnitude G<20.5. We present the spectroscopic classification of five GMP candidates using VLT/ERIS and Keck/OSIRIS, and compare them with the classifications obtained from: a) the near-IR colors of 7 systems obtained with LBT/LUCI, and b) the analysis of the total, spatially-unresolved spectra. We conclude that colors and integrated spectra can already provide reliable classifications of many systems. Finally, we summarize the 14 confirmed dual AGNs at z>0.5 selected by the GMP technique, and compare this sample with other such systems from the literature, concluding that GMP can provide a large number of confirmed dual AGNs at separations below 7 kpc.

Astrophysical observations of neutron stars probe the structure of dense nuclear matter and have the potential to reveal phase transitions at high densities. Most recent analyses are based on parametrized models of the equation of state with a finite number of parameters and occasionally include extra parameters intended to capture phase transition phenomenology. However, such models restrict the types of behavior allowed and may not match the true equation of state. We introduce a complementary approach that extracts phase transitions directly from the equation of state without relying on, and thus being restricted by, an underlying parametrization. We then constrain the presence of phase transitions in neutron stars with astrophysical data. Current pulsar mass, tidal deformability, and mass-radius measurements disfavor only the strongest of possible phase transitions (latent energy per particle $\gtrsim 100\,\mathrm{MeV}$). Weaker phase transitions are consistent with observations. We further investigate the prospects for measuring phase transitions with future gravitational-wave observations and find that catalogs of \result{$O(100)$} events will (at best) yield Bayes factors of $\sim 10:1$ in favor of phase transitions even when the true equation of state contains very strong phase transitions. Our results reinforce the idea that neutron star observations will primarily constrain trends in macroscopic properties rather than detailed microscopic behavior. Fine-tuned equation of state models will likely remain unconstrained in the near future.

The effect of interplanetary plasma on pulsed pulsar radiation passing through is considered. The pulses of two rotating radio transients (J0609+16, J1132+25) and a pulsar (B0320+39) detected on the Large Phased Array (Pushchino observatory) were analyzed. It is shown that in observations at the frequency of 111 MHz, on elongations of 20o-40o, both an increase and a decrease in the number of received pulses are observed. The change in the number of pulses is explained by the distortion of the energy distribution of pulses due to interplanetary scintillation. These changes in the number of observed pulses are in qualitative agreement with the expected dependence of the scintillation index on the observed sources elongation. Analytical expressions are obtained that allow estimating the effective modulation index from observations of individual pulses for the power distribution of pulses by energy.

Phase-space data, chemistry, and ages together reveal a complex structure in the outer low-${\alpha}$ disc of the Milky Way. The age-vertical velocity dispersion profiles beyond the Solar Neighbourhood show a significant jump at 6 Gyr for stars beyond the Galactic plane. Stars older than 6 Gyr are significantly hotter than younger stars. The chemistry and age histograms reveal a bump at [Fe/H] = -0.5, [${\alpha}$/Fe] = 0.1, and an age of 7.2 Gyr in the outer disc. Finally, viewing the stars beyond 13.5 kpc in the age-metallicity plane reveals a faint streak just below this bump, towards lower metallicities at the same age. Given the uncertainty in age, we believe these features are linked and suggest a pericentric passage of a massive satellite 6 Gyr ago that heated pre-existing stars, led to a starburst in existing gas. New stars also formed from the metal-poorer infalling gas. The impulse approximation was used to characterise the interaction with a satellite, finding a mass of ~1e11 M$_{\odot}$, and a pericentric position between 12 and 16 kpc. The evidence points to an interaction with the Sagittarius dwarf galaxy, likely its first pericentric passage.

The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to perform maximum a posteriori estimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.

$\gamma$-ray-emitting narrow-line Seyfert 1 galaxies ($\gamma$-NLS1) constitute an intriguing small population of Active Galactic Nuclei with $\gamma$-ray emission resembling low power flat-spectrum radio quasars (FSRQ), but with differing physical properties. They are jetted, $\gamma$/radio-loud Seyfert galaxies, with relatively low black hole masses, accreting at exceptionally high, near-Eddington rates. Certain of these sources exhibit highly variable emission states on relatively short time scales, the physical origin of which remains elusive. In this work, varying emission states of two bona-fide NLS1s, 1H 0323+342 and PMN J0948+0022, and one little studied FSRQ/intermediate object, B2 0954+25A, are examined. For each source, we analyzed quasi-simultaneous multiwavelength data for different states of $\gamma$-ray activity and present the results of their broad-band emission modelling, taking into account all available physical constraints to limit the range of the model parameters. Two different scenarios are discussed, in the framework of a one-zone leptonic model, where the high energy emission is due to the inverse Compton scattering of the disc and broad line region (BLR) or torus photons by relativistic electrons within the jet. The transition from low to high state is well described by variations of the jet parameters, leaving the external photon fields unchanged. The parameterisation favours an emission scenario with particle injection on a stationary shock inside the jet. When considering all physical constraints, the disc & BLR scenario is preferred for all three sources. We use the multi-epoch modelling to characterize total jet powers and discuss the intrinsic nature of $\gamma$-NLS1 galaxies and FSRQs.

We have measured the angular auto-correlation function of near-infrared galaxies in SERVS+DeepDrill, the Spitzer Extragalactic Representative Volume Survey and its follow-up survey of the Deep Drilling Fields, in three large fields totalling over 20 sq. deg on the sky, observed in two bands centred on 3.6 and 4.5 micron. We performed this analysis on the full sample as well as on sources selected by [3.6]-[4.5] colour in order to probe clustering for different redshift regimes. We estimated the spatial correlation strength as well, using the redshift distribution from S-COSMOS with the same source selection. The strongest clustering was found for our bluest subsample, with z~0.7, which has the narrowest redshift distribution of all our subsamples. We compare these estimates to previous results from the literature, but also to estimates derived from mock samples, selected in the same way as the observational data, using deep light-cones generated from the SHARK semi-analytical model of galaxy formation. For all simulated (sub)samples we find a slightly steeper slope than for the corresponding observed ones, but the spatial clustering length is comparable in most cases.

We present a magnetic-field surface map for both stellar components of the young visual binary ksi Boo AB (A: G8V, B: K5V). Employed are high resolution Stokes-V spectra obtained with the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) at the Large Binocular Telescope (LBT). Stokes V line profiles are inverted with our iMAP software and compared to previous inversions. We employed an iterative regularization scheme without the need of a penalty function and incorporated a three-component description of the surface magnetic-field vector. The spectral resolution of our data is 130,000 (0.040-0.055A) and have signal-to-noise ratios (S/N) of up to three thousand per pixel depending on wavelength. A singular-value decomposition (SVD) of a total of 1811 spectral lines is employed for averaging Stokes-V profiles. Our mapping is accompanied by a residual bootstrap error analysis. Magnetic flux densities of the radial field component of up to plus/minus 115 +/- 5 G were reconstructed for ksi Boo A while up to plus/minus 55 +/- 3G were reconstructed for ksi Boo B. ksi Boo A's magnetic morphology is characterized by a very high latitude, nearly polar, spot of negative polarity and three low-to-mid latitude spots of positive polarity while ksi Boo B's morphology is characterized by four low-to-mid latitude spots of mixed polarity. No polar magnetic field is reconstructed for the cooler ksi Boo B star. Both our maps are dominated by the radial field component, containing 86 and 89 percent of the magnetic energy of ksi Boo A and B, respectively. We found only weak azimuthal and meridional field densities on both stars (plus/minus 15-30 G), about a factor two weaker than what was seen previously for ksi Boo A. The phase averaged longitudinal field component and dispersion is +4.5 +/- 1.5G for ksi Boo A and -5.0 +/- 3.0 G for ksi Boo B.

Post-starburst galaxies (PSBs) are transition galaxies showing evidence of recent rapid star formation quenching. To understand the role of galaxy mergers in triggering quenching, we investigate the incidence of PSBs and resolved PSB properties in post-merger galaxies using both SDSS single-fiber spectra and MaNGA resolved IFU spectra. We find post-mergers have a PSB excess of 10 - 20 times that relative to their control galaxies using single-fiber PSB diagnostics. A similar excess of ~ 19 times is also found in the fraction of central (C)PSBs and ring-like (R)PSBs in post-mergers using the resolved PSB diagnostic. However, 60% of the CPSBs + RPSBs in both post-mergers and control galaxies are missed by the single-fiber data. By visually inspecting the resolved PSB distribution, we find that the fraction of outside-in quenching is 7 times higher than inside-out quenching in PSBs in post-mergers while PSBs in control galaxies do not show large differences in these quenching directions. In addition, we find a marginal deficit of HI gas in PSBs relative to non-PSBs in post-mergers using the MaNGA-HI data. The excesses of PSBs in post-mergers suggest that mergers play an important role in triggering quenching. Resolved IFU spectra are important to recover the PSBs missed by single-fiber spectra. The excess of outside-in quenching relative to inside-out quenching in post-mergers suggests that AGN are not the dominant quenching mechanism in these galaxies, but that processes from the disk (gas inflows/consumption and stellar feedback) play a more important role.

Torsional Alfv\'en waves in coronal plasma loops are usually considered to be non-collective, i.e. consist of cylindrical surfaces evolving independently, which significantly complicates their detection in observations. This non-collective nature, however, can get modified in the nonlinear regime. To address this question, the propagation of nonlinear torsional Alfv\'en waves in straight magnetic flux tubes has been investigated numerically using the astrophysical MHD code Athena++ and analytically, to support numerical results, using the perturbation theory up to the second order. Numerical results have revealed that there is radially uniform induced density perturbation whose uniformity does not depend on the radial structure of the mother Alfv\'en wave. Our analysis showed that the ponderomotive force leads to the induction of the radial and axial velocity perturbations, while the mechanism for the density perturbation is provided by a non-equal elasticity of a magnetic flux tube in the radial and axial directions. The latter can be qualitatively understood by the interplay between the Alfv\'en wave perturbations, external medium, and the flux tube boundary conditions. The amplitude of these nonlinearly induced density perturbations is found to be determined by the amplitude of the Alfv\'en driver squared and the plasma parameter $\beta$. The existence of the collective and radially uniform density perturbation accompanying nonlinear torsional Alfv\'en waves could be considered as an additional observational signature of Alfv\'en waves in the upper layers of the solar atmosphere.

We report a 325(-7, +8) day quasi-periodic oscillation (QPO) in the X-ray emission of the blazar Mkn 421, based on data obtained with the Rossi X-ray Timing Explorer (RXTE). The QPO is seen prominently in the ASM data (at least 15 cycles), due to the fact that it has had near-continuous sampling for more than a decade. The PCA data, where the sampling is not uniform and shows many large gaps, provide supporting evidence at lower significance. This QPO is present in both the Proportional Counter Array (PCA) and All-Sky Monitor (ASM) light curves, however it is far more secure (32 sigma significance) in the ASM data since much of the PCA data are from target-of-opportunity flare observations and thus have substantial gaps. QPOs are an important observable in accretion disks, can be modulated by various orbital timescales, and may be generated by a number of mechanisms. They have been studied extensively in X-ray binaries, and should be present in active galactic nuclei (AGN) if they are governed by a common set of physical principles. In jetted sources, QPOs can probe jet-disk interactions or helical oscillations. This QPO previously has been claimed intermittently in X-ray, radio and gamma-ray data, but the continuous, 15-year extent (1996-2011) of the ASM observations (in which Mkn 421 is the brightest AGN observed) provides a unique window. The QPO appears present for nearly the entire extent of the ASM observations. We explore various physical origins and modulating mechanisms, particularly interpretations of the QPO as a result of disk-jet interactions, either due to an accretion disk limit cycle, jet instabilities or helical motions. Limit-cycle related oscillations would not interact with either Keplerian or Lense-Thirring modulated oscillations, however those associated with jet instabilities or helical motions in the jet would likely be modulated by Lense-Thirring precession.

We present the Differentiable Lensing Lightcone (DLL), a fully differentiable physical model designed for being used as a forward model in Bayesian inference algorithms requiring access to derivatives of lensing observables with respect to cosmological parameters. We extend the public FlowPM N-body code, a particle-mesh N-body solver, simulating lensing lightcones and implementing the Born approximation in the Tensorflow framework. Furthermore, DLL is aimed at achieving high accuracy with low computational costs. As such, it integrates a novel Hybrid Physical-Neural parameterisation able to compensate for the small-scale approximations resulting from particle-mesh schemes for cosmological N-body simulations. We validate our simulations in an LSST setting against high-resolution $\kappa$TNG simulations by comparing both the lensing angular power spectrum and multiscale peak counts. We demonstrate an ability to recover lensing $C_\ell$ up to a 10% accuracy at $\ell=1000$ for sources at redshift 1, with as few as $\sim 0.6$ particles per Mpc/h. As a first use case, we use this tool to investigate the relative constraining power of the angular power spectrum and peak counts statistic in an LSST setting. Such comparisons are typically very costly as they require a large number of simulations, and do not scale well with the increasing number of cosmological parameters. As opposed to forecasts based on finite differences, these statistics can be analytically differentiated with respect to cosmology, or any systematics included in the simulations at the same computational cost of the forward simulation. We find that the peak counts outperform the power spectrum on the cold dark matter parameter $\Omega_c$, on the amplitude of density fluctuations $\sigma_8$, and on the amplitude of the intrinsic alignment signal $A_{IA}$.

We present the first X-ray observation at sub-arcsecond resolution of the high-redshift ($z=6.18$) quasar CFHQS J142952+544717 (J1429). The ~100 net-count 0.3-7 keV spectrum obtained from $\sim 30$ ksec Chandra exposure is best fit by a single power-law model with a photon index $\Gamma=2.0\pm0.2$ and no indication of an intrinsic absorber, implying a 3.6-72 keV rest-frame luminosity $L_{\rm X}=(2.3^{+0.6}_{-0.5})\times10^{46}$ erg s$^{-1}$. We identify a second X-ray source at 30 arcsec, distance from J1429 position, with a soft ($\Gamma\simeq 2.8$) and absorbed (equivalent hydrogen column density $N_{\rm H} <13.4\times 10^{20}$ cm$^{-2}$) spectrum, which likely contaminated J1429 spectra obtained in lower angular resolution observations. Based on the analysis of the Chandra image, the bulk of the X-ray luminosity is produced within the central $\sim 3$ kpc region, either by the disk/corona system, or by a moderately aligned jet. In this context, we discuss the source properties in comparison with samples of low- and high-redshift quasars. We find indication of a possible excess of counts over the expectations for a point-like source in a 0.5 arcsec-1.5 arcsec ($\sim 3-8$ kpc) annular region. The corresponding X-ray luminosity at J1429 redshift is $4\times 10^{45}$ erg s$^{-1}$. If confirmed, this emission could be related to either a large-scale X-ray jet, or a separate X-ray source.

Jets of energetic particles, as seen in FR type-I and FR type-II sources, ejected from the center of Radio-Loud AGN affect the sources surrounding intracluster medium/intergalactic medium. Placing constraints on the age of such sources is important in order to measure the jet powers and determine the effects on feedback. To evaluate the age of these sources using spectral age models, we require high-resolution multi-wavelength data. The new sensitive and high-resolution MIGHTEE survey of the XMM-LSS field along with data from the Low Frequency Array (LOFAR) and the Giant Metrewave Radio Telescope (GMRT) provide data taken at different frequencies with similar resolution, which enables us to determine the spectral age distribution for radio loud AGN in the survey field. In this study we present a sample of 28 radio galaxies with their best fitting spectral age distribution analyzed using the Jaffe-Perola (JP) model on a pixel-by-pixel basis. Fits are generally good and objects in our sample show maximum ages within the range of 2.8 Myr to 115 Myr with a median of 8.71 Myr. High-resolution maps over a range of frequencies are required to observe detailed age distributions for small sources and high-sensitivity maps will be needed in order to observe fainter extended emission. We do not observe any correlation between the total physical size of the sources and their age and we speculate both dynamical models and the approach to spectral age analysis may need some modification to account for our observations.

Hot DA white dwarfs have fully radiative pure hydrogen atmospheres that are the least complicated to model. Pulsationally stable, they are fully characterized by their effective temperature Teff, and surface gravity log g, which can be deduced from their optical spectra and used in model atmospheres to predict their spectral energy distribution (SED). Based on this, three bright DAWDs have defined the spectrophotometric flux scale of the CALSPEC system of HST. In this paper we add 32 new fainter (16.5 < V < 19.5) DAWDs spread over the whole sky and within the dynamic range of large telescopes. Using ground based spectra and panchromatic photometry with HST/WFC3, a new hierarchical analysis process demonstrates consistency between model and observed fluxes above the terrestrial atmosphere to < 0.004 mag rms from 2700 {\AA} to 7750 {\AA} and to 0.008 mag rms at 1.6{\mu}m for the total set of 35 DAWDs. These DAWDs are thus established as spectrophotometric standards with unprecedented accuracy from the near ultraviolet to the near-infrared, suitable for both ground and space based observatories. They are embedded in existing surveys like SDSS, PanSTARRS and GAIA, and will be naturally included in the LSST survey by Rubin Observatory. With additional data and analysis to extend the validity of their SEDs further into the IR, these spectrophotometric standard stars could be used for JWST, as well as for the Roman and Euclid observatories.

The presence of minor bodies in exoplanetary systems is in most cases inferred through infra-red excesses, with the exception of exocomets. Even if over 35 years have passed since the first detection of exocomets around beta Pic, only ~ 25 systems are known to show evidence of evaporating bodies, and most of them have only been observed in spectroscopy. With the appearance of new high-precision photometric missions designed to search for exoplanets, such as CHEOPS, a new opportunity to detect exocomets is available. Combining data from CHEOPS and TESS we investigate the lightcurve of 5 Vul, an A-type star with detected variability in spectroscopy, to search for non periodic transits that could indicate the presence of dusty cometary tails in the system. While we did not find any evidence of minor bodies, the high precision of the data, along with the combination with previous spectroscopic results and models, allows for an estimation of the sizes and spatial distribution of the exocomets.

While hot ICM in galaxy clusters makes these objects powerful X-ray sources, the cluster's outskirts and overdense gaseous filaments might give rise to much fainter sub-keV emission. Cosmological simulations show a prominent "focusing" effect of rich clusters on the space density of the Warm-Hot Intergalactic Medium (WHIM) filaments up to a distance of $\sim 10\,{\rm Mpc}$ ($\sim$ turnaround radius, $r_{ta}$) and beyond. Here, we use Magneticum simulations to characterize their properties in terms of integrated emission measure for a given temperature and overdensity cut and the level of contamination by the more dense gas. We suggest that the annuli $(\sim 0.5-1)\times \,r_{ta}$ around massive clusters might be the most promising sites for the search of the gas with overdensity $\lesssim 50$. We model spectral signatures of the WHIM in the X-ray band and identify two distinct regimes for the gas at temperatures below and above $\sim 10^6\,{\rm K}$. Using this model, we estimate the sensitivity of X-ray telescopes to the WHIM spectral signatures. We found that the WHIM structures are within reach of future high spectral resolution missions, provided that the low-density gas is not extremely metal-poor. We then consider the Coma cluster observed by SRG/eROSITA during the CalPV phase as an example of a nearby massive object. We found that beyond the central $r\sim 40'$ ($\sim 1100\,{\rm kpc}$) circle, where calibration uncertainties preclude clean separation of the extremely bright cluster emission from a possible softer component, the conservative upper limits are about an order of magnitude larger than the levels expected from simulations.

We present a semi-analytic model for predicting kilonova light curves from the mergers of neutron stars with black holes (NSBH). The model is integrated into the MOSFiT platform, and can generate light curves from input binary properties and nuclear equation-of-state considerations, or incorporate measurements from gravitational wave (GW) detectors to perform multi-messenger parameter estimation. The rapid framework enables the generation of NSBH kilonova distributions from binary populations, light curve predictions from GW data, and statistically meaningful comparisons with an equivalent BNS model in MOSFiT. We investigate a sample of kilonova candidates associated with cosmological short gamma-ray bursts, and demonstrate that they are broadly consistent with being driven by NSBH systems, though most have limited data. We also perform fits to the very well sampled GW170817, and show that the inability of an NSBH merger to produce lanthanide-poor ejecta results in a significant underestimate of the early ($\lesssim 2$ days) optical emission. Our model indicates that NSBH-driven kilonovae may peak up to a week after merger at optical wavelengths for some observer angles. This demonstrates the need for early coverage of emergent kilonovae in cases where the GW signal is either ambiguous or absent; they likely cannot be distinguished from BNS mergers by the light curves alone from $\sim 2$ days after the merger. We also discuss the detectability of our model kilonovae with the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST).

The goal of helioseismology is to provide accurate information about the Sun's interior from the observations of the wave field at its surface. In the last three decades, both global and local helioseismology studies have made significant advances and breakthroughs in solar physics. However, 3-d mapping of the structure and dynamics of sunspots and active regions below the surface has been a challenging task and are among the longest standing and intriguing puzzles of solar physics due to the complexity of the turbulent and dynamic nature of sunspots. Thus the key problems that need to be addressed during the next decade are: (i) Understanding the wave excitation mechanisms in the quiet Sun and magnetic regions, (ii) Characterizing the wave propagation and transformation in strong and inclined magnetic field regions and understanding the magnetic portals in the chromosphere, (iii) Improving helioseismology techniques and investigating the whole life cycle of active regions, from magnetic flux emergence to dissipation, and (iv) Detecting helioseismic signature of the magnetic flux of active regions before it becomes visible on the surface so as to provide warnings several days before the emergence. For a transformative progress on these problems require full disk, simultaneous Doppler and vector magnetic field measurements of the photosphere up to the chromosphere with a spatial resolution of about 2 arc-sec as well as large-scale radiative MHD simulations of the plasma dynamics from the sub-photosphere to the chromosphere.

The present work analyzes various aspects of M31 gamma-ray halo emission in its relation to annihilating dark matter (DM). The main aspect is the predicted effect of asymmetry of the intensity of emission due to inverse Compton scattering (ICS) of a possible population of relativistic electrons and positrons ($e^\pm$) in the galactic halo on starlight photons. This asymmetry is expected to exist around the major galactic axis, and arises due to anisotropy of the interstellar radiation field and the inclination of M31. ICS emission and its asymmetry were modeled by GALPROP code for the trial case of $e^\pm$ generated by annihilating weakly interacting massive particles (WIMPs) with various properties. The asymmetry was obtained to appear at photon energies above $\sim$ 0.1 MeV. Morphological and spectral properties of the asymmetry were studied in detail. Potential observational detection of the asymmetry may allow to infer the leptonic fraction in the emission generation mechanism, thus providing valuable inferences for understanding the nature of M31 gamma-ray halo emission. Specific asymmetry predictions were made for the recently claimed DM interpretation of the outer halo emission. The paper also studied the role of secondary -- ICS and bremsstrahlung -- emissions due to DM annihilation for that interpretation. And, finally, the latter was shown to be in significant tension with the recently derived WIMP constraints by radio data on M31.

We detect the cross-correlation between 2.7 million DESI quasar targets across 14,700 deg$^2$ (180 quasars deg$^{-2}$) and Planck 2018 CMB lensing at $\sim$30$\sigma$. We use the cross-correlation on very large scales to constrain local primordial non-Gaussianity via the scale dependence of quasar bias. The DESI quasar targets lie at an effective redshift of 1.51 and are separated into four imaging regions of varying depth and image quality. We select quasar targets from Legacy Survey DR9 imaging, apply additional flux and photometric redshift cuts to improve the purity and reduce the fraction of unclassified redshifts, and use early DESI spectroscopy of 194,000 quasar targets to determine their redshift distribution and stellar contamination fraction (2.6%). Due to significant excess large-scale power in the quasar autocorrelation, we apply weights to mitigate contamination from imaging systematics such as depth, extinction, and stellar density. We use realistic contaminated mocks to determine the greatest number of systematic modes that we can fit, before we are biased by overfitting and spuriously remove real power. We find that linear regression with one to seven imaging templates removed per region accurately recovers the input cross-power, $f_{\textrm{NL}}$ and linear bias. As in previous analyses, our $f_{\textrm{NL}}$ constraint depends on the linear primordial non-Gaussianity bias parameter, $b_{\phi} = 2(b - p)\delta_c$ assuming universality of the halo mass function. We measure $f_{\textrm{NL}} = -26^{+45}_{-40}$ with $p=1.6$ $(f_{\textrm{NL}} = -18^{+29}_{-27}$ with $p=1.0$), and find that this result is robust under several systematics tests. Future spectroscopic quasar cross-correlations with Planck lensing lensing can tighten the $f_{\textrm{NL}}$ constraints by a factor of 2 if they can remove the excess power on large scales in the quasar auto power spectrum.

Supernovae (SNe) inject $\sim 10^{51}$ erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the FLASH code to calculate the cooling radiation from shock-heated gas using collisional excitation data from MAPPINGS V. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ($[$O III$]$, $[$N II$]$, $[$S II$]$, H${\alpha}$, H${\beta}$) are usually best observable. Our shock detection scheme shows that [S II] and [N II] emissions arise from the thin shell surrounding the SNR, while [O III], H$\rm \alpha$, and H$\rm \beta$ originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR's complex structure and its projection onto the plane of the sky because the escaping line luminosity can be reduced by 10 -- 80\% due to absorption along the line-of-sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity.

Seven years after the first direct detection of gravitational waves, from the collision of two black holes, the field of gravitational wave astronomy is firmly established. A first detection of continuous gravitational waves from rapidly-spinning neutron stars could be the field's next big discovery. I review the last twenty years of efforts to detect continuous gravitational waves using the LIGO and Virgo gravitational wave detectors. I summarise the model of a continuous gravitational wave signal, the challenges to finding such signals in noisy data, and the data analysis algorithms that have been developed to address those challenges. I present a quantitative analysis of 291 continuous wave searches from 78 papers, published from 2003 to 2022, and compare their sensitivities and coverage of the signal model parameter space.

When studied at finite temperature, Yang-Mills theories in $3+1$ dimensions display the presence of confinement/deconfinement phase transitions, which are known to be of first order -- the $SU(2)$ gauge theory being the exception. Theoretical as well as phenomenological considerations indicate that it is essential to establish a precise characterisation of these physical systems in proximity of such phase transitions. We present and test a new method to study the critical region of parameter space in non-Abelian quantum field theories on the lattice, based upon the Logarithmic Linear Relaxation (LLR) algorithm. We apply this method to the $SU(3)$ Yang Mills lattice gauge theory, and perform extensive calculations with one fixed choice of lattice size. We identify the critical temperature, and measure interesting physical quantities near the transition. Among them, we determine the free energy of the model in the critical region, exposing for the first time its multi-valued nature with a numerical calculation from first principles, providing this novel evidence in support of a first order phase transition. This study sets the stage for future high precision measurements, by demonstrating the potential of the method.

In this article we investigate the cumulative stochastic gravitational wave spectra as a tool to gain insight on the creation mechanism of primordial black holes. We consider gravitational waves from the production mechanism of primordial black holes and from the gravitational interactions of those primordial black holes among themselves and other astrophysical black holes. We specifically focus on asynchronous bubble nucleation during a first order phase transition as the creation mechanism. We have used two benchmark phase transitions through which the primordial black holes and the primary gravitational wave spectra have been generated. We have considered binary systems and close hyperbolic interactions of primordial black holes with other primordial and astrophysical black holes as the source of the secondary part of the spectra. We have shown that this unique cumulative spectra have features which directly and indirectly depend on the specifics of the production mechanism.

Nonlocal gravity (NLG), a classical extension of Einstein's theory of gravitation, has been studied mainly in linearized form. In particular, nonlinearities have thus far prevented the treatment of cosmological models in NLG. In this essay, we discuss the local limit of NLG and apply this limit to the expanding homogenous and isotropic universe. The theory only allows spatially flat cosmological models; furthermore, de Sitter spacetime is forbidden. The components of the model will have different dynamics with respect to cosmic time as compared to the standard $\Lambda$CDM model; specifically, instead of the cosmological constant, the modified flat model of cosmology involves a dynamic dark energy component in order to account for the accelerated phase of the expansion of the universe.