Papers for Thursday, Aug 12 2021

The LIGO Scientific Collaboration and Virgo Collaboration (LVC) have recently reported in GWTC-2.1 eight additional high-significance candidates with a probability of astrophysical origin greater than 0.5 in the LVC search on O3a running. In GWTC-2.1, the majority of the effective inspiral spins ($\chi_{\rm eff}$) show magnitudes consistent with zero, while two (GW190403$_{-}$051519 and GW190805$_{-}$211137) of the eight new events have $\chi_{\rm eff}$ $> 0$(at 90\% credibility). We note that GW190403$_{-}$051519 was reported with $\chi_{\rm eff}$ = $0.70^{+0.15}_{-0.27}$ and mass ratio $q$ = $0.25^{+0.54}_{-0.11}$, respectively. For binaries with very unequal masses, a stronger statements of spins can be made. In the case of GW190403$_{-}$051519, its primary dimensionless spin, $\chi_1$ = $0.92^{+0.07}_{-0.22}$, marks the fastest first-born black hole ever measured in all current gravitational-wave events. If we believe that GW190403$_{-}$051519 formed through isolated binary evolution channel, this extremely high spin thus challenges the existence of efficient angular momentum transport mechanism, i.e., Spruit-Tayler (ST) dynamo \citep{2002A&A...381..923S}, between the stellar core and its radiative envelopes of massive sta

Recently, several experiments reported a muon deficit in air shower simulations with respect to the data. This problem can be studied using an estimator that quantifies the relative muon content of the data with respect to those of proton and iron Monte Carlo air shower simulations. We analyze two estimators. The first one, based on the logarithm of the mean of the muon content, is built from experimental considerations. It is ideal for comparing results from different experiments as it is independent of the detector resolution. The second estimator is based on the mean of the logarithm of the muon content, which implies that it depends on shower-to-shower fluctuations. It is linked to the mean-logarithmic mass $\left \langle \ln A \right \rangle$ through the Heitler-Matthews model. We study the properties of the estimators and their biases considering the knowns and unknowns of typical experiments. Furthermore, we study these effects in measurements of the muon density at $1000\,$m from the shower axis obtained by the Akeno Giant Air Shower Array (AGASA). Finally, we report the estimates of the relative muon content of the AGASA data, which support a muon deficit in simulations. These estimates constitute valuable additional information of the muon content of air showers at the highest energies.

We present a field-based approach to the analysis of cosmic shear data to infer jointly cosmological parameters and the dark matter distribution. This forward modelling approach samples the cosmological parameters and the initial matter fluctuations, using a physical gravity model to link the primordial fluctuations to the non-linear matter distribution. Cosmological parameters are sampled and updated consistently through the forward model, varying (1) the initial matter power spectrum, (2) the geometry through the distance-redshift relationship, and (3) the growth of structure and light-cone effects. Our approach extracts more information from the data than methods based on two-point statistics. We find that this field-based approach lifts the strong degeneracy between the cosmological matter density, $\Omega_\mathrm{m}$, and the fluctuation amplitude, $\sigma_8$, providing tight constraints on these parameters from weak lensing data alone. In the simulated four-bin tomographic experiment we consider, the field-based likelihood yields marginal uncertainties on $\sigma_8$ and $\Omega_\mathrm{m}$ that are, respectively, a factor of 3 and 5 smaller than those from a two-point power spectrum analysis applied to the same underlying data.

In this work we revisit five different point sources within or behind galaxy clusters in order to constrain the coupling constant between axion-like particles (ALPs) and photons. We use three distinct machine learning (ML) techniques and compare our results with a standard $\chi^2$ analysis. For the first time we apply approximate Bayesian computation to searches for ALPs and find consistently good performance across ML classifiers. Further, we apply more realistic 3D magnetic field simulations of galaxy clusters and compare our results with previously used 1D simulations. We find constraints on the ALP-photon coupling at the level of state-of-the-art bounds with $g_{a\gamma\gamma} \lesssim 0.6 \times 10^{-12}$ GeV${}^{-1}$, hence improving on previous constraints obtained from the same observations.

Multiple experiments reported evidences of a muon deficit in air shower simulations with respect to data, which increases with the primary energy. In this work, we study the muon deficit using measurements of the muon density at $1000\,$m from the shower axis obtained by the Akeno Giant Air Shower Array (AGASA). The selected events have reconstructed energies in the range $18.83\,\leq\,\log_{10}(E_{R}/\textrm{eV})\,\leq\,19.46$ and zenith angles $\theta\leq 36^{\circ}$. We compare these muon density measurements to proton, iron, and mixed composition scenarios, obtained by using the high-energy hadronic interaction models EPOS-LHC, QGSJetII-04, and Sibyll2.3c. We find that AGASA data are compatible with a heavier composition, lying above the predictions of the mixed composition scenarios. The average muon density divided by the energy in AGASA data is greater than in the mixed composition scenarios by a factor of $1.49\pm0.11\,\textrm{(stat)}\pm^{0.49}_{0.30}\,\textrm{(syst)}$, $1.54\pm0.12\,\textrm{(stat)}\pm^{0.50}_{0.31}\,\textrm{(syst)}$, and $1.66\pm0.13\,\textrm{(stat)} \pm ^{0.54}_{0.34}\,\textrm{(syst)}$ for EPOS-LHC, Sibyll2.3c, and QGSJetII-04, respectively. We interpret this as further evidence of a muon deficit in air shower simulations at the highest energies.

SCExAO at the Subaru telescope is a visible and near-infrared high-contrast imaging instrument employing extreme adaptive optics and coronagraphy. The instrument feeds the near-infrared light (JHK) to the integral-field spectrograph CHARIS. The spectropolarimetric capability of CHARIS is enabled by a Wollaston prism and is unique among high-contrast imagers. We present a detailed Mueller matrix model describing the instrumental polarization effects of the complete optical path, thus the telescope and instrument. From measurements with the internal light source, we find that the image derotator (K-mirror) produces strongly wavelength-dependent crosstalk, in the worst case converting ~95% of the incident linear polarization to circularly polarized light that cannot be measured. Observations of an unpolarized star show that the magnitude of the instrumental polarization of the telescope varies with wavelength between 0.5% and 1%, and that its angle is exactly equal to the altitude angle of the telescope. Using physical models of the fold mirror of the telescope, the half-wave plate, and the derotator, we simultaneously fit the instrumental polarization effects in the 22 wavelength bins. Over the full wavelength range, our model currently reaches a total polarimetric accuracy between 0.08% and 0.24% in the degree of linear polarization. We propose additional calibration measurements to improve the polarimetric accuracy to <0.1% and plan to integrate the complete Mueller matrix model into the existing CHARIS post-processing pipeline. Our calibrations of CHARIS' spectropolarimetric mode will enable unique quantitative polarimetric studies of circumstellar disks and planetary and brown dwarf companions.

Becker et al. 2021 measured the mean free path of Lyman limit photons in the IGM at $z=6$. The short value suggests that absorptions may have played a prominent role in reionization. Here we study physical properties of ionizing photon sinks in the wake of ionization fronts (I-fronts) using radiative hydrodynamic simulations. We quantify the contributions of gaseous structures to the Lyman limit opacity by tracking the column density distributions in our simulations. Within $\Delta t = 10$ Myr of I-front passage, we find that self-shielding systems ($N_{\rm HI} > 10^{17.2}$ cm$^{-2}$) are comprised of two distinct populations: (1) over-density $\Delta \sim 50$ structures in photo-ionization equilibrium with the ionizing background; (2) $\Delta \gtrsim 100$ density peaks with fully neutral cores. The self-shielding systems contribute more than half of the opacity at these times, but the IGM evolves considerably in $\Delta t \sim 100$ Myr as structures are flattened by pressure smoothing and photoevaporation. By $\Delta t = 300$ Myr, they contribute $\lesssim 10 \%$ to the opacity in an average 1 Mpc$^3$ patch of the Universe. The percentage can be a factor of a few larger in over-dense patches, where more self-shielding systems survive. We quantify the characteristic masses and sizes of self-shielding structures. Shortly after I-front passage, we find $M=10^{4} - 10^8$ M$_\odot$ and effective diameters $d_{\rm eff} = 1 - 20$ ckpc$/h$. These scales increase as the gas relaxes. The picture herein presented may be different in dark matter models with suppressed small-scale power.

