Papers for Thursday, Jun 16 2022

The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ($E$, $L_z$, $J_r$, $J_z$) is one method of identifying accreted structure. We produce 32 simulated stellar halos using particle tagging with the \textit{Caterpillar} simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover $6\%$ of UFDs remnants at best. These are remnants that accreted recently, $z_{\text{accretion}}\lesssim 0.5$. Based on these results, we make recommendations to help guide the search for dynamically-linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and $J_r$, providing a way to help identify real vs. fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.

Using the largest sample of galaxies observed with an optical integral field unit (IFU, the SDSS-IV MaNGA survey, $\sim$10000 targets), we derive the radial distribution of the physical properties obtained from the stellar continuum and the ionized-gas emission lines. Given the large sample, we are able to explore the impact of the total stellar mass and morphology by averaging those radial distributions for different bins of both global properties. We use a piece-wise analysis to characterize the slopes of the gradients from those properties at different galactocentric distances. In general we find that most of the properties -- derived from both the stellar continuum and the ionized gas emission lines -- exhibit a negative gradient with a secondary impact by global properties such as the total stellar mass or morphology. Our results confirm the intimate interplay between the properties of the stellar component and those of the ionized gas at local (kpc) scales in order to set the observed gradients. Furthermore, the resemblance of the gradients for similar global properties (in particular for the stellar parameters) indicates statistical similar histories of star formation and chemical enrichment with an initial radial gas distribution following the potential of the galaxy.

Low-mass early-type galaxies, including dwarf spheroidals (dSph) and brighter dwarf ellipticals (dE), dominate the galaxy population in groups and clusters. Recently, an additional early-type population of more extended ultra-diffuse galaxies (UDGs) has been identified, sparking a discussion on the potential morphological and evolutionary connections between the three classifications. Here we present the first measurements of spatially resolved stellar kinematics from deep integrated-light spectra of KDG 64 (UGC 5442), a large dSph galaxy in the M 81 group. From these data we infer stellar population properties and dark matter halo parameters using Jeans dynamical modelling. We find an old, metal-poor stellar population with no young stars and a dark matter mass fraction of ~ 90% within the half-light radius. These properties and the position of KDG 64 on the Fundamental Plane indicate that it is a local analogue of smaller UDGs in the Coma and Virgo clusters and is probably a transitional dSph-UDG object. Its evolutionary path cannot be uniquely established from the existing data, but we argue that supernovae feedback and tidal heating played key roles in shaping KDG 64.

We study the structure of atomic hydrogen (HI) in the host galaxy of GRB 171205A / SN 2017iuk at z=0.037 through HI 21cm emission line observations with the Karl G. Jansky Very Large Array. These observations reveal unusual morphology and kinematics of the HI in this otherwise apparently normal galaxy. High column density, cold HI is absent from an extended North-South region passing by the optical centre of the galaxy, but instead is extended towards the South, on both sides of the galaxy. Moreover, the HI kinematics do not show a continuous change along the major axis of the galaxy as expected in a classical rotating disk. We explore several scenarios to explain the HI structure and kinematics in the galaxy: feedback from a central starburst and/or an active galactic nucleus, ram pressure stripping, accretion, and tidal interaction from a companion galaxy. All of these options are ruled out. The most viable remaining explanation is the penetrating passage of a satellite through the disk only a few Myr ago, redistributing the HI in the GRB host without yet affecting its stellar distribution. It can also lead to the rapid formation of peculiar stars due to a violent induced shock. The location of GRB 171205A in the vicinity of the distorted area suggests that its progenitor star(s) originated in extreme conditions that share the same origin as the peculiarities in HI. This could explain the atypical location of GRB 171205A in its host galaxy.

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e. kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic Particle-in-Cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e. MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron-ion simulations, on which we will focus in future work.

We present here the analysis performed using the pyPipe3D pipeline for the final MaNGA dataset included in the SDSS seventeenth data-release. This dataset comprises more than 10,000 individual datacubes, being the integral field spectroscopy galaxy survey with the largest number of galaxies. pyPipe3D processes the IFS datacubes to extract spatially-resolved spectroscopic properties of both the stellar population and the ionized-gas emission lines. A brief summary of the properties of the sample and the characteristics of the analyzed data are included. The article provides details on (i) the performed analysis, (ii) a description of the pipeline, (iii) the adopted stellar population library, (iv) the morphological and photometric analysis, (v) the adopted datamodel for the derived spatially resolved properties and (vi) the individual integrated and characteristic galaxy properties included in a final catalog. Comparisons with results from a previous version of the pipeline for earlier data releases and from other tools using this dataset are included. A practical example on how to use of the full dataset, and final catalog illustrates how to handle the delivered product. Our full analysis can be accessed and downloaded from the webpage this http URL

A thermonuclear explosion triggered by a helium-shell detonation on a carbon-oxygen white dwarf core has been predicted to have strong UV line blanketing at early times due to the iron-group elements produced during helium-shell burning. We present the photometric and spectroscopic observations of SN 2016dsg, a sub-luminous peculiar Type I SN consistent with a thermonuclear explosion involving a thick He shell. With a redshift of 0.04, the $i$-band peak absolute magnitude is derived to be around -17.5. The object is located far away from its host, an early-type galaxy, suggesting it originated from an old stellar population. The spectra collected after the peak are unusually red, show strong UV line blanketing and weak O I $\lambda$7773 absorption lines, and do not evolve significantly over 30 days. An absorption line around 9700-10500 \AA is detected in the near-infrared spectrum and is likely from the unburnt helium in the ejecta. The spectroscopic evolution is consistent with the thermonuclear explosion models for a sub-Chandrasekhar mass white dwarf with a thick helium shell, while the photometric evolution is not well described by existing models.

The optical properties of the second generation dust that we observe in debris disks remain quite elusive, whether it is the absorption efficiencies at millimeter wavelengths or the (un)polarized phase function at near-infrared wavelengths. Thankfully the same particles are experiencing forces that are size dependent (e.g., radiation pressure), and with high angular resolution observations we can take advantage of this natural spatial segregation. Observations at different wavelengths probe different ranges of sizes, and there is therefore a great synergy in multi-wavelength observations to better constrain the optical properties of the particles. We present a new approach to simultaneously model SPHERE and ALMA observations and apply it to the debris disk around HD\,32297, putting the emphasis on the spatial distribution of the grains with different $\beta$ values. This modeling approach requires few assumptions on the actual sizes of the particles and the interpretation can therefore be done a posteriori. We find that the ALMA observations are best reproduced with a combination of small and large $\beta$ values ($0.03$ and $0.42$) while the SPHERE observations require several intervals of $\beta$ values. We discuss the nature of the halo previously reported in ALMA observations, and hypothesize it could be caused by over-abundant $\mu$m-sized particles (the over-abundance being the consequence of their extended lifetime). We model the polarized phase function at near-infrared wavelengths and fluffy aggregates larger than a few $\mu$m provide the best solution. Comparing our results with comets of the solar system, we postulate that the particles released in the disk originate from rather pristine cometary bodies (to avoid compaction of the fluffy aggregates) and are then set on highly eccentric orbits, which could explain the halo detected at long wavelengths.

