Papers for Monday, Sep 19 2022

RAPOC (Rosseland and Planck Opacity Converter) is a Python 3 code that calculates Rosseland and Planck mean opacities (RPMs) from wavelength-dependent opacities for a given temperature, pressure, and wavelength range. In addition to being user-friendly and rapid, RAPOC can interpolate between discrete data points, making it flexible and widely applicable to the astrophysical and Earth-sciences fields, as well as in engineering. For the input data, RAPOC can use ExoMol and DACE data, or any user-defined data, provided that it is in a readable format. In this paper, we present the RAPOC code and compare its calculated Rosseland and Planck mean opacities with other values found in the literature. The RAPOC code is open-source and available on Pypi and GitHub.

Achernar, the closest and brightest classical Be star, presents rotational flattening, gravity darkening, occasional emission lines due to a gaseous disk, and an extended polar wind. It is also a member of a close binary system with an early A-type dwarf companion. We aim to determine the orbital parameters of the Achernar system and to estimate the physical properties of the components. We monitored the relative position of Achernar B using a broad range of high angular resolution instruments of the VLT/VLTI (VISIR, NACO, SPHERE, AMBER, PIONIER, GRAVITY, and MATISSE) over a period of 13 years (2006-2019). These astrometric observations are complemented with a series of more than 700 optical spectra for the period from 2003 to 2016. We determine that Achernar B orbits the Be star on a seven-year period, eccentric orbit (e = 0.7255 +/- 0.0014) which brings the two stars within 2 au at periastron. The mass of the Be star is found to be mA = 6.0 +/- 0.6 Msun for a secondary mass of mB = 2.0 +/- 0.1 Msun. We find a good agreement of the parameters of Achernar A with the evolutionary model of a critically rotating star of 6.4 Msun at an age of 63 million years. We also identify a resolved comoving low-mass star, which leads us to propose that Achernar is a member of the Tucana-Horologium moving group. Achernar A is presently in a short-lived phase of its evolution following the turn-off, during which its geometrical flattening ratio is the most extreme. Considering the orbital parameters, no significant interaction occurred between the two components, demonstrating that Be stars may form through a direct, single-star evolution path without mass transfer. Since component A will enter the instability strip in a few hundred thousand years, Achernar appears to be a promising progenitor of the Cepheid binary systems.

Stars that interact with supermassive black holes (SMBHs) can either be completely or partially destroyed by tides. In a partial tidal disruption event (TDE) the high-density core of the star remains intact, and the low-density, outer envelope of the star is stripped and feeds a luminous accretion episode. The TDE AT2018fyk, with an inferred black hole mass of $10^{7.7\pm0.4}$ M$_{\odot}$, experienced an extreme dimming event at X-ray (factor of $>$6000) and UV (factor $\sim$15) wavelengths $\sim$500--600 days after discovery. Here we report on the re-emergence of these emission components roughly 1200 days after discovery. We find that the source properties are similar to those of the pre-dimming accretion state, suggesting that the accretion flow was rejuvenated to a similar state. We propose that a repeating partial TDE, where the partially disrupted star is on a $\sim 600$ day orbit about the SMBH and is periodically stripped of mass during each pericenter passage, powers its unique lightcurve. This scenario provides a plausible explanation for AT2018fyk's overall properties, including the rapid dimming event and the rebrightening at late times. We also provide testable predictions for the behavior of the accretion flow in the future: if the second encounter was also a partial disruption then we predict another strong dimming event around day 1800 (August 2023), and a subsequent rebrightening around day 2400 (March 2025). This source provides strong evidence of the partial disruption of a star by a SMBH.

We study the formation of ultra-diffuse galaxies (UDGs) using the cosmological hydrodynamical simulation TNG50 of the Illustris-TNG suite. We define UDGs as dwarf galaxies in the stellar mass range $\rm{7.5 \leq log (M_{\star} / M_{\odot}) \leq 9 }$ that are in the $5\%$ most extended tail of the simulated mass-size relation. This results in a sample of UDGs with half-mass radii $\rm{r_{h \star } \gtrsim 2 \ kpc}$ and surface brightness between $\rm{24.5}$ and $\rm{28 \ mag \ arcsec^{-2}}$, similar to definitions of UDGs in observations. The large cosmological volume in TNG50 allows for a comparison of UDGs properties in different environments, from the field to galaxy clusters with virial mass $\rm{M_{200} \sim 2 \times 10^{14} ~ M_{\odot}}$. All UDGs in our sample have dwarf-mass haloes ($\rm{M_{200}\sim 10^{11} ~ M_{\odot} }$) and show the same environmental trends as normal dwarfs: field UDGs are star-forming and blue while satellite UDGs are typically quiescent and red. The TNG50 simulation predicts UDGs that populate preferentially higher spin haloes and more massive haloes at fixed $\rm{M_{\star}}$ compared to non-UDG dwarfs. This applies also to most satellite UDGs, which are actually ``born" UDGs in the field and infall into groups and clusters without significant changes to their size. We find, however, a small subset of satellite UDGs ($\lesssim 10 \%$) with present-day stellar size a factor $\geq 1.5$ larger than at infall, confirming that tidal effects, particularly in the lower mass dwarfs, are also a viable formation mechanism for some of these dwarfs, although subdominant in this simulation.

The Multi Unit Spectroscopic Explorer (MUSE) is an integral field spectrograph on the Very Large Telescope Unit Telescope 4, capable of laser guide star assisted and tomographic adaptive optics using the GALACSI module. Its observing capabilities include a wide field (1 square arcmin), ground layer AO mode (WFM-AO) and a narrow field (7.5"x7.5"), laser tomography AO mode (NFM-AO). The latter has had several upgrades in the 4 years since commissioning, including an optimisation of the control matrices for the AO system and a new sub-electron noise detector for its infra-red low order wavefront sensor. We set out to quantify the NFM-AO system performance by analysing $\sim$230 spectrophotometric standard star observations taken over the last 3 years. To this end we expand upon previous work, designed to facilitate analysis of the WFM-AO system performance. We briefly describe the framework that will provide a user friendly, semi-automated way for system performance monitoring during science operations. We provide the results of our performance analysis, chiefly through the measured Strehl ratio and full width at half maximum (FWHM) of the core of the point spread function (PSF) using two PSF models, and correlations with atmospheric conditions. These results will feed into a range of applications, including providing a more accurate prediction of the system performance as implemented in the exposure time calculator, and the associated optimization of the scientific output for a given set of limiting atmospheric conditions.

