Papers for Friday, Nov 05 2021

We review two initial orbit determination methods for too short arcs (TSAs) of optical observations of a solar system body. These methods employ the conservation laws of Kepler's problem, and allow to attempt the linkage of TSAs referring to quite far epochs, differing by even more than one orbital period of the observed object. The first method ({\tt Link2}) concerns the linkage of 2 TSAs, and leads to a univariate polynomial equation of degree 9. An optimal property of this polynomial is proved using Gr\"obner bases theory. The second method ({\tt Link3}) is thought for the linkage of 3 TSAs, and leads to a univariate polynomial equation of degree 8. A numerical test is shown for both algorithms.

Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

Population III stars were the first generation of stars, formed in minihalos of roughly primordial element abundances, and therefore metal-free. They are thought to have formed at the cores of dense dark matter clouds. Interactions between baryons and dark matter can therefore have had an important impact on their evolution. In this paper we consider the capture of non- or weakly-annihilating dark matter by these early massive stars. In a wide region of parameter space, interactions of dark matter with baryons lead to premature death of the star as a black hole. We sketch how this modification of the standard evolutionary history of Population III stars might impact the epoch of reionisation, by modifying the amount of UV emission, the transition to Population II star formation, and the X-ray and radio emission from accretion onto the black hole remnants. Signals of massive black holes originating from Population III stars could be observed through gravitational waves from their mergers. Finally, the observation of pair-instability supernovae could effectively preclude premature black hole death across a wide range of parameter space, ranging in mass from $m_{DM} \sim 0.1\text{ GeV}$ to $m_{DM} \sim m_{\rm Pl}$.

We present the first results from Citizen ASAS-SN, a citizen science project for the All-Sky Automated Survey for Supernovae (ASAS-SN) hosted on the Zooniverse platform. Citizen ASAS-SN utilizes the newer, deeper, higher cadence ASAS-SN $g$-band data and tasks volunteers to classify periodic variable star candidates based on their phased light curves. We started from 40,640 new variable candidates from an input list of ${\sim} 7.4$ million stars with $\delta < -60^\circ$ and the volunteers identified 10,420 new discoveries which they classified as 4,234 pulsating variables, 3,132 rotational variables, 2,923 eclipsing binaries, and 131 variables flagged as Unknown. They classified known variable stars with an accuracy of 89% for pulsating variables, 81% for eclipsing binaries, and 49% for rotational variables. We examine user performance, agreement between users, and compare the citizen science classifications with our machine learning classifier updated for the $g$-band light curves. In general, user activity correlates with higher classification accuracy and higher user agreement. We used the user's "Junk" classifications to develop an effective machine learning classifier to separate real from false variables, and there is a clear path for using this "Junk" training set to significantly improve our primary machine learning classifier. We also illustrate the value of Citizen ASAS-SN for identifying unusual variables with several examples.

We analyse some open debates in cosmology in the light of the most updated quasar (QSO) sample, covering a redshift range up to $z\sim7.5$, combined with type Ia supernovae (SNe) and baryon acoustic oscillations (BAO) data. Indeed, extending the cosmological analyses with high-redshift data is the key to distinguishing between different cosmological models and allowing a better constraint on dark energy (DE) evolution. Also, we discuss different combinations of the BAO, SNe, and QSO data to understand their compatibility and their implications for the non-flat Universe and extensions of the standard cosmological model. Specifically, we consider a flat and non-flat $\Lambda \mathrm{CDM}$ cosmology, a flat and non-flat DE model with a constant DE EoS parameter ($w$), and four flat DE models with variable $w$: Chevallier-Polarski-Linder, Jassal-Bagla-Padmanabhan, an "exponential" and "rational" parameterisations. We find that a joint analysis of QSO+SNe with BAO is only possible in the context of a flat Universe. Indeed BAO confirms the flatness condition assuming a curved geometry, whilst SNe+QSO show evidence of a closed space. We also find that the matter component is fully consistent with $\Omega_{M,0}=0.3$ in all the data sets assuming a flat $\Lambda \mathrm{CDM}$ model. Yet, all the other analysed models show a statistically significant deviation at $>3\sigma$ from this prediction by making use of the SNe+QSO sample, which remains still significant (2-3$\sigma$) with the combined SNe+QSO+BAO data set. In the models where the DE density evolves with time, the SNe+QSO+BAO data always prefer $\Omega_{M,0}>0.3$, $w_{0}<-1$ and $w_{a}$ greater, but statistically consistent, than $w_{a}=0$. This DE phantom behaviour is mainly driven by the contribution of SNe+QSO, while BAO are closer to the prediction of the flat $\Lambda \mathrm{CDM}$ model (also due to their large uncertainties).

To extract information from the clustering of galaxies on non-linear scales, we need to model the connection between galaxies and halos accurately and in a flexible manner. Standard halo occupation distribution (HOD) models make the assumption that the galaxy occupation in a halo is a function of only its mass, however, in reality, the occupation can depend on various other parameters including halo concentration, assembly history, environment, spin, etc. Using the IllustrisTNG hydrodynamic simulation as our target, we show that machine learning tools can be used to capture this high-dimensional dependence and provide more accurate galaxy occupation models. Specifically, we use a random forest regressor to identify which secondary halo parameters best model the galaxy-halo connection and symbolic regression to augment the standard HOD model with simple equations capturing the dependence on those parameters, namely the local environmental overdensity and shear, at the location of a halo. This not only provides insights into the galaxy-formation relationship but, more importantly, improves the clustering statistics of the modeled galaxies significantly. Our approach demonstrates that machine learning tools can help us better understand and model the galaxy-halo connection, and are therefore useful for galaxy formation and cosmology studies from upcoming galaxy surveys.

The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.

Many massive early-type galaxies host central radio sources and hot X-ray atmospheres indicating the presence of radio-mechanical active galactic nucleus (AGN) feedback. The duty cycle and detailed physics of the radio-mode AGN feedback is still a matter of debate. To address these questions, we present 1-2 GHz Karl G. Jansky Very Large Array (VLA) radio observations of a sample of the 42 nearest optically and X-ray brightest early-type galaxies. We detect radio emission in 41/42 galaxies. However, the galaxy without a radio source, NGC 499, has recently been detected at lower frequencies by the Low-Frequency Array (LOFAR). Furthermore, 27/42 galaxies in our sample host extended radio structures and 34/42 sources show environmental interactions in the form of X-ray cavities. We find a significant correlation between the radio flux density and the largest linear size of the radio emission and between the radio power and the luminosity of the central X-ray point-source. The central radio spectral indices of the galaxies span a wide range of values, with the majority of the systems having steep spectra and the rest flat spectra. These results are consistent with AGN activity, where the central radio sources are mostly switched on, thus the duty cycle is very high. 7/14 galaxies with point-like radio emission (Fanaroff-Riley Class 0; FR 0) also show X-ray cavities indicating that, despite the lack of extended radio structures at 1-2 GHz, these AGN do launch jets capable of inflating lobes and cavities.

We discuss new consistency relations for single field models of inflation capable of generating primordial black holes (PBH), and their observational implications for CMB $\mu$-space distortions. These inflationary models include a short period of non-attractor evolution: the scale-dependent profile of curvature perturbation is characterized by a pronounced dip, followed by a rapid growth leading to a peak responsible for PBH formation. We investigate the squeezed and the collapsed limits of three and four point functions of curvature perturbation around the dip, showing that they satisfy consistency relations connecting their values to the total amplification of the curvature spectrum, and to the duration of the non-attractor era. Moreover, the corresponding non-Gaussian parameters are scale-dependent in proximity of the dip, with features that again depend on the amplification of the spectrum. For typical PBH scenarios requiring an order ${\cal O}(10^7)$ enhancement of the spectrum from large towards small scales, we generally find values $f_{\rm NL}^{\rm sq}\,=\,{\cal O}(10)$ and $\tau_{\rm NL}^{\rm col}\,=\,{\cal O}(10^3)$ in a range of scales that can be probed by CMB $\mu$-space distortions. Using these consistency relations, we carefully analyze how the scale-dependence of non-Gaussian parameters leads to characteristic features in $\langle \mu T \rangle$ and $\langle \mu \mu \rangle$ correlators, providing distinctive probes of inflationary PBH scenarios that can be tested using well-understood CMB physics.