We propose an analytic model, CuspCore II, for the response of dark matter (DM) haloes to central gas ejection, as a mechanism for generating DM-deficient cores in dwarfs and high-z massive galaxies. We test this model and three other methods using idealized N-body simulations. The current model is physically justified and provides more accurate predictions than the earlier version, CuspCore I (Freundlich et al. 2020). The CuspCore model assumes an instantaneous change of potential, followed by a relaxation to a new Jeans equilibrium. The relaxation turns out to be violent relaxation during the first orbital period, followed by phase mixing. By tracing the energy diffusion dE=dU(r) iteratively, the model reproduces the simulated DM profiles with ~10% accuracy or better. A method based on adiabatic invariants shows similar precision for moderate mass change but underestimates the DM expansion for strong gas ejection. A method based on a simple empirical relation between DM and total mass ratios makes slightly inferior predictions. The crude assumption used in CuspCore I, of energy conservation for shells that encompass a fixed DM mass, turns out to underestimate the DM response, which can be partially remedied by introducing an alternative "energy" definition. Our model is being generalized to address the differential response of a multi-component system of stars and DM in the formation of DM-deficient galaxies.

Luminous red novae (LRNe) are astrophysical transients believed to be caused by the partial ejection of a binary star's common envelope (CE) and the merger of its components. The formation of the CE is likely to occur during unstable mass transfer, initiated by a primary star which is evolving off the main sequence (a Hertzsprung gap star) and a lower mass companion. In agreement with observations, theoretical studies have shown that outflows from the pre-CE phase produce a detectable brightening of the progenitor system a few years before the ejection event. Based on these assumptions, we present a method to identify Galactic LRNe precursors, the resulting precursor candidates, and our follow-up analysis to uncover their nature. We begin by constructing a sample of progenitor systems, i.e. Hertzsprung gap stars, by statistically modelling the density of a colour magnitude diagram formed from "well behaved" Gaia DR2 sources. Their time-domain evolution from the Zwicky Transient Facility (ZTF) survey is used to search for slowly brightening events, as pre-CE precursor candidates. The nature of the resulting candidates is further investigated using archival data and our own spectroscopic follow-up. Overall, we constructed a sample of $\sim5.4\times{10^4}$ progenitor sources, from which 21 were identified as candidate LRNe precursors. Further analysis revealed 16 of our candidates to be H$\alpha$ emitters, with their spectra often suggesting hotter (albeit moderately extincted) A-type or B-type stars. Because of their long-term variability in optical and mid-infrared wavelengths, we propose that many of our candidates are mass-transferring binaries with compact companions surrounded by dusty circumstellar disks or alternatively magnetically active stellar merger remnants.

One outstanding problem in extrasolar planet studies is why no co-orbital exoplanets have been found, despite numerous searches among the many known planetary systems, many of them in other mean-motion resonances. Here we examine the hypothesis that dissipation of energy by tides in Trojan planets is preventing their survival. The Appendix of this paper generalizes the conventional theory of tides to include tidal forces independent of dissipation, as well as the effects of one body on tides raised by another. The main text applies this theory to a model system consisting of a primary of stellar mass, a secondary of sub-stellar mass in a circular orbit about the primary, and a much lighter Trojan planet librating with small amplitude about an equilateral point of the system. Next, we linearize the equations of motion about the Trojan points, including the tidal forces, and solve for the motion of the Trojan. The results indicate that tides damp out the Trojan's motion perpendicular to the orbital plane of the primary and secondary, as well as its epicycles due to its eccentricity; but they pump up the amplitude of its tadpole librations exponentially. We then verify our analytic solutions by integrating the non-linearized equations of motion numerically for several sample cases. In each case, we find that the librations grow until the Trojan escapes its libration, which leads to a close encounter with either the primary or the secondary.

Using observations of X-ray pulsar Her X-1 by the Imaging X-ray Polarimetry Explorer, we report on a highly significant detection of the polarization signal from an accreting neutron star. The observed degree of the polarization of $\sim10$% is found to be far below theoretical expectations for this object, and stays low throughout the spin cycle of the pulsar. Both the polarization degree and the angle exhibit variability with pulse phase, which allowed us to measure the pulsar spin position angle and magnetic obliquity of the neutron star, which is an essential step towards detailed modeling of intrinsic emission of X-ray pulsars. Combining our results with the optical polarimetric data, we find that the spin axis of the neutron star and the angular momentum of the binary orbit are misaligned by at least $\sim$20 deg, which is a strong argument in support of the models explaining stability of the observed super-orbital variability with the precession of the neutron star.

A new method for reconstruction of coronal magnetic fields as force-free fields (FFFs) is presented. Our method employs poloidal and toroidal functions to describe divergence-free magnetic fields. This magnetic field representation naturally enables us to implement the boundary conditions at the photospheric boundary, i.e., the normal magnetic field and the normal current density there, in a straightforward manner. At the upper boundary of the corona, a source-surface condition can be employed, which accommodates magnetic flux imbalance at the bottom boundary. Although our iteration algorithm is inspired by extant variational methods, it is non-variational and requires much less iteration steps than most of them. The computational code based on our new method is tested against the analytical FFF solutions by Titov & D\'{e}moulin (1999). It is found to excel in reproducing a tightly wound flux rope, a bald patch and quasi-separatrix layers with a hyperbolic flux tube.

We characterize the optical counterparts to the compact X-ray source population within the nearby spiral galaxy M81 using multi-band Hubble Space Telescope (HST) imaging data. By comparing the optical luminosities and colors measured for candidate donor stars and host clusters to stellar and cluster evolutionary models, respectively, we estimate the likely masses and upper age limits of the field and cluster X-ray binaries. We identify 15 low-mass X-ray binaries (i.e. donor star mass $\leq$ 3 solar masses) within ancient globular clusters, as well as 42 candidate high-mass X-ray binaries (i.e. donor star mass $\geq$ 8 solar masses). To estimate the likelihood of misclassifications, we inject 4,000 artificial sources into the HST mosaic image and conclude that our classifications of globular clusters and high-mass X-ray binaries are reliable at the >90% level. We find that globular clusters that host X-ray binaries are on average more massive and more compact than globular clusters that do not. However, there is no apparent correlation between the X-ray brightness of the clusters and their masses or densities, nor are X-ray binary hosts more X-ray luminous than the general field population of low-mass X-ray binaries. This work represents one of the first in-depth analyses of the population of X-ray binaries within globular clusters in a spiral galaxy.

The spontaneous evolution of magnetic reconnection in generalized situations (with thermodynamic asymmetry regarding the current sheet and magnetic shear) is investigated using two-dimensional magnetohydrodynamic simulation. We focus on the asymptotic state of temporal evolution, i.e., the self-similarly expanding phase (Nitta et al. 2001). 1) A long fast-mode shock is generated in front of the shorter plasmoid as in the shear-less thermodynamically asymmetric case; however, the sheared magnetic component weakens the shock. This fast shock may work as a particle acceleration site. 2) The shorter plasmoid-side plasma infiltrates the longer plasmoid across the current sheet. Then, the plasmas from both sides of the current sheet coexist on the same magnetic field lines in the longer plasmoid. This may result in efficient plasma mixing. 3) The thermodynamic asymmetry and magnetic shear drastically decrease the reconnection rate in many orders of magnitude.