The modern time-domain photometric surveys collect a lot of observations of various astronomical objects, and the coming era of large-scale surveys will provide even more information. Most of the objects have never received a spectroscopic follow-up, which is especially crucial for transients e.g. supernovae. In such cases, observed light curves could present an affordable alternative. Time series are actively used for photometric classification and characterization, such as peak and luminosity decline estimation. However, the collected time series are multidimensional, irregularly sampled, contain outliers, and do not have well-defined systematic uncertainties. Machine learning methods help extract useful information from available data in the most efficient way. We consider several light curve approximation methods based on neural networks: Multilayer Perceptrons, Bayesian Neural Networks, and Normalizing Flows, to approximate observations of a single light curve. Tests using both the simulated PLAsTiCC and real Zwicky Transient Facility data samples demonstrate that even few observations are enough to fit networks and achieve better approximation quality than other state-of-the-art methods. We show that the methods described in this work have better computational complexity and work faster than Gaussian Processes. We analyze the performance of the approximation techniques aiming to fill the gaps in the observations of the light curves, and show that the use of appropriate technique increases the accuracy of peak finding and supernova classification. In addition, the study results are organized in a Fulu Python library available on GitHub, which can be easily used by the community.

The existence of eV-mass sterile neutrinos is not ruled out because of persistent experimental anomalies. Upcoming multi-messenger detections of neutron-star merger remnants could provide indirect constraints on the existence of these particles. We explore the active-sterile flavor conversion phenomenology in a two-flavor scenario (1 active + 1 sterile species) as a function of the sterile neutrino mixing parameters, neutrino emission angle from the accretion torus, and temporal evolution of the merger remnant. The torus geometry and the neutron richness of the remnant are responsible for the occurrence of multiple resonant active-sterile conversions. The number of resonances strongly depends on the neutrino emission direction above or inside the remnant torus and leads to large production of sterile neutrinos (and no antineutrinos) in the proximity of the polar axis as well as more sterile antineutrinos than neutrinos in the equatorial region. As the black hole torus evolves in time, the shallower baryon density is responsible for more adiabatic flavor conversion, leading to larger regions of the mass-mixing parameter space being affected by flavor mixing. Our findings imply that the production of sterile states can have indirect implications on the disk cooling rate, its outflows, and related electromagnetic observables which remain to be assessed.

Simple monomial inflationary scenarios have been ruled out by recent observations. In this work we revisit the next simplest scenario, a single--field model where the scalar potential is a polynomial of degree four which features a concave ``almost'' saddle point. We focus on trans--Planckian field values. We reparametrize the potential, which greatly simplifies the procedure for finding acceptbale model parameters. This allows for the first comprehensive scan of parameter space consistent with recent Planck and BICEP/Keck 2018 measurements. Even for trans--Planckian field values the tensor--to--scalar ratio $r$ can be as small as $\mathcal{O}(10^{-8})$, but the model can also saturate the current upper bound. In contrast to the small--field version of this model, radiative stability does not lead to strong constraints on the parameters of the inflaton potential. For very large field values the potential can be approximated by the quartic term; as well known, this allows eternal inflation even for field energy well below the reduced Planck mass $M_{\rm Pl}$, with Hubble parameter $H \sim 10^{-2} M_{\rm Pl}$. More interestingly, we find a region of parameter space that even supports {\em two phases of eternal inflation}. The second epoch only occurs if the slope at the would--be saddle point is very small, and has $H \sim 10^{-5} M_{\rm Pl}$; it can only be realized if $r \sim 10^{-2}$, within the sensitivity range of next--generation CMB observations.

The Great Oxidation Event was a period during which Earth's atmospheric oxygen (O$_2$) concentrations increased from $\sim 10^{-5}$ times its present atmospheric level (PAL) to near modern levels, marking the start of the Proterozoic geological eon 2.4 billion years ago. Using WACCM6, an Earth System Model, we simulate the atmosphere of Earth-analogue exoplanets with O$_2$ mixing ratios between 0.1% and 150% PAL. Using these simulations, we calculate the reflection/emission spectra over multiple orbits using the Planetary Spectrum Generator. We highlight how observer angle, albedo, chemistry, and clouds affect the simulated observations. We show that inter-annual climate variations, as well short-term variations due to clouds, can be observed in our simulated atmospheres with a telescope concept such as LUVOIR or HabEx. Annual variability and seasonal variability can change the planet's reflected flux (including the reflected flux of key spectral features such as O$_2$ and H$_2$O) by up to factors of 5 and 20, respectively, for the same orbital phase. This variability is best observed with a high-throughput coronagraph. For example, HabEx (4 m) with a starshade performs up to a factor of two times better than a LUVOIR B (6 m) style telescope. The variability and signal-to-noise ratio of some spectral features depends non-linearly on atmospheric O$_2$ concentration. This is caused by temperature and chemical column depth variations, as well as generally increased liquid and ice cloud content for atmospheres with O$_2$ concentrations of $<$1% PAL.

Little is known about the radio astronomical universe at frequencies below 10 MHz because such radiation does not penetrate the ionosphere. A cubesat-based telescope for the 1--10 MHz band could be rapidly and economically deployed in low earth orbit. We consider possible transient and steady sources, and application to study of the ionosphere itself.

We present a database of solar flares registered by the Konus-Wind instrument during more than 27 years of operation, from 1994 November to now (2022 June). The constantly updated database (hereafter KW-Sun) contains over 1000 events detected in the instrument's triggered mode and is accessible online at this http URL For each flare, the database provides time-resolved energy spectra in energy range from ~20 keV to ~15 MeV in FITS format along with count rate light curves in three wide energy bands G1 (~20-80 keV), G2 (~80-300 keV), and G3 (~300-1200 keV) with high time resolution (down to 16 ms) in ASCII and IDL SAV formats. This article focuses on the instrument capabilities in the context of solar observations, the structure of the KW-Sun data and their intended usage. The presented homogeneous data set obtained in the broad energy range with high temporal resolution during more than two full solar cycles is beneficial for both statistical and case studies as well as a source of context data for solar flare research.

