Papers for Tuesday, May 31 2022

We present a citation pattern analysis between astronomical papers and 13 other disciplines, based on the arXiv database over the past decade ($2010 - 2020$). We analyze 12,600 astronomical papers citing over 14,531 unique publications outside astronomy. Two striking patterns are unraveled. First, general relativity recently became the most cited field by astronomers, a trend highly correlated with the discovery of gravitational waves. Secondly, the fast growth of referenced papers in computer science and statistics, the first with a notable 15-fold increase since 2015. Such findings confirm the critical role of interdisciplinary efforts involving astronomy, statistics, and computer science in recent astronomical research.

Gravitational wave observations of quasicircular compact binary mergers imply complicated posterior measurements of their parameters. Though Gaussian approximations to the pertinent likelihoods have decades of history in the field, the relative generality and practical utility of these approximations hasn't been appreciated, given focus on careful, comprehensive generic Bayesian parameter inference. Building on our previous work in three dimensions, we demonstrate by example that bounded multivariate normal likelihood approximations are accurate, provide useful insight into individual sources and populations, and enable powerful fast calculations which would otherwise be inaccessible for population and low-latency parameter inference. We provide Normal Approximate Likelihood (NAL) fits for each event published in the Gravitational-Wave Transient Catalogs at https://gitlab.com/xevra/nal-data, with public code releases in the near future.

Collapsar disks have been proposed to be rich factories of heavy elements, but the major question of whether their outflows are neutron-rich, and could therefore represent significant sites of the rapid neutron-capture (r-) process, or dominated by iron-group elements, remains unresolved. We present global models of collapsars that start from a stellar progenitor and self-consistently describe the evolution of the disk, its composition, and its outflows in response to the imploding stellar mantle using energy-dependent M1 neutrino transport and an alpha-viscosity to capture turbulent angular-momentum transport. We find that a neutron-rich, neutrino-dominated accretion flow (NDAF) is established only marginally--either for very short times or very low viscosities--because the disk tends to disintegrate into an advective disk (ADAF) soon after its formation, launching powerful outflows but preventing it from developing a hot and dense, and therefore neutron-rich core. Viscous outflows disrupt the star within ~100 s with explosion energies close to that of hypernovae. If viscosity is neglected, a stable NDAF with disk mass of about 1 Msun is formed but is unable to release neutron-rich ejecta, while it produces a relatively mild explosion powered by a neutrino-driven wind blown off its surface. With ejecta electron fractions close to 0.5, all models presumably produce large amounts of 56Ni. Our results suggest that neutron-rich disks, and correspondingly r-process viable outflows, do not occur as readily as in remnant disks of neutron-star mergers. A weak effective viscosity generated by magnetohydrodynamic turbulence would improve the prospects for obtaining neutron-rich ejecta.

A non-negligible fraction of binary neutron star mergers are expected to form long-lived neutron star remnants, dramatically altering the multi-messenger signatures of a merger. Here, we extend existing models for magnetar-driven kilonovae and explore the diversity of kilonovae and kilonova afterglows. Focusing on the role of the (uncertain) magnetic field strength, we study the resulting electromagnetic signatures as a function of the external dipolar and internal toroidal fields. These two parameters govern, respectively, the competition between magnetic-dipole spindown and gravitational-wave spindown (due to magnetic-field deformation) of the rapidly-rotating remnant. We find that even in the parameter space where gravitational-wave emission is dominant, a kilonova with a magnetar central engine will be significantly brighter than one without an engine, as this parameter space is where more of the spin-down luminosity is thermalised. In contrast, a system with minimal gravitational-wave emission will produce a kilonova that may be difficult to distinguish from ordinary kilonovae unless early-epoch observations are available. However, as the bulk of the energy in this parameter space goes into accelerating the ejecta, such a system will produce a brighter kilonova afterglow that will peak on shorter times. To effectively hide the presence of the magnetar from the kilonova and kilonova afterglow, the rotational energy inputted into the ejecta must be $\lesssim 10^{-3}-10^{-2} E_{\rm rot}$. We discuss the different diagnostics available to identify magnetar-driven kilonovae in serendipitous observations and draw parallels to other potential magnetar-driven explosions, such as superluminous supernovae and broad-line supernovae Ic.

Recent analyses of Gaia data have resulted in the identification of new stellar structures, including a new class of extended stellar filaments called stellar "strings", first proposed by Kounkel & Covey 2019. We explore the spatial, kinematic, and chemical composition of strings to demonstrate that these newfound structures are largely inconsistent with being physical objects whose members share a common origin. Examining the 3D spatial distribution of string members, we find that the spatial dispersion around the claimed string spine does not improve in the latest Gaia EDR3 data release -- despite tangible gains in the signal-to-noise of the parallax measurements -- counter to expectations of a bona fide structure. Using the radial velocity dispersion of the strings (averaging $\rm\sigma_{V_r} \approx 15\;km\;s^{-1}$) to estimate their virial masses, we find that all strings are gravitationally unbound. Given the finding that the strings are dispersing, the reported stellar ages of the strings are typically $100\times$ larger than their measured dispersal times. Finally, we validate prior work that stellar strings are more chemically homogeneous than their local field stars, but show it is possible to obtain the same signatures of chemical homogeneity by drawing random samples of stars from spatially, temporally, and kinematically unrelated open clusters. Our results show that while some strings may be composed of real sub-structures, there is no consistent evidence for larger string-like connections over the sample. These results underline the need for caution in over-interpreting the significance of these strings and their role in understanding the star formation history of the Milky Way.

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, allows galaxies to be fit with a range of IMFs. Three surprising new features appear: (1) most star-forming galaxies are best fit with a bottom-lighter IMF than the Milky Way; (2) most star-forming galaxies at fixed redshift are fit with a very similar IMF; and (3) the most massive star-forming galaxies at fixed redshift instead exhibit a less bottom-light IMF, similar to that measured in quiescent galaxies. Additionally, since stellar masses and star formation rates both depend on the IMF, these results slightly modify the resulting relationship, while yielding similar qualitative characteristics to previous studies.

Observational signatures of accretion disc winds have been found in a significant number of low-mass X-ray binaries at either X-ray or optical wavelengths. The 2015 outburst of the black hole transient V404 Cygni provided a unique opportunity to study both types of outflows in the same system. Using contemporaneous X-ray (Chandra Observatory) and optical (Gran Telescopio Canarias, GTC) spectroscopy, in addition to hard X-ray light-curves (INTEGRAL), we show that the kinetic properties of the wind, as derived from P-Cyg profiles detected in the optical at low hard X-ray fluxes and in a number of X-ray transitions during luminous flares, are remarkably similar. Furthermore, strictly simultaneous data taken at intermediate hard X-ray fluxes show consistent emission line properties between the optical and the X-ray emission lines, which most likely arise in the same accretion disc wind. We discuss several scenarios to explain the properties of the wind, favouring the presence of a dynamic, multi-phase outflow during the entire outburst of the system. This study, together with the growing number of wind detections with fairly similar characteristic velocities at different wavelengths, suggest that wind-type X-ray binary outflows might predominantly have a multi-phase nature.

Direct-collapse black holes (DCBHs) may be the seeds of the first quasars, over 200 of which have now been detected at $z > 6$ . The James Webb Space Telescope (JWST) could detect DCBHs in the near infrared (NIR) at $z \lesssim 20$ and probe the evolution of primordial quasars at their earliest stages, but only in narrow fields that may not capture many of them. Wide-field NIR surveys by Euclid and the Nancy Grace Roman Survey Telescope (RST) would enclose far greater numbers of DCBHs but only directly detect them at $z \lesssim 6 - 8$ because of their lower sensitivities. However, their large survey areas will cover thousands of galaxy clusters and massive galaxies that could gravitationally lense flux from DCBHs, boosting them above current Euclid and RST detection limits and revealing more of them than could otherwise be detected. Here, we estimate the minimum number density of strongly lensed DCBHs and supermassive primordial stars required for detection in surveys by Euclid, RST and JWST at $z \lesssim 20$. We find that for reasonable estimates of host halo numbers RST, Euclid, and JWST could potentially find hundreds of strongly-lensed DCBHs at $z = 7 - 20$. RST would detect the most objects at $z \lesssim 10$ and JWST would find the most at higher redshifts. Lensed supermassive primordial stars could potentially also be found, but in fewer numbers because of their short lifetimes.

The behavior of an adaptive optics (AO) system for ground-based high contrast imaging (HCI) dictates the achievable contrast of the instrument. In conditions where the coherence time of the atmosphere is short compared to the speed of the AO system, the servo-lag error can become the dominant error term of the AO system. While the AO system measures the wavefront error and subsequently applies a correction (typically taking a total of one or a few milliseconds), the atmospheric turbulence above the telescope has changed resulting in the servo-lag error. In addition to reducing the Strehl ratio, the servo-lag error causes a build-up of speckles along the direction of the dominant wind vector in the coronagraphic image, severely limiting the contrast at small angular separations. One strategy to mitigate this problem is to predict the evolution of the turbulence over the delay time. Our predictive wavefront control algorithm minimizes, in a mean square sense, the wavefront error over the delay and has been implemented on the Keck II AO bench. In this paper, we report on the latest results of our algorithm and discuss updates to the algorithm itself. We explore how to tune various filter parameters based on both daytime laboratory tests and on-sky tests. We show a reduction in the residual-mean-square wavefront error for the predictor compared to the leaky integrator (the standard controller for Keck) implemented on Keck for three separate nights. Finally, we present contrast improvements for daytime and on-sky tests for the first time. Using the L-band vortex coronagraph for Keck's NIRC2 instrument, we find a contrast gain of up to 2 at a separation of 3 lambda/D and up to 3 for larger separations (3-7 lambda/D).

To combine information from measurements of the redshift-space power spectrum from spectroscopic data with angular weak lensing, galaxy clustering and galaxy-galaxy lensing power spectra from photometric surveys (i.e. the $3 \times 2$ point statistics), we must account for the covariance between the two probes. Currently any covariance between the two types of measurements is neglected as existing photometric and spectroscopic surveys largely probe different cosmological volumes. This will cease to be the case as data arrives from Stage-IV surveys. In this paper we derive an analytic expression for the covariance between photometric 2D angular power spectra and the 3D redshift-space power spectrum for Gaussian fields under the plane-parallel approximation. We find that the two probes are covariant on large radial scales, but because the information content of these modes is extremely low due to sample variance, we forecast that it is safe to neglect this covariance when performing cosmological parameter inference.