The Deep Atmosphere Venus Investigation of Noble gases, Chemistry, and Imaging (DAVINCI) mission described herein has been selected for flight to Venus as part of the NASA Discovery Program. DAVINCI will be the first mission to Venus to incorporate science-driven flybys and an instrumented descent sphere into a unified architecture. The anticipated scientific outcome will be a new understanding of the atmosphere, surface, and evolutionary path of Venus as a possibly once-habitable planet and analog to hot terrestrial exoplanets. The primary mission design for DAVINCI as selected features a preferred launch in summer/fall 2029, two flybys in 2030, and descent sphere atmospheric entry by the end of 2031. The in situ atmospheric descent phase subsequently delivers definitive chemical and isotopic composition of the Venus atmosphere during a cloud-top to surface transect above Alpha Regio. These in situ investigations of the atmosphere and near infrared descent imaging of the surface will complement remote flyby observations of the dynamic atmosphere, cloud deck, and surface near infrared emissivity. The overall mission yield will be at least 60 Gbits (compressed) new data about the atmosphere and near surface, as well as first unique characterization of the deep atmosphere environment and chemistry, including trace gases, key stable isotopes, oxygen fugacity, constraints on local rock compositions, and topography of a tessera.

The activation energy for desorption (Edes) and that for surface diffusion (Esd) of adsorbed molecules on dust grains are two of the most important parameters for the chemistry in the interstellar medium. Although Edes is often measured by laboratory experiments, the measurement of Esd is sparse. Due to the lack of data, astrochemical models usually assume a simple scaling relation, Esd = fEdes, where f is a constant, irrespective of adsorbed species. Here, we experimentally measure Esd for CH4, H2S, OCS, CH3OH, and CH3CN on water-ice surfaces using an ultra-high-vacuum transmission electron microscope (UHV-TEM). Compiling the measured Esd values and Edes values from the literature, we find that the value of f ranges from ~0.2 to ~0.7, depending on the species. Unless f (or Esd) for the majority of species is available, a natural alternative approach for astrochemical models is running multiple simulations, varying f for each species randomly. In this approach, ranges of molecular abundances predicted by multiple simulations, rather than abundances predicted by each simulation, are important. We here run 10,000 simulations of astrochemical models of molecular clouds and protostellar envelopes, randomly assigning a value of f for each species. In the former case, we identify several key species whose Esd most strongly affects the uncertainties of the model predictions; Esd for those species should be investigated in future laboratory and quantum chemical studies. In the latter case, uncertainties in the Esd of many species contribute to the uncertainties in the model predictions.

Using the high resolution Milky Way-like model from Auriga simulation we study the chemical properties of the Galactic bulge, focusing on the metallicity difference between stars on the near side (in front of the Galactic center) and the far side (behind the Galactic center). In general, along certain sight lines the near side is more metal-rich than the far side, consistent with the negative vertical metallicity gradient of the disk, since the far side is located higher above the disk plane than the near side. However, at the region $l<0^\circ$ and $|b|\le6^\circ$, the near side is even more metal-poor than the far side, and their difference changes with the Galactic longitude. This is mainly due to the fact that stars around the minor axis of the bar are more metal-poor than those around the major axis. Since the bar is tilted, in the negative longitude region, the near side is mainly contributed by stars close to the minor axis region than the far side to result in such metallicity difference. We extract stars in the X-shape structure by identifying the overdensities in the near and far sides. Their metallicity properties are consistent with the results of the whole Galactic bulge. The boxy/peanut-shaped bulge can naturally explain the metallicity difference of the double red clump stars in observation. There is no need to involve a classical bulge component with different stellar populations.

In this study, we examine the driving mechanism for the atmospheric overturning circulation on dry, tidally-locked rocky planets without the condensation of water vapor or other species. We find that the main driving process is the radiative cooling of CO2 (or other non-condensable greenhouse gases) rather than CO2 greenhouse warming or stellar radiation. Stellar radiation is the ultimate mechanism but not the direct mechanism. Due to the combination of the uneven distribution in the stellar radiation and effective horizontal energy transports in the free troposphere, there is strong temperature inversion in the area away from the substellar region. This inversion makes CO2 to have a radiative cooling effect rather than a radiative warming effect for the atmosphere, same as that in the stratosphere of Earth's atmosphere. This cooling effect produces negative buoyancy and drives large-scale downwelling, supporting the formation of a global-scale overturning circulation. If CO2 is excluded from the atmosphere, the overturning circulation becomes very weak, regardless the level of stellar radiation. This mechanism is completely different from that for the atmospheric overturning circulation on Earth or on moist, tidally-locked rocky planets, where latent heat release and/or baroclinic instability are the dominated mechanisms. Our study improves the understanding of the atmospheric circulation on tidally-locked exoplanets and also on other dry planets, such as Venus and Mars in the solar system.

We apply the variational autoencoder (VAE) to the LAMOST-K2 low-resolution spectra to detect the magnetic activity of the stars in the K2 field. After the training on the spectra of the selected inactive stars, the VAE model can efficiently generate the synthetic reference templates needed by the spectral subtraction procedure, without knowing any stellar parameters. Then we detect the peculiar spectral features, such as chromospheric emissions, strong nebular emissions and lithium absorptions, in our sample. We measure the emissions of the chromospheric activity indicators, H$\alpha$ and Ca$~{\rm {\small II}}$ infrared triplet (IRT) lines, to quantify the stellar magnetic activity. The excess emissions of H$\alpha$ and Ca$~{\rm {\small II}}$ IRT lines of the active stars are correlated well to the rotational periods and the amplitudes of light curves derived from the K2 photometry. We degrade the LAMOST spectra to simulate the slitless spectra of the planned China Space Station Telescope (CSST) and apply the VAE to the simulated data. For cool active stars, we reveal a good agreement between the equivalent widths (EWs) of H$\alpha$ line derived from the spectra with two resolutions. The result indicates the ability of identifying the magnetically active stars in the future CSST survey, which will deliver an unprecedented large database of low-resolution spectra as well as simultaneous multi-band photometry of stars.

Precise radial velocity (PRV) surveys are important for the search of Earth analogs around nearby bright stars. Such planets induce a small stellar reflex motion with RV amplitude of $\sim$10 cm/s. Detecting such a small RV signal poses important challenges to instrumentation, data analysis, and the precision of astrophysical models to mitigate stellar jitter. In this work, we investigate an important component in the PRV error budget - the spectral contamination from the Earth's atmosphere (tellurics). We characterize the effects of telluric absorption on the RV precision and quantify its contribution to the RV budget over time and across a wavelength range of 350 nm - 2.5$\mu$m. We investigate the effectiveness in mitigating tellurics using simulated spectra of a solar twin star with telluric contamination over a year's worth of observations, and we extracted the RVs using two commonly adopted algorithms: dividing out a telluric model before performing cross-correlation or Forward Modeling the observed spectrum incorporating a telluric model. We assume various degrees of cleanness in removing the tellurics, including mimicking the lack of accurate knowledge of the telluric lines by using a mismatched line profile to model the "observed" tellurics. We conclude that the RV errors caused by telluric absorption can be suppressed to close to or even below the photon-limited precision in the optical region, especially in the blue, around 1-10 cm/s. At red through near-infrared wavelengths, however, the residuals of tellurics can induce an RV error on the m/s level even under the most favorable assumptions for telluric removal, leading to significant systematic noise in the RV time series and periodograms. If the red-optical or near-infrared becomes critical in the mitigation of stellar activity, systematic errors from tellurics can be eliminated with a space mission such as EarthFinder.

