Papers for Tuesday, Jan 11 2022

Understanding the formation and growth of supermassive black holes (SMBHs) at high redshift represents a major challenge for theoretical models. In this work we investigate the early evolution of the first SMBHs by constraining their distribution in mass and luminosity at $z > 4$. In particular, we focus on the poorly explored low-mass end of the nuclear black hole (BH) distribution down to $z \simeq 4$, and explore its connection with the nature of the first BH seeds and the processes governing their mass growth. To this aim, we have developed CAT (Cosmic Archaeology Tool), a new semi-analytic model that describes the formation of the first stars and black holes in a self-consistent way and follows the co-evolution of nuclear BHs and their host galaxies for a representative population at $z > 4$. We find that current observational constraints favour models where the growth of BH seeds is Eddington limited and occurs at the Bondi-Hoyle-Lyttleton rate or where super-Eddington accretion occurs via a slim disk during gas rich galaxy mergers. The main difference between these two model variants lies at the low-end of the predicted mass and luminosity functions at $4 \le z \le 6$, where a clear gap appears in the first model, reflecting the stunted growth of light BH seeds formed as remnants of the first stars. Detecting this signature will be extremely challenging even for the future generation of space observatories, such as JWST, Athena and Lynx.

We present an absorption line study of the physical and chemical properties of the Leo HI Ring and the Leo I Group as traced by 11 quasar sightlines spread over a 600 kpc X 800 kpc region. Using HST/COS G130/G160 archival observations as constraints, we couple cloud-by-cloud, multiphase, Bayesian ionization modeling with galaxy property information to determine the plausible origin of the absorbing gas along these sightlines. We search for absorption in the range 600 km/s - 1400 km/s consistent with the kinematics of the Leo Ring/Group. We find absorption plausibly associated with the Leo Ring towards five sightlines. Along three other sightlines, we find absorption likely to be associated with individual galaxies, intragroup gas, and/or large-scale filamentary structure. The absorption along these five sightlines is stronger in metal lines than expected from individual galaxies, indicative of multiple contributions, and of the complex kinematics of the region. We also identify three sightlines within a 7-degree X 6-degree field around the Leo Ring, along which we do not find any absorption. We find that the metallicities associated with the Leo Ring are generally high, with values between solar and several times solar. The inferred high metallicities are consistent with the origin of the ring as tidal debris from a major galaxy merger.

We present a fast methodology to produce mock observations of the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects based on dark matter only $N$-body simulations coupled with the analytic intra-cluster medium model. The methods employ two different approaches: halo-based pasting (HP) and particle-based pasting (PP). The former pastes gas density and pressure onto halos and requires only halo catalogue, and the latter also considers contribution from field particles as well, i.e., particles which do not belong to any halos, and thus utilises the full particle information. Therefore, the PP algorithm incorporates secondary effects beyond the HP algorithm: asphericity of halos and contribution from diffuse gas. In particular, such diffuse component is the dominant source of the kinetic SZ effect. As validation of our methods, we have produced 108 all-sky maps with HP and 108 flat-sky maps which cover $5 \times 5 \, \mathrm{deg}^2$ with both of HP and PP. Our method can produce a mock map within a few hours even for all-sky coverage with parallel computational environment. The resultant power spectra of these maps are consistent with theoretical predictions. We discuss the utility of baryon pasted mock SZ maps for estimating the covariance matrix of cross-correlation between SZ effects and other large-scale structure probes as well as modelling the selection and projection effects for cluster cosmology.

We analyze the angular power spectrum (APS) of the unresolved gamma-ray background (UGRB) emission and combine it with the measured properties of the resolved gamma-ray sources of the \Fermi-LAT 4FGL catalog. Our goals are to dissect the composition of the gamma-ray sky and to establish the relevance of different classes of source populations of active galactic nuclei in determining the observed size of the UGRB anisotropy, especially at low energies. We find that, under physical assumptions for the spectral energy dispersion, two populations are required to fit APS data, namely flat spectrum radio quasars (FSRQs) at low energies and BL Lacs (BLLs) at higher energies. The inferred luminosity functions agree well with the extrapolation of the FSRQ and BLL ones obtained from the 4FLG catalog. We use these luminosity functions to calculate the UGRB intensity from blazars, finding a contribution of 20\% at 1GeV and 30\% above 10 GeV. Finally, bounds on an additional gamma-ray emission due to annihilating dark matter are derived.

We present the $\textit{AUTOmated Photometry Of Transients}$ (AutoPhOT) package, a novel automated pipeline that is designed for rapid, publication-quality photometry of transients. AutoPhOT is built from the ground up using Python 3 - with no dependencies on legacy software. Capabilities of AutoPhOT include aperture and PSF-fitting photometry, template subtraction, and calculation of limiting magnitudes through artificial source injection. AutoPhOT is also capable of calibrating photometry against either survey catalogs (e.g. SDSS, PanSTARRS), or using a custom set of local photometric standards. We demonstrate the ability of AutoPhOT to reproduce lightcurves found in the published literature. AutoPhOT's ability to recover source fluxes is consistent with commonly used software e.g. DAOPHOT, using both aperture and PSF photometry. We also demonstrate that AutoPhOT can reproduce published lightcurves for a selection of transients with minimal human intervention.

We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 207 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $\kappa_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios are reproduced without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$).

Upcoming ground and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21-cm signal of the IGM. In the early universe, when the relative abundance of low-mass DM halos is important, measuring the global signal would place constraints on the damping of structure formation caused by DM having a higher relic velocity (warm dark matter, or WDM) than in cold dark matter (CDM). Such damping, however, can be mimicked by altering the star formation efficiency (SFE) and difficult to detect because of the presence of Pop III stars with unknown properties. We study these various cases and their degeneracies with the WDM mass parameter $m_X$ using a Fisher matrix analysis. We study the $m_X = 7$ keV case and a star-formation model that parametrizes the SFE as a strong function of halo mass and include several variations of this model along with three different input noise levels for the likelihood. We find that when the likelihood includes only Pop II stars, $m_X$ is constrained to $\sim 0.4$ keV for all models and noise levels at 68$\%$ CI. When the likelihood includes weak Pop III stars, $m_X \sim 0.3$ keV, and if Pop III star formation is relatively efficient, $m_X \sim 0.1$ keV, with tight Pop III star-formation parameter constraints. Our results show that the global 21-cm signal is a promising test-bed for WDM models, even in the presence of strong degeneracies with astrophysical parameters.

Spectral-timing techniques have proven valuable in studying the interplay between the X-ray corona and the accretion disc in variable active galactic nuclei (AGN). Under certain conditions, photo-ionized outflows emerging from central AGN regions also play a role in the observable spectral-timing properties of the nuclear components. The variable ionizing flux causes the intervening gas to ionize or recombine, resulting in a time-dependent absorption spectrum. Understanding the spectral-timing properties of these outflows is critical for the determination of their role in the AGN environment, but also the correct interpretation of timing signatures of other AGN components. In this paper, we test the capabilities of the Athena X-IFU instrument in studying the spectral and spectral-timing properties of a black hole system displaying a variable outflow. We take the narrow-line Seyfert 1 IRAS 13224-3809 as a test case. Our findings show that while the non-linear response of the absorbing medium can result in complex behaviour of time lags, the resulting decrease in the coherence can be used to constrain gas density and distance to the central source. Ultimately, modelling the coherence spectra of AGN outflows may constitute a valuable tool in studying the physical properties of the outflowing gas.

Ram pressure stripping is one of the most efficient mechanisms able to affect the gas reservoir in cluster galaxies and in the last decades many studies have characterized the properties of stripped galaxies. A definite census of the importance of this process in local clusters is though still missing. Here we characterize the fraction of galaxies showing signs of stripping at optical wavelengths, using the data of 66 clusters from the WINGS and OMEGAWINGS surveys. We focus on the infalling galaxy population and hence only consider blue, bright (B<18.2) late-type spectroscopically confirmed cluster members within 2 virial radii. In addition to "traditional" stripping candidates (SC) -- i.e. galaxies showing unilateral debris and tails -- we also consider unwinding galaxies (UG) as potentially stripped galaxies. Recent work has indeed unveiled a connection between unwinding features and ram pressure stripping and even though only integral field studies can inform on how often these features are indeed due to ram pressure, it is important to include them in the global census. We performed a visual inspection of B-band images and here we release a catalog of 143 UG. SC and UG each represent ~15-20% of the inspected sample. If we make the assumption that they both are undergoing ram pressure stripping, we can conclude that at any given time in the low-z universe about 35% of the infalling cluster population show signs of stripping in their morphology at optical wavelengths. These fractions depend on color, mass, morphology, and little on clustercentric distance. Making some rough assumptions on the duration of the tail visibility and on the time cluster galaxies can maintain blue colors, we infer that almost all bright blue late-type cluster galaxies undergo a stripping phase during their life, boosting the importance of ram pressure stripping in cluster galaxy evolution.

Some contracting or expanding stars are thought to host a large-scale magnetic field in their radiative interior. By interacting with the contraction-induced flows, such fields may significantly alter the rotational history of the star. They thus constitute a promising way to address the problem of angular momentum transport during the rapid phases of stellar evolution. In this work, we aim at studying the interplay between flows and magnetic fields in a contracting radiative zone. We propose a scenario that may account for the post-main sequence evolution of solar-like stars, in which a quasi-solid rotation can be maintained by a large-scale magnetic field during a contraction timescale. Then, an axisymmetric instability would destroy this large-scale structure and enables the differential rotation to set in. Such a contraction driven instability could also be at the origin of the observed dichotomy between strongly and weakly magnetic intermediate-mass stars.

With the sample of observed tidal disruption events (TDEs) now reaching several tens, distinct spectroscopic classes have emerged: TDEs with only hydrogen lines (TDE-H), only helium lines (TDE-He), or hydrogen in combination with He II and often N III/O III (TDE-H+He). Here we model the light curves of 32 optically-bright TDEs using the Modular Open Source Fitter for Transients (MOSFiT) to estimate physical and orbital properties, and look for statistical differences between the spectroscopic classes. For all types, we find a shallow distribution of star masses, compared to a typical initial mass function, between $\sim 0.1-1$ M$_\odot$, and no TDEs with very deep ($\beta \gg 1$) encounters. Our main result is that TDE-H events appear to come from less complete disruptions (and possibly lower SMBH masses) than TDE-H+He, with TDE-He events fully disrupted. We also find that TDE-H events have more extended photospheres, in agreement with recent literature, and argue that this could be a consequence of differences in the self-intersection radii of the debris streams. Finally, we identify an approximately linear correlation between black hole mass and radiative efficiency. We suggest that TDE-H may be powered by collision-induced outflows at relatively large radii, while TDE-H+He could result from prompt accretion disks, formed more efficiently in closer encounters around more massive SMBHs.

Tidal disruption events (TDEs) provide a means to probe the low end of the supermassive black hole (SMBH) mass distribution, as they are only observable below the Hills mass ($\lesssim 10^8$ M$_\odot$). Here we attempt to calibrate the scaling of SMBH mass with host galaxy bulge mass, enabling SMBH masses to be estimated for large TDE samples without the need for follow-up observations or extrapolations of relations based on high-mass samples. We derive host galaxy masses using Prospector fits to the UV-MIR spectral energy distributions for the hosts of 29 well-observed TDEs with BH mass estimates from MOSFiT. We then conduct detailed bulge/disk decomposition using SDSS and PanSTARRS imaging, and provide a catalog of bulge masses. We measure a positive correlation between SMBH and bulge mass for the TDE sample, with a power-law slope of 0.28 and significance $p=0.06$ (Spearmans) and $p=0.05$ (Pearsons), and an intrinsic scatter of 0.2 dex. Applying MC resampling and bootstrapping, we find a more conservative estimate of the slope is $0.18\pm0.11$, dominated by the systematic errors from Prospector and MOSFiT. This is shallower than the slope at high SMBH mass, which may be due to a bias in the TDE sample towards lower mass BHs that can more easily disrupt low-mass stars outside of the event horizon. When combining the TDE sample with that of the high mass regime, we find that TDEs are successful in extending the SMBH - stellar mass relationship further down the mass spectrum and provide a relationship across the full range of SMBH masses.