The measurement of the large scale distribution of neutral hydrogen in the late Universe, obtained with radio telescopes through the hydrogen 21cm line emission, has the potential to become a key cosmological probe in the upcoming years. We explore the constraining power of 21cm intensity mapping observations on the full set of cosmological parameters that describe the $\Lambda$CDM model. We assume a single-dish survey for the SKA Observatory and simulate the 21cm linear power spectrum monopole and quadrupole within six redshift bins in the range $z=0.25-3$. Forecasted constraints are computed numerically through Markov Chain Monte Carlo techniques. We extend the sampler \texttt{CosmoMC} by implementing the likelihood function for the 21cm power spectrum multipoles. We assess the constraining power of the mock data set alone and combined with Planck 2018 CMB observations. We include a discussion on the impact of extending measurements to non-linear scales in our analysis. We find that 21cm multipoles observations alone are enough to obtain constraints on the cosmological parameters comparable with other probes. Combining the 21cm data set with CMB observations results in significantly reduced errors on all the cosmological parameters. The strongest effect is on $\Omega_ch^2$ and $H_0$, for which the error is reduced by almost a factor four. The percentage errors we estimate are $\sigma_{\Omega_ch^2} = 0.25\%$ and $\sigma_{H_0} = 0.16\%$, to be compared with the Planck only results $\sigma_{\Omega_ch^2} = 0.99\%$ and $\sigma_{H_0} = 0.79\%$. We conclude that 21cm SKAO observations will provide a competitive cosmological probe, complementary to CMB and, thus, pivotal for gaining statistical significance on the cosmological parameters constraints, allowing a stress test for the current cosmological model.

In a companion paper, a faceted wideband imaging technique for radio interferometry, dubbed Faceted HyperSARA, has been introduced and validated on synthetic data. Building on the recent HyperSARA approach, Faceted HyperSARA leverages the splitting functionality inherent to the underlying primal-dual forward-backward algorithm to decompose the image reconstruction over multiple spatio-spectral facets. The approach allows complex regularization to be injected into the imaging process while providing additional parallelization flexibility compared to HyperSARA. The present paper introduces new algorithm functionalities to address real datasets, implemented as part of a fully fledged MATLAB imaging library made available on Github. A large scale proof-of-concept is proposed to validate Faceted HyperSARA in a new data and parameter scale regime, compared to the state-of-the-art. The reconstruction of a 15 GB wideband image of Cyg A from 7.4 GB of VLA data is considered, utilizing 1440 CPU cores on a HPC system for about 9 hours. The conducted experiments illustrate the reconstruction performance of the proposed approach on real data, exploiting new functionalities to set, both an accurate model of the measurement operator accounting for known direction-dependent effects (DDEs), and an effective noise level accounting for imperfect calibration. They also demonstrate that, when combined with a further dimensionality reduction functionality, Faceted HyperSARA enables the recovery of a 3.6 GB image of Cyg A from the same data using only 91 CPU cores for 39 hours. In this setting, the proposed approach is shown to provide a superior reconstruction quality compared to the state-of-the-art wideband CLEAN-based algorithm of the WSClean software.

Galaxy mergers trigger both star formation and accretion onto the central supermassive black hole. As a result of subsequent energetic feedback processes, it has long been proposed that star formation may be promptly extinguished in galaxy merger remnants. However, this prediction of widespread, rapid quenching in late stage mergers has been recently called into question with modern simulations and has never been tested observationally. Here we perform the first empirical assessment of the long-predicted end phase in the merger sequence. Based on a sample of ~500 post-mergers identified from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), we show that the frequency of post-merger galaxies that have rapidly shutdown their star formation following a previous starburst is 30-60 times higher than expected from a control sample of non-merging galaxies. No such excess is found in a sample of close galaxy pairs, demonstrating that mergers can indeed lead to a rapid halt to star formation, but that this process only manifests after coalescence.

Coronagraphs allow for faint off-axis exoplanets to be observed, but are limited to angular separations greater than a few beam widths. Accessing closer-in separations would greatly increase the expected number of detectable planets, which scales inversely with the inner working angle. The Vortex Fiber Nuller (VFN) is an instrument concept designed to characterize exoplanets within a single beam-width. It requires few optical elements and is compatible with many coronagraph designs as a complementary characterization tool. However, the peak throughput for planet light is limited to about 20%, and the measurement places poor constraints on the planet location and flux ratio. We propose to augment the VFN design by replacing its single-mode fiber with a six-port mode-selective photonic lantern, retaining the original functionality while providing several additional ports that reject starlight but couple planet light. We show that the photonic lantern can also be used as a nuller without a vortex. We present monochromatic simulations characterizing the response of the Photonic Lantern Nuller (PLN) to astrophysical signals and wavefront errors, and show that combining exoplanet flux from the nulled ports significantly increases the overall throughput of the instrument. We show using synthetically generated data that the PLN detects exoplanets more effectively than the VFN. Furthermore, with the PLN, the exoplanet can be partially localized, and its flux ratio constrained. The PLN has the potential to be a powerful characterization tool complementary to traditional coronagraphs in future high-contrast instruments.

In this article, we provide an alternative up-sampling and PSF deconvolution method for the iterative multi-exposure coaddition. Different from the previous works, the new method has a ratio-correction term, which allows the iterations to converge more rapidly to an accurate representation of the underlying image than those with difference-correction terms. By employing this method, one can coadd the under-sampled multi-exposures to a super-resolution and obtain a higher peak signal-to-noise ratio. A set of simulations show that we can take many advantages of the new method, e.g. in the signal-to-noise ratio, the average deviation of all source fluxes, super-resolution, and source distortion ratio, etc., which are friendly to astronomical photometry and morphology, and benefits faint source detection and shear measurement of weak gravitational lensing. It provides an improvement in fidelity over the previous works tested in this paper.