Most galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic flux density variations in the blazar PKS 2131-021, which strongly suggest a SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our well-sampled, 45.1-year, 14.5-15.5 GHz, light curve shows two epochs of strong sinusoidal variation with the same period to within <2%, and the same phase to within $\sim 10\%$, straddling a 22-year period when this periodic variation was absent. We generated simulated light curves that accurately reproduce the "red noise" character of the radio light curve of this object, and carried out Lomb-Scargle, weighted wavelet Z transform, and least-squares sine wave analyses that demonstrate conclusively, at the $4.6\sigma$ significance level, that the periodicity in this object is not due to random fluctuations, but to a physical periodicity in the source. The observed period measured across the whole 45.1-year span of our observations is $4.758\pm 0.007$ years (i.e., $\delta P/P \sim 1.5 \times 10^{-3}$), translating to $2.082\pm 0.003$ years in the rest frame of PKS 2131-021. While periodic variations are expected from a SMBHB, the periodic variation in PKS 2131-021 is remarkably sinusoidal, which should provide important insights into this system. We present a model for the putative SMBHB, in which the orbital motion, when combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily accounts for the amplitude of the observed modulation. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS 2131-021 and searches using the techniques we have developed are strongly motivated.

Magnetically arrested accretion discs (MADs) around black holes (BH) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here we propose instead that it is mainly a dynamically important toroidal field -- the result of dynamo action triggered by the significant but still relatively weak vertical field -- that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotatonal instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (strongly magnetized but not "arrested") discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.

Know thy star, know thy planet,... especially in the ultraviolet (UV). Over the past decade, that motto has grown from mere wish to necessity in the M dwarf regime, given that the intense and highly variable UV radiation from these stars is suspected of strongly impacting their planets' habitability and atmospheric loss. This has led to the development of the Star-Planet Activity Research CubeSat (SPARCS), a NASA-funded 6U CubeSat observatory fully devoted to the photometric monitoring of the UV flaring of M dwarfs hosting potentially habitable planets. The SPARCS science imaging system uses a 9-cm telescope that feeds two delta-doped UV-optimized CCDs through a dichroic beam splitter, enabling simultaneous monitoring of a target field in the near-UV and far-UV. A dedicated onboard payload processor manages science observations and performs near-real time image processing to sustain an autonomous dynamic exposure control algorithm needed to mitigate pixel saturation during flaring events. The mission is currently half-way into its development phase. We present an overview of the mission's science drivers and its expected contribution to our understanding of star-planet interactions. We also present the expected performance of the autonomous dynamic exposure control algorithm, a first-of-its-kind on board a space-based stellar astrophysics observatory.

A multiwavelength observing campaign of the T Tauri star (TTS) GM Aur was undertaken in 2019 December that obtained Swift X-ray and NUV fluxes, HST NUV spectra, LCOGT ugri and TESS photometry, CHIRON Halpha spectra, ALMA 13CO and C18O line fluxes, and VLA 3 cm continuum fluxes taken contemporaneously over one month. The X-ray to optical observations were presented previously. Here we present the ALMA and VLA data and make comparisons to GM Aur's accretion and X-ray properties. We report no variability in the observed millimeter CO emission. Using disk chemistry models, we show that the magnitude of the changes seen in the FUV luminosity of GM Aur could lead to variation of up to ~6% in CO line emission and changes in the X-ray luminosity could lead to larger changes of ~25%. However, the FUV and X-ray luminosity increases must last at least 100 years in order to induce changes, which seems implausible in the TTS stage; also, these changes would be too small to be detectable by ALMA. We report no variability in the 3 cm emission observed by the VLA, showing that changes of less than a factor of ~3 in the accretion rates of TTSs do not lead to detectable changes in the mass-loss rate traced by the jet at centimeter wavelengths. We conclude that typically seen changes in the FUV and X-ray luminosities of TTSs do not lead to observable changes in millimeter CO line emission or jet centimeter continuum emission.

Context: The Nova V1405 Cas, which erupted on 18 March 2021, showed a pre-maximum phase of almost two months. After the main maximum, during several month the nova only weakened slightly and showed clear fluctuations in brightness. Aims: The aim of this work is to investigate the relationship between brightness and half-width (FWHM) and equivalent width (EW) of the H{\alpha} emission line, and from this to infer possible causes for the abnormal behavior of the Nova. Methods: Magnitudes and spectra from online databases were used, the H{\alpha} line of the spectra was measured and 2-day normal values were formed in order to establish a temporal comparison. Results: After the main maximum, during the period investigated by JD 2459292 to 2459305 at intervals of 2-5 weeks secondary maxima with different characteristics. The maxima of the equivalent width were 6-14 days after the respective brightness maxima. The maxima are interpreted as subsequent outbreaks of the nova. Possibly these are subsequent weaker eruptions of the nova or other mechanisms that lead to an increase in density in the envelope.

The current experimentally measured parameters of the Standard Model (SM) suggest that our Universe lies in a metastable electroweak vacuum, where the Higgs field is prone to vacuum decay to a lower state with catastrophic consequences. Our measurements dictate that such an event has not taken place yet, despite the many different mechanisms that could have triggered it in our past light-cone. The focus of our work has been to calculate the probability of the false vacuum to decay during the period of inflation and use it to constrain the last unknown renormalisable SM parameter $\xi$, which couples the Higgs field with space-time curvature. More specifically, we derived lower $\xi$-bounds from vacuum stability in three inflationary models: quadratic and quartic chaotic inflation, and Starobinsky-like power-law inflation. We also took the time-dependence of the Hubble rate into account both in the geometry of our past light-cone and in the Higgs effective potential, which is approximated with three-loop renormalisation group improvement supplemented with one-loop curvature corrections.

Globular clusters are prone to lose stars while moving around the Milky Way. These stars escape the clusters and are distributed throughout extended envelopes or tidal tails. However, such extra-tidal structures are not observed in all globular clusters, and yet there is no structural or dynamical parameters that can predict their presence or absence. NGC\,6864 is an outer halo globular cluster with reported no observed tidal tails. We used Dark Energy Camera (DECam) photometry reaching $\sim$ 4 mag underneath its main sequence turnoff to confidently detect an extra-tidal envelope, and stellar debris spread across the cluster outskirts. These features emerged once robust field star filtering techniques were applied to the fainter end of the observed cluster main sequence. NGC\,6864 is associated to the {\it Gaia}-Enceladus dwarf galaxy, among others 28 globular clusters. Up-to-date, nearly 64$\%$ of them have been targeted looking for tidal tails and most of them have been confirmed to exhibit tidal tails. Thus, the present outcomes allow us to speculate on the possibility that {\it Gaia}-Enceladus globular clusters share a common pattern of mass loss by tidal disruption.

Our current measurements of the Standard Model parameters imply that the Higgs field resides in a metastable electroweak vacuum, where the vacuum can decay to a lower ground state, with cataclysmic repercussions for our Universe. According to our observations, no such event has happened in the observable universe, in spite of the various energetic processes that could have triggered it. This work serves as an overview of a method that uses the metastability of the false vacuum during cosmological inflation to provide constraints on the Higgs curvature coupling $\xi$. Considering also the effects of the time-dependent Hubble rate on the effective Higgs potential and on our past light-cone space-time geometry, results in state-of-the-art lower $\xi$-bounds from quadratic and quartic chaotic inflation, and Starobinsky-like power-law inflation.

M87 has been the target of numerous astronomical observations across the electromagnetic spectrum and Very Long Baseline Interferometry (VLBI) resolved an edge-brightened jet. However, the origin and formation of its jets remain unclear. In our current understand black holes (BH) are the driving engine of jet formation, and indeed the recent Event Horizon Telescope (EHT) observations revealed a ring-like structure in agreement with theoretical models of accretion onto a rotating Kerr BH. In addition to the spin of the BH being a potential source of energy for the launching mechanism, magnetic fields are believed to play a key role in the formation of relativistic jets. A priori, the spin, $a_\star$, of BH in M87* is unknown, however, when accounting for the estimates on the X-ray luminosity and jet power, values $\left |a_\star \right| \gtrsim 0.5$ appear favoured. Besides the properties of the accretion flow and the BH spin, the radiation microphysics including the particle distribution (thermal and non-thermal) as well as the particle acceleration mechanism play a crucial role. We show that general-relativistic magnetohydrodynamics simulations and general-relativistic radiative transfer calculations can reproduce the broadband spectrum from the radio to the near-infrared regime and simultaneously match the observed collimation profile of M87, thus allowing us to set rough constraints on the dimensionless spin of M87* to be $0.5\lesssim a_{\star}\lesssim 1.0$, with higher spins being possibly favoured.

The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal emission. In order to bridge the gap between these scales and to provide a theoretical interpretation of these observations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching. M87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique details on the collimation profile of the jet down to several gravitational radii. In this work we aim to model the observed broad-band spectrum of M87 from the radio to the NIR regime and at the same time fit the jet structure as observed with Global mm-VLBI at 86 GHz. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The radiative signatures of these simulations are computed taking different electron distribution functions into account and a detailed parameter survey is performed in order to match the observations. The results of our simulations show that magnetically arrested disks around fast spinning black holes ($a_\star\geq0.5$) together with a mixture of thermal and non-thermal particle distributions are able to model simultaneously the broad-band spectrum and the innermost jet structure of M87