Studies of future space- and ground-based exoplanet surveys often rely on models of planetary systems to simulate instrument response, estimate scientific yields, perform trade analyses, and study efficient observation strategies. Until now, no planetary system models contained all of the basic physics necessary to enable study with all of the major exoplanet detection methods. Here we introduce a suite of such models generated by a new tool, exoVista. The exoVista tool quickly generates thousands of models of quasi-self-consistent planetary systems around known nearby stars at scattered light wavelengths and efficiently records the position, velocity, spectrum, and physical parameters of all bodies as functions of time. The modeled planetary systems can be used to simulate surveys using the direct imaging, transit, astrometric, and radial velocity techniques, as well as the overlap of these different methods.

We present late-time images of the site of the peculiar jet-driven TypeIIn supernova SN2010jp, including HST images taken 2-5 yr post explosion and deep ground-based images over a similar time. These are used to characterise its unusually remote environment and to constrain the progenitor's initial mass and age. The position of SN2010jp is found to reside along a chain of diffuse starlight that is probably an outer spiral arm or tidal tail of the interacting galaxy pair NGC2207/IC2163. There is one bright HII region projected within 1 kpc, and there is faint extended Halpha emission immediately surrounding the continuum source at the position of SN2010jp, which has $M_{F555W} = -7.7 (\pm 0.2)$ mag. In principle, the lingering light could arise from late-time circumstellar material (CSM) interaction, an evolved supergiant, a host star cluster, or some combination of these. Steady flux over 3 yr and a lack of strong, spatially unresolved Halpha emission make ongoing CSM interaction unlikely. If an evolved supergiant dominates, its observed luminosity implies an initial mass of roughly 22 $M_{\odot}$ and an age older than roughly 8 Myr. If the source is a star cluster, then its colour and absolute magnitude imply an age of 8-13 Myr and a modest cluster initial mass of log($M/M_{\odot}$) = 3.6-3.8. Extended Halpha emission out to a radius of 30 pc reveals a faint evolved HII region, pointing to recent star formation with at least one late O-type star. Based on these various clues, we conclude that the progenitor of SN2010jp had a likely initial mass of 18-22 $M_{\odot}$.

Within an imaging instrument's field of view, there may be many observational targets of interest. Similarly, within a spectrograph's bandpass, there may be many emission lines of interest. The brightness of these targets and lines can be orders of magnitude different, which poses a challenge to instrument and mission design. A single exposure can saturate the bright emission and/or have a low signal to noise ratio (SNR) for faint emission. Traditional high dynamic range (HDR) techniques solve this problem by either combining multiple sequential exposures of varied duration or splitting the light to different sensors. These methods, however, can result in the loss of science capability, reduced observational efficiency, or increased complexity and cost. The simultaneous HDR method described in this paper avoids these issues by utilizing a special type of detector whose rows can be read independently to define zones that are then composited, resulting in areas with short or long exposure measured simultaneously. We demonstrate this technique for the sun, which is bright on disk and faint off disk. We emulated these conditions in the lab to validate the method. We built an instrument simulator to demonstrate the method for a realistic solar imager and input. We then calculated SNRs, finding a value of 45 for a faint coronal mass ejection (CME) and 200 for a bright CME, both at 3.5 $R_{\odot}$ -- meeting or far exceeding the international standard for digital photography that defines a SNR of 10 as acceptable and 40 as excellent. Future missions should consider this type of hardware and technique in their trade studies for instrument design.

Observing the ultraviolet (UV) sky for time-variable phenomena is one of the many exciting science goals that can be achieved by a relatively small aperture telescope in space. The Near Ultraviolet Transient Surveyor (NUTS) is a wide-field ($3^\circ$) imager with a photon-counting detector in the near-UV (NUV, 200-300 nm), to be flown on an upcoming small satellite mission. It has a Ritchey-Chretien (RC) telescope design with correction optics to enable wide-field observations while minimizing optical aberrations. We have used an intensified CMOS detector with a solar blind photocathode, to be operated in photon-counting mode. The main science goal of the instrument is the observation of transient sources in the UV, including flare stars, supernovae, and active galactic nuclei. NUTS's aperture size and effective area enable observation of relatively unexplored, brighter parts of the UV sky which are usually not accessible to larger missions. We have designed, fabricated, and assembled the instrument, and the final calibrations and environmental tests are being carried out. In this paper, we provide the scientific motivation and technical overview of the instrument and describe the assembly and calibration steps.

This paper describes the Habitable Energy balance model for eXoplaneT ObseRvations (HEXTOR), which is a model for calculating latitudinal temperature profiles on Earth and other rapidly rotating planets. HEXTOR includes a lookup table method for calculating the outgoing infrared radiative flux and planetary albedo, which provides improvements over other approaches at parameterizing radiative transfer in an energy balance model. Validation cases are presented for present-day Earth and other Earth-sized planets with aquaplanet and land planet conditions from 0 to 45 degrees obliquity. A tidally locked coordinate system is also implemented in the energy balance model, which enables calculation of the horizontal temperature profile for planets in synchronous rotation around low mass stars. This coordinate transformed model is applied to cases for TRAPPIST-1e as defined by the TRAPPIST Habitable Atmosphere Intercomparison protocol, which demonstrates better agreement with general circulation models compared to the latitudinal energy balance model. Advances in applying energy balance models to exoplanets can be made by using general circulation models as a benchmark for tuning as well as by conducting intercomparisions between energy balance models with different physical parameterizations.

We present a near-infrared transmission spectrum of the long period (P=542 days), temperate ($T_{eq}$=294 K) giant planet HIP 41378 f obtained with the Wide-Field Camera 3 (WFC3) instrument aboard the Hubble Space Telescope (HST). With a measured mass of 12 $\pm$ 3 $M_{\oplus}$ and a radius of 9.2 $\pm$ 0.1 $R_{\oplus}$, HIP 41378 f has an extremely low bulk density (0.09 $\pm$ 0.02 g/cm$^{3}$). We measure the transit depth with a median precision of 84 ppm in 30 spectrophotometric channels with uniformly-sized widths of 0.018 microns. Within this level of precision, the spectrum shows no evidence of absorption from gaseous molecular features between 1.1-1.7 microns. Comparing the observed transmission spectrum to a suite of 1D radiative-convective-thermochemical-equilibrium forward models, we rule out clear, low-metallicity atmospheres and find that the data prefer high-metallicity atmospheres or models with an additional opacity source such as high-altitude hazes and/or circumplanetary rings. We explore the ringed scenario for this planet further by jointly fitting the K2 and HST light curves to constrain the properties of putative rings. We also assess the possibility of distinguishing between hazy, ringed, and high-metallicity scenarios at longer wavelengths with JWST. HIP 41378 f provides a rare opportunity to probe the atmospheric composition of a cool giant planet spanning the gap between the Solar System giants, directly imaged planets, and the highly-irradiated hot Jupiters traditionally studied via transit spectroscopy.

Transit spectroscopy is a powerful tool to decode the chemical composition of the atmospheres of extrasolar planets. In this paper we focus on unsupervised techniques for analyzing spectral data from transiting exoplanets. We demonstrate methods for i) cleaning and validating the data, ii) initial exploratory data analysis based on summary statistics (estimates of location and variability), iii) exploring and quantifying the existing correlations in the data, iv) pre-processing and linearly transforming the data to its principal components, v) dimensionality reduction and manifold learning, vi) clustering and anomaly detection, vii) visualization and interpretation of the data. To illustrate the proposed unsupervised methodology, we use a well-known public benchmark data set of synthetic transit spectra. We show that there is a high degree of correlation in the spectral data, which calls for appropriate low-dimensional representations. We explore a number of different techniques for such dimensionality reduction and identify several suitable options in terms of summary statistics, principal components, etc. We uncover interesting structures in the principal component basis, namely, well-defined branches corresponding to different chemical regimes of the underlying atmospheres. We demonstrate that those branches can be successfully recovered with a K-means clustering algorithm in fully unsupervised fashion. We advocate for a three-dimensional representation of the spectroscopic data in terms of the first three principal components, in order to reveal the existing structure in the data and quickly characterize the chemical class of a planet.

The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope.

After providing an overview of solar activity as measured by the sunspot number (SSN) and space weather events during solar cycles (SCs) 21-24, we focus on the weak solar activity in SC 24. The weak solar activity reduces the number of energetic eruptions from the Sun and hence the number of space weather events. The speeds of coronal mass ejections (CMEs), interplanetary (IP) shocks, and the background solar wind all declined in SC 24. One of the main heliospheric consequences of weak solar activity is the reduced total (magnetic + gas) pressure, magnetic field strength, and Alfv\'en speed. There are three groups of phenomena that decline to different degrees in SC 24 relative to the corresponding ones in SC 23: (i) those that decline more than SSN does, (ii) those that decline like SSN, and (iii) those that decline less than SSN does. The decrease in the number of severe space weather events such as high-energy solar energetic particle (SEP) events and intense geomagnetic storms is deeper than the decline in SSN. CMEs expand anomalously and hence their magnetic content is diluted resulting in weaker geomagnetic storms. The reduction in the number of intense geomagnetic storms caused by corotating interaction regions is also drastic. The diminished heliospheric magnetic field in SC 24 reduces the efficiency of particle acceleration, resulting in fewer high-energy SEP events. The numbers of IP type II radio bursts, IP socks, and high-intensity energetic storm particle events closely follow the number of fast and wide CMEs (and approximately SSN). The number of halo CMEs in SC 24 declines less than SSN does, mainly due to the weak heliospheric state. Phenomena such as IP CMEs and magnetic clouds related to frontside halos also do not decline significantly. The mild space weather is likely to continue in SC 25, whose strength has been predicted to be not too different from that of SC 24.