Understanding the inner structure of the clumpy molecular torus surrounding the active galactic nucleus is essential in revealing the forming mechanism. However, spatially resolving the torus is difficult because of its size of a few parsecs. Thus, to probe the clump conditions in the torus, we performed the velocity decomposition of the CO ro-vibrational absorption lines ($\Delta{v}=0\to 1,\ \Delta{J}=\pm 1$) at $\lambda\sim 4.67\ \mathrm{\mu{m}}$ observed toward an ultra-luminous infrared galaxy IRAS 08572+3915 NW with the high-resolution spectroscopy ($R\sim 10000$) of Subaru Telescope. Consequently, we found that each transition had two outflowing components, i.e., (a) and (b), both at approximately $\sim -160\ \mathrm{km\ s^{-1}}$, but with broad and narrow widths, and an inflowing component, i.e., (c), at approximately $\sim +100\ \mathrm{km\ s^{-1}}$, which were attributed to the torus. The ratios of the velocity dispersions of each component lead to those of the rotating radii around the black hole of $R_\mathrm{rot,a}:R_\mathrm{rot,b}:R_\mathrm{rot,c}\approx 1:5:17$, indicating the torus where clumps are outflowing in the inner regions and inflowing in the outer regions if a hydrostatic disk with $\sigma_V\propto R_\mathrm{rot}^{-0.5}$ is assumed. Based on the kinetic temperature of components (a) and (b) of $\sim 720\ \mathrm{K}$ and $\sim 25\ \mathrm{K}$ estimated from the level population, the temperature gradient is $T_\mathrm{kin}\propto R_\mathrm{rot}^{-2.1}$. Magnetohydrodynamic models with large density fluctuations of two orders of magnitude or more are necessary to reproduce this gradient.

We report on the search for electromagnetic counterparts to the nine gravitational-wave events with a $>$60\% probability of containing a neutron star during the third (O3) LIGO-Virgo Collaboration (LVC) observing run with the All-Sky Automated Survey for SuperNovae (ASAS-SN). No optical counterparts associated with a gravitational wave event was found. However, thanks to its network of telescopes, the average area visible to at least one ASAS-SN site during the first 10 hours after the trigger contained $\sim$30\% of the integrated source location probability. Through a combination of normal operations and target-of-opportunity observations, ASAS-SN observations of the highest probability fields began within one hour of the trigger for four of the events. After 24 hours, ASAS-SN observed $>$60\% of total probability for three events and $>$40\% for all but one of the events. This is the largest area coverage to a depth of $g = 18.5$ mag from any survey with published coverage statistics for seven of the nine events. With its observing strategy, five sites around the world, and a large field of view, ASAS-SN will be one of the leading surveys to optically search for nearby neutron star mergers during LVC O4.

We estimate the detectability of X-ray metal-line emission from the circumgalactic medium (CGM) of galaxies over a large halo mass range ($\mathrm{M}_{\mathrm{200c}} =10^{11.5}$-$10^{14.5}\,\mathrm{M}_{\odot}$) using the EAGLE simulations. With the XRISM Resolve instrument, a few bright (K-$\alpha$ or Fe L-shell) lines from $\mathrm{M}_{\mathrm{200c}} \gtrsim 10^{13}\,\mathrm{M}_{\odot}$ haloes should be detectable. Using the Athena X-IFU or the Lynx Main Array, emission lines (especially from O$\,$VII and O$\,$VIII) from the inner CGM of $\mathrm{M}_{\mathrm{200c}} \gtrsim10^{12.5}\,\mathrm{M}_{\odot}$ haloes become detectable, and intragroup and intracluster gas will be detectable out to the virial radius. With the Lynx Ultra-high Resolution Array, the inner CGM of haloes hosting $\mathrm{L}_{*}$ galaxies is accessible. These estimates do assume long exposure times ($\sim 1\,$Ms) and large spatial bins ($\sim1$-$10\,\mathrm{arcmin}^{2}$). We also investigate the properties of the gas producing this emission. CGM emission is dominated by collisionally ionized (CI) gas, and tends to come from halo centres. The gas is typically close to the maximum emissivity temperature for CI gas ($\mathrm{T}_\mathrm{peak}$), and denser and more metal-rich than the bulk of the CGM at a given distance from the central galaxy. However, for the K-$\alpha$ lines, emission can come from hotter gas in haloes where the virialized, volume-filling gas is hotter than $\mathrm{T}_\mathrm{peak}$. Trends of emission with halo mass can largely be explained by differences in virial temperature. Differences between lines generally result from the different behaviour of the emissivity as a function of temperature of the K-$\alpha$, He-$\alpha$-like, and Fe~L-shell lines. We conclude that upcoming X-ray missions will open up a new window onto the hot CGM.

A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a new source that can generate outward migration of low-mass planets. While the Lindblad and corotation torques depend mostly on the properties of the protoplanetary disc and on the planet mass, the thermal torque depends also on the luminosity of the planet, arising mainly from the accretion of solids. Thus, the accretion of solids plays an important role not only in the formation of the planet but also in its migration process. In a previous work, we evaluated the thermal torque effects on planetary growth and migration mainly in the planetesimal accretion paradigm. In this new work, we study the role of the thermal torque within the pebble accretion paradigm. Computations are carried out consistently in the framework of a global model of planet formation that includes disc evolution, dust growth and evolution, and pebble formation. We also incorporate updated prescriptions of the thermal torque derived from high resolution hydrodynamical simulations. Our simulations show that the thermal torque generates extended regions of outward migration in low viscosity discs. This has a significant impact in the formation of the planets.

Aims: To provide the first demonstration of STIX Fourier-transform X-ray imaging using semi-calibrated (amplitude-only) visibility data acquired during the Solar Orbiter's cruise phase. Methods: We use a parametric imaging approach by which STIX visibility amplitudes are fitted by means of two non-linear optimization methods: a fast meta-heuristic technique inspired by social behavior, and a Bayesian Monte Carlo sampling method, which, although slower, provides better quantification of uncertainties. Results: When applied to a set of solar flare visibility amplitudes recorded by STIX on November 18, 2020 the two parametric methods provide very coherent results. The analysis also demonstrates the ability of STIX to reconstruct high time resolution information and, from a spectral viewpoint, shows the reliability of a double-source scenario consistent with a thermal versus nonthermal interpretation. Conclusions: In this preliminary analysis of STIX imaging based only on visibility amplitudes, we formulate the imaging problem as a non-linear parametric issue we addressed by means of two high-performance optimization techniques that both showed the ability to sample the parametric space in an effective fashion, thus avoiding local minima.

In the last two decades, very-high-energy gamma-ray astronomy has reached maturity: over 200 sources have been detected, both Galactic and extragalactic, by ground-based experiments. At present, Active Galactic Nuclei (AGN) make up about 40% of the more than 200 sources detected at very high energies with ground-based telescopes, the majority of which are blazars, i.e. their jets are closely aligned with the line of sight to Earth and three quarters of which are classified as high-frequency peaked BL Lac objects. One challenge to studies of the cosmological evolution of BL Lacs is the difficulty of obtaining redshifts from their nearly featureless, continuum- dominated spectra. It is expected that a significant fraction of the AGN to be detected with the future Cherenkov Telescope Array (CTA) observatory will have no spectroscopic redshifts, compromising the reliability of BL Lac population studies, particularly of their cosmic evolution. We started an effort in 2019 to measure the redshifts of a large fraction of the AGN that are likely to be detected with CTA, using the Southern African Large Telescope (SALT). In this contribution, we present two results from an on-going SALT program focused on the determination of BL Lac object redshifts that will be relevant for the CTA observatory.

The Astrophysical Multimessenger Observatory Network (AMON), has developed a real-time multi-messenger alert system. The system performs coincidence analyses of datasets from gamma-ray and neutrino detectors, making the Neutrino-Electromagnetic (NuEM) alert channel. For these analyses, AMON takes advantage of sub-threshold events, i.e., events that by themselves are not significant in the individual detectors. The main purpose of this channel is to search for gamma-ray counterparts of neutrino events. We will describe the different analyses that make up this channel and present a selection of recent results.

