Papers for Wednesday, Feb 22 2023

We announce FluxCT, a web tool for identifying contaminating flux in Kepler and TESS target pixel files due to secondary visual sources. We demonstrate the usage of this tool and discuss the benefits of this tool over a simple Gaia radius search. FluxCT focuses on clarity and simplicity, where the only input needed from the user is a KIC or TIC ID. By more appropriately accounting for the actual shape of the photometric pixel apertures, FluxCT can produce much more accurate estimates of contaminating flux than simple radial cone searches.

We report the discovery of two apparently isolated stellar remnants that exhibit rotationally modulated magnetic Balmer emission, adding to the emerging DAHe class of white dwarf stars. While the previously discovered members of this class show Zeeman-split triplet emission features corresponding to single magnetic field strengths, these two new objects exhibit significant fluctuations in their apparent magnetic field strengths with variability phase. The Zeeman-split hydrogen emission lines in LP $705{-}64$ broaden from $9.4$ MG to $22.2$ MG over an apparent spin period of $72.629$ minutes. Similarly, WD J$143019.29{-}562358.33$ varies from $5.8$ MG to $8.9$ MG over its apparent $86.394$-minute rotation period. This brings the DAHe class of white dwarfs to at least five objects, all with effective temperatures within $500$ K of $8000$ K and masses ranging from $0.65{-}0.83\,M_{\odot}$.

The origin of Carbon Enhanced Metal-Poor (CEMP-no) stars with low abundances of neutron-capture elements is still unclear. These stars are ubiquitous, found primarily in the Milky Way halo and ultra-faint dwarf galaxies (UFDs). To make a major step forward, we developed a data-calibrated model for B\"ootes I that simultaneously includes all carbon sources: supernovae and asymptotic giant branch (AGB) stars both from first (Pop III) stars, and subsequent normal star formation (Pop II). We demonstrate that each of these sources leave a specific chemical signature in the gas, allowing us to identify the origin of present day CEMP-no stars through their location in the A(C)-[Fe/H] diagram. The CEMP stars with A(C)>6 are predominantly enriched by AGB Pop II stars. We identify a new class of 'moderate CEMP-s' stars with A(C)~ 7 and 0<[Ba/Fe]<+1 , imprinted by winds from AGB stars. True Pop III descendants are predicted to have A(C)<6 and a constant [C/Mg] with [Fe/H], in perfect agreement with observations in B\"ootes I and the Milky Way halo. For the first time we now have a complete picture of the origins of CEMP-no stars which can and will be verified with future observations.

We present an analytic model for measuring the jet core angle ($\theta_c$) and viewing angle ($\theta_{obs}$) of off-axis gamma-ray bursts independently of the jet angular structure outside of the core. We model the images of off-axis jets and using this model we show that $\theta_{obs}$ and $\theta_c$ can be measured using any two of the three following observables: the afterglow light curve, the flux-centroid motion, and the image width. The model is calibrated using 2D relativistic hydrodynamic simulations with a broad range of jet angular structures. We study the systematic errors due to the uncertainty in the jet structure and find that when using the light curve and centroid motion to determine $\theta_{obs}$ and $\theta_c$, our formulae can be accurate to a level of $5-10\%$ and $30\%$, respectively. In light of the Hubble tension, the systematic error in $\cos\theta_{obs}$ in GRBs originating in a binary compact object merger is of special interest. We find that the systematic uncertainty on the measurement of $\cos\theta_{obs}$ due to the unknown jet structure is smaller than $1.5\%$ for well-observed events. A similar error is expected if the microphysical parameters evolve at a level that is not easily detected by the light curve. Our result implies that this type of systematic uncertainty will not prevent measurement of $H_0$ to a level of $2\%$ with a sample of well-observed GW events with resolved afterglow image motion. Applying our model to the light curve and centroid motion observations of GW170817 we find $\theta_{obs}=19.2\pm 2~\deg$ (1$\sigma$) and $\theta_c=1.5-4~\deg$.

We use the large spectroscopic dataset of the MOSFIRE Deep Evolution Field (MOSDEF) survey to investigate some of the key factors responsible for the elevated ionization parameters (U) inferred for high-redshift galaxies, focusing in particular on the role of star-formation-rate surface density (Sigma_SFR). Using a sample of 317 galaxies with spectroscopic redshifts z~1.9-3.7, we construct composite rest-frame optical spectra in bins of Sigma_SFR and infer electron densities, n_e, using the ratio of the [OII] 3727, 3730 doublet. Our analysis suggests a significant (~3 sigma) correlation between n_e and Sigma_SFR. We further find significant correlations between U and Sigma_SFR for composite spectra of a subsample of 113 galaxies, and for a smaller sample of 25 individual galaxies with inferences of U. The increase in n_e -- and possibly also the volume filling factor of dense clumps in HII regions -- with Sigma_SFR appear to be important factors in explaining the relationship between U and Sigma_SFR. Further, the increase in n_e and SFR with redshift at a fixed stellar mass can account for most of the redshift evolution of U. These results suggest that the gas density, which sets n_e and the overall level of star-formation activity, may play a more important role than metallicity evolution in explaining the elevated ionization parameters of high-redshift galaxies.

A small subset of extragalactic double radio sources, termed HYMORS (HYbrid MOrpholgy Radio Sources), is distinguished by a very unusual, hybrid morphology in terms of the Fanaroff-Riley (FR) classification. In HYMORS, one radio lobe appears edge-darkened (FR I), while the other shows a well-defined emission peak near its outer edge (edge-brightened, FR II). Such sources are rare but critical for constraining the mechanism responsible for the FR dichotomy, a widely debated issue in extragalactic astrophysics. Here we highlight the need for caution in assigning FR type, in view of some upcoming observational campaigns to confirm HYMORS among the candidates. To illustrate this we highlight the cases of 3 radio sources which have been perceived to be HYMORS, including the radio galaxy 0500+630 (4C +63.07) which has been claimed to be a good, original example of a HYMORS, with a FR I western lobe and a FR II eastern lobe marked by a prominent terminal hot spot. However, its recent VLASS map at 3 GHz has revealed that the western lobe actually extends much farther out than reported and terminates in a well-defined emission peak. This implies that the source is a regular FR II radio galaxy and not a HYMORS. We also provide a brief perspective of the HYMORS phenomenon and underscore the need to confirm a FR I classification by ruling out additional FR II characteristics, such as an inward lobe-widening and spectral steepening, as well as a lack of prominent radio jet within the lobe.

We report the discovery of an extreme galaxy overdensity at $z = 5.4$ in the GOODS-S field using JWST/NIRCam imaging from JADES and JEMS alongside JWST/NIRCam wide field slitless spectroscopy from FRESCO. We identified potential members of the overdensity using HST+JWST photometry spanning $\lambda = 0.4-5.0\ \mu\mathrm{m}$. These data provide accurate and well-constrained photometric redshifts down to $m \approx 29-30\,\mathrm{mag}$. We subsequently confirmed $N = 96$ potential members of the overdensity using JWST slitless spectroscopy over $\lambda = 3.9-5.0\ \mu\mathrm{m}$ through a targeted line search for $\mathrm{H} \alpha$ around the best-fit photometric redshift. We verified that $N = 53$ galaxies reside in the field at $z = 5.2-5.5$ while $N = 43$ galaxies reside in an overdensity at $z = 5.4$ around $\sim 10-12$ times that of a random volume. Stellar populations for these galaxies were inferred from the photometry and used to construct the star-forming main sequence, where protocluster members appeared more massive and exhibited earlier star formation (and thus older stellar populations) when compared to their field galaxy counterparts. We estimate the total halo mass of this large-scale structure to be $13.0 \lesssim \mathrm{log}_{10} \left( M_{\mathrm{halo}}/M_{\odot} \right) \lesssim 13.5$ using an empirical stellar mass to halo mass relation, although this is likely an underestimate as a result of incompleteness. Our discovery demonstrates the power of JWST at constraining dark matter halo assembly and galaxy formation at very early cosmic times.

Accretion discs properties should deviate from standard theory when magnetic pressure exceeds the thermal pressure. To quantify these deviations, we present a systematic study of the dynamical properties of magnetically arrested discs (MADs), the most magnetized type of accretion disc. Using an artificial cooling function to regulate the gas temperature, we study MADs of three different thermal thicknesses, $h_\mathrm{th}/r=0.3, 0.1$ and $0.03$. We find that the radial structure of the disc is never mostly supported by the magnetic field. In fact, thin MADs are very near Keplerian. However, as discs gets colder, they become more magnetized and the largest deviations from standard theory appear in our thinnest disc with $h_\mathrm{th}/r=0.03$. In this case, the disc is much more extended vertically and much less dense than in standard theory because of vertical support from the turbulent magnetic pressure and wind-driven angular momentum transport that enhances the inflow speed. The thin disc also dissipates a lot of thermal energy outside of $z/r = \pm 0.03$ and a significant fraction of this dissipation happens in mildly relativistic winds. The enhanced dissipation in low-density regions could possibly feed coronae in X-ray binaries (XRBs) and active galactic nuclei (AGN). Wind-driven accretion will also impact the dynamical evolution of accretion discs and could provide a mechanism to explain the rapid evolution of changing-look AGN and the secular evolution of XRBs. Finally, our MAD winds have terminal velocities and mass loss rates in good agreement with the properties of ultra-fast outflows observed in AGN.