Colors and binarity provide important constraints on the Kuiper belt formation. The cold classical objects at radial distance r=42-47 au from the Sun are predominantly very red (spectral slope s>17%) and often exist as equal-size binaries (~30% observed binary fraction). This has been taken as evidence for the in-situ formation of cold classicals. Interestingly, a small fraction (~10%) of cold classicals is less red with s<17%, and these "blue" bodies are often found in wide binaries. Here we study the dynamical implantation of blue binaries from r<42 au. We find that they can be implanted into the cold classical belt from a wide range of initial radial distances, but the survival of the widest blue binaries -- 2001 QW322 and 2003 UN284 -- implies formation at r>30 au. This would be consistent with the hypothesized less-red to very-red transition at 30<r<40 au. For any reasonable choice of parameters (Neptune's migration history, initial disk profile, etc.), however, our model predicts a predominance of blue singles, rather than blue binaries, which contradicts existing observations. We suggest that wide blue binaries formed in situ at r=42-47 au and their color reflects early formation in a protoplanetary gas disk. The predominantly VR colors of cold classicals may be related to the production of methanol and other hydrocarbons during the late stages of the disk, when the temperature at 45 au dropped to 20 K and carbon monoxide was hydrogenated.

We study the optical $gri$ photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 200$ nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10-15 yr baselines and $\lesssim 100$ epochs. We find that the average damping timescale $\tau_{\rm DRW}$ continues to rise with increased baseline, reaching a median value of $\sim 750$ days ($g$ band) in the rest-frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained $\tau_{\rm DRW}$ (less than 20\% of the baseline), we confirm a weak wavelength dependence of $\tau_{\rm DRW}\propto \lambda^{0.51\pm0.20}$. We further quantify optical variability of these quasars over days to decades timescales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping timescales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on timescales less than $\sim$ a month for these luminous quasars, and this second break point correlates with the longer DRW damping timescale.

On the basis of GWTC-3, we discuss the detection prospect of extra-Galactic binary black holes (BBHs) by space gravitational-wave interferometers. In particular, targeting BBHs with component masses around 5-100$M_\odot$, we directly incorporate the chirp mass distribution of the 62 BBHs detected at high significance. We find that, due to the reduction of both the comoving merger rate and a weighted average of chirp masses, the expected detection numbers are generally much smaller than the results obtained by the same authors immediately after the report of GW150914. For LISA, the total BBH detections are estimated to be $N_{\rm tot}\sim 2 (T/4{\rm yr})^{3/2}(\rho_{\rm thr}/10)^{-3}$, dominated by nearly monochromatic BBHs ($\rho_{\rm thr}$: the detection threshold, $T$: the observational period). TianQin will have a total detection number $N_{\rm tot}$ similar to LISA. Meanwhile, TianQin has potential to find $N_{\rm mer}\sim0.6 (T/4{\rm yr})^{7/4}(\rho_{\rm thr}/10)^{-3}$ BBHs that merge in the observational period. This number for merging BBHs is 4-5 times larger than that of LISA, because of the difference between the optimal bands. We also investigate prospects for joint operations of multiple detectors, finding that concurrent observations will be more advantageous than sequential ones.

1SWASP J162545.15-043027.9 (WASP 1625-04) has been announced as one of EL CVn candidates showing total primary eclipses and ellipsoidal variations. This paper presents the absolute properties of the binary star, based on our high-resolution spectroscopy conducted from 2015 through 2020. From the spectral analysis, the radial velocities (RVs) for both components were obtained with the effective temperature $T_{\rm eff,1} = 8990 \pm 200$ K and the rotational rate $v_1\sin$$i=53\pm5$ km s$^{-1}$ for the more massive primary. The RV measurements were analyzed with archival WASP photometry. From the modelling we obtained: $M_1 = 1.745 \pm 0.013$ M$_\odot$, $M_2 = 0.187 \pm 0.002$ M$_\odot$, $R_1 = 1.626 \pm 0.008$ M$_\odot$, $R_2 = 0.290 \pm 0.003$ M$_\odot$, $L_1 = 15.5 \pm 1.4$ L$_\odot$, and $L_2 = 1.84 \pm 0.16$ L$_\odot$. In the Hertzsprung-Russell diagram, WASP 1625-04 A lies on the zero-age main sequence and its companion accords well with the helium-core white dwarf models of 0.19 M$_\odot$ in the constant luminosity phase. Our improved results demonstrate that WASP 1625-04 is a typical EL CVn-type binary with a low mass ratio and $M_2$ combination in the thin-disk population and is the product of the stable, non-conservative mass transfer of the precursor binary.

Recent observations revealed that several extremely low-density exoplanets show featureless transmission spectra. While atmospheric aerosols are a promising explanation for both the low density and featureless spectra, there is another attractive possibility: the presence of circumplanetary rings. Previous studies suggested that rings cause anomalously large transit radii. However, it remains poorly understood how rings affect the transmission spectrum. Here, we provide a framework to characterize the transmission spectra of ringed exoplanets. We develop an analytical prescription to include rings in the transmission spectra for arbitrarily viewing geometries. We also establish a simple post-processing model that can include the ring's effects on precomputed ring-free spectra. The ring flattens the transmission spectrum for a wide range of viewing geometries, consistent with the featureless spectra of extremely low-density exoplanets. Near-future observations by JWST at longer wavelengths would be able to distinguish the aerosol and ring scenarios. We also find that rocky rings might cause a silicate feature at $\sim$10 $\mu$m if the ring's optical depth is around unity. Thus, the ring's spectral features, if detected, would provide tight constrains on the physical properties of exoplanetary rings. We also discuss the ring's stability and suggest that thick rings are sustainable only at the equilibrium temperature of $\lesssim$300 K for the ring's age comparable to Kepler planets. This might indicate the intrinsic deficit of thick rings in the Kepler samples, unless rings are much younger than the planets as suggested for Saturn.

It is expected since the early 1970s that tenuous dust rings are formed by grains ejected from the Martian moons Phobos and Deimos by impacts of hypervelocity interplanetary projectiles. In this paper, we perform direct numerical integrations of a large number of dust particles originating from Phobos and Deimos. In the numerical simulations, the most relevant forces acting on dust are included: Martian gravity with spherical harmonics up to 5th degree and 5th order, gravitational perturbations from the Sun, Phobos, and Deimos, solar radiation pressure, as well as the Poynting-Robertson drag. In order to obtain the ring configuration, simulation results of various grain sizes ranging from submicron to 100 microns are averaged over a specified initial mass distribution of ejecta. We find that for the Phobos ring grains smaller than about 2 microns are dominant; while the Deimos ring is dominated by dust in the size range of about 5-20 microns. The asymmetries, number densities and geometrical optical depths of the rings are quantified from simulations. The results are compared with the upper limits of the optical depth inferred from Hubble observations. We compare to previous work and discuss the uncertainties of the models.

The technetium-rich (Tc-rich) M stars reported in the literature (Little-Marenin & Little (1979); Uttenthaler et al. (2013)) are puzzling objects since no isotope of technetium has a half-life longer than a few million years, and 99Tc, the longest-lived isotope along the s-process path, is expected to be detected only in thermally-pulsing stars enriched with other s-process elements (like zirconium). Carbon should also be enriched, since it is dredged up at the same time, after each thermal pulse on the asymptotic giant branch (AGB). However, these Tc-enriched objects are classified as M stars, meaning that they neither have any significant zirconium enhancement (otherwise they would be tagged as S-type stars) nor any large carbon overabundance (in which case they would be carbon stars). Here we present the first detailed chemical analysis of a Tc-rich M-type star, namely S Her. We first confirm the detection of the Tc lines, and then analyze its carbon and s-process abundances, and draw conclusions on its evolutionary status. Understanding these Tc-rich M stars is an important step to constrain the threshold luminosity for the first occurrence of the third dredge-up and the composition of s-process ejecta during the very first thermal pulses on the AGB.

In this paper, which is the first in a series of papers devoted to a detailed analysis of the death line of radio pulsars, we consider a possibility of producing secondary particles at a sufficiently long pulsar period P. To this end, we reconsidered the potential drop necessary for secondary plasma generation in the inner gap over magnetic polar caps. Our research made it possible to refine the conditions for generating secondary plasma, such as the multiplicity of the production of secondary particles and their spatial distribution. Our research also made it possible to further quantitatively analyse the dependence of the possibility of secondary plasma generation on all parameters, including the inclination angle of the magnetic axis to the rotation axis, the polar cap size and the magnetic field geometry.

Global warming is one of the problems of human civilization and decarbonization policy is the main solution to this problem. In this work, we propose using the gravity-assist by the asteroids to increase the orbital distance of the Earth from the Sun. We can manipulate the orbit of asteroids in the asteroid belt by solar sailing and propulsion engines to guide them towards the Mars orbit and a gravitational scattering can put asteroids in a favorable direction to provide an energy loss scattering from the Earth. The result would be increasing the orbital distance of the earth and consequently cooling down the Earth's temperature. We calculate the increase in the orbital distance of the earth for each scattering and investigate the feasibility of performing this project.

The surface brightness profiles (SBPs) of star clusters hold invaluable information on the dynamical state of clusters. The observed SBPs of star clusters, especially that of globular clusters, are in good agreement with the SBPs expected for isothermal spheres containing stars of reduced kinetic energies. However, the SBPs of configurations that satisfy these theoretical criteria cannot be uniquely expressed by analytical formulae, which had hindered the analysis of dynamical state of observed clusters in external galaxies. To counter this shortcoming, it has become a practice to use empirical fitting formulae that best represent the core and halo characteristics of theoretical models. We here present a general purpose code, named nProFit, that allows fitting of the surface brightness profiles of extragalactic star clusters to theoretical star clusters, defined by dynamical models of King (1966) and Wilson (1975). In addition, we also incorporated theoretical models that result in power-law surface brightness profiles represented by Elson et al. 1987. The code returns the basic size parameters such as core radius, half-light radius and tidal radius, as well as dynamically relevant parameters, such as the volume and surface density profiles, velocity dispersion profile, total mass and the binding energy for a user-fixed mass-to-light ratio. The usefulness of the code in the dynamical study of extragalactic clusters has been already illustrated in Cuevas-Otahola et al. 2020. The code, which is python-based at the user end, but makes calls to advanced routines in Pyraf and Fortran, is now available for public use. We provide example scripts and mock clusters in the installation package as guide to users.

Atomic species in the interstellar medium (ISM) transition out of their gas phase mainly by depletion onto dust. In this study, we examine if there is any change to the spectral line ratio predictions from a photoionization model of the Orion H II region when the degree of dust depletions is altered according to the most recently published model. We use equations and parameters published by previous works, in order to streamline the calculation of depleted abundances within CLOUDY. Our aim is for CLOUDY users to be able to vary the level of depletion using a single parameter in the input file. This makes it possible to explore predictions for a large range of depletions more efficiently. Finally, we discuss the results obtained for a model of the Orion Nebula when the degree of depletions are manipulated in this way. We found that the intensity of line ratios are significantly affected by depletions onto dust grains. Further, we found that adjusting dust abundances along with depletion affects the structure and the overall temperature of the H$^+$ layer across the H II region.

A family of unidentified infrared emission (UIE) bands has been observed throughout the Universe. The current observed spectral properties of the UIE bands are summarized. These properties are discussed in the frameworks of different models of the chemical carriers of these bands. The UIE carriers represent a large reservoir of carbon in the Universe, and play a significant role in the physical and chemical processes in the interstellar medium and galactic environment. A correct identification of the carrier of the UIE bands is needed to use these bands as probes of galactic evolution.

The effects of the heliospheric current sheet (HCS) on the evolution of Alfv\'enic turbulence in the solar wind are studied using MHD simulations incorporating the expanding-box-model (EBM). The simulations show that near the HCS, the Alfv\'enicity of the turbulence decreases as manifested by lower normalized cross helicity and larger excess of magnetic energy. The numerical results are supported by a superposed-epoch analysis using OMNI data, which shows that the normalized cross helicity decreases inside the plasma sheet surrounding HCS, and the excess of magnetic energy is significantly enhanced at the center of HCS. Our simulation results indicate that the decrease of Alfv\'enicity around the HCS is due to the weakening of radial magnetic field and the effects of the transverse gradient in the background magnetic field. The magnetic energy excess in the turbulence may be a result of the loss of Alfv\'enic correlation between velocity and magnetic field and the faster decay of transverse kinetic energy with respect to magnetic energy in a spherically expanding solar wind.

The Sun often produces coronal mass ejections with similar structure repeatedly from the same source region, and how these homologous eruptions are initiated remains an open question. Here, by using a new magnetohydrodynamic simulation, we show that homologous solar eruptions can be efficiently produced by recurring formation and disruption of coronal current sheet as driven by continuously shearing of the same polarity inversion line within a single bipolar configuration. These eruptions are initiated by the same mechanism, in which an internal current sheet forms slowly in a gradually sheared bipolar field and reconnection of the current sheet triggers and drives the eruption. Each of the eruptions does not release all the free energy but with a large amount left in the post-flare arcade below the erupting flux rope. Thus, a new current sheet can be more easily formed by further shearing of the post-flare arcade than by shearing a potential field arcade, and this is favorable for producing the next eruption. Furthermore, it is found that the new eruption is stronger since the newly formed current sheet has a larger current density and a lower height. In addition, our results also indicate the existence of a magnetic energy threshold for a given flux distribution, and eruption occurs once this threshold is approached.

Popular Fast Radio Burst models involve rotating magnetized neutron stars, yet no rotational periodicities have been found, even in observations of over 1500 bursts from each of FRB 121102 and FRB 20201124A. Periodograms of events with cosine-distributed random offsets as large as $\pm 0.6 P$ from a strict period $P$ reveal the underlying periodicity. Models of repeating FRB without intrinsic periodicity are considered, as are models of apparently non-repeating FRB.

In 2013, the method of determining mass ratios of dwarf novae using stage A superhumps (growing superhumps) was established. This method is a dynamical one in that it relies only on celestial mechanics. It is not dependent on an experimental calibration. Since then, more than 100 objects have been measured by this method. In this paper, I provide an updated description of the method. Comparisons with the results of the modern eclipse modeling method, which is considered to be the golden standard, have shown that these two methods agree very well and the stage A superhump method has been confirmed to be as accurate and as reliable as the modern eclipse modeling method. The number of the object by the former methods is now a few time those by the latter method, and the former is indispensable to study the terminal evolution of cataclysmic variables. I also showed that past studies by the other groups assumed incorrect fractional superhump excess relations, causing biases in discussing the evolution. I also derived a new experimental relation for stage B superhumps. The updated evolutionary track around the period minimum suggests that the angular momentum loss is 1.9 times larger than expected by gravitational wave radiation. The measurements of stage A superhumps greatly owe to international collaborations with amateurs and professionals. There is similarity with the world of ornithology in that both play a role in uniting the world via international exchanges of observations. I describe a summary of these collaborations and describe my thoughts about the relation between astronomy and ornithology, and give prospects how multidisciplinary works can be made possible between these seemingly distant fields of science (abridged).

The interest in the structure of ice-rich planetary bodies, in particular the differentiation between ice and rock, has grown due to the discovery of Kuiper belt objects and exoplanets. We thus carry out a parameter study for a range of planetary masses $M$, yielding radii $50 \aplt R \aplt 3000$~km, and for rock/ice mass ratios between 0.25 and 4, evolving them for 4.5~Gyr in a cold environment, to obtain the present structure. We use a thermal evolution model that allows for liquid and vapor flow in a porous medium, solving mass and energy conservation equations under hydrostatic equilibrium for a spherical body in orbit around a central star. The model includes the effect of pressure on porosity and on the melting temperature, heating by long-lived radioactive isotopes, and temperature-dependent serpentinization and dehydration. We obtain the boundary in parameter space [size, rock-content] between bodies that differentiate, forming a rocky core, and those which remain undifferentiated: small bodies, bodies with a low rock content, and the largest bodies considered, which develop high internal pressures and barely attain the melting temperature. The final differentiated structure comprises a rocky core, an ice-rich mantle, and a thin dense crust below the surface. We obtain and discuss the bulk density-radius relationship. The effect of a very cold environment is investigated and we find that at an ambient temperature of $\sim$20~K, small bodies preserve the ice in amorphous form to the present.

During the early phases of low-mass star formation, episodic accretion causes the ejection of high-velocity outflow bullets, which carry a fossil record of the driving protostar's accretion history. We present 44 SPH simulations of $1\,\mathrm{M}_{\odot}$ cores, covering a wide range of initial conditions, and follow the cores for five free-fall times. Individual protostars are represented by sink particles, and the sink particles launch episodic outflows using a subgrid model. The Optics algorithm is used to identify individual episodic bullets within the outflows. The parameters of the overall outflow and the individual bullets are then used to estimate the age and energetics of the outflow, and the accretion events that triggered it; and to evaluate how reliable these estimates are, if observational uncertainties and selection effects (like inclination) are neglected. Of the commonly used methods for estimating outflow ages, it appears that those based on the length and speed of advance of the lobe are the most reliable in the early phases of evolution, and those based on the width of the outflow cavity and the speed of advance are most reliable during the later phases. We describe a new method that is almost as accurate as these methods, and reliable throughout the evolution. In addition we show how the accretion history of the protostar can be accurately reconstructed from the dynamics of the bullets if each lobe contains at least two bullets. The outflows entrain about ten times more mass than originally ejected by the protostar.

Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions.

A highly misaligned gas disk around one component of a binary star system can undergo global Kozai-Lidov (KL) oscillations for which the disk inclination and eccentricity are exchanged. With hydrodynamical simulations of a gas and dust disk we explore the effects of these oscillations on the dust density distribution. For dust that is marginally coupled to the gas (${\rm St} \approx 1$), we find that the dust undergoes similar dynamical behaviour to the gas disk but the radial distribution of dust may be very different to the gas. The inward radial drift of the dust is faster in an eccentric disk leading to a smaller outer dust disk radius. The dust breaks into multiple narrow eccentric rings during the highly eccentric disk phase. Eccentric dust ring formation may have significant implications for the formation of planets in misaligned disks. We suggest that multiple dust rings may generally occur within gas disks that have sufficiently strong eccentricity peaks at intermediate radii.

Black hole X-ray binaries (BHXRBs) display a wide range of variability phenomena, from long duration spectral state changes to short-term broadband variability and quasi-periodic oscillations (QPOs). A particularly puzzling aspect is the production of QPOs, which -- if properly understood -- could be used as a powerful diagnostic tool of black hole accretion and evolution. In this work we analyse a high resolution three-dimensional general relativistic magnetohydrodynamic simulation of a geometrically thin accretion disk which is tilted by $65^{\circ}$ with respect to the black hole spin axis. We find that the Lense-Thirring torque from the rapidly spinning 10 $M_\odot$ black hole causes several sub-disks to tear off within $\sim 10-20$ gravitational radii. Tearing occurs in cycles on timescales of seconds. During each tearing cycle the inner sub-disk precesses for 1-5 periods before it falls into the black hole. We find a precession frequency of $\sim 3\rm Hz$, consistent with observed low-frequency QPOs. In addition, we find a high frequency QPO (HFQPO) with centroid frequency of $\sim55$Hz in the power spectra of the mass-weighted radius of the inner disk. This signal is caused by radial epicyclic oscillations of a dense ring of gas at the tearing radius, which strongly suggests a corresponding modulation of the X-ray lightcurve and may thus explain some of the observed HFQPOs.

MAGRATHEA is an open-source planet structure code that considers the case of fully differentiated spherically symmetric interiors. Given the mass of each layer, the code iterates the hydrostatic equations in order to shoot for the correct planet radius. The density may be discontinuous at a layer's boundary whose location is unknown. Therefore, in our case, shooting methods, which do not require predefined grid points, are preferred over relaxation methods in solving the two-point boundary value problem. The first version of MAGRATHEA supports a maximum of four layers of iron, silicates, water, and ideal gas. The user has many options for the phase diagram and equation of state in each layer and we document how to change/add additional equations of state. In this work, we present MAGRATHEA capabilities and discuss its applications. We encourage the community to participate in the development of MAGRATHEA at https://github.com/Huang-CL/Magrathea.

We present a new catalogue of distances and peculiar velocities (PVs) of $34,059$ early-type galaxies derived from Fundamental Plane (FP) measurements using data from the Sloan Digital Sky Survey (SDSS). This $7016\,\mathrm{deg}^{2}$ sample comprises the largest set of peculiar velocities produced to date and extends the reach of PV surveys up to a redshift limit of $z=0.1$. Alongside the data, we produce an ensemble of $2,048$ mock galaxy catalogues that reproduce the data selection function, and are used to validate our fitting pipelines and check for systematics. We uncover a significant trend between group richness and mean surface brightness within the sample, which hints at an environmental dependence within the FP and can result in biased peculiar velocities. This is removed using multiple FP fits as function of group richness, a procedure made tractable through a new analytic derivation for the integral of a 3D Gaussian over non-trivial limits. Our catalogue is calibrated to the zero-point of the CosmicFlows-III sample with an uncertainty of $0.004$ dex, which is cross-validated using the independent, predicted zero-point from the 2M++ reconstruction of our local velocity field. We achieve a mean uncertainty on the PVs of $\sim23\%$ compared to the CMB-frame redshift. Finally, as an example of what is possible with our new catalogue, we obtain preliminary bulk flow measurements up to a depth of $135\,h^{-1}\mathrm{Mpc}$. We find a slightly larger-than-expected bulk flow, although this could be caused by the presence of the Shapley supercluster lying just beyond the bounds of our data.

Deep convolutional neural networks (DCNNs) have become the most common solution for automatic image annotation due to their non-parametric nature, good performance, and their accessibility through libraries such as TensorFlow. Among other fields, DCNNs are also a common approach to the annotation of large astronomical image databases acquired by digital sky surveys. One of the main downsides of DCNNs is the complex non-intuitive rules that make DCNNs act as a ``black box", providing annotations in a manner that is unclear to the user. Therefore, the user is often not able to know what information is used by the DCNNs for the classification. Here we demonstrate that the training of a DCNN is sensitive to the context of the training data such as the location of the objects in the sky. We show that for basic classification of elliptical and spiral galaxies, the sky location of the galaxies used for training affects the behavior of the algorithm, and leads to a small but consistent and statistically significant bias. That bias exhibits itself in the form of cosmological-scale anisotropy in the distribution of basic galaxy morphology. Therefore, while DCNNs are powerful tools for annotating images of extended sources, the construction of training sets for galaxy morphology should take into consideration more aspects than the visual appearance of the object. In any case, catalogs created with deep neural networks that exhibit signs of cosmological anisotropy should be interpreted with the possibility of consistent bias.

Many disc galaxies in clusters have been found with bulges of similar age or younger than their surrounding discs, at odds with field galaxies of similar morphology and their expected inside-out formation. We use the EAGLE simulations to test potential origins for this difference in field and cluster galaxies. We find, in agreement with observations, that on average disc-dominated field galaxies in the simulations have older inner regions, while similar galaxies in groups and clusters have similarly aged or younger inner regions. This environmental difference is a result of outside-in quenching of the cluster galaxies. Prior to group/cluster infall, galaxies of a given present-day mass and morphology exhibit a similar evolution in their specific star formation rate (sSFR) profiles. Post-infall, the outer sSFRs of group and cluster galaxies significantly decrease due to interstellar medium stripping, while the central sSFR remains similar to field galaxies. Field disc galaxies instead generally retain radially increasing sSFR profiles. Thus, field galaxies continue to develop negative age gradients (younger discs), while cluster galaxies instead develop positive age gradients (younger bulges).

Chandra observations of the nearby, candidate Lyman-continuum (LyC) emitting galaxy Tol 0440-381 show brightening of an X-ray source by at least a factor of 4 to a luminosity of 1.6E40 erg/s over 3.8~days. The X-ray emission likely arises from either a low-luminosity AGN or an ultraluminous X-ray source. The properties of the X-ray source are similar to those found in Haro 11 and Tololo 1247-232, the only other LyC-emitting galaxies that have been resolved in X-rays. All three galaxies host luminous, variable, and hard spectrum X-ray sources that are likely accretion-powered. Accretion onto compact objects produces powerful outflows and ionizing radiation that could help enable LyC escape.

To better understand the physical processes associated with Jovian decametric (DAM) radio emissions, we present the statistical study of DAMs and inferred characteristics of DAM sources based on multi-view observation from Wind and STEREO spacecraft. The distribution of the apparent rotation speed of DAMs derived from multiple spacecraft suggests that the rotation speed of Io-related DAMs is in range of 0.15-0.6{\Omega}_J and that of non-Io-DAMs is between 0.7-1.2{\Omega}_J. Based on the method of Wang et al. (2020), we locate the sources of the DAMs and infer their emission angles and associated electron energies. The statistical results show that the DAM source locations have three preferred regions, two in the southern hemisphere and one in the northern hemisphere, which is probably caused by the non-symmetrical topology of Jupiter's magnetic field. The difference between Io-DAM source footprints and Io auroral spots changes with the Io's position in longitude, consistent with the previous results from Hess et al. (2010), Bonfond et al. (2017) and Hinton et al. (2019). In addition, the emission angles for non-Io-DAMs are smaller than that for Io-DAMs from the same source regions and all the emission angles range from 60{\deg} to 85{\deg}. Correspondingly, the electron energy is mainly distributed between 0.5 and 20 keV.

In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR due to the lack of direct magnetic field measurement in the corona. In this study, we investigate the MFR structures, originating from two active regions (ARs), AR 11719 and AR 12158, and estimate their magnetic properties quantitatively. We perform the nonlinear force-free field extrapolations with preprocessed photospheric vector magnetograms. In addition, remote-sensing observations are employed to find indirect evidence of MFRs on the Sun and to analyze the time evolution of magnetic reconnection flux associated with the flare ribbons during the eruption. A coherent "pre-existing" MFR structure prior to the flare eruption is identified quantitatively for one event from the combined analysis of the extrapolation and observation. Then the characteristics of MFRs for two events on the Sun before and during the eruption, forming the CME-MFR, including the axial magnetic flux, field-line twist, and reconnection flux, are estimated and compared with the corresponding in situ modeling results. We find that the magnetic reconnection associated with the accompanying flares for both events injects significant amount of flux into the erupted CME-MFRs.

We present high-resolution maps of the dust reddening in the Magellanic Clouds (MCs). The maps cover the Large and Small Magellanic Cloud (LMC and SMC) area and have a spatial angular resolution between $\sim$ 26 arcsec and 55 arcmin. Based on the data from the optical and near-infrared (IR) photometric surveys, including the Gaia Survey, the SkyMapper Southern Survey (SMSS), the Survey of the Magellanic Stellar History (SMASH), the Two Micron All Sky Survey (2MASS) and the near-infrared $YJK_{\rm{S}}$ VISTA survey of the Magellanic Clouds system (VMC), we have obtained multi-band photometric stellar samples containing over 6 million stars in the LMC and SMC area. Based on the measurements of the proper motions and parallaxes of the individual stars from Gaia Early Data Release 3 (Gaia EDR3), we have built clean samples that contain stars from the LMC, SMC and Milky Way (MW), respectively. We apply the spectral energy distribution (SED) fitting to the individual sample stars to estimate their reddening values. As a result, we have derived the best-fitting reddening values of ~ 1.9 million stars in the LMC, 1.5 million stars in the SMC and 0.6 million stars in the MW, which are used to construct dust reddening maps in the MCs. Our maps are consistent with those from the literature. The resultant high-resolution dust maps in the MCs are not only important tools for reddening correction of sources in the MCs, but also fundamental for the studies of the distribution and properties of dust in the two galaxies.

It has been suggested that there is possibly a class of stellar-mass black holes (BHs) residing near (distance $\le 10^3 M$) the galactic center massive black hole, Sgr A*. Possible formation scenarios include the mass segregation of massive stellar-mass black holes and/or the disk migration if there was an active accretion flow near Sgr A* within $\mathcal{O}(10)$ Myr. In this work, we explore the application of this type of objects as sources of space-borne gravitational wave detectors, such as Laser Interferometer Space Antenna (LISA). We find it is possible to probe the spin of Sgr A* based on the precession of the orbital planes of these stellar-mass black holes moving around Sgr A*. We also show that the dynamical friction produced by accumulated cold dark matter near Sgr A* generally produces small measurable phase shift in the gravitational waveform. In the case that there is an axion cloud near Sgr A*, the dynamical friction induced modification to gravitational waveform is measurable only if the mass of the axion field is in a narrow range of the mass spectrum. Gravitational interaction between the axion cloud and the stellar-mass black holes may introduce additional precession around the spin of Sgr A*. This additional precession rate is generally weaker than the spin-induced Lense-Thirring precession rate, but nevertheless may contaminate the spin measurement in a certain parameter regime. At last, we point out that the multi-body gravitational interaction between these stellar-mass black holes generally causes negligible phase shift during the LISA lifetime.

Astronomical observations in the X-ray band are subject to atmospheric attenuation and have to be performed in the space. CubeSats offer a cost effective means for space-based X-ray astrophysics but allow only limited mass and volume. In this article, we describe two successful CubeSat-based missions, HaloSat and PolarLight, both sensitive in the keV energy range. HaloSat was a 6U CubeSat equipped with silicon drift detectors. It conducted an all-sky survey of oxygen line emission and revealed the clumpy nature of the circumgalactic medium surrounding the Milky Way. PolarLight is a dedicated X-ray polarimeter performing photoelectron tracking using a gas pixel detector in a 1U payload. It observed the brightest X-ray objects and helped constrain their magnetic field or accretion geometry. On-orbit operation of both missions for multiple years demonstrates the capability of CubeSats as an effective astronomical platforms. The rapid time scales for development and construction of the missions makes them particularly attractive for student training.

We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxies, are still underdense compared to the mean density from the entire field. In addition, we classify galaxies into star-forming, quiescent, and green valley populations, and find that void galaxies tend to have slightly higher fractions of star-forming galaxies under the mass and redshift control, although the significance of this result is only moderate (2$\sigma$). However, when we focus on the star-forming population, the distribution of the specific star formation rate (sSFR) of void galaxies shows little difference from that of the control galaxies. Similarly, the median sSFR of star-forming void galaxies is also in good agreement with that of the star-forming control galaxies. Moreover, the effective green valley fraction of void galaxies, defined as the number of green valley galaxies over the number of nonquiescent galaxies, is comparable to that of the control ones, supporting the suggestion that void and control galaxies evolve under similar physical processes and quenching frequencies. Our results thus favor a scenario of the galaxy assembly bias.

Radio relics are extended radio emission features which trace shock waves in the periphery of galaxy clusters originating from cluster mergers. Some radio relics show a highly polarised emission, which make relics an excellent probe for the magnetisation of the intra-cluster medium. The origin of the relic polarisation is still debated. It could be a result of tangentially stretching the magnetic field at the shock surface. This scenario would naturally explain the alignment of the polarisation (E-vectors) with the shock normal. We have implemented a toy model for the relic polarisation according to this scenario. We find that the magnetic field strength itself crucially affects the fractional polarisation. Moreover, we find that the shock strength has surprisingly little effect on the overall polarisation fraction. Finally, we find that the fractional polarisation may decrease downstream depending on the magnetic field strength. Our results demonstrates that the shock compression scenario provides a very plausible explanation for the radio relic polarisation which specific features permitting to test the origin of radio relic polarisation.

We report here a new result extracted from the Fermi Large Area Telescope observation of the classical nova ASASSN-16ma that exhibits coherent gamma-ray pulsations at 544.84(7) seconds during its outburst in 2016. Considering the number of independent trials, the significance of the evidence is 4.0 sigma, equivalent to a false alarm probability of 5.9e-5. The periodicity was steady during the 4 days of its appearance, indicating its origin as the spinning signal of the white dwarf. Given that the optical and gamma-ray light curves of some shock-powered gamma-ray novae have been recently shown closely correlated to each other, the gamma-ray pulsation phenomenon likely implies an existence of the associated optical pulsations, which would provide detailed ephemerides for these extreme white dwarf binaries for further investigations in the near future.

Collecting and analysing X-ray photons over either spatial or temporal scales encompassing varying optical depth values requires knowledge about the optical depth distribution. In the case of sufficiently broad optical depth distribution, assuming a single column density value leads to a misleading interpretation of the source emission properties, nominally its spectral model. We present a model description for the interstellar medium absorption in X-ray spectra at moderate energy resolution, extracted over spatial or temporal regions encompassing a set of independent column densities. The absorption model (named disnht) approximates the distribution with a lognormal one and is presented in table format. The solution table and source code are made available and can be further generalized or tailored for arbitrary optical depth distributions encompassed by the extraction region. The disnht absorption model presented and its generalized solution are expected to be relevant for present and upcoming large angular scale analyses of diffuse X-ray emission, such as the ones from the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) and the future Athena missions.