We present a new HI mass estimator which relates the HI-to-stellar mass ratio to four galaxy properties: stellar surface mass density, color index $u-r$, stellar mass and concentration index, with the scatter of individual galaxies around the mean HI mass modeled with a Gaussian distribution. We calibrate the estimator using the xGASS sample, including both HI detection and non-detection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased HI masses for optical samples like the SDSS, thus suitable for statistical studies of HI gas contents in galaxies and dark matter halos. We apply our estimator to the SDSS spectroscopic sample to estimate the local HI mass function (HIMF), the conditional HI mass function (CHIMF) in galaxy groups and the HI-halo mass (HIHM) relation. Our HIMF agrees with the ALFALFA measurements at $M_{HI}\gtrsim 5\times 10^9M_{\odot}$, but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance which is corrected for the SDSS sample but not for the ALFALFA. The CHIMFs for all halo masses can be described by a single Schechter function, and this is true for red, blue and satellite galaxies. For central galaxies the CHIMFs show a double-Gaussian profile, with the two components contributed by the red and blue galaxies, respectively. The total HI mass in a group increases monotonically with halo mass. The HI mass of central galaxies in galaxy groups increases rapidly with halo mass only at $M_h\lesssim10^{12}M_{\odot}$, while the mass dependence becomes much weaker at higher halo masses. The observed HI-halo mass relation is not reproduced by current hydrodynamic simulations and semi-analytic models of galaxy formation.

We investigate the second order gravitational waves induced by the primordial scalar and tensor perturbations during radiation-dominated era. The explicit expressions of the power spectra of the second order GWs are presented. We calculate the energy density spectra of the second order GWs for a monochromatic primordial power spectra. For large $k$ $\left( k>k_* \right)$, the effects of the primordial tensor perturbation with tensor-to-scalar ratio $r=A_{h}/A_{\zeta}=0.2$ lead to an around $50\% $ increase of the signal-to-noise ratio (SNR) for LISA observations.

Analyzing the effects of space weather on aviation is a new and developing topic. It has been commonly accepted that the flight time of the polar flights may increase during solar proton events because the flights have to change their route to avoid the high-energy particles. However, apart from such phenomenon, researches related to the flight time during space weather events is very rare. Based on the analyses of 39 representative international air routes around westerlies, it is found that 97.44% (94.87%) of the commercial airplanes on the westbound (eastbound) air routes reveal shorter (longer) flight time during solar proton events compared to those during quiet periods, and the averaged magnitude of change in flight time is ~10 min or 0.21%-4.17% of the total flight durations. Comparative investigations reassure the certainty of such phenomenon that the directional differences in flight time are still incontrovertible regardless of over-land routes (China-Europe) or over-sea routes (China-Western America). Further analyses suggest that the solar proton events associated atmospheric heating will change the flight durations by weakening certain atmospheric circulations, such as the polar jet stream. While the polar jet stream will not be obviously altered during solar flares so that the directional differences in flight time are not found. Besides the conventional space weather effects already known, this paper is the first report that indicates a distinct new scenario of how the solar proton events affect flight time. These analyses are also important for aviation since our discoveries could help the airways optimize the air routes to save passenger time costs, reduce fuel costs and even contribute to the global warming issues.

Using the IllustrisTNG simulation, we study the interaction of large-scale shocks with the circumgalactic medium (CGM) and interstellar medium (ISM) of star-forming (SF) satellite galaxies in galaxy clusters. These shocks are usually produced by mergers and massive accretion. Our visual inspection shows that approximately half of SF satellites have encountered shocks in their host clusters at $z\leq0.11$. After a satellite crosses a shock front and enters the postshock region, the ram pressure on it is boosted significantly. Both the CGM and ISM can be severely impacted, either by striping or compression. The stripping of the ISM is particularly important for low-mass galaxies with $\log (M_{*}/M_{\odot})<10$ and can occur even in the outskirts of galaxy clusters. In comparison, satellites that do not interact with shocks lose their ISM only in the inner regions of clusters. About half of the ISM is stripped within about 0.6 Gyr after it crosses the shock front. Our results show that shock-induced stripping plays an important role in quenching satellite galaxies in clusters.

BP Psc is an active late-type (sp:G9) star with unclear evolutionary status lying at high galactic latitude $b=-57^{\circ}$. It is also the source of the well collimated bipolar jet. We present results of the proper motion and radial velocity study of BP Psc outflow based on the archival $H\alpha$ imaging with the GMOS camera at 8.1-m Gemini-North telescope as well as recent imaging and long-slit spectroscopy with the SCORPIO multi-mode focal reducer at 6-m BTA telescope of SAO RAS. The 3D kinematics of the jet revealed the full spatial velocity up to $\sim$140 km$\cdot$s$^{-1}$ and allows us to estimate the distance to BP Psc system as $D=135\pm40$ pc. This distance leads to an estimation of the central source luminosity $L_*\approx1.2L_{\odot}$, indicating that it is the $\approx$1.3$M_{\odot}$ T Tauri star with an age $t\lesssim$ 7 Myr. We measured the electron density of order $N_e\sim10^2$ cm$^{-3}$ and mean ionization fraction $f\approx0.04$ within the jet knots and estimated upper limit of the mass-loss rate in NE lobe as $\dot{M}_{out}\approx1.2\cdot10^{-8}M_{\odot}\cdot yr^{-1}$. The physical characteristics of the outflow are typical for the low-excitation YSO jets and consistent with the magnetocentrifugal mechanism of its launching and collimation. Prominent wiggling pattern revealed in $H\alpha$ images allowed us to suppose the existence of a secondary substellar companion in a non-coplanar orbit and estimate its most plausible mass as $M_p\approx 30M_{Jup}$. We conclude that BP Psc is one of the closest to the Sun young jet-driving systems and its origin is possibly related to the episode of star formation triggered by expanding supershells in Second Galactic quadrant.