We use the observed cumulative statistics of CIV absorbers and dark matter halos to infer the distribution of CIV-absorbing gas relative to galaxies at redshifts $0\!\leq\!z\!\leq\!5$. We compare the cosmic incidence $dN/dX$ of CIV absorber populations and galaxy halos, finding that massive $L \geq L_{\star}$ halos alone cannot account for all the observed $W_r \geq 0.05$~{\AA} absorbers. However, the $dN/dX$ of lower mass halos exceeds that of $W_r \geq 0.05$~{\AA} absorbers. We also estimate the characteristic gas radius of absorbing structures required for the observed CIV $dN/dX$, assuming each absorber is associated with a single galaxy halo. The $W_r \geq 0.3$~{\AA} and $W_r \geq 0.6$~{\AA} CIV gas radii are $\sim30-70\%$ ($\sim20-40\%$) of the virial radius of $L_{\star}$ ($0.1L_{\star}$) galaxies, and the $W_r \geq 0.05$~{\AA} gas radius is $\sim100-150\%$($\sim60-100\%$) of the virial radius of $L_{\star}$ ($0.1L_{\star}$) galaxies. For stronger absorbers, the gas radius relative to virial radius rises across Cosmic Noon and falls afterwards, while for weaker absorbers, the relative gas radius declines across Cosmic Noon and then dramatically rises at $z\!<\!1$. A strong luminosity-dependence of gas radius implies highly extended CIV envelopes around massive galaxies before Cosmic Noon, while a luminosity-independent gas radius implies highly extended envelopes around dwarf galaxies after Cosmic Noon. From available absorber-galaxy and CIV evolution data, we favor a scenario in which low-mass galaxies enrich the volume around massive galaxies at early epochs and propose that the outer halo gas ($>0.5R_v$) was produced primarily in ancient satellite dwarf galaxy outflows, while the inner halo gas ($<0.5R_v$) originated from the central galaxy and persists as recycled accreting gas.

We present polynomial conversions between each of 11 different strong line gas-phase metallicity calibrations, each based on $\sim$ 1.1 million star-forming spaxels in the public Sloan Digital Sky Survey (SDSS) Data Release 15 (DR15) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. For this sample, which is $\sim$ 20 times larger than previous works, we present 5th order polynomial fits for each of 110 possible calibration conversions, for both Small Magellanic Cloud (SMC)-type and Milky Way (MW)-type dust corrections. The typical $2\sigma$ scatter around our polynomial fits is 0.1 dex; we present the range over which the metallicities are valid. Conversions between metallicities which rely on the same set of line ratios, or a heavily shared set of emission lines, have reduced scatter in their conversions relative to those conversions with little overlap in required emission lines. Calibration conversions with less consistent sets of emission lines also have increased galaxy-to-galaxy variability, and this variability can account for up to 35% of the total scatter. We also compare our conversions to previous work with the single fibre SDSS DR7 spectra along with higher spatial resolution data from the TYPHOON Integral Field Spectroscopy survey, resulting in comparison samples with spatial resolutions from several kpc down to $\sim$100 pc. Our metallicity conversions, obtained with the large sample of MaNGA, are robust against the influence of diffuse ionized gas, redshift, effective radius and spatial blurring, and are therefore consistent across both integrated spectra and the high resolution integral field spectroscopy data.

The STRings for Absorption length in Water (STRAW) experiment is the first in a series of pathfinder missions for the Pacific Ocean Neutrino Experiment (P-ONE), a future large-scale neutrino telescope in the north-eastern Pacific Ocean. STRAW consists of two 150 m long mooring lines instrumented with optical emitters and detectors. The experiment is designed to measure the attenuation length of the water and perform a long-term assessment of the optical background at the future P-ONE site. After two years of continuous operation, measurements from STRAW show an optical attenuation length of about 28 metres at 450 nm. Additionally, the data allows a study of the ambient undersea background. The overall optical environment reported here is comparable to other deep-water neutrino telescopes and qualifies the site for the deployment of P-ONE.

The IceCube Neutrino Observatory detects neutrinos at energies orders of magnitude higher than those available to current accelerators. Above 40 TeV, neutrinos traveling through the Earth will be absorbed as they interact via charged current interactions with nuclei, creating a deficit of Earth-crossing neutrinos detected at IceCube. The previous published results showed the cross section to be consistent with Standard Model predictions for 1 year of IceCube data. We present a new analysis that uses 8 years of IceCube data to fit the $\nu_\mu$ absorption in the Earth, with statistics an order of magnitude better than previous analyses, and with an improved treatment of systematic uncertainties. It will measure the cross section in three energy bins that span the range 1 TeV to 100 PeV. We will present Monte Carlo studies that demonstrate its sensitivity.

Magnetic fields are dynamically important in the diffuse interstellar medium. Understanding how gravitationally bound, star-forming clouds form requires modeling of the fields in a self-consistent, supernova-driven, turbulent, magnetized, stratified disk. We employ the FLASH magnetohydrodynamics code to follow the formation and early evolution of clouds with final masses of 3-8 $\times 10^3 M_{\odot}$ within such a simulation. We use the code's adaptive mesh refinement capabilities to concentrate numerical resolution in zoom-in regions covering single clouds, allowing us to investigate the detailed dynamics and field structure of individual self-gravitating clouds in a consistent background medium. Our goal is to test the hypothesis that dense clouds are dynamically evolving objects far from magnetohydrostatic equilibrium. We find that the cloud envelopes are magnetically supported with field lines parallel to density gradients and flow velocity, as indicated by the histogram of relative orientations and other statistical measures. In contrast, the dense cores of the clouds are gravitationally dominated, with gravitational energy exceeding internal, kinetic, or magnetic energy and accelerations due to gravity exceeding those due to magnetic or thermal pressure gradients. In these regions field directions vary strongly, with a slight preference towards being perpendicular to density gradients, as shown by three-dimensional histograms of relative orientation.

Peptide bonds, as the molecular bridges that connect amino acids, are crucial to the formation of proteins. Searches and studies of molecules with embedded peptide-like bonds are thus important for the understanding of protein formation in space. Here we report the first tentative detection of propionamide (C2H5CONH2), the largest peptide-like molecule detected in space toward Sagittarius B2(N1) at a position called N1E that is slightly offset from the continuum peak. A new laboratory measurements of the propionamide spectrum were carried out in the 9-461 GHz, which provide good opportunity to check directly for the transition frequencies of detected interstellar lines of propionamide. Our observing result indicates that propionamide emission comes from the warm, compact cores in Sagittarius B2, in which massive protostellars are forming. The column density of propionamide toward Sgr B2(N1E) was derived to be 1.5\times 10^{16} cm^-2, which is three fifths of that of acetamide, and one nineteenth of that of formamide. This detection suggests that large peptide-like molecules can form and survive during star-forming process and may form more complex molecules in the interstellar medium. The detection of propionamide bodes well for the presence of polypeptides, as well as other complex prebiotic molecules in the interstellar medium.

The observed large-scale alignment of polarization angles and galaxy axis have been challenging the fundamental assumption of homogeneity and isotropy in standard cosmology since more than two decades. The intergalactic magnetic field, and its correlations in real space, potentially seems as a viable candidate for explaining this phenomenon. It has been shown earlier that the large-scale intergalactic magnetic field correlations can explain the alignment signal of quasars over Gpc scale, interestingly they can also explain the radio polarization alignment observed in JVAS/CLASS data over 100 Mpc. Motivated with recent observations of galaxy axis alignment over several tens of Mpc, and Mpc scale, i.e., the cluster scale, we further explore the correlations of background magnetic field to explain these relatively small scale alignment observations. In particular, we explore two recently claimed signals of alignment in the radio sources in the FIRST catalog and in the ACO clusters. We find that both of these can be explained in terms of the intergalactic magnetic field with a spectral index of $-2.62\pm 0.03$. The large-scale magnetic field correlations potentially seem to explain the polarization and galaxy axis alignment from Gpc to Mpc scales.

We investigate the stellar population properties of bulges and disks separately for 34 S0s using integral field spectroscopy from the Calar Alto Legacy Integral Field Area survey. The spatially resolved stellar age and metallicity of bulge and disk components have been simultaneously estimated using the penalized pixel fitting method with photometrically defined weights for the two components. We find a tight correlation between age and metallicity for bulges, while the relation for disks has a larger scatter than that of bulges. This implies that the star formation histories of the disks are more complicated than those of the bulges. Bulges of the high-mass S0s are mostly comparable in metallicity, while bulges appear to be systematically more metal-rich than disks for the low-mass S0s. The ages of bulges and disks in the high-mass S0s appear to increase with local density. The bulge ages of the low-mass S0s also increases with local density, but such a trend is not clear in the disk ages of low-mass S0s. In addition, the age difference between bulge and disk components (delta Age) tends to increase with local density, both for the high-mass and low-mass S0s. The high-mass S0s have systematically higher delta Age than the low-mass S0s at given local density. Our results indicate that the stellar mass significantly influences the evolution of S0 galaxies, but the environment also plays an important role in determining the evolution of bulges and disks at given stellar mass.