General classical arguments on the time evolution of the phase-space density can be used to derive constraints on the mass of particle candidates for the cosmological dark matter (DM). The resulting Tremaine-Gunn limit is extremely useful in constraining particle DM models. In certain models, however, the DM particle mass varies appreciably over time. In this work, we generalize the phase-space limits on possible DM particle masses to these scenarios. We then examine the ensuing cosmological implications on the effective DM equation of state and indirect DM detection.

Variability in the far ultraviolet (FUV) emission produced by stellar activity affects photochemistry and heating in orbiting planetary atmospheres. We present a comprehensive analysis of the FUV variability of GJ 436, a field-age, M2.5V star ($P_\mathrm{rot}\approx44$ d) orbited by a warm, Neptune-size planet ($M \approx 25\ M_\oplus$, $R \approx 4.1\ R_\oplus$, $P_\mathrm{orb}\approx2.6$ d). Observations at three epochs from 2012 to 2018 span nearly a full activity cycle, sample two rotations of the star and two orbital periods of the planet, and reveal a multitude of brief flares. Over 2012-2018, the star's $7.75\pm0.10$ yr activity cycle produced the largest observed variations, $38\pm3$% in the summed flux of major FUV emission lines. In 2018, variability due to rotation was $8\pm2$%. An additional $11\pm1$% scatter at 10 min cadence, treated as white noise in fits, likely has both instrumental and astrophysical origins. Flares increased time-averaged emission by 15% over the 0.88 d of cumulative exposure, peaking as high as 25$\times$ quiescence. We interpret these flare values as lower limits given that flares too weak or too infrequent to have been observed likely exist. GJ 436's flare frequency distribution (FFD) at FUV wavelengths is unusual compared to other field-age M dwarfs, exhibiting a statistically-significant dearth of high energy ($>4\times 10^{28}$ erg) events that we hypothesize to be the result of a magnetic star-planet interaction (SPI) triggering premature flares. If an SPI is present, GJ 436 b's magnetic field strength must be $\lesssim$100 G to explain the statistically insignificant increase in orbit-phased FUV emission.

We present the first results from an imaging and a spectroscopic survey of the optical emission associated with supernova remnant (SNR) G107.5$-$1.5. We discovered optical diffuse and filamentary emission from G107.5$-$1.5 using the 1.5-m and 1-m telescopes. The optical emission from the North East (NE) and North West (NW) regions show the diffuse structure, while the South East (SE) and East (E) regions show filamentary structure. From long-slit spectra, we found [SII]/H$\alpha$ $>$ 0.5 for the SE and E regions, which supports the origin of the emission being from shock-heated gas. The average [SII]/H$\alpha$ ratio is found to be $\sim$0.4 and $\sim$0.3 for the NW and NE regions, respectively, indicating that the optical emission is coming from ionized gas in an HII region. From the ratios of [SII]$\lambda\lambda$ 6716/6731, we estimate an average electron density of $\sim$2400 cm$^{-3}$ for the NW region, which can be attributed to the dense ionized gas. The average [SII]$\lambda\lambda$ 6716/6731 ratios are $\sim$1.25 and $\sim$1.15 for the SE and E regions, respectively, which are indicative of low electron density.

Superfluid dark matter, consisting of self-interacting light particles that thermalize and condense to form a superfluid in galaxies, provides a novel theory that matches the success of the standard $\Lambda$CDM model on cosmological scales while simultaneously offering a rich phenomenology on galactic scales. Within galaxies, the dark matter density profile consists of a nearly homogeneous superfluid core surrounded by an isothermal envelope. In this work we compute the density profile of superfluid dark matter around supermassive black holes at the center of galaxies. We show that, depending on the fluid equation of state, the dark matter profile presents distinct power-law behaviors, which can be used to distinguish it from the standard results for collisionless dark matter.

Massive stars have far-reaching feedback effects that alter the surrounding environment on local, global, and cosmic scales. Spectral analyses of massive stars with adequate stellar-atmosphere models are important to study massive star feedback in detail. We discuss the most recent UV and optical studies of massive metal-poor stars, including those with metallicities ranging from half to one twentieth of solar, connected with large-scale ISM structures in the Magellanic Clouds and the tidal Magellanic Bridge. We present ionizing fluxes from massive stars with low metallicity along with mechanical energy, and we further compare these to the observed energetics in the ISM. The results give hints on the leakage of hot gas and ionizing photons in the Magellanic Clouds. The paper outlines feedback from individual massive stars to population-level collective feedback, the significance of various feedback mechanisms (radiation, wind, supernova), and the influence by the physical conditions of the ISM

Directly imaging temperate rocky planets orbiting nearby, Sun-like stars with a 6-m-class IR/O/UV space telescope, recently dubbed the Habitable Worlds Observatory, is a high priority goal of the Astro2020 Decadal Survey. To prepare for future direct imaging surveys, the list of potential targets should be thoroughly vetted to maximize efficiency and scientific yield. We present an analysis of archival radial velocity data for southern stars from the NASA/NSF Extreme Precision Radial Velocity Working Group's list of high priority target stars for future direct imaging missions (drawn from the HabEx, LUVOIR, and Starshade studies). For each star, we constrain the region of companion mass and period parameter space we are already sensitive to based on the observational baseline, sampling, and precision of the archival RV data. Additionally, for some of the targets we report new estimates of magnetic activity cycle periods, rotation periods, improved orbital parameters for previously known exoplanets, and new candidate planet signals that require further vetting or observations to confirm. Our results show that for many of these stars we are not yet sensitive to even Saturn-mass planets in the habitable zone, let alone smaller planets, highlighting the need for future EPRV vetting efforts before the launch of a direct imaging mission. We present evidence that the candidate temperate super-Earth exoplanet HD 85512 b is most likely due to the star's rotation, and report an RV acceleration for delta Pav which supports the existence of a distant giant planet previously inferred from astrometry.

There is controversy whether dusty starbursts selected at submillimeter wavelengths can trace galaxy overdensities. We perform the first systematic search for protoclusters around a homogeneously selected sample of 12 spectroscopically confirmed submillimeter galaxies (SMGs) at $z\sim1.2-5.3$ in the GOODS-N field. We applied the Poisson Probability Method (PPM) to search for Mpc scale overdensities around these SMGs using three photometric redshift catalogs. We detect galaxy overdensities for 11 out of the 12 SMGs ($92\%\pm8$\%), distributed over eight protoclusters. We confirm three previously discovered protoclusters, and we detect five new ones around the SMGs SMMJ123634 ($z=1.225$), ID.19 ($z=2.047$), SMMJ123607 ($z=2.487$), SMMJ123606 ($z=2.505$), and GN10 ($z=5.303$). A wavelet-based analysis shows that the SMGs live in protocluster cores with a complex morphology (compact, filamentary, or clumpy) and an average size of $\sim(0.4-1)$Mpc. By comparing the PPM results obtained using independently the three redshift catalogs, we possibly witness a transitioning phase at $z\gtrsim4$ for the galaxy populations. While $z\lesssim4$ protoclusters appear to be populated by dusty galaxies, those at highest redshifts $z\sim5$ are detected as overdensities of Lyman$\alpha$ emitters or Lyman break galaxies. We also find a good correlation between the molecular (H$_2$) gas mass of the SMG and the overdensity significance. To explain the overall phenomenology, we suggest that galaxy interactions in dense environments likely triggered the starburst and gas-rich phase of the SMGs. Altogether, we support the scenario that SMGs are excellent tracers of distant protoclusters. Those presented in this work are excellent targets for the {\it James Webb Space Telescope.} Surveys with forthcoming facilities (e.g., {\it Euclid}, LSST) can be tuned to detect even larger samples of distant protoclusters.

Improved observational constraints on the orbital parameters of the low-mass X-ray binary Scorpius~X-1 were recently published in Killestein et al (2023). In the process, errors were corrected in previous orbital ephemerides, which have been used in searches for continuous gravitational waves from Sco~X-1 using data from the Advanced LIGO detectors. We present the results of a re-analysis of LIGO detector data from the third observing run of Advanced LIGO and Advanced Virgo using a model-based cross-correlation search. The corrected region of parameter space, which was not covered by previous searches, was about 1/3 as large as the region searched in the original O3 analysis, reducing the required computing time. We have confirmed that no detectable signal is present over a range of gravitational-wave frequencies from $25\textrm{Hz}$ to $1600\textrm{Hz}$, analogous to the null result of Abbott et al (2022). Our search sensitivity is comparable to that of Abbott et al (2022), who set upper limits corresponding, between $100\textrm{Hz}$ and $200\textrm{Hz}$, to an amplitude $h_0$ of about $10^{-25}$ when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than $4\times 10^{-26}$ assuming the optimal orientation.