Stellar physical and dynamical properties are essential knowledge to understanding the structure, formation, and evolution of our Galaxy. We produced an all-sky uniformly derived catalog of stellar astrophysical parameters (APs; age, mass, temperature, bolometric luminosity, distance, dust extinction) to give insight into the physical properties of Milky-Way stars. Exploiting the power of multi-wavelength and multi-survey observations from Gaia DR2 parallaxes and integrated photometry along with 2MASS and AllWISE photometry, we introduce an all-sky uniformly derived catalog of stellar astrophysical parameters, including dust extinction (A0) and average grain size (R0) along the line of sight, for 123,097,070 stars. In contrast with previous works, we do not use a Galactic model as prior in our analysis. We validate our results against other literature (e.g., benchmark stars, interferometry, Bayestar, StarHorse). The limited optical information in the Gaia photometric bands or the lack of ultraviolet or spectroscopic information renders the chemistry inference prior dominated. We demonstrate that Gaia parallaxes bring sufficient leverage to explore the detailed structures of the interstellar medium in our Milky Way. In Gaia DR3, we will obtain the dispersed optical light information to break through some limitations of this analysis, allowing us to infer stellar chemistry in particular. Gaia promises us data to construct the most detailed view of the chemo-dynamics of field star populations in our Galaxy. Our catalog is available from GAVO at this http URL (soon Gaia Archive and VizieR)