We present the first determination of the abundance ratios of $^{13}$C substitutions of cyanoacetylene (HC$_{3}$N), [H$^{13}$CCCN]:[HC$^{13}$CCN]:[HCC$^{13}$CN] in Titan's atmosphere measured using millimeter-wave spectra obtained by the Atacama Large Millimeter-submillimeter Array. To compare the line intensities precisely, datasets which include multiple molecular lines were extracted to suppress effects of Titan's environmental conditions and observation settings. The [HC$^{13}$CCN]:[HCC$^{13}$CN] and [H$^{13}$CCCN]:[HCC$^{13}$CN}] ratios were obtained from 12 and 1 selected datasets, respectively. As a result, nearly the uniform [H$^{13}$CCCN]:[HC$^{13}$CCN]:[HCC$^{13}$CN] abundance ratios as 1.17 ($\pm$0.20) : 1.09 ($\pm$0.25) : 1 (1$\sigma$) were derived, whereas previously reported ratios for interstellar medium (ISM) have shown large anomalies that may be caused by $^{13}$C concentrations in precursors. The result obtained here suggests that $^{13}$C concentration processes suggested in the ISM studies do not work effectively on precursors of HC$_{3}$N and HC$_{3}$N itself due to Titan's high atmospheric temperature and/or depletion of both $^{13}$C and $^{13}$C$^+$.

The Cherenkov Telescope Array (CTA) is the next-generation gamma-ray observatory that is expected to reach one order of magnitude better sensitivity than that of current telescope arrays. The Large-Sized Telescopes (LSTs) have an essential role in extending the energy range down to 20 GeV. The prototype LST (LST-1) proposed for CTA was built in La Palma, the northern site of CTA, in 2018. LST-1 is currently in its commissioning phase and moving towards scientific observations. The LST-1 camera consists of 1855 photomultiplier tubes (PMTs) which are sensitive to Cherenkov light. PMT signals are recorded as waveforms sampled at 1 GHz rate with Domino Ring Sampler version 4 (DRS4) chips. Fast sampling is essential to achieve a low energy threshold by minimizing the integration of background light from the night sky. Absolute charge calibration can be performed by the so-called F-factor method, which allows calibration constants to be monitored even during observations. A calibration pipeline of the camera readout has been developed as part of the LST analysis chain. The pipeline performs DRS4 pedestal and timing corrections, as well as the extraction and calibration of charge and time of pulses for subsequent higher-level analysis. The performance of each calibration step is examined, and especially charge and time resolution of the camera readout are evaluated and compared to CTA requirements. We report on the current status of the calibration pipeline, including the performance of each step through to signal reconstruction, and the consistency with Monte Carlo simulations.

Globally-propagating shocks in the solar corona have long been studied to quantify their involvement in the acceleration of energetic particles. However, this work has tended to focus on large events associated with strong solar flares and fast coronal mass ejections (CMEs), where the waves are sufficiently fast to easily accelerate particles to high energies. Here we present observations of particle acceleration associated with a global wave event which occurred on 1 October 2011. Using differential emission measure analysis, the global shock wave was found to be incredibly weak, with an Alfv\'en Mach number of ~1.008-1.013. Despite this, spatially-resolved type III radio emission was observed by the Nan\c{c}ay RadioHeliograph at distinct locations near the shock front, suggesting localised acceleration of energetic electrons. Further investigation using a magnetic field extrapolation identified a fan structure beneath a magnetic null located above the source active region, with the erupting CME contained within this topological feature. We propose that a reconfiguration of the coronal magnetic field driven by the erupting CME enabled the weak shock to accelerate particles along field lines initially contained within the fan and subsequently opened into the heliosphere, producing the observed type III emission. These results suggest that even weak global shocks in the solar corona can accelerate energetic particles via reconfiguration of the surrounding magnetic field.

We present the spectroscopic and photometric study of five intermediate-luminosity red transients (ILRTs), namely AT 2010dn, AT 2012jc, AT 2013la, AT 2013lb, and AT 2018aes. They share common observational properties and belong to a family of objects similar to the prototypical ILRT SN~2008S. These events have a rise time that is less than 15 days and absolute peak magnitudes of between $-11.5$ and $-14.5$ mag. Their pseudo-bolometric light curves peak in the range $0.5$ - $9.0 \times10^{40}~\mathrm{erg~s}^{-1}$ and their total radiated energies are on the order of (0.3 - 3) $\times$~10$^{47}$~erg. After maximum brightness, the light curves show a monotonic decline or a plateau, resembling those of faint supernovae IIL or IIP, respectively. At late phases, the light curves flatten, roughly following the slope of the $^{56}$Co decay. If the late-time power source is indeed radioactive decay, these transients produce $^{56}$Ni masses on the order of $10^{-4}$ to $10^{-3}$~\msun. The spectral energy distribution of our ILRT sample, extending from the optical to the mid-infrared (MIR) domain, reveals a clear IR excess soon after explosion and non-negligible MIR emission at very late phases. The spectra show prominent H lines in emission with a typical velocity of a few hundred km~s$^{-1}$, along with Ca~II features. In particular, the [Ca~II] $\lambda$7291,7324 doublet is visible at all times, which is a characteristic feature for this family of transients. The identified progenitor of SN~2008S, which is luminous in archival Spitzer MIR images, suggests an intermediate-mass precursor star embedded in a dusty cocoon. We propose the explosion of a super-asymptotic giant branch star forming an electron-capture supernova as a plausible explanation for these events.

We present a general relativistic, ray-tracing radiative transfer code RAIKOU for multi-wavlength studies of spectra and images including the black hole shadows around Kerr black holes. Important radiative processes in hot plasmas around black holes, i.e., (cyclo-)synchrotron, bremsstrahlung emission/absorption and Compton/inverse-Compton scattering, are incorporated. The Maxwell-J\"uttner and single/broken power-law electron distribution functions are implemented to calculate the radiative transfer via both of the thermal and the nonthermal electrons. Two calculation algorithms are implemented for studies of both the images and broadband spectra. An observer-to-emitter algorithm, which inversely solve the radiative transfer equation from the observer screen to emitting plasmas, is suitable for efficient calculations of the images, e.g., the black hole shadows, and spectra without the Compton effects. On the other hand, an emitter-to-observer algorithm, by which photons are transported with a Monte-Carlo method including the effects of Compton/inverse-Compton scatterings, enables us to compute multi-wavelength spectra with their energy bands broadly ranging from radio to very-high-energy gamma-ray. The code is generally applicable to accretion flows around Kerr black holes with relativistic jets and winds/coronae with various mass accretion rate (i.e., radiatively inefficient accretion flows, super-Eddington accretion flows, and others). We demonstrate an application of the code to a radiatively innefficent accretion flow onto a supermassive black hole.

The IceCube Neutrino Observatory is highly sensitive to neutrino bursts of $\mathcal{O}$(10) MeV energy that are would be generated by core collapse supernovae in our Galaxy. It will resolve temporal structures in supernova light curves particularly well. In the light of an improved understanding of the ice properties and the detector response, the effective area and the corresponding uncertainties were newly determined with a Geant4-based Monte Carlo. Uncertainties due to cross sections and oscillation effects in the Earth were also investigated. This analysis has been extended by simulating a very large sample to determine the small coincidence probability between optical modules that bears information on the average neutrino energy. These simulation results were then used to interpret the data in time and frequency space. While the availability to record data for low energy neutrinos from supernovae is close to perfect (99.2$\%$ between 2013-2020), the analysis requires that the detector works faultlessly and artifacts do not mimic the signal in the 13 years of data taken so far. An effort has been made to keep the uptime after all analysis steps similarly high. The frequency space can be studied in a range between 1 Hz and 1/year to test the detector stability with high accuracy, to study the influence of cosmic rays, and to search for periodic phenomena that lead to sub-threshold increases in the count rates. Here we discuss the results of the simulations and the corresponding systematic limitations, the method to reconstruct the mean neutrino energy for a recorded supernova, as well as aspects of the analyses of continuously taken optical module rate data in the time and frequency domain.