Extragalactic fast radio bursts (FRBs) are a new class of astrophysical transients with unknown origins that have become a main focus of radio observatories worldwide. FRBs are highly energetic ($\sim 10^{36}$-$10^{42}$ ergs) flashes that last for about a millisecond. Thanks to its broad bandwidth (400-800 MHz), large field of view ($\sim$200 sq. deg.), and massive data rate (1500 TB of coherently beamformed data per day), the Canadian Hydrogen Intensity Mapping Experiment / Fast Radio Burst (CHIME/FRB) project has increased the total number of discovered FRBs by over a factor 10 in 3 years of operation. CHIME/FRB observations are hampered by the constant exposure to radio frequency interference (RFI) from artificial devices (e.g., cellular phones, aircraft), resulting in $\sim$20% loss of bandwidth. In this work, we describe our novel technique for mitigating RFI in CHIME/FRB real-time intensity data. We mitigate RFI through a sequence of iterative operations, which mask out statistical outliers from frequency-channelized intensity data that have been effectively high-pass filtered. Keeping false positive and false negative rates at very low levels, our approach is useful for any high-performance surveys of radio transients in the future.

We present a new method for probabilistic generative modelling of stellar colour-magnitude diagrams (CMDs) to infer the frequency of binary stars and their mass-ratio distribution. The method invokes a mixture model to account for overlapping populations of single stars, binaries and outliers in the CMD. We apply the model to Gaia observations of the old open cluster, M67, and find a frequency $f_B(q > 0.5) = 0.258 \pm 0.019$ for binary stars with mass ratio greater than 0.5. The form of the mass-ratio distribution function rises towards higher mass ratios for $q > 0.3$.

The physics of early stellar evolution (e.g. accretion processes) is often not properly included in the calculations of pre-main-sequence models, leading to insufficient model grids and hence systematic errors in the results. We aim to investigate current and improved approaches for the asteroseismic modelling of pre-main-sequence delta Scuti stars. We calculated an extensive grid of pre-main-sequence models including the early accretion phase and used the resulting equilibrium models as input to calculate theoretical frequency spectra. These spectra were used to investigate different approaches in modelling echelle diagrams to find the most reliable methods. By applying Petersen diagrams, we present a simple algorithm to extract echelle diagrams from observed pulsation frequencies. We show that model grids with insufficient input physics and imperfect modelling approaches lead to underestimated uncertainties and systematic errors in the extracted stellar parameters. Our re-discussion of HD 139614 leads to different stellar parameters than the ones derived by Murphy et al. (2021). We performed a model comparison between this previous investigation and our results by applying the Akaike and Bayesian information criteria. While the results with regard to our 10-d model are inconclusive, they show (very) strong evidence of a 6-d model with fixed accretion parameters (leading to almost identical stellar parameters to those of the 10-d model) to be preferred over the model applied by Murphy et al. (2021). In general, our modelling approach can provide narrow constraints on the stellar parameters (\Delta R ~ 0.05 R_\odot, \Delta log g <~ 0.01, and \Delta M_\star ~ 0.1 M_\odot). The extensively tested modelling approaches and automatic extraction of echelle diagrams should allow us to study many more pre-main-sequence delta Scuti stars in the future and lead to reliable stellar parameters.

The X-ray emission from the knots of the kilo-parsec scale jet of active galactic nuclei (AGN) suggests the high energy emission process is different from the radio/optical counterpart. Interpretation based on the Inverse Compton scattering of cosmic microwave photons has been ruled out through Fermi gamma-ray observations for low redshift sources. As an alternate explanation, synchrotron emission from a different electron population is suggested. We propose a model considering the advected electron distribution from the sites of particle acceleration in AGN knots. This advected electron distribution is significantly different from the accelerated electron distribution and satisfies the requirement of the second electron population. The synchrotron emission from the accelerated and the advected electron distribution can successfully reproduce the observed radio to X-ray fluxes of the knots of 3C 273. For the chosen combination of the model parameters, the spectrum due to inverse Compton scattering of cosmic microwave photons falls within the Fermi gamma-ray upper limits.

This is the second paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the (super)massive regime. In the present work we focus on (super)massive BHs and provide an ab-initio computation of their mass function across cosmic times. We consider two main mechanisms to grow the central BH, that are expected to cooperate in the high-redshift star-forming progenitors of local massive galaxies. The first is the gaseous dynamical friction process, that can cause the migration toward the nuclear regions of stellar-mass BHs originated during the intense bursts of star formation in the gas-rich host progenitor galaxy, and the buildup of a central heavy BH seed $M_\bullet\sim 10^{3-5}\, M_\odot$ within short timescales $\lesssim$ some $10^7$ yr. The second mechanism is the standard Eddington-type gas disk accretion onto the heavy BH seed, through which the central BH can become (super)massive $M_\bullet\sim 10^{6-10}\, M_\odot$ within the typical star-formation duration $\lesssim 1$ Gyr of the host. We validate our semi-empirical approach by reproducing the observed redshift-dependent bolometric AGN luminosity functions and Eddington ratio distributions, and the relationship between the star-formation and the bolometric luminosity of the accreting central BH. We then derive the relic (super)massive BH mass function at different redshifts via a generalized continuity equation approach, and compare it with present observational estimates. Finally, we reconstruct the overall BH mass function from the stellar to the (super)massive regime, over more than ten orders of magnitudes in BH mass.

In this work, we extend the catalog of low-surface brightness (LSB) galaxies, including Ultra-Diffuse Galaxy (UDG) candidates, within $\approx 0.4R_{vir}$ of the Hydra I cluster of galaxies, based on deep images from the VST Early-type GAlaxy Survey (VEGAS). The new galaxies are found by applying an automatic detection tool and carrying out additional visual inspections of $g$ and $r$ band images. This led to the detection of 11 UDGs and 8 more LSB galaxies. For all of them, the cluster membership has been assessed using the color-magnitude relation derived for early-type giant and dwarf galaxies in Hydra I. The UDGs and new LSB galaxies found in Hydra I span a wide range of central surface brightness ($22.7 \lesssim \mu_{0,g} \lesssim 26.5$ mag/arcsec$^2$), effective radius ($0.6 \lesssim R_e \lesssim 4.0$ kpc) and color ($0.4 \leq g-r \leq 0.9$ mag), and have stellar masses in the range $\sim 5\times 10^6 - 2\times 10^8$M$_{\odot}$. The 2D projected distribution of both galaxy types is similar to the spatial distribution of dwarf galaxies, with over-densities in the cluster core and north of the cluster centre. They have similar color distribution and comparable stellar masses to the red dwarf galaxies. Based on photometric selection, we identify a total of 9 globular cluster candidates associated to the UDGs and 4 to the LSB galaxies, with the highest number of candidates in an individual UDG being three. We find that there are no relevant differences between dwarfs, LSB galaxies and UDGs: the structural parameters (that is surface brightness, size, colors, n-index) and GCs content of the three classes have similar properties and trends. This finding is consistent with UDGs being the extreme LSB tail of the size-luminosity distribution of dwarfs in this environment.