We report results of optical imaging and low-resolution spectroscopic monitoring of supernova (SN) 2017iro that occurred in the nearby ($\sim$\,31 Mpc) galaxy NGC 5480. The \ion{He}{1} 5876 \AA\, feature present in the earliest spectrum (--\,7 d) classified it as a Type Ib SN. The follow-up observations span from --\,7 to +\,266 d with respect to the $B$-band maximum. With a peak absolute magnitude in $V$-band, ($M_{V}$)\,=\,$-17.76\pm0.15$ mag and bolometric luminosity (log$_{10}$\,L)\,=\,42.39\,$\pm$\, 0.09 erg s$^{-1}$, SN 2017iro is a moderately luminous Type Ib SN. The overall light curve evolution of SN 2017iro is similar to SN 2012au and SN 2009jf during the early (up to $\sim$100 d) and late phases ($>$150 d), respectively. The line velocities of both \ion{Fe}{2} 5169 \AA\, and \ion{He}{1} 5876 \AA\, are $\sim$\,9000 km s$^{-1}$ near the peak. The analysis of the nebular phase spectrum ($\sim$\,+209 d) indicates an oxygen mass of $\sim$\,0.35 M$_{\odot}$. The smaller [\ion{O}{1}]/[\ion{Ca}{2}] flux ratio of $\sim$\,1 favours a progenitor with a zero-age main-sequence mass in the range $\sim$\,13--15 M$_{\odot}$, most likely in a binary system, similar to the case of iPTF13bvn. The explosion parameters are estimated by applying different analytical models to the quasi-bolometric light curve of SN 2017iro. $^{56}$Ni mass synthesized in the explosion has a range of $\sim$\,0.05\,--\,0.10 M$_{\odot}$, the ejecta mass $\sim$1.4\,--\,4.3 M$_{\odot}$ and the kinetic energy $\sim$\,0.8\,--\,1.9\,$\times$ 10$^{51}$ erg.