Neutron star mergers have been proposed as the main source of heavy $r$-process nucleosynthesis in the Universe. However, the mergers' significant expected delay after binary formation is in tension with observed very early $r$-process enrichment, e.g., in the dwarf galaxy Reticulum II. The LIGO and Virgo gravitational-wave observatories discovered two binary mergers with lighter companion masses ($\sim 2.6$ M$_\odot$) similar to the total mass of many binary neutron star systems in the Galaxy. The progenitor of such mergers could be a neutron star binary orbiting a black hole. Here we show that a significant fraction of neutron star binaries in hierarchical triples merge rapidly ($\gtrsim3\%$ within $\lesssim10$ Myr after neutron star formation) and could explain the observed very early $r$-process enrichment. The neutron star binary can become eccentric via von Zeipel-Kozai-Lidov oscillations, promoting a fast coalescence followed later by a merger of the low-mass black hole with the higher-mass black hole in the system. We show that this scenario is also consistent with an overall binary neutron star merger rate density of $\sim100$ Gpc$^{-3}$yr$^{-1}$ in such triples. Using hydrodynamic simulations we show that highly eccentric neutron star mergers dynamically eject several times more mass than standard mergers, with exceptionally bright kilonovae with an "early blue bump" as unique observational signatures.

Magnetars, neutron stars thought to be with ultra-strong magnetic fields of $10^{14 - 15}$ G, are observed to be much hotter than ordinary pulsars with $\sim 10^{12}$ G, and additional heating sources are required. One possibility is heating by the ambipolar diffusion in the stellar core. This scenario is examined by calculating the models using the relativistic thermal evolutionary code without making the isothermal approximation. The results show that this scenario can be consistent with most of the observed magnetar temperature data.

UV-FIR SED modeling is an effective way to disentangle emission between star formation (SF) and active galactic nuclei (AGN) in galaxies; however, this approach becomes uncertain for composite AGN/SF galaxies that comprise 50-70% of IR-samples. Cosmic X-ray background (XRB) models require a large fraction of obscured AGN to reproduce the observed XRB peak, motivating reliable SED analyses in objects where the AGN may be ``buried" in the galaxy and in the mid-IR to far-IR SED. In this paper, we study a 24$\mu$m-selected ($S_{24}$ > 100$\mu$Jy) sample of 95 galaxies with $0 \% < f_{MIR,AGN} < 100 \%$, 0.4 < z < 2.7, and $10^{11}$L$_{\odot}$ < L$_{IR}$ < $10^{13}$L$_{\odot}$. We test the performance of AGN models ranging in torus optical depth via SED fitting, comparing results with Spitzer MIR spectroscopy and X-ray observations. Best-fit torus optical depth can shed light on whether these galaxies host a luminous obscured AGN population. We find that permitting a broader AGN SED parameter space results in improved fit quality with higher optical depths, higher FIR AGN contributions, and higher $L_{Bol}$, impacting the bright-end of the $L_{Bol}$ luminosity function. Our results suggest there may be a population of dust-obscured composites that are bolometrically significant but have their AGN mostly hidden in the mid-IR SED. If so, literature applications of SED fitting that often simplify AGN models or omit optically thick tori may largely underestimate AGN contribution from composite sources, as these sources are both numerous and have solutions sensitive to the assumed range of AGN models.

Kepler Space Telescope and Transiting Exoplanet Survey Satellite unveiled that Sun-like stars frequently host exoplanets. These exoplanets are subject to fluxes of ionizing radiation in the form of X-ray and extreme-ultraviolet (EUV) radiation that may cause changes in their atmospheric dynamics and chemistry. While X-ray fluxes can be observed directly, EUV fluxes cannot be observed because of severe interstellar medium absorption. Here, we present a new empirical method to estimate the whole stellar XUV (X-ray plus EUV) and FUV spectra as a function of total unsigned magnetic fluxes of stars. The response of the solar XUV and FUV spectrum (0.1-180 nm) to the solar total unsigned magnetic flux is investigated by using the long-term Sun-as-a-star dataset over 10 yrs, and the power-law relation is obtained for each wavelength with a spectral resolution of 0.1-1 nm. We applied the scaling relations to active young Sun-like stars (G-dwarfs), EK Dra (G1.5V), $\pi^1$ Uma (G1.5V) and $\kappa^1$ Ceti (G5V), and found that the observed spectra (except for the unobservable longward EUV wavelength) are roughly consistent with the extension of the derived power-law relations with errors of an order of magnitude. This suggests that our model is a valuable method to derive the XUV/FUV fluxes of Sun-like stars including the EUV band mostly absorbed at wavelengths longward of 36 nm. We also discuss differences between the solar extensions and stellar observations at the wavelength in the 2-30 nm band and concluded that simultaneous observations of magnetic and XUV/FUV fluxes are necessary for further validations.

We develop an accretion disc (AD) fitting method, utilising thin and slim disc models and Bayesian inference with the Markov-Chain Monte-Carlo approach, testing it on the most luminous known quasar, SMSS J215728.21-360215.1, at redshift $z=4.692$. With a spectral energy distribution constructed from near-infrared spectra and broadband photometry, the AD models find a black hole mass of $\log(M_{\rm{AD}}/M_{\odot}) = 10.31^{+0.17}_{-0.14}$ with an anisotropy-corrected bolometric luminosity of $\log{(L_{\rm{bol}}/\rm{erg\,s^{-1}})} = 47.87 \pm 0.10$, and derive an Eddington ratio of $0.29^{+0.11}_{-0.10}$ as well as a radiative efficiency of $0.09^{+0.05}_{-0.03}$. Using the near-infrared spectra, we estimate the single-epoch virial black hole mass estimate to be $\log(M_{\rm{SE}}/M_{\odot}) = 10.33 \pm 0.08$, with a monochromatic luminosity at 3000\AA\ of $\log{(L(\rm{3000\text{\AA}})/\rm{erg\,s^{-1}})} = 47.66 \pm 0.01$. As an independent approach, AD fitting has the potential to complement the single-epoch virial mass method in obtaining stronger constraints on properties of massive quasar black holes across a wide range of redshifts.

How does the appearance of a strongly lensed system change if a gravitational wave is produced by the lens? In this work we address this question by considering a supermassive black hole binary at the center of the lens emitting gravitational waves propagating either colinearly or orthogonally to the line of sight. Specializing to an Einstein ring configuration (where the source, the lens and the observer are aligned), we show that the gravitational wave induces changes on the ring's angular size and on the optical path of photons. The changes are the same for a given pair of antipodal points on the ring, but maximally different for any pair separated by $90^{\circ}$. For realistic lenses and binaries, we find that the change in the angular size of the Einstein ring is dozens of orders of magnitude smaller than the precision of current experiments. On the other hand, the difference in the optical path induced on a photon by a gravitational wave propagating \textit{orthogonally} to the line of sight triggers, at peak strain, time delays in the range $\sim 0.01 - 1$ seconds, making the chance of their detection (and thus the use of Einstein rings as gravitational wave detectors) less hopeless.

Ba-enhanced stars are interesting probes of stellar astrophysics and Galactic formation history. In this work, we investigate the chemistry and kinematics for a large sample of Ba-enhanced ([Ba/Fe]$>$1.0) dwarf and subgiant stars with $5000 < T_{\rm eff }< 6700$\,K from LAMOST. We find that both stellar internal evolution process and external mass exchange due to binary evolution are responsible for the origins of the Ba-enhancement of our sample stars. About one third of them exhibit C and N enhancement and ultraviolet brightness excess, indicating they are products of binary evolution. The remaining Ba-enhanced stars with normal C and N abundances are mostly warm stars with $T_{\rm eff} > 6000$\,K. They are likely consequences of stellar internal elemental transport processes, but they show very different elemental patterns to the hotter Am/Fm stars. Our results reveal a substantially lack of high-[$\alpha$/Fe] Ba-enhanced stars in the [Fe/H]--[$\alpha$/Fe] plane, which we dub as a {\em high-$\alpha$ desert}. We suggest it is due to a lower efficiency for producing Ba-enhanced stars by low-mass AGB progenitors in binary systems. Our results call for detailed modellings for these Ba-enhanced stellar peculiars, in the context of both stellar internal elemental transport and external mass accretion.

The Chang'E 5 (CE-5) samples represent the youngest mare basalt ever known and provide an access into the late lunar evolution. Recent studies have revealed that CE-5 basalts are the most evolved lunar basalt, yet controversy remains over the nature of their mantle sources. Here we combine Fe and Mg isotope analyses with a comprehensive study of petrology and mineralogy on two CE-5 basalt clasts. These two clasts have a very low Mg# (~29) and show similar Mg isotope compositions with Apollo low-Ti mare basalts as well as intermediate TiO2 and Fe isotope compositions between low-Ti and high-Ti mare basalts. Fractional crystallization or evaporation during impact cannot produce such geochemical signatures which otherwise indicate a hybrid mantle source that incorporates both early- and late-stage lunar magma ocean (LMO) cumulates. Such a hybrid mantle source would be also compatible with the KREEP-like REE pattern of CE-5 basalts. Overall, our new Fe-Mg isotope data highlight the role of late LMO cumulate for the generation of young lunar volcanism.

This study examines interacting quintessence dark energy models and their observational constraints for a general parameterization of the quintessence potential, which encompasses a broad range of popular potentials. Four different forms of interactions are considered. The analysis is done by expressing the system as a set of autonomous equations for each interaction. The Bayesian Model Comparison has been used to compare these models with the standard Lambda Cold Dark Matter ($\Lambda$CDM) model. Our analysis shows that, compared to all of the considered interacting models, the $\Lambda$CDM model is still preferred.