MAXI J1421-613 is an X-ray burster discovered by Monitor of All-sky X-ray Image (MAXI) on 9 January 2014 and is considered to be a low-mass X-ray binary. A previous study analyzing follow-up observation data obtained by Suzaku on 31 January to 3 February 2014 reported that an annular emission of ~3'-9' radius was found around the transient source. The most plausible origin of the annular emission is a dust scattering echo by the outburst of MAXI J1421-613. In this paper, we confirm the annular emission by analyzing the data of the Swift follow-up observation which was conducted by the photon counting mode on 18 January 2014. In a radial profile, we found an annular emission at ~2'.5-4'.5. Its spectrum was well explained by an absorbed power law, and the photon index was higher than that of MAXI J1421-613 itself by delta Gamma~2. The flux and radius of the annular emission observed by Swift are explained by dust scattering of the same outburst as is responsible for the annular emission observed by Suzaku. Assuming that the dust layer causing the annular emission found by Swift is located at the same position as the CO cloud in front of MAXI J1421-613, the distance to the transient source was estimated to be ~3 kpc, which is consistent with the value estimated by the previous study of Suzaku.

One of the major challenges of mid-infrared astronomical heterodyne interferometry is its sensitivity limitations. Detectors capable of handling several 10 GHz bandwidths have been identified as key building blocks of future instruments. Intersubband detectors based on heterostructures have recently demonstrated their ability to provide such performances. In this work we characterize a Quantum Well Infrared Photodetector in terms of noise, dynamic range and bandwidth in a non-interferometric heterodyne set-up. We discuss the possibility to use them on astronomical systems to measure the beating between the local oscillator and the astronomical signal.

Statistical studies of X-ray selected Active Galactic Nuclei (AGN) indicate that the fraction of obscured AGN increases with increasing redshift, and the results suggest that a significant part of the accretion growth occurs behind obscuring material in the early universe. We investigate the obscured fraction of highly accreting X-ray AGN at around the peak epoch of supermassive black hole growth utilizing the wide and deep X-ray and optical/IR imaging datasets. A unique sample of luminous X-ray selected AGNs above $z>2$ was constructed by matching the XMM-SERVS X-ray point-source catalog with a PSF-convolved photometric catalog covering from $u^*$ to 4.5$\mu \mathrm{m}$ bands. Photometric redshift, hydrogen column density, and 2-10 keV AGN luminosity of the X-ray selected AGN candidates were estimated. Using the sample of 306 2-10 keV detected AGN at above redshift 2, we estimate the fraction of AGN with $\log N_{\rm H}\ (\rm cm^{-2})>22$, assuming parametric X-ray luminosity and absorption functions. The results suggest that $76_{-3}^{+4}\%$ of luminous quasars ($\log L_X\ (\rm erg\ s^{-1}) >44.5$) above redshift 2 are obscured. The fraction indicates an increased contribution of obscured accretion at high redshift than that in the local universe. We discuss the implications of the increasing obscured fraction with increasing redshift based on the AGN obscuration scenarios, which describe obscuration properties in the local universe. Both the obscured and unobscured $z>2$ AGN show a broad range of SEDs and morphology, which may reflect the broad variety of host galaxy properties and physical processes associated with the obscuration.

A tantalizing hint of isotropic cosmic birefringence has been found in the $E B$ cross-power spectrum of the cosmic microwave background (CMB) polarization data with a statistical significance of $3\sigma$. A pseudoscalar field coupled to the CMB photons via the Chern-Simons term can explain this observation. The same field may also be responsible for early dark energy (EDE), which alleviates the so-called Hubble tension. Since the EDE field evolves significantly during the recombination epoch, the conventional formula that relates $E B$ to the difference between the $E$- and $B$-mode auto-power spectra is no longer valid. Solving the Boltzmann equation for polarized photons and the dynamics of the EDE field consistently, we find that currently favored parameter space of the EDE model yields a variety of shapes of the $EB$ spectrum, which can be tested by CMB experiments.

Icy grains in the interstellar medium and star-formation regions consist of a variety of materials. Such composite grains interact differently with cosmic-ray (CR) particles compared to simple single-material grains. We aim to calculate the spectra of energies and temperatures of mixed-composition grains undergoing whole-grain heating by CRs. The grains were assumed to consist of a mixture of carbon and olivine, covered by ices consisting of carbon oxides and water. The energy and temperature spectra for grains with radii 0.05; 0.1, and 0.2 microns impacted by CRs were calculated for eight values of column density, relevant to molecular clouds and star-forming cores. The approach takes into account changes in ice thickness and composition with increasing column density. These detailed data for CR interaction with interstellar grains are intended for applications in astrochemical models. The main finding is that the a more accurate approach on grain heat capacity and other factors prevent a frequent heating of 0.1 micron or larger icy grains to high temperatures.

We numerically explore on galaxy scales the Dipolar dark matter (DM) model based on the concept of gravitational polarization. This DM model has been proposed as a natural way to reproduce observed tight galactic scaling relations such as the baryonic Tully-Fisher relation and the Radial Acceleration Relation. We present a customized version of the \texttt{RAMSES} code including for the first time the dynamics of this Dipolar DM in $N$-body simulations. As a first application of this code, we check that we recover an equilibrium configuration that had been found analytically, where a low density Dipolar DM halo is at rest with respect to its central galaxy, recovering the aforementioned scaling relations. A characteristic signature of this equilibrium model is that it harbours a dynamical instability with a characteristic time depending on the Dipolar DM halo density, which we recover numerically. This represents a first step towards more involved simulations needed to test this framework, ranging from galaxy interactions to structure formation.

We present the results of our studies of deuterated molecules (DCN, DNC, DCO$^+$, N$_2$D$^+$ and NH$_2$D) in regions of high-mass star formation, which include a survey of such regions with the 20-m Onsala radio telescope and mapping of several objects in various lines with the 30-m IRAM and 100-m MPIfR radio telescopes. The deuteration degree reaches $\sim$10$^{-2}$ in these objects. We discuss its dependencies on the gas temperature and velocity dispersion, as well as spatial distributions of deuterated molecules. We show that the H$^{13}$CN/HN$^{13}$C intensity ratio may be a good indicator of the gas kinetic temperature and estimate densities of the investigated objects.

In interferometry, the quality of the reconstructed image depends on the algorithm used and its parameters, and users often need to compare the results of several algorithms to disentangle artifacts from actual features of the astrophysical object. Such comparisons can rapidly become cumbersome, as these software packages are very different. OImaging is a graphical interface intended to be a common frontend to image reconstruction software packages. With OImaging, the user can now perform multiple reconstructions within a single interface. From a given dataset, OImaging allows benchmarking of different image reconstruction algorithms and assessment of the reliability of the image reconstruction process. To that end, OImaging uses the IMAGE-OI OIFITS extension proposed to standardize communication with image reconstruction algorithms.