The Cherenkov Telescope Array (CTA) will be the next generation gamma-ray observatory, which will consist of three kinds of telescopes of different sizes. Among those, the Large Size Telescope (LST) will be the most sensitive in the low energy range starting from 20 GeV. The prototype LST (LST-1) proposed for CTA was inaugurated in October 2018 in the northern hemisphere site, La Palma (Spain), and is currently in its commissioning phase. MAGIC is a system of two gamma-ray Cherenkov telescopes of the current generation, located approximately 100 m away from LST-1, that have been operating in stereoscopic mode since 2009. Since LST-1 and MAGIC can observe the same air shower events, we can compare the brightness of showers, estimated energies of gamma rays, and other parameters event by event, which can be used to cross-calibrate the telescopes. Ultimately, by performing combined analyses of the events triggering the three telescopes, we can reconstruct the shower geometry more accurately, leading to better energy and angular resolutions, and a better discrimination of the background showers initiated by cosmic rays. For that purpose, as part of the commissioning of LST-1, we performed joint observations of established gamma-ray sources with LST-1 and MAGIC. Also, we have developed Monte Carlo simulations for such joint observations and an analysis pipeline which finds event coincidence in the offline analysis based on their timestamps. In this work, we present the first detection of an astronomical source, the Crab Nebula, with combined observation of LST-1 and MAGIC. Moreover, we show results of the inter-telescope cross-calibration obtained using Crab Nebula data taken during joint observations with LST-1 and MAGIC.

In these proceedings we summarize the current status of the study of the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. Gamma-ray emission is expected in galaxy clusters both from interactions of cosmic rays (CR) with the intra-cluster medium, or as a product of annihilation or decay of dark matter (DM) particles in case they are weakly interactive massive particles (WIMPs). The observation of Perseus has been proposed as one of the CTA Key Science Projects. In this contribution, we focus on the DM-induced component of the flux. Our DM modelling includes the substructures we expect in the main halo which will boost the annihilation signal significantly. We adopt an ON/OFF observation strategy and simulate the expected gamma-ray signals. Finally we compute the expected CTA sensitivity using a likelihood maximization analysis including the most recent CTA instrument response functions. In absence of signal, we show that CTA will allow us to provide stringent and competitive constraints on TeV DM, especially for the case of DM decay.

Stellar systems are often formed through the collapse of dense molecular clouds which, in turn, return copious amounts of atomic and molecular material to the interstellar medium. An in-depth understanding of chemical evolution during this cyclic interaction between the stars and the interstellar medium is at the heart of astrochemistry. Systematic chemical composition changes as interstellar clouds evolve from the diffuse stage to dense, quiescent molecular clouds to star-forming regions and proto-planetary disks further enrich the molecular diversity leading to the evolution of ever more complex molecules. In particular, the icy mantles formed on interstellar dust grains and their irradiation are thought to be the origin of many of the observed molecules, including those that are deemed to be prebiotic; that is those molecules necessary for the origin of life. This review will discuss both observational (e.g., ALMA, SOFIA, Herschel) and laboratory investigations using millimeter, submillimeter, and terahertz and far-IR (THz/F-IR) spectroscopies and the role that they play in contributing to our understanding of the formation of prebiotic molecules. Mid-IR spectroscopy has typically been the primary tool used in laboratory studies. However, THz/F-IR spectroscopy offers an additional and complementary approach in that it provides the ability to investigate intermolecular interactions compared to the intramolecular modes available in the mid-IR. THz/F-IR spectroscopy is still somewhat under-utilized, but with the additional capability it brings, its popularity is likely to significantly increase in the near future. This review will discuss the strengths and limitations of such methods, and will also provide some suggestions on future research areas that should be pursued in the coming decade exploiting both space-borne and laboratory facilities.

Context: Blue horizontal-branch stars are very old objects that can be used as markers in studies of the Galactic structure and formation history. To create a clean sky catalogue of blue horizontal-branch stars, we cross-matched the Gaia data release 2 (DR2) dataset with existing reference catalogues to define selection criteria based on Gaia DR2 parameters. Following the publication of Gaia early data release 3 (EDR3), these methods were verified and subsequently applied to this latest release. Aims: The purpose of this catalogue is to identify a set of blue horizontal-branch star candidates that have been selected using photometric and astrometric observations and exhibits a low contamination rate. Methods: We cross-matched reference blue horizontal-branch datasets with the Gaia DR2 database and defined two sets of selection criteria. Firstly, in Gaia DR2 - colour and absolute G magnitude space, and secondly, in Gaia DR2 - colour and reduced proper motion space. The main-sequence contamination in both subsets of the catalogue was reduced, at the expense of completeness, by concentrating on the Milky Way's Galactic halo, where relatively young main-sequence stars were not expected. Results: We present a catalogue, based on Gaia EDR3, of 57,377 blue horizontal-branch stars. The Gaia EDR3 parallax was used in selecting 16,794 candidates and the proper motions were used to identify a further 40,583 candidates.

We numerically investigate the motion of stars on the meridional plane of an axially symmetric disk galaxy model, containing a central supermassive black hole, represented by the Paczynski-Wiita potential. By using this pseudo-Newtonian potential we can replicate important relativistic properties, such as the existence of the Schwarzschild radius. After classifying extensive samples of initial conditions of trajectories, we manage to distinguish between collisional, ordered, and chaotic motion. Besides, all starting conditions of regular orbits are further classified into families of regular orbits. Our results are presented through modern color-coded basin diagrams on several types of two-dimensional planes. Our analysis reveals that both the mass of the black hole (in direct relation with the Schwarzschild radius) as well as the angular momentum play an important role in the character of orbits of stars. More specifically, the trajectories of low angular momentum stars are highly affected by the mass of the black hole, while high angular momentum stars seem to be unaffected by the central black hole. Comparison with previous related outcomes, using Newtonian potentials for the central region of the galaxy, is also made.

We aim to use signatures of microlensing induced by stars in the foreground lens galaxy to infer the size of the accretion disk in the gravitationally lensed quasar Q 0957+561. The long-term photometric monitoring of this system (which so far has provided the longest available light curves of a gravitational lens system) permits us to evaluate the impact of uncertainties on our recently developed method (controlled by the distance between the modeled and the experimental magnitude difference histograms between two lensed images), and thus to test the robustness of microlensing-based disk-size estimates. We analyzed the well-sampled 21-year GLENDAMA optical light curves of the double-lensed quasar and studied the intrinsic and extrinsic continuum variations. Using accurate measurements for the time delay between the images A and B, we modeled and removed the intrinsic quasar variability, and from the statistics of microlensing magnifications we used a Bayesian method to derive the size of the region emitting the continuum at 2558 angstroms. Analyses of the Q 0957+561 R-band light curves show a slow but systematic increase in the brightness of the B relative to the A component during the past ten years. The relatively low strength of the magnitude differences between the images indicates that the quasar has an unusually big optical accretion disk of half-light radius $R_{1/2} = 17.6 \pm 6.1 \sqrt{M/0.3M_\odot}$ lt-days.

We investigate the possibility of detection of the VHE gamma-ray counterparts to the neutrino astrophysical sources within the Neutrino Target of Opportunity (NToO) program of CTA using the populations simulated by the FIRESONG software to resemble the diffuse astrophysical neutrino flux measured by IceCube. We derive the detection probability for different zenith angles and geomagnetic field configurations. The difference in detectability of sources between CTA-North and CTA-South for the average geomagnetic field is not substantial. We investigate the effect of a higher night-sky background and the preliminary CTA Alpha layout on the detection probability.

We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that in excess of a million SMBH shadows meeting these criteria are potentially accessible to interferometric observations with sufficient angular resolution and sensitivity. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Our model predicts that with modest improvements to sensitivity, as many as $\sim$5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experiments that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings.