The CARMENES instrument was conceived to deliver high-accuracy radial velocity (RV) measurements with long-term stability to search for temperate rocky planets around a sample of nearby cool stars. The broad wavelength coverage was designed to provide a range of stellar activity indicators to assess the nature of potential RV signals and to provide valuable spectral information to help characterise the stellar targets. The CARMENES Data Release 1 (DR1) makes public all observations obtained during the CARMENES guaranteed time observations, which ran from 2016 to 2020 and collected 19,633 spectra for a sample of 362 targets. The CARMENES survey target selection was aimed at minimising biases, and about 70% of all known M dwarfs within 10 pc and accessible from Calar Alto were included. The data were pipeline-processed, and high-level data products, including 18,642 precise RVs for 345 targets, were derived. Time series data of spectroscopic activity indicators were also obtained. We discuss the characteristics of the CARMENES data, the statistical properties of the stellar sample, and the spectroscopic measurements. We show examples of the use of CARMENES data and provide a contextual view of the exoplanet population revealed by the survey, including 33 new planets, 17 re-analysed planets, and 26 confirmed planets from transiting candidate follow-up. A subsample of 238 targets was used to derive updated planet occurrence rates, yielding an overall average of 1.44+/-0.20 planets with 1 M_Earth < M sin i < 1000 M_Earth and 1 d < P_orb < 1000 d per star, and indicating that nearly every M dwarf hosts at least one planet. CARMENES data have proven very useful for identifying and measuring planetary companions as well as for additional applications, such as the determination of stellar properties, the characterisation of stellar activity, and the study of exoplanet atmospheres.

Quasi-periodic Oscillations (QPOs) from magnetars bursts are rare, some high-confidence detections are mostly from magnetar Gaint Flares. Following a systematically research of bursts data from Fermi/GBM, here we report a 4.98$\sigma$ ($p$-value is 6.51e-7) detection of high frequency QPO $\sim$ 80 Hz in SGR 150228213. This burst has initially attributed to 4U 0142+61 by Fermi/GBM, however, it shows a different spectral features with other magnetar bursts, including burts from 4U 0142+61. More interestingly, we find the location of SGR 150228213 almost coincides with the periodic Fast Radio Burst source FRB 180916. Finally, we discuss the possible physical origins of this burst.

We study environmental quenching in the Eagle}/C-Eagle cosmological hydrodynamic simulations over the last 11 Gyr (i.e. $z=0-2$). The simulations are compared with observations from the SAMI Galaxy Survey at $z=0$. We focus on satellite galaxies in galaxy groups and clusters ($10^{12}\,\rm M_{\odot}$ $\lesssim$ $M_{200}$ < $3 \times 10^{15}\, \rm M_{\odot}$). A star-formation concentration index [$C$-index $= \log_{10}(r_\mathrm{50,SFR} / r_\mathrm{50,rband})$] is defined, which measures how concentrated star formation is relative to the stellar distribution. Both Eagle/C-Eagle and SAMI show a higher fraction of galaxies with low $C$-index in denser environments at $z=0-0.5$. Low $C$-index galaxies are found below the SFR-$M_{\star}$ main sequence (MS), and display a declining specific star formation rate (sSFR) with increasing radii, consistent with ``outside-in'' environmental quenching. Additionally, we show that $C$-index can be used as a proxy for how long galaxies have been satellites. These trends become weaker at increasing redshift and are absent by $z=1-2$. We define a quenching timescale $t_{\rm quench}$ as how long it takes satellites to transition from the MS to the quenched population. We find that simulated galaxies experiencing ``outside-in'' environmental quenching at low redshift ($z=0\sim0.5$) have a long quenching timescale (median $t_{\rm quench}$ > 2 Gyr). The simulated galaxies at higher redshift ($z=0.7\sim2$) experience faster quenching (median $t_{\rm quench}$ < 2Gyr). At $z\gtrsim 1-2$ galaxies undergoing environmental quenching have decreased sSFR across the entire galaxy with no ``outside-in'' quenching signatures and a narrow range of $C$-index, showing that on average environmental quenching acts differently than at $z\lesssim 1$.

The European Space Agency's Rosetta mission to comet 67P/Churyumov-Gerasimenko has offered scientists the opportunity to study a comet in unprecedented detail. Four instruments of the Rosetta orbiter, namely, the Micro-Imaging Dust Analysis System (MIDAS), the Grain Impact Analyzer and Dust Accumulator (GIADA), the COmetary Secondary Ion Mass Analyser (COSIMA), and the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) have provided information on cometary dust particles. Cross-instrument comparisons are crucial to characterize cometary dust particles beyond the capabilities of individual sensors, as they are sensitive to different dust components. We present the first comparison between detections of the ROSINA COmet Pressure Sensor (COPS) and GIADA. These two instruments are complementary as the former is sensitive solely to volatiles of icy particles, while the latter measured the dust particle as a whole, including refractories and condensed (semi)volatiles. Our goal is to correlate the particles detected by COPS and GIADA and to assess whether they belong to a common group. We statistically analyzed the in situ data of COPS and GIADA by calculating Pearson correlation coefficients. Among the several types of particles detected by GIADA, we find that COPS particles are significantly correlated solely with GIADA fluffy agglomerates (Pearson correlation coefficient of 0.55 and p-value of $4.6\cdot 10^{-3}$). This suggests that fluffy particles are composed of both refractories and volatiles. COPS volatile volumes, which may be represented by equivalent spheres with a diameter in the range between 0.06 $\mu$m and 0.8 $\mu$m, are similar to the sizes of the fractal particle's subunits identified by MIDAS (i.e., 0.05-0.18 $\mu$m).

We report the first measurement of polarized thermal dust emission toward the entire early and massive Infrared Dark Cloud G11.11$-$0.12 taken by the polarimeter SOFIA/HAWC+ at 214 $\mu m$ wavelength. Magnetic fields (B-fields) obtained from the polarized emission tend to be perpendicular to the filament's spine. We produce a map of B-field strengths for the center region of the filament. The strengths vary in the range of 100-600 $\mu\rm{G}$ and are strongest along the filament's spine. The central region is sub-Alfv\'enic and mostly sub-critical meaning that B-fields dominate over turbulence and are strong enough to resist gravitational collapse. The alignment and properties of dust grains are studied in the filament using the RAdiative Torque (RAT) theory. We find the decrease of polarization degree $P$ with emission intensity $I$, i.e., depolarization effect, of the form $P\propto I^{-\alpha}$ with $\alpha\sim$0.8-0.9, implying a significant loss of grain alignment in the filament's spine. The depolarization can be explained by the decrease in RAT alignment efficiency toward the denser regions with lower dust temperature, and cannot be explained by the B-field tangling. We study the effect of the enhanced magnetic relaxation by embedded iron inclusions on RAT alignment and find that the high polarization fraction $P\sim$20-30\% in the outer layer of the filament is potential evidence for the enhanced RAT alignment by magnetic relaxation. This is the first time this effect is evaluated in a filament. Based on the polarization fraction and RAT alignment theory, we find evidence for grain growth in the filament.

Evidence for the presence of ion cyclotron waves, driven by turbulence, at the boundaries of the current sheet is reported in this paper. By exploiting the full potential of the joint observations performed by Parker Solar Probe and the Metis coronagraph on board Solar Orbiter, local measurements of the solar wind can be linked with the large-scale structures of the solar corona. The results suggest that the dynamics of the current sheet layers generates turbulence, which in turn creates a sufficiently strong temperature anisotropy to make the solar-wind plasma unstable to anisotropy-driven instabilities such as the Alfv\'en ion-cyclotron, mirror-mode, and firehose instabilities. The study of the polarization state of high-frequency magnetic fluctuations reveals that ion cyclotron waves are indeed present along the current sheet, thus linking the magnetic topology of the remotely imaged coronal source regions with the wave bursts observed in situ. The present results may allow improvement of state-of-the-art models based on the ion cyclotron mechanism, providing new insights into the processes involved in coronal heating.

CCDs characterization is the preliminary step to perform before the CCD can be properly used at the telescope. Most of the scientific instrumentation at the Italian National Telescope Galileo use CCDs as detectors. In particular the optical imager (OIG) and the high resolution spectrograph (SARG) use a mosaic of two 2k X 4k CCD manufactured by EEV (EEV 4280). The technical characteristics of the EEV4280 can be found in Cosentino et al (these proceedings).

AGILE is a space mission launched in 2007 to study X-ray and gamma-ray astronomy. The AGILE team developed real-time analysis pipelines to detect transient phenomena such as Gamma-Ray Bursts (GRBs) and to react to external science alerts received by other facilities. The AGILE anti-coincidence system (ACS) comprises five panels (four lateral and one on the top) that surround the AGILE detectors to reject background charged particles. It can also detect hard X-ray photons in the energy range 50 - 200 KeV. The acquisition of the ACS data produces a time series for each panel. These time series can be merged in a single multivariate time series (MTS). We present in this work a new Deep Learning model for GRBs detection in the MTSs, generated by the ACS, using an anomaly detection technique. The model is implemented with a Deep Convolutional Neural Network autoencoder architecture. We trained the model with an unsupervised learning algorithm using a dataset of MTSs randomly extracted from the AGILE ACS data. The reconstruction error of the autoencoder is used as the anomaly score to classify the MTS. If the anomaly score is higher than a predefined threshold, the MTS is flagged as a GRB. The trained model is evaluated using a list of MTSs containing GRBs. The tests confirmed the model's ability to detect transient events, providing a new promising technique to identify GRBs in the ACS data that can be implemented in the AGILE real-time analysis pipeline.