A key unknown of the Milky Way (MW) satellites is their orbital history, and, in particular, the time they were accreted onto the MW system since it marks the point where they experience a multitude of environmental processes. We present a new methodology for determining infall times, namely using a neural network (NN) algorithm. The NN is trained on MW-analogues in the EAGLE hydrodynamical simulation to predict if a dwarf galaxy is at first infall or a backsplash galaxy and to infer its infall time. The resulting NN predicts with 85\% accuracy if a galaxy currently outside the virial radius is a backsplash satellite and determines the infall times with a typical 68\% confidence interval of 4.4 Gyrs. Applying the NN to MW dwarfs with Gaia EDR3 proper motions, we find that all of the dwarfs within 300 kpc had been inside the Galactic halo. The overall MW satellite accretion rate agrees well with the theoretical prediction except for late times when the MW shows a second peak at a lookback time of 1.5 Gyrs corresponding to the infall of the LMC and its satellites. We also find that the quenching times for ultrafaint dwarfs show no significant correlation with infall time and thus supporting the hypothesis that they were quenched during reionisation. In contrast, dwarfs with stellar masses above $10^5~M_\odot$ are found to be consistent with environmental quenching inside the Galactic halo, with star-formation ceasing on average at $0.5^{+0.9}_{-1.2}$ Gyrs after infall.

Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make detections in massive sources hard to achieve. We present an application of the radiative transfer code POLARIS, with the goal to test the observability of the ortho-H$_2$D$^+$ transition $1_{10}$-$1_{11}$ (372.42 GHz) using simulations of high-mass collapsing cores that include deuteration chemistry. We analyzed an early and a late stage of the collapse of a 60 M$_{\odot}$ core, testing different source distances. For all cases, we generated synthetic single-dish and interferometric observations and studied the differences in both techniques. The column densities we derive are comparable to values reported for similar sources. These estimates depend on the extent over which they are averaged, and sources with compact emission they can be highly affected by beam dilution. Combined ALMA-ACA observations improve in signal-to-noise ratio and lead to better column density estimates as compared to ALMA alone. We confirm the feasibility to study ortho-H$_2$D$^+$ emission up to distances of 7 kpc. We provide a proof-of-concept of our framework for synthetic observations and highlight its importance when comparing numerical simulations with real observations. This work also proves how relevant it is to combine single-dish and interferometric measurements to derive appropriate source column densities.

In a magnetised plasma on scales well above ion kinetic scales, any constant-magnitude magnetic field, accompanied by parallel Alfv\'enic flows, forms a nonlinear solution in an isobaric, constant-density background. These structures, which are also known as spherically polarised Alfv\'en waves, are observed ubiquitously in the solar wind, presumably created by the growth of small-amplitude fluctuations as they propagate outwards in the corona. Here, we present a computational method to construct such solutions of arbitrary amplitude with general multi-dimensional structure, and explore some of their properties. The difficulty lies in computing a zero-divergence, constant-magnitude magnetic field, which leaves a single, quasi-free function to define the solution, while requiring strong constraints on any individual component of the field. Motivated by the physical process of wave growth in the solar wind, our method circumvents this issue by starting from low-amplitude Alfv\'enic fluctuations dominated by a strong mean field, then "growing" magnetic perturbations into the large-amplitude regime. We present example solutions with nontrivial structure in one, two, and three dimensions, demonstrating a clear tendency to form very sharp gradients or discontinuities, unless the solution is one dimensional. As well as being useful as an input for other calculations, particularly the study of parametric decay, our results provide a natural explanation for the extremely sharp field discontinuities observed across magnetic-field switchbacks in the low solar wind.

We report the results of an observation of low mass X-ray binary GX 3+1 with {\it AstroSat}'s Large Area X-ray Proportional Counter (LAXPC) and Soft X-ray Telescope (SXT) instruments on-board for the first time. We have detected one Type-1 thermonuclear burst ($\sim$ 15 s) present in the LAXPC 20 light curve, with a double peak feature at higher energies and our study of the hardness-intensity diagram reveals that the source was in a soft banana state. The pre-burst emission could be described well by a thermally Comptonised model component. The burst spectra is modelled adopting a time-resolved spectroscopic method using a single color blackbody model added to the pre-burst model, to monitor the parametric changes as the burst decays. Based on our time-resolved spectroscopy, we claim that the detected burst is a photospheric radius expansion (PRE) burst. During the PRE phase, the blackbody flux is found to be approximately constant at an averaged value $\sim$ 2.56 in $10^{-8}$ ergs s$^{-1}$ cm$^{-2}$ units. On the basis of literature survey, we infer that \textit{AstroSat}/LAXPC 20 has detected a burst from GX 3+1 after more than a decade which is also a PRE one. Utilising the burst parameters obtained, we provide a new estimation to the source distance, which is $\sim$ 9.3 $\pm$ 0.4 kpc, calculated for an isotropic burst emission. Finally, we discuss and compare our findings with the published literature reports.

We present results of the comparative analysis of the two wide binary systems -- 16 Cyg, with a giant gas planet orbiting around 16 Cyg B, and HD 219542 without planet detected. Atmospheric parameters of the binary components and the Sun were determined using their high-resolution spectra and the SME tools for automatic spectral analysis. By applying the synthetic spectrum method, we derived abundances of 29 and 23 chemical elements in 16 Cyg and HD 219542, respectively. For 19 of these elements, our results are based on the non-local thermodynamic equilibrium (NLTE) line formation. For both 16 Cyg and HD 219542, we obtained a small abundance difference between the A and B components: +0.019$\pm$0.012 and -0.014$\pm$0.019, respectively, suggesting only a weak influence of the giant gas planet formation on chemical composition of the host star atmosphere. For HD 219542 A and B, trends of the relative-to-solar abundances with the dust condensation temperature are similar to the literature data for the solar analogues without detected planets. The components of 16 Cyg reveal very similar behaviour of [X/H] with the condensation temperature, however, it is different from that for HD 219542. This indicates a specific chemical composition of the cloud from which the 16 Cyg binary system formed.

Cosmological probes at any redshift are necessary to reconstruct consistently the cosmic history. Studying properly the tension on the Hubble constant, $H_0$, obtained by Supernovae Type Ia (SNe Ia) and the Planck measurements of the Cosmic Microwave Background Radiation would require complete samples of distance indicators at any epoch. Gamma-Ray Bursts (GRBs) are necessary for the aforementioned task because of their huge luminosity that allows us to extend the cosmic ladder at very high redshifts. However, using GRBs alone as standard candles is challenging because their luminosity varies widely. To this end, we choose a reliable correlation for GRBs with a very small intrinsic scatter: the so-called fundamental plane correlation for GRB afterglows corrected for selection biases and redshift evolution. We choose a well-defined sample: the platinum sample, composed of 50 Long GRBs. To further constrain cosmological parameters, we use Baryon Acoustic Oscillations (BAOs) given their reliability as standard rulers. Thus, we have applied GRBs, SNe Ia, and BAOs in a binned analysis in redshifts so that GRBs' contribution is fully included in the last redshift bin, which reaches $z=5$. We use the fundamental plane correlation together with SNe Ia and BAOs, to constrain $H_0$ and the density matter today, $\Omega_{M}$. This methodology allows us to assess the role of GRBs combined with SNe Ia and BAOs. We have obtained results for $H_0$ and $\Omega_{M}$ using GRBs+ SNe Ia+BAOs with better precision than the SNe Ia alone for every bin, thus confirming the beneficial role of BAOs and GRBs added together. In addition, consistent results between GRBs+ SNe Ia +BAOs are obtained when compared with the SNe Ia +BAOs, showing the importance of GRBs since the distance ladder is extended up to $z=5$ with a similar precision obtained with other probes without including the GRBs.