Young massive stars influence their surroundings from local to galactic scales, but the observational challenges associated with their distance and embedded nature has, until the recent decade, made high-resolution studies of these objects difficult. In particular, comparative analyses of massive young stellar object (MYSO) discs are currently lacking and our understanding of their evolution is limited. Here, we combine the results of two studies with the aim to attribute geometrical features to an evolutionary sequence for a sample of seven MYSOs. The time evolution is based on a near-IR spectral features, while the geometry is determined from a multi size-scale study of MYSOs. We find that MYSO discs with determined geometrical substructure turn out to be also spectroscopically more evolved. This implies that disc evolution and dispersal are occurring within MYSOs, similar to low-mass YSO disc evolution, despite their faster formation timescales.

The high stellar density in the central parsecs around the Galactic Centre makes it a seemingly favourable environment for finding relativistic binary pulsars. These include pulsars orbiting other neutron stars, stellar-mass black holes or the central supermassive black hole, Sagittarius A*. Here we present multi-epoch pulsar searches of the Galactic Centre at four observing frequencies, (4.85, 8.35, 14.6 18.95) GHz, using the Effelsberg 100-m radio telescope. Observations were conducted one year prior to the discovery of, and during monitoring observations of, the Galactic Centre magnetar PSR J1745-2900. Our data analysis features acceleration searches on progressively shorter time series to maintain sensitivity to relativistic binary pulsars. The multi-epoch observations increase the likelihood of discovering transient or nulling pulsars, or ensure orbital phases are observed at which acceleration search methods work optimally. In ~147 h of separate observations, no previously undiscovered pulsars have been detected. Through calibration observations, we conclude this might be due to insufficient instantaneous sensitivity; caused by the intense continuum emission from the Galactic Centre, its large distance and, at higher frequencies, the aggregate effect of steep pulsar spectral indices and atmospheric contributions to the system temperature. Additionally we find that for millisecond pulsars in wide circular orbits ~<800 d around Sagittarius A*, linear acceleration effects cannot be corrected in deep observations (9 h) with existing software tools. Pulsar searches of the Galactic Centre with the next generation of radio telescopes - such as MeerKat, ngVLA and SKA1-mid - will have improved chances of uncovering this elusive population.

We present a multi-wavelength study of the gaseous medium surrounding the nearby active galactic nucleus (AGN) Fornax A. Using MeerKAT, ALMA and MUSE observations we reveal a complex distribution of the atomic (HI), molecular (CO), and ionised gas in its centre and along the radio jets. By studying the multi-scale kinematics of the multi-phase gas, we reveal the presence of concurrent AGN feeding and feedback phenomena. Several clouds and an extended 3 kpc filament -- perpendicular to the radio jets and the inner disk ($r\lesssim 4.5$ kpc) -- show highly-turbulent kinematics, which likely induces nonlinear condensation and subsequent Chaotic Cold Accretion (CCA) onto the AGN. In the wake of the radio jets and in an external ($r\gtrsim 4.5$ kpc) ring, we identify an entrained massive ($\sim$ $10^7$ M$_\odot$) multi-phase outflow ($v_{\rm OUT}\sim 2000$ km s$^{-1}$). The rapid flickering of the nuclear activity of Fornax A ($\sim$ 3 Myr) and the gas experiencing turbulent condensation raining onto the AGN provide quantitative evidence that a recurrent, tight feeding and feedback cycle may be self-regulating the activity of Fornax A, in agreement with CCA simulations. To date, this is one of the most in-depth probes of such a mechanism, paving the way to apply these precise diagnostics to a larger sample of nearby AGN hosts and their multi-phase ISM.

This is the index of all contributions of the H.E.S.S. Collaboration to the 37th International Cosmic-Ray Conference, held virtually, July 12 - 23, 2021.

Asteroid families with ages younger than $1$ Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes. We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on the characteristic dispersal velocity at subkilometer sizes corresponding to the smallest visible Hobson members. We identified the Hobson family using an up-to-date asteroid catalog. A significant increase in the number of its known members allowed us to study their size distribution and compare it with computer simulations of catastrophic disruptions. Backward orbital integrations of the heliocentric orbits allowed us to confirm the previously suggested age of Hobson and helped to estimate limits of the ejection speed. The Hobson family has an unusual size distribution: two nearly equal-size bodies, followed by a population of smaller asteroids, whose distribution takes a characteristic power law. There are two possibilities to explain these data. Either a canonical impact onto a single parent body, requiring fine-tuned impact conditions that have not been studied so far, or an unconventional model for the parent object of the Hobson family, namely a binary with $\simeq 7-9$ km primary and a $\simeq 2.5$ km secondary. In the latter case, the primary was disrupted, leaving behind the largest remnant (18777) Hobson and a suite of subkilometer asteroids. The second largest asteroid, (57738) 2001 UZ160, is the nearly intact satellite of the parent binary. The excellent convergence of nominal orbits of Hobson members sets an upper limit of $\simeq (10-20)$ m s$^{-1}$ for the initial dispersal velocity of the known members, which is consistent with both formation models.

We present for the first time the timing and spectral analyses for a narrow-line Seyfert 1 galaxy, SBS 1353+564, using \it{XMM-Newton} and \it{Swift} multi-band observations from 2007 to 2019. Our main results are as follows: 1) The temporal variability of SBS 1353+564 is random, while the hardness ratio is relatively constant over a time span of 13 years; 2) We find a prominent soft X-ray excess feature below 2 keV, which cannot be well described by a simple blackbody component; 3) After comparing the two most prevailing models for interpreting the origin of the soft X-ray excess, we find that the relativistically smeared reflection model is unable to fit the data above 5 keV well and the X-ray spectra do not show any reflection features, such as the Fe K\alpha emission line. However, the warm corona model can obtain a good fitting result. For the warm corona model, we try to use three different sets of spin values to fit the data and derive different best-fitting parameter sets; 4) We compare the UV/optical spectral data with the extrapolated values of the warm corona model to determine which spin value is more appropriate for this source, and we find that the warm corona model with non-spin can sufficiently account for the soft X-ray excess in SBS 1353+564.

The IceCube Neutrino Observatory opened the window on neutrino astronomy by discovering high-energy astrophysical neutrinos in 2013 and identifying the first compelling astrophysical neutrino source, the blazar TXS0506+056, in 2017. In this talk, we will discuss the science reach and ongoing development of the IceCube-Gen2 facility---a planned extension to IceCube. IceCube-Gen2 will increase the rate of observed cosmic neutrinos by an order of magnitude, be able to detect five-times fainter neutrino sources, and extend the measurement of astrophysical neutrinos several orders of magnitude higher in energy. We will discuss the envisioned design of the instrument, which will include an enlarged in-ice optical array, a surface array for the study of cosmic-rays, and a shallow radio array to detect ultra-high energy (>100 PeV) neutrinos. we will also highlight ongoing efforts to develop and test new instrumentation for IceCube-Gen2.

Galaxy cluster mergers are a powerful laboratory for testing cosmological and astrophysical models. However, interpreting individual merging clusters depends crucially on their merger configuration, defined by the masses, velocities, impact parameters, and orientation of the merger axis with respect to the plane of the sky. In this work, we investigate the impact of merger parameters on the X-ray emitting intracluster medium and gravitational lensing maps using a suite of idealised simulations of binary cluster mergers performed using the GAMER-2 code. As a test case, we focus on modeling the Bullet Cluster-like merging system Abell 2146, in which deep \textit{Chandra} X-ray and lensing observations revealed prominent merger shocks as well as the mass distribution and substructures associated with this merging cluster. We identify the most interesting parameter combinations, and evaluate the effects of various parameters on the properties of merger shocks observed by deep \textit{Chandra} and lensing observations. We show that due gravitational compression of the cluster halos during the merger, previous mass estimates from weak lensing are too high. The plane of the merger is tilted further from the plane of the sky than estimated previously, up to $30^\circ$ from the plane of the sky. We discuss the applicability of our results to multi-wavelength observations of merging galaxy clusters and their use as probes of cosmology and plasma physics.

The IceCube Neutrino Observatory at the geographic South Pole instruments a gigaton of glacial Antarctic ice with over 5000 photosensors. The detector, by now running for over a decade, will be upgraded with seven new densely instrumented strings. The project focuses on the improvement of low-energy and oscillation physics sensitivities as well as re-calibration of the existing detector. Over the last few years we developed a Precision Optical Calibration Module (POCAM) providing self-monitored, isotropic, nanosecond, light pulses for optical calibration of large-volume detectors. Over 20 next-generation POCAMs will be calibrated and deployed in the IceCube Upgrade in order to reduce existing detector systematics. We report a general overview of the POCAM instrument, its performance and calibration procedures.