One of the main objectives of the CREDO project is to search for so-called Cosmic-Ray Ensembles (CRE) \cite{homola2020cosmic}. To confirm the existence of such phenomena a massive scale observation of even relatively low energy Extensive Air Showers (EAS) and an analysis of their correlations in time must be performed. To make such observations possible, an infrastructure of widely spread detectors connected in a global network should be developed using low-cost devices capable of collecting data for a long period of time. For each of these detectors or small detector systems the probability of detection of an EAS has to be determined. Such information is crucial in the analysis and interpretation of collected data. In the case of large number of systems with different properties the standard approach based on detailed and extensive simulations is not possible, thus a faster method is developed. Knowing the characteristics of EAS from more general simulations any required probability is calculated using appropriate parameterization taking into account EAS spectrum, energy dependence of particle density and zenith angle dependence. This allows to estimate expected number of EAS events measured by a set of small detectors \cite{Karbowiak_2020} and compare results of calculations with these measurements.

We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 \lambda$. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularised expectation-maximisation algorithm with consensus optimisation (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data.

Ram-pressure stripping by the intracluster medium (ICM) is one of the most advocated mechanisms that affect the properties of cluster galaxies. A recent study based on a small sample has found that many galaxies showing strong signatures of ram-pressure stripping also possess an active galactic nucleus (AGN), suggesting a possible correlation between the two phenomena. This result has not been confirmed by a subsequent study. Building upon previous findings, here we combine MUSE observations conducted within the GASP program and a general survey of the literature to robustly measure the AGN fraction in cluster's ram pressure stripped galaxies using BPT emission line diagrams. Considering a sample of 115 ram pressure stripped galaxies with stellar masses $\geq 10^9 \, M_{\odot}$, we find an AGN fraction of $\sim27\%$. This fraction strongly depends on stellar mass: it raises to 51\% when only ram-pressure stripped galaxies of masses $M_* \geq 10^{10} \, M_{\odot}$ are considered. We then investigate whether the AGN incidence is in excess in ram pressure stripped galaxies compared to non-stripped galaxies, using as comparison a sample of non-cluster galaxies observed by the survey MaNGA. Considering mass-matched samples, we find that the incidence of AGN activity is significantly higher (at a confidence level $>99.95\%$) when ram-pressure stripping is on act, supporting the hypothesis of an AGN-ram pressure connection.

We study the host galaxy properties of active galactic nuclei (AGN) that have been detected in X-rays in the nearby Universe ($\rm z<0.2$). For that purpose, we use the catalogue provided by the ROSAT-2RXS in the 0.1-2.4\,keV energy band, one of the largest X-ray datasets with spectroscopic observations. Our sample consists of $\sim 900$ X-ray AGN. The catalogue provides classification of the sources into type 1 and 2, based on optical spectra. $\sim 25\%$ of the AGN are type 2. We use the available optical, near-IR and mid-IR photometry to construct SEDs. We measure the stellar mass ($M_*$) and star formation rate (SFR) of the AGN, by fitting these SEDs with the X-CIGALE code. We compare the $M_*$ and SFR of the two AGN populations, taking into account their different redshift and luminosity distributions. Based on our results, type 2 AGN tend to live in more massive galaxies compared to their type 1 counterparts ($\rm log\,[M_*(M_\odot)]=10.49^{+0.16}_{-0.10}$ vs. $10.23^{+0.05}_{-0.08}$), in agreement with previous studies at higher redshifts. In terms of SFR, our analysis shows, that in the nearby Universe, the number of X-ray AGN that live in quiescent systems is increased compared to that at higher redshifts, in accordance with previous studies in the local universe. However, the majority of AGN ($\sim 75\%$) live inside or above the main sequence.

We report on gamma-ray orbital modulation of the transitioning millisecond pulsar binary XSS J12270-4859 detected in the Fermi Large Area Telescope (LAT) data. We use long-term optical data taken with the XMM-Newton Optical Monitor and the Swift UltraViolet Optical Telescope to inspect radio timing solutions that are limited to relatively short time intervals, and find that extrapolation of the solutions aligns well the phasing of the optical data over 15 years. The Fermi-LAT data folded on the timing solutions exhibit significant modulation (p = 5x10^-6) with a gamma-ray minimum at the inferior conjunction of the pulsar. Intriguingly, the source seems to show similar modulation in both the low-mass X-ray binary (LMXB) and the millisecond pulsar (MSP) states, implying that mechanisms for gamma-ray emission in the two states are similar. We discuss these findings and their implications using an intra-binary shock scenario.

It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multi-wavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multi-wavelength flares. Nevertheless, previous works have not directly evaluated $\gamma$-ray signatures from first-principle simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multi-wavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and {\it Fermi-LAT} $\gamma$-ray bands as well as closely correlated optical and $\gamma$-ray flares. The optical polarization angle swings are also accompanied by $\gamma$-ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast $\gamma$-ray flares or orphan $\gamma$-ray flares under the leptonic scenario, complementary to the frequently considered mini-jet scenario. It may also infer neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

We present mass and radius measurements of K2-79b and K2-222b, two transiting exoplanets orbiting active G-type stars. Their respective 10.99d and 15.39d orbital periods fall near periods of signals induced by stellar magnetic activity. The two signals might therefore interfere and lead to an inaccurate estimate of exoplanet mass. We present a method to mitigate these effects when radial velocity and activity indicator observations are available over multiple observing seasons and the orbital period of the exoplanet is known. We perform correlation and periodogram analyses on sub-sets composed of each target's two observing seasons, in addition to the full data sets. For both targets, these analyses reveal an optimal season with little to no interference at the orbital period of the known exoplanet. We make a confident mass detection of each exoplanet by confirming agreement between fits to the full radial velocity set and the optimal season. For K2-79b, we measure a mass of 11.8 $\pm$ 3.6 $M_{Earth}$ and a radius of 4.09 $\pm$ 0.17 $R_{Earth}$. For K2-222b, we measure a mass of 8.0 $\pm$ 1.8 $M_{Earth}$ and a radius of 2.35 $\pm$ 0.08 $R_{Earth}$. According to model predictions, K2-79b is a highly irradiated Uranus-analog and K2-222b hosts significant amounts of water ice. We also present an RV solution for a candidate second companion orbiting K2-222 at 147.5d.

We present a model for the broadband radio-to-$\gamma$-ray spectral energy distribution of the compact radio source, PKS 1718-649. Because of its young age (100 years) and proximity ($z=0.014$), PKS 1718-649 offers a unique opportunity to study nuclear conditions and the jet/host galaxy feedback process at the time of an initial radio jet expansion. PKS 1718-649 is one of a handful of young radio jets with $\gamma$-ray emission confirmed with the Fermi/LAT detector. We show that this $\gamma$-ray emission can be successfully explained by Inverse Compton scattering of the ultraviolet photons, presumably from an accretion flow, off non-thermal electrons in the expanding radio lobes. The origin of the X-ray emission in PKS 1718-649 is more elusive. While Inverse Compton scattering of the infrared photons emitted by a cold gas in the vicinity of the expanding radio lobes contributes significantly to the X-ray band, the data require that an additional X-ray emission mechanism is at work, e.g. a weak X-ray corona or a radiatively inefficient accretion flow, expected from a LINER type nucleus such as that of PKS 1718-649. We find that the jet in PKS 1718-649 has low power, $L_j \simeq 2.2 \times 10^{42}$ erg s$^{-1}$, and expands in an environment with density $n_0 \simeq 20$ cm$^{-3}$. The inferred mass accretion rate and gas mass reservoir within 50-100 pc are consistent with estimates from the literature obtained by tracing molecular gas in the innermost region of the host galaxy with SINFONI and ALMA.

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SNe) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion amongst metal poor stars, which suggests that the interstellar medium (ISM) was not very well mixed at these early times. The purpose of this Letter is to help isolate via a series of numerical experiments some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at parsec scales. By investigating the statistics of variations of $\alpha$ elements in these simulations we are able to derive meaningful constraints on the star formation history (SFH) of the Milky Way (MW). We argue that the observed dispersion in stellar abundances of metal poor halo stars are naturally explained by star formation conditions expected in dwarf satellites or in an early low starforming MW progenitor. Metal variations in phase-mixed stars give further credence to the idea that the MW halo could have been assembled primarily by disrupted dwarf satellites.