In 2021 BL Lacertae underwent an extraordinary activity phase, which was intensively followed by the Whole Earth Blazar Telescope (WEBT) Collaboration. We present the WEBT optical data in the BVRI bands acquired at 36 observatories around the world. In mid 2021 the source showed its historical maximum, with R = 11.14. The light curves display many episodes of intraday variability, whose amplitude increases with source brightness, in agreement with a geometrical interpretation of the long-term flux behaviour. This is also supported by the long-term spectral variability, with an almost achromatic trend with brightness. In contrast, short-term variations are found to be strongly chromatic and are ascribed to energetic processes in the jet. We also analyse the optical polarimetric behaviour, finding evidence of a strong correlation between the intrinsic fast variations in flux density and those in polarisation degree, with a time delay of about 13 h. This suggests a common physical origin. The overall behaviour of the source can be interpreted as the result of two mechanisms: variability on time scales greater than several days is likely produced by orientation effects, while either shock waves propagating in the jet, or magnetic reconnection, possibly induced by kink instabilities in the jet, can explain variability on shorter time scales. The latter scenario could also account for the appearance of quasi-periodic oscillations, with periods from a few days to a few hours, during outbursts, when the jet is more closely aligned with our line of sight and the time scales are shortened by relativistic effects.

AGILE is an ASI space mission launched in 2007 to study X-ray and gamma-ray phenomena in the energy range from $\sim20$ keV to $\sim10$ GeV. The AGILE Team developed a real-time analysis pipeline for the fast detection of transient sources, and the follow-up of external science alerts received through networks such as the General Coordinates Network. We developed a new Deep Learning method for detecting and localizing Gamma-Ray Bursts (GRB) in the AGILE/GRID sky maps. We trained the model using sky maps with GRBs simulated in a radius of 20 degrees from the center of the map, which is larger than 99.5 \% of the error region present in the GRBWeb catalog. We also plan to apply this method to search for counterparts of gravitational wave events, which typically have a wider localization error region. The method comprises two Deep Learning models implemented with two Convolutional Neural Networks. The first model detects and filters sky maps containing a GRB, while the second model localizes its position. We trained and tested the models using simulated data. The detection model achieves an accuracy of 95.7 \%, and the localization model has a mean error lower than 0.8 degrees. We configured a Docker container with all the required software for data simulation and deployed it using the Amazon Web Service to calculate the p-value distribution under different conditions. With the p-value distribution, we can calculate the statistical significance of a detection.

In May 1997 a consistent part of the services and structures committed to the industry had already been released to the commissioning group. The telescope itself was, with the exception of the Nasmyth derotators, motors and all the optics groups, basically ready in its mechanical parts to accept the integration of all services and control equipment. Also the verification of the cabling (interlocks, data-nets, power and controls) already mounted was started in the same period. Starting from June 1998 (telescope first-light date) the telescope went gradually in use, several nights per week, in order to test and tune the tracking and pointing system, the optics and the first derotator system (Nasmyth A station). At the end of the commissioning period and with the first scientific instruments mounted (April 1999) also the first routinely observations started. In this moment the telescope is doing astronomy 80% of time and the complete first-light instrumentation is mounted.

SARG is the high resolution spectrograph of TNG. It has been in operation since late spring 2000. SARG is a cross dispersed echelle spectrograph; it offers both single object and long slit (up to 26 arcsec) observing modes covering a spectral range from {\lambda}=0.37 up to 1 {\mu}m, with resolution ranging from R=29,000 up to R=164,000. Cross dispersion is provided by means of a selection of four grisms; interference filters may be used for the long slit mode (up to 26 arcsec). A dioptric camera images the cross dispersed spectra onto a mosaic of two 2048x4096 EEV CCDs (pixel size: 13.5 {\mu}m) allowing complete spectral coverage at all resolving power for {\AA} <0.8 {\mu}m. An iodine-absorbing cell allows to obtain high precision radial velocities. A Distributed Active Temperature Control System (DATCS) maintains constant the temperature of all spectrograph components at a preset value. Early results show that SARG works according to original specifications in terms of wavelength coverage, efficiency (measured peak efficiency is about 13%), resolution (maximum resolution R~164,000 using a 0.3 arcsec slit, R~144,000 using an image slicer), and stability (preliminary estimates of the radial velocity accuracy is ~5 m/s using the iodine cell and ~ 150 m/s without the cell).

We investigate the fraction of baryon mass in intergalactic medium ($f_\mathrm{IGM}$), using 18 well-localized FRBs in the redshift range $z\in (0.0039,0.66)$. We construct a five-parameter Bayesian inference model, with the probability distributions of dispersion measures (DM) of IGM and host galaxy properly taken into account. To check the possible redshift evolution, we parameterize $f_\mathrm{IGM}$ as a mildly evolving function of redshift, $f_\mathrm{IGM}=f_\mathrm{IGM,0}[1+\alpha z/(1+z)]$. By simultaneously constraining five parameters, we get $f_\mathrm{IGM,0} = 0.92^{+0.06}_{-0.12}$ and $\alpha = 0.49^{+0.59}_{-0.47}$, and the median value of DM of host galaxy is $\exp(\mu)=72.49^{+33.31}_{-25.62}~{\rm pc ~ cm ^ {-3}}$. By fixing two parameters which can be constrained independently with other observations, we obtain $\alpha =0.11^{+0.24}_{-0.27}$ in the three-parameter fit, which is consistent with zero within $1\sigma$ uncertainty. Monte Carlo simulations show that even 300 FRBs are not enough to tightly constrain five parameters simultaneously. This is mainly caused by the correlation between parameters. Only if two parameters are fixed, 100 FRBs are necessary to achieve unbiased constraints on the remaining parameters.

The spectral shape of the gamma-ray emission observed for dynamically old supernova remnants that interact with molecular clouds triggered an exciting scenario of adiabatic compression and farther re-acceleration of Galactic cosmic rays (GCRs) in radiative shells of the remnants, which was extensively discussed and applied to various sources over recent years. Indeed, the observed gamma-ray spectrum from a number of remnants strongly resembles the expected spectrum of the gamma-ray emission from the compressed population of Galactic cosmic rays. In the following we discuss the feasibility of this scenario and show that it is very unlikely that compressed GCRs could produce sufficient amount of gamma-rays and that the observed spectral shape is putting strong limits on the allowed compression factors. Further, absence of curvature in featureless power-law spectra of evolved supernova remnants at radio wavelengths is strongly disfavoring the compression scenario for electrons and hence for hadrons. Our calculations show that the contribution of compressed electrons to the observed radio-flux could reach at most ~10%.

Recently, it has been proposed that the magnetic-field-induced transition (MIT) in Fe X can be used to measure coronal magnetic field strengths. Several techniques, the direct line ratio technique and the weak and strong magnetic field techniques, are developed to apply the MIT theory to spectroscopic observations taken by EUV Imaging Spectrometer (EIS) onboard Hinode. However, the suitability of coronal magnetic field measurements based on the weak and strong magnetic field techniques has not been evaluated. Besides, temperature diagnostics is also important for measuring coronal magnetic field based on the MIT theory, but how to determine the accurate formation temperature of the Fe X lines from EIS observations still needs investigation. In this study, we synthesized emissions of several spectral lines from a 3D radiation magnetohydrodynamic model of a solar active region, and then derived magnetic field strengths using different methods. We first compared the magnetic field strengths derived from the weak and strong magnetic field techniques to the values in the model. Our study suggests that both weak and strong magnetic field techniques underestimate the coronal magnetic field strength. Then we developed two methods to calculate the formation temperature of the Fe X lines. One is based on differential emission measure analyses, and the other is deriving temperature from the Fe IX and Fe XI line pairs. However, neither of the two methods can provide temperature determination for accurate coronal magnetic field measurements as those derived from the Fe X 174/175 and 184/345 {\AA} line ratios. More efforts are still needed for accurate coronal magnetic field measurements using EIS observations.

The study of ultraluminous X-ray sources (ULXs) has changed dramatically over the last decade. In this review we first describe the most important observations of ULXs in various wavebands, and across multiple scales in space and time. We discuss recent progress and current unanswered questions. We consider the range of current theories of ULX properties in the light of this observational progress. Applying these models to neutron-star ULXs offers particularly stringent tests, as this is the unique case where the mass of the accretor is effectively fixed.