Detecting the 21-cm signal from the Epoch of Reionisation (EoR) is challenging due to the strong astrophysical foregrounds, ionospheric effects, radio frequency interference and instrumental effects. Understanding and calibrating these effects are crucial for the detection. In this work, we introduce a newly developed direction-dependent (DD) calibration algorithm DDECAL and compare its performance with an existing algorithm, SAGECAL, in the context of the LOFAR-EoR 21-cm power spectrum experiment. In our data set, the North Celestial Pole (NCP) and its flanking fields were observed simultaneously. We analyse the NCP and one of its flanking fields. The NCP field is calibrated by the standard pipeline, using SAGECAL with an extensive sky model and 122 directions, and the flanking field is calibrated by DDECAL and SAGECAL with a simpler sky model and 22 directions. Additionally, two strategies are used for subtracting Cassiopeia A and Cygnus A. The results show that DDECAL performs better at subtracting sources in the primary beam region due to the application of a beam model, while SAGECAL performs better at subtracting Cassiopeia A and Cygnus A. This indicates that including a beam model during DD calibration significantly improves the performance. The benefit is obvious in the primary beam region. We also compare the 21-cm power spectra on two different fields. The results show that the flanking field produces better upper limits compared to the NCP in this particular observation. Despite the minor differences between DDECAL and SAGECAL due to the beam application, we find that the two algorithms yield comparable 21-cm power spectra on the LOFAR-EoR data after foreground removal. Hence, the current LOFAR-EoR 21-cm power spectrum limits are not likely to depend on the DD calibration method.

Based on a large spectroscopic sample of $\sim$ 4,300 RR Lyrae stars with metallicity, systemic radial velocity and distance measurements, we present a detailed analysis of the chemical and kinematic properties of the Galactic halo. Using this sample, the metallicity distribution function (MDF) as a function of $r$ and the velocity anisotropy parameter $\beta$ profiles (for different metallicity populations) are derived for the stellar halo. Both the chemical and kinematic results suggest that the Galactic halo is composed of two distinct parts, the inner halo and outer halo. The cutoff radius ($\sim$ 30 kpc) is similar to the previous break radius found in the density distribution of the stellar halo. We find that the inner part is dominated by a metal-rich population with extremely radial anisotropy ($\beta \sim 0.9$). These features are in accordance with those of ``{\it Gaia}-Enceladus-Sausage'' (GES) and we attribute this inner halo component as being dominantly composed of stars deposited from this ancient merged satellite. We find that GES probably has a slightly negative metallicity gradient. The metal-poor populations in the inner halo are characterized as a long-tail in MDF with an anisotropy of $\beta \sim 0.5$, which is similar to that of the outer part. The MDF for the outer halo is very broad with several weak peaks and the value of $\beta$ is around 0.5 for all metallicities.

Over the last years, large (sub-)millimetre surveys of protoplanetary disks have well constrained the demographics of disks, such as their millimetre luminosities, spectral indices, and disk radii. Additionally, several high-resolution observations have revealed an abundance of substructures in the disks dust continuum. The most prominent are ring like structures, likely due to pressure bumps trapping dust particles. The origins and characteristics of these bumps, nevertheless, need to be further investigated. The purpose of this work is to study how dynamic pressure bumps affect observational properties of protoplanetary disks. We further aim to differentiate between the planetary- versus zonal flow-origin of pressure bumps. We perform one-dimensional gas and dust evolution simulations, setting up models with varying pressure bump features. We subsequently run radiative transfer calculations to obtain synthetic images and the different quantities of observations. We find that the outermost pressure bump determines the disks dust size across different millimetre wavelengths. Our modelled dust traps need to form early (< 0.1 Myr), fast (on viscous timescales), and must be long lived (> Myr) to obtain the observed high millimetre luminosities and low spectral indices of disks. While the planetary bump models can reproduce these observables irrespectively of the opacity prescription, the highest opacities are needed for the zonal flow bump model to be in line with observations. Our findings favour the planetary- over the zonal flow-origin of pressure bumps and support the idea that planet formation already occurs in early class 0-1 stages of circumstellar disks. The determination of the disks effective size through its outermost pressure bump also delivers a possible answer to why disks in recent low-resolution surveys appear to have the same sizes across different millimetre wavelengths.

The changes in brightness of an astronomical source as a function of time are key probes into that source's physics. Periodic and quasi-periodic signals are indicators of fundamental time (and length) scales in the system, while stochastic processes help uncover the nature of turbulent accretion processes. A key method of studying time variability is through Fourier methods, the decomposition of the signal into sine waves, which yields a representation of the data in frequency space. With the extension into \textit{spectral timing} the methods built on the Fourier transform can not only help us characterize (quasi-)periodicities and stochastic processes, but also uncover the complex relationships between time, photon energy and flux in order to help build better models of accretion processes and other high-energy dynamical physics. In this Chapter, we provide a broad, but practical overview of the most important relevant methods.

Direct seismic imaging of sub-surface flow, sound-speed and magnetic field is crucial for predicting flux tube emergence on the solar surface, an important ingredient for space weather. The sensitivity of helioseismic mode-amplitude cross-correlation to $p$- and $f$-mode oscillations enable formal inversion of such sub-photospheric perturbations. It is well-known that such problems are written in the form of an integral equation that connects the perturbations to the observations via ``sensitivity kernels". While the sensitivity kernels for flow and sound-speed have been known for decades and have been used extensively, formulating kernels for general magnetic perturbations had been elusive. A recent study proposed sensitivity kernels for Lorentz-stresses corresponding to global magnetic fields of general geometry. The present study is devoted to proposing kernels for inferring Lorentz-stresses as well as the solenoidal magnetic field in a local patch on the Sun via Cartesian mode-coupling. Moreover, for the first time in solar physics, Slepian functions are employed to parameterize perturbations in the horizontal dimension. This is shown to increase the number of data constraints in the inverse problem, implying an increase in the precision of inferred parameters. This paves the path to reliably imaging sub-surface solar magnetic features in, e.g., supergranules, sunspots and (emerging) active regions.