In this paper we emphasize the non-Lorentzian behavior of the Balmer series in helium-dominated DBA white dwarf stars for which the decades-old problem exists for the determination of the hydrogen abundance. In a very recent work, we have shown that quasi-molecular line satellites due to H-He and H-H collisions are responsible for the asymmetrical shape of the Lyman-alpha lines observed with the Cosmic Origin Spectrograph (COS) and that a similar asymmetry exists for the Balmer-alpha line profiles. In continuation with very recent work, where the n=2, 3 potential energies and transition dipole moments from the ground state were determined, here, we present accurate H-He potential energies and electronic transition dipole moments concerning the molecular states correlated with H(n=4)+He and their transition dipole moments with the states correlated with H(n=2)+He. Those new data are used to provide a theoretical investigation of the collisional effects in the blue wing of the Balmer-beta line of H perturbed by He. Because of the general trend characterizing the repulsive Sigma states of the potential energies involved in the Balmer series, the amplitude in the core of the line is decreasing very fast with the order of the series when the helium density gets as large as 10^21 cm^-3. This study is undertaken by applying a unified theory of spectral line broadening that is valid at very high helium densities found in DZA white dwarf stars. The treatment includes collision-induced (CI) line satellites due to asymptotically forbidden transitions, and it explains the asymmetry observed in their spectra.

Radial metallicity trends provide a key indicator of physical processes such as star formation and radial gas migration within a galaxy. Large IFU surveys allow for detailed studies of these radial variations, with recent observations detecting central dips in the metallicity, which may trace the impact of various evolutionary processes. However, the origin of these dips has not been conclusively determined, with suggestions that they may be diagnostic dependent. In this paper, we use the SDSS-IV MaNGA survey to investigate whether the observed dips represent genuine decreases in the central metallicity, or if they could be an artefact of the diagnostic used. Using a sub-sample of 758 local star-forming galaxies at low inclinations, we investigate in detail the impact of using different strong line diagnostics on the shapes of the returned profiles, and the prevalence of dips. We find no clear evidence of the dips being caused by changing values of the ionisation parameter within galaxies. To investigate physical causes, we explore both global and spatially-resolved parameters, finding that galaxies exhibiting central dips in the O3N2 metallicity profile have on average lower H$\alpha$EW values out to $R/R_\rm{e} \sim 1.5$, and higher values of D$_N$(4000) in the central regions. We additionally find a higher prevalence of dips in galaxies with high stellar mass, and lower values of global specific star formation rate, suggesting a possible link to central quenching. Nevertheless, these results are dependent on the diagnostic used, suggesting caution should be taken when interpreting observed features in galaxy metallicity gradients.

WASP-103b is the exoplanet with the highest expected deformation signature in its transit light curve and one of the shortest expected spiral-in times. Measuring the tidal deformation of the planet would allow us to estimate the second degree fluid Love number and gain insight into the planet's internal structure. Moreover, measuring the tidal decay timescale would allow us to estimate the stellar tidal quality factor, which is key to constraining stellar physics. We obtained 12 transit light curves of WASP-103b with the CHEOPS to estimate the tidal deformation and tidal decay of this extreme system. We modelled the high-precision CHEOPS transit light curves together with systematic instrumental noise using multi-dimensional Gaussian process regression informed by a set of instrumental parameters. To model the tidal deformation, we used a parametrisation model which allowed us to determine the second degree fluid Love number of the planet. We combined our light curves with previously observed transits of WASP-103b with HST and Spitzer. We estimate the radial Love number of WASP-103b to be $h_f = 1.59^{+0.45}_{-0.53}$. This is the first time that the tidal deformation is directly detected (at $3\, \sigma$) from the transit light curve of an exoplanet. Combining the transit times derived from CHEOPS, with the other transit times available in the literature, we find no significant orbital period variation for WASP-103b. However, the data show a hint of an orbital period increase instead of a decrease, as is expected for tidal decay. This could be either due to a visual companion star if this star is bound, the Applegate effect, or a statistical artefact. The estimated Love number of WASP-103b is similar to Jupiter's. This will allow us to constrain the internal structure and composition of WASP-103b, which could provide clues on the inflation of hot Jupiters. (Abridged)

We present the first search for signatures of brane-world extra-dimensional dark matter (DM) in the very-high-energy gamma-ray band by scrutinizing observations of the dwarf spheroidal galaxy Segue 1 with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope system. Branons are new degrees of freedom that appear within flexible brane-world models: they are weakly interacting massive particles and natural DM candidates. The ground-based gamma-ray telescopes MAGIC could indirectly detect branon DM in the multi-TeV mass range by observing secondary products of DM annihilation into Standard Model particles. In the absence of a signal, we place constraints on the branon DM parameter space by using a binned likelihood analysis of almost 160-hours deep exposure on the Segue 1 dwarf spheroidal galaxy by the MAGIC telescopes. Our most stringent limit to the thermally-averaged annihilation cross-section (at $95\%$ confidence level) corresponds to $ \langle \sigma v \rangle \simeq 1.4 \times 10^{-23}~\text{cm}^{3}\text{s}^{-1} $ at a branon mass of $ \sim 0.7~\text{TeV}$.

In the fuzzy dark matter (FDM) model, dark matter is composed of ultra-light particles with a de Broglie wavelength of $\sim$kpc, above which it behaves like cold dark matter (CDM). Due to this, FDM suppresses the growth of structure on small scales, which delays the onset of the cosmic dawn (CD) and the subsequent epoch of reionization (EoR). This leaves potential signatures in the sky averaged 21-cm signal (global), as well as in the 21-cm fluctuations, which can be sought for with ongoing and future 21-cm global and intensity mapping experiments. To do so reliably, it is crucial to include effects such as the dark-matter/baryon relative velocity and Lyman-Werner star-formation feedback, which also act as delaying mechanisms, as well as CMB and \lya heating effects, which can significantly change the amplitude and timing of the signal, depending on the strength of X-ray heating sourced by the remnants of the first stars. Here we model the 21-cm signal in FDM cosmologies across CD and EoR using a modified version of the public code 21cmvFAST that accounts for all these additional effects, and is directly interfaced with the Boltzmann code CLASS so that degeneracies between cosmological and astrophysical parameters can be fully explored. We examine the prospects to distinguish between the CDM and FDM models and forecast joint astrophysical, cosmological and FDM parameter constraints achievable with intensity mapping experiments such as HERA and global signal experiments like EDGES. We find that HERA will be able to detect FDM particle masses up to $m_{\rm FDM}\! \sim \!10^{-19}\,{\rm eV}\!-\!10^{-18}\,{\rm eV}$, depending on foreground assumptions, despite the mitigating effect of the delaying and heating mechanisms included in the analysis.

New ALMA observations of protoplanetary disks allow us to probe planet formation in other systems, giving us new constraints on planet formation processes. Meanwhile, studies of our own Solar System rely on constraints derived in a completely different way. However, it is still unclear what features the Solar System's disk could have produced during its gas phase. By running 2D isothermal hydro-simulations and a dust evolution model, we derive synthetic images at 1.3 mm wavelength using the radiative transfer code RADMC3D. We find that the embedded multiple giant planets strongly perturb the radial gas velocities of the disk, creating traffic jams in the dust. They produce over-densities different from the ones created by pressure traps and located away from the planets' positions in the disk. By deriving the images at 1.3mm from these dust distributions, we show that the traffic jams, observable with a high resolution, further blur the link between the number of gaps and rings in disks and the number of embedded planets. We additionally show that a system of 3 compact giant planets does not automatically produce bright outer rings at large radii in the disk. This means that high resolution observations of disks of various sizes are needed to distinguish between different giant planet formation scenarios during the disk phase, where the giants form either in the outer regions of the disks or in the inner regions. Finally, we find that, even when the dust temperature is determined self-consistently, the dust masses derived observationally might be off by up to a factor of ten compared to the dust contained in our simulations due to the creation of optically thick regions. Our study clearly shows that in addition to the constraints from exoplanets and the Solar System, ALMA has the power to constrain different stages of planet formation already during the first few million years.

Studying the cosmological sources at their cosmological rest-frames is crucial to track the cosmic history and properties of compact objects. In view of the increasing data volume of existing and upcoming telescopes/detectors, we here construct a 1--dimensional convolutional neural network (CNN) with a residual neural network (ResNet) structure to estimate the redshift of quasars in Sloan Digital Sky Survey IV (SDSS-IV) catalog from DR16 quasar-only (DR16Q) of eBOSS on a broad range of signal-to-noise ratios, named \code{FNet}. Owing to its $24$ convolutional layers and the ResNet structure with different kernel sizes of $500$, $200$ and $15$, FNet is able to discover the "\textit{local}" and "\textit{global}" patterns in the whole sample of spectra by a self-learning procedure. It reaches the accuracy of 97.0$\%$ for the velocity difference for redshift, $|\Delta\nu|< 6000~ \rm km/s$ and 98.0$\%$ for $|\Delta\nu|< 12000~ \rm km/s$. While \code{QuasarNET}, which is a standard CNN adopted in the SDSS routine and is constructed by 4 convolutional layers (no ResNet structure), with kernel sizes of $10$, to measure the redshift via identifying seven emission lines (\textit{local} patterns), fails in estimating redshift of $\sim 1.3\%$ of visually inspected quasars in DR16Q catalog, and it gives 97.8$\%$ for $|\Delta\nu|< 6000~ \rm km/s$ and 97.9$\%$ for $|\Delta\nu|< 12000~ \rm km/s$. Hence, FNet provides similar accuracy to \code{QuasarNET}, but it is applicable for a wider range of SDSS spectra, especially for those missing the clear emission lines exploited by \code{QuasarNET}. These properties of \code{FNet}, together with the fast predictive power of machine learning, allow \code{FNet} to be a more accurate alternative for the pipeline redshift estimator and can make it practical in the upcoming catalogs to reduce the number of spectra to visually inspect.

Classical T Tauri stars are surrounded by a circumstellar disk from which they are accreting material. This process is essential in the formation of Sun-like stars. Although often described with simple and static models, the accretion process is inherently time variable. Our aim is to examine the accretion process of the low-mass young stellar object CR Cha on a wide range of timescales from minutes to a decade by analyzing both photometric and spectroscopic observations from 2006, 2018, and 2019. We carried out period analysis on the light curves of CR Cha from the TESS mission and the ASAS-SN and the ASAS-3 databases. We studied the color variations of the system using $I,J,H,K$-band photometry obtained contemporaneously with the TESS observing window. We analyzed the amplitude, timescale, and the morphology of the accretion tracers found in a series of high-resolution spectra obtained in 2006 with the AAT/UCLES, in 2018 with the HARPS, and in 2019 with the ESPRESSO and the FEROS spectrographs. All photometric data reveal periodic variations compatible with a 2.327 days rotational period, which is stable in the system over decades. Moreover, the ASAS-SN and ASAS-3 data hint at a long-term brightening by 0.2 mag, between 2001 and 2008, and of slightly less than 0.1 mag in the 2015 - 2018 period. The near-infrared color variations can be explained by either changing accretion rate or changes in the inner disk structure. Our results show that the amplitude of the variations in the H$\alpha$ emission increases on timescales from hours to days/weeks, after which it stays similar even when looking at decadal timescales. On the other hand, we found significant morphological variations on yearly/decadal timescales, indicating that the different physical mechanisms responsible for the line profile changes, such as accretion or wind, are present to varying degrees at different times.