Exploiting broad- and narrow-band images of the \textit{Hubble Space Telescope} from near-UV to I-band restframe, we study the star-forming clumps of six galaxies of the GASP sample undergoing strong ram-pressure stripping (RPS). Clumps are detected in H$\alpha$ and near-UV, tracing star formation on different timescales. We consider clumps located in galaxy disks, in the stripped tails and those formed in stripped gas but still close to the disk, called extraplanar. We detect 2406 H$\alpha$-selected clumps (1708 in disks, 375 in extraplanar regions, and 323 in tails) and 3750 UV-selected clumps (2026 disk clumps, 825 extraplanar clumps and 899 tail clumps). Only $\sim15\%$ of star-forming clumps are spatially resolved, meaning that most are smaller than $\sim 140$ pc. We study the luminosity and size distribution functions (LDFs and SDFs, respectively) and the luminosity-size relation. The average LDF slope is $1.79\pm 0.09$, while the average SDF slope is $3.1\pm 0.5$. Results suggest the star formation to be turbulence driven and scale-free, as in main-sequence galaxies. All the clumps, whether they are in the disks or in the tails, have an enhanced H$\alpha$ luminosity at a given size, compared to the clumps in main-sequence galaxies. Indeed, their H$\alpha$ luminosity is closer to that of clumps in starburst galaxies, indicating that ram pressure is able to enhance the luminosity. No striking differences are found among disk and tail clumps, suggesting that the different environments in which they are embedded play a minor role in influencing the star formation.

Dark sirens, i.e., gravitational-wave (GW) sources without electromagnetic counterparts, are new probes of the expansion of the universe. The efficacy of this method relies on correctly localizing the host galaxies. However, recent theoretical studies have shown that astrophysical environments could mislead the spatial localization by distorting the GW signals. It is unclear whether and to what degree the incorrect spatial localizations of dark sirens would impair the accuracy of the measurement of the cosmological parameters. To address this issue, we consider the future observations of dark sirens using the Cosmic Explore and the Einstein Telescope, and we design a Bayesian framework to access the precision of measuring the Hubble-Lema\^itre constant $H_0$. Interestingly, we find that the precision is not compromised when the number of well-localized dark sirens is significantly below $300$, even in the extreme scenario that all the dark sirens are localized incorrectly. As the number exceeds $300$, the incorrect spatial localizations start to produce statistically noticeable effects, such as a slow convergence of the posterior distribution of $H_0$. We propose several tests that can be used in future observations to verify the spatial localizations of dark sirens. Simulations of these tests suggest that incorrect spatial localizations will dominate a systematic error of $H_0$ if as much as $10\%$ of a sample of $300$ well-localized dark sirens are affected by their environments. Our results have important implications for the long-term goal of measuring $H_0$ to a precision of $<1\%$ using dark sirens.

We consider dynamical environments of (486958) Arrokoth, focusing on both their present state and their long-term evolution, starting from the KBO's formation. Both analytical (based on an upgraded Kepler-map formalism) and numerical (based on massive simulations and construction of stability diagrams in the 3D setting of the problem) approaches to the problem are used. The debris removal is due to either absorption by the KBO or by leaving the Hill sphere; the interplay of these processes is considered. The clearing mechanisms are explored, and the debris removal timescales are estimated. We assess survival opportunities for any debris orbiting around Arrokoth. The generic chaotization of Arrokoth's circumbinary debris disk's inner zone and generic cloudization of the disk's periphery, which is shown to be essential in the general 3D case, naturally explains the current absence of any debris in its vicinities.

The coronal magnetic field, despite its overwhelming importance to the physics and dynamics of the corona, has only rarely been measured. Here, the electron density maps derived from images acquired during the total solar eclipse of August 21st, 2017 are employed to demonstrate a new technique to measure the coronal magnetic fields. The strength of the coronal magnetic fields is derived with a semiempirical formula relating the plasma magnetic energy density to the gravitational potential energy. The resulting values are compared with those provided by more advanced coronal field reconstruction methods based on MHD simulations of the whole corona starting from photospheric field measurements, finding a very good agreement. Other parameters such as the plasma-$\beta$ and Alfv\'en velocity are also derived and compared with those of MHD simulations. Moreover, the plane-of-sky (POS) orientation of the coronal magnetic fields is derived from the observed inclination of the coronal features in the filtered images, also finding a close agreement with magnetic field reconstructions. Hence, this work demonstrates for the first time that the 2D distribution of coronal electron densities measured during total solar eclipses is sufficient to provide the coronal magnetic field strengths and inclinations projected on the POS. These are among the main missing pieces of information that limited so far our understanding of physical phenomena going on in the solar corona.

The distribution of the spin frequencies of neutron stars in low-mass X-ray binaries exhibits a cut-off at 730 Hz, below the break-up frequency (mass-shedding limit) of neutron stars. The absence of the sub-millisecond pulsars presents a problem, given that these systems are older than the spin-up timescale. We confront models of disc-magnetosphere interaction near torque equilibrium balanced by the torque due to gravitational wave emission. We note that field lines penetrating the disc beyond the inner radius reduce the maximum rotation frequency of the star, a result well-known since the seminal work of Ghosh & Lamb. We show that the polar cap area corresponds to about half the neutron star surface area at the cut-off frequency if the inner radius is slightly smaller than the corotation radius. We then include the change in the moment of inertia of the star due to the accretion of mass and find that this effect further reduces the maximum rotation frequency of the star. Finally, we include the torque due to gravitational wave emission and calculate its contribution to the torque equilibrium. Our results suggest that all three processes are significant at the cut-off frequency, and all of them must be considered in addressing the absence of millisecond pulsars.

TOI-1055 is a Sun-like star known to host a transiting Neptune-sized planet on a 17.5-day orbit (TOI-1055 b). Radial velocity (RV) analyses carried out by two independent groups using almost the same set of HARPS spectra provide planetary masses that differ by $\sim$2$\sigma$. Our aim was to solve the inconsistency in the published planetary masses by significantly extending the set of HARPS RV measurements and employing a new analysis tool able to account and correct for stellar activity. Our further aim was to improve the precision on the planetary radius by observing two transits of the planet with the CHEOPS space telescope. We fitted a skew normal function to each cross correlation function extracted from the HARPS spectra to obtain RV measurements and hyperparameters to be used for the detrending. We evaluated the correlation changes of the hyperparameters along the RV time series using the breakpoint technique, further performing a joint photometric and RV analysis using a Markov chain Monte Carlo scheme simultaneously detrending the light curves and the RV time series. We firmly detected the Keplerian signal of TOI-1055 b, deriving a planetary mass of $M_b=20.4_{-2.5}^{+2.6}\,M_{\oplus}$ ($\sim$12%). This value agrees with one of the two estimates in the literature, but it is significantly more precise. Thanks to TESS transit light curve combined with the exquisite CHEOPS photometry, we also derived a planetary radius of $R_b=3.490_{-0.064}^{+0.070}\,R_{\oplus}$ ($\sim$1.9%). Our mass and radius measurements imply a mean density of $\rho_b=2.65_{-0.35}^{+0.37}$ g cm$^{-3}$ ($\sim$14%). We further inferred the planetary structure and found that TOI-1055 b very likely hosts a substantial gas envelope with a mass of $0.41^{+0.34}_{-0.20}$ M$_\oplus$ and a thickness of $1.05^{+0.30}_{-0.29}$ R$_\oplus$.

The number counts of homogeneous samples of radio sources are a tried and true method of probing the large scale structure of the Universe, as most radio sources outside the galactic plane are at cosmological distances. As such they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although the cosmic radio dipole matches the direction of the CMB dipole, it has a significantly larger amplitude. This result challenges our assumption of the Universe being isotropic, which can have large repercussions for the current cosmological paradigm. Though significant measurements have been made, sensitivity to the radio dipole is generally hampered by systematic effects that can cause large biases in the measurement. Here we assess these systematics with data from the MeerKAT Absorption Line Survey (MALS). We present the analysis of ten MALS pointings, focusing on systematic effects that could lead to an inhomogeneous catalogue. We describe the calibration and creation of full band continuum images and catalogues, producing a combined catalogue containing 16,313 sources and covering 37.5 square degrees of sky down to a sensitivity of 10 $\mu$Jy/beam. We measure the completeness, purity, and flux recovery statistics for these catalogues using simulated data. We investigate different source populations in the catalogues by looking at flux densities and spectral indices, and how they might influence source counts. Using the noise characteristics of the pointings, we find global measures that can be used to correct for the incompleteness of the catalogue, producing corrected number counts down to 100 - 200 $\mu$Jy. We show that we can homogenise the catalogues and properly account for systematic effects. We determine that we can measure the dipole to $3\sigma$ significance with 100 MALS pointings.

Earth is the only known habitable planet and it serves as a testbed to benchmark the observations of temperate and more Earth-like exoplanets. It is required to observe the disc-integrated signatures of Earth for a large range of phase angles, resembling the observations of an exoplanet. In this work, an AOTF (Acousto-Optic Tunable Filter) based experiment is designed to observe the spectro-polarimetric signatures of Earth. The results of spectroscopic and polarimetric laboratory calibration are presented here along with a brief overview of a possible instrument configuration. Based on the results of the spectro-polarimetric calibration, simulations are carried out to optimize the instrument design for the expected signal levels for various observing conditions. The usefulness of an AOTF based spectro-polarimeter is established from this study and it is found that, in the present configuration, the instrument can achieve a polarimetric accuracy of $<0.3$\% for linear polarization for an integration time of 100 ms or larger. The design configuration of the instrument and the planning of conducting such observations from Lunar orbit are discussed.