X-ray SOI pixel sensors, "XRPIX", are being developed for the next-generation X-ray astronomical satellite, "FORCE". The XRPIX are fabricated with the SOI technology, which makes it possible to integrate a high-resistivity Si sensor and a low-resistivity Si CMOS circuit. The CMOS circuit in each pixel is equipped with a trigger function, allowing us to read out outputs only from the pixels with X-ray signals at the timing of X-ray detection. This function thus realizes high throughput and high time resolution, which enables to employ anti-coincidence technique for background rejection. A new series of XRPIX named XRPIX6E developed with a pinned depleted diode (PDD) structure improves spectral performance by suppressing the interference between the sensor and circuit layers. When semiconductor X-ray sensors are used in space, their spectral performance is generally degraded owing to the radiation damage caused by high-energy protons. Therefore, before using an XRPIX in space, it is necessary to evaluate the extent of degradation of its spectral performance by radiation damage. Thus, we performed a proton irradiation experiment for XRPIX6E for the first time at HIMAC in the NIRS. We irradiated XRPIX6E with high-energy protons with a total dose of up to 40 krad, equivalent to 400 years of irradiation in orbit. The 40-krad irradiation degraded the energy resolution of XRPIX6E by 25 $\pm$ 3%, yielding an energy resolution of 260.1 $\pm$ 5.6 eV at the full width half maximum for 5.9 keV X-rays. However, the value satisfies the requirement for FORCE, 300 eV at 6 keV, even after the irradiation. It was also found that the PDD XRPIX has enhanced radiation hardness compared to previous XRPIX devices. In addition, we investigated the degradation of the energy resolution; it was shown that the degradation would be due to increasing energy-independent components, e.g., readout noise.

Through the Backyard Worlds: Planet 9 citizen science project, we have identified a wide-separation ($\sim$10\arcmin, $\sim$9900 au projected) substellar companion to the nearby ($\sim$17.5 pc), mid-M dwarf Ross 19. We have developed a new formalism for determining chance alignment probabilities based on the BANYAN $\Sigma$ tool, and find a 100% probability that this is a physically associated pair. Through a detailed examination of Ross 19A, we find that the system is metal-poor ([Fe/H]=$-$0.40$\pm$0.12) with an age of 7.2$^{+3.8}_{-3.6}$ Gyr. Combining new and existing photometry and astrometry, we find that Ross 19B is one of the coldest known wide-separation companions, with a spectral type on the T/Y boundary, an effective temperature of 500$^{+115}_{-100}$ K, and a mass in the range 15--40 $M_{\rm Jup}$. This new, extremely cold benchmark companion is a compelling target for detailed characterization with future spectroscopic observations using facilities such as the Hubble Space Telescope or James Webb Space Telescope.

TOI-2202 b is a transiting warm Jovian-mass planet with an orbital period of P=11.91 days identified from the Full Frame Images data of five different sectors of the TESS mission. Ten TESS transits of TOI-2202 b combined with three follow-up light curves obtained with the CHAT robotic telescope show strong transit timing variations (TTVs) with an amplitude of about 1.2 hours. Radial velocity follow-up with FEROS, HARPS and PFS confirms the planetary nature of the transiting candidate (a$_{\rm b}$ = 0.096 $\pm$ 0.002 au, m$_{\rm b}$ = 0.98 $\pm$ 0.06 M$_{\rm Jup}$), and dynamical analysis of RVs, transit data, and TTVs points to an outer Saturn-mass companion (a$_{\rm c}$ = 0.155 $\pm$ 0.003 au, m$_{\rm c}$= $0.37 \pm 0.10$ M$_{\rm Jup}$) near the 2:1 mean motion resonance. Our stellar modeling indicates that TOI-2202 is an early K-type star with a mass of 0.82 M$_\odot$, a radius of 0.79 R$_\odot$, and solar-like metallicity. The TOI-2202 system is very interesting because of the two warm Jovian-mass planets near the 2:1 MMR, which is a rare configuration, and their formation and dynamical evolution are still not well understood.

Our recently reported lack of Intra-Night Optical Variability (INOV) among Broad-Absorption-Line (BAL) quasars exhibiting some blazar-like radio properties, either questions polar ejection of BAL clouds, and/or hints at a physical state of the relativistic jet modified due to interaction with the thermal BAL wind on the innermost sub-parsec scale. As a robust check on this unexpected finding for the BAL_blazar candidates, we report here the INOV study of a new and much more rigorously defined comparison sample consisting of 9 normal (non-BAL) blazars, matched in both magnitude and redshift to the aforementioned sample of BAL_blazar candidates. The present campaign spanning 27 sessions yields an INOV duty cycle of ~23% for the comparison sample of normal blazars, employing the enhanced F-test. However, even this more sensitive test does not detect INOV for the sample of BAL_blazar candidates. This stark INOV contrast found here between the BAL_blazar candidates and normal blazars can probably be traced to a physical interaction of the relativistic jet with the thermal wind, within sub-parsec range from the nucleus. The consequent enfeebling of the jet would additionally explain the striking deficiency among BAL quasars of powerful FR II radio lobes on the much larger scale of 10-100 kpc.

The highest-energy blazars exhibit non-thermal radiation extending beyond 1 TeV with high luminosities and strong variabilities, indicating extreme particle acceleration in their relativistic jets. The gamma-ray spectra of blazars contain information about the distribution and cooling processes of high-energy particles in jets, the extragalactic background light between the source and the observer, and potentially, the environment of the gamma-ray emitting region and exotic physics that may modify the opacity of the universe to gamma rays. We use data from Fermi-LAT and VERITAS to study the variability and spectra of a sample of TeV blazars across a wide range of gamma-ray energies, taking advantage of more than ten years of data from both instruments. The variability in both the GeV and TeV gamma-ray bands is investigated using a Bayesian blocks method to identify periods with a steady flux, during which the average gamma-ray spectra, after correcting for the pair absorption effect from propagation, can be parameterized without the risk of mixing different flux states. We report on the search for intrinsic spectral curvature and spectral variability in these blazars, in an effort to understand the physical mechanisms behind the high-energy gamma-ray spectra of TeV blazars.

We explore the potential of Pulsar Timing Arrays (PTAs) such as NANOGrav, EPTA, and PPTA to detect the Stochastic Gravitational Wave Background (SGWB) in theories of massive gravity. In General Relativity, the function describing the dependence of the correlation between the arrival times of signals from two pulsars on the angle between them is known as the Hellings-Downs curve. We compute the analogous overlap reduction function for massive gravity, including the additional polarization states and the correction due to the mass of the graviton, and compare the result with the Hellings-Downs curve. The primary result is a complete analytical form for the analog Hellings-Downs curve, providing a starting point for future numerical studies aimed at a detailed comparison between PTA data and the predictions of massive gravity. We study both the massless limit and the stationary limit as checks on our calculation, and discuss how our formalism also allows us to study the impact of massive spin-2 dark matter candidates on data from PTAs.

The recently introduced swampland criterion for de Sitter (arXiv:1806.08362) can be viewed as a (hierarchically large) bound on the smallness of the slow roll parameter $\epsilon_V$. This leads us to consider the other slow roll parameter $\eta_V$ more closely, and we are lead to conjecture that the bound is not necessarily on $\epsilon_V$, but on slow roll itself. A natural refinement of the de Sitter swampland conjecture is therefore that slow roll is violated at ${\cal O}(1)$ in Planck units in any UV complete theory. A corollary is that $\epsilon_V$ need not necessarily be ${\cal O}(1)$, if $\eta_V \lesssim -{\cal O}(1)$ holds. We consider various tachyonic tree level constructions of de Sitter in IIA/IIB string theory (as well as closely related models of inflation), which superficially violate arXiv:1806.08362, and show that they are consistent with this refined version of the bound. The phrasing in terms of slow roll makes it plausible why both versions of the conjecture run into trouble when the number of e-folds during inflation is high. We speculate that one way to evade the bound could be to have a large number of fields, like in $N$-flation.