In the present paper, we use Si IV 1393.755 \AA\ spectral line observed by the Interface Region Imaging Spectrograph (IRIS) in the quiet-Sun to determine physical nature of the solar transition region (TR) oscillations. We analyze the properties of these oscillations using wavelet tools (e.g., power, cross-power, coherence, and phase difference) along with the stringent noise model (i.e., power-law + constant). We estimate the period of the intensity and Doppler velocity oscillations at each chosen location in the quiet-Sun (QS) and quantify the distribution of the statistically significant power and associated periods in one bright and two dark regions. In the bright TR region, the mean periods in intensity and velocity are 7 min, and 8 min respectively. In the dark region, the mean periods in intensity and velocity are 7 min, and 5.4 min respectively. We also estimate the phase difference between the intensity and Doppler velocity oscillations at each location. The statistical distribution of phase difference is estimated, which peaks at -119\degree $\pm$ 13\degree, 33\degree $\pm$ 10\degree, 102\degree $\pm$ 10\degree\ in the bright region, while at -153\degree $\pm$ 13\degree, 6\degree $\pm$ 20\degree, 151\degree $\pm$ 10\degree\ in the dark region. The statistical distribution reveals that the oscillations are caused by propagating slow magnetoacoustic waves encountered with the TR. Some of these locations may also be associated with the standing slow waves. Even, in the given time domain, several locations exhibit presence of both propagating and standing oscillations at different frequencies.

1991T-like supernovae are the luminous, slow-declining extreme of the Branch shallow-silicon (SS) subclass of Type Ia supernovae. They are distinguished by extremely weak Ca II H & K and Si II $\lambda6355$ and strong Fe III absorption features in their optical spectra at pre-maximum phases, and have long been suspected to be over-luminous compared to normal Type Ia supernovae. In this paper, the pseudo equivalent width of the Si II $\lambda$6355 absorption obtained at light curve phases from $\leq+10$ days is combined with the morphology of the $i$-band light curve to identify a sample of 1991T-like supernovae in the Carnegie Supernova Project-II. Hubble diagram residuals show that, at optical as well as near-infrared wavelengths, these events are over-luminous by $\sim$0.1-0.5 mag with respect to the less extreme Branch SS (1999aa-like) and Branch core-normal supernovae with similar $B$-band light curve decline rates.

Cosmological observations in the new millennium have dramatically increased our understanding of the Universe, but several fundamental questions remain unanswered. This topical group report describes the best opportunities to address these questions over the coming decades by extending observations to the $z<6$ universe. The greatest opportunity to revolutionize our understanding of cosmic acceleration both in the modern universe and the inflationary epoch would be provided by a new Stage V Spectroscopic Facility (Spec-S5) which would combine a large telescope aperture, wide field of view, and high multiplexing. Such a facility could simultaneously provide a dense sample of galaxies at lower redshifts to provide robust measurements of the growth of structure at small scales, as well as a sample at redshifts $2<z<5$ to measure cosmic structure at the largest scales, spanning a sufficient volume to probe primordial non-Gaussianity from inflation, to search for features in the inflationary power spectrum on a broad range of scales, to test dark energy models in poorly-explored regimes, and to determine the total neutrino mass and effective number of light relics. A number of compelling opportunities at smaller scales should also be pursued alongside Spec-S5. The science collaborations analyzing DESI and LSST data will need funding for a variety of activities, including cross-survey simulations and combined analyses. The results from these experiments can be greatly improved by smaller programs to obtain complementary data, including follow-up studies of supernovae and spectroscopy to improve photometric redshift measurements. The best future use of the Vera C. Rubin Observatory should be evaluated later this decade after the first LSST analyses have been done. Finally, investments in pathfinder projects could enable powerful new probes of cosmology to come online in future decades.

Interactions of ultrahigh energy cosmic rays in the surroundings of their accelerators can naturally explain the observed spectrum and composition of UHECRs, including the abundance of protons below the ankle. We show that astrophysical properties of the UHECR source environment such as the temperature, size, and magnetic field can be constrained by UHECR and neutrino data. Applying this to candidate sources with a simple structure shows that starburst galaxies are consistent with these constraints, but galaxy clusters may be in tension with them. For multi-component systems like AGNs and GRBs the results are indicative but customized analysis is needed for definitive conclusions.

Cosmology may give rise to appreciable populations of both particle dark matter and primordial black holes (PBH) with the combined mass density providing the observationally inferred value $\Omega_{\rm DM}\approx0.26$. However, previous studies have highlighted that scenarios with both particle dark matter and PBH are strongly excluded by $\gamma$-ray limits for particle dark matter with a velocity independent thermal cross section $\langle\sigma v\rangle\sim3\times10^{-26}{\rm cm}^3/{\rm s}$, as is the case for classic WIMP dark matter. Here we extend these existing studies on $s$-wave annihilating particle dark matter to ascertain the limits from diffuse $\gamma$-rays on velocity dependent annihilations which are $p$-wave with $\langle\sigma v \rangle\propto v^2$ or $d$-wave with $\langle\sigma v \rangle\propto v^4$, which we find to be considerably less constraining. Furthermore, we highlight that even if the freeze-out process is $p$-wave it is relatively common for (loop/phase-space) suppressed $s$-wave processes to actually provide the leading contributions to the experimentally constrained $\gamma$-ray flux from the PBH halo. This work also utilyses a refined treatment of the PBH dark matter density profile and outlines an improved application of extra-galactic $\gamma$-ray bounds.

Weakly Interacting Massive Particles (WIMPs) are among the best-motivated dark matter candidates. In the standard scenario where the freeze-out happens well after the end of inflationary reheating, they are in tension with the severe experimental constraints. Here, we investigate the thermal freeze-out of WIMPs occurring {\it during} reheating, while the inflaton $\phi$ coherently oscillates in a generic potential $\propto \phi^n$. Depending on the value of $n$ and the spin of the inflaton decaying products, the evolution of the radiation and inflaton energy densities can show distinct features, therefore, having a considerable impact on the freeze-out behavior of WIMPs. As a result of the injection of entropy during reheating, the parameter space compatible with the observed DM relic abundance is enlarged. In particular, the WIMP thermally averaged annihilation cross-section can be several magnitudes lower than that in the standard case. Finally, we discuss the current bounds from dark matter indirect detection experiments, and explore future challenges and opportunities.