Purpose: The Low-Energy X-ray telescope (LE) is a main instrument of the Insight-HXMT mission and consists of 96 Swept Charge Devices (SCD) covering the 1-10 keV energy band. The energy gain and resolution are continuously calibrated by analysing Cassiopeia A (Cas A) and blank sky data, while the effective areas are also calibrated with the observations of the Crab Nebula. In this paper, we present the evolution of the in-orbit performances of LE in the first 5 years since launch. Methods: The Insight-HXMT Data Analysis Software package (HXMTDAS) is utilized to extract the spectra of Cas A, blank sky, and Crab Nebula using different Good Time Interval (GTI) selections. We fit a model with a power-law continuum and several Gaussian lines to different ranges of Cas A and blank sky spectra to get peak energies of their lines through xspec. After updating the energy gain calibration in CALibration DataBase (CALDB), we rerun the Cas A data to obtain the energy resolution. An empirical function is used to modify the simulated effective areas so that the background-subtracted spectrum of the Crab Nebula can best match the standard model of the Crab Nebula. Results: The energy gain, resolution, and effective areas are calibrated every month. The corresponding calibration results are duly updated in CALDB, which can be downloaded and used for the analysis of Insight-HXMT data. Simultaneous observations with NuSTAR and NICER can also be used to verify our derived results. Conclusion: LE is a well calibrated X-ray telescope working in 1-10 keV band. The uncertainty of LE gain is less than 20 eV in 2-9 keV band and the uncertainty of LE resolution is less than 15eV. The systematic errors of LE, compared to the model of the Crab Nebula, are lower than 1.5% in 1-10 keV.

We aim to catalog all dust particles collected and analyzed by MIDAS, together with their main statistical properties such as size, height, basic shape descriptors, and collection time. Furthermore, we aim to present the scientific results that can be extracted from the catalog (e.g., the size distribution and statistical characteristics of cometary dust particles). The existing MIDAS particle catalog has been greatly improved by a careful re-analysis of the AFM images, leading to the addition of more dust particles and a detailed description of the particle properties. The catalog documents all images of identified dust particles and includes a variety of derived information tabulated one record per particle. Furthermore, the best image of each particle was chosen for subsequent studies. Finally, we created dust coverage maps and clustering maps of the MIDAS collection targets and traced any possible fragmentation of collected particles with a detailed algorithm. The revised MIDAS catalog includes 3523 MIDAS particles in total, where 1857 particles are expected to be usable for further analysis (418 scans of particles before perihelion + 1439 scans of particles after perihelion, both after the removal of duplicates), ranging from about 40 nm to about 8 {\mu}m in size. The mean value of the equivalent radius derived from the 2D projection of the particles is 0.91 {\pm} 0.79 {\mu}m. A slightly improved equivalent radius based on the particle's volume coincides in the range of uncertainties with a value of 0.56 {\pm} 0.45 {\mu}m. We note that those sizes and all following MIDAS particle size distributions are expected to be influenced by the fragmentation of MIDAS particles upon impact on the collection targets. Furthermore, fitting the slope of the MIDAS particle size distribution with a power law of a r {^b} yields an index b of {\sim} -1.67 to -1.88.

We built Galactic open star cluster mass functions (CMFs) for different age sub-samples in the wider solar neighbourhood. We present a simple cluster formation and evolution model to reproduce the main features of the CMFs. We used an unbiased sample of 2227 clusters of the Milky Way Star Cluster (MWSC) catalogue, which occupy the heliocentric cylinders with magnitude-dependent completeness radii of 1-5 kpc. We derived tidal masses of clusters with an accuracy of 70%. Our cluster formation and evolution model is based on the cluster initial mass function, the cluster formation rate, cluster mass loss due to stellar evolution and the clusters' dynamical evolution in the Galactic tidal field. The obtained tidal masses have been added to the MWSC catalogue. A general CMF (GCMF), built for all cluster ages around the Sun, extends over four decades in mass. The high-mass slope is +1.14. The CMFs for different age groups show the same high-mass slopes, while the low-mass slope is flat for the youngest sub-sample and about -0.7 for the others. The sub-samples inside and outside the solar Galactocentric radius are consistent with the GCMF, once the exponential decline of the Galactic disc density is taken into account. The model suggests star formation with low efficiency of 15%, where 10% of stars remain bound in a cluster after gas expulsion and violent relaxation. The cluster formation rate required to reproduce the observed age-mass-distribution is 0.4 solar masses per square pc and Gyr. The obtained high-mass slope of the GCMF for the wide solar neighbourhood is similar to slopes determined in nearby galaxies. The MWSC catalogue supports models with low star-formation efficiency, where 90% of stars are lost quickly after gas expulsion. The cluster formation rate corresponds to open clusters' contribution to the stellar content of the thin disc of 30%.

Mutually misaligned circumbinary planets may form in a warped or broken gas disc or from later planet-planet interactions. With numerical simulations and analytic estimates we explore the dynamics of two circumbinary planets with a large mutual inclination. A coplanar inner planet causes prograde apsidal precession of the binary and the stationary inclination for the outer planet is higher for larger outer planet orbital radius. In this case a coplanar outer planet always remains coplanar. On the other hand, a polar inner planet causes retrograde apsidal precession of the binary orbit and the stationary inclination is smaller for larger outer planet orbital radius. For a range of outer planet semi-major axes, an initially coplanar orbit is librating meaning that the outer planet undergoes large tilt oscillations. Circumbinary planets that are highly inclined to the binary are difficult to detect -- it is unlikely for a planet to have an inclination below the transit detection limit in the presence of a polar inner planet. These results suggest that there could be a population of circumbinary planets that are undergoing large tilt oscillations.

The production site of gamma rays in blazars is closely related to their interaction with the photon fields surrounding the active galactic nucleus. In this work we discuss an indirect method that may help to unveil the presence of ambient structures in BL Lacs through the analysis of their gamma-ray spectrum. Passing through structures at different distances from the black hole, gamma rays interact with the corresponding photon fields via gamma-gamma pair production, producing absorption features in their spectral energy distribution. An interaction of the gamma-ray photons with a putative broad-line region may reduce the gamma-ray flux only if its production site were very close to the central engine. On the other hand, if jet photons interact with optical-UV seed photons produced by a pc-scale narrow-line region, the consequent gamma-gamma process may cause absorption features at a few hundreds GeV. Sources with spectra reaching TeV energies, such as HBLs and EHBLs (extreme blazars), may represent exceptional probes to investigate this topic. In this regard, we discuss recent observations of sources which may show evidence of such absorption features in their gamma-ray spectra. Finally, we discuss how sub-TeV absorption features in the spectra of BL Lacs may affect their broadband modeling, and eventually represent a powerful diagnostic tool to constrain the gamma-ray production site and the jet environment.

The control system and the entire architecture of the High Resolution Spectrograph (SARG) for the Italian National Telescope "Galileo" (TNG) are here described. The concept of SARG instrument controls is similar to that of the other TNG instruments, in particular the CCD detector driving and the image acquisition use the same TNG standard boards and the same selected bus: the VME. The link between the SARG VME and the other telescope components is based on the same GATE software that guarantees the compatibility with the entire distributed TNG software. The control of the moving parts as well as the other parts of instrument that is the lamp controller and the temperature sensors, is based on a commercial controller connected to the system through a serial link. Furthermore a specialized software running on a PC has been realized to test the rotating tables independently of the VME system. Test of accuracy and repeatability of the positioning were done and some results are presented.

The amount and distribution of water on the lunar surface are related to the input and production of water by solar wind and meteoroid bombardment, balanced by photodestruction and mobility across the surface. Using the Stratospheric Observatory for Infrared Astronomy (SOFIA), we imaged the 6.1 micron feature that uniquely traces molecular water, covering 1/4 of the lunar nearside surface south of -60 degrees latitude with 5 km resolution on 2022 Feb 17 UTC. The water feature strength varies significantly across the region, being drier at +28 degrees longitude to more wet (~170 ppm) at -7 degrees longitude, and also decreasing toward the pole. Significant local enhancements are found, associated with south-facing, high-altitude topographic features. This includes relatively high H2O concentration in a "wet ridge" just north of Curtius crater; the south-facing, northern, inner rims of most prominent craters; the south face of the central peak of Moretus crater; and permanently-shadowed polar regions.

To explain Venus' young surface age and lack of plate tectonics, Venus' tectonic regime has often been proposed to be either an episodic-lid regime with global lithospheric overturns, or an equilibrium resurfacing regime with numerous volcanic and tectonic activities. Here, we use global 2-D thermochemical convection models with realistic parameters, including rheology (dislocation creep, diffusion creep, and plastic yielding), an experiment-based plagioclase (An$_{75}$) crustal rheology, and intrusive magmatism, to investigate the tectonics and mantle evolution of Venus. We find that surface tectonics is strongly affected by crustal rheology. With a ''weak'' plagioclase-rheology crust, models exhibit episodic overturns but with continuously high surface mobility and high distributed surface strain rates between overturns, leading to a new tectonic regime that we name ''deformable episodic lid''. On the other hand, olivine-crustal-rheology models exhibit either standard episodic-lid tectonics, i.e. with mobility that is high during overturns and near zero between overturns, or stagnant-lid tectonics, i.e. with near-zero mobility over the entire model time. Also, a combination of plagioclase crustal rheology and dislocation creep can weaken the lithosphere sufficiently to facilitate lithospheric overturns without applying plastic yielding. Internally, the composition-dependent density profile results in a ''basalt barrier'' at the mantle transition zone, which strongly affects Venus' mantle evolution. Only strong plumes can penetrate this basalt barrier and cause global lithospheric overturns. This basalt barrier also causes global internal episodic overturns that generate global volcanic resurfacing in stagnant-lid models, which suggests a new resurfacing mechanism (we name it ''stagnant episodic-volcanic-resurfacing'') that does not involve lithospheric overturns.

