Papers for Monday, Jan 17 2022

Hot atmospheres pervading galaxy clusters, groups, and early-type galaxies are rich in metals, produced during epochs and diffused via processes that are still to be determined. While this enrichment has been routinely investigated in clusters, metals in lower mass systems are more challenging to probe with standard X-ray exposures and spectroscopy. In this paper, we focus on very deep XMM-Newton ($\sim$350 ks) observations of NGC 1404, a massive elliptical galaxy experiencing ram-pressure stripping of its hot atmosphere while infalling toward the centre of the Fornax cluster, with the aim to derive abundances through its hot gas extent. Importantly, we report the existence of a new fitting bias - the "double Fe bias" - leading to an underestimate of the Fe abundance when two thermal components cannot realistically model the complex temperature structure present in the outer atmosphere of the galaxy. Contrasting with the ''metal conundrum'' seen in clusters, the Fe and Mg masses of NGC 1404 are measured 1-2 orders of magnitude below what stars and supernovae could have reasonably produced and released. In addition, we note the remarkable Solar abundance ratios of the galaxy's halo, different from its stellar counterpart but similar to the chemical composition of the ICM of rich clusters. Completing the clusters regime, all these findings provide additional support toward a scenario of early enrichment, at play over two orders of magnitude in mass. A few peculiar and intriguing features, such as a possible double metal peak as well as an apparent ring of enhanced Si near the galaxy core, are also discussed.

More than 50 per cent of present-day massive disc galaxies show a rotating stellar bar. Their formation and dynamics have been widely studied both numerically and observationally. Although numerical simulations in the $\Lambda$CDM cosmological framework predict the formation of such stellar components, there seems to be a tension between theoretical and observational results. Simulated bars are typically larger in size and have slower pattern speed than observed ones. We study the formation and evolution of barred galaxies, using two $\Lambda$CDM zoom-in hydrodynamical simulations of the CLUES project that follow the evolution of a cosmological Local Group-like volume. We found that our simulated bars, at $z = 0$, are both shorter and faster rotators than previous ones found in other studies on cosmological simulations alleviating the tension mentioned above. These bars match the short tail-end of the observed bar length distribution. In agreement with previous numerical works, we find that bars form in those systems where the disc self-gravity is dominant over the dark matter halo, making them unstable against bar formation. Our bars developed in the last 3-4 Gyr until they achieve their current length and strength; as bars grow, their lengths increase while their rotation speeds decrease. Despite this slowdown, at redshift $z = 0$ their rotation speeds and size match well the observational data.

We present new thermophysical model (TPM) fits of 1,847 asteroids, deriving thermal inertia, diameter, and Bond and visible geometric albedo. We use thermal flux measurements obtained by the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010; Mainzer et al. 2011) during its fully cryogenic phase, when both the 12$\mu$m (W3) and 22$\mu$m (W4) bands were available. We take shape models and spin information from the Database of Asteroid Models from Inversion Techniques (DAMIT; \v{D}urech et al. 2010) and derive new shape models through lightcurve inversion and combining WISE photometry with existing DAMIT lightcurves. When we limit our sample to the asteroids with the most reliable shape models and thermal flux measurements, we find broadly consistent thermal inertia relations with recent studies. We apply fits to the diameters $D$ (km) and thermal inertia $\Gamma$ (J m$^{-2}$ s$^{-0.5}$ K$^{-1}$) normalized to 1 au with a linear relation of the form $\log[\Gamma]=\alpha+\beta\log[D]$, where we find $\alpha = 2.667 \pm 0.059$ and $\beta = -0.467 \pm 0.044$ for our sample alone and $\alpha = 2.509 \pm 0.017$ and $\beta = -0.352 \pm 0.012$ when combined with other literature estimates. We find little evidence of any correlation between rotation period and thermal inertia, owing to the small number of slow rotators to consider in our sample. While the large uncertainties on the majority of our derived thermal inertia only allow us to identify broad trends between thermal inertia and other physical parameters, we can expect a significant increase in high-quality thermal flux measurements and asteroid shape models with upcoming infrared and wide-field surveys, enabling even more thermophysical modeling of higher precision in the future.

We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities.

Alfv\'{e}nic component of MHD turbulence damps Alfv\'{e}nic waves. The consequences of this effect are important for many processes, from cosmic ray (CR) propagation to launching outflows and winds in galaxies and other magnetized systems. We discuss the differences in the damping of the streaming instability by turbulence and the damping of a plane parallel wave. The former takes place in the system of reference aligned with the local direction of magnetic field along which CRs stream. The latter is in the reference frame of the mean magnetic field and traditionally considered in plasma studies. We also compare the turbulent damping of streaming instability with ion-neutral collisional damping, which becomes the dominant damping effect at a sufficiently low ionization fraction. Numerical testing and astrophysical implications are also discussed.

We aim at characterizing the hot phase of the Circum-Galactic Medium in a large sample of galaxies. We stack X-ray events from the SRG/eROSITA eFEDS survey around central galaxies in the GAMA 9hr field to construct radially projected soft X-ray luminosity profiles as a function of their stellar mass and specific star formation rate. We consider samples of quiescent (star-forming) galaxies in the stellar mass range $2\times 10^{10}$ -- $10^{12}$ M$_\odot$ ($3\times 10^9$ -- $6\times10^{11}$ M$_\odot$). For quiescent galaxies, the X-ray profiles are clearly extended throughout the available mass range; however, the measured profile is likely biased high due to projection effects, as these galaxies tend to live in dense and hot environments. For the most massive star forming samples ($\geq10^{11}$ M$_\odot$), there is a hint of detection of extended emission. For star-forming galaxies with $< 10^{11}$ M$_\odot$ the X-ray stacked profiles are compatible with unresolved sources and consistent with the expected emission from faint Active Galactic Nuclei and X-ray binaries. We measure for the first time the mean relation between average X-ray luminosity and stellar mass separately for quiescent and star-forming galaxies. High-mass ($\geq 10^{11}$ M$_\odot$) star-forming or quiescent galaxies follow the expected scaling of virialized hot haloes, while lower mass star-forming galaxies show a less prominent luminosity and a weaker dependence on stellar mass, consistent with empirical models of the weak AGN population. When comparing our results with state-of-the art numerical simulations, we find an overall consistency on large ($>80$ kpc) scales at masses $\geq 10^{11}$ M$_\odot$, but disagreement on the small scales, where brighter than observed compact cores are predicted. Simulations also do not predict the clear differentiation that we observe between quiescent and star-forming galaxies.

Large excursion of the inflaton field can trigger interesting dynamics. One important example is a first-order phase transition in a spectator sector which couples to the inflaton. Gravitational waves (GWs) from such a first-order phase transition during inflation, an example of an instantaneous source, have an oscillatory feature. In this work, we show that this feature is generic for a source in an era of accelerated expansion. We also demonstrate that the shape of the GW signal contains information about the evolution of the early universe following the phase transition. In particular, the slope of the infrared part of the GW spectrum is sensitive to the evolution of the Hubble parameter when the GW modes reenter the horizon after inflation. The slope of the profile of the intermediate oscillatory part and the ultraviolet part of the GW spectrum depend on the evolution of the Hubble parameter when the modes exit horizon during the inflation and when they reenter the horizon during the reheating. The ultraviolet spectrum also depends on the details of the dynamics of the phase transition. We consider the GW signal in several models of evolution during and after inflation, and compare them with the minimal scenario of quasi-de Sitter inflation followed by radiation domination after a fast reheating, and demonstrate that the shape of the GW can be used to distinguish them. In this way, the GW signal considered in this paper offers a powerful probe to the dynamics of the early universe which is otherwise difficult to explore directly through CMB, large scale structure, big bang nucleosynthesis (BBN), and other well-studied cosmological observables.