Gravitational wave echoes may provide a smoking gun signal for new physics in the immediate vicinity of black holes. As a quasi-periodic signal in time, echoes are characterized by the nearly constant time delay, and its precise measurement can help reveal a Planck scale deviation outside of the would-be horizon. Different search methods have been developed for this quasi-periodic signal, while the searches suffer from large theoretical uncertainties of the echo waveform associated with the near-horizon physics. On the other hand, a coherent combine of a large number of pulses gives rise to a generic narrow resonance structure for the echo amplitude in frequency. The quasi-periodic resonance structure sets a complementary search target for echoes, and the time delay is inversely related to the average resonance spacing. A uniform comb has been proposed to look for the resonance structure in a rather model independent way. In this paper, we develop a Bayesian algorithm to search for the resonance structure based on combs, where a phase-marginalized likelihood plays an essential role. The algorithm is validated with signal injections in detector noise from Advanced LIGO. With special treatments of the non-Gaussian artifacts, the noise outliers of the log Bayes factor distribution are properly removed. An echo signal not significantly below noise is detectable, and the time delay can be determined to very high precision. We perform the proposed search on real gravitational wave strain data of the first observing run of Advanced LIGO. We find no clear evidence of a comb-like structure for GW150914 and GW151012.

One contribution to any dark sector's abundance comes from its gravitational production during inflation. If the dark sector is weakly coupled to the inflaton and the Standard Model, this can be its only production mechanism. For non-interacting dark sectors, such as a free massive fermion or a free massive vector field, this mechanism has been studied extensively. In this paper we show, via the example of dark massive QED, that the presence of interactions can result in a vastly different mass for the dark matter (DM) particle, which may well coincide with the range probed by upcoming experiments. In the context of dark QED we study the evolution of the energy density in the dark sector after inflation. Inflation produces a cold vector condensate consisting of an enormous number of bosons, which via interesting processes - Schwinger pair production, strong field electromagnetic cascades, and plasma dynamics - transfers its energy to a small number of "dark electrons" and triggers thermalization of the dark sector. The resulting dark electron DM mass range is from 50 MeV to 30 TeV, far different from both the $10^{-5}$ eV mass of the massive photon dark matter in the absence of dark electrons, and from the $10^9$ GeV dark electron mass in the absence of dark photons. This can significantly impact the search strategies for dark QED and, more generally, theories with a self-interacting DM sector. In the presence of kinetic mixing, a dark electron in this mass range can be searched for with upcoming direct detection experiments, such as SENSEI-100g and OSCURA.

The Horizon-10T collaboration \cite{Beisembaev:2016cyg,2017EPJWC.14514001B,Beznosko:2019cI,2019EPJWC.20806002B,Beisembaev:2019nzd} have reported observation of Multi-Modal Events (MME) containing multiple peaks suggesting their clustering origin. These events are proven to be hard to explain in terms of conventional cosmic rays (CR). We propose that these MMEs might be result of the dark matter annihilation events within the so-called axion quark nugget (AQN) dark matter model, which was originally invented for completely different purpose to explain the observed similarity between the dark and the visible components in the Universe, i.e. $\Omega_{\rm DM}\sim \Omega_{\rm visible}$ without any fitting parameters. We support this proposal by demonstrating that the observations \cite{Beisembaev:2016cyg,2017EPJWC.14514001B,Beznosko:2019cI,2019EPJWC.20806002B,Beisembaev:2019nzd}, including the frequency of appearance, intensity, the spatial distribution, the time duration, the clustering features, and many other properties nicely match the emission characteristics of the AQN annihilation events in atmosphere. We list a number of features of the AQN events which are very distinct from conventional CR air showers. The observation (non-observation) of these features may substantiate (refute) our proposal.

Secondary school teachers often lack the necessary content background in astronomy to teach such a course confidently. Our theory of change postits that an increased confidence level will increase student retention in astronomy and related STEM fields. Beyond the science content knowledge though, teachers need opportunities to embed the content in pedagogically sound practices, and with appropriate technology tools. We report on our interdisciplinary approach to designing, developing, fielding, and iteratively improving the San Antonio Teacher Training Astronomy Academy (SATTAA), an annually offered Teacher Professional Development program. In particular, we present how our separate areas of expertise, in content and in STEM pedagogy, led to a synergistic process of teacher professional development that has now resulted in three cohorts of alumni. In this paper, we share our interdisciplinary processes and lessons learned; program metrics are described elsewhere in detail.

We review the ab initio symmetry-adapted (SA) framework for determining the structure of stable and unstable nuclei, along with related electroweak, decay and reaction processes. This framework utilizes the dominant symmetry of nuclear dynamics, the shape-related symplectic Sp(3,R) symmetry, which has been shown to emerge from first principles and to expose dominant degrees of freedom that are collective in nature, even in the lightest species or seemingly spherical states. This feature is illustrated for a broad scope of nuclei ranging from helium to titanium isotopes, enabled by recent developments of the ab initio symmetry-adapted no-core shell model expanded to the continuum through the use of the SA basis and that of the resonating group method. The review focuses on energies, electromagnetic transitions, quadrupole and magnetic moments, radii, form factors, and response function moments, for ground-state rotational bands and giant resonances. The method also determines the structure of reaction fragments that is used to calculate decay widths and alpha-capture reactions for simulated x-ray burst abundance patterns, as well as nucleon-nucleus interactions for cross sections and other reaction observables.

We propose a new Harten-Lax-van Leer discontinuities (HLLD) approximate Riemann solver to improve the stability of shocks and the accuracy of low-speed flows in multidimensional magnetohydrodynamic (MHD) simulations. Stringent benchmark tests verify that the new solver is more robust against a numerical shock instability and is more accurate for low-speed, nearly incompressible flows than the original solver, whereas additional computational costs are quite low. The novel ability of the new solver enables us to tackle MHD systems, including both high and low Mach number flows.

We analyze the IGETS based gravitational acceleration measurements to search for a variability in Newton's constant $G$, as a complement to the analysis carried out in arXiv:2103.11157. To understand the variation in gravitational acceleration, we fit the aforementioned data to a constant model, two sinusoidal models, and a linear model. We find that none of the four models provide a good fit to the data, showing that there is no evidence for periodicity or linear temporal variation in the acceleration measurements. Based on the sinusoidal models, we infer that the relative variance of $G$ is $\leqslant 2 \times 10^{-9}$, which is about four orders of magnitude more stringent than the amplitude of periodic variations inferred from previous $G$ measurements. From the linear model, we obtain $\frac{\dot{G}}{G} < 2.14 \times 10^{-10} \rm{yr^{-1}}$.

After scattering off nuclei in the Sun, dark matter particles can be gravitationally captured by the Sun, accumulate in the Sun's core and annihilate into Standard Model particles. Neutrinos originating from these annihilations can be detected by the IceCube Neutrino Observatory, located at the South Pole. Due to the non-observation of these neutrinos, constraints on the standard spin-dependent and spin-independent dark matter-nucleon scattering cross sections have been placed. Based on these constraints, we present upper limits on the coupling constants of the non-relativistic effective theory of dark matter-nucleon interactions, including velocity and momentum dependent interactions.

We propose the use of trapped ions for detection of millicharged dark matter. Millicharged particles will scatter off the ions, giving a signal either in individual events or in the overall heating rate of the ions. Ion traps have several properties which make them ideal detectors for such a signal. First, ion traps have demonstrated significant isolation of the ions from the environment, greatly reducing the background heating and event rates. Second, ion traps can have low thresholds for detection of energy deposition, down to $\sim \text{neV}$. Third, since the ions are charged, they naturally have large cross sections for scattering with the millicharged particles, further enhanced by the low velocities of the thermalized millicharges. Despite ion-trap setups being optimized for other goals, we find that existing measurements put new constraints on millicharged dark matter which are many orders of magnitude beyond previous bounds. For example, for a millicharge dark matter mass $m_Q=10~\textrm{GeV}$ and charge $10^{-3}$ of the electron charge, ion traps limit the local density to be $n_Q \lesssim 1 \, \textrm{cm}^{-3}$, a factor $\sim 10^8$ better than current constraints. Future dedicated ion trap experiments could reach even further into unexplored parameter space.

High energy neutrinos from SN1987A were searched for using upward-going muons recorded by the Kamiokande-II experiment and the IMB experiment. Between August 11 and October 20, 1987, and from an angular window of 10 degree radius, 2 upward-going muon events were recorded by Kamiokande-II, and also 2 events were recorded by IMB. The probability that these upward-going muons were explained by a chance coincidence of atmospheric neutrinos was calculated to be 0.27%. This shows a possible evidence of high energy neutrinos from SN1987A.