New, high-resolution MARCS synthetic spectra have been calculated for more than a dozen mixtures of the metals allowing, in turn, for variations in C:N:O, [CNO/Fe], and enhanced abundances of C, O, Mg, and Si. Bolometric Corrections (BCs) for many of the broad-band filters currently in use have been generated from these spectra. Due to improved treatments of molecules that involve atoms of C, N, and O, the BCs for UV and blue passbands, in particular, differ substantially from those derived from previous MARCS models. These differences, and the effects on BCs of varying the abundances of the metals, are shown in a number of instructive plots. Stellar evolutionary grids for -2.5 <= [Fe/H] <= -0.5 have also been computed for the different mixtures. Isochrones based on these tracks are intercompared on the theoretical H-R diagram and on a few of the colour-magnitude diagrams that can be constructed from HST Wide Field Camera 3 (WFC3) F336W, F438W, F606W, F814W, F110W, and F160W observations. For the first application of these models, isochrones have been fitted to WFC3 photometry of the globular cluster NGC 6496 from the HST UV Legacy Survey, with very encouraging results.

Stellar models for -2.5 <= [Fe/H] <= -0.5 that have been computed for variations in the C:N:O abundance ratio (for two different values of [CNO/Fe]) are compared with HST Wide Field Camera 3 (WFC3) observations of the globular clusters (GCs) 47 Tuc, NGC 6362, M5, M5, M55, and M92. The bolometric corrections used to transpose the models to the observed planes are based on new MARCS synthetic spectra that incorporate improved treatments of molecules that involve atoms of C, N, and O. On the assumption of well-supported distance moduli and reddenings, isochrones for [O/Fe] = 0.6 and [m/Fe] = 0.4 for the other alpha elements, which are favoured by binary stars in GCs, generally reproduce the main features of observed colour-magnitude diagrams (CMDs) to within ~ 0.03 mag. In particular, they appear to match the spreads in the observed (M_F336W - M_F438W)_0 colours that are spanned by CN-weak and CN-strong stars along the lower giant branch quite well, but not the bluest giants, which are suspected to be N-poor ([N/Fe] <~ -0.5). Both the absolute (M_F438W - M_F606W)_0 colours and the variations in these colours at a given M_F606W magnitude on the giant branch are difficult to explain unless the reddest stars are C-rich ([C/Fe] >~ +0.5). Allowing for moderate He abundance variations (delta Y ~ 0.05) improves the fits to the observations.

The stellar initial mass function (IMF) is expressed by $\phi(m) \propto m^{-\alpha}$ with the slope $\alpha$, and known as the poorly-constrained but very important function in studies of star and galaxy formation. There are no sensible observational constraints on the IMF slopes beyond Milky Way and nearby galaxies. Here we combine two sets of observational results, 1) cosmic densities of core-collapse supernova explosion (CCSNe) rates and 2) cosmic far ultraviolet radiation (and infrared re-radiation) densities, which are sensitive to massive ($\simeq 8-50 \,{\rm M}_\odot$) and moderately massive ($\simeq 2.5-7 \,{\rm M}_\odot$) stars, respectively, and constrain the IMF slope at $m>1\,{\rm M}_\odot$ with a freedom of redshift evolution. Although no redshift evolution is identified beyond the uncertainties, we find that the cosmic average IMF slope at $z=0$ is $\alpha=1.8-3.2$ at the 95\% confidence level that is comparable with the Salpeter IMF, $\alpha=2.35$, which marks the first constraint on the cosmic average IMF. We show a forecast for the Nancy Grace Roman Space Telescope supernova survey that will provide significantly strong constraints on the IMF slope with $\delta \alpha\simeq 0.5$ over $z=0-2$. Moreover, as for an independent IMF probe instead of 1), we suggest to use diffuse supernovae neutrino background (DSNB), relic neutrinos from CCSNe. We expect that the Hyper-Kamiokande neutrino observations over 20 years will improve the constraints on the IMF slope and the redshift evolution significantly better than those obtained today, if the systematic uncertainties of DSNB production physics are reduced in the future numerical simulations.

Supermassive black holes and supernovae explosions at the centres of active galaxies power cycles of outflowing and inflowing gas that affect galactic evolution and the overall structure of the Universe. While simulations and observations show that this must be the case, the range of physical scales (over ten orders of magnitude) and paucity of available tracers, make both the simulation and observation of these effects difficult. By serendipity, there lies an active galaxy, Centaurus A (NGC 5128), at such a close proximity as to allow its observation over this entire range of scales and across the entire electromagnetic spectrum. In the radio band, however, details on scales of 10-100 kpc from the supermassive black hole have so far been obscured by instrumental limitations. Here we report low-frequency radio observations that overcome these limitations and show evidence for a broad, bipolar outflow with velocity 1100 km per s and mass outflow rate of 2.9 solar masses per year on these scales. We combine our data with the plethora of multi-scale, multi-wavelength historical observations of Centaurus A to probe a unified view of feeding and feedback, which we show to be consistent with the Chaotic Cold Accretion self-regulation scenario.

Current cosmological observations point to a serious discrepancy between the observed Hubble parameter obtained using direct and cosmic microwave background radiation (CMBR) measurements. Besides this, the so called Hubble--Lema\^itre tension, we also find considerable evidence in diverse cosmological observables that indicate violation of the cosmological principle. In this paper, we suggest that both these discrepancies are related and can be explained by invoking superhorizon perturbations in the Universe. We implement this by considering a single superhorizon mode and showing that it leads to both a dipole in large scale structures and a shift in the Hubble--Lema\^itre parameter. Furthermore, the shift is found to be independent of redshift up to a certain distance. This is nicely consistent with the data.