The polarization of photons emitted by astrophysical sources might be altered as they travel through a dark matter medium composed of ultra light axion-like particles (ALPs). In particular, the coherent oscillations of the ALP background in the galactic halo induce a periodic change on the polarization of the electromagnetic radiation emitted by local sources such as pulsars. Building up on previous works, we develop a new, more robust, analysis based on the generalised Lomb-Scargle periodogram to search for this periodic signal in the emission of the Crab supernova remnant observed by the QUIJOTE MFI instrument and 20 galactic pulsars from the Parkes Pulsar Timing Array (PPTA) project. We also carefully take into account the stochastic nature of the axion field, an effect often overlooked in previous works. This refined analysis leads to the strongest limits on the axion-photon coupling for a wide range of dark matter masses spanning $10^{-23}\text{ eV}\lesssim m_a\lesssim10^{-19} \text{ eV}$. Finally, we survey possible optimal targets and the potential sensitivity to axionic dark-matter in this mass range that could be achieved using pulsar polarimetry in the future.

There is evidence that the star formation process is linked to the intricate net of filaments in molecular clouds, which may be also due to gas compression from external triggers. We studied the southern region of the Perseus NGC 1333 molecular cloud, known to be heavily shaped by similar external triggers, to shed light on the process that perturbed the filament where the Class 0 IRAS4 protostars lie. We use new IRAM-NOEMA observations of SiO and CH3OH, both known to trace violent events as shocks, toward IRAS 4A as part of the Large Program Seeds Of Life in Space (SOLIS). We detected three parallel elongated ($>$6000 au) structures, called fingers, with narrow line profiles (~1.5 $km s^{-1}$) peaked at the cloud systemic velocity, tracing gas with high density (5-20 $10^5 cm^{-3}$) and high temperature (80-160 K). They are chemically different, with the northern finger traced by both SiO and CH3OH ([CH3OH]/[SiO]~160-300), while the other two only by SiO ([CH3OH]/[SiO]$<$ 40). Among various possibilities, a train of three shocks, distanced by $>$5000 yr, would be consistent with the observations if a substantial fraction of silicon, frozen onto the grain mantles, is released by the shocks.We suggest that the shock train is due to an expanding gas bubble, coming behind NGC 1333 from the southwest and clashing against the filament, where IRAS 4A lies. Finally, we propose a solution to the two-decades long debate on the nature and origin of the widespread narrow SiO emission observed in the south part of NGC 1333, namely that it is due to unresolved trains of shocks.

The stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA) is public and is widely used by the community. It includes the possibility of taking several non-standard processes such as atomic diffusion into account. Even if the effect of gravitational settling is considered a standard ingredient in stellar modelling today, this is not the case for radiative accelerations. The specific treatment of atomic diffusion along with the radiative accelerations has never been compared with other stellar evolution codes. Benchmarking these codes is important because improved accuracy is required in order to analyse data from present and future space missions, such as the \textit{Kepler}, Transiting Exoplanet Survey Satellite (TESS), and PLAnetary Transits and Oscillations of stars (PLATO) missions. The aim of this paper is to compare MESA models including atomic diffusion (with radiative accelerations) with models computed with the Montreal/Montpellier stellar evolution code and with the Code d'Evolution Stellaire Adaptatif et Modulaire (CESTAM). Additionally, we assess the impact of some MESA options related to atomic diffusion. We calculated atomic diffusion, including radiative accelerations, following the abundance profiles of 14 elements with MESA models. This was then compared with 1.1 and 1.4~$M_{\odot}$ models computed with the Montreal/Montpellier and CESTAM codes. Various tests of MESA options for atomic diffusion were also carried out by varying only one of them at a time. We find that the abundance profiles of the considered elements in the MESA models compare rather well with the models computed with the two other codes when atomic diffusion options are carefully set. We also show that some options in MESA are crucial for a proper treatment of atomic diffusion.

The chemistry of phosphorus in star- and planet-forming regions is poorly understood, despite the central role of phosphorus in terrestrial biochemistry. We present ALMA Band 3 and 4 observations of PO and PN towards the Class I protostar B1-a, representing the first spatially resolved observations of phosphorus carriers towards a Solar-type star forming region. The phosphorus molecules emit from two distinct clumps, which coincide with regions where the protostellar outflow (traced by SiO) interacts with a filament of dense gas (traced by CCS). Thus, the gas-phase phosphorus seems to originate from the shocking of dense interstellar clumps. Based on the observed emission patterns, PO and PN appear to be daughter products of a solid phosphorus carrier with an intermediate volatility between ices and silicate grains. Interstellar shocks may therefore play an important role in converting semi-refractory phosphorus to a more volatile form prior to incorporation into cometary ices. Indeed, the (PO+PN)/CH3OH ratio is similar in B1-a and comet 67P, implying a comparable reservoir of volatile phosphorus. The PO/PN ratio ranges from ~1-8 across B1-a. The northern emission clump exhibits a lower PO/PN ratio and weaker 13CH3OH emission than southern clump, indicating distinct shock physics and chemistry at the two positions. Resolved observations of P carriers towards additional sources are needed to better understand what regulates such variations in the PO/PN ratio in protostellar environments.

The morphology of haloes inform about both cosmological and galaxy formation models. We use the Minkowski Functionals (MFs) to characterize the actual morphology of haloes, only partially captured by smooth density profile, going beyond the spherical or ellipsoidal symmetry. We employ semi-analytical haloes with NFW and $\alpha\beta\gamma$-profile and spherical or ellipsoidal shape to obtain a clear interpretation of MFs as function of inner and outer slope, concentration and sphericity parameters. We use the same models to mimic the density profile of $N$-body haloes, showing that their MFs clearly differ as sensitive to internal substructures. This highlights the benefit of MFs at the halo scales as promising statistics to improve the spatial modeling of dark matter, crucial for future lensing, Sunyaev-Zel'dovich, and X-ray mass maps as well as dark matter detection based on high-accuracy data.

Binary driven hypernova (BdHN) models long gamma-ray burst (GRBs) as occurring in the binary systems involving a carbon-oxygen core (CO$_{\rm core}$) and a companion neutron star (NS) or a black hole (BH). This model, first proposed in 2012, succeeds and improves upon the fireshell model and the induced gravitational collapse (IGC) paradigm. After nearly a decade of development, the BdHN model has reached a nearly complete structure, explaining all the observables of long bursts into its theoretical framework, and has given a refined classification of long GRBs according to the original properties of the progenitors. In this article, we present a summary of the BdHN model and the physical processes at work in each of the envisaged Episodes during its occurrence and lifetime, duly contextualized in the framework of GRB observations.

With this Paper we complete a comprehensive study of substructure in dark matter haloes. In Paper I we derived the radial distribution and mass function (MF) of accreted subhaloes (scaled to the radius and mass of the host halo) and showed they are essentially universal. This is not the case, however, for those of stripped subhaloes, which depend on halo mass and assembly history. In Paper II we derived these latter properties in the simplest case of purely accreting haloes. Here we extend the study to ordinary haloes having suffered major mergers. After showing that all the properties of substructure are encoded in the mean truncated-to-original subhalo mass ratio profile, we demonstrate that the dependence of the subhalo MF on halo mass arises from their mass-dependent concentration, while the shape of the subhalo radial distribution depends on the time of the last major merger of the host halo. In this sense, the latter property is a better probe of halo formation time than the former. Unfortunately, this is not the case for the radial distribution of satellites as this profile is essentially disconnected from subhalo stripping and the properties of accreted subhaloes are independent of the halo formation time.

Multi-wavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nuclei (AGN) transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole. However, the formation, structure and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at 35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically-thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semi-opaque corona reach, respectively, Ti=10^10 K and Te=5x10^8 K. The system produces radiation, powerful collimated jets and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning black hole by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium can naturally lead to broadened iron line emission observed in the hard state.

Atmospheric retrievals of exoplanet transmission spectra allow constraints on the composition and structure of the day-night terminator region. Such retrievals in the past have typically assumed one-dimensional temperature structures which were adequate to explain extant observations. However, the increasing data quality expected from exoplanet spectroscopy with JWST motivates considerations of multidimensional atmospheric retrievals. We present AURA-3D, a three-dimensional atmospheric retrieval framework for exoplanet transmission spectra. AURA-3D includes a forward model that enables rapid computation of transmission spectra in 3D geometry for a given atmospheric structure and can, therefore, be used for atmospheric retrievals as well as for computing spectra from General Circulation Models (GCMs). In order to efficiently explore the space of possible 3D temperature structures in retrievals, we develop a parametric 3D pressure-temperature profile which can accurately represent azimuthally-averaged temperature structures of a range of hot Jupiter GCMs. We apply our retrieval framework to simulated JWST observations of hot Jupiter transmission spectra, obtaining accurate estimates of the day-night temperature variation across the terminator as well as the abundances of chemical species. We demonstrate an example of a model hot Jupiter transmission spectrum for which a traditional 1D retrieval of JWST-quality data returns biased abundance estimates, whereas a retrieval including a day-night temperature gradient can accurately retrieve the true abundances. Our forward model also has the capability to include inhomogeneous chemistry as well as variable clouds/hazes. This new retrieval framework opens the field to detailed multidimensional atmospheric characterisation using transmission spectra of exoplanets in the JWST era.

Plasma streaming instabilities play an important role in magnetic field amplification and particle acceleration in relativistic shocks and their environments. However, in the far shock precursor region where accelerated particles constitute a highly relativistic and dilute beam, streaming instabilities typically become inefficient and operate at very small scales when compared to the gyroradii of the beam particles. We report on a plasma cavitation instability that is driven by dilute relativistic beams and can increase both the magnetic field strength and coherence scale by orders of magnitude to reach near-equipartition values with the beam energy density. This instability grows after the development of the Weibel instability and is associated with the asymmetric response of background leptons and ions to the beam current. The resulting net inductive electric field drives a strong energy asymmetry between positively and negatively charged beam species. Large-scale particle-in-cell simulations are used to verify analytical predictions for the growth and saturation level of the instability and indicate that it is robust over a wide range of conditions, including those associated with pair-loaded plasmas. These results can have important implications for the magnetization and structure of shocks in gamma-ray bursts, and more generally for magnetic field amplification and asymmetric scattering of relativistic charged particles in plasma astrophysical environments.

We extend various recent results regarding the derivation of effective cosmological Friedmann equations from the dynamics of group field theory (GFT). Restricting ourselves to a single GFT field mode (or fixed values of Peter-Weyl representation labels), we first consider dynamics given by a quadratic Hamiltonian, which takes the form of a squeezing operator, and then add a quartic interaction that can be seen as a toy model for interactions in full GFT. Our derivation of effective Friedmann equations does not require a mean-field approximation; we mostly follow a general approach in which these equations in fact hold for any state. The resulting cosmological equations exhibit corrections to classical Friedmann dynamics similar to those of loop quantum cosmology, leading to generic singularity resolution, but also involve further state-dependent terms. We then specify these equations to various types of coherent states, such as Fock coherent states or Perelomov-Gilmore states based on the su(1,1) structure of harmonic quantum cosmology. We compute relative uncertainties of volume and energy in these states, clarifying whether they can be interpreted as semiclassical. In the interacting case, both analytical and numerical approximations are used to obtain modified cosmological dynamics. Our results clarify how effective cosmological equations derived from GFT can provide reliable approximations to the full dynamics.