Exoplanetary science is a very active field of astronomy nowadays, with questions still opened such as how planetary systems form and evolve (occurrence, process), why such a diversity of exoplanets is observed (mass, radius, orbital parameters, temperature, composition), and what are the interactions between planets, circumstellar disk and their host star. Several complementary methods are used for the detection of exoplanets. Among these, imaging aims at the direct detection of the light reflected, scattered or emitted by exoplanets and circumstellar disks. This allows their spectral and polarimetric characterization. Such imaging remains challenging because of the large luminosity ratio (1e4-1e10$) and the small angular separation (fraction of an arcsecond) between the star and its environment. Over the past two decades, numerous techniques, including coronagraphy, have been developed to make exoplanet imaging a reality. This review gives a broad overview of the subsystems that make up a coronagraphic instrument for imaging exoplanetary systems. It is especially intended for non-specialists or newcomers in the field. We explain the principle of coronagraphy and propose a formalism to understand their behavior. We discuss the impact of wavefront aberrations on the performance of coronagraphs and how they induce stellar speckles in the scientific image. Finally, we present instrumental and signal processing techniques used for on-sky minimization or a posteriori calibration of these speckles in order to improve the performance of coronagraphs.

Astronomy is entering the era of large surveys of the variable sky such as Zwicky Transient Facility (ZTF) and the forthcoming Legacy Survey of Space and Time (LSST) which are intended to produce up to a million alerts per night. Such an amount of photometric data requires efficient light-curve pre-processing algorithms for the purposes of subsequent data quality cuts, classification, and characterization analysis. In this work, we present the new library "light-curve" for Python and Rust, which is intended for feature extraction from light curves of variable astronomical sources. The library is suitable for machine learning classification problems: it provides a fast implementation of feature extractors, which outperforms other public available codes, and consists of dozens features describing shape, magnitude distribution, and periodic properties of light curves. It includes not only features which had been shown to provide a high performance in classification tasks, but also new features we developed to improve classification quality of selected types of objects. The "light-curve" library is currently used by the ANTARES, AMPEL, and Fink broker systems for analyzing the ZTF alert stream, and has been selected for use with the LSST.

Decline and recovery timescales surrounding eclipse are indicative of the controlling physical processes in Io's atmosphere. Recent studies have established that the majority of Io's molecular atmosphere, SO2 and SO, condenses during its passage through Jupiter's shadow. The eclipse response of Io's atomic atmosphere is less certain, having been characterized solely by ultraviolet aurorae. Here we explore the response of optical aurorae for the first time. We find oxygen to be indifferent to the changing illumination, with [O I] brightness merely tracking the plasma density at Io's position in the torus. In shadow, line ratios confirm sparse SO2 coverage relative to O, since their collisions would otherwise quench the emission. Io's sodium aurora mostly disappears in eclipse and e-folding timescales, for decline and recovery differ sharply: ~10 minutes at ingress and nearly 2 hr at egress. Only ion chemistry can produce such a disparity; Io's molecular ionosphere is weaker at egress due to rapid recombination. Interruption of a NaCl+ photochemical pathway best explains Na behavior surrounding eclipse, implying that the role of electron impact ionization is minor relative to photons. Auroral emission is also evident from potassium, confirming K as the major source of far red emissions seen with spacecraft imaging at Jupiter. In all cases, direct electron impact on atomic gas is sufficient to explain the brightness without invoking significant dissociative excitation of molecules. Surprisingly, the nonresponse of O and rapid depletion of Na is opposite the temporal behavior of their SO2 and NaCl parent molecules during Io's eclipse phase.

Properties of the to-date-observed binary black hole (BBH) merger events suggest a preference towards spin-orbit aligned mergers. Naturally, this has caused widespread interest and speculations regrading implications on various merger formation channels. Here we show that (i) not only the BBH-merger population from isolated binaries, but also (ii) BBH population formed in young massive clusters (YMC) would possess an asymmetry in favour of aligned mergers, in the distribution of the events' effective spin parameter ($\chi_{\rm eff}$). In our analysis, we utilize BBH-merger outcomes from state-of-the-art N-body evolutionary models of YMCs and isolated binary population synthesis. We incorporate, for the first time in such an analysis, misalignments due to both natal kicks and dynamical encounters. The YMC $\chi_{\rm eff}$ distribution has a mean (an anti-aligned merger fraction) of $\langle\chi_{\rm eff}\rangle\leq0.05$ ($f_X-\approx40\%$) which is smaller (larger) than but consistent with the observed asymmetry of $\langle\chi_{\rm eff}\rangle\approx0.06$ ($f_X-\approx28\%$). In contrast, isolated binaries alone tend to produce a much stronger asymmetry; for the tested physical models, $\langle\chi_{\rm eff}\rangle\approx0.25$ and $f_X-\lesssim7\%$. Although the YMC $\chi_{\rm eff}$ distribution is more similar to the observed counterpart, none of the channels correctly reproduce the observed distribution. Our results suggest that further extensive model explorations for both isolated-binary and dynamical channels as well as better observational constraints are necessary to understand the physics of 'the symmetry breaking' of the BBH-merger population.

Gravitons radiated from light, evaporating black holes contribute to the stochastic background of gravitational waves. The spectrum of such emission depends on both the mass and the spin of the black holes, as well as on the redshifting that occurs between the black hole formation and today. This, in turn, depends on the expansion history of the universe, which is largely unknown and unconstrained at times prior to the synthesis of light elements. Here, we explore the features of the stochastic background of gravitational waves from black hole evaporation under a broad range of possible early cosmological histories. The resulting gravitational wave signals typically peak at very high frequencies, and offer opportunities for proposed ultra-high frequency gravitational wave detectors. Lower-frequency peaks are also possible, albeit with a suppressed intensity that is likely well below the threshold of detectability. We find that the largest intensity peaks correspond to cosmologies dominated by fluids with equations of state that have large pressure-to-density ratios. Such scenarios can be constrained on the basis of violation of $\Delta N_{\rm eff}$ bounds.

If dark matter is composed of axions, then axion stars form in the cores of dark matter halos. These stars are unstable above a critical mass, decaying to radio photons that heat the intergalactic medium, offering a new channel for axion indirect detection. We recently provided the first accurate calculation of the axion decay rate due to axion star mergers. In this work we show how existing data concerning the CMB optical depth leads to strong constraints on the axion photon coupling in the mass range $10^{-14}\,{\rm eV}\lesssim m_a\lesssim 10^{-8}\,{\rm eV}$. Axion star decays lead to efficient reionization of the intergalactic medium during the dark ages. By comparing this non-standard reionization with Planck legacy measurements of the Thompson optical width, we show that couplings in the range $10^{-14}\,{\rm GeV}^{-1} \lesssim g_{a\gamma\gamma} \lesssim 10^{-10}\,{\rm GeV}^{-1}$ are excluded for our benchmark model of axion star abundance. Future measurements of the 21cm emission of neutral hydrogen at high redshift could improve this limit by an order of magnitude or more, providing complementary indirect constraints on axion dark matter in parameter space also targeted by direct detection haloscopes.

We present a new search for weakly interacting massive particles utilizing ten years of public IceCube data, setting more stringent bounds than previous IceCube analysis on massive dark matter to neutrino annihilation. We also predict the future potential of the new neutrino observatory, P-ONE, showing that it may even exceed the sensitivities of gamma-ray searches by about 1-2 orders of magnitude in 1-10 TeV regions. This analysis considers the diffuse dark matter self-annihilation to neutrinos via direct and indirect channels, from the galactic dark matter halo and extra-galactic sources. We also predict that P-ONE will be capable of pushing these bounds further than IceCube, even reaching the thermal relic abundance utilizing a galactic center search for extended run-time.

Observational astronomy is plagued with selection effects that must be taken into account when interpreting data from astronomical surveys. Because of the physical limitations of observing time and instrument sensitivity, datasets are rarely complete. However, determining specifically what is missing from any sample is not always straightforward. For example, there are always more faint objects (such as galaxies) than bright ones in any brightness-limited sample, but faint objects may not be of the same kind as bright ones. Assuming they are can lead to mischaracterizing the population of objects near the boundary of what can be detected. Similarly, starting with nearby objects that can be well observed and assuming that objects much farther away (and sampled from a younger universe) are of the same kind can lead us astray. Demographic models of galaxy populations can be used as inputs to observing system simulations to create ``mock'' catalogues that can be used to characterize and account for multiple, interacting selection effects. The use of simulations for this purpose is common practice in astronomy, and blurs the line between observations and simulations; the observational data cannot be interpreted independent of the simulations. We will describe this methodology and argue that astrophysicists have developed effective ways to establish the reliability of simulation-dependent observational programs. The reliability depends on how well the physical and demographic properties of the simulated population can be constrained through independent observations. We also identify a new challenge raised by the use of simulations, which we call the ``problem of uncomputed alternatives.'' Sometimes the simulations themselves create unintended selection effects when the limits of what can be simulated lead astronomers to only consider a limited space of alternative proposals.