In the context of the most general scalar-vector-tensor theory, we study the stability of static spherically symmetric black holes under linear odd-parity perturbations. We calculate the action to second order in the linear perturbations to derive a master equation for these perturbations. For this general class of models, we obtain the conditions of no-ghost and Laplacian instability. Then, we study in detail the generalized Regge-Wheeler potential of particular cases to find their stability conditions.

We study the non-equilibrium dynamics of a pseudoscalar axion-like particle (ALP) weakly coupled to degrees of freedom in thermal equilibrium by obtaining its reduced density matrix. Its time evolution is determined by the in-in effective action which we obtain to leading order in the (ALP) coupling but to \emph{all orders} in the couplings of the bath to other fields within or beyond the standard model. The effective equation of motion for the (ALP) is a Langevin equation with noise and friction kernels obeying the fluctuation dissipation relation. A ``misaligned'' initial condition yields damped coherent oscillations, however, the (ALP) population increases towards thermalization with the bath. As a result, the energy density features a mixture of a cold component from misalignment and a hot component from thermalization with proportions that vary in time $(cold)\,e^{-\Gamma t}+(hot)\,(1-e^{-\Gamma t})$, providing a scenario wherein the ``warmth'' of the dark matter evolves in time from colder to hotter. As a specific example we consider the (ALP)-photon coupling $g a \vec{E}\cdot \vec{B}$ to lowest order, valid from recombination onwards. For $T \gg m_a$ the long-wavelength relaxation rate is substantially enhanced $\Gamma_T = \frac{g^2\,m^2_a\,T}{16\pi} $. The ultraviolet divergences of the (ALP) self-energy require higher order derivative terms in the effective action. We find that at high temperature, the finite temperature effective mass of the (ALP) is $m^2_a(T) = m^2_a(0)\Big[ 1-(T/T_c)^4\Big]$, with $T_c \propto \sqrt{m_a(0)/g}$, \emph{suggesting} the possibility of an inverted phase transition, which when combined with higher derivatives may possibly indicate exotic new phases. We discuss possible cosmological consequences on structure formation and the effective number of relativistic species.

We show that dark-matter candidates with large masses and large nuclear interaction cross sections are detectable with terrestrial radar systems. We develop our results in close comparison to successful radar searches for tiny meteoroids, aggregates of ordinary matter. The path of a meteoroid (or suitable dark-matter particle) through the atmosphere produces ionization deposits that reflect incident radio waves. We calculate the equivalent radar echoing area or `radar cross section' for dark matter. By comparing the expected number of dark-matter-induced echoes with observations, we set new limits in the plane of dark-matter mass and cross section, complementary to pre-existing cosmological limits. Our results are valuable because (A) they open a new detection technique for which the reach can be greatly improved and (B) in case of a detection, the radar technique provides differential sensitivity to the mass and cross section, unlike cosmological probes.

We use the quantum unimodular theory of gravity to relate the value of the cosmological constant, $\Lambda$, and the energy scale for the emergence of cosmological classicality. The fact that $\Lambda$ and unimodular time are complementary quantum variables implies a perennially quantum Universe should $\Lambda$ be zero (or, indeed, fixed at any value). Likewise, the smallness of $\Lambda$ puts an upper bound on its uncertainty, and so a lower bound on the unimodular clock's uncertainty or the cosmic time for the emergence of classicality. Far from being the Planck scale, classicality arises at around $7 \times 10^{11}$ GeV for the observed $\Lambda$, and taking the region of classicality to be our Hubble volume. We confirm this argument with a direct evaluation of the wavefunction of the Universe in the connection representation for unimodular theory. Our argument is robust, with the only leeway being in the comoving volume of our cosmological classical patch, which should be bigger than that of the observed last scattering surface. Should it be taken to be the whole of a closed Universe, then the constraint depends weakly on $\Omega_k$: for $-\Omega_k < 10^{-3}$ classicality is reached at $ > 4 \times 10^{12}$ GeV. If it is infinite, then this energy scale is infinite, and the Universe is always classical within the minisuperspace approximation. It is a remarkable coincidence that the only way to render the Universe classical just below the Planck scale is to define the size of the classical patch as the scale of non-linearity for a red spectrum with the observed spectral index $n_s = 0.967(4)$ (about $10^{11}$ times the size of the current Hubble volume). In the context holographic cosmology, we may interpret this size as the scale of confinement in the dual 3D quantum field theory, which may be probed (directly or indirectly) with future cosmological surveys.

This is the report for the topical group Theory of Neutrino Physics (TF11/NF08) for Snowmass 2021. This report summarizes the progress in the field of theoretical neutrino physics in the past decade, the current status of the field, and the prospects for the upcoming decade.

We study theoretically mode conversion and resonant overreflection of magnetohydrodynamic waves in an inhomogeneous plane-stratified plasma in the presence of a nonuniform shear flow, using precise numerical calculations of the reflection and transmission coefficients and the field distributions based on the invariant imbedding method. The cases where the flow velocity and the external magnetic field are directed perpendicularly to the inhomogeneity direction and both the flow velocity and the plasma density vary arbitrarily along it are considered. When there is a shear flow, the wave frequency is modulated locally by the Doppler shift and resonant amplification and overreflection occur where the modulated frequency is negative and its absolute value matches the local Alfv\'en or slow frequency. For many different types of the density and flow velocity profiles, we find that, especially when the parameters are such that the incident waves are totally reflected, there arises a giant overreflection where the reflectance is much larger than 10 in a fairly broad range of the incident angles, the frequency, and the plasma $\beta$ and its maximum attains values larger than $10^5$. In a finite $\beta$ plasma, both incident fast and slow magnetosonic waves are found to cause strong overreflection and there appear multiple positions exhibiting both Alfv\'en and slow resonances inside the plasma. We explain the mechanism of overreflection in terms of the formation of inhomogeneous and open cavities close to the resonances and the strong enhancement of the wave energy due to the occurrence of semi-bound states there. We give discussions of the observational consequences in magnetized terrestrial and solar plasmas.