We continue our analysis of a quantum cosmology model describing a flat Friedmann--Lema\^itre--Robertson--Walker universe filled with a (free) massless scalar field and an arbitrary perfect fluid. For positive energy density in the scalar and fluid, each classical solution has a singularity and expands to infinite volume. When quantising we view the cosmological dynamics in relational terms, using one degree of freedom as a clock for the others. Three natural candidates for this clock are the volume, a time variable conjugate to the perfect fluid, and the scalar field. We have previously shown that requiring unitary evolution in the "fluid" time leads to a boundary condition at the singularity and generic singularity resolution, while in the volume time semiclassical states follow the classical singular trajectories. Here we analyse the third option of using the scalar field as a clock, finding further dramatic differences to the previous cases: the boundary condition arising from unitarity is now at infinity. Rather than singularity resolution, this theory features a quantum recollapse of the universe at large volume, as was shown in a similar context by Paw{\l}owski and Ashtekar. We illustrate the properties of the theory analytically and numerically, showing that the ways in which the different quantum theories do or do not depart from classical behaviour directly arise from demanding unitarity with respect to different clocks. We argue that using a Dirac quantisation would not resolve the issue. Our results further illustrate the problem of time in quantum gravity.

The neutron lifetime anomaly has been used to motivate the introduction of new physics with hidden-sector particles coupled to baryon number, and on which neutron stars provide powerful constraints. Although the neutron lifetime anomaly may eventually prove to be of mundane origin, we use it as motivation for a broader review of the ways that baryon number violation, be it real or apparent, and dark sectors can intertwine and how neutron star observables, both present and future, can constrain them.

In the Penrose process and the Blandford-Znajek mechanism, the rotational energy of a black hole (BH) is extracted via particle fission and magnetic tension, respectively. Recently, inspired by a fundamental trait in plasma astrophysics known as magnetic reconnection (MR), a new energy extraction mechanism based on the fast reconnection of the magnetic field lines inside the ergosphere has been proposed by Comisso and Asenjo. In this paper, we investigate energy extraction caused by MR in the ergosphere of a rapidly spinning BH with broken Lorentz symmetry by a background bumblebee vector field. The desired rotating BH solution differentiates from the standard Kerr BH via the Lorentz symmetry breaking (LSB) parameter $l$, which comes from non-minimal coupling between the bumblebee field with non-zero vacuum expectation value and gravity. We find that incorporating $l<0$ in the background is in the interest of the energy extraction via MR for the fast-spinning BH surrounded by the plasma with weak magnetization, below what is expected from the scenario by Comisso and Asenjo. Our analysis robustly indicates that the power of energy extraction and efficiency of the plasma energization process through fast MR is more efficient than the Comisso-Asenjo solution, provided that the LSB parameter is negative, $l<0$. Compared to the Blandford-Znajek mechanism arising from the underlying background, we also show MR is a more efficient energy extraction mechanism if $l<0$.

We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea level production rate of $^{37}$Ar in natural xenon is estimated to be 0.024~atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1~tonne/month, the average $^{37}$Ar activity after 10~tonnes are purified and transported underground is 0.058--0.090~$\mu$Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic $^{37}$Ar will appear as a noticeable background in the early science data, while decaying with a 35~day half-life. This newly-noticed production mechanism of $^{37}$Ar should be considered when planning for future liquid xenon-based experiments.

In this paper we discuss on recent attempts aimed at demonstrating that the second clock effect (SCE) does not take place in Weyl spaces, which is contrary to well-known results. These attempts include Weyl gauge theories of gravity, as well as the symmetric teleparallel theories (STTs). Our approach to this issue is based on the power of Weyl gauge symmetry (WGS) which is a manifest symmetry of the basic laws of Weyl geometry. Through proper consideration of WGS we shall show that the SCE, being an effect of purely geometric nature, does not depend on the chosen theory of gravity and matter. Quite the contrary, the SCE singles out those matter couplings which are phenomenologically compatible with the underlying geometric laws. Here we consider both, spacetimes based in Weyl geometry with arbitrary nonmetricity (generalized Weyl geometry), as well as, standard Weyl spaces where the nonmetricity is proportional to the product of a Weyl gauge vector by the metric. This issue is of special relevance for the fate of the STTs which are being intensively applied in the cosmological framework. As we shall show, if realize that WGS is a manifest symmetry of generalized Weyl spaces, neither the Weyl gauge theories nor the STTs are free of the second clock effect, unless Weyl integrable geometry (WIG) spaces are considered.

We present a supersymmetric model where energy scales of a discrete $R$-symmetry breaking ($Z_{6R}$) and cosmic inflation are commonly attributed to the confinement scale of a hidden $Sp(2)$ strong dynamics. Apart from these, SUSY-breaking scale, the Higgsino mass and the right-handed neutrino masses are all shown to stem from $Z_{6R}$ breaking scale inferred from CMB observables. We will show that the model is characterized by the SUSY-breaking soft mass $m_{\rm soft}\simeq100-1000{\rm TeV}$ and the reheating temperature $T_{\rm rh}\simeq10^{9}{\rm GeV}$. Then we discuss how these predictions of the model can be tested with the help of the spectrum of the gravitational wave induced by the short-lived cosmic string present during the reheating era.

An inflationary scenario is expected to be embedded into an ultraviolet (UV) complete theory such as string theory. Quasi-heavy fields are ubiquitous in UV complete theory. The effect of these heavy fields may appear as nontrivial kinetic terms in the low energy effective field theory, which provides a nontrivial geometry in field space. In this paper, we study the effect of the geometry of multi-form-field space on an inflationary scenario. In particular, we focus on the geometric destabilization mechanism which induces the phase transition from the conventional slow-roll inflation to a novel inflationary scenario. Anisotropic inflation is a typical example of the new phase. To conform to observations, we restrict us to isotropic configuration of form fields. We clarify the conditions for the onset of the destabilization and reveal the geometric structure of attractors after the destabilization. We classify the viable models from the observational point of view. We also investigate the features of the primordial fluctuations and find the similarity to hyperbolic inflation. By calculating the power spectrum, we make several phenomenological predictions which are useful to discriminate our models from others inflation models.

Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding energies, two-nucleon separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. The rms deviation from the 637 mass data is 1.518 MeV, providing one of the best microscopic descriptions for nuclear masses. The drip lines obtained from DRHBc calculations are compared with other calculations, including the spherical relativistic continuum Hartree-Bogoliubov (RCHB) and triaxial relativistic Hartree-Bogoliubov (TRHB) calculations with PC-PK1. The deformation and continuum effects on the limits of the nuclear landscape are discussed. Possible peninsulas consisting of bound nuclei beyond the two-neutron drip line are predicted. The systematics of the two-nucleon separation energies, two-nucleon gaps, rms radii, quadrupole deformations, potential energy curves, neutron densities, neutron mean-field potentials, and pairing energies in the DRHBc calculations are also discussed. In addition, the $\alpha$ decay energies extracted are in good agreement with available data.

White dwarfs (WDs) and neutron stars (NSs) are among the most magnetized astrophysical objects in the universe, with magnetic fields that can reach up to $10^9\,\mathrm{G}$ for WDs and up to $10^{15}\,\mathrm{G}$ for NSs. The galaxy is expected to be populated with approximately one hundred million of double WD and millions of NS-WD binaries. Throughout the duration of the mission, the Laser Interferometer Space Antenna (LISA) will observe gravitational waves (GWs) emitted simultaneously by more than ten thousand of such galactic binaries. In this paper, we investigate the effect of the magnetic dipole-dipole interaction on the GW signal emitted by magnetic galactic binaries. We derive the secular equations governing the orbital and rotational motion of these objects. Then, we integrate these equations both numerically and analytically. We conclude that the overall visible effect is an additional secular drift of the mean longitude. This drift is proportional to the product of the magnetic moments and is inversely proportional to the $7/2$ power of the semi-major axis. Finally, we show that, at zeroth-order in eccentricity, the magnetic dipole-dipole interaction shifts the main frequency of the gravitational strain measured by LISA.

We highlight here several solutions developed to make high-level Cherenkov data FAIR: Findable, Accessible, Interoperable and Reusable. The first three FAIR principles may be ensured by properly indexing the data and using community standards, protocols and services, for example provided by the International Virtual Observatory Alliance (IVOA). However, the reusability principle is particularly subtle as the question of trust is raised. Provenance information, that describes the data origin and all transformations performed, is essential to ensure this trust, and it should come with the proper granularity and level of details. We developed a prototype platform to make the first H.E.S.S. public test data findable and accessible through the Virtual Observatory (VO). The exposed high-level data follows the gamma-ray astronomy data format (GADF) proposed as a community standard to ensure wider interoperability. We also designed a provenance management system in connection with the development of pipelines and analysis tools for CTA (ctapipe and gammapy), in order to collect rich and detailed provenance information, as recommended by the FAIR reusability principle. The prototype platform thus implements the main functionalities of a science gateway, including data search and access, online processing, and traceability of the various actions performed by a user.

A plasma whose Coulomb-collision rate is very small may relax on a shorter time scale to non-Maxwellian quasi-equilibria, which, nevertheless, have a universal form, with dependence on initial conditions retained only via an infinite set of Casimir invariants enforcing phase-volume conservation. These are distributions derived by Lynden-Bell (1967) via a statistical-mechanical entropy-maximisation procedure, assuming perfect mixing of phase-space elements. To show that these equilibria are reached dynamically, one must derive an effective `collisionless collision integral' for which they are fixed points -- unique and inevitable provided the integral has an appropriate H-theorem. We describe how such collision integrals are derived and what assumptions are required for them to have a closed form, how to prove the H-theorems for them, and why, for a system carrying sufficiently large electric-fluctuation energy, collisionless relaxation should be fast. It is suggested that collisionless dynamics may favour maximising entropy locally in phase space before converging to global maximum-entropy states.

Entropic-force cosmology provides, in contrast with dark energy descriptions, a concrete physical understanding of the accelerated expansion of the universe. The acceleration appears to be a consequence of the entropy associated with the information storage in the universe. We study the effects of including a subdominant power-law term within a thermodynamically admissible entropic-force model. The temperature of the universe horizon is obtained by requiring that the Legendre structure of thermodynamics is preserved. The correction term is introduced to explain different periods of acceleration and deceleration in the late-time universe. We analyze the various types of behaviors, and we satisfactorily compare them with the observational red-shift dependencies of the Hubble parameter $H$ and of the luminosity distance data available from supernovae.

[Background] Abundance anomalies in some globular clusters, such as the enhancement of potassium and the depletion of magnesium, can be explained in terms of an earlier generation of stars polluting the presently observed ones. It was shown that the potential range of temperatures and densities of the polluting sites depends on the strength of a few number of critical reaction rates. The reaction has been identified as one of these important reactions. [Purpose] The key ingredient for evaluating the thermonuclear reaction rate is the strength of the resonances which, at low energy, are proportional to their proton width. Therefore the goal of this work is to determine the proton widths of unbound 31P states. [Method] States in 31P were studied at the Maier-Leibnitz-Laboratorium using the one-proton transfer reaction. Deuterons were detected with the Q3D magnetic spectrometer. Angular distribution and spectroscopic factors were extracted for 27 states, and proton widths and resonance strengths were calculated for the unbound states. [Results] Several unbound states have been observed for the first time in a one-proton transfer reaction. Above 20 MK, the reaction rate is now entirely estimated from the observed properties of states. The reaction rate uncertainty from all resonances other than the resonance has been reduced down to less than a factor of two above that temperature. The unknown spin and parity of the resonance dominates the uncertainty in the rate in the relevant temperature range. [Conclusion] The remaining source of uncertainty on the reaction rate comes from the unknown spin and parity of the resonance which can change the reaction rate by a factor of ten in the temperature range of interest.