Isolated black holes and neutron stars can be revealed through the observation of long duration gravitational microlensing events. A few candidates have been found in surveys of stars in the direction of the Galactic bulge. Recently, thanks to the addition of astrometric information at milliarcsecond level, it has been possible to reduce the uncertainties in the masses and distances for some of these "dark" gravitational lenses and select the most promising candidates. These isolated compact objects might emit X-rays powered by accretion from the interstellar medium. Using data of the Chandra, XMM-Newton, and INTEGRAL satellites, we searched for X-ray emission in the isolated black hole candidate OGLE-2011-BLG-0462, and in several other putative collapsed objects found with gravitational microlensing. OGLE-2011-BLG-0462 has been recently interpreted as a 7.1 M_sun black hole at a distance of 1.6 kpc, although a different group obtained a mass range (1.6-4.4 M_sun) that cannot exclude a massive neutron star. We have derived upper limits on the flux from OGLE-2011-BLG-0462 of 9$\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ in the 0.5-7 keV range and $\sim2\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$ in the 17-60 keV range. The implied X-ray luminosity is consistent with the small radiative efficiency expected for a black hole and disfavours a neutron star interpretation. Limits down to a factor about five lower are obtained for the soft X-ray flux of other candidates, but their interpretation is affected by larger uncertainties in the masses, distances and spatial velocities.