The reionization of the intergalactic medium at redshifts $z\gtrsim 6$ is expected to have a lasting impact on galaxies residing in low-mass dark matter halos. Unable to accrete or retain gas photo-heated to temperatures $T \gtrsim 10^4$ K, the star formation histories of faint galaxies in the early Universe are expected to decline as they exhaust their gas supply, resulting in a turn-over in the galaxy luminosity function (LF) and a potential solution to the missing satellites problem in the local group. Unfortunately, there are several challenges to constraining `reionization feedback' empirically, most notably that galaxies in low-mass halos are intrinsically faint, and that there are other physical mechanisms capable of inducing a turn-over in the LF. In this work, we investigate a new signature of reionization feedback that is in principle distinct from other processes: as faint galaxies are quenched by reionization, their stellar populations passively age and grow redder while the brighter galaxies nearby continue to form stars at an increasing rate and so remain relatively blue. We find that this contrast, between quenched and un-quenched galaxies induces a scale and colour-dependent signature in the present-day near-infrared background comparable to the expected sensitivity of NASA's upcoming SPHEREx mission. Whereas models with pure mass suppression largely affect the signal at wavelengths $\lesssim 2 \mu\rm{m}$, $\sim 5$%-level differences in the background persist out to $\sim 5 \mu\rm{m}$ for reionization feedback models. Finally, the power spectra of intensity ratio maps exhibit larger variations, and may thus be a promising target for future analyses.

The active galactic nucleus (AGN) phenomena results from a supermassive black hole accreting its surrounding gaseous and dusty material. The infrared (IR) regime provides most of the information to characterize the dusty structures that bridge from the galaxy to the black hole, providing clues to the black hole growth and host galaxy evolution. Over the past several decades, with the commissioning of various ground, airborne and space IR observing facilities, our interpretations of the AGN circumnuclear structures have advanced significantly through improved understanding of how their dust emission changes as a function of wavelength and how the heating of the dusty structures responds to variations of the energy released from the central engine. In this review, we summarize the current observational knowledge of the AGN IR broad-band spectral energy distributions (SEDs) and the IR time variability behavior covering large ranges of AGN luminosity and redshift, and discuss some first-order insights into the obscuring structures and host galaxy IR properties that can be obtained by integrating the relevant observations into a coherent picture.

The variability of blazars in the X-ray and optical regions informs both the physics of their emitting region and places demands on the observer if a program requires that the object be bright or faint. The extensive simultaneous X-ray and optical observation by Swift provides the best insight into the variable nature of these objects. This program uses \textit{Swift} data for 19 X-ray-bright blazars, generally at $z > 0.1$, to determine their variability properties. The analysis is based on structure functions and provides insight into the nature of the variability and how it depends on time, luminosity, and redshift. We also consider strategies for observing blazars at or above average brightness, given a time delay between planning an observation and obtaining the data. This is critical to observations with orbiting X-ray telescopes, current or future. The variability in the soft X-ray band is typically three to eight times larger than at UV-optical wavelengths, at fixed time differences (i.e., 30 or 100 days). There is almost no difference in the amplitude of variation (X-ray or UV-optical) as a function of redshift (time delay of 30 days) and a modest positive correlation with luminosity. In the X-ray band, blazars that become brighter than normal typically remain bright for at least 2-3 months, although with significant flickering. One can avoid observing objects that are significantly below the average X-ray flux by scheduling the observation when the $F_X > 0.9 F_{X,avg}$, which requires monitoring observations near the time of the scheduling activity.

We present 370 candidate eclipsing binaries (EBs), identified from ~510,000 short cadence TESS light curves. Our statistical criteria identify 5,105 light curves with features consistent with eclipses (~1% of the initial sample). After visual confirmation of the light curves, we have a final sample of 2,288 EB candidates. Among these, we find 370 sources that were not included in the catalog recently published by Prsa et al. We publish our full sample of 370 new EB candidates, and statistical features used for their identification, reported per observation sector.

We report the AGILE observations of GRB 220101A, which took place at the beginning of 1st January 2022 and was recognized as one of the most energetic gamma-ray bursts (GRBs) ever detected since their discovery. The AGILE satellite acquired interesting data concerning the prompt phase of this burst, providing an overall temporal and spectral description of the event in a wide energy range, from tens of keV to tens of MeV. Dividing the prompt emission into three main intervals, we notice an interesting spectral evolution, featuring a notable hardening of the spectrum in the central part of the burst. The average fluxes encountered in the different time intervals are relatively moderate, with respect to those of other remarkable bursts, and the overall fluence exhibits a quite ordinary value among the GRBs detected by MCAL. However, GRB 220101A is the second farthest event detected by AGILE, and the burst with the highest isotropic equivalent energy of the whole MCAL GRB sample, releasing E_iso=2.54x10^54 erg and exhibiting an isotropic luminosity of L_iso=2.34x10^52 erg/s (both in the 400 keV - 10 MeV energy range). We also analyzed the first 10^6 s of the afterglow phase, using the publicly available Swift XRT data, carrying out a theoretical analysis of the afterglow, based on the forward shock model. We notice that GRB 220101A is with high probability surrounded with a wind-like density medium, and that the energy carried by the initial shock shall be a fraction of the total E_iso, presumably near 50%.

Type 2 active galactic nuclei (AGN) show signatures of accretion onto a supermassive black hole through strong, high-ionization, narrow emission lines extended on scales of 100s to 1000s of parsecs, but they lack the broad emission lines from close in to the black hole that characterize type 1 AGN. The lack of broad emission could indicate obscuration of the innermost nuclear regions, or could indicate that the black hole is no longer strongly accreting. Since high-energy X-rays can penetrate thick obscuring columns, they have the power to distinguish these two scenarios. We present high-energy NuSTAR observations of 9 Seyfert 2 AGN from the IRAS 12 micron survey, supplemented with low-energy X-ray observations from Chandra, XMM-Newton, and Swift. The galaxies were selected to have anomalously low observed 2-10 keV luminosities compared to their [O III] optical luminosities, a traditional diagnostic of heavily obscured AGN, reaching into the Compton-thick regime for the highest hydrogen column densities. Based on updated [O III] luminosities and intrinsic X-ray luminosities based on physical modeling of the hard X-ray spectra, we find that one galaxy was misclassified as type 2 (NGC 5005) and most of the remaining AGN are obscured, including three confirmed as Compton-thick (IC 3639, NGC 1386, and NGC 3982). One galaxy, NGC 3627, appears to be recently deactivated. We also find a new X-ray changing-look AGN in NGC 6890.

The planetary nebula (PN) IC 4997 is one of a few rapidly evolving objects with variable brightness and nebular emission around a hydrogen-deficient star. In this study, we have determined the physical conditions and chemical abundances of this object using the collisionally excited lines (CELs) and optical recombination lines (ORLs) measured from the medium-resolution spectra taken in July 2014 with the FIbre-fed \'Echelle Spectrograph on the Nordic Optical Telescope at La Palma Observatory. We derived electron densities of $\gtrsim 3 \times 10^4$ cm$^{-3}$ and electron temperatures of $\gtrsim 14,000$ K from CELs, whereas cooler temperatures of $\sim 11,000$ and $\sim 7,000$ K were obtained from helium and heavy element ORLs, respectively. The elemental abundances deduced from CELs point to a metal-poor progenitor with [O/H]$\lesssim-0.75$, whereas the ORL abundances are slightly above the solar metallicity, [O/H]$\approx 0.15$. Our abundance analysis indicates that the abundance discrepancy factors (ADFs$\equiv$ORLs/CELs) of this PN are relatively large: ADF(O$^{2+})\gtrsim 8$ and ADF(N$^{2+})\gtrsim 7$. Further research is needed to find out how the ADFs and variable emissions are formed in this object and whether they are associated with a binary companion or a very late thermal pulse.

Common envelope evolution (CEE) physics plays a fundamental role in the formation of binary systems, such as mergering stellar gravitational wave sources, pulsar binaries and type Ia supernovae. A precisely constrained CEE has become more important in the age of large surveys and gravitational wave detectors. We use an adiabatic mass loss model to explore how the total envelope energy changes as a function of the remnant mass. This provides a more self-consistent way to calculate the binding energy of the donor. For comparison, we also calculate the binding energy through integrating the total energy from the core to the surface. We apply our results to 142 hot subdwarf binaries. For shorter orbital period sdBs, the binding energy is highly consistent. For longer orbital period sdBs in our samples, the binding energy can differ by up to a factor of 2. The CE efficiency parameter $\beta_\mathrm{CE}$ becomes smaller than $\alpha_\mathrm{CE}$ for the final orbital period $\log_{10} P_{\mathrm{orb}}/\mathrm{d} > -0.5$. We also find the mass ratios $\log_{10} q$ and CE efficiency parameters $\log_{10} \alpha_{\mathrm{CE}}$ and $\log_{10} \beta_{\mathrm{CE}}$ linearly correlate in sdBs, similarly to De Marco et al. (2010) for post-AGB binaries.

Both the lack of observation of ultra-high energy (UHE) photons and the limitations of the state-of-the-art methodology being applied for their identification motivate studies on alternative approaches to the relevant simulations and the related observational strategies. One of such new approaches is proposed in this report and it concerns new observables allowing indirect identification of UHE photons through cosmic ray phenomena composed of many spatially correlated extensive air showers or primary cosmic rays observed at one time. The study is based on simulations of interactions of UHE photons with the magnetic field of the Sun using the PRESHOWER program with some essential modifications. One of the expected results of such interactions is a generation of cosmic ray ensembles (CRE) in the form of very thin and very elongated cascades of secondary photons of energies spanning the whole cosmic ray energy spectrum. Upon entering the Earth's atmosphere, these cascades or their parts may generate uniquely characteristic walls of spatially correlated extensive air showers, and the effect is expected also in cases when primary UHE photons are not directed towards the Earth. Particle distributions in these multi-primary UHE photon footprints are expected to have thicknesses of the order of meters and elongations reaching even hundreds of millions kilometers, making them potentially observable with a global, multi-experiment approach, including re-exploring of the historical data, with the expected event rate exceeding the capabilities of even very large cosmic ray observatories. The methods described in this report allow for simulating potentially observable quantities related to CRE induced by UHE photons: densities, energy spectra and geographical orientations of particles at the top of the Earth's atmosphere.

Refereeing is a crucial component of publishing astronomical research, but few professional astronomers receive formal training on how to effectively referee a manuscript. In this article, we lay out considerations and best practices for referees. This document is intended as a tool for early career researchers to develop a fair, effective, and efficient approach to refereeing.

Some of the small inner moons of Uranus have very closely-spaced orbits. Multiple numerical studies have found that the moons Cressida and Desdemona, within the Portia sub-group, are likely to collide in less than 100 Myr. The subsequent discovery of three new moons (Cupid, Perdita, and Mab) made the system even more crowded. In particular, it has been suggested that the Belinda group (Cupid, Belinda, and Perdita) will become unstable in as little as 10$^5$ years. Here we revisit the issue of the stability of the inner moons of Uranus using updated orbital elements and considering tidal dissipation. We find that the Belinda group can be stable on $10^8$-year timescales due to an orbital resonance between Belinda and Perdita. We find that tidal evolution cannot form the Belinda-Perdita resonance, but convergent migration could contribute to the long-term instability of the Portia group. We propose that Belinda captured Perdita into the resonance during the last episode of disruption and re-accretion among the inner moons, possibly hundreds of Myr ago.

Stellar flares are characterized by sudden enhancement of electromagnetic radiation from the atmospheres of stars. Compared to their solar counterparts, our knowledge on the coronal plasma dynamics of stellar flares and their connection to coronal mass ejections (CMEs) remains very limited. With time-resolved high-resolution spectroscopic observations from the \textit{Chandra} X-ray observatory, we detected noticeable coronal plasma flows during several stellar flares on a nearby dMe star EV Lac. In the observed spectra of O~{\sc{viii}} (3 MK), Fe~{\sc{xvii}} (6 MK), Mg~{\sc{xii}} (10 MK), and Si~{\sc{xiv}} (16 MK) lines, these flare-induced upflows/downflows appear as significant Doppler shifts of several tens to \speed{130}, and the upflow velocity generally increases with temperature. Variable line ratios of the Si~{\sc{xiii}} triplet reveal that these plasma flows in most flares are accompanied by an increase of the coronal plasma density and temperature. We interpret these results as X-ray evidences for chromospheric evaporation on EV Lac. In two successive flares, the plasma flow pattern and a sharp increase of the measured coronal density are highly suggestive of explosive evaporation. The transition from redshifts to blueshifts in such an explosive evaporation occurs at a temperature of at least 10 MK, much higher than that observed in solar flares ($\sim$1 MK). However, in one flare the cool and warm upflows appear to be accompanied by a decreasing plasma density, which might be explained by a stellar filament/prominence eruption coupled to this flare. These results provide important clues to understand the coronal plasma dynamics during flares on M dwarfs.

RS Oph has persistently displayed flickering at optical wavelengths when observed away from its repeating nova outbursts. During the 2006 eruption the flickering disappeared, and this repeated during the recent 2021 event. We have been monitoring RS Oph looking for the reappearance of flickering at B-band following the 2021 outburst. The flickering was still absent (sigma(B)<0.002 mag) on day +210 (counted from nova optical maximum), appeared at sigma(B)=0.008 mag on day +224, and raised to sigma(B)=0.029 mag on day +250. On following dates the amplitude remained large, although fluctuating. The recovery of B-band quiescence brightness by RS Oph begun around day +225 and was completed by day +260. The parallel patterns followed by the rise in system brightness and the reappearance of flickering confirm the central role played in RS Oph by the return to pre-outburst conditions of the accretion disk and the refilling by the RG wind of the immediate circumstellar space.

We conduct a wide-band X-ray spectral analysis in the energy range of 1.5-100 keV to study the time evolution of the M7.6 class flare of 2016 July 23, with the Miniature X-ray Solar Spectrometer (MinXSS) CubeSat and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) spacecraft. With the combination of MinXSS for soft X-rays and RHESSI for hard X-rays, a non-thermal component and three-temperature multi-thermal component -- "cool" ($T \approx$ 3 MK), "hot" ($T \approx$ 15 MK), and "super-hot" ($T \approx$ 30 MK) -- were measured simultaneously. In addition, we successfully obtained the spectral evolution of the multi-thermal and non-thermal components with a 10 s cadence, which corresponds to the Alfv\'en time scale in the solar corona. We find that the emission measures of the cool and hot thermal components are drastically increasing more than hundreds of times and the super-hot thermal component is gradually appearing after the peak of the non-thermal emission. We also study the microwave spectra obtained by the Nobeyama Radio Polarimeters (NoRP), and we find that there is continuous gyro-synchrotron emission from mildly relativistic non-thermal electrons. In addition, we conducted a differential emission measure (DEM) analysis by using Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) and determine that the DEM of cool plasma increases within the flaring loop. We find that the cool and hot plasma components are associated with chromospheric evaporation. The super-hot plasma component could be explained by the thermalization of the non-thermal electrons trapped in the flaring loop.

We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter $D \ge 100\,{\rm km}$ bodies) should have been created, but the observed number is only 20. Therefore we computed a standard set of 125 simulations of collisions between representative $D = 100\,{\rm km}$ asteroids and high-speed icy projectiles (comets). According to our results, the largest remnant mass $M_{\rm lr}$ is similar as in low-speed collisions, due to appropriate scaling with the effective strength $Q_{\rm eff}$, but the largest fragment mass $M_{\rm lf}$ exhibits systematic differences - it is typically smaller for craterings and bigger for super-catastrophic events. This trend does not, however, explain the non-existence of old families. The respective parametric relations can be used in other statistical (Monte-Carlo) models to better understand collisions between asteroidal and cometary populations.

We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$\mu$m emission is detected at 10$\sigma$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $\Delta V_{\rm obs}/2\sigma_{\rm tot}=0.84\pm0.23$, where $\Delta V_{\rm obs}$ is the observed maximum velocity difference and $\sigma_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/\sigma_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $\sigma_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$.

Understanding the physical structures of the accreated matter very close to the black hole in quasars and active galactic nucleus (AGNs) is an important milestone to constrain the activities occurring in their centers. In this paper, we numerically investigate the effects of the asymptotic velocities on the physical structures of the accretion disk around the Kerr and Einstein-Gauss-Bonnet (EGB) rapidly rotating black holes. The Bondi-Hoyle accretion is considered with a falling gas towards the black hole in upstream region of the computational domain. The shock cones are naturally produced in the downstream part of the flow around both black holes. It is found that the structure of the cones and the amount of the accreated matter depend on asymptotic velocity $V_{\infty}$ (Mach number) and the types of the gravities (Kerr or EGB). Increasing the Mach number of the inflowing matter in the supersonic region causes the shock opening angle and accretion rates getting smaller because of the rapidly falling gas towards the black hole. The EGB gravity leads to an increase in the shock opening angle of the shock cones while the mass accretion rates $\dot{M}$ are decreasing in EGB gravity with a Gauss-Bonnet (GB) coupling constant $\alpha$. It is also confirmed that accretion rates and drag forces are significantly altered in the EGB gravity. Our numerical simulation results could be used to identify the accreation mechanism and physical properties of the accretion disk and black hole in the observed $X-$ rays such as NGC $1313$ $X-1$ and $1313$ $X-2$ and MAXI $J1803-298$.

We update constraints on the amplitude of the primordial trispectrum, using the final Planck mission temperature and polarization data. In the squeezed limit, a cosmological local trispectrum would be observed as a spatial modulation of small-scale power on the CMB sky. We reconstruct this signal as a source of statistical anisotropy via quadratic estimator techniques. We systematically demonstrate how the estimated power spectrum of a reconstructed modulation field can be translated into a constraint on $\tau_\mathrm{NL}$ via likelihood methods, demonstrating the procedures effectiveness by inferring known $\tau_\mathrm{NL}$ signal(s) from simulations. Our baseline results constrain $\tau_\mathrm{NL} < 1700$ at the 95\% confidence level, providing the most stringent constraints to date.

In strongly compressible magnetohydrodynamic turbulence, obliqueness between the large-scale density gradient and magnetic field gives an electromotive force mediated by density variance (intensity of density fluctuation). This effect is named ``magnetoclinicity'', and is expected to play an important role in large-scale magnetic-field generation in astrophysical compressible turbulent flows. Analysis of large-scale instability due to the magnetoclinicity effect shows that the mean magnetic-field perturbation is destabilised at large scales in the vicinity of strong mean density gradient in the presence of density variance.

The recent ALMA maps together with observations of H$_2$O maser emission seem to suggest the presence of a counter-rotation in the obscuring torus of NGC1068. We propose to explain this phenomenon as due to the influence of a wind and the effects of orientation. In order to test this idea: 1. we use N-body simulation of a clumpy torus in the gravitational field of a supermassive black hole taking into account mutual forces between clouds; 2. we consider the influence of the asymmetrical wind on the clouds in the torus throat. An axis of such a wind is tilted with respect to the torus symmetry axis; 3. we apply ray-tracing algorithm to take into account the effects of obscuration and of the torus orientation relative to an observer. By choosing the AGN parameters corresponding to NGC1068 we obtain the model velocity maps that emulate the effect of an apparent counter-rotation and can explain the discovery made by ALMA.

Previous Ku-band (15 GHz) imaging with data obtained from the Very Long Baseline Array (VLBA) had shown two compact, sub-pc components at the location of a presumed kpc-scale radio core in the Seyfert galaxy NGC 7674. It was then presumed that these two unresolved and compact components were dual radio cores corresponding to two supermassive black holes (SMBHs) accreting surrounding gas and launching radio-bright relativistic jets. However, utilizing the original VLBA dataset used to claim the detection of a binary SMBH, in addition to later multi-epoch/multi-frequency datatsets obtained from both the VLBA and the European VLBI Network, we find no evidence to support the presence of a binary SMBH. We place stringent upper limits to the flux densities of any sub-pc-scale radio cores which are at least an order of magnitude lower than the original VLBI radio-core detections, directly challenging the original binary SMBH detection claim. With this in mind, we discuss the possible reasons for the non-detection of any VLBI radio cores in our imaging, the possibility of a binary SMBH still residing in NGC 7674, and the prospect of future observations shedding further light on the true nature of this active galactic nucleus.

RR Lyrae stars are useful chemical tracers thanks to the empirical relationship between their heavy-element abundance and the shape of their light curves. However, the consistent and accurate calibration of this relation across multiple photometric wavebands has been lacking. We have devised a new method for the metallicity estimation of fundamental-mode RR Lyrae stars in the Gaia optical $G$ and near-infrared VISTA $K_s$ wavebands by deep learning. First, an existing metallicity prediction method is applied to large photometric data sets, which are then used to train long short-term memory recurrent neural networks for the regression of the [Fe/H] to the light curves in other wavebands. This approach allows an unbiased transfer of our accurate, spectroscopically calibrated $I$-band formula to additional bands at the expense of minimal additional noise. We achieve a low mean absolute error of $0.1$ dex and high $R^2$ regression performance of $0.84$ and $0.93$ for the $K_s$ and $G$ bands, respectively, measured by cross-validation. The resulting predictive models are deployed on the Gaia DR2 and VVV inner-bulge RR Lyrae catalogs. We estimate mean metallicities of $-1.35$ dex for the inner bulge and $-1.7$ for the halo, which are significantly less than values obtained by earlier photometric prediction methods. Using our results, we establish a public catalog of photometric metallicities of over 60,000 Galactic RR Lyrae stars, and provide an all-sky map of the resulting RR Lyrae metallicity distribution. The software code used for training and deploying our recurrent neural networks is made publicly available in the open-source domain.

Much effort has been done in order to better understand the active galactic nuclei mechanisms behind the relativistic jets observed in radio-loud sources. These phenomena are commonly seen in luminous objects with intermediate/high redshift such as quasars, so that the analysis of the spectroscopic properties of these sources may be a way to clarify this issue. Measurements are presented and contextualized taking advantage of the set of correlations associated with the quasar Main Sequence (MS), a parameter space that allows to connect observed properties to the relative relevance of radiative and gravitational forces. In the redshift range we consider, the low-ionization HI Balmer line H\b{eta} is shifted into the near infrared. Here we present first results of a sample of 22 high-luminosity quasars with redshift between 1.4 and 3.8. Observations covering the H\b{eta} spectral region were collected with the IR spectrometer ISAAC at ESO-VLT. Additional data were collected from SDSS in order to cover the UV region of our sources. The comparison between the strong C IV {\lambda}1549 high-ionization line and H\b{eta} in terms of line widths and shifts with respect to the rest-frame leads to an evaluation of the role of radiative forces in driving an accretion disk wind. While for non-jetted quasars the wind properties have been extensively characterized as a function of luminosity and other physical parameters, the situation is by far less clear for jetted sources. The overarching issue is the effect of the relativistic jets on the wind, and on the structure of the emitting region in general. We present results from our analysis of the optical and UV line profiles aimed to identifying the wind contribution to the line emission.

Luminous Blue Variables (LBVs) are massive stars that show strong spectral and photometric variability. The question of what evolutionary stages they represent and what exactly drives their instability is still open, and thus is it important to understand whether LBVs without significant ongoing activity exist, and for how long such dormant LBVs may ``sleep''. In this article we investigate the long-term variability properties of the LBV candidate MN112, by combining its optical and infrared spectral data covering 12 years with photometric data covering nearly a century acquired by both modern time-domain sky surveys and historical photographic plates. We analyze the spectra, derive physical properties of the star by modelling its atmosphere and use a new distance estimate from Gaia data release 3 (DR3) to determine the position of MN112 both inside the Galaxy and in the Hertzsprung--Russell diagram. Distance estimation increased in almost two times in comparison with Gaia DR2. Due to that MN112 moved in upper part of the diagram and according to our modeling it lies on evolutionary track for star with initial mass $M_*=70 M_\odot$ near Humphreys--Davidson limit. Given the absence of any significant variability we conclude that the star is a dormant LBV that has been inactive for at least a century now.

Infrared (IR) multiple-angle incidence resolution spectrometry (IR-MAIRS) is a recently developed spectroscopic technique that combines oblique incidence transmission measurements and chemometrics (multivariate analysis) to obtain both pure in-plane (IP) and out-of-plane (OP) vibration spectra for a thin sample. IR-MAIRS is established for analyzing the molecular orientation of organic thin films at atmospheric pressure, but it should also be powerful for the structural characterization of vapor-deposited thin samples prepared in a vacuum. The application of IR-MAIRS to vapor-deposited amorphous water is particularly interesting in the fields of physical and interstellar chemistry, because it is a representative model material for interstellar icy dust grains. We recently developed an experimental setup for in situ IR-MAIRS under low-temperature, ultra-high-vacuum conditions, which thus facilitates measurements of interstellar ice analogues such as vapor-deposited amorphous water. This review considers the theoretical framework of IR-MAIRS and our recent experimental results for vapor-deposited amorphous water. We present spectroscopic signatures for the perpendicular orientation of dangling OH bonds for three-coordinated water molecules at the surface of amorphous water at 90 K. The absolute absorption cross-section of the three-coordinated dangling OH bonds is quantitatively measured. As IR-MAIRS can essentially be conducted using only a Fourier-transform IR spectrometer and an angle-controllable linear polarizer, it is a useful, low-cost, and simple spectroscopic technique for studying laboratory analogues of interstellar ices including vapor-deposited amorphous water.

The particle diffusion coefficients of three TeV pulsar halos observed so far are inferred to be significantly smaller than the typical value of the interstellar medium (ISM). The anisotropic diffusion model ascribes the slow diffusion to the cross-field diffusion assuming sub-Alfv{\'e}nic turbulence in the ISM around the pulsar if the viewing angle between the observer's line-of-sight (LOS) to the pulsar and the local mean field direction is small. In general, the TeV halo's morphology under this model highly depends on the viewing angle, and an elongated, asymmetric morphology is predicted if the LOS is not approximately aligned with the local mean field direction. While the specific requirement of a small viewing angle is supposedly established only for a small fraction of TeV halos, TeV halos with apparent asymmetric morphology has not been detected. In this paper we will study the expectation of TeV halos measured by the TeV-PeV gamma-ray detector LHAASO in the framework of anisotropic diffusion model, with a particular focus on the influence of the viewing angle on the detectability. We show that a TeV halo is more detectable with a smaller viewing angle and this selection effect may explain why the morphologies of all three detected TeV halos so far are consistent being spherical. We also demonstrate that LHAASO is capable of detecting asymmetric TeV halos after several-year operation with reasonable source parameters. This can serve as a critical test of the anisotropic diffusion model.

We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1'$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy's measured properties are consistent with other FRB hosts with secure associations.

We study the impact of the loss of axial symmetry around the optical axis on the polarimetric properties of a telescope with segmented primary mirror when each segment is present in a different aging stage. The different oxidation stage of each segment as they are substituted in time leads to non-negligible crosstalk terms. This effect is wavelength dependent and it is mainly determined by the properties of the reflecting material. For an aluminum coating, the worst polarimetric behavior due to oxidation is found for the blue part of the visible. Contrarily, dust -- as modeled in this work -- does not significantly change the polarimetric behavior of the optical system . Depending on the telescope, there might be segment substitution sequences that strongly attenuate this instrumental polarization.

Astronomical data of the several recent years, which present an evidence in favour of abundant antimatter population in our Galaxy, Milky Way, are analysed. The data include: registration of gamma-rays with energy 0.511 MeV, which surely originate from electron-positron annihilation at rest, very large flux of anti-helium nuclei, discovered at AMS, and 14 stars which produce excessive gamma-rays with energies of several hundred MeV which may be interpreted as indication that these stars consist of antimatter. Theoretical predictions of these phenomena, made much earlier ago are described

Fluctuations during a prolonged maximum have been observed in several nova eruptions, although it is not clear, and can not be deduced directly from observations, if the phenomenon is an actual physical reaction to some mechanism originating in the erupting white dwarf, if it is occurring in the expanding ejected shell or if it is a form of interaction with the red dwarf companion. A handful of erupting nova models are investigated in this work, in order to assess the possibility of this sort of feature being an actual part of the eruption itself. The results explain the conditions that may possibly cause the formation of such fluctuations, favoring low mass WDs.

Primordial black holes (PBHs) are a fascinating candidate for being the dark matter, albeit one which has been heavily constrained. This review presents an in depth look at those observational constraints, particularly at their nuances and uncertainties. Despite their varied origins, the standard PBH formation path is assumed to be collapse of perturbations after inflation, which should leave signals visible in the CMB at certain scales. Other constraints come from microlensing surveys, which severely limit PBHs as dark matter in the solar to satellite range, but there are diminishing results in regards to lower mass ranges. Gravitational waves signals and PBH evaporation from Hawking radiation also make for useful probes, but the former requires the next generation of experiments before making constraints beyond the solar mass range, and the later is severely limited above $10^{-16} \rm M_{\odot}$. Other dynamical and accretion constraints exist for PBH of large masses. Care also has to be given, as all these constraints can carry different implications coming from differences between monochromatic and extended mass distributions, and their degree of clustering. Beyond all these issues, a window still exists for primordial black holes to be all of the dark matter between $10^{-16}$ and $10^{-11} \rm M_{\odot}$

Aims. This study aims to explore the relationship between the physical properties of different galaxy subclasses and their environment based on the analysis of 31 631 VIMOS Public Extragalactic Redshift Survey (VIPERS) galaxies observed at 0.5 < z < 0.9. Methods. We use the results of an unsupervised clustering algorithm to distinguish 11 subclasses of VIPERS galaxies based on the multi-dimensional feature space defined by rest-frame UV to NIR colours presented in Siudek et al (2018a). We investigate the relationship between the properties of these subclasses of galaxies and their local environment, defined as the galaxy density contrast derived from the 5th nearest neighbour technique. Results. We confirm that the galaxy population-density relation is already in place at z ~ 0.9, with the blue galaxy fraction decreasing with density, compensated by an increase of the red fraction. On average red galaxies in the high-density environment are larger by 28% than the ones in low-density environments. In particular, we find one group of galaxies, subclass C3, whose increase of size with time can be explained mainly as the result of mergers; for other red subclasses, mergers would not seem to play a major role (subclass C2) or play a negligible role (subclass C1). The properties of the green galaxies (subclasses C4-6) depend on whether their stellar mass is above or below a transition mass. Low-mass green galaxies appear to have grown through secular processes, while in high-mass green galaxies mass assembly appears to be dominated by mergers. When it comes to blue galaxies, the trend of decreasing fraction with denser environments seen for the group as a whole (subclasses C7-11) is found to be driven mostly by one group of galaxies, subclass C10. These are compact low-mass galaxies with high sSFRs, that are preferentially found in low-density environments.

The leading contribution to the uncertainties of atmospheric neutrino flux calculations arise from the cosmic ray nucleon flux and the production cross sections of secondary particles in hadron-air interactions. The data-driven model (DDM) developed in this work parametrizes particle yields from fixed-target accelerator data. The propagation of errors from the accelerator data to the inclusive muon and neutrino flux predictions results in smaller uncertainties than in previous estimates, and the description of atmospheric flux data is good. The model is implemented as part of the MCEq package, and hence can be flexibly employed for theoretical flux error estimation at neutrino telescopes.

A dynamical encounter between a stellar binary and Sgr A* in the Galactic Centre (GC) can tidally separate the binary and eject one member with a velocity beyond the escape speed of the Milky Way. These hypervelocity stars (HVSs) can offer insight into the stellar population in the inner parsecs of the Milky Way. In a previous work, our simulations showed that the lack of main sequence HVS candidates in current data releases from the Gaia space mission with precise astrometry and measured radial velocities places a robust upper limit on the ejection rate of HVSs from the GC of $3\times10^{-2} \, \mathrm{yr^{-1}}$. We improve this constraint in this work by additionally considering the absence of post main sequence HVSs in Gaia Early Data Release 3 as well the existence of the HVS candidate S5-HVS1. This evidence offers degenerate joint constraints on the HVS ejection rate and the stellar initial mass function (IMF) in the GC. For a top-heavy GC IMF as suggested by recent works, our modelling motivates an HVS ejection rate of $\eta=0.7_{-0.5}^{+1.5} \times10^{-4} \, \mathrm{yr^{-1}}$. This preferred ejection rate can be as large as $10^{-2} \, \mathrm{yr^{-1}}$ for a very top-light IMF and as low as 10$^{-4.5} \, \mathrm{yr^{-1}}$ if the IMF is extremely top-heavy. Constraints will improve further with future Gaia data releases, regardless of how many HVS candidates are found therewithin.

I obtained time-series photometry of the compact binary candidate for the Fermi source 4FGL J0935.3+0901. Superposed on the 2.44 hr orbital modulation are day-to-day variations and frequent flaring as seen in several redback and black widow millisecond pulsars (MSPs). The short orbital period favors a black widow. While the modulation of $\leq 1$ mag is smaller than that of most black widows, it could indicate a low orbital inclination. Although a published optical spectrum shows strong emission lines, the light curve evinces pulsar heating of the companion star rather than accretion-disk emission of a transitional MSP. Emission lines and flaring occur in the same objects, probably powered by shocks between the relativistic pulsar wind and a wind driven off the companion star. I also recovered the period in photometry from the Zwicky Transient Facility. A phase-connected ephemeris derived from MDM Observatory and ZTF data spanning 4 years yields a period of 0.10153276(36) days and an epoch for the ascending node of the putative pulsar.

GW170817 is the first binary neutron star (NS) merger detected in gravitational waves (GWs) and photons, and so far remains the only GW event of its class with a definitive electromagnetic (EM) counterpart. Radio emission from the structured jet associated with GW170817 has faded below the sensitivity achievable via deep radio observations with the most sensitive radio arrays currently in operation. Hence, we now have the opportunity to probe the radio re-brightening that some models predict, should emerge at late times from the interaction of the dynamically-stripped merger ejecta with the interstellar medium. Here we present the latest results from our deep radio observations of the GW170817 field with the Karl G. Jansky Very Large Array (VLA), 4.5 years after the merger. Our new data at $3\,$GHz do not show any compelling evidence for emission in excess to the tail of the jet afterglow ($<3.3\,\mu$Jy), confirming our previous results. We thus set new constraints on the dynamical ejecta afterglow models. These constraints favor single-speed ejecta with energy $\lesssim 10^{50}\,$erg (for an ejecta speed of $\beta_0=0.5$), or steeper energy-speed distributions of the kilonova ejecta. Our results also suggest larger values of the cold, non-rotating maximum NS mass in equal mass scenarios. However, without a detection of the dynamical ejecta afterglow, obtaining precise constraints on the NS equation of state remains challenging.

The dissipation of intense crustal electric currents produces high Joule heating rates in cooling neutron stars. Here it is shown that Joule heating can counterbalance fast cooling, making it difficult to infer the presence of hyperons (which accelerate cooling) from measurements of the observed thermal luminosity $L_\gamma$. Models with and without hyperon cores match $L_{\gamma}$ of young magnetars (with poloidal-dipolar field $B_{\textrm{dip}} \gtrsim 10^{14}$ G at the polar surface and $L_{\gamma} \gtrsim 10^{34}$ erg s$^{-1}$ at $t \lesssim 10^5$ yr) as well as mature, moderately magnetized stars (with $B_{\textrm{dip}} \lesssim 10^{14}$ G and $10^{31} \ \textrm{erg s}^{-1} \lesssim L_{\gamma} \lesssim 10^{32}$ erg s$^{-1}$ at $t \gtrsim 10^5$ yr). In magnetars, the crustal temperature is almost independent of hyperon direct Urca cooling in the core, regardless of whether the latter is suppressed or not by hyperon superfluidity. The thermal luminosities of light magnetars without hyperons and heavy magnetars with hyperons have $L_{\gamma}$ in the same range and are almost indistinguishable. Likewise, $L_{\gamma}$ data of neutron stars with $B_{\textrm{dip}} \lesssim 10^{14}$ G but with strong internal fields are not suitable to extract information about the equation of state as long as hyperons are superfluid, with maximum amplitude of the energy gaps of the order $\approx 1$ MeV.

AliCPT is the first Chinese cosmic microwave background (CMB) experiment which will make the most precise measurements of the CMB polarization in the northern hemisphere. The key science goal for AliCPT is the detection of primordial gravitational waves (PGWs). It is well known that an epoch of cosmic inflation, in the very early universe, can produce PGWs, which leave an imprint on the CMB in form of odd parity $B$-mode polarization. In this work, we study the performance of the component separation and parameter estimation pipelines in context of constraining the value of the tensor-to-scalar ratio. Based on the simulated data for one observation season, we compare five different pipelines with different working principles. Three pipelines perform component separation at map or spectra level before estimating $r$ from the cleaned spectra, while the other two pipelines performs a global fit for both foreground parameters and $r$. We also test different methods to account for the effects of time stream filtering systematics. This work shows that our pipelines provide consistent and robust constraints on the tensor-to-scalar ratio and a consistent sensitivity $\sigma(r) \sim 0.02$. This showcases the potential of precise $B$-mode polarization measurement with AliCPT-1. AliCPT will provide a powerful opportunity to detect PGWs, which is complementary with various ground-based CMB experiments in the southern hemisphere.

We present a study of the influence of the Gaia-Sausage-Enceladus (GSE) on the density shape of the Galactic stellar halo using 11624 K giants from the LAMOST survey. Every star is assigned a probability of being a member of the GSE based on its spherical velocities and metallicity by a Gaussian Mixture Model. We divide the stellar halo into two parts by the obtained probabilities, of which one is composed of the GSE members and defined as the GSE-related halo, and the other one is referred to as the GSE-removed halo. Using a non-parametric method, the radial number density profiles of the two stellar halos can be well described by a single power law with a variable flattening $q$ ($r = \sqrt{R^2+[(Z/q(r))]^2}, \nu = {\nu_0}r^{-\alpha}$). The index $\alpha$ is $4.92\pm0.12$ for the GSE-related halo and $4.25\pm0.14$ for the GSE-removed halo. Both the two stellar halos are vertically flattened at smaller radii but become more spherical at larger radii. We find that the GSE-related halo is less vertically flattened than the GSE-removed halo, and the difference of $q$ between the two stellar halos ranges from 0.07 to 0.15. However, after the consideration of the bias, it is thought to be within 0.08 at most of the radii. Finally, we compare our results with two Milky Way analogues which experience a significant major merger in the TNG50 simulation. The study of the two analogues also shows that the major merger-related stellar halo has a smaller ellipticity than the major merger-removed stellar halo.

The spectral observations and analysis for the W80 Region are presented by using the data of Medium-Resolution Spectroscopic Survey of Nebulae (MRS-N) with the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). A total of 2982 high-quality nebular spectra have been obtained in the 20 square degree field of view (FoV) which covers the W80 complex, and the largest sample of spectral data have been established for the first time. The relative intensities, radial velocities (RVs), and Full Widths at Half Maximum (FWHMs) are measured with the high spectral resolution of LAMOST MRS, for H$\alpha$ $\lambda$ 6563 \AA, [\ion{N}{ii}] $\lambda$$\lambda$ 6548 \AA, 6584 \AA \ , and [\ion{S}{ii}] $\lambda$$\lambda$ 6716 \AA, 6731 \AA \ emission lines. In the field of view of whole W80 Region, the strongest line emissions are found to be consistent with the bright nebulae, NGC 7000, IC 5070, and LBN 391, and weak line emissions also truly exist in the Middle Region, where no bright nebulae are detected by the wide-band optical observations. The large-scale spectral observations to the W80 Region reveal the systematic spatial variations of RVs and FWHMs, and several unique structural features. A 'curved feature' to the east of the NGC 7000, and a 'jet feature' to the west of the LBN 391 are detected to be showing with larger radial velocities. A 'wider FWHM region' is identified in the eastern part of the NGC 7000. The variations of [\ion{S}{ii}] / H${\alpha}$ ratios display a gradient from southwest to northeast in the NGC 7000 region, and manifest a ring shape around the 'W80 bubble' ionized by an O-type star in the L935. Further spectral and multi-band observations are guaranteed to investigate in detail the structural features.

After the morphological classification for the 18,190 $^{12}$CO molecular clouds, we further investigate the properties of their internal molecular gas structures traced by the $^{13}$CO($J=$ 1$-$0) line emissions. Using three different methods to extract the $^{13}$CO gas structures within each $^{12}$CO cloud, we find that $\sim$ 15$\%$ of $^{12}$CO clouds (2,851) have $^{13}$CO gas structures and these $^{12}$CO clouds contribute about 93$\%$ of the total integrated fluxes of $^{12}$CO emission. In each of 2,851 $^{12}$CO clouds with $^{13}$CO gas structures, the $^{13}$CO emission area generally does not exceed 70$\%$ of the $^{12}$CO emission area, and the $^{13}$CO integrated flux does not exceed 20$\%$ of its $^{12}$CO integrated flux. We reveal the strong correlation between the velocity-integrated intensities of $^{12}$CO lines and that of $^{13}$CO lines in both $^{12}$CO and $^{13}$CO emission regions. This indicates the H$_{2}$ column densities of molecular clouds are crucial for the $^{13}$CO lines emission. After linking the $^{13}$CO structures detection rates of the 18,190 $^{12}$CO molecular clouds to their morphologies, i.e. nonfilaments and filaments, we find that the $^{13}$CO gas structures are primarily detected in the $^{12}$CO clouds with filamentary morphologies. Moreover, these filaments tend to harbor more than one $^{13}$CO structure. That demonstrates filaments not only have larger spatial scales, they also have more molecular gas structures traced by $^{13}$CO lines, i.e. the local gas density enhancements. Our results favor the turbulent compression scenario for filaments formation, in which dynamical compression of turbulent flows induces the local density enhancements. The nonfilaments tend to be in the low-pressure and quiescent turbulent environments of the diffuse ISM.

Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding $2{\times}10^{17}$ eV using about 5.5 years of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between $10^{17}$ and $10^{18}$ eV.

We study the capabilities of present and future radio very-long-baseline-interferometry arrays to distinguish black holes from horizonless spacetimes. We consider an example of a horizonless spacetime, obtained by overspinning a regular black hole. Its image is distinct from the image of a Kerr spacetime due to a second set of photon rings interior to the shadow. These photon rings cannot be directly resolved by present and even next-generation Event Horizon telescope arrays, but instead imprint themselves in horizon-scale images as excess central brightness relative to that of a black hole. We demonstrate that future arrays can detect such indirect imprints.

We present a study of the bulge-to-total ratio (B/T) of a Ks-band-selected sample of 88 close major-merger pairs of galaxies (H-KPAIR) based on 2-D decomposition of SDSS r-band images with \textsc{galfit}. We investigate the dependence of the interaction-induced specific star formation rate enhancement ($\rm sSFR_{enh}$) on the B/T ratio, and the effects of this dependence on the differences between star-forming galaxies (SFGs) in spiral+spiral (S+S) and spiral+elliptical (S+E) pairs. Of all 132 spiral galaxies in H-KPAIR, the 44 in S+E pairs show higher B/T than those in the 44 S+S pairs, with means of $\rm B/T = 0.35 \pm 0.05$ and $\rm B/T = 0.26 \pm 0.03$, respectively. There is a strong negative dependence of $\rm sSFR_{enh}$ on the B/T ratio and only paired SFGs with $\rm B/T<0.3$ show significant ($>5\sigma$) enhancement. Paired SFGs in S+S pairs show a similar trend, and many disky SFGs ($\rm B/T<0.1$) in S+S have strong sSFR enhancements ($\rm sSFR_{enh} > 0.7$~dex). For SFGs in S+E, the sSFR has no clear B/T dependence, nor any significant enhancement in any B/T bin. Disky SFGs in S+S show significant ($>4\sigma$) enhancement in the molecular gas content ($\rm M_{H_2}/M_{star}$), while SFGs in S+E have no such enhancement in any B/T bin. No significant enhancement on total gas content ($\rm M_{gas}/M_{star}$) is found in any B/T bin for paired galaxies. The star formation efficiency of either the total gas ($\rm SFE_{gas} = SFR/M_{gas}$) or the molecular gas ($\rm SFE_{H_2} = SFR/M_{H_2}$) does not depend on the B/T ratio. The only significant ($>4\sigma$) SFE enhancement found for paired SFGs is the $\rm SFE_{gas}$ for disky SFGs in S+S pairs.

We investigate ways to produce the bifurcation observed in the stellar stream of the Sagittarius dwarf galaxy (Sgr). Our method consists in running $N$-body simulations of Sgr falling into the Milky Way for the last 3~Gyr, with added test particles on disk orbits that span a wide range of initial positions, energies, and angular momenta. We find that particles that end up in the faint branch are predominantly high angular momentum particles that can all originate from a single plane within the progenitor, nearly perpendicular both to the orbital plane of the progenitor and to the Milky Way stellar disk. Their original configuration at the start of the simulation corresponds to spiral features already present 3~Gyr ago, which could be, e.g., the result of a disk-like component being tidally perturbed, or the tidal tails of a satellite being disrupted within Sgr. We then run a simulation including the self-gravity of this disky component. Despite the remaining ambiguity of its origin, this disk component of the Sgr dwarf with spiral over-densities provides a first step towards a working model to reproduce the observed faint branch of the bifurcated Sgr stream.

We propose a new coupled three-form dark energy model to relieve the Hubble tension in this paper. Firstly, by performing a dynamical analysis with the coupled three-form dark energy model, we obtain four fixed points, including a saddle point representing a radiation dominated Universe, a saddle point representing a matter dominated Universe, and two attractors representing two saturated de Sitter Universes. Secondly, by confronting the coupled three-form dark energy model and the $\Lambda$ cold dark matter model (the $\Lambda$CDM model) with cosmic microwave background (CMB), baryonic acoustic oscillations (BAO), Type Ia supernovae (SN Ia) observations, we obtain $H_0= 67.8_{-0.6}^{+0.7}$($1\sigma$ level) km/s/Mpc for the coupled three-form dark energy model and $H_0=67.6_{-0.5}^{+0.5}$($1\sigma$ level) km/s/Mpc for the $\Lambda$CDM model, the former is in strong tension with the latest local measured $H_0$ value at $4.3\sigma$ confidence level, while the latter is in strong tension with the latest local measured $H_0$ value at $5.1\sigma$ level.

We carry out a thorough numerical examination of field theory monodromy inflation at strong coupling. We perform an MCMC analysis using a Gaussian likelihood, fitting multiparameter models using CMB constraints on the spectral index and the tensor to scalar ratio. We show that models with uniquely positive Wilson coefficients are ruled out. If there are coefficients that can take on both signs, there can be a cancellation of terms that flattens the potentials and allows one to satisfy current data, and forecasts with strong constraints on the tensor to scalar ratio. Models of field theory monodromy are naturally enhanced to include a mechanism for canceling off radiative corrections to vacuum energy, via vacuum energy sequestering (VES). Although they include a much larger parameter space, we find that a similar numerical examination yields no significant change in the Bayesian evidence for VES enhanced models, with naturalness considerations making them more attractive from a theoretical perspective.

Context: In the centre of pre-stellar cores, the deuterium fractionation is enhanced due to the cold temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to probe the evolution and the kinematics in the earliest stages of star formation. Aims: We analyse emission maps of cyclopropenylidene, c-C$_3$H$_2$, to study the distribution of the deuteration throughout the prototypical pre-stellar core L1544. Methods: We use single-dish observations of c-C$_3$H$_2$, c-H$^{13}$CC$_2$H, c-C$_3$HD, and c-C$_3$D$_2$ towards the pre-stellar core L1544, performed at the IRAM 30m telescope. We derive the column density and deuterium fraction maps, and compare these observations with non-LTE radiative transfer simulations. Results: The highest deuterium fractions are found close to the dust peak at the centre of L1544, where the increased abundance of H$_2$D$^+$ ions drives the deuteration process. The peak values are N(c-C$_3$HD)/N(c-C$_3$H$_2)=0.17\pm0.01$, N(c-C$_3$D$_2$)/N(c-C$_3$H$_2)=0.025\pm0.003$ and N(c-C$_3$D$_2$)/N(c-C$_3$HD$)=0.16\pm0.03$, which is consistent with previous single point observations. The distributions of c-C$_3$HD and c-C$_3$D$_2$ indicate that the deuterated forms of c-C$_3$H$_2$ in fact trace the dust peak and not the c-C$_3$H$_2$ peak. Conclusions: The N(c-C$_3$D$_2$)/N(c-C$_3$HD) map confirms that the process of deuteration is more efficient towards the centre of the core and demonstrates that carbon-chain molecules are still present at high densities. This is likely caused by an increased abundance of He$^+$ ions destroying CO, which increases the amount of carbon atoms in the gas phase.

Characterising the atmospheres of hot Jupiters is important in understanding the formation and migration of these exotic planets. However, there are still many open questions about the chemical and physical properties of these atmospheres. Here, we confirm the detection of water vapour in thermal emission from the non-transiting hot Jupiter {\tau} Bootis Ab with the high resolution NIR CARMENES spectrograph. Combining over 17 h of observations (560 spectra) and using a Bayesian cross-correlation to log-likelihood approach, we measure a systemic velocity of $V_{sys} = -11.51^{+0.59}_{-0.60}$ km/s and a radial velocity semi-amplitude of $K_{P} = 106.21^{+1.76}_{-1.71}$ km/s for the planet, which results in an absolute mass of $M_{P} = 6.24^{+0.17}_{-0.18} M_{J}$ and an orbital inclination of $41.6^{+1.0}_{-0.9}$ degrees. Our retrieved $V_{sys}$ shows a significant shift (+5 km/s) from the literature value, which could be caused by an inaccurate time of periastron. Within the explored model grid, we measure a preference for solar water abundance (VMR = $10^{-3}$) and no evidence for additional minor species in the atmosphere. Given the extensive orbital coverage of the data, we searched for a phase dependency in the water signal but found no strong evidence of variation with orbital phase. This detection is at odds with recent observations from SPIRou/CFHT and their tight upper limit on water abundance. We recommend further observations of the atmosphere {\tau} Bootis Ab to try and resolve these discrepancies.

The elastic crust of a neutron star fractures repeatedly as it spins down electromagnetically. An idealised, macroscopic model of inhomogeneous crustal failure is presented based on a cellular automaton with nearest-neighbour tectonic interactions involving strain redistribution and thermal dissipation. Predictions are made of the size and waiting-time distributions of failure events, as well as the rate of failure as the star spins down. The last failure event typically occurs when the star spins down to approximately 1% of its birth frequency with implications for rotational glitch activity. Neutron stars are commonly suggested as sources of continuous gravitational waves. The output of the automaton is converted into predictions of the star's mass ellipticity and gravitational wave strain as functions of its age, with implications for future observations with instruments such as the Laser Interferometer Gravitational Wave Observatory (LIGO).

The third observing run of the LIGO/Virgo/KARGA collaboration reported a few neutron star - black hole (NSBH) merger events. While NSBH mergers have yet to receive extensive theoretical attention, they may have a promising electromagnetic signature in the form of short gamma - ray bursts. Here we show that NSBH dynamical mergers can naturally form from ultra - wide binaries in the field. Flyby gravitational interactions with other neighbors in the galaxy in these ultra - wide systems may result in high eccentricity that drives the binary into a merger. We show that this process can result in a merger rate at the order of $\sim 10$~Gpc$^{-3}$~yr$^{-1}$ ($\sim 5$~Gpc$^{-3}$~yr$^{-1}$) for elliptical (spiral) galaxies. This channel predicts higher merger rate with higher velocity dispersion of the host - galaxy, delay time distribution which shallower than uniform but steeper that $1/t$, higher merger rate for lower BH to NS mass ratio.

The distribution of the muon content of highly inclined Monte Carlo cosmic ray showers is affected by the influence of Earth's geomagnetic field. It is found that the shapes of the positive and negative muon distributions get affected/modified by the influence of the Earth's geomagnetic field. Such a correlation between the earth's geomagnetic activity and the cosmic ray (CR) air shower muons is found sensitive to the primary cosmic ray mass composition.

Coronal Mass Ejections (CMEs) influence the interplanetary environment over vast distances in the solar system by injecting huge clouds of fast solar plasma and energetic particles (SEPs). A number of fundamental questions remain about how SEPs are produced, but current understanding points to CME-driven shocks and compressions in the solar corona. At the same time, unprecedented remote and in situ (Parker Solar Probe, Solar Orbiter) solar observations are becoming available to constrain existing theories. Here we present a general method for recognition and tracking on solar images of objects such as CME shock waves and filaments. The calculation scheme is based on a multi-scale data representation concept a trous wavelet transform, and a set of image filtering techniques. We showcase its performance on a small set of CME-related phenomena observed with the SDO/AIA telescope. With the data represented hierarchically on different decomposition and intensity levels, our method allows to extract certain objects and their masks from the imaging observations, in order to track their evolution in time. The method presented here is general and applicable to detecting and tracking various solar and heliospheric phenomena in imaging observations. It holds potential to prepare large training data sets for deep learning. We have implemented this method into a freely available Python library.

We present ALMA band 6/7 (1.3 mm/0.87 mm) and VLA Ka band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star forming region. We characterize the continuum and associated molecular line emission towards the most luminous protostars, i.e., IRS1 and IRS3, on ~100 au (0. 2") scales. IRS1 is partly resolved in millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well resolved disk appearance in millimeter continuum and is further resolved into a close binary system separated by ~40 au at 9 mm. Both sources exhibit clear velocity gradients across their disk major axes in multiple spectral lines including C18O, H2CO, SO, SO2, and complex organic molecules like CH3OH, 13CH3OH and CH3OCHO. We use an analytic method to fit the Keplerian rotation of the disks, and give constraints on physical parameters with a MCMC routine. The IRS3 binary system is estimated to have a total mass of 1.4-1.5$M_\odot$. IRS1 has a central mass of 3-5$M_\odot$ based on both kinematic modeling and its spectral energy distribution, assuming that it is dominated by a single protostar. For both IRS1 and IRS3, the inferred ejection directions from different tracers, including radio jet, water maser, molecular outflow, and H2 emission, are not always consistent, and for IRS1, these can be misaligned by ~50$^{\circ}$. IRS3 is better explained by a single precessing jet. A similar mechanism may be present in IRS1 as well but an unresolved multiple system in IRS1 is also possible.

FUV spectra of Eta Car, recorded across two decades with HST/STIS, document multiple changes in resonant lines caused by dissipating extinction in our line of sight. The FUV flux has increased nearly ten-fold which has led to increased ionization of the multiple shells within the Homunculus and photo-destruction of molecular hydrogen. Comparison of observed resonant line profiles with CMFGEN model profiles allows separation of wind-wind collision and shell absorptions from the primary wind, P Cygni profiles.The dissipating occulter preferentially obscured the central binary and interacting winds relative to the very extended primary wind. We are now able to monitor changes in the colliding winds with orbital phase. High velocity transient absorptions occurred across the most recent periastron passage, indicating acceleration of the primary wind by the secondary wind which leads to a downstream, high velocity bowshock that is newly generated every orbital period. There is no evidence of changes in the properties of the binary winds.

Context: Population III (Pop III) stars may be fast rotating. An expected consequence of fast rotation is strong internal mixing that deeply affects their evolutionary tracks in the Hertzsprung-Russell diagram and hence their ionising power. aims: We investigate the impact on the ionising power of Pop III stars in an extreme case of internal mixing, the one leading to chemically homogeneous evolution (CHE). In that situation, during the main sequence phase, the star keeps the same chemical composition from its center to its surface. Homogeneous stars have larger effective temperatures and luminosities than stars evolving non-homogeneously and thus are stronger ionising sources. Methods: The stellar evolution models are based on $n=3$ polytropes with a time varying hydrogen mass fraction. The ionisation model employs the self-similar champagne flow solution from Shu et al. (2002), as well as numerical simulations for the stochastic treatment of star clusters over a grid of redshifts and halo masses. Results: We find that haloes containing chemically homogeneous stars have an escape fraction of ionising photons about twice that of haloes containing classical Pop III stars. By extrapolating the high-$z$ ionisation history powered by Pop III stars (at $z\gtrsim 15$) to the post-reionisation epoch, we derive the Thomson scattering optical depth $\tau$, which is compared with the value measured by $\textit{Planck}$. We find that $\tau$ is overproduced by $\sim1.5- 5\sigma$, when all Pop III stars evolve homogeneously. This indicates that CHE is unlikely to be realised in the majority of Pop III stars, although the present study cannot exclude that a fraction of them undergoes CHE. Conclusions: Fast rotation might have a significant impact on the ionising budget of Pop III stars, and thus on early cosmic reionisation.

We derived radial velocities and CaT metallicity of more than 150 red giants stars in six SMC star clusters and their surrounding fields, with the instrument GMOS on GEMINI-S. The mean cluster radial velocity and metallicity were obtained with mean errors of 2.2 km\,s$^{-1}$ and 0.03 dex, while the mean field metallicities have a mean error of 0.13 dex. We add this information to that available for another 51 clusters and 30 fields with CaT metallicities on the same scale. Using this expanded sample we analize the chemical properties of the SMC Main Body, defined as the inner 3.4 degrees in semimajor axis. We found a high probability that the metallicity distribution of the Main Body clusters is bimodal with a metal-rich and a metal-poor cluster group, having mean metallicities with a dispersion of $\mu = -0.80$, $\sigma = 0.06$ and $\mu = -1.15$, $\sigma = 0.10$ dex, respectively. On the other hand, Main Body field stars show a unimodal metallicity distribution peaking at $[Fe/H] \sim -1$ and dispersion of $0.3$. Neither metal-rich nor metal-poor clusters present a metallicity gradient. However the full Main Body cluster sample and field stars have a negative metallicity gradient consistent with each other, but the one corresponding to clusters has a large error due to the large metallicity dispersion present in the clusters studied in that region. Metal-rich clusters present a clear age-metallicity relation, while metal-poor clusters present no chemical enrichment throughout the life of the galaxy. We present observational evidence that the chemical enrichment is complex in the SMC Main Body. Two cluster groups with potential different origins could be coexisting in the Main Body. More data with precise and homogeneous metallicities and distances are needed and dynamical simulations are required to understand possible different origins for the two possible cluster groups.

We present accurate and deep multi-band ($g,r,i$) photometry of the Local Group dwarf irregular galaxy NGC 6822. The images were collected with wide field cameras at 2m/4m- (INT,CTIO,CFHT) and 8m-class telescopes (SUBARU) covering a 2 square degrees FoV across the center of the galaxy. We performed PSF photometry of $\approx$7,000 CCD images and the final catalog includes more than 1 million objects. We developed a new approach to identify candidate field and galaxy stars, and performed a new estimate of the galaxy center by using old stellar tracers finding that it differs by 1.15 (RA) and 1.53 (DEC) arcmin from previous estimates. We also found that young (Main Sequence, Red Supergiants), intermediate (Red Clump, Asymptotic Giant Branch [AGB]) and old (Red Giant Branch [RGB]) stars display different radial distributions. Old stellar population is spherically distributed and extends to radial distances larger than previously estimated ($\sim$1 degree). The young population shows a well defined bar and a disk-like distribution, as suggested by radio measurements, that is off-center compared with old population. We discuss pros and cons of the different diagnostics adopted to identify AGB stars and develop new ones based on optical-NIR-MIR color-color diagrams (CCDs) to characterize Oxygen and Carbon (C) rich stars. We found a mean population ratio between Carbon and M-type (C/M) stars of 0.67$\pm$0.08 (optical/NIR/MIR) and we used the observed C/M ratio with empirical C/M-metallicity relations to estimate a mean iron abundance of [Fe/H]$\sim$-1.25 ($\sigma$=0.04 dex) that agrees quite well with literature estimates.

Dust emission at 8 micron has been extensively calibrated as an indicator of current star formation rate for galaxies and ~kpc-size regions within galaxies. Yet, the exact link between the 8 micron emission and the young stellar populations in galaxies is still under question, as dust grains can be stochastically heated also by older field stars. In order to investigate this link, we have combined mid-infrared images from the Spitzer Space Telescope with a published star cluster candidates catalog for the Local Group galaxy M33. M33 is sufficiently close that the Spitzer's 8 micron images resolve individual regions of star formation. Star clusters represent almost-single-age stellar populations, which are significantly easier to model than more complex mixtures of stars. We find a decrease in the 8 micron luminosity per unit stellar mass as a function of age of the star clusters, with a large scatter that is consistent with varying fractions of stellar light absorbed by dust. The decrease and scatter both confirm findings based on more distant galaxies and are well described by simple models for the dust emission of a young stellar population. We conclude that the dust emission at 8 micron depends sensitively on the age of the stellar population, out to at least the oldest age analyzed here, ~400 Myr. This dependence complicates the use of the 8 micron emission as a star formation rate indicator, at least for small galactic regions and individual star forming regions. By leveraging the Spitzer legacy, this investigation paves the way for future explorations with the James Webb Space Telescope.

The $^{26}$Al short-lived radioactive nuclide is the source of the observed galactic diffuse $\gamma$-ray emission at 1.8 MeV. While different sources of $^{26}$Al have been explored, such as AGB stars, massive stars winds, and supernovae, the contribution of very massive stars has never been studied. We study the stellar wind contribution of very massive stars, i.e stars with initial masses between 150 and 300 M$_\odot$, to the enrichment in $^{26}$Al of the galactic interstellar medium. We discuss the production of $^{26}$Al by studying rotating and non-rotating very massive stellar models with initial masses between 150 and 300 M$_\odot$ for metallicities Z=0.006, 0.014, and 0.020. We confront this result to a simple Milky Way model taking into account both the metallicity and the star formation rate gradients. We obtain that very massive stars in the Z=0.006-0.020 metallicity range might be very significant contributors to the $^{26}$Al enrichment of the interstellar medium. Typically, the contribution of the winds of massive stars to the total quantity of $^{26}$Al in the Galaxy increases by 150\% when very massive stars are considered. Very massive stars, despite their rarity, might be important contributors to $^{26}$Al and overall very important actors for nucleosynthesis in the Galaxy.

In this paper, we carry out a semi-analytic general relativistic study of a Gamma-Ray Bursts (GRB) jet that is breaking out of a cocoon or stellar envelope. We solve hydrodynamic equations with the relativistic equation of state that takes care of fluid composition. In short GRBs, a general relativistic approach is required to account for curved spacetime in strong gravity. The piercing of the jet through the cocoon resembles a de Laval nozzle and the jet may go through recollimation shock transitions. We show that the possibility of shock transition and the shock properties are sensitive to the matter composition and the cocoon strength. Obtained Lorentz factors in thermally driven jets comfortably reach a few $\times$10.

We report observations of a unique, large prominence eruption that was observed in the He II 304 {\AA} passband of the the Extreme Ultraviolet Imager/Full Sun Imager telescope aboard Solar Orbiter on 15-16 February 2022. Observations from several vantage points (Solar Orbiter, the Solar-Terrestrial Relations Observatory, the Solar and Heliospheric Observatory, and Earth-orbiting satellites) were used to measure the kinematics of the erupting prominence and the associated coronal mass ejection. Three-dimensional reconstruction was used to calculate the deprojected positions and speeds of different parts of the prominence. Observations in several passbands allowed us to analyse the radiative properties of the erupting prominence. The leading parts of the erupting prominence and the leading edge of the corresponding coronal mass ejection propagate at speeds of around 1700 km/s and 2200 km/s, respectively, while the trailing parts of the prominence are significantly slower (around 500 km/s). Parts of the prominence are tracked up to heights of over $6 R_\sun$. The He II emission is probably produced via collisional excitation rather than scattering. Surprisingly, the brightness of a trailing feature increases with height. The reported prominence is the first observed in He II 304 {\AA} emission at such a great height (above 6 $R_\sun$).

An intense type III radio storm has been disrupted by a fast halo coronal mass ejection (CME) on 2000 April 4. The CME is also associated with a large solar energetic particle (SEP) event. The storm recovers after about10 hrs. We identified another CME that occurs on 2003 November 11 with similar CME properties but there is no type III storm in progress. The 2003 November 11 CME is also not associated with an SEP event above the background (less than 2 pfu), whereas the one with type III storm has an intense SEP event (about 56 pfu). One of the factors affecting the intensity of SEP events is the presence of seed particles that are accelerated by CME-driven shocks. We suggest that the type III storm source, which accelerates electrons to produce the storm, also accelerates ions that serve as seed particles to the CME shock.

A pseudo-Nambu Goldstone Boson (pNGB) arising from the breaking of a global symmetry ($G\rightarrow H$) can be one of the most sound candidate to describe the thawing quintessence model, which explains the late time acceleration of our universe. Motivated from the Composite Higgs scenario, we have investigated the case where the pNGB associated with $SO(N)/ SO(N-1)$ develops a potential through its couplings with the particles that do not form the complete representations of $G$. The Coleman Weinberg (CW) potential is generated via the external particles in the loop which are linked with the strongly interacting dynamics and can be computed predicatively. The model of Dark Energy (DE) is tested against several latest cosmological observations such as supernovae data of Pantheon, Baryon Acoustic Oscillation (BAO) data, Redshift-space distortion (RSD) data etc. We have found that the fit prefers sub-Planckian value of the pNGB field decay constant. Moreover, we have found that the model predicts cosmological parameters well within the allowed range of the observation and thus gives a well motivated model of quintessence.

We study two of the most theoretically promising models of inflation, namely the Natural inflation and the Mutated Hilltop inflation, in the Einstein-Gauss Bonnet(EGB) gravity framework. In this work, we try to explore these models in EGB framework, keeping the observational constraint from $GW170817$ on the speed of gravitational wave to be equal to the speed of light. This has direct implication on the non-minimal coupling to the Gauss-Bonnet invariant in the action. Thus, the effective potential gets new features. We have not only analysed the inflationary dynamics, but also the reheating dynamics and finally the corresponding energy spectrum of the gravitational wave.

Globular clusters contain a unique pulsar population, with many exotic systems that can form only in their dense stellar environments. The leap in sensitivity of the upgraded Giant Metrewave Radio Telescope (uGMRT) in India, especially at low radio frequencies ($<$ 1 GHz) has motivated a new search for radio pulsars in a group of eight Southern globular clusters. We discovered PSR J1835$-$3259B, a 1.83-ms pulsar in NGC 6652; this is in a near-circular wide orbit of 28.7 hr with a low-mass ($ \sim 0.2 \, M_{\rm \odot}$) companion, likely a Helium white dwarf. We derived a 10-year timing solution for this system. We also present measurements of scattering, flux densities and spectral indices for some of the previously known pulsars in these GCs. A significant fraction of the pulsars in these clusters have steep spectral indices. Additionally, we detected eight radio point sources not associated with any known pulsar positions in the radio images. There are four newly identified sources, three in NGC 6652 and one in NGC 6539, and one previously identified source each in NGC 1851, NGC 6440, NGC 6544, and Terzan 5. Surprisingly, our images show that our newly discovered pulsar, PSR J1835$-$3259B, is the brightest pulsar in all GCs we have imaged; like other pulsars with broad profiles (Ter 5 C and O), its flux density in the radio images is much larger than in its pulsations. This indicates that their pulsed emission is only a fraction of their total emission. The detection of radio sources outside the core radii but well within the tidal radii of these clusters show that future GC surveys should complement the search analysis by using the imaging capability of interferometers, and preferentially synthesize large number of search beams in order to obtain a larger field of view.

The motion of a scalar field that interacts with a hot plasma, like the inflaton during reheating, is damped, which is a dissipative process. At high temperatures the damping can be described by a local term in the effective equation of motion. The damping coefficient is sensitive to multiple scattering. In the loop expansion its computation would require an all-order resummation. Instead we solve an effective Boltzmann equation, similarly to the computation of transport coefficients. For an interaction with another scalar field we obtain a simple relation between the damping coefficient and the bulk viscosity, so that one can make use of known results for the latter. The numerical prefactor of the damping coefficient turns out to be rather large, of order $ 10 ^ 4 $.

The sensitivities of ground-based gravitational-wave (GW) detectors are limited by quantum shot noise at a few hundred Hertz and above. Nonetheless, one can use a quantum-correlation technique proposed by Martynov, et al. [Phys. Rev. A 95, 043831 (2017)] to remove the expectation value of the shot noise, thereby exposing underlying classical signals in the cross spectrum formed by cross-correlating the two outputs in a GW interferometer's anti-symmetric port. We explore here the prospects and analyze the sensitivity of using quantum correlation to detect astrophysical GW signals. Conceptually, this technique is similar to the correlation of two different GW detectors as it utilizes the fact that a GW signal will be correlated in the two outputs but the shot noise will be uncorrelated. Quantum correlation also has its unique advantages as it requires only a single interferometer to make a detection. Therefore, quantum correlation could increase the duty cycle, enhance the search efficiency, and enable the detection of highly polarized signals. In particular, we show that quantum correlation could be especially useful for detecting post-merger remnants of binary neutron stars with both short ($< 1\,{\rm s}$) and intermediate ($\sim 10-10^4\,{\rm s}$) durations and setting upper limits on continuous emissions from unknown pulsars.

Uncertainty quantification is crucial for assessing the predictive ability of AI algorithms. A large body of work (including normalizing flows and Bayesian neural networks) has been devoted to describing the entire predictive distribution (PD) of a target variable Y given input features $\mathbf{X}$. However, off-the-shelf PDs are usually far from being conditionally calibrated; i.e., the probability of occurrence of an event given input $\mathbf{X}$ can be significantly different from the predicted probability. Most current research on predictive inference (such as conformal prediction) concerns constructing prediction sets, that do not only provide correct uncertainties on average over the entire population (that is, averaging over $\mathbf{X}$), but that are also approximately conditionally calibrated with accurate uncertainties for individual instances. It is often believed that the problem of obtaining and assessing entire conditionally calibrated PDs is too challenging to approach. In this work, we show that recalibration as well as validation are indeed attainable goals in practice. Our proposed method relies on the idea of regressing probability integral transform (PIT) scores against $\mathbf{X}$. This regression gives full diagnostics of conditional coverage across the entire feature space and can be used to recalibrate misspecified PDs. We benchmark our corrected prediction bands against oracle bands and state-of-the-art predictive inference algorithms for synthetic data, including settings with distributional shift and dependent high-dimensional sequence data. Finally, we demonstrate an application to the physical sciences in which we assess and produce calibrated PDs for measurements of galaxy distances using imaging data (i.e., photometric redshifts).

In this paper we use a Modified Newton's method based on the Continuous analog of Newton's method and high precision arithmetic for a general numerical search of periodic orbits for the planar three-body problem. We consider relatively short periods and a relatively coarse search-grid. As a result, we found 123 periodic solutions belonging to 105 new topological families that are not included in the database in [Science China Physics, Mechanics and Astronomy 60.12 (2017)]. The extensive numerical search is achieved by a parallel solving of many independent tasks using many cores in a computational cluster.

The running of the Higgs self coupling may lead to numerous phenomena in early universe cosmology. In this paper we introduce a scenario where the Higgs running induces turns in the trajectory passing a region with tachyonic mass, leading to a temporal tachyonic growth in the curvature power spectrum. This effect induced by the Higgs leaves phenomena in the form of primordial black holes and stochastic gravitational waves, where proposed GW observatories will be able to probe in the near future.

The first measurement on the temperature of hydrogen 21-cm signal reported by EDGES strongly favors Coulomb-like interaction between WIMPless dark matter and baryon fluid. We investigate such dark matter both in one- and two-component context, with the light force carrier(s) essential for the Coulomb-like interaction not being photon. Using a conversion of cross sections used by relevant experiments and Boltzmann equations to encode effects of the dark matter-baryon interaction, we show that both cases are robustly excluded by the stringent stellar cooling bounds in the sub-GeV dark matter mass range. The exclusion of one-component case applies to simplified WIMPless dark matter with the light force carrier as dark photon, gauged $B-L$, $L_{e}-L_{\mu}$ or $L_{e}-L_{\tau}$, or axion-like particle, while the exclusion of two-component case applies to simplified WIMPless dark matter with the two light force carriers as two axion-like particles coupled to standard model quarks and leptons respectively.

The ringdown signal emitted during a binary black hole coalescence can be modeled as a linear superposition of the characteristic damped modes of the remnant black hole that get excited during the merger phase. While checking the consistency of the measured frequencies and damping times against the Kerr BH spectrum predicted by General Relativity~(GR) is a cornerstone of strong-field tests of gravity, the consistency of measured excitation amplitudes and phases have been largely left unexplored. For a nonprecessing, quasi-circular binary black hole merger, we find that GR predicts a narrow region in the space of mode amplitude ratio and phase difference, independently of the spin of the binary components. % Using this unexpected result, we develop a new null test of strong-field gravity which demands that the measured amplitudes and phases of different ringdown modes should lie within this narrow region predicted by GR. We call this the \emph{amplitude-phase consistency test} and introduce a procedure for performing it using information from the ringdown signal. Lastly, we apply this test to the GW190521 event, using the multimodal ringdown parameters inferred by Capano et al.~(2021)~\cite{Capano:2021etf}. While ringdown measurements errors for this event are large, we show that GW190521 is consistent with the amplitude-phase consistency test. Our test is particularly well suited for accommodating multiple loud ringdown detections as those expected in the near future, and can be used complementarily to standard black-hole spectroscopy as a proxy for modified gravity, compact objects other than black holes, and binary precession.

We study the production of heavy, $\mu \gtrsim 1$ TeV, bosonic spin $s=0,1$ dark matter (DM) via the simultaneous processes of Hawking evaporation and superradiance (SR) from an initial population of small, $\lesssim 10^6$ kg, primordial black holes (PBHs). Even for small initial PBH spins the SR process can produce extremely dense gravitationally-bound DM Bose or Proca soliton "stars" of radius $\lesssim {\rm pm}$ and mass $\sim 10^{\rm few}$ kg that can survive to today, well after PBH decay. These solitons can constitute a significant fraction of the DM density, rising to $\gtrsim 50\%$ in the vector DM case.