The amplitude of the primordial magnetic field (PMF) is constrained from observational limits on primordial nuclear abundances. Within this constraint, it is possible that nuclear motion is regulated by Coulomb scattering with electrons and positrons ($e^\pm$'s), while $e^\pm$'s are affected by a PMF rather than collisions. For example, at a temperature of $10^9$ K, thermal nuclei typically experience $\sim 10^{21}$ scatterings per second that are dominated by very small angle scattering leading to minuscule changes in the nuclear kinetic energy of order $\mathcal{O}$(1) eV. In this paper the upper limit on the effects of a possible discretization of the $e^\pm$ momenta by the PMF on the nuclear momentum distribution is estimated under the extreme assumptions that the momentum of the $e^\pm$ is relaxed before and after Coulomb scattering to Landau levels, and that during Coulomb scattering the PMF is neglected. This assumption explicitly breaks the time reversal invariance of Coulomb scattering, and the Maxwell-Boltzmann distribution is not a trivial steady state solution of the Boltzmann equation under these assumptions. We numerically evaluate the collision terms in the Boltzmann equation, and show that the introduction of a special direction in the $e^\pm$ distribution by the PMF generates no directional dependence of the collisional destruction term of nuclei. Large anisotropies in the nuclear distribution function are then constrained from big bang nucleosynthesis. Ultimately, we conclude that a PMF does not significantly affect the isotropy or BBN.

The study of spiral structures in protoplanetary disks is of great importance for understanding of processes in the disks, including planet formation. Bright spiral arms were detected in the disk of young star CQ Tau by Uyama et al. (2020) in H and L bands. The spiral arms are located inside the gap in millimeter size dust, recovered earlier using ALMA observations (Ubeira Gabellini et al. 2019). To explain the gap, Ubeira Gabellini et al. (2019) proposed the existence of a planet with the semimajor axis of 20 AU. We obtained multiepoch observations of a spiral feature in the circumstellar envelope of CQ Tau in Ic band using a novel technique of differential speckle polarimetry. The observations covering period from 2015 to 2021 allow us to estimate the pattern speed of spiral: $-0.2\pm1.1^{\circ}$/yr (68% credible interval, positive value indicates counter-clockwise rotation), assuming face-on orientation of the disk. This speed is significantly smaller than expected for a companion-induced spiral, if the perturbing body has the semimajor axis of 20 AU. We emphasize that the morphology of the spiral structure is likely to be strongly affected by shadows of a misaligned inner disk detected by Eisner et al. (2004).

The ultraviolet (UV) bright accretion disc in active galactic nuclei (AGN) should give rise to line driving, producing a powerful wind which may play an important role in AGN feedback as well as in producing structures like the broad line region. However, coupled radiation-hydrodynamics codes are complex and expensive, so we calculate the winds instead using a non-hydrodynamical approach (the Qwind framework). The original Qwind model assumed the initial conditions in the wind, and had only simple radiation transport. Here, we present an improved version which derives the wind initial conditions and has significantly improved ray-tracing to calculate the wind absorption self consistently given the extended nature of the UV emission. We also correct the radiation flux for relativistic effects, and assess the impact of this on the wind velocity. These changes mean the model is more physical, so its predictions are more robust. We find that, even when accounting for relativistic effects, winds can regularly achieve velocities $\simeq$ (0.1-0.5) $c$, and carry mass loss rates which can be up to 30% of the accreted mass for black hole masses of $10^{7-9}$ $\mathrm{M}_\odot$, and mass accretion rates of 50% of the Eddington rate. Overall, the wind power scales as a power law with the black hole mass accretion rate, unlike the weaker scaling generally assumed in current cosmological simulations that include AGN feedback. The updated code, Qwind3, is publicly available in GitHub

A peculiar aspect of Cherenkov telescopes is that they are designed to detect atmospheric light flashes on the time scale of nanoseconds, being almost blind to stellar sources. As a consequence, the pointing calibration of these instruments cannot be done in general exploiting the standard astrometry of the focal plane. In this paper we validate a procedure to overcome this problem for the case of the innovative ASTRI telescope, developed by INAF, exploiting sky images produced as an ancillary output by its novel Cherenkov camera. In fact, this instrument implements a statistical technique called "Variance method" (VAR) owning the potentiality to image the star field (angular resolution $\sim 11'$). We demonstrate here that VAR images can be exploited to assess the alignment of the Cherenkov camera with the optical axis of the telescope down to $\sim 1''$. To this end, we evaluate the position of the stars with sub-pixel precision thanks to a deep investigation of the convolution between the point spread function and the pixel distribution of the camera, resulting in a transformation matrix that we validated with simulations. After that, we considered the rotation of the field of view during long observing runs, obtaining light arcs that we exploited to investigate the alignment of the Cherenkov camera with high precision, in a procedure that we have already tested on real data. The strategy we have adopted, inherited from optical astronomy, has never been performed on Variance images from a Cherenkov telescope until now, and it can be crucial to optimize the scientific accuracy of the incoming MiniArray of ASTRI telescopes.

We present a photometric study of the newly discovered eclipsing IP $J183221.56-162724.25$ (in short $J1832$) with an orbital period of $8.87hr$. The system features a box-like deep eclipse with a full width at 50% depth of $1970\pm2s$ and a large-amplitude coherent pulsation with $P_\mathrm{obs}\!=\!65.18min$, which represents either the synodic (beat) period or the spin period of the white dwarf (WD). The period ratio is either $P_\mathrm{spin}\!/P_\mathrm{orb} = 0.1091$ or $0.1225$, respectively. The eclipsed light originates almost entirely from the two accretion spots and columns on the WD, with characteristics indicative of pole flipping. There is no evidence for an accretion disk, and we identify J1832 as the first deeply eclipsing stream-fed intermediate polar. Our $grizy$ photometry in eclipse yielded an $i$-band AB magnitude of the Roche-lobe-filling secondary star of 18.98(3), an extinction $E_\mathrm{B-V}\!=\!0.54\!\pm\!0.17$, and a spectral type $\sim\,K6$. Dynamic models, fitting the photometry, limit the distance to between 1270 and 2500pc for masses of the secondary star, $M_2$, between $0.16$ and $1.0M_\mathrm{\odot}$, well within the Gaia EDR3 confidence limits. Employing a luminosity selection inspired by binary population studies yields a mean $M_2\!=\!0.32 M_\mathrm{\odot}$ with a 2$\sigma$ upper limit of $0.60 M_\mathrm{\odot}$ and a mean distance d = 1596pc with a 2$\sigma$ upper limit of 1980pc. The secondary star is located in its Hertzsprung-Russell diagram at a mean $T_\mathrm{eff,2}\!=\!4120K$ and $log(L_2/L_\mathrm{\odot})\!=\!-0.92$, from where the binary can evolve into either a polar or an ultracompact binary with a highly magnetic primary. The system displays a variable accretion rate, lapses repeatedly into short-lived low states of negligible accretion, and currently displays an orbital period that decreases on a timescale of $\tau\!\sim\!3*10^5yr$.

Filaments are a ubiquitous morphological feature of the molecular interstellar medium and are identified as sites of star formation. In recent years, more than 100 large-scale filaments (with a length $>10$\,pc) have been observed in the inner Milky Way. As they appear linked to Galactic dynamics, studying those structures represents an opportunity to link kiloparsec-scale phenomena to the physics of star formation, which operates on much smaller scales. In this letter, we use newly acquired Outer Galaxy High Resolution Survey (OGHReS) $^{12}$CO(2-1) data to demonstrate that a significant number of large-scale filaments are present in the outer Galaxy as well. The 37 filaments identified appear tightly associated with inter-arm regions. In addition, their masses and linear masses are, on average, one order of magnitude lower than similar-sized molecular filaments located in the inner Galaxy, showing that Milky Way dynamics is able to create very elongated features in spite of the lower gas supply in the Galactic outskirts.

The CHEX-MATE sample was built to provide an overview of the statistical properties of the underlying cluster population and to set the stage for future X-ray missions. In this work, we perform a morphological analysis of the 118 clusters included in the sample with the aim to provide a classification of their dynamical state which will be useful for future studies of the collaboration.

Context: The ionosphere of a comet is known to deflect the solar wind through mass loading, but the interaction is dependent on cometary activity. We investigate the details of this process at comet 67P using the Rosetta Ion Composition Analyzer. Aims: This study aims to compare the interaction of the solar wind and cometary ions during two different time periods in the Rosetta mission. Methods: We compare both the integrated ion moments (density, velocity, and momentum flux) and the velocity distribution functions for two days four months apart. The velocity distribution functions are projected into a coordinate system dependent on the magnetic field direction and averaged over seven hours. Results: The first case shows highly scattered H+ in both ion moments and velocity distribution function. He2+ are somewhat scattered, but less so, and look more like those of H2O+ pickup ions. The second case shows characteristic evidence of mass-loading, where the solar wind species are deflected but the velocity distribution function is not significantly changed. Conclusions: The distributions of H+ in the first case, when compared to He2+ and H2O+ pickup ions, are indicative of a narrow cometosheath on the scale of the H+ gyroradius. He2+ and H2O+, with larger gyroradii, are largely able to pass through this cometosheath. Examination of the momentum flux tensor suggests that all species in the first case have a significant non-gyrotropic momentum flux component, higher than for the second mass loaded case. Mass loading is not a sufficient explanation for the distribution functions and momentum flux tensor in the first case, and so it is evidence of bow shock formation.

To reduce and analyze astronomical images, astronomers can rely on a wide range of libraries providing low-level implementations of legacy algorithms. However, combining these routines into robust and functional pipelines requires a major effort which often ends up in instrument-specific and poorly maintainable tools, yielding products that suffer from a low-level of reproducibility and portability. In this context, we present prose, a Python framework to build modular and maintainable image processing pipelines. Built for astronomy, it is instrument-agnostic and allows the construction of pipelines using a wide range of building blocks, pre-implemented or user-defined. With this architecture, our package provides basic tools to deal with common tasks such as automatic reduction and photometric extraction. To demonstrate its potential, we use its default photometric pipeline to process 26 TESS candidates follow-up observations and compare their products to the ones obtained with AstroImageJ, the reference software for such endeavors. We show that prose produces light curves with lower white and red noise while requiring less user interactions and offering richer functionalities for reporting.

We present results of a long-term optical, UV and X-ray study of variability of the nearby changing-look (CL) Seyfert NGC 1566 which was observed with the Swift Observatory from 2007 to 2020. We summarize our previously published spectroscopic and photometric results and present new observations. We reported on the alteration in the spectral type of NGC 1566 in 2018 (REF1). Moreover, we focused on the exceptional postmaximum behavior after 2018 July, when all bands dropped with some fluctuations (REF2). We observed four significant re-brightenings in the post-maximum period. We have found differences in X-ray and UV/Optical variability. The L_uv/L_x-ray ratio was decreased during 2018-2020. New post-maximum spectra covering the period 2018 November 31 to 2019 September 23 show dramatic changes compared to 2018 August 2, with fading of the broad lines and [Fe X] 6374 until 2019 March (REF2). Effectively, two CL states were observed for this object: changing to type 1.2 and then returning to the low state as a type 1.8 Sy. We suggest that the changes are mostly due to fluctuations in the energy generation. Variability properties of NGC1566 are compared with our results for other CL AGNs.

Standard cosmological data analyses typically constrain simple phenomenological dark-energy parameters, for example the present-day value of the equation of state parameter, $w_0$, and its variation with scale factor, $w_a$. However, results from such an analysis cannot easily indicate the presence of modified gravity. Even if general relativity does not hold, experimental data could still be fit sufficiently well by a phenomenological $w_0w_a$CDM, unmodified-gravity model. Hence, it would be very useful to know if there are generic signatures of modified gravity in standard analyses. Here we present, for the first time to our knowledge, a quantitative mapping showing how modified gravity models look when (mis)interpreted within the standard unmodified-gravity analysis. Scanning through a broad space of modified-gravity (Horndeski) models, and assuming a near-future survey consisting of CMB, BAO, and SNIa observations, we report values of the best-fit set of cosmological parameters including $(w_0, w_a)$ that would be inferred if modified gravity were at work. We find that modified gravity models that can masquerade as standard gravity lead to very specific biases in standard-parameter spaces. We also comment on implications for measurements of the amplitude of mass fluctuations described by the parameter $S_8$.

We model the multiwavelength properties of binaries accreting at super-critical rates with the aim to better understand the observational properties of Ultra-luminous X-ray Sources (ULXs). We calculate an extended grid of binary systems which evolve through Roche Lobe Overflow and undergo case A mass transfer from massive donors (up to 25 Msol) onto massive Black Holes (BH) (up to 100 Msol). Angular momentum loss with the ejection of mass through an outflow is incorporated. We apply our super-Eddington accretion model to these systems, computing their evolutionary tracks on the color-magnitude diagram (CMD) for the Johnson and HST photometric systems. We found that the tracks occupy specific positions on the CMD depending on the evolutionary stage of the donor and of the binary. Moreover, their shapes are similar, regardless the BH mass. More massive BHs lead to more luminous tracks. We additionally compute their optical-through-X-ray Spectral Energy Distribution (SED) considering the effects of a Comptonizing corona which surrounds the innermost regions of the disc. We apply our model to four ULXs: NGC4559 X-7, NGC 5204 X-1, Holmberg II X-1 and NGC 5907 ULX-2. We found that accretion onto BHs with mass in the range 35-55 Msol is consistent with to the observational properties of these sources. We finally explore and discuss the possibility to extend our model also to ULXs powered by accreting Pulsars (PULXs).

The supernova remnant (SNR) B0532$-$67.5 in the Large Magellanic Cloud (LMC) was first diagnosed by its nonthermal radio emission and its SNR nature was confirmed by diffuse X-ray emission; however, no optical SNR shell is detected. The OB association LH75, or NGC 2011, is projected within the boundary of this SNR. We have analyzed the massive star population in and around SNR B0532$-$67.5: using optical photometric data to construct color-magnitude diagrams (CMDs), using stellar evolutionary tracks to estimate stellar masses, and using isochrones to assess the stellar ages. From these analyses, we find a 20-25 Myr population in LH75 and a younger population less than 10 Myr old to the southwest of LH75. The center of SNR B0532$-$67.5 is located closer to the core of LH75 than the massive stars to its southwest. We conclude that the SN progenitor was probably a member of LH75 with an initial mass $\sim$15 $M_\odot$. The SN exploded in an H I cavity excavated by the energy feedback of LH75. The low density of the ambient medium prohibits the formation of a visible nebular shell. Despite the low density in the ambient medium, physical properties of the hot gas within the SNR interior do not differ from SNRs with a visible shell by more than a factor of 2-3. The large-scale H I map shows that SNR B0532$-$67.5 is projected in a cavity that appears to be connected with the much larger cavity of the supergiant shell LMC-4.

Context. IRC+10216, the carbon-rich envelope of the asymptotic giant branch (AGB) star CW Leo, is one of the richest molecular sources in the sky. Available spectral surveys below 51 GHz are more than 25 years old and new work is needed. Aims. Characterizing the rich molecular content of this source, specially for heavy species, requires to carry out very sensitive spectral surveys at low frequencies. In particular in this work we have achieved an rms in the range 0.2-0.6 mK per MHz. Methods. long Q-band (31.0-50.3 GHz) single dish integrations were carried out with the Yebes 40m telescope using specifically built receivers. State of the art line catalogs are used for line identification. Results. A total of 652 spectral features corresponding to 713 transitions from 81 species (we count as different the isomers, isotopologues and ortho/para species) are present in the data. Only 57 unidentified lines remain with signal to noise ratios >3. Some new species and/or vibrational modes have been discovered for the first time with this survey. Conclusions. This IRC+10216 spectral survey is, by far, the most sensitive one carried out to this date in the Q-band. It therefore provides the most complete view of IRC+10216 from 31.0 to 50.3 GHz, giving unique information on its molecular content, specially for heavy species. Rotational diagrams built from the data provide valuable information on the physical conditions and chemical content of this circumstellar envelope.

Starburst-driven galactic outflows in star-forming galaxies have been observed to contain complex thermal structures and emission line features that are difficult to explain by adiabatic fluid models and plasmas in photoionization equilibrium (PIE) and collisional ionization equilibrium (CIE). We previously performed hydrodynamic simulations of starburst-driven outflows using the MAIHEM module for non-equilibrium atomic chemistry and radiative cooling functions in the hydrodynamics code FLASH, and calculated emission lines in combined CIE and PIE conditions. In the present study, we consider time-dependent non-equilibrium ionization (NEI) states produced by our MAIHEM simulations. Through extensive CLOUDY calculations made with the NEI states from our hydrodynamic simulations, we predict the UV and optical emission line profiles for starburst-driven outflows in time-evolving non-equilibrium photoionization conditions. Our hydrodynamic results demonstrate applications of non-equilibrium radiative cooling for H II regions in starburst galaxies hosting cool outflows.

This study aims to determine the main physical parameters (N(H2) hydrogen column density and Td dust temperature) of the Interstellar medium, and their distribution in the extended star-forming region, which includes IRAS 05156+3643, 05162+3639, 05168+3634, 05177+3636, and 05184+3635 sources. We also provide a comparative analysis of the properties of the Interstellar medium and young stellar objects. Analysis of the results revealed that Interstellar medium forms relativity dense condensations around IRAS sources, which are interconnected by a filament structure. In general, in sub-regions Td varies from 11 to 24 K, and N(H2) - from 1.0 to 4.0 x 10^23 cm^(-2). The masses of the ISM vary from 1.7 x 10^4 to 2.1 x 10^5 Msol. All BGPSv2 objects identified in this star-forming region are located at the N(H2) maximum. The direction of the outflows, which were found in two sub-regions, IRAS 05168+3634 and 05184+3635, correlates well with the isodenses direction. The sub-regions with the highest N(H2) and Interstellar medium mass have the largest percentage of young stellar objects with Class I evolutionary stage. The wide spread of the evolutionary ages of stars in all sub-regions (from 10^5 to 10^7 years) suggests that the process of star formation in the considered region is sequential. In those sub-regions where the mass of the initial, parent molecular cloud is larger, this process is likely to proceed more actively. On the Gaia EDR3 database, it can be assumed that all sub-regions are embedded in the single molecular cloud and belong to the same star-forming region, which is located at a distance of about 1.9 kpc.

We present a full set of numerical tools to extract Galaxy Cluster pressure profiles from the joint analysis of Planck and South Pole Telescope (SPT) observations. Pressure profiles are powerful tracers of the thermodynamic properties and the internal structure of the clusters. Tracing the pressure over the cosmic times allows to constraints the evolution of the cluster structure and the contribution of astrophysical phenomena. SPT and Planck are complementary to constrain the cluster structure at various spatial scales. The SPT cluster catalogue counts 677 cluster candidates up to redshift 1.7, it is a nearly mass limited sample, an ideal benchmark to test cluster evolution. We developed a pipeline to first separate the cluster signal from the background and foreground components and then jointly fit a parametric profile model on a combination of Planck and SPT data. We validate our algorithm on a sub-sample of six clusters, common to the SPT and the CHEX-MATE catalogues, comparing the results with the profiles obtained from X-ray observations with XMM-Newton.

We present results constraining the multiplicity of the very low mass stars and sub-stellar objects in the Orion Nebula Cluster (ONC). Our sample covers primary masses 0.012-0.1M$_{\odot}$ using archival Hubble Space Telescope data obtained with the Advanced Camera for Surveys using multiple filters. Studying the binary populations of clusters provides valuable constraints of how the birth environment affects binary formation and evolution. Prior surveys have shown that the binary populations of high-mass, high-density star clusters like the ONC may substantially differ from those in low-mass associations. Very low mass stellar and sub-stellar binaries at wide separations, $>$20AU, are statistically rare in the Galactic field and have been identified in stellar associations like Taurus-Auriga and Ophiuchus. They also may be susceptible to dynamical interactions, and their formation may be suppressed by feedback from on-going star formation. We implement a double point-spread function (PSF) fitting algorithm using empirical, position dependent PSF models to search for binary companions at projected separations $>$ 10 AU (25 mas). With this technique, we identify 7 very low mass binaries, 5 of which are new detections, resulting in a binary frequency of 12$^{+6}_{-3.2}$% over mass ratios of 0.5 - 1.0 and projected separations of 20 - 200 AU. We find an excess of very low mass binaries in the ONC compared to the Galactic field, with a probability of 10$^{-6}$ that the populations are statistically consistent. The sub-stellar population of the ONC may require further dynamical processing of the lowest binding energy binaries to resemble the field population.

Solar variability occurs over a broad range of spatial and temporal scales, from the Sun's brightening over its lifetime to the fluctuations commonly associated with magnetic activity over minutes to years. The latter activity includes most prominently the 11-year sunspot solar cycle and its modulations. Space weather events, in the form of solar flares, solar energetic particles, coronal mass ejections, and geomagnetic storms, have long been known to approximately follow the solar cycle occurring more frequently at solar maximum than solar minimum. These events can significantly impact our advanced technologies and critical infrastructures, making the prediction for the strength of future solar cycles particularly important. Several methods have been proposed to predict the strength of the next solar cycle, cycle 25, with results that are generally not always consistent. Most of these methods are based on the international sunspot number time series, or other indicators of solar activity. We present here a new approach that uses more than 100 years of measured fractional areas of the visible solar disk covered by sunspots and plages and an empirical relationship for each of these two indices of solar activity in even-odd cycles. We anticipate that cycle 25 will peak in 2024 and will last for about 12 years, slightly longer than cycle 24. We also found that, in terms of sunspot and plage areas coverage, the amplitude of cycle 25 will be substantially similar or slightly higher than cycle 24.

The component black holes (BHs) observed in gravitational-wave (GW) binary black hole (BBH) events tend to be more massive and slower spinning than those observed in black hole X-ray binaries (BH-XRBs). Without the use of astrophysical models, we investigate whether these apparent tensions can be explained by GW observational selection effects alone. We find that this is indeed the case for the discrepancy between BH masses in BBHs and the observed high-mass X-ray binaries (HMXBs), when we account for statistical uncertainty from the small sample size of just three HMXBs. On the other hand, the BHs in observed low-mass X-ray binaries (LMXBs) are significantly lighter than the underlying BBH population, but this may just be due to a correlation between component masses in a binary system. Given their light stellar companions, we expect light BHs in LMXBs. The observed spins in HMXBs and LMXBs, however, are in tension with the inferred BBH spin distribution at the $>99.9\%$ level. We discuss possible scenarios behind the significantly larger spins in observed BH-XRBs. One possibility is that a small subpopulation (conservatively $< 30\%$) of BBHs have rapidly spinning primary components, indicating that they may have followed a similar evolutionary pathway to the observed HMXBs. In LMXBs, it has been suggested that BHs can spin up by accretion. If LMXB natal spins follow the BBH spin distribution, we find LMXBs must gain an average dimensionless spin of $0.47^{+0.10}_{-0.11}$, but if their natal spins follow the observed HMXB spins, the average spinup must be $< 0.03$

The extended stellar halos of galaxies contain important clues for investigating their assembly history and evolution. We investigate the resolved stellar content and the extended halo of NGC 5128 as a function of galactocentric distance. We used HST images to resolve individual red giant branch (RGB) stars in 28 independent pointings. Star counts from deep VI color-magnitude diagrams reaching at least 1.5 mag below the tip of the RGB are used to derive the surface density distribution of the halo. The contamination by Milky Way stars is assessed with a new control field, with models, and by combining optical and near-IR photometry. We present a new calibration of the WFC3 F606W+F814W photometry to the ground-based VI photometric system. The photometry shows that the stellar halo of NGC 5128 is dominated by old RGB stars that are present in all fields. The V-band surface brightness changes from 23 to 32 mag/arcsec$^2$ between 8.3 kpc from the galaxy center to our outermost halo fields located 140 kpc away from the center along the major axis and 92 kpc along the minor axis. Within ~30 kpc, we also find evidence for a 2-3 Gyr old population traced by bright asymptotic giant branch stars. This population contributes only up to 10% in total stellar mass if it is 2 Gyr old, but a larger fraction of 30-40% is required if its age is 3 Gyr. The stellar surface density profile is well fit by a r$^{1/4}$ curve or a power-law $\sim r^{-3.1}$ over the full radial range, with no obvious break in the slope, but with large field-to-field scatter. The ellipticity measured from integrated-light photometry in the inner parts, $e=(b/a)=0.77$, flattens to $e=0.54 \pm 0.02$ beyond 30 kpc. Considering the flattening of the outer halo, the projection of the elliptical isophote on the semimajor axis for our most distant field reaches nearly 30 effective radii. [abridged]

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma ray data, and so on. It has been claimed that only the DM annihilation to a muon-pair, namely muonphilic DM, is compatible with the null detection of all the corresponding astrophysical measurements [Di Mauro and Winkle (2021)]. On the other hand, a muonphilic DM model may also lead to a signal in the recent muon g-2 measurement or the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including different DM and mediator spins. In agreement with GCE (not only $2\mu$ but also $4\mu$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the mediator resonance region can explain GCE and relic density simultaneously. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon g-2 excess on top of GCE, relic density, and other DM and mediator constraints.

The orbital International Gamma-Ray Astrophysics Laboratory (INTEGRAL), launched in 2002, continues its successful work in observing the sky at energies above 20 keV. The growing INTEGRAL data archive allows one to conduct a hard X-ray all-sky survey including a number of deep extragalactic fields and the deepest ever hard X-ray survey of the Galaxy. Taking advantage of the data gathered over 17 years with the IBIS coded-mask telescope on board INTEGRAL, we conducted a survey of hard X-ray sources in the 17-60 keV band, providing flux information in different energy bands up to 290 keV. The catalog of sources includes 929 objects, 890 of which exceed a detection threshold of 4.5 sigma and the rest are detected at 4.0-4.5 sigma and belong to known cataloged INTEGRAL sources and sources from the on-going all-sky survey by the BAT telescope of the Neil Gehrels Swift Observatory. Among the identified sources with known or suspected nature, 376 are associated with Galaxy and Magellanic clouds, including 145 low-mass X-ray binaries, 115 high-mass X-ray binaries, 79 cataclysmic variables, and 37 of other types; and 440 are extragalactic, including 429 active galactic nuclei (AGNs), 2 ultra-luminous sources (ULXs), supernova remnant AT2018cow and 8 galaxy clusters. 113 sources from the catalog remain unclassified. 46 objects are detected in the hard X-ray band for the first time. The cumulative LogN-LogS distribution of non-blazar AGNs, based on 356 sources detected at S/N>4.5 sigma, is measured down to flux 2E-12 erg/s/cm^2 and can be described by a power law with a slope of 1.44 +/- 0.09 and a normalization of 8E-3 sources per deg^2 at fluxes >1E-11 erg/s/cm^2. The LogN-LogS distribution of unclassified sources indicates that the majority of them are of extragalactic origin.

We explore the connection between the $\gamma$-ray and radio emission in the jet of the blazar 0716$+$714 by using 15, 37, and 230 GHz radio and 0.1$-$200 GeV $\gamma$-ray light curves spanning 10.5 years (2008$-$2019). We find significant positive and negative correlations between radio and $\gamma$-ray fluxes in different time ranges. The time delays between radio and $\gamma$-ray emission suggest that the observed $\gamma$-ray flares originated from multiple regions upstream of the radio core, within a few parsecs from the central engine. Using time-resolved 43 GHz VLBA maps we identified 14 jet components moving downstream along the jet. Their apparent speeds range from 6 to 26 $c$, showing notable variations in their position angles upstream the stationary component ($\sim$0.53 mas from the core). The brightness temperature declines as function of distance from the core according to a power-law which becomes shallower at the location of the stationary component. We also find that the periods at which significant correlations between radio and $\gamma$-ray emission occur overlap with the times when the jet was oriented to the north. Our results indicate that the passage of a propagating disturbance (or shock) through the radio core and the orientation of the jet might be responsible for the observed correlation between the radio and $\gamma$-ray variability. We present a scenario that connects the positive correlation and the unusual anti-correlation by combining the production of a flare and a dip at $\gamma$-rays by a strong moving shock at different distances from the jet apex.

X-ray polarimetry will soon open a new window on the high energy universe with the launch of NASA's Imaging X-ray Polarimetry Explorer (IXPE). Polarimeters are currently limited by their track reconstruction algorithms, which typically use linear estimators and do not consider individual event quality. We present a modern deep learning method for maximizing the sensitivity of X-ray telescopic observations with imaging polarimeters, with a focus on the gas pixel detectors (GPDs) to be flown on IXPE. We use a weighted maximum likelihood combination of predictions from a deep ensemble of ResNets, trained on Monte Carlo event simulations. We derive and apply the optimal event weighting for maximizing the polarization signal-to-noise ratio (SNR) in track reconstruction algorithms. For typical power-law source spectra, our method improves on the current state of the art, providing a ~40% decrease in required exposure times for a given SNR.

We have investigated the structure of the Hercules supercluster (SCL160) based on data originally extracted from the Sloan Digital Sky Survey SDSS-DR7. We have traced the mass distribution in the field through the numerical density-weighted by the $r^\prime$-luminosity of the galaxies and classified them based on their spatial position and redshift. This has allowed us not only to address the kinematics of the supercluster as a whole, but also the internal kinematic of each cluster, which was no further explored before. We have confirmed that the Hercules supercluster is composed of the galaxy clusters A2147, A2151, and A2152. A2151 consists of five subclusters, A2147 on two and A2152 on at least two. They form the heart of the Hercules supercluster. We also have found two other gravitationally bond clusters, increasing, therefore, the known members of the supercluster. We have estimated a total mass of $2.1\pm0.2 \times 10^{15}$ M$_\odot$ for the Hercules supercluster. To determine the dynamical masses in this work, we have resorted to the $M_{200}-\sigma$ scaling relation and the caustic technique. Comparing both methods with simulated data of bimodal merging clusters, we found the caustic, as well as the $\sigma$-based masses, are biased through the merger age, showing a boost just after the pericentric passage. This is not in line with the principle of the caustic method that affirms it is not depending on the cluster dynamical state.

We determine the proper motion of the Solar system from the Pantheon sample of supernovae (SNe) of type Ia. The posterior distribution of the Solar system proper velocity, its direction and relevant cosmological parameters are obtained based on the observed distance moduli, heliocentric redshifts, and positions of SNe by means of a Markov Chain Monte Carlo method. We account for the unknown peculiar motion of SNe by including their expected covariance from linear theory. We find that the Solar system moves with $v_o = 249 \pm 51$ km/s towards $RA = 166 \pm 16$ deg, $Dec = 10 \pm 19$ deg (J2000) (all at 68\% C.L.). The direction of motion agrees with the direction of the dipole observed in the cosmic microwave background (CMB) ($RA = 166$ deg, $Dec = -7$ deg). The inferred velocity is $2.4 \sigma$ smaller than the value inferred from a purely kinematic interpretation of the CMB dipole ($370$ km/s). Assuming a flat $\Lambda$ cold dark matter model, we find no degeneracy of Solar proper motion with other cosmological parameters. The dimensionless matter density is $\Omega_M = 0.305 \pm 0.022$, in excellent agreement with CMB measurements. We also find no degeneracy of the Solar proper motion with the SN calibration nuisance parameter. We conclude that a larger sample of SNe will allow an independent and robust test of the kinematic nature of the CMB dipole.

The recent observational evidence for cosmic filament spin on megaparsec scales (Wang et al, Nature Astronomy 5, 839-845 (2021)) demands an explanation in the physics of dark matter. Conventional collisionless cold particle dark matter is conjectured to generate cosmic filament spin through tidal torquing, but this explanation requires extrapolating from the quasi-linear regime to the non-linear regime. Meanwhile no alternative explanation exists in the context of ultra-light (e.g., axion) dark matter, and indeed these models would naively predict zero spin for cosmic filaments. In this Letter we study cosmic filament spin in theories of ultra-light dark matter, such as ultra-light axions, and bosonic and fermionic condensates, such as superfluids and superconductors. These models are distinguished from conventional particle dark matter models by the possibility of dark matter vortices. We take a model agnostic approach, and demonstrate that a collection of dark vortices can explain the data reported in Wang et al. Modeling a collection of vortices with a simple two-parameter analytic model, corresponding to an averaging of the velocity field, we find an excellent fit to the data. We perform a Markov Chain Monte Carlo analysis and find constraints on the number of vortices, the dark matter mass, and the radius of the inner core region where the vortices are distributed, in order for ultra-light dark matter to explain spinning cosmic filaments.

The effect of gravitational wave of extremely low frequency on time delays between different locations on the Einstein ring in a lens system with an aligned source-deflector-observer configuration is investigated. The observer will observe an Einstein ring from the lens system aligned in a highly symmetric configuration. Time delays between different locations on the Einstein ring cannot emerge without gravitational wave of cosmological wavelength propagating through the lens system. Otherwise, different locations on the Einstein ring will show time delays. Our previous studies showed that time delays from different locations on the Einstein ring in the presence of gravitational wave of cosmological wavelength present a special relationship. But this result is limited to a specific aligned source-deflector-observer configuration where the source and the observer are equidistant from the deflector and the gravitational wave has special direction of propagation and polarization. In order to investigate the property in the general case, we extend to the general condition that the source and the observer aren't equidistant from the deflector and the direction and polarization of the gravitational wave is arbitrary in the present work. Results in this work show that time delays between different locations on the Einstein ring with the general condition present the same relationship as that in our previous studies.

We introduce NRPyElliptic, an elliptic solver for numerical relativity (NR) built within the NRPy+ framework. As its first application, NRPyElliptic sets up conformally flat, binary black hole (BBH) puncture initial data (ID) on a single numerical domain, similar to the widely used TwoPunctures code. Unlike TwoPunctures, NRPyElliptic employs a hyperbolic relaxation scheme, whereby arbitrary elliptic PDEs are trivially transformed into a hyperbolic system of PDEs. As consumers of NR ID generally already possess expertise in solving hyperbolic PDEs, they will generally find NRPyElliptic easier to tweak and extend than other NR elliptic solvers. When evolved forward in (pseudo)time, the hyperbolic system exponentially reaches a steady state that solves the elliptic PDEs. Notably NRPyElliptic accelerates the relaxation waves, which makes it many orders of magnitude faster than the usual constant-wavespeed approach. While it is still ${\sim}12$x slower than TwoPunctures at setting up full-3D BBH ID, NRPyElliptic requires only ${\approx}0.3\%$ of the runtime for a full BBH simulation in the Einstein Toolkit. Future work will focus on improving performance and generating other types of ID, such as binary neutron star.

Phonon-mediated particle detectors employing Kinetic Inductance Detectors (KIDs) on Silicon substrates have demonstrated both O(10) eV energy resolution and mm position resolution, making them strong candidates for instrumenting next generation rare-event experiments such as in looking for dark matter or in neutrino measurements. Previous work has demonstrated the performance of an 80-KID array on a Si wafer, however current energy resolution measurements show a 25x difference between otherwise identical KIDs on the same wafer - between 5 to 125 eV on energy absorbed by the KID. Here, we use a first principles approach and attempt to identify the drivers behind the resolution variation. In particular, we analyze a subset of 8 KIDs using the unique approach of pulsing neighboring KIDs to generate signals in the target. We tentatively identify differences in quality factor terms as the likely culprit for the observed variation.

The NEWS-G direct dark matter search experiment uses spherical proportional counters (SPC) with light noble gases to explore low WIMP masses. The first results obtained with an SPC prototype operated with Ne gas at the Laboratoire Souterrain de Modane (LSM) have already set competitive results for low-mass WIMPs. The forthcoming next phase of the experiment consists of a large 140 cm diameter SPC installed at SNOLAB with a new sensor design, with improved detector performance and data quality. Before its installation at SNOLAB, the detector was commissioned with pure methane gas at the LSM, with a temporary water shield, offering a hydrogen-rich target and reduced backgrounds. After giving an overview of the improvements of the detector, preliminary results of this campaign will be discussed, including UV laser and Ar-37 calibration data.

The abundance and fractionation of the stable strontium (Sr) isotope system are being increasingly utilized to move forward our understanding in geological and cosmological processes. Two analytical techniques are commonly used to measure stable Sr isotopes: 1) double-spike thermal ionization mass spectrometry (DS-TIMS) and 2) Zr-doped sample-standard bracketing multi-collector inductively coupled plasma mass spectrometry (Zr-doped SSB via MC-ICP-MS). Relative to DS-TIMS, Zr-doped SSB via MC-ICP-MS allows simultaneous determinations of both 87Sr/86Sr and 88Sr/86Sr ratios, increasing measurement efficiency and sample throughput. However, this technique is currently associated with greater uncertainties in measurement precision and accuracy. In this study, we evaluated potential factors that can affect the quality of Sr isotope measurements during Zr-doped SSB. Our tests show that incomplete Sr recovery during chromatographic separation, mismatches of Sr and Zr concentrations and acid molarity between sample and bracketing standard, and cation contamination could all affect the precision and accuracy of Sr isotope measurements. We present evidence that, with updated preparation procedures and diligent concentration checks, a long-term reproducibility (2{\sigma}SD: 87Sr/86Sr = +/-0.000015 and {\delta}88/86Sr = +/-0.03 permil) comparable to that of DS-TIMS is achievable when using the Zr-doped SSB method via MC-ICP-MC.