The number of type Ia supernova observations will see a significant growth within the next decade, especially thanks to the Legacy Survey of Space and Time undertaken by the Vera Rubin Observatory in Chile. With this rise, the statistical uncertainties will decrease and the flux calibration will become the main uncertainty for the characterization of dark energy. The uncertainty over the telescope transmission is a major systematic when measuring SNe Ia colors. Here we introduce the Collimated Beam Projector (CBP), a device that can measure the transmission of a telescope and its filters. Composed of a tunable monochromatic light source and optics to provide a parallel output beam this device is able to measure with high precision the filter transmissions. In the following, we will show how measuring precisely a telescope transmission can also improve the precision of the dark energy parameters. As an example, we will present the first results of the CBP in the context of the StarDice experiment.

Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the "Feedback In Realistic Environment" (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66% completeness, depending upon merger phase and the field-of-view of the observation. We detect fewer mergers in the pair phase (40%) and many more in the merging and post-coalescence phases (97%). We find that merger detectability decreases with field-of-view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combination of photometric and kinematic asymmetry, we increase these detections to 89% overall, 79% in pairs, and close to 100% in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields-of-view subtended by massively multiplexed integral field spectroscopy programmes.

Context. Gliese 832 (GJ 832) is an M2V star hosting a massive planet on a decade-long orbit, GJ 832b, discovered by radial velocity (RV). Later, a super Earth or mini-Neptune orbiting within the stellar habitable zone was reported (GJ 832c). The recently determined stellar rotation period (45.7 $\pm$ 9.3 days) is close to the orbital period of putative planet c (35.68 $\pm$ 0.03 days). Aims. We aim to confirm or dismiss the planetary nature of the RV signature attributed to GJ 832c, by adding 119 new RV data points, new photometric data, and an analysis of the spectroscopic stellar activity indicators. Additionally, we update the orbital parameters of the planetary system and search for additional signals. Methods. We performed a frequency content analysis of the RVs to search for periodic and stable signals. Radial velocity time series were modelled with Keplerians and Gaussian process (GP) regressions alongside activity indicators to subsequently compare them within a Bayesian framework. Results. We updated the stellar rotational period of GJ 832 from activity indicators, obtaining $37.5^{+1.4}_{-1.5}$ days, improving the precision by a factor of 6. The new photometric data are in agreement with this value. We detected an RV signal near 18 days (FAP < 4.6%), which is half of the stellar rotation period. Two Keplerians alone fail at modelling GJ 832b and a second planet with a 35-day orbital period. Moreover, the Bayesian evidence from the GP analysis of the RV data with simultaneous activity indices prefers a model without a second Keplerian, therefore negating the existence of planet c.

Only a handful of late-T brown dwarfs have been monitored for spectro-photometric variability, leaving incomplete the study of the atmospheric cloud structures of the coldest brown dwarfs, that share temperatures with some cold, directly-imaged exoplanets. 2MASSJ00501994-332240 is a T7.0 rapidly rotating, field brown dwarf that showed low-level photometric variability in data obtained with the Spitzer Space telescope. We monitored 2MASSJ00501994-332240 during ~2.6 hr with MOSFIRE, installed at the Keck I telescope, with the aim of constraining its near-infrared spectro-photometric variability. We measured fluctuations with a peak-to-peak amplitude of 1.48+\-0.75% in the J-band photometric light curve, an amplitude of 0.62+/-0.18% in the J-band spectro-photometric light curve, and an amplitude of 1.26+/-0.93% in the H-band light curve, and an amplitude of 5.33+/-2.02% in the CH_4-H_2O band light curve. Nevertheless, the Bayesian Information Criterion does not detect significant variability in any of the light curves. Thus, given the detection limitations due to the MOSFIRE sensitivity, we can only claim tentative low-level variability for 2M0050-3322 in the best-case scenario. The amplitudes of the peak-to-peak fluctuations measured for 2MASSJ00501994-332240 agree with the variability amplitude predictions of General Circulation Models for a T7.0 brown dwarf for an edge on object. Radiative-transfer models predict that the Na_2S and KCl clouds condense at pressures lower than that traced by the CH_4-H_2O band, which might explain the higher peak-to-peak fluctuations measured for this light curve. Finally, we provide a visual recreation of the map provided by General Circulation Models and the vertical structure of 2MASSJ00501994-332240 provided by radiative-transfer models.

IXPE (Imaging X-ray Polarimetry Explorer) is a NASA Small Explorer mission -- in partnership with the Italian Space Agency (ASI) -- dedicated to X-ray polarimetry in the 2--8 keV energy band. The IXPE telescope comprises three grazing incidence mirror modules coupled to three detector units hosting each one a Gas Pixel Detector (GPD), a gas detector that allows measuring the polarization degree by using the photoelectric effect. A wide and accurate ground calibration was carried out on the IXPE Detector Units (DUs) at INAF-IAPS, in Italy, where a dedicated facility was set-up at this aim. In this paper, we present the results obtained from this calibration campaign to study the IXPE focal plane detector response to polarized radiation. In particular, we report on the modulation factor, which is the main parameter to estimate the sensitivity of a polarimeter.

This paper investigates a new methodology to search for periods in light-curves of high-energy gamma-ray sources such as Active Galactic Nuclei (AGNs). High-energy light curves have significant stochastic components, making period detection somewhat challenging. In our model, periodic terms, drifts of the light-curves and random walk with correlation between flux points due to colored noise are taken into account independently. The parameters of the model are obtained directly from a Markov Chain Monte-Carlo minimization. The time periods found are compared to the output of the publicly available Agatha program. The search method is applied to high-energy periodic AGN candidates from the Fermi-LAT catalogue. The significance of periodic models over pure noise models is discussed. Finally, the variability of the period and amplitude of oscillating terms is studied on the most significant candidates.

The origin of the positron excess is one of the most intriguing mysteries in astroparticle physics. The recent discovery of extended $\gamma$-ray halos around the pulsars Geminga, Monogem and PSR J0621+3755 have brought indirect evidence that pulsar wind nebulae accelerate $e^{\pm}$ up to very-high-energy. While the precision of previous data does not permit precise evaluation of the parameters for the pulsars, we are able to find the more precise shape of the injection spectrum using new data released by HAWC and LHAASO in 2020 and 2021. We find that this is well fitten by a power-law with an exponential cutoff. The spectral index is quite hard with values around 1 while the cutoff energy is roughly 100 TeV. We also derive the strength of the diffusion coefficient around the pulsars finding that it is two orders of magnitude lower than the average of the Galaxy. Finally, we use the above mentioned results to estimate the contribution of Geminga to the positron excess. This source alone can contribute to the entire positron excess at around 1 TeV.

Galaxy protoclusters, which will eventually grow into the massive clusters we see in the local universe, are usually traced by locating overdensities of galaxies. Large spectroscopic surveys of distant galaxies now exist, but their sensitivity depends mainly on a galaxy's star formation activity and dust content rather than its mass. Tracers of massive protoclusters that do not rely on their galaxy constituents are therefore needed. Here we report observations of Lyman-$\alpha$ absorption in the spectra of a dense grid of background galaxies, which we use to locate a substantial number of candidate protoclusters at redshifts 2.2-2.8 via their intergalactic gas. We find that the structures producing the most absorption, most of which were previously unknown, contain surprisingly few galaxies compared to the dark matter content of their analogs in cosmological simulations. Nearly all are expected to be protoclusters, and we infer that half of their expected galaxy members are missing from our survey because they are unusually dim at rest-frame ultraviolet wavelengths. We attribute this to an unexpectedly strong and early influence of the protocluster environment on the evolution of these galaxies that reduced their star formation or increased their dust content.

We present a sample of 135,873 RR Lyrae stars (RRLs) with precise photometric metallicity and distance estimates from our newly calibrated $P$-$\phi_{31}$-[Fe/H] and $G$-band absolute magnitude-metallicity relations. The $P$-$\phi_{31}$-[Fe/H] relations for RRab and RRc are obtained from nearly 3000 {\it Gaia}-identified RRLs with precise $\phi_{31}$ measurements from the light curves and metallicity estimates from spectroscopy. Using over one thousand nearby RRLs with accurate distances estimated from the parallax measurements with Gaia EDR3, new $G$-band absolute magnitude metallicity relations and near-infrared period-absolute magnitude-metallicity relations for $K_{\rm s}$ and $W1$ bands are constructed. External checks, using other spectroscopic samples of field RRLs and RRL members of globular clusters, show that the typical uncertainties in photometric metallicity are 0.27/0.17 dex for RRab/RRc stars, respectively, without significant systematic offsets. The accuracies of these metallicity estimates are much improved, especially for RRab stars, when compared to those provided by the Gaia DR3. Validations of our distance estimates, again by using RRL members of globular clusters, show that the typical distance errors are 4.6%/3.4% for RRab/RRc stars, respectively. The distance scale from this study is consistent with that of the globular clusters. The distance modulus $\mu_{0}=18.444\pm0.135$ mag for the LMC and $\mu_{0}=18.940\pm0.147$ mag for the SMC are estimated from our RRab star sample, respectively, and are in excellent agreement with previous measurements. The mean metallicity of the LMC and SMC derived in this work are also consistent with the previous determinations. Using our sample, a steep metallicity gradient of $-0.021\pm0.001$ dex/kpc is found for the LMC, while a negligible metallicity gradient of $-0.005\pm0.003$ dex/kpc is obtained for the SMC.

The search for the primordial B-modes of the cosmic microwave background (CMB) relies on the separation from the brighter foreground dust signal. In this context, the characterisation of the spectral energy distribution (SED) of thermal dust in polarization has become a critical subject of study. We present a power-spectra analysis of Planck data, which improves on previous studies by using the newly released SRoll2 maps that correct residual data systematics, and by extending the analysis to regions near the Galactic plane. Our analysis focuses on the lowest multipoles between l=4 and 32, and three sky areas with sky fractions of fsky = 80%, 90%, and 97%. The mean dust SED for polarization and the 353 GHz Q and U maps are used to compute residual maps at 100, 143 and 217 GHz, highlighting spatial variations of the dust polarization SED. Residuals are detected at the three frequencies for the three sky areas. We show that models based on total intensity data are underestimating by a significant factor the complexity of dust polarized CMB foreground. Our analysis emphasizes the need to include variations of polarization angles of the dust polarized CMB foreground. The frequency dependence of the EE and BB power spectra of the residual maps yields further insight. We find that the moments expansion to the first order of the modified black-body (MBB) spectrum provides a good fit to the EE power- spectra. This result suggests that the residuals could follow mainly from variations of dust MBB spectral parameters. However, this conclusion is challenged by cross-spectra showing that the residuals maps at the three frequencies are not fully correlated, and the fact that the BB power-spectra do not match the first order moment expansion of a MBB SED. This work sets new requirements for simulations of the dust polarized foreground and component separation methods (abridged)

We introduce deep-field METACALIBRATION, a new technique that reduces the pixel noise in METACALIBRATION estimators of weak lensing shear signals by using a deeper imaging survey for calibration. In standard METACALIBRATION, when estimating the object's shear response, extra noise is added to correct the effect of shearing the noise in the image, increasing the uncertainty on shear estimates by ~ 20%. Our new deep-field METACALIBRATION technique leverages a separate, deeper imaging survey to calculate calibrations with less degradation in image noise. We demonstrate that weak lensing shear measurement with deep-field METACALIBRATION is unbiased up to second-order shear effects. We provide algorithms to apply this technique to imaging surveys and describe how to generalize it to shear estimators that rely explicitly on object detection (e.g., METACALIBRATION). For the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the improvement in weak lensing precision will depend on the somewhat unknown details of the LSST Deep Drilling Field (DDF) observations in terms of area and depth, the relative point-spread function properties of the DDF and main LSST surveys, and the relative contribution of pixel noise vs. intrinsic shape noise to the total shape noise in the survey. We conservatively estimate that the degradation in precision is reduced from 20% for METACALIBRATION to ~ 5% or less for deep-field METACALIBRATION, which we attribute primarily to the increased source density and reduced pixel noise contributions to the overall shape noise. Finally, we show that the technique is robust to sample variance in the LSST DDFs due to their large area, with the equivalent calibration error being ~ 0.1%. The deep-field METACALIBRATION technique provides higher signal-to-noise weak lensing measurements while still meeting the stringent systematic error requirements of future surveys.

A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, Earth-like planet discoveries are expected from the planned infrared astrometry space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars on the astrometric signal. This effect interferes with the astrometric detections of Earth-mass planets. While the first two papers considered stars rotating at the solar rotation rate, this paper focuses on stars having solar effective temperature and metallicity but rotating faster than the Sun, and consequently more active. By simulating the distribution of active regions on such stars using the Flux Emergence And Transport model, we show that the contribution of magnetic activity to the astrometric measurements becomes increasingly significant with increasing rotation rates. We further show that the jitter for the most variable periodic Kepler stars is high enough to be detected by Gaia. Furthermore, due to a decrease in the facula-to-spot area ratio for more active stars, the magnetic jitter is found to be spot-dominated for rapid rotators. Our simulations of the astrometric jitter has the potential to aid the interpretation of data from Gaia and upcoming space astrometry missions.

We present Hubble Space Telescope observations of active asteroid P/2020 O1 taken to examine its development for a year after perihelion. We find that the mass loss peaks <~1 kg/s in 2020 August and then declines to nearly zero over four months. Long-duration mass loss (~180 days) is consistent with a sublimation origin, indicating that this object is likely an ice-bearing main-belt comet. Equilibrium sublimation of water ice from an area as small as 1580 m^2 can supply the observed mass loss. Time-series photometry shows tentative evidence for extremely rapid rotation (double-peaked period < 2 hr) of the small nucleus (effective radius ~420 m). Ejection velocities of 0.1 mm particles are comparable to the 0.3 m/s gravitational escape speed from the nucleus, while larger particles are ejected at speeds less than the escape velocity. These properties are consistent with the sublimation of near-surface ice aided by centripetal acceleration. If water ice sublimation is confirmed, P/2020 O1 would be the icy asteroid with the smallest semimajor axis (highest temperature), setting new bounds on the distribution of ice in the asteroid belt.

Very recently, it was suggested that combining the Swampland program with the smallness of the dark energy and confronting these ideas to experiment lead to the prediction of the existence of a single extra-dimension (dubbed the dark dimension) with characteristic length-scale in the micron range. We show that the rate of Hawking radiation slows down for black holes perceiving the dark dimension and discuss the impact of our findings in assessing the dark matter fraction that could be composed of primordial black holes. We demonstrate that for a species scale of ${\cal O}(10^{10}~{\rm GeV})$, an all-dark-matter interpretation in terms of primordial black holes should be feasible for masses in the range $10^{15} \alt \mbh/{\rm g} \alt 10^{21}$. This range is extended compared to that in the 4D theory by 3 orders of magnitude in the low mass region.

We systematically study the top-down model of loop quantum black holes (LQBHs), recently derived by Alesci, Bahrami and Pranzetti (ABP). To understand the structure of the model, we first derive several well-known LQBH solutions by taking proper limits. These include the B\"ohmer-Vandersloot and Ashtekar-Olmedo-Singh models, which were all obtained by the so-called bottom-up polymerizations within the framework of the minisuperspace quantizations. Then, we study the ABP model, and find that the inverse volume corrections become important only when the radius of the two-sphere is of the Planck size. For macroscopic black holes, the minimal radius obtained at the transition surface is always much larger than the Planck scale, and hence these corrections are always sub-leading. The transition surface divides the whole spacetime into two regions, and in one of them the spacetime is asymptotically Schwarzschild-like, while in the other region, the asymptotical behavior sensitively depends on the ratio of two spin numbers involved in the model, and can be divided into three different classes. In one class, the spacetime in the 2-planes orthogonal to the two spheres is asymptotically flat, and in the second one it is not even conformally flat, while in the third one it can be asymptotically conformally flat by properly choosing the free parameters of the model. In the latter, it is asymptotically de Sitter. However, in any of these three classes, sharply in contrast to the models obtained by the bottom-up approach, the spacetime is already geodesically complete, and no additional extensions are needed in both sides of the transition surface. In particular, identical multiple black hole and white hole structures do not exist.

Solar flares not only pose risks to outer space technologies and astronauts' well being, but also cause disruptions on earth to our hight-tech, interconnected infrastructure our lives highly depend on. While a number of machine-learning methods have been proposed to improve flare prediction, none of them, to the best of our knowledge, have investigated the impact of outliers on the reliability and those models' performance. In this study, we investigate the impact of outliers in a multivariate time series benchmark dataset, namely SWAN-SF, on flare prediction models, and test our hypothesis. That is, there exist outliers in SWAN-SF, removal of which enhances the performance of the prediction models on unseen datasets. We employ Isolation Forest to detect the outliers among the weaker flare instances. Several experiments are carried out using a large range of contamination rates which determine the percentage of present outliers. We asses the quality of each dataset in terms of its actual contamination using TimeSeriesSVC. In our best finding, we achieve a 279% increase in True Skill Statistic and 68% increase in Heidke Skill Score. The results show that overall a significant improvement can be achieved to flare prediction if outliers are detected and removed properly.

The Maxwell equations imply that, under the background of non-zero $\boldsymbol{B}$, varying $\theta$ term produces $\boldsymbol{E} \cdot \boldsymbol{B}$. An interesting example is the Witten effect where a magnetic monopole becomes a dyon which, however, should disappear in the exact massless limit of the fermion. Underlying mechanism of this phenomenon has been understood by Callan by the presence of an effective axion-like degree of freedom around the monopole, which is roughly the phase of the fermions. The configuration of this axion cancels the effect of the $\theta$ term. Now, the chiral anomaly implies that non-vanishing $\boldsymbol{E} \cdot \boldsymbol{B}$ induces the chiral charge in the system. The question is whether the chiral charge is generated in the massless limit when we take into account the axion-like degree of freedom in the discussion. The discussion is relevant for the mechanism of baryogenesis under the background of time-dependent $\theta$. We solve the system of the massless QED with time dependent $\theta$ by reducing it to the two-dimensional QED. We demonstrate the occurrence of chiral charge generation in the background of static magnetic field for two cases: a magnetic monopole and a uniform magnetic flux. For the monopole case, the chiral charge comes out from the monopole while canceling the Witten effect. For the case of the uniform flux, on the other hand, the effect of the backreaction cannot be ignored, giving a more non-trivial time dependence. We also discuss their implications on baryogenesis.

The Saffman helicity integral of Hosking and Schekochihin (2021, PRX 11, 041005) has emerged as an important quantity that may govern the decay properties of magnetically dominated turbulence. Using a range of different computational methods, we confirm that this quantity is indeed gauge-invariant and nearly perfectly conserved in the limit of large Lundquist numbers. For direct numerical simulations with ordinary viscosity and magnetic diffusivity operators, we find that the solution develops in a nearly self-similar fashion. In a diagram quantifying the instantaneous decay coefficients of magnetic energy and integral scale, we find that the solution evolves along a line that is indeed suggestive of the governing role of the Saffman helicity integral. The solution settles near a line in this diagram that is expected for a self-similar evolution of the magnetic energy spectrum. However, the solution may settle in a slightly different position when the magnetic diffusivity decreases with time, which would be compatible with the decay being governed by the reconnection time scale rather than the Alfv\'en time.

We propose a new paradigm in the field of black hole physics, i.e., with more and more high precision observations, one must make a metric selection to determine which black hole is observationally preferred. In light of the shadow imaging data from the Event Horizon Telescope and gravitational wave measurements from the LIGO-VIRGO collaboration, we attempt to address this issue. Although not finding any obvious preference for a specific black hole metric based on current data, as a presentation of this new research paradigm, we rank these black hole metrics using the Bayesian information criterion. We find that Kerr and Reissner-Nordstr\"{o}m tie the first place in the black hole Olympics Games. Interestingly, we give the $2\,\sigma$ upper bound of the average electronic charge $Q<2.82\times10^{18}$ C for the Reissner-Nordstr\"{o}m spacetime, and the $2\,\sigma$ constraint on the average spin parameter $a=-0.02^{+4.07}_{-4.04}$ m for the Kerr spacetime, which is consistent with the prediction of zero rotation for a distant observer. Future multi-messenger and multi-wavelength observations will enhance the application of this new paradigm into tests of gravitational theories and black holes.

Finsler geometry is a natural and fundamental generalization of Riemann geometry. The Finsler structure depends on both coordinates and velocities. We present the arrival time delay of astroparticles subject to Lorentz violation in the framework of Finsler geometry, and the result corresponds to that derived by Jacob and Piran in the standard model of cosmology.