The redshift of $z\sim1.5$ is the transition epoch of protoclusters (PCs) from the star-forming phase into the quenching phase, and hence an appropriate era to investigate the build up of the quenched population. We define a `core' as the most massive halo in a given PC, where environmental effects are likely to work most effectively, and search for cores at $1<z<1.5$. We use a photometric redshift catalogue of a wide (effective area of $\sim22.2\,\mathrm{deg}^{2}$) and deep ($i\sim26.8\,\mathrm{mag}$) optical survey with Subaru Hyper-Suprime Cam. Regarding galaxies with $\log(M_{*}/M_{\odot})>11.3$ as the central galaxies of PC cores, we estimate their average halo mass by clustering analysis and find it to be $\log(M_\mathrm{h}/M_{\odot})\sim13.7$. An expected mass growth by the Illustris TNG simulation and the observed overdensities around them suggest that the PC cores we find are progenitors of present-day clusters. Classifying our galaxy sample into red and blue galaxies, we calculate the stellar mass function (SMF) and the red galaxy fraction. The SMFs in the PC cores are more-top heavy than field, implying early high-mass galaxy formation and disruption of low-mass galaxies. We also find that the red fraction increases with stellar mass, consistent with stellar-mass dependent environmental quenching recently found at $z>1$. Interestingly, although the cores with red and blue centrals have similar halo masses, only those with red centrals show a significant red fraction excess compared to the field, suggesting a conformity effect. Some observational features of PC cores may imply that the conformity is caused by assembly bias.

We describe the use of Multi Order Coverage (MOC) maps as a practical way to manage complex regions of the sky for the planning of multi-messenger observations. MOC maps are a data structure that provides a multi-resolution representation of irregularly shaped and fragmentary regions over the sky based on the HEALPix (Hierarchical Equal Area isoLatitude Pixelization) tessellation. We present a new application of MOC, in combination with the \texttt{astroplan} observation planning package, to enable the efficient computation of sky regions and the visibility of these regions from a specific location on the Earth at a particular time. Using the example of the low-latency gravitational-wave alerts, and a simulated observational campaign with three observatories, we show that the use of MOC maps allows a high level of interoperability to support observing schedule plans. Gravitational-wave detections have an associated credible region localization on the sky. We demonstrate that these localizations can be encoded as MOC maps, and how they can be used in visualisation tools, and processed (filtered, combined) and also their utility for access to Virtual Observatory services which can be queried 'by MOC' for data within the region of interest. The ease of generating the MOC maps and the fast access to data means that the whole system can be very efficient, so that any updates on the gravitational-wave sky localization can be quickly taken into account and the corresponding adjustments to observing schedule plans can be rapidly implemented. We provide example python code as a practical example of these methods. In addition, a video demonstration of the entire workflow is available.

We investigate a cosmological scenario in which the dark matter particles can be created during the evolution of the Universe. By regarding the Universe as an open thermodynamic system and using non-equilibrium thermodynamics, we examine the mechanism of gravitational particle production. In this setup, we study the large-scale structure (LSS) formation of the Universe in the Newtonian regime of perturbations and derive the equations governing the evolution of the dark matter overdensities. Then, we implement the cosmological data from Planck 2018 CMB measurements, SNe Ia and BAO observations, as well as the SH0ES local measurement for $H_0$ to provide some cosmological constraints for the parameters of our model. We see that the best case of our scenario ($\chi_{{\rm tot}}^{2}=3834.40$) fits the observational data better than the baseline $\Lambda$CDM model ($\chi_{{\rm tot}}^{2} = 3838.00$) at the background level. We also see that this case results in the Hubble constant as $H_0 = 68.79\pm 0.59\,{\rm km\,s^{-1}\,Mpc^{-1}}$ which is greater than $H_0 = 68.20^{+0.42}_{-0.38}\,{\rm km\,s^{-1}\,Mpc^{-1}}$ given by the $\Lambda$CDM model, and hence we can alleviate the $H_0$ tension to some extent in our framework. Furthermore, the best case of our scenario gives a lower value for the best-fit of the $S_8$ parameter than the $\Lambda$CDM result, and therefore it also reduces the LSS tension slightly. We moreover estimate the growth factor of linear perturbations and show that the best case of our model ($\chi_{f\sigma_{8}}^{2}=40.84$) fits the LSS data significantly better than the $\Lambda$CDM model ($\chi_{f\sigma_{8}}^{2}=44.29$). Consequently, our model also makes a better performance at the level of the linear perturbations compared to the standard cosmological model.

We present high signal-to-noise measurements of three-point shear correlations and the third moment of the mass aperture statistic using the first 3 years of data from the Dark Energy Survey. We additionally obtain the first measurements of the configuration and scale dependence of the four three-point shear correlations which carry cosmological information. With the third-order mass aperture statistic, we present tomographic measurements over angular scales of 4 to 60 arcminutes with a combined statistical significance of 15.0$\sigma$. Using the tomographic information and measuring also the second-order mass aperture, we additionally obtain a skewness parameter and its redshift evolution. We find that the amplitudes and scale-dependence of these shear 3pt functions are in qualitative agreement with measurements in a mock galaxy catalog based on N-body simulations, indicating promise for including them in future cosmological analyses. We validate our measurements by showing that B-modes, parity-violating contributions and PSF modeling uncertainties are negligible, and determine that the measured signals are likely to be of astrophysical and gravitational origin.

Dating stars is a major challenge with a deep impact on many astrophysical fields. One of the most promising techniques for this is using chemical abundances. Recent space- and ground-based facilities have improved the quantity of stars with accurate observations. This has opened the door for using Bayesian inference tools to maximise the information we can extract from them. Our aim is to present accurate and reliable stellar age estimates of FGK stars using chemical abundances and stellar parameters. We used one of the most flexible Bayesian inference techniques (hierarchical Bayesian models) to exceed current possibilities in the use of chemical abundances for stellar dating. Our model is a data-driven model. We used a training set that has been presented in the literature with ages estimated with isochrones and accurate stellar abundances and general characteristics. The core of the model is a prescription of certain abundance ratios as linear combinations of stellar properties including age. We gathered four different testing sets to assess the accuracy, precision, and limits of our model. We also trained a model using chemical abundances alone. We found that our age estimates and those coming from asteroseismology, other accurate sources, and also with ten Gaia benchmark stars agree well. The mean absolute difference of our estimates compared with those used as reference is 0.9 Ga, with a mean difference of 0.01 Ga. When using open clusters, we reached a very good agreement for Hyades, NGC 2632, Ruprecht 147, and IC4651. We also found outliers that are a reflection of chemical peculiarities and/or stars at the limit of the validity ranges of the training set. The model that only uses chemical abundances shows slightly worse mean absolute difference (1.18 Ga) and mean difference (-0.12 Ga).

SXP 15.3 and SXP 305 are two Be X-ray binaries in the Small Magellanic Cloud that are spatially separated by ~7 arcsec. The small separation between these sources has, in the past, resulted in confusion about the origin of the emission from the combined region. We present long-term optical and X-ray monitoring results of both sources, where we study the historic and recent behaviour. In particular, from data collected as part of the S-CUBED project we see repeating X-ray outbursts from the combined region of the two sources in the recent lightcurve from the Neil Gehrels Swift Observatory, and we investigate the origin of this emission. Using the H-alpha emission line from the Southern African Large Telescope (SALT) and photometric flux from the Optical Gravitational Lensing Experiment (OGLE) to study the changes in the size and structure of the Be disc, we demonstrate that the X-ray emission likely originates from SXP 15.3. Timing analysis reveals unusual behaviour, where the optical outburst profile shows modulation at twice the frequency of the X-ray outbursts. We consider either of these periodicities being the true orbital period in SXP 15.3 and propose models based on the geometric orientations of the Be disc and neutron star to explain the physical origin of the outbursts.

We show that the Platinum gamma-ray burst (GRB) data compilation, probing the redshift range $0.553 \leq z \leq 5.0$, obeys a cosmological-model-independent three-parameter fundamental plane (Dainotti) correlation and so is standardizable. While they probe the largely unexplored $z \sim 2.3-5$ part of cosmological redshift space, the GRB cosmological parameter constraints are consistent with, but less precise than, those from a combination of baryon acoustic oscillation (BAO) and Hubble parameter [$H(z)$] data. In order to increase the precision of GRB-only cosmological constraints, we exclude common GRBs from the larger Amati-correlated A118 data set composed of 118 GRBs and jointly analyze the remaining 101 Amati-correlated GRBs with the 50 Platinum GRBs. This joint 151 GRB data set probes the largely unexplored $z \sim 2.3-8.2$ region; the resulting GRB-only cosmological constraints are more restrictive, and consistent with, but less precise than, those from $H(z)$ + BAO data.

We investigate the combined evolution of the dipolar surface magnetic field (B$_{s}$) and the spin-period (P$_s$) of known magnetars and high magnetic field (B$_s$ $ \gtrsim 10^{13}$~G) radio pulsars. We study the long term behaviour of these objects assuming a simple ohmic dissipation of the magnetic field. Identifying the regions (in the P$_s$-B$_s$ plane) in which these neutron stars would likely move into, before crossing the death-line to enter the pulsar graveyard, we comment upon the possible connection between the magnetars and other classes of neutron stars.

The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small ($\sim$10 per cent) reference sample of high-quality redshifts. Within pencil-beam-like sub-volumes surrounding each galaxy, we use the precise dN/d$z$ distribution of the reference sample to recover new redshifts and assign them one-to-one to galaxies such that the original rank order of redshifts is preserved. Preserving the rank order is motivated by the fact that random variables drawn from Gaussian probability density functions with different means but equal standard deviations satisfy stochastic ordering. The process is repeated for sub-volumes surrounding each galaxy in the survey. This results in every galaxy with an uncertain redshift being assigned multiple "recovered" redshifts from which a new redshift estimate can be determined. An earlier paper applied SORT to a mock Sloan Digital Sky Survey at $z \lesssim$ 0.2 and accurately recovered the two-point correlation function on scales $\gtrsim$4 $h^{-1}$Mpc. In this paper, we test the performance of SORT in surveys spanning the redshift range 0.75$<z<$2.25. We used two mock surveys extracted from the Small MultiDark-Planck and Bolshoi-Planck N-body simulations with dark matter haloes that were populated by the Santa Cruz semi-analytic model. We find that SORT is able to improve redshift estimates and recover distinctive large-scale features of the cosmic web. Further, it provides unbiased estimates of the redshift-space two-point correlation function $\xi(s)$ on scales $\gtrsim$2.5 $h^{-1}$Mpc, as well as local densities in regions of average or higher density. This may allow improved understanding of how galaxy properties relate to their local environments.

Derivation of physical properties of galaxies using spectral energy distribution (SED) fitting is a powerful method, but can suffer from various systematics arising from model assumptions. Previously, such biases were mostly studied for individual galaxies. Here we specifically study a large sample (9000) of composite (merged) galaxies extracted from GALEX-SDSS-WISE Legacy Catalog (GSWLC) in order to check if SED fitting accurately determines the properties (stellar mass and star formation rate) of an unresolved merged galaxy. We carry out the assessment by artificially merging the galaxies that could merge in reality and running UV/optical SED fitting on them as well as on pre-merger individual galaxies (to establish ground truth). With this simple approach we know what mass and SFR to ideally expect for the merger. No biases are found when comparing the stellar mass, star formation rate (SFR), and specific star formation rate (sSFR) values from the SED fitting to the combined values from the individual galaxies, although there are tails in the SFR and sSFR merger comparisons with low (s)SFR. Furthermore, no significant residuals were found when focusing on galaxies with different mass ratios (major vs.\ minor mergers), different contrasts in SFR, or dust content. We conclude that the SED fitting of post-mergers does not suffer from significant biases that would arise from the fact that the merger is a composite of two galaxies with potentially very different star formation histories and with different dust properties.

Feedback from active galactic nuclei (AGN) is believed to be the most promising solution to the cooling flow problem in cool-core clusters, though how exactly the jet energy is transformed into heat is a subject of debate. Dissipation of sound waves is considered as one of the possible heating mechanisms; however, its relative contribution to heating remains unclear. To estimate the energy budget for heating, we perform 3D hydrodynamic simulations of AGN jet injections in a Perseus-like cluster and quantify the amount of energy stored in the forms of weak shocks and waves. We find that, for a single jet injection with typical parameters in cool-core clusters, $\sim9\%$ of the total jet energy is stored in compressional waves (including both shocks and waves). However, due to the destructive effects among randomly phased waves as well as the dissipation of shock energies, in our simulations including self-regulated AGN feedback, no more than $3\%$ of the total injected energy goes into compressional waves. We further separate the energy contribution from shocks and waves and find that, for a single outburst, the shocks can only contribute to $\sim20-30\%$ of the total compressional energy in the inner radii, and they dissipate very quickly as they travel outward. However, because of the repeated generation of shocks by multiple AGN outbursts, in the self-regulated case shocks completely dominate over sound waves in the inner region and can still provide $\sim40-50\%$ of the total compressional energy even at outer radii. Our results suggest that the production of sound waves is not as efficient as what was found in previous single-outburst simulations, and thus sound wave dissipation may be a subdominant source of heating in cool-core clusters.

In this paper we study scenarios of the super-Eddington accretion onto black holes at high redshifts z > 10, which are expected to be seeds to evolve to supermassive black holes until redshift z ~ 7. Such an accretion disk inevitably emitted high-energy photons which had heated the cosmological plasma of the inter-galactic medium continuously from high redshifts. In this case, cosmic history of cosmological gas temperature is modified, by which the absorption feature of the cosmological 21 cm lines are suppressed. By comparing theoretical predictions of the 21cm line absorption with the observational data, we can obtain a cosmological upper bound on the mass-accretion rate as a function of the seed BH masses. In order to realize $M_{\rm BH} \sim 10^9 M_{\odot}$ at $z \sim 7$ by a continuous mass-accretion on to a seed BH, to be consistent with the cosmological 21cm line absorption at $z \sim 17$, we obtained an upper bound on the initial mass of the seed BH to be $M_{\rm BH, ini} \lesssim 10^2 M_{\odot}$ ($M_{\rm BH, ini} \lesssim 10^6 M_{\odot}$) for a seed BH with its comoving number density $n_{\rm seed,0} \sim 10^{-3} {\rm Mpc}^{-3}$ ($n_{\rm seed,0} \sim 10^{-7} {\rm Mpc}^{-3}$).

In the optically thin regime, the intensity ratio of the two Si IV resonance lines (1394 and 1403 \AA\ ) are theoretically the same as the ratio of their oscillator strengths, which is exactly 2. Here, we study the ratio of the integrated intensity of the Si IV lines ($R=\int I_{1394}(\lambda)\mathrm{d}\lambda/\int I_{1403}(\lambda)\mathrm{d}\lambda$) and the ratio of intensity at each wavelength point ($r(\Delta\lambda)=I_{1394}(\Delta\lambda)/I_{1403}(\Delta\lambda)$) in two solar flares observed by the Interface Region Imaging Spectrograph. We find that at flare ribbons, the ratio $R$ ranges from 1.8 to 2.3 and would generally decrease when the ribbons sweep across the slit position. Besides, the distribution of $r(\Delta\lambda)$ shows a descending trend from the blue wing to the red wing. In loop cases, the Si IV line presents a wide profile with a central reversal. The ratio $R$ deviates little from 2, but the ratio $r(\Delta\lambda)$ can vary from 1.3 near the line center to greater than 2 in the line wings. Hence we conclude that in flare conditions, the ratio $r(\Delta\lambda)$ varies across the line, due to the variation of the opacity at the line center and line wings. We notice that, although the ratio $r(\Delta\lambda)$ could present a value which deviates from 2 as a result of the opacity effect near the line center, the ratio $R$ is still close to 2. Therefore, caution should be taken when using the ratio of the integrated intensity of the Si IV lines to diagnose the opacity effect.

Recently it was found from Cassini data that the mean recession speed of Titan from Saturn is $v=11.3\pm 2.0$ cm/yr which corresponds to a tidal quality factor of Saturn $Q\cong 100$ while the standard estimate yields $Q\ge 6\cdot 10^4$. It was assumed that such a large speed $v$ is due to a resonance locking mechanism of five inner mid-sized moons of Saturn. In this paper, we show that an essential part of $v$ may come from a local Hubble expansion, where the Hubble-Lema\^{\i}tre constant $H_0$ recalculated to the Saturn-Titan distance $D$ is 8.15 cm/(yr $D$). Our hypothesis is based on many other observations showing a slight expansion of the Solar system and also of our Galaxy at a rate comparable with $H_0$. We demonstrate that the large disproportion in estimating the $Q$ factor can be just caused by the local expansion effect.

We utilize the ALMA-MaNGA QUEnch and STar formation (ALMaQUEST) survey to investigate the kpc-scale scaling relations, presented as the resolved star forming main sequence (rSFMS: $\Sigma_{\rm SFR}$ vs. $\Sigma_{*}$), the resolved Schmidt-Kennicutt relation (rSK: $\Sigma_{\rm SFR}$ vs. $\Sigma_{\rm H_{2}}$), and the resolved molecular gas main sequence (rMGMS: $\Sigma_{\rm H_{2}}$ vs. $\Sigma_{*}$), for 11478 star-forming and 1414 retired spaxels (oversampled by a factor of $\sim20$) located in 22 green valley (GV) and 12 main sequence (MS) galaxies. For a given galaxy type (MS or GV), the retired spaxels are found to be offset from the sequences formed by the star-forming spaxels on the rSFMS, rSK, and rMGMS planes, toward lower absolute values of sSFR, SFE, and $f_{\rm H_{2}}$ by $\sim$ 1.1, 0.6, and 0.5 dex. The scaling relations for GV galaxies are found to be distinct from that of the MS galaxies, even if the analyses are restricted to the star-forming spaxels only. It is found that for star-forming spaxels, sSFR, SFE, and $f_{\rm H_{2}}$ in GV galaxies are reduced by $\sim$0.36, 0.14, and 0.21 dex, respectively, compared to those in MS galaxies. Therefore, the suppressed sSFR/SFE/$f_{\rm H_{2}}$ in GV galaxies are associated with not only an increased proportion of retired regions in GV galaxies but also a depletion of these quantities in star-forming regions. Finally, the reduction of SFE and $f_{\rm H_{2}}$ in GV galaxies relative to MS galaxies is seen in both bulge and disk regions (albeit with larger uncertainties), suggesting that statistically, quenching in the GV population may persist from the inner to the outer regions.

There is a growing concern about an impact of low-Earth-orbit (LEO) satellite constellations on ground-based astronomical observations, in particular, on wide-field surveys in the optical and infrared. The Zwicky Transient Facility (ZTF), thanks to the large field of view of its camera, provides an ideal setup to study the effects of LEO megaconstellations - such as SpaceX's Starlink - on astronomical surveys. Here, we analyze the archival ZTF observations collected between 2019 November and 2021 September and find 5301 satellite streaks that can be attributed to Starlink satellites. We find that the number of affected images is increasing with time as SpaceX deploys more and more satellites. Twilight observations are particularly affected - a fraction of streaked images taken during twilight has increased from less than 0.5% in late 2019 to 18% in 2021 August. We estimate that once the size of the Starlink constellation reaches 10,000, essentially all ZTF images taken during twilight may be affected. However, despite the increase in satellite streaks observed during the analyzed period, the current science operations of ZTF are not yet strongly affected. We also find that redesigning Starlink satellites (by installing visors intended to block sunlight from reaching the satellite antennas to prevent reflection) reduces their brightness by a factor of 4.6 +/- 0.1 with respect to the original design in g, r, and i bands.

Context. Many well-known bright stars have been observed by the ongoing Transiting Exoplanet Survey Satellite (TESS) space mission. For several of them, these new data reveal previously unobserved variability, such as tidally perturbed pulsations in close binary stars. Aims. Using newly detected gravity-mode (g-mode) pulsations in V456 Cyg, we aim to determine the global stellar properties of this short-period eclipsing binary and evaluate the interaction between these pulsations and the tides. Methods. We model the binary orbit and determine the physical properties of the component stars using the TESS photometry and published spectroscopy. We then measure the pulsation frequencies from the residuals of the light curve fit using iterative prewhitening, and analyse them to determine the global asteroseismic stellar parameters. We evaluate the pulsation parameters as a function of the orbital phase. Results. We find that the pulsations belong to the secondary component of V456 Cyg and that this star likely has a uniform radial rotation profile, synchronous ($\nu_{\rm rot} = 1.113(14) \rm d^{-1}$) with the binary orbit ($\nu_{\rm orb} = 1.122091(8) \rm d^{-1}$). The observed g modes are amplified by almost a factor three in the stellar hemisphere facing the primary. We present evidence that this is caused by tidal perturbation of the pulsations, with the mode coupling being strongly affected. Conclusions. V456 Cyg is only the second object for which tidally perturbed high-order g-mode pulsations are identified, after $\pi^5$ Ori. This opens up new opportunities for tidal g-mode asteroseismology, as it demonstrates another avenue in which g modes and tides can interact with each other.

The recently discovered near-Earth asteroid 2021 PH27 has the shortest orbital period of all known asteroids. It cannot be excluded that 2021 PH27 is also an active asteroid, as (3200) Phaethon. We intend to estimate the consequences of this hypothesis, although testing is difficult with ground-based observations during perihelion passages, due to low solar elongation. Assuming a surface activity similar to that of Phaethon, an increase in brightness of about 1.4 mag can be estimated. Since it is an asteroid with a MOID of $0.014660 \pm 0.000034$ AU with Venus, 2021 PH27 could be the equivalent of Phaethon for the Earth and be the progenitor body of a venusian meteor shower. A good opportunity to observe the hypothetical fireballs in Venus's atmosphere will take place on the days around Jun 07, 2023, when Venus will pass at the minimum distance from the nominal orbit of 2021 PH27. Another favorable date will be Jul 05, 2026. Finally, on Mar 28, 2022 the asteroid will also be at the maximum Sun elongation of about $52.3^{\circ}$ and at the aphelion of its orbit, the most favorable configuration to characterize it from the physical point of view with photometric, polarimetric and spectroscopic observations.

Measuring the redshift of active galactic nuclei (AGNs) requires the use of time-consuming and expensive spectroscopic analysis. However, obtaining redshift measurements of AGNs is crucial as it can enable AGN population studies, provide insight into the star formation rate, the luminosity function, and the density rate evolution. Hence, there is a requirement for alternative redshift measurement techniques. In this project, we aim to use the Fermi gamma-ray space telescope's 4LAC Data Release (DR2) catalog to train a machine learning model capable of predicting the redshift reliably. In addition, this project aims at improving and extending with the new 4LAC Catalog the predictive capabilities of the machine learning (ML) methodology published in Dainotti et al. (2021). Furthermore, we implement feature engineering to expand the parameter space and a bias correction technique to our final results. This study uses additional machine learning techniques inside the ensemble method, the SuperLearner, previously used in Dainotti et al.(2021). Additionally, we also test a novel ML model called Sorted L-One Penalized Estimation (SLOPE). Using these methods we provide a catalog of estimated redshift values for those AGNs that do not have a spectroscopic redshift measurement. These estimates can serve as a redshift reference for the community to verify as updated Fermi catalogs are released with more redshift measurements.

We report the discoveries of a nuclear ring of diameter 10$\arcsec$ ($\sim$1.5 kpc) and a potential low luminosity active galactic nucleus (LLAGN) in the radio continuum emission map of the edge-on barred spiral galaxy NGC~5792. These discoveries are based on the Continuum Halos in Nearby Galaxies - an Expanded Very Large Array (VLA) Survey, as well as subsequent VLA observations of sub-arcsecond resolution. Using a mixture of H$\alpha$ and 24 $\mu$m calibration, we disentangle the thermal and non-thermal radio emission of the nuclear region, and derive a star formation rate (SFR) of $\sim 0.4~M_{\sun}$ yr$^{-1}$. We find that the nuclear ring is dominated by non-thermal synchrotron emission. The synchrotron-based SFR is about three times of the mixture-based SFR. This result indicates that the nuclear ring underwent more intense star-forming activity in the past, and now its star formation is in the low state. The sub-arcsecond VLA images resolve six individual knots on the nuclear ring. The equipartition magnetic field strength $B_{\rm eq}$ of the knots varies from 77 to 88 $\mu$G. The radio ring surrounds a point-like faint radio core of $S_{\rm 6GHz}=(16\pm4)$ $\mu$Jy with polarized lobes at the center of NGC~5792, which suggests an LLAGN with an Eddington ratio $\sim10^{-5}$. This radio nuclear ring is reminiscent of the Central Molecular Zone (CMZ) of the Galaxy. Both of them consist of a nuclear ring and LLAGN.

Metal-deficient stars are important tracers for understanding the early formation of the Galaxy. Recent large-scale surveys with both photometric and spectroscopic data have reported an increasing number of metal-deficient stars whose kinematic features are consistent with those of the disk stellar populations. We report the discovery of an RR~Lyrae variable (hereafter RRL) that is located within the thick disk and has an orbit consistent with the thick-disk kinematics. Our target RRL (HD 331986) is located at around 1 kpc from the Sun and, with V=11.3, is among the 130 brightest RRLs known so far. However, this object was scarcely studied because it is in the midplane of the Galaxy, the Galactic latitude around -1 deg. Its near-infrared spectrum (0.91-1.32 micron) shows no absorption line except hydrogen lines of the Paschen series, suggesting [Fe/H] less than -2.5. It is the most metal-deficient RRL, at least, among the RRLs whose orbits are consistent with the disk kinematics, although we cannot determine to which of the disk and the halo it belongs. This unique RRL would provide us with essential clues for studying the early formation of stars in the inner Galaxy with further investigations, including high-resolution optical spectroscopy.

In this paper, we discuss the influence of the gravitational darkening effect on the emergent spectrum of a fast-rotating, flattened neutron star. Model atmosphere codes always calculate spectra of emergent intensities and fluxes emitted from the unit surface on the star in plane-parallel geometry. Here we took a step beyond that and calculated a small sample grid of theoretical spectra integrated over the distorted surface of a sample rotating neutron star seen by a distant observer at various inclination angles. We assumed parameters like two dimensionless angular velocities $\bar{\Omega}^2=0.30$ and 0.60, the effective temperature of a nonrotating star $T_{\rm eff}=2.20\times 10^7\,$K, the logarithm of the surface gravity of a spherical star $\log(g)=14.40$ (cgs), and inclination angles from $i=0^\circ$ to $i=90^\circ$ with step $\Delta i=10^\circ$. We assumed that the atmosphere consists of a mixture of hydrogen and helium with $M_{\rm H}=0.70$ and $M_{\rm He}=0.30$. At each point on the neutron star surface, we calculated true intensities for local values of parameters ($T_{\rm eff}$ and $\log(g)$), and these monochromatic intensities are next integrated over the whole surface to obtain the emergent spectrum. In this paper, we compute for the first time theoretical spectra of the fast-rotating neutron star. Our work clearly shows that the gravitational darkening effect strongly influences the spectrum and should be included in realistic models of the atmospheres of rotating neutron stars.

This article is an attempt to review 50 years of high-energy cosmic particle physics at the Institute for Nuclear Research of the Russian Academy of Sciences. It is written by an outsider whose scientific career, to a large part, was formed by collaborating with INR scientists in the late 1980s and 1990s. The review covers the fields of cosmic-ray, gamma-ray and high-energy neutrino physics. The main focus will be on INR's large infrastructures in the Baksan Valley and at Lake Baikal. Research at these facilities is flanked by participation in top experiments at different places around the world, recently the Telescope Array in the USA and the LHAASO detector in China.

While the slope of the dust attenuation curve ($\delta$) is found to correlate with effective dust attenuation ($A_V$) as obtained through spectral energy distribution (SED) fitting, it remains unknown how the fitting degeneracies shape this relation. We examine the degeneracy effects by fitting SEDs of a sample of local star-forming galaxies (SFGs) selected from the Galaxy And Mass Assembly survey, in conjunction with mock galaxy SEDs of known attenuation parameters. A well-designed declining starburst star formation history is adopted to generate model SED templates with intrinsic UV slope ($\beta_0$) spanning over a reasonably wide range. The best-fitting $\beta_0$ for our sample SFGs shows a wide coverage, dramatically differing from the limited range of $\beta_0<-2.2$ for a starburst of constant star formation. Our results show that strong degeneracies between $\beta_0$, $\delta$, and $A_V$ in the SED fitting induce systematic biases leading to a false $A_V$--$\delta$ correlation. Our simulation tests reveal that this relationship can be well reproduced even when a flat $A_V$--$\delta$ relation is taken to build the input model galaxy SEDs. The variations in best-fitting $\delta$ are dominated by the fitting errors. We show that assuming a starburst with constant star formation in SED fitting will result in a steeper attenuation curve, smaller degeneracy errors, and a stronger $A_V$--$\delta$ relation. Our findings confirm that the $A_V$--$\delta$ relation obtained through SED fitting is likely driven by the systematic biases induced by the fitting degeneracies between $\beta_0$, $\delta$, and $A_V$.

We propose a new, efficient multi-scale method to decompose a map (or signal in general) into components maps that contain structures of different sizes. In the widely-used wave transform, artifacts containing negative values arise around regions with sharp transitions due to the application of band-limited filters. In our approach, the decomposition is achieved by solving a modified, non-linear version of the diffusion equation. This is inspired by the anisotropic diffusion methods, which establish the link between image filtering and partial differential equations. In our case, the artifact issue is addressed where the positivity of the decomposed images is guaranteed. Our new method is particularly suitable for signals which contain localized, non-linear features, as typical of astronomical observations. It can be used to study the multi-scale structures of astronomical maps quantitatively and should be useful in observation-related tasks such as background removal. We thus propose a new measure called the ''scale spectrum'', which describes how the image values distribute among different components in the scale space, to describe maps. The method allows for input arrays of an arbitrary number of dimensions, and a python3 implementation of the algorithms is included in the Appendix and available at https://gxli.github.io/Constrained-Diffusion-Decomposition/.

Different relationships between the H$\alpha$ and Ca II chromospheric emissions have been reported in solar-type stars. In particular, the time-series of emissions in these two lines are clearly anti-correlated for a few percent of the stars, contrary to what is observed on the Sun. Our objective is to characterise these relationships in more detail using complementary criteria, and to constrain the properties of filaments and plages that are necessary to explain the observations. We analysed the average level and variability of the H$\alpha$ and Ca II emission for 441 F-G-K stars, paying particular attention to their (anti-)correlations on both short and long timescales. We also computed synthetic H$\alpha$ and Ca II time-series for different assumptions of plage and filament properties and compared them with the observations. We were not able to find plage properties that, alone, are sufficient to reproduce the observations at all timescales simultaneously, even when allowing different H$\alpha$ and Ca II emission relationships for different stars. We also specified the complex and surprising relationship between the average activity levels of both indexes, in particular for low-activity stars.We conclude that plages alone are unlikely to explain the observed variety of relationships between Ca II and H$\alpha$ emission, and that the presence of other phenomena like filaments may help to reconcile the models with observations.

Relativistic radiation mediated shocks (RRMS) likely form in prodigious cosmic explosions. The structure and emission of such shocks is regulated by copious production of electron-positron pairs inside the shock transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the nonlinear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced transport equation for the particles via pitch angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to nonadiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of nonthermal distributions, arising from the microturbulence, can affect the observed shock breakout signal in ways unaccounted for by current single-fluid models.

We report on NICER X-ray monitoring of the magnetar SGR 1830-0645 covering 223 days following its October 2020 outburst, as well as Chandra and radio observations. We present the most accurate spin ephemerides of the source so far: $\nu=0.096008680(2)$~Hz, $\dot{\nu}=-6.2(1)\times10^{-14}$~Hz~s$^{-1}$, and a significant second and third frequency derivative terms indicative of non-negligible timing noise. The phase-averaged 0.8--7~keV spectrum is well fit with a double-blackbody (BB) model throughout the campaign. The BB temperatures remain constant at 0.46 and 1.2 keV. The areas and flux of each component decreased by a factor of 6, initially through a steep decay trend lasting about 46 days followed by a shallow long-term one. The pulse shape in the same energy range is initially complex, exhibiting three distinct peaks, yet with clear continuous evolution throughout the outburst towards a simpler, single-pulse shape. The rms pulsed fraction is high and increases from about 40% to 50%. We find no dependence of pulse shape or fraction on energy. These results suggest that multiple hotspots, possibly possessing temperature gradients, emerged at outburst-onset, and shrank as the outburst decayed. We detect 84 faint bursts with \nicer, having a strong preference for occurring close to the surface emission pulse maximum the first time this phenomenon is detected in such a large burst sample. This likely implies a very low altitude for the burst emission region, and a triggering mechanism connected to the surface active zone. Finally, our radio observations at several epochs and multiple frequencies reveal no evidence of pulsed or burst-like radio emission.

Observationally inferred crystalline abundance in silicates in comets, which should have been formed in the outer region of a protoplanetary disk, is relatively high (~ 10-60%), although crystalline silicates would be formed by annealing of amorphous precursors in the disk inner region. In order to quantitatively address this puzzle, we have performed Monte Carlo simulation of advection/diffusion of silicate particles in a turbulent disk, in the setting based on pebble accretion model: pebbles consisting of many small amorphous silicates embedded in icy mantle are formed in the disk outer region, silicate particles are released at the snow line, crystalline silicate particles are produced at the annealing line, the silicate particles diffused beyond the snow line, and they eventually stick to drifting pebbles to come back to the snow line. In a simple case without the sticking and with a steady pebble flux, we show through the simulations and analytical arguments that crystalline components in silicate materials beyond the snow line is robustly and uniformly ~ 5%. On the other hand, in a more realistic case with the sticking and with a decaying pebble flux, the crystalline abundance is raised up to ~ 20-25%, depending on the ratio of decay and diffusion timescales. This abundance is consistent with the observations. In this investigation, we assume a simple steady accretion disk. The simulations coupled with the disk evolution is needed for more detailed comparison with observed data.

Numerous magnetic hot stars exhibit gyrosynchrotron radio emission. The source electrons were previously thought to be accelerated to relativistic velocities in the current sheet formed in the middle magnetosphere by the wind opening magnetic field lines. However, a lack of dependence of radio luminosity on the wind power, and a strong dependence on rotation, has recently challenged this paradigm. We have collected all radio measurements of magnetic early-type stars available in the literature. When constraints on the magnetic field and/or the rotational period are not available, we have determined these using previously unpublished spectropolarimetric and photometric data. The result is the largest sample of magnetic stars with radio observations that has yet been analyzed: 131 stars with rotational and magnetic constraints, of which 50 are radio-bright. We confirm an obvious dependence of gyrosynchrotron radiation on rotation, and furthermore find that accounting for rotation neatly separates stars with and without detected radio emission. There is a close correlation between H$\alpha$ emission strength and radio luminosity. These factors suggest that radio emission may be explained by the same mechanism responsible for H$\alpha$ emission from centrifugal magnetospheres, i.e. centrifugal breakout (CBO), however, whereas the H$\alpha$-emitting magnetosphere probes the cool plasma before breakout, radio emission is a consequence of electrons accelerated in centrifugally-driven magnetic reconnection.

Magnetars, isolated neutron stars with magnetic field strengths typically $\gtrsim10^{14}$~G, exhibit distinctive months-long outburst epochs during which strong evolution of soft X-ray pulse profiles, along with nonthermal magnetospheric emission components, is often observed. Using near-daily NICER observations of the magnetar SGR 1830-0645 during the first 37 days of a recent outburst decay, a pulse peak migration in phase is clearly observed, transforming the pulse shape from an initially triple-peaked to a single-peaked profile. Such peak merging has not been seen before for a magnetar. Our high-resolution phase-resolved spectroscopic analysis reveals no significant evolution of temperature despite the complex initial pulse shape. Yet the inferred surface hot spots shrink during the peak migration and outburst decay. We suggest two possible origins for this evolution. For internal heating of the surface, tectonic motion of the crust may be its underlying cause. The inferred speed of this crustal motion is $\lesssim100$~m~day$^{-1}$, constraining the density of the driving region to $\rho\sim10^{10}$~g~cm$^{-3}$, at a depth of $\sim200$~m. Alternatively, the hot spots could be heated by particle bombardment from a twisted magnetosphere possessing flux tubes or ropes, somewhat resembling solar coronal loops, that untwist and dissipate on the 30-40~day timescale. The peak migration may then be due to a combination of field-line footpoint motion (necessarily driven by crustal motion) and evolving surface radiation beaming. These novel dataset paints a vivid picture of the dynamics associated with magnetar outbursts, yet it also highlights the need for a more generic theoretical picture where magnetosphere and crust are considered in tandem.

Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as well as archival NIR/FIR, cm continuum, CO(3-2) and HI data are folded into the analysis. Results: While the overall [CII] morphology follows the general ionized gas, the channel maps show multiple bubble-like structures with sizes on the order of ~80-100" (~1.0-1.28pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the region, the other bubble-morphologies may also be due to the intrinsically porous structure of the HII region. An analysis of the expansion velocities around 10km s^{-1} indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The most blue-shifted [CII] component has barely any molecular or atomic counterparts. At the interface to the molecular cloud, we find a typical photon-dominated region (PDR) with a bar-shape. Ionized, atomic and molecular carbon show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized form with atomic and molecular masses of ~0.45+-0.1M_{\odot} and ~1.2+-0.1M_{\odot}, respectively, compared to the ionized carbon in the range of 3.6-9.7M_{\odot}. Conclusions: The NGC7538 HII region exhibits a diverse set of sub-structures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of HII regions (e.g., Orion Veil, RCW120, RCW49), bubble-shape morphologies revealed in [CII] emission, indicative of expanding shells, are recurring structures of PDRs.

In the present work we utilize ACE/SWICS in-situ measurements of the properties of the solar wind outside ICMEs in order to determine whether, and to what extent are the solar wind properties affected by the solar cycle. We focus on proton temperatures and densities, ion temperatures and differential speeds, charge state distributions and both relative and absolute elemental abundances. We carry out this work dividing the wind in velocity bins to investigate how winds at different speeds react to the solar cycle. We also repeat this study, when possible, to the subset of SWICS measurements less affected by Coulomb collisions. We find that with the only exception of differential speeds (for which we do not have enough measurements) all wind properties change as a function of the solar cycle. Our results point towards a scenario where both the slow and fast solar wind are accelerated by waves, but originate from different sources (open/closed magnetic structures for the fast/slow wind, respectively) whose relative contribution changes along the solar cycle. We also find that the signatures of heating and acceleration on one side, and of the FIP effect on the other, indicate that wave-based plasma heating, acceleration and fractionation remain active throughout the solar cycle, but decrease their effectiveness in all winds, although the slow wind is much affected than the fast one.

We present a first-principles model of pitch-angle and energy distribution function evolution as particles are sequentially accelerated by multiple flare magnetic islands. Data from magnetohydrodynamic (MHD) simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. Magnetic islands, which are created by sporadic reconnection at the self-consistently formed flare current sheet, contract and accelerate the particles. The particle distributions are evolved using rules derived in our previous work. In this investigation, we assume that a prescribed fraction of particles sequentially "hops" to another accelerator and receives an additional boost in energy and anisotropy. This sequential process generates particle number spectra that obey an approximate power law at mid-range energies and presents low- and high-energy breaks. We analyze these spectral regions as functions of the model parameters. We also present a fully analytic method for forming and interpreting such spectra, independent of the sequential acceleration model. The method requires only a few constrained physical parameters, such as the percentage of particles transferred between accelerators, the energy gain in each accelerator, and the number of accelerators visited. Our investigation seeks to bridge the gap between MHD and kinetic regimes by combining global simulations and analytic kinetic theory. The model reproduces and explains key characteristics of observed flare hard X-ray spectra as well as the underlying properties of the accelerated particles. Our analytic model provides tools to interpret high-energy observations for missions and telescopes, such as RHESSI, FOXSI, NuSTAR, Solar Orbiter, EOVSA, and future high-energy missions.

The Kepler supernova remnant (SNR) is the only historic supernova remnant lacking a detection at GeV and TeV energies which probe particle acceleration. A recent analysis of Fermi-LAT data reported a likely GeV gamma-ray candidate in the direction of the SNR. Using approximately the same dataset but with an optimized analysis configuration, we confirm the gamma-ray candidate to a solid $>6\sigma$ detection and report a spectral index of $2.14 \pm 0.12_{\rm stat} \pm 0.15_{\rm syst}$ for an energy flux above 100 MeV of $(3.1 \pm 0.6_{\rm stat} \pm 0.3_{\rm syst}) \times 10^{-12}$ erg~cm$^{-2}$~s$^{-1}$. The gamma-ray excess is not significantly extended and is fully compatible with the radio, infrared or X-ray spatial distribution of the SNR. We successfully characterized this multi-wavelength emission with a model in which accelerated particles interact with the dense circumstellar material in the North-West portion of the SNR and radiate GeV gamma-rays through $\pi^{o}$ decay. The X-ray synchrotron and inverse-Compton (IC) emission mostly stem from the fast shocks in the southern regions with a magnetic field B$\sim$100 $\mu$G or higher. Depending on the exact magnetic field amplitude, the TeV emission could arise from either the South region (IC dominated) or the interaction region ($\pi^{o}$ decay dominated).

Context. Several methods have been proposed to build 3D extinction maps of the Milky Way (MW), most often based on Bayesian approaches. Although some studies employed machine learning (ML) methods in part of their procedure, or to specific targets, no 3D extinction map of a large volume of the MW solely based on a Neural Network method has been reported so far. Aims. We aim to apply deep learning as a solution to build 3D extinction maps of the MW. Methods. We built a convolutional neural network (CNN) using the CIANNA framework, and trained it with synthetic 2MASS data. We used the Besan\c{c}on Galaxy model to generate mock star catalogs, and 1D Gaussian random fields to simulate the extinction profiles. From these data we computed color-magnitude diagrams (CMDs) to train the network, using the corresponding extinction profiles as targets. A forward pass with observed 2MASS CMDs provided extinction profile estimates for a grid of lines of sight. Results. We trained our network with data simulating lines of sight in the area of the Carina spiral arm tangent and obtained a 3D extinction map for a large sector in this region ($l = 257 - 303$ deg, $|b| \le 5$ deg), with distance and angular resolutions of $100$ pc and $30$ arcmin, respectively, and reaching up to $\sim 10$ kpc. Although each sightline is computed independently in the forward phase, the so-called fingers-of-God artifacts are weaker than in many other 3D extinction maps. We found that our CNN was efficient in taking advantage of redundancy across lines of sight, enabling us to train it with only 9 sightlines simultaneously to build the whole map. Conclusions. We found deep learning to be a reliable approach to produce 3D extinction maps from large surveys. With this methodology, we expect to easily combine heterogeneous surveys without cross-matching, and therefore to exploit several surveys in a complementary fashion.

We model the response of spherical, non-rotating Milky Way (MW) dark matter and stellar halos to the Large Magellanic Cloud (LMC) using the matrix method of linear response theory. Our computations reproduce the main features of the dark halo response from simulations. We show that these features can be well separated by a harmonic decomposition: the large scale over/underdensity in the halo (associated with its reflex motion) corresponds to the $\ell=1$ terms, and the local overdensity to the $\ell\geq2$ multipoles. Moreover, the dark halo response is largely dominated by the first order 'forcing' term, with little influence from self-gravity. This makes it difficult to constrain the underlying velocity distribution of the dark halo using the observed response of the stellar halo, but it allows us to investigate the response of stellar halo models with various velocity anisotropies: a tangential (respectively radial) halo produces a shallower (respectively stronger) response. We also show that only the local wake is responsible for these variations, the reflex motion being solely dependent on the MW potential. Therefore, we identify the structure (orientation and winding) of the in-plane quadrupolar ($m=2$) response as a potentially good probe of the stellar halo anisotropy. Finally, our method allows us to tentatively relate the wake strength and shape to resonant effects: the strong radial response could be associated with the inner Lindblad resonance, and the weak tangential one with corotation.

We investigate the possible entropic nature of the force responsible for the discrepancy between the observed galactic rotation curves and those expected from the distribution of visible matter in the galaxy. Observations from the Spitzer Photometry and Accurate Rotation Curves (SPARC) data base are used to study the adequacy of the proposed models. A concrete model derived from a simple solution of the Fokker-Planck equation is used to fit SPARC data, resulting in strong agreement with observations compared to the popular dark matter NFW mass profile. We also show that correlations exist between a parameter of the proposed model and the flat velocity as well as the Luminosity at 3.6 $\mu$m of the sample of galaxies.

We determine the low-redshift X-ray luminosity function (XLF), active black hole mass function (BHMF), and Eddington-ratio distribution function (ERDF) for both unobscured (Type 1) and obscured (Type 2) active galactic nuclei (AGN) using the unprecedented spectroscopic completeness of the BAT AGN Spectroscopic Survey (BASS) data release 2. In addition to a straightforward 1/Vmax approach, we also compute the intrinsic distributions, accounting for sample truncation by employing a forward modeling approach to recover the observed BHMF and ERDF. As previous BHMFs and ERDFs have been robustly determined only for samples of bright, broad-line (Type 1) AGNs and/or quasars, ours is the first directly observationally constrained BHMF and ERDF of Type 2 AGN. We find that after accounting for all observational biases, the intrinsic ERDF of Type 2 AGN is significantly skewed towards lower Eddington ratios than the intrinsic ERDF of Type 1 AGN. This result supports the radiation-regulated unification scenario, in which radiation pressure dictates the geometry of the dusty obscuring structure around an AGN. Calculating the ERDFs in two separate mass bins, we verify that the derived shape is consistent, validating the assumption that the ERDF (shape) is mass independent. We report the local AGN duty cycle as a function of mass and Eddington ratio, by comparing the BASS active BHMF with the local mass function for all SMBH. We also present the log N-log S of Swift-BAT 70-month sources.

The integrated shear 3-point correlation function $\zeta_{\pm}$ is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function $\xi_{\pm}$ by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate $\zeta_{\pm}$ that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale $\zeta_{\pm}$ is dominated by the nonlinear matter bispectrum in the squeezed limit, which can be evaluated accurately using the nonlinear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale $\zeta_{\pm}$ measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the nonlinear matter power spectrum, thereby permitting to account for these astrophysical effects on $\zeta_{\pm}$ similarly to how they are currently accounted for on $\xi_{\pm}$. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to $\xi_{\pm}$, a combined $\xi_{\pm}\ \&\ \zeta_{\pm}$ analysis can lead to improvements of order $20-40\%$ on the constraints of cosmological parameters such as $\sigma_8$ or the dark energy equation of state parameter $w_0$. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These are encouraging results that motivate future works on the integrated shear 3-point correlation function towards applications to real survey data.

Cosmic domain walls are harmless, provided that their tension decreases with expansion of the Universe. This setup can be realized, if the scale of spontaneous symmetry breaking is induced dynamically through the interaction with hot primordial plasma. In that case, the domain wall tension can attain large values in the early Universe without any conflict with observations. Owing to the large initial tension, these topological defects may serve as a powerful source of gravitational waves. We make a preliminary estimate of the gravitational wave spectrum and argue that it is distinct from the spectrum produced by other sources, in particular by domain walls of a constant tension. The resulting gravitational wave signal is in the range accessible by Einstein Telescope, DECIGO, TianQin, LISA, IPTA, or SKA, if the field constituting the domain walls is very feebly coupled with hot primordial plasma and has tiny self-interactions. In particular, one can consider this field for the role of Dark Matter. We discuss various Dark Matter production mechanisms and properties of the emitted gravitational waves associated with them. We find that the conventional freeze-out and freeze-in mechanisms lead to large and perhaps unobservable frequency of gravitational waves. However, the Dark Matter production is also possible at the second order phase transition leading to the domain wall formation or at the inverse phase transition, when the domain walls get dissolved eventually. In both cases, there is essentially no lower bound on the frequency of emitted gravitational waves.

We study the coupled dark energy model constructed from the general conformal transformation in which the coefficient of the conformal transformation depends on both the scalar field and its kinetic term. Under this conformal transformation, the action for subclass of Degenerate Higher-Order Scalar-Tensor (DHOST) theories is related to the Einstein-Hilbert action. The evolution of the background universe has the scaling fixed point which corresponds to acceleration of the universe at late time. For the choices of parameters that make the late-time scaling point stable, the fixed point corresponding to $\phi$-matter-dominated-era ($\phi$MDE) is a saddle point, and the universe can evolve from radiation dominated epoch through $\phi$MDE before reaching the scaling point at late time with the cosmological parameters which satisfy the observational bound. During the $\phi$MDE, the effective equation of state parameter is slightly positive, so that one of possible mechanisms for alleviating the $H_0$ tension can be achieved. In this coupled dark energy model, the effective gravitational coupling for dark matter perturbations on small scales can be smaller than that in the $\Lambda$CDM model. Therefore the growth rate of the dark matter perturbations is suppressed compare with the $\Lambda$CDM model, which implies that the $\sigma_8$ tension could be alleviated.

Now that we know that Earth-like planets are ubiquitous in the universe, as well as that most of them are much older than the Earth, it is justified to ask to what extent evolutionary outcomes on other such planets are similar, or indeed commensurable, to the outcomes we perceive around us. In order to assess the degree of specialty or mediocrity of our trajectory of biospheric evolution, we need to take into account recent advances in theoretical astrobiology, in particular (i) establishing the history of habitable planets' formation in the Galaxy, and (ii) understanding the crucial importance of "Gaian" feedback loops and temporal windows for the interaction of early life with its physical environment. Hereby we consider an alternative macroevolutionary pathway that may result in tight functional integration of all sub-planetary ecosystems, eventually giving rise to a true superorganism at the biospheric level. The blueprint for a possible outcome of this scenario has been masterfully provided by the great Polish novelist Stanis{\l}aw Lem in his 1961 novel Solaris. In fact, Solaris offers such a persuasive and powerful case for an "extremely strong" Gaia hypothesis that it is, arguably, high time to investigate it in a discursive astrobiological and philosophical context. In addition to novel predictions in the domain of potentially detectable biosignatures, some additional cognitive and heuristic benefits of studying such extreme cases of functional integration are briefly discussed.

Recent studies have outlined the interest for the evaluation of transport coefficients in space plasmas, where the observed velocity distributions of plasma particles are conditioned not only by the binary collisions, e.g., at low energies, but also by the energisation of particles from their interaction with wave turbulence and fluctuations, generating the suprathermal Kappa-distributed populations. This paper provides a first estimate of the main transport coefficients based on regularised Kappa distributions (RKDs), which, unlike standard Kappa distributions (SKDs), enable macroscopic parameterisation without mathematical divergences or physical inconsistencies. All transport coefficients derived here, i.e., the diffusion and mobility coefficients, electric conductivity, thermoelectric coefficient and thermal conductivity, are finite and well defined for all values of $\kappa > 0$. Moreover, for low values of $\kappa$ (i.e., below the SKD poles), the transport coefficients can be orders of magnitudes higher than the corresponding Maxwellian limits, meaning that significant underestimations can be made if suprathermal electrons are ignored.

It is well-known that stars have the potential to be excellent dark matter detectors. Infalling dark matter that scatters within stars could lead to a range of observational signatures, including stellar heating, black hole formation, and modified heat transport. To make robust predictions for such phenomena, it is necessary to calculate the scattering rate for dark matter inside the star. As we show in this paper, for small enough momentum transfers, this requires taking into account collective effects within the dense stellar medium. These effects have been neglected in many previous treatments; we demonstrate how to incorporate them systematically, and show that they can parametrically enhance or suppress dark matter scattering rates depending on how dark matter couples to the Standard Model. We show that, as a result, collective effects can significantly modify the potential discovery or exclusion reach for observations of compact objects such as white dwarfs and neutron stars. While the effects are more pronounced for dark matter coupling through a light mediator, we show that even for dark matter coupling via a heavy mediator, scattering rates can differ by orders of magnitude from their naive values for dark matter masses <~ 100 MeV. We also illustrate how collective effects can be important for dark matter scattering in more dilute media, such as the Solar core. Our results demonstrate the need to systematically incorporate collective effects in a wide range of astroparticle contexts; to facilitate this, we provide expressions for in-medium self-energies for a variety of different media, which are applicable to many other processes of interest (such as particle production).

Since the initial discovery of gravitational waves in 2015, significant developments have been made towards waveform interpretation and estimation of compact binary source parameters. We present herein an implementation of the generalized precession parameter $\langle \chi_p \rangle$, which averages over all angular variations on the precession timescale, within the RIFT parameter estimation framework. Relative to the originally-proposed precession parameter $\chi_p$, which characterizes the single largest dynamical spin in a binary, this new parameter $\langle \chi_p \rangle$ has a unique domain $1 < \langle \chi_p \rangle < 2$, which is exclusive to binaries with two precessing spins. After reviewing the physical differences between these two parameters, we describe how $\langle \chi_p \rangle$ was implemented in RIFT and apply it to all 36 events from the second half of LIGO's third operating run (O3b). In O3b, ten events show significant amounts of precession $\langle \chi_p \rangle > 0.5$. Of particular interest is GW191109_010717; we show it has a $\sim28\%$ probability that the originating system necessarily contains two misaligned spins.

Parametric Instability (PI) is a phenomenon that results from resonant interactions between optical and acoustic modes of a laser cavity. This is problematic in gravitational wave interferometers where the high intra-cavity power and low mechanical loss mirror suspension systems create an environment where three mode PI will occur without intervention. We demonstrate a technique for real time imaging of the amplitude and phase of the optical modes of PI yielding the first ever images of this phenomenon which could form part of active control strategies for future detectors.

Friedmann-Lema\^itre-Robertson-Walker cosmology is examined from the point of view of gravitoelectromagnetism, in the approximation of spacetime regions small in comparison with the Hubble radius. The usual Lorentz gauge is not appropriate for this situation, while the Painlev\'e-Gullstrand gauge is rather natural. Several non-trivial features and differences with respect to ``standard'' asymptotically flat gravitoelectromagnetism are discussed.

The thermal production mechanism of dark matter is attractive and well-motivated by predictivity. A representative of this type of dark matter candidate is the canonical, weakly interacting massive particles. An alternative is semi-annihilating dark matter, which exhibits different phenomenological aspects from the former example. In this study, we constructed a model of dark matter semi-annihilating into a pair of anti-dark matter and a Majoron based on a global $U(1)_{B-L}$ symmetry, and show that semi-annihilation induces the core formation of dark matter halos, which can alleviate the so-called small-scale problems. In addition, the box-shaped spectrum of neutrinos was produced by the subsequent decay of the Majoron. This can be a distinctive signature of the dark matter in the model. We find a parameter space where the produced neutrinos can be detected by the future large-volume neutrino detector Hyper-Kamiokande. We also compared the dark matter scenario with the case of halo core formation by the strongly self-interacting dark matter.

The IBEX-Lo instrument on board the Interstellar Boundary Explorer (IBEX) mission samples interstellar neutral (ISN) helium atoms penetrating the heliosphere from the very local interstellar medium (VLISM). In this study, we analyze the IBEX-Lo ISN helium observations covering a complete solar cycle, from 2009 through 2020 using a comprehensive uncertainty analysis including statistical and systematic sources.W e employ the Warsaw Test Particle Model to simulate ISN helium fluxes at IBEX, which are subsequently compared with the observed count rate in the three lowest energy steps of IBEX-Lo. The $\chi^2$ analysis shows that the ISN helium flows from ecliptic $(\lambda,\beta)=(255.59^{\circ}\pm0.23^{\circ}, 5.14^{\circ}\pm0.08^{\circ})$, with speed $v_\text{HP}=25.86\pm0.21$ km s$^{-1}$ and temperature $T_\text{HP}=7450\pm140$ K at the heliopause. Accounting for gravitational attraction and elastic collisions, the ISN helium speed and temperature in the pristine VLISM far from the heliopause are $v_\text{VLISM}=25.9$ km s$^{-1}$ and $T_\text{VLISM}=6150$ K, respectively. The time evolution of the ISN helium fluxes at 1 au over 12 years suggests significant changes in the IBEX-Lo detection efficiency, higher ionization rates of ISN helium atoms in the heliosphere than assumed in the model, or an additional unaccounted signal source in the analyzed observations. Nevertheless, we do not find any indication of the evolution of the derived parameters of ISN helium over the period analyzed. Finally, we argue that the continued operation of IBEX-Lo to overlap with the Interstellar Mapping and Acceleration Probe (IMAP) will be pivotal in tracking possible physical changes in the VLISM.

We introduce a new approach to renormalize physical quantities in curved space-time by adiabatic subtraction. We use a comoving infrared cut-off in defining the adiabatic counterpart of the physical quantity under consideration, building on the fact that the adiabatic approximation is ill-defined in the infrared tail of the spectrum. We show how this infrared cut-off should be used to obtain a completely well-defined renormalization scheme and how it is fundamental to avoid unphysical divergences that can be generated by a pathological behavior of the adiabatic subtraction extended to the infrared tail. The infrared cut-off appears as a new degree of freedom introduced in the theory and its actual value has to be consistently fixed by a physical prescription. As an example, we show how such degree of freedom can be set to obtain the correct value of the conformal anomaly in the symptomatic case of an inflationary model with gauge fields coupled to a pseudo-scalar inflaton.