Papers for Thursday, Jan 28 2021

We report the detection of 23 OH+(1-0) absorption, emission, or P-Cygni-shaped lines and CO(9-8) emission lines in 18 Herschel-selected z=2-6 starburst galaxies with ALMA and NOEMA, taken as part of the Gas And Dust Over cosmic Time (GADOT) Galaxy Survey. We find that the CO(9-8) luminosity is higher than expected based on the far-infrared luminosity when compared to nearby star-forming galaxies. Together with the strength of the OH+ emission components, this may suggest that shock excitation of warm, dense molecular gas is more prevalent in distant massive dusty starbursts than in nearby star-forming galaxies on average, perhaps due to an impact of galactic winds on the gas. OH+ absorption is found to be ubiquitous in massive high-redshift starbursts, and is detected toward 89% of the sample. The majority of the sample shows evidence for outflows or inflows based on the velocity shifts of the OH+ absorption/emission, with a comparable occurrence rate of both at the resolution of our observations. A small subsample appears to show outflow velocities in excess of their escape velocities. Thus, starburst-driven feedback appears to be important in the evolution of massive galaxies in their most active phases. We find a correlation between the OH+ absorption optical depth and the dust temperature, which may suggest that warmer starbursts are more compact and have higher cosmic ray energy densities, leading to more efficient OH+ ion production. This is in agreement with a picture in which these high-redshift galaxies are "scaled-up" versions of the most intense nearby starbursts.

To date, about two dozen low-mass embedded protostars exhibit rich spectra with complex organic molecule (COM) lines. These protostars seem to possess different enrichment in COMs. However, the statistics of COM abundance in low-mass protostars are limited by the scarcity of observations. This study introduces the Perseus ALMA Chemistry Survey (PEACHES), which aims at unbiasedly characterizing the chemistry of COMs toward the embedded (Class 0/I) protostars in the Perseus molecular cloud. Among the 50 embedded protostars surveyed, 58% of them have emission from COMs. A 56%, 32%, and 40% of protostars have CH$_3$OH, CH$_3$OCHO, and N-bearing COMs, respectively. The detectability of COMs depends on neither the averaged continuum brightness temperature, a proxy of the H$_2$ column density, nor the bolometric luminosity and the bolometric temperature. For the protostars with detected COMs, CH$_3$OH has a tight correlation with CH$_3$CN, spanning more than two orders of magnitude in column densities normalized by the continuum brightness temperature, suggesting a chemical relation between CH$_3$OH and CH$_3$CN and large chemical diversity among the PEACHES samples at the same time. A similar trend with more scatter is also found between all identified COMs, hinting at common chemistry for the sources with COMs. The correlation between COMs is insensitive to the protostellar properties, such as the bolometric luminosity and the bolometric temperature. The abundance of larger COMs (CH$_3$OCHO and CH$_3$OCH$_3$) relative to that of smaller COMs (CH$_3$OH and CH$_3$CN) increases with the inferred gas column density, hinting at an efficient production of complex species in denser envelopes.

This document describes the general process of setting up, running, and analysing disc galaxy simulations using the freely available program Phantom of RAMSES (PoR). This implements Milgromian Dynamics (MOND) with a patch to the RAMSES grid-based $N$-body and hydrodynamical code that uses adaptive mesh refinement. We discuss the procedure of setting up isolated and interacting disc galaxy initial conditions for PoR, running the simulations, and analysing the results. This manual also concisely documents all previously developed MOND simulation codes and the results obtained with them.

We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model from D\'iaz-Rodr\'guez et al. (2018). For this inference, we use the local star formation histories to estimate the age for each supernova (SN) as inferred in Williams et al. (2018). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are random bursts of star formation. The primary inference is the progenitor age distribution. Due to the limited number of historic SNe and highly uncertain star formation at young ages, we restrict our inference to the slope of the age distribution and the maximum age for CCSNe. Using single-star evolutionary models, we transform the progenitor age distribution into a progenitor mass distribution. Under these assumptions, the minimum mass for CCSNe is ${M_\textrm{min}}~=~8.60^{+0.37}_{-0.41} \ M_\odot$ and the slope of the progenitor mass distribution is $\alpha = -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($\alpha = -2.35$). These values are consistent with previous estimates using precursor imaging and the age-dating technique, further confirming that using stellar populations around SN and supernova remnants is a reliable way to infer the progenitor masses.

High-redshift star-forming galaxies have very different morphologies compared to nearby ones. Indeed, they are often dominated by bright star-forming structures of masses up to $10^{8-9}$ $\mathrm{M}_\odot$ dubbed {\guillemotleft}giant clumps{\guillemotright}. However, recent observations questioned this result by showing only low-mass structures or no structure at all. We use Adaptative Mesh Refinement hydrodynamical simulations of galaxies with parsec-scale resolution to study the formation of structures inside clumpy high-redshift galaxies. We show that in very gas-rich galaxies star formation occurs in small gas clusters with masses below $10^{7-8}$ $\mathrm{M}_\odot$ that are themselves located inside giant complexes with masses up to $10^8$ and sometimes $10^9$ $\mathrm{M}_\odot$. Those massive structures are similar in mass and size to the giant clumps observed in imaging surveys, in particular with the Hubble Space Telescope. Using mock observations of simulated galaxies, we show that at very high resolution with instruments like the Atacama Large Millimeter Array or through gravitational lensing, only low-mass structures are likely to be detected, and their gathering into giant complexes might be missed. This leads to the non-detection of the giant clumps and therefore introduces a bias in the detection of these structures. We show that the simulated giant clumps can be gravitationally bound even when undetected in mocks representative for ALMA observations and HST observations of lensed galaxies. We then compare the top-down fragmentation of an initially warm disc and the bottom-up fragmentation of an initially cold disc to show that the process of formation of the clumps does not impact their physical properties.

We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from $N$-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of $N$-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these nonlinear spectra from redshifts $z = 0$ to $z=2$. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h {\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift $z\simeq 0.4$ to scales of $k_{\rm max} \approx 0.6\, h\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy--galaxy lensing analyses with data drawn from an independent $N$-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming datasets from galaxy surveys.

Many observations in recent times have shown evidence against the standard assumption of isotropy in the Big Bang model. Introducing a superhorizon scalar metric perturbation has been able to explain some of these anomalies. In this work, we probe the net velocity arising due to the perturbation, which does not cancel out for large scale structure, unlike in the case of CMB. Thus, within this model's framework, our velocity with respect to the CMB is different from the velocity with respect to the large scale structure. Taking this extra velocity component into account, we study the superhorizon mode's implications for the excess dipole observed in the NRAO VLA Sky Survey (NVSS). We find that the mode can consistently explain both the CMB and NVSS observations. We also find that the model is consistent with the observed Hubble constant dipole and the Hubble bulk flow velocity. The model leads to several predictions which can be tested in future surveys. In particular, it implies that the observed dipole in large scale structure should be redshift dependent and should show an increase in amplitude with redshift. We also find that the Hubble parameter should show a dipole anisotropy whose amplitude must increase with redshift in the CMB frame. Similar anisotropic behaviour is expected for the observed redshift as a function of the luminosity distance.

(Abridged) Several phenomena in astrophysics generate light curves with time delays. Among these are reverberation mapping, and lensed quasars. In some systems, the measurement of the time-delay is complicated by the fact that the delayed components are unresolved and that the light curves are generated from a red-noise process. We derive the likelihood function of the observations given a model of either a combination of time-delayed light curves or a single light curve. This likelihood function is different from the auto-correlation function. We demonstrate that given a single-band light curve that is a combination of two (or more) time-shifted copies of an original light curve, generated from a red-noise probability distribution, we can test if the total-flux light curve is a composition of time-delayed copies or, alternatively, is consistent with being the original light curve. Furthermorew, in some realistic cases, it is possible to measure the time delays and flux ratios between these unresolved components even when the flux ratio is about 1/10. This method is useful for identifying lensed quasars and simultaneously measuring their time delays, and for estimating the reverberation time scales of active galactic nuclei. In a companion paper, we derive a method that uses the center-of-light position (e.g., of a lensed quasar) along with the combined flux. This allow us to identify lensed quasars and supernovae and measure their time delays, with higher fidelity compared to the flux-only method. The astrometry + flux method, however, is not suitable for quasar reverberation mapping. We also comment on the commonly used method of fitting a power-law model to a power spectrum, and present the proper likelihood function for such a fit. We test the new method on simulations and provide Python and MATLAB implementations.

The physical origin of the scatter in the relation between galaxy stellar mass and the metallicity of the interstellar medium, i.e. the Mass-Metallicity Relation (MZR), reflects the relative importance of key processes in galaxy evolution. The \eagle cosmological hydrodynamical simulation is used to investigate the correlations between the residuals of the MZR and the residuals of the relations between stellar mass and, respectively, specific inflow, outflow and star formation rate as well as the gas fraction for central galaxies. At low redshift, all these residuals are found to be anti-correlated with the residuals of the MZR for $M_\star/\mathrm{M}_\odot \lesssim 10^{10}$. The correlations between the residuals of the MZR and the residuals of the other relations with mass are interrelated, but we find that gas fraction, specific inflow rate and specific outflow rate all have at least some independent influence on the scatter of the MZR. We find that, while for $M_\star/\mathrm{M}_\odot > 10^{10.4}$ the specific mass of the nuclear black hole is most important, for $M_\star/\mathrm{M}_\odot \lesssim 10^{10.3}$ gas fraction and specific inflow rate are the variables that correlate most strongly with the MZR scatter. The timescales involved in the residual correlations and the time that galaxies stay above the MZR are revealed to be a few Gyr. However, most galaxies that are below the MZR at $z=0$ have been below the MZR throughout their lifetimes.

The levels of heavy elements in stars are the product of enhancement by previous stellar generations, and the distribution of this metallicity among the population contains clues to the process by which a galaxy formed. Most famously, the "G-dwarf problem" highlighted the small number of low-metallicity G-dwarf stars in the Milky Way, which is inconsistent with the simplest picture of a galaxy formed from a "closed box" of gas. It can be resolved by treating the Galaxy as an open system that accretes gas throughout its life. This observation has classically only been made in the Milky Way, but the availability of high-quality spectral data from SDSS-IV MaNGA and the development of new analysis techniques mean that we can now make equivalent measurements for a large sample of spiral galaxies. Our analysis shows that high-mass spirals generically show a similar deficit of low-metallicity stars, implying that the Milky Way's history of gas accretion is common. By contrast, low-mass spirals show little sign of a G-dwarf problem, presenting the metallicity distribution that would be expected if such systems evolved as pretty much closed boxes. This distinction can be understood from the differing timescales for star formation in galaxies of differing masses.

Lensed quasars and supernovae can be used to study galaxies' gravitational potential and measure cosmological parameters. The typical image separation of objects lensed by galaxies is of the order of 0.5". Therefore, finding the ones with small separations, and measuring their time-delays using ground-based observations is challenging. We suggest a new method to identify lensed quasars and simultaneously measure their time-delays, using seeing-limited synoptic observations in which the lensed quasar images and the lensing galaxy are unresolved. We show that using the light curve of the combined flux, and the astrometric measurements of the center-of-light position of the lensed images, the lensed nature of a quasar can be identified, and its time-delay can be measured. We provide the analytic formalism to do so, taking into account the measurement errors and the fact that the power spectra of quasar light curves is red (i.e., the light curve is highly correlated). We demonstrate our method on simulated data, while its implementation to real data will be presented in future papers. Our simulations suggest that, under reasonable assumptions, the new method can detect unresolved lensed quasars and measure their time delays, even when the image separation is below 0.1", or the flux ratio between the faintest and brightest images is as low as 0.03. Python and MATLAB implementations are provided. In a companion paper, we present a method for measuring the time delay using the combined flux observations. Although the flux-only method is less powerful, it may be useful in cases in which the astrometric information is not relevant (e.g., reverberation mapping).

The size of convective cores remains uncertain, despite its substantial influence on stellar evolution, and thus on stellar ages. The seismic modeling of young subgiants can be used to obtain indirect constraints on the core structure during main sequence, thanks to the high probing potential of mixed modes. We selected the young subgiant KIC10273246, observed by Kepler, based on its mixed-mode properties. We thoroughly modeled this star, with the aim of placing constraints on the size of its main sequence convective core. We first extracted the parameters of the oscillation modes of the star using the full Kepler data set. To overcome the challenges posed by the seismic modeling of subgiants, we proposed a method which is specifically tailored for subgiants with mixed modes and consists in a nested optimization. We then applied this method to perform a detailed seismic modeling of KIC10273246. We obtained models that show good statistical agreements with the observations, both seismic and non-seismic. We showed that including core overshooting in the models significantly improves the quality of the seismic fit, optimal models being found for $\alpha_{\mathrm{ov}} = 0.15$. Higher amounts of core overshooting strongly worsen the agreement with the observations and are thus firmly ruled out. We also found that having access to two g-dominated mixed modes in young subgiants allows us to place stronger constraints on the gradient of molecular weight in the core and on the central density. This study confirms the high potential of young subgiants with mixed modes to investigate the size of main-sequence convective cores. It paves the way for a more general study including the subgiants observed with Kepler, TESS, and eventually PLATO.

The HAWC Collaboration has observed gamma rays at energies above 56 TeV from a collection of nine sources. It has been suggested that this emission could be hadronic in nature, requiring that these systems accelerate cosmic-ray protons or nuclei up to PeV-scale energies. In this paper, we instead show that the spectra of these objects favors a leptonic (inverse Compton) origin for their emission. More specifically, the gamma-ray emission from these objects can be straightforwardly accommodated within a model in which $\sim$\,10-20\% of the host pulsar's spindown power is transferred into the acceleration of electrons and positrons with a power-law spectrum that extends to several hundred TeV or higher. The spectral break that is observed among these sources is naturally explained within the context of this simple model, and occurs at the energy where the timescale for energy losses matches the age of the pulsar. In contrast, this spectral feature cannot be straightforwardly accommodated in hadronic scenarios. Furthermore, hadronic models predict that these sources should produce more emission at GeV-scale energies than is observed. In light of these considerations, we conclude that HAWC's highest energy sources should be interpreted as TeV halos or pulsar wind nebulae, which produce their emission through inverse Compton scattering, and are powered by the rotational kinetic energy of their host pulsar.

An excess of $\gamma$ rays in the data measured by the Fermi Large Area Telescope in the direction of the Galactic center has been reported in several publications. This excess, labeled as the Galactic center excess (GCE), is detected analyzing the data with different interstellar emission models, point source catalogs and analysis techniques. The characteristics of the GCE, recently measured with unprecedented precision, are all compatible with dark matter particles (DM) annihilating in the main halo of our Galaxy, even if other interpretations are still not excluded. We investigate the DM candidates that fit the observed GCE spectrum and spatial morphology. We assume a simple scenario with DM annihilating into a single channel but we inspect also more complicated models with two and three channels. We perform a search for a $\gamma$-ray flux from a list of 48 Milky Way dwarf spheroidal galaxies (dSphs) using state-of-the-art estimation of the DM density in these objects. Since we do not find any significant signal from the dSphs, we put upper limits on the annihilation cross section that result to be compatible with the DM candidate that fits the GCE. However, we find that the GCE DM signal is excluded by the AMS-02 $\bar{p}$ flux data for all hadronic and semi-hadronic annihilation channels unless the vertical size of the diffusion halo is smaller than 2 kpc -- which is in tension with radioactive cosmic ray fluxes and radio data. Furthermore, AMS-02 $e^+$ data rule out pure or mixed channels with a component of $e^+ e^-$. The only DM candidate that fits the GCE spectrum and is compatible with constraints obtained with the combined dSphs analysis and the AMS-02 $\bar{p}$ and $e^+$ data annihilates purely into $\mu^+\mu^-$, has a mass of 60 GeV and roughly a thermal cross section.

It is widely believed that super-Eddington accretion flow can produce powerful outflow, but where it originates from and how much mass and energy are carried away to which directions? To answer to these questions, we newly perform a large-box, two-dimensional radiation hydrodynamic simulation, paying special attention lest the results should depend on adopted initial and boundary conditions. We could achieve a quasi-steady state in an unprecedentedly large range, $r=2~r_{\rm S}$-$600~r_{\rm S}$ (with $r_{\rm S}$ being the Schwarzschild radius) from the black hole. The accretion rate onto the central $10 ~M_{\odot}$ black hole is $\dot{M}_{\rm BH} \sim 180 ~L_{\rm Edd}/c^{2}$, whereas the mass outflow rate is ${\dot M}_{\rm outflow} \sim 24 ~L_{\rm Edd}/c^2$ (where $L_{\rm Edd}$ and $c$ are the Eddington luminosity and the speed of light, respectively). The ratio (${\dot M}_{\rm outflow}/{\dot M}_{\rm BH} \sim 0.14$) is much less than those reported previously. By careful inspection we find that most of outflowing gas which reach the outer boundary originates from the region at $R\lesssim140~r_{\rm S}$, while gas at $140~r_{\rm S}$-$230 ~r_{\rm S}$ forms failed outflow. Therefore, significant outflow occurs inside the trapping radius $\sim 450 ~r_{\rm S}$. The mechanical energy flux (or mass flux) reaches its maximum in the direction of $\sim 15^\circ$ ($\sim 80^\circ$) from the rotation axis. The total mechanical luminosity is $L_{\rm mec}\sim 0.16~L_{\rm Edd}$, while the isotropic X-ray luminosity varies from $L_{\rm X}^{\rm ISO}\sim 2.9~L_{\rm Edd}$, (for a face-on observer) to $\sim 2.1~L_{\rm Edd}$ (for a nearly edge-on observer). The power ratio is $L_{\rm mec}/L_{\rm X}^{\rm ISO}\sim 0.05$-$0.08$, in good agreement with the observations of Ultra-Luminous X-ray sources surrounded by optical nebulae.

We present an analysis of the $R\lesssim 1.5$ kpc core regions of seven simulated Milky Way mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period ($\Delta t = 2.2$ Myr) of 22 Myr at $z \approx 0$, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (CMZ). Despite not being tuned to reproduce the detailed structure of the CMZ, we find that four of these galaxies are consistent with CMZ observations at some point during this 22 Myr period. The galaxies presented here are not homogeneous in their central structures, roughly dividing into two morphological classes; (a) several of the galaxies have very asymmetric gas and SFR distributions, with intense (compact) starbursts occurring over a period of roughly 10 Myr, and structures on highly eccentric orbits through the CMZ, whereas (b) others have smoother gas and SFR distributions, with only slowly varying SFRs over the period analyzed. In class (a) centers, the orbital motion of gas and star-forming complexes across small apertures ($R \lesssim 150$pc, analogously $|l|<1^\circ$ in the CMZ observations) contributes as much to tracers of star formation/dense gas appearing in those apertures, as the internal evolution of those structures does. These asymmetric/bursty galactic centers can simultaneously match CMZ gas and SFR observations, demonstrating that time-varying star formation can explain the CMZ's low star formation efficiency.

We use photometry and proper motions from Gaia DR2 to determine the blue straggler star (BSS) populations of 16 old (1-10 Gyr), nearby ($d< 3500$ pc) open clusters. We find that the fractional number of BSS compared to RGB stars increases with age, starting near zero at 1 Gyr and flattening to $\sim 0.35$ by 4 Gyr. Fitting stellar evolutionary tracks to these BSS, we find that their mass distribution peaks at a few tenths of a solar mass above the main-sequence turnoff. BSS more than 0.5 $M_\odot$ above the turnoff make up only $\sim25$\% of the sample, and BSS more than 1.0 $M_\odot$ above the turnoff are rare. We compare this to population synthesis models of BSS formed via mass transfer using the Compact Object Synthesis and Monte Carlo Investigation Code (COSMIC). We find that standard population synthesis assumptions dramatically under-produce the number of BSS in old open clusters. We also find that these models overproduce high mass BSS relative to lower mass BSS. Expected numbers of BSS formed through dynamics do not fully account for this discrepancy. We conclude that in order to explain the observed BSS populations from Roche lobe overflow, mass-transfer from giant donors must be more stable than assumed in canonical mass-transfer prescriptions, and including non-conservative mass transfer is important in producing realistic BSS masses. Even with these modifications, it is difficult to achieve the large number of BSS observed in the oldest open clusters. We discuss some additional physics that may explain the large number of observed blue stragglers among old stellar populations.

We measured $^{35}$Cl abundances in 52 M giants with metallicities between -0.5 $<$ [Fe/H] $<$ 0.12. Abundances and atmospheric parameters were derived using infrared spectra from CSHELL on the IRTF and from optical echelle spectra. We measured Cl abundances by fitting a H$^{35}$Cl molecular feature at 3.6985 $\mu$m with synthetic spectra. We also measured the abundances of O, Ca, Ti, and Fe using atomic absorption lines. We find that the [Cl/Fe] ratio for our stars agrees with chemical evolution models of Cl and the [Cl/Ca] ratio is broadly consistent with the solar ratio over our metallicity range. Both indicate that Cl is primarily made in core-collapse supernovae with some contributions from Type Ia SN. We suggest other potential nucleosynthesis processes, such as the $\nu$-process, are not significant producers of Cl. Finally, we also find our Cl abundances are consistent with H II and planetary nebular abundances at a given oxygen abundance, although there is scatter in the data.

We model the Neupert effect that relates flare heating energy with the observed SXR emission. The traditional form of the Neupert effect refers to the correlation between the time-integrated HXR or microwave light curve and the SXR light curve. In this paper, instead, we use as the proxy for heating energy the ultraviolet (UV) emission at the foot-points of flare loops, and modify the model of the Neupert effect by taking into account the discrete nature of flare heating as well as cooling. In the modified empirical model, spatially resolved UV lightcurves from the transition region or upper chromosphere are each convolved with a kernel function characterizing the decay of the flare loop emission. Contributions by all loops are summed to compare with the observed total SXR emission. The model has successfully reproduced the observed SXR emission from its rise to decay. To estimate heating energies in flare loops, we also employ the UV Foot-point Calorimeter (UFC) method that infers heating rates in flare loops from these UV light curves and models evolution of flare loops with a zero-dimensional hydrodynamic code. The experiments show that a multitude of impulsive heating events do not well reproduce the observed flare SXR light curve, but a two-phase heating model leads to better agreement with observations. Comparison of the two models of the Neupert effect further allows us to calibrate the UFC method, and improve the estimate of heating rates in flare loops continuously formed by magnetic reconnection throughout the flare evolution.

We present a Markov Chain Monte Carlo (MCMC)-based parameter estimation package, CosmoReionMC, to jointly constrain cosmological parameters of the $\Lambda$CDM model and the astrophysical parameters related to hydrogen reionization. The package is based on a previously developed physically motivated semi-analytical model for reionization, a similar semi-analytical model for computing the global 21~cm signal during the cosmic dawn and using an appropriately modified version of the publicly available CAMB for computing the CMB anisotropies. These calculations are then coupled to an MCMC ensemble sampler \texttt{emcee} to compute the posterior distributions of the model parameter. The model has twelve free parameters in total: five cosmological and seven related to the stellar populations. We constrain the parameters by matching the theoretical predictions with CMB data from Planck, observations related to the quasar absorption spectra and, for the first time, the global 21~cm signal from EDGES. We find that incorporating the quasar spectra data in the analysis tightens the bounds on the electron scattering optical depth $\tau$ and consequently the normalization $A_s$ of the primordial matter power spectrum (or equivalently $\sigma_8$). Furthermore, when we include the EDGES data in the analysis, we find that an early population of metal-free stars with efficient radio emission is necessary to match the absorption amplitude. The CosmoReionMC package should have interesting future applications, e.g., probing non-standard extensions to the $\Lambda$CDM model.

The LIGO-Virgo Collaboration has so far detected around 90 black holes, some of which have masses larger than what were expected from the collapse of stars. The mass distribution of LIGO-Virgo black holes appears to have a peak at $\sim30M_{\odot}$ and two tails on the ends. By assuming that they all have a primordial origin, we analyze the GWTC-1 (O1\&O2) and GWTC-2 (O3a) datasets by performing maximum likelihood estimation on a broken power law mass function $f(m)$, with the result $f\propto m^{1.2}$ for $m<35M_{\odot}$ and $f\propto m^{-4}$ for $m>35M_{\odot}$. This appears to behave better than the popular log-normal mass function. Surprisingly, such a simple and unique distribution can be realized in our previously proposed mechanism of PBH formation, where the black holes are formed by vacuum bubbles that nucleate during inflation via quantum tunneling. Moreover, this mass distribution can also provide an explanation to supermassive black holes formed at high redshifts.

The formation of giant planets is best studied through direct imaging by observing planets both during and after formation. Giant planets are expected to form either by core accretion, which is typically associated with low initial entropy (cold-start models) or by gravitational instability, which corresponds to a high initial entropy of the gas (hot-start models). Thus, constraining the initial entropy provides insight into the planet formation mechanism and determines the resultant brightness evolution. We find that, by observing planets in nearby moving groups of known age both through direct imaging and astrometry with Gaia, it will be possible to constrain the initial entropy of giant planets. We simulate a set of planetary systems in stars in nearby moving groups identified by BANYAN $\Sigma$ and assume a model for planet distribution consistent with radial velocity detections. We find that Gaia should be able to detect approximately 50% of planets in nearby moving groups greater than ~0.3 M$_\text{J}$. Using 5$\sigma$ contrast limits of current and future instruments, we calculate the flux uncertainty, and using models for the evolution of the planet brightness, we convert this to an initial entropy uncertainty. We find that, for future instruments such as MICADO and METIS on E-ELT and VIKiNG with VLTI, the entropy uncertainty is less than 0.5 $k_{B}$/baryon, showing that these instruments should be able to distinguish between hot and cold-start models.

Exoplanetary systems are prime targets for the Search for Extraterrestrial Intelligence (SETI). With the recent uptick in the identification of candidate and confirmed exoplanets through the work of missions like the Transiting Exoplanet Survey Satellite (TESS), we are beginning to understand that Earth-like planets are common. In this work, we extend the Breakthrough Listen (BL) search for extraterrestrial intelligence to include targeted searches of stars identified by TESS as potential exoplanet hosts. We report on 113 30-min cadence observations collected for 28 targets selected from the TESS Input Catalog (TIC) from among those identified as containing signatures of transiting planets. The targets were searched for narrowband signals from 1-11 GHz using the turboSETI pipeline architecture modified for compatibility with the Google Cloud environment. Data were searched for drift rates of +/-4 Hz/s above a minimum signal-to-noise threshold of 10, following the parameters of previous searches conducted by Price et al. (2020) and Enriquez et al. (2017). The observations presented in this work establish some of the deepest limits to date over such a wide band (1-11 GHz) for life beyond Earth. We determine that fewer than 12.72% of the observed targets possess transmitters operating at these frequencies with an Equivalent Isotropic Radiated Power greater than our derived threshold of 4.9*10^(14) W.

Double peaked light curves are observed for some Type Ic supernovae (SNe Ic) including LSQ14efd, iPTF15dtg and SN 2020bvc. One possible explanation of the first peak would be shock-cooling emission from massive extended material around the progenitor, which is produced by mass eruption or rapid expansion of the outermost layers of the progenitor shortly before the supernova explosion. We investigate the effects of such circumstellar matter (CSM) on the multi-band optical light curves of SNe Ic using the radiation hydrodynamics code STELLA. Two different SNe Ic progenitor masses at the pre-SN stage (3.93$M_\odot$ and 8.26$M_\odot$) are considered in the SN models. The adopted parameter space consists of the CSM mass of $M_\mathrm{CSM} = 0.05 - 0.3 M_\odot$, the CSM radius of $R_\mathrm{CSM} = 10^{13} - 10^{15}$cm and the explosion energy of $E_\mathrm{burst} = (1.0 - 12.0)\times10^{51}$erg. We also investigate the effects of the radioactive nickel distribution on the overall shape of the light curve and the color evolution. Comparison of our SN models with the double peaked SNe Ic LSQ14efd, iPTF15dtg and SN 2020bvc indicate that these three SNe Ic had a similar CSM structure (i.e., $M_\mathrm{CSM} \approx 0.1 - 0.2 M_\odot$ and $R_\mathrm{CSM} = 10^{13} - 10^{14}~\mathrm{cm}$), which might imply a common mechanism for the CSM formation. The implied mass loss rate of $\dot{M} \gtrsim 1.0~M_\odot~\mathrm{yr^{-1}}$ is too high to be explained by the previously suggested scenarios for pre-SN eruption, which calls for a novel mechanism.

Radio astronomy has provided unique and plentiful information about our universe. Wide--field continuum radio surveys have made vital contributions to our understanding of the formation and evolution of galaxies, clusters and AGNs (active galactic nuclei) over cosmic time. Using upgraded and proposed radio interferometers, increasingly wide and sensitive radio continuum surveys are planned, or being conducted, which require wide-field imaging methods to produce wide-survey images. In the past the dynamic range of images produced by existing methods is limited, and has often failed to achieve the dynamic range objective of the survey. To deal with this problem, a high dynamic range wide-field imaging method for radio interferometers is proposed. This method improves the widely-used W-stacking method, enabling radio astronomers to produce high--fidelity wide-field interferometric images with a dynamic range (peak:rms) exceeding 10^6:1. This method has been implemented in WSCLEAN and NIFTY by their respective authors.

Emission from the interstellar medium can be a significant contaminant of measurements of the intensity and polarization of the cosmic microwave background (CMB). For planning CMB observations, and for optimizing foreground-cleaning algorithms, a description of the statistical properties of such emission can be helpful. Here we examine a machine learning approach to inferring the statistical properties of dust from either observational data or physics-based simulations. In particular, we apply a type of neural network called a Variational Auto Encoder (VAE) to maps of the intensity of emission from interstellar dust as inferred from Planck sky maps and demonstrate its ability to a) simulate new samples with similar summary statistics as the training set, b) provide fits to emission maps withheld from the training set, and c) produce constrained realizations. We find VAEs are easier to train than another popular architecture: that of Generative Adversarial Networks (GANs), and are better-suited for use in Bayesian inference.

Using a density dependent quark mass (QMDD) model for strange quark matter, we investigate the effects of non-Newtonian gravity on the properties of strange stars and constrain the parameters of the QMDD model by employing the mass of PSR J0740+6620 and the tidal deformability of GW170817. We find that for QMDD model these mass and tidal deformability observations would rule out the existenceof str ange stars if non-Newtonian gravity effects are ignored. For the current quark masses of $m_{u0}=2.16$ MeV, $m_{d0}=4.67$ MeV, and $m_{s0}=93$ MeV, we find that a strange star can exist for values of the non-Newtonian gravity parameter $g^{2}/\mu^{2}$ in the range of 4.58 GeV$^{-2}\leq g^{2}/\mu^{2}\leq$ 9.32 GeV$^{-2}$, and that the parameters $D$ and $C$ of the QMDD modelare restricted to 158.3 MeV$\leq D^{1/2}\leq$ 181.2 MeV and $-0.65\leq C \leq -0.12$. It is found that the largest possible maximum mass of a strange star obtained with the QMDD model is $2.42 \, M_{\odot}$, and that the secondary componentof GW190814 with a mass of 2.59_{-0.09}^{+0.08}\, M_{\dot} could not be a static strange star. We also find that forthe mass and radius of PSR J0030+0451 given by Riley et al. through the analysis of observational data of NICER, there exists a very tiny allowed parameter space for which strange stars computed for the QMDD model agree with the observations of PSR J0740+6620, GW17 0817 and PSR J0030+0451 simultaneously. However, for the mass and radius given by Miller et al., no such parameter space exist.

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multi-band image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all the change points (in time domain), as well as all the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length (MDL) principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real datasets show that the proposed method enjoys very promising empirical properties. Applications to two real datasets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

The recent report of NANOGrav is gathering attention since its signal can be explained by the stochastic gravitational waves (GWs) with $\Omega_{\rm GW}\sim 10^{-9}$ at $f\sim 10^{-8}$Hz. The PBH formation scenario is one of the candidates for the NANOGrav signal, which can simultaneously explain the observed $30 M_\odot$ black holes in the binary merger events in LIGO-Virgo collaboration. We focus on the type II axion-like curvaton model of the PBH formation. In type II model the complex field whose phase part is the axion rolls down from the origin of the potential. It is found that type II model achieves the broad power spectrum of the density perturbations and can simultaneously explain the LIGO-Virgo events and the NANOGrav signal. We also improve the treatment of the non-Gaussianity of perturbations in our model to estimate the amplitude of the induced GWs precisely.

Xova is a software package that implements baseline-dependent time and channel averaging on Measurement Set data. The uv-samples along a baseline track are aggregated into a bin until a specified decorrelation tolerance is exceeded. The degree of decorrelation in the bin correspondingly determines the amount of channel and timeslot averaging that is suitable for samples in the bin. This necessarily implies that the number of channels and timeslots varies per bin and the output data loses the rectilinear input shape of the input data.

We constrain the coupling between axionlike particles (ALPs) and photons, measured with the superconducting resonant detection circuit of a cryogenic Penning trap. By searching the noise spectrum of our fixed-frequency resonant circuit for peaks caused by dark matter ALPs converting into photons in the strong magnetic field of the Penning-trap magnet, we are able to constrain the coupling of ALPs with masses around $2.7906-2.7914\,\textrm{neV/c}^2$ to $g_{a\gamma}< 1 \times 10^{-11}\,\textrm{GeV}^{-1}$. This is more than one order of magnitude lower than the best laboratory haloscope and approximately 5 times lower than the CERN axion solar telescope (CAST), setting limits in a mass and coupling range which is not constrained by astrophysical observations. Our approach can be extended to many other Penning-trap experiments and has the potential to provide broad limits in the low ALP mass range.

We perform a comprehensive demographic study of the CO extent relative to dust of the disk population in the Lupus clouds, in order to find indications of dust evolution and possible correlations with other properties. We increase up to 42 the number of disks of the region with measured CO and dust sizes ($R_{\mathrm{CO}}$, $R_{\mathrm{dust}}$) from observations with the Atacama Large Millimeter/submillimeter Array (ALMA). The sizes are obtained from modeling the ${^{12}}$CO $J = 2-1$ line emission and continuum emission at $\sim 0.89$ mm with an empirical function (Nuker profile or Gaussian function). The CO emission is more extended than the dust continuum, with a $R_{68\%}^{\mathrm{CO}}$/$R_{68\%}^{\mathrm{dust}}$ median value of 2.5, for the entire population and for a sub-sample with high completeness. 6 disks, around $15\%$ of the Lupus disk population have a size ratio above 4. Based on thermo-chemical modeling, this value can only be explained if the disk has undergone grain growth and radial drift. These disks do not have unusual properties in terms of stellar mass ($M_{\star}$), disk mass ($M_{\mathrm{disk}}$), CO and dust sizes ($R_{\mathrm{CO}}$, $R_{\mathrm{dust}}$), and mass accretion. We search for correlations between the size ratio and $M_{\star}$, $M_{\mathrm{disk}}$, $R_{\mathrm{CO}}$ and $R_{\mathrm{dust}}$: only a weak monotonic anti-correlation with the $R_{\mathrm{dust}}$ is found. The lack of strong correlations is remarkable and suggests that the bulk of the population may be in a similar evolutionary stage, independent of the stellar and disk properties. These results should be further investigated, since the optical depth difference between CO and dust continuum may play a role in the inferred size ratios. Lastly, the CO emission for the majority of the disks is consistent with optically thick emission and an average CO temperature of around 30 K.

Star formation in galaxies is a hierarchical process with a wide range of scales from smaller clusters to larger stellar complexes. Here, we present an ultra-violet imaging study of the nearby flocculent spiral galaxy NGC 7793, observed with the Ultra-Violet Imaging Telescope (UVIT). We find that the disk scale-length estimated in Far-UV (2.64$\pm$0.16 kpc) is larger than that in Near-UV (2.21$\pm$0.21 kpc) and optical (1.08 kpc), which supports the inside-out growth scenario of the galaxy disk. The star-forming UV disk is also found to be contained within the extent of H~I gas of column density greater than $10^{21}$cm$^{-2}$. With the spatial resolution of UVIT (1 pixel $\sim$ 6.8 pc), we identified 2046 young star-forming clumps in the galaxy with radii between $\sim$ 12 - 70 pc, which matches well with the size of GMCs detected in the galaxy. Around 61\% of the regions identified in our study have age younger than 20 Myr, which points to a recent enhancement of star formation across the galaxy. We also noticed that the youngest star-forming regions, with age $<$ 10 Myr, distinctly trace the flocculent arms of the galaxy. The estimated mass of the clumps cover a range between $10^3 - 10 ^6 M_{\odot}$. We noticed a gradient in the mass distribution of identified clumps along the spiral arms. We have also studied the nuclear star cluster of the galaxy and found that the stellar populations in the cluster outskirts are younger than the inner part.

Ellerman bombs are regions with enhanced Balmer line wing emission and mark magnetic reconnection in the deep solar atmosphere in active regions and quiet Sun. They are often found in regions where opposite magnetic polarities are in close proximity. Recent high resolution observations suggest that Ellerman bombs are more prevalent than thought before. We aim to determine the occurrence of Ellerman bombs in the penumbra of sunspots. We analyze high spatial resolution observations of sunspots in the Balmer H-alpha and H-beta lines as well as auxiliary continuum channels obtained with the Swedish 1-m Solar Telescope and apply the k-means clustering technique to systematically detect and characterize Ellerman Bombs. Features with all the defining characteristics of Ellerman bombs are found in large numbers over the entire penumbra. The true prevalence of these events is only fully appreciated in the H-beta line due to highest spatial resolution and lower chromospheric opacity. We find that the penumbra hosts some of the highest Ellerman bomb densities, only surpassed by the moat in the immediate surroundings of the sunspot. Some penumbral Ellerman bombs show flame morphology and rapid dynamical evolution. Many penumbral Ellerman bombs are fast moving with typical speed of 3.7 km/s and sometimes more than 10 km/s. Many penumbral Ellerman bombs migrate from the inner to the outer penumbra over hundreds of km and some continue moving beyond the outer penumbral boundary into the moat. Many penumbral Ellerman bombs are found in the vicinity of regions with opposite magnetic polarity. We conclude that reconnection is a near continuous process in the low atmosphere of the penumbra of sunspots as manifest in the form of penumbral Ellerman bombs. These are so prevalent that they may be a major sink of sunspot magnetic energy.

SN 2017ein is a narrow-lined Type Ic SN that was found to share a location with a point-like source in the face on spiral galaxy NGC 3938 in pre-supernova images, making SN 2017ein the first credible detection of a Type Ic progenitor. Results in the literature suggest this point-like source is likely a massive progenitor of 60-80 M$_{\odot}$, depending on if the source is a binary, a single star, or a compact cluster. Using new photometric and spectral data collected for 200 days, including several nebular spectra, we generate a consistent model covering the photospheric and nebular phase using a Monte Carlo radiation transport code. Photospheric phase modelling finds an ejected mass 1.2-2.0 M$_{\odot}$ with an $E_\mathrm{k}$ of $\sim(0.9 \pm0.2)\times 10^{51}$ erg, with approximately 1 M$_{\odot}$ of material below 5000 km s$^{-1}$ found from the nebular spectra. Both photospheric and nebular phase modelling suggests a $^{56}$Ni mass of 0.08-0.1 M$_{\odot}$. Modelling the [\OI] emission feature in the nebular spectra suggests the innermost ejecta is asymmetric. The modelling results favour a low mass progenitor of to 16-20 M$_{\odot}$, which is in disagreement with the pre-supernova derived high mass progenitor. This contradiction is likely due to the pre-supernova source not representing the actual progenitor.

MAVIS (MCAO-Assisted Visible Imager and Spectrograph) is an instrument proposed for the VLT Adaptive Optics Facility (AOF), which is currently in the phase-A conceptual design study. It will be the first instrument performing Multi-conjugate adaptive optics at visible wavelengths, enabling a new set of science observations. MAVIS will be installed at the Nasmyth platform of VLT UT-4 taking advantage of the already operational Adaptive Optics Facility that consists of 4 LGS and an adaptive secondary mirror with 1170 actuators. In addition, two post-focal deformable mirrors and 3 Natural Guide Stars (NGS) are foreseen for the tomographic reconstruction and correction of atmospheric turbulence. The MAVIS AO module is intended to feed both an imager and a spectrograph that will take advantage of the increased resolution and depth with respect to current instrumentation. In this paper we present the trade-off study for the optical design of the MAVIS AO module, highlighting the peculiarities of the system and the requirements imposed by AO. We propose a set of possible optical solutions able to provide a compact and efficient implementation of the different subsystems and we compare them in terms of delivered optical quality, overall throughput, encumbrance, ease of alignment and residual distortion.

We investigate the effect of strong emission line galaxies on the performance of empirical photometric redshift estimation methods. In order to artificially control the contribution of photometric error and emission lines to total flux, we develop a PCA-based stochastic mock catalogue generation technique that allows for generating infinite signal-to-noise ratio model spectra with realistic emission lines on top of theoretical stellar continua. Instead of running the computationally expensive stellar population synthesis and nebular emission codes, our algorithm generates realistic spectra with a statistical approach, and - as an alternative to attempting to constrain the priors on input model parameters - works by matching output observational parameters. Hence, it can be used to match the luminosity, colour, emission line and photometric error distribution of any photometric sample with sufficient flux-calibrated spectroscopic follow-up. We test three simple empirical photometric estimation methods and compare the results with and without photometric noise and strong emission lines. While photometric noise clearly dominates the uncertainty of photometric redshift estimates, the key findings are that emission lines play a significant role in resolving colour space degeneracies and good spectroscopic coverage of the entire colour space is necessary to achieve good results with empirical photo-z methods. Template fitting methods, on the other hand, must use a template set with sufficient variation in emission line strengths and ratios, or even better, first estimate the redshift empirically and fit the colours with templates at the best-fit redshift to calculate the K-correction and various physical parameters.

We report the 888 MHz radio detection in the Rapid ASKAP Continuum Survey (RACS) of VIK J2318$-$3113, a z=6.44 quasar. Its radio luminosity (1.2 $\times 10^{26}$ W Hz$^{-1}$ at 5 GHz) compared to the optical one (1.8 $\times 10^{24}$ W Hz$^{-1}$ at 4400 A) makes it the most distant radio-loud quasar observed so far, with a radio loudness R$\sim$70 (R$=L_{5GHz}/L_{4400A}$). Moreover, the large bolometric luminosity associated to the accretion disk (L$_{bol}$=7.5 $\times 10^{46}$ erg s$^{-1}$) suggests the presence of a supermassive black hole with a large mass ($\sim$10$^9$ M$_\odot$) when the Universe was less than a billion years old. Combining the new radio observation from RACS with archival radio data at the same frequency, we found that the flux density of the source may have varied by a factor of $\sim$2, which could suggest the presence of a relativistic jet oriented towards the line of sight, i.e. a blazar nature. However, currently available radio data do not allow us to firmly characterise the orientation of the source. Further radio and X-ray observations are needed.

Intensity maps of the 21cm emission line of neutral hydrogen are lensed by intervening large-scale structure, similar to the lensing of the cosmic microwave background temperature map. We extend previous work by calculating the lensing contribution to the full-sky 21cm bispectrum in redshift space. The lensing contribution tends to peak when equal-redshift fluctuations are lensed by a lower redshift fluctuation. At high redshift, lensing effects can become comparable to the contributions from density and redshift-space distortions.

We present a background model for dark matter searches using an array of NaI(Tl) crystals in the COSINE-100 experiment that is located in the Yangyang underground laboratory. The model includes background contributions from both internal and external sources, including cosmogenic radionuclides and surface $^{210}$Pb contamination. To improve the model in the low energy region, with the threshold lowered to 1 keV, we used a depth profile of $^{210}$Pb contamination in the surface of the NaI(Tl) crystals determined in a comparison between measured and simulated spectra. We also considered the effect of the energy scale errors propagated from the statistical uncertainties and the nonlinear detector response at low energies. The 1.7 years COSINE-100 data taken between October 21, 2016 and July 18, 2018 were used for this analysis. The Geant4 toolkit version 10.4.2 was utilized throughout the Monte Carlo simulations for the possible internal and external origins. In particular, the version provides a non-Gaussian peak around 50 keV originating from beta decays of $^{210}$Pb in a good agreement with the measured background. This improved model estimates that the activities of $^{210}$Pb and $^{3}$H are the dominant sources of the backgrounds with an average level of 2.73$\pm$0.14 counts/day/keV/kg in the energy region of 1-6 keV, using COSINE-100 data with a total exposure of 97.7 kg$\cdot$years.

The plethora of photometric data collected by the Kepler space telescope has promoted the detection of tens of thousands of stellar rotation periods. However, these periods are not found to an equal extent among different spectral types. Interestingly, early G-type stars with near-solar rotation periods are strongly underrepresented among those stars with known rotation periods. In this study we investigate whether the small number of such stars can be explained by difficulties in the period determination from photometric time series. For that purpose, we generate model light curves of early G-type stars with solar rotation periods for different inclination angles, metallicities and (magnitude-dependent) noise levels. We find that the detectability is determined by the predominant type of activity (i.e. spot or faculae domination) on the surface, which defines the degree of irregularity of the light curve, and further depends on the level of photometric noise. These two effects significantly complicate the period detection and explain the lack of solar-like stars with known near-solar rotation periods. We conclude that the rotation periods of the majority of solar-like stars with near-solar rotation periods remain undetected to date. Finally, we promote the use of new techniques to recover more periods of near-solar rotators.

We have developed an inversion procedure designed for high-resolution solar spectro-polarimeters, such as Hinode/SP or DKIST/ViSP. The procedure is based on artificial neural networks trained with profiles generated from random atmospheric stratifications for a high generalization capability. When applied to Hinode data we find a hot fine-scale network structure whose morphology changes with height. In the middle layers this network resembles what is observed in G-band filtergrams but it is not identical. Surprisingly, the temperature enhancements in the middle and upper photosphere have a reversed pattern. Hot pixels in the middle photosphere, possibly associated to small-scale magnetic elements, appear cool at the log(tau_500)=-3 and -4 level, and viceversa. Finally, we find hot arcs on the limb side of magnetic pores, which we interpret as the first direct observational evidence of the "hot wall" effect in temperature.

Large-scale magnetic field is believed to play a key role in launching and collimating jets/outflows. It was found that advection of external field by a geometrically thin disk is rather inefficient, while the external weak field may be dragged inwards by fast radially moving tenuous or/and hot gas above the thin disk. We investigate the field advection in a thin (cold) accretion disk covered with hot corona, in which turbulence is responsible for the angular momentum transfer of the gas in the disk and corona. The radial velocity of the gas in the corona is significantly higher than that in the thin disk. Our calculations show that the external magnetic flux is efficiently transported inwards by the corona, and the field line is strongly inclined towards the disk surface, which help launching outflows. The field configurations are consistent with those observed in the numerical simulations. The strength of the field is substantially enhanced in the inner region of the disk (usually several orders of magnitude higher than the external field strength), which is able to drive a fraction of gas in the corona into outflows. This mechanism may be useful in explaining the observational features in X-ray binaries and active galactic nuclei. Our results may help understanding the physics of the magneto-hydrodynamic (MHD) simulations.

We present a method to mitigate the atmospheric effects (residual atmospheric lines) in single-dish radio spectroscopy caused by the elevation difference between the target and reference positions. The method is developed as a script using the Atmospheric Transmission at Microwaves (ATM) library built into the Common Astronomy Software Applications (CASA) package. We apply the method to the data taken with the Total Power Array of the Atacama Large Millimeter/submillimeter Array. The intensities of the residual atmospheric (mostly O3) lines are suppressed by, typically, an order of magnitude for the tested cases. The parameters for the ATM model can be optimized to minimize the residual line and, for a specific O3 line at 231.28 GHz, a seasonal dependence of a best-fitting model parameter is demonstrated. The method will be provided as a task within the CASA package in the near future. The atmospheric removal method we developed can be used by any radio/millimeter/submillimeter observatory to improve the quality of its spectroscopic measurements.

We present the analysis of XMM-Newton European Photon Imaging Camera (EPIC) observations of the nova shell IPHASX J210204.7$+$471015. We detect X-ray emission from the progenitor binary star with properties that resemble those of underluminous intermediate polars such as DQ Her: an X-ray-emitting plasma with temperature of $T_\mathrm{X}=(6.4\pm3.1)\times10^{6}$ K, a non-thermal X-ray component, and an estimated X-ray luminosity of $L_\mathrm{X}=10^{30}$ erg s$^{-1}$. Time series analyses unveil the presence of two periods, the dominant with a period of $2.9\pm0.2$ hr, which might be attributed to the spin of the white dwarf, and a secondary of $4.5\pm0.6$ hr that is in line with the orbital period of the binary system derived from optical observations. We do not detect extended X-ray emission as in other nova shells probably due to its relatively old age (130-170 yr) or to its asymmetric disrupted morphology which is suggestive of explosion scenarios different to the symmetric ones assumed in available numerical simulations of nova explosions.

We have used archival infrared images obtained with the Wide Field Camera 3 on board the Hubble Space Telescope to constrain the initial mass function of low-mass stars and brown dwarfs in the W3 star-forming region. The images cover 438 arcmin$^2$, which encompasses the entire complex, and were taken in the filters F110W, F139M, and F160W. We have estimated extinctions for individual sources in these data from their colors and have dereddened their photometry accordingly. By comparing an area of the images that contains the richest concentration of previously identified W3 members to an area that has few members and is dominated by background stars, we have estimated the luminosity function for members of W3 with masses of 0.03-0.4 $M_\odot$. That luminosity function closely resembles data in typical nearby star-forming regions that have much smaller stellar populations than W3 ($\lesssim$500 vs. several thousand objects). Thus, we do not find evidence of significant variations in the initial mass function of low-mass stars and brown dwarfs with star forming conditions, which is consistent with recent studies of other distant massive star-forming regions.

We study the behaviour and properties of the solar wind using a 2.5D Alfv\'en wave driven wind model. We first systematically compare the results of an Alfv\'en wave (AW) driven wind model with a polytropic approach. Polytropic magnetohydrodynamic wind models are thermally driven, while Alfv\'en waves act as additional acceleration and heating mechanisms in the Alfv\'en wave driven model. We confirm that an AW-driven model is required to reproduce the observed bimodality of slow and fast solar winds. We are also able to reproduce the observed anti-correlation between the terminal wind velocity and the coronal source temperature with the AW-driven wind model. We also show that the wind properties along an eleven-year cycle differ significantly from one model to the other. The AW-driven model again shows the best agreement with observational data. Indeed, solar surface magnetic field topology plays an important role in the Alfv\'en wave driven wind model, as it enters directly into the input energy sources via the Poynting flux. On the other hand, the polytropic wind model is driven by an assumed pressure gradient; thus it is relatively less sensitive to the surface magnetic field topology. Finally, we note that the net torque spinning down the Sun exhibits the same trends in the two models, showing that the polytropic approach still captures correctly the essence of stellar winds.

The evidence for benzonitrile (C$_6$H$_5$CN}) in the starless cloud core TMC-1 makes high-resolution studies of other aromatic nitriles and their ring-chain derivatives especially timely. One such species is phenylpropiolonitrile (3-phenyl-2-propynenitrile, C$_6$H$_5$C$_3$N), whose spectroscopic characterization is reported here for the first time. The low resolution (0.5 cm$^{-1}$) vibrational spectrum of C$_6$H$_5$C$_3$N} has been recorded at far- and mid-infrared wavelengths (50 - 3500 cm$^{-1}$) using a Fourier Transform interferometer, allowing for the assignment of band centers of 14 fundamental vibrational bands. The pure rotational spectrum of the species has been investigated using a chirped-pulse Fourier transform microwave (FTMW) spectrometer (6 - 18 GHz), a cavity enhanced FTMW instrument (6 - 20 GHz), and a millimeter-wave one (75 - 100 GHz, 140 - 214 GHz). Through the assignment of more than 6200 lines, accurate ground state spectroscopic constants (rotational, centrifugal distortion up to octics, and nuclear quadrupole hyperfine constants) have been derived from our measurements, with a plausible prediction of the weaker bands through calculations. Interstellar searches for this highly polar species can now be undertaken with confidence since the astronomically most interesting radio lines have either been measured or can be calculated to very high accuracy below 300 GHz.

We carried out a detailed study of the temporal and broadband spectral behaviour of one of the brightest misaligned active galaxies in gamma-rays, NGC 1275 utilising 11 years of Fermi, and available Swift and AstroSat observations. Based on the cumulative flux distribution of the gamma-ray lightcurve, we identified four distinct activity states and noticed an increase in the baseline flux during the first three states. Similar nature of the increase in the average flux was also noticed in X-ray and UV bands. A large flaring activity in gamma-rays was noticed in the fourth state. The source was observed twice by AstroSat for shorter intervals (~days) during the longer observing periods (~years) state 3 and 4. During AstroSat observing periods, the source gamma-ray flux was higher than the average flux observed during longer duration states. The increase in the average baseline flux from state 1 to state 3 can be explained considering a corresponding increase of jet particle normalisation. The inverse Comptonisation of synchrotron photons explained the average X-ray and gamma-ray emission by jet electrons during the first three longer duration states. However, during the shorter duration AstroSat observing periods, a shift of the synchrotron peak frequency was noticed, and the synchrotron emission of jet electrons well explained the observed X-ray flux.

The 21 cm linear polarization due to Thomson scattering off free electrons can probe the distribution of neutral hydrogen in the intergalactic medium during the epoch of reionization, complementary to the 21 cm temperature fluctuations. Previous study (Babich & Loeb 2005) estimated the strength of polarization with a toy model and claimed that it can be detected with 1-month observation of the Square Kilometre Array (SKA). Here we revisit this investigation with account of nonlinear terms due to inhomogeneous reionization, using seminumerical reionization simulations to provide the realistic estimation of the 21 cm TE and EE angular power spectra ($C^{\rm TE}_\ell$ and $C^{\rm EE}_\ell$). We find that (1) both power spectra are enhanced on sub-bubble scales but suppressed on super-bubble scales, compared with previous results; (2) $C^{\rm TE}_\ell$ flips its sign at $\ell \simeq 50$, which can probe the H I bias at large scales; (3) the ratios of the power spectrum to its maximum value during reionization at a given $\ell$, i.e. $C^{\rm TE}_\ell / C^{\rm TE}_{\ell,{\rm max}} $ and $C^{\rm EE}_{\ell}/C^{\rm EE}_{\ell,{\rm max}}$, show robust correlations with the global ionized fraction. However, measurement of this signal will be very challenging not only because the overall strength is weaker than the sensitivity of SKA, but also because Faraday rotation due to Galactic and extragalactic magnetic fields significantly modifies the observed polarization. Nevertheless, it may still be possible that the 21 cm linear polarization signal may be detected through other approaches, e.g. its cross-correlation with other probes.

We present a concept for a machine-learning classification of hard X-ray (HXR) emissions from solar flares observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), identifying flares that are either occulted by the solar limb or located on the solar disk. Although HXR observations of occulted flares are important for particle-acceleration studies, HXR data analyses for past observations were time consuming and required specialized expertise. Machine-learning techniques are promising for this situation, and we constructed a sample model to demonstrate the concept using a deep-learning technique. Input data to the model are HXR spectrograms that are easily produced from RHESSI data. The model can detect occulted flares without the need for image reconstruction nor for visual inspection by experts. A technique of convolutional neural networks was used in this model by regarding the input data as images. Our model achieved a classification accuracy better than 90 %, and the ability for the application of the method to either event screening or for an event alert for occulted flares was successfully demonstrated.

The recently commissioned Dark Energy Spectroscopic Instrument (DESI) will measure the expansion historyof the universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasarsover 14000 sq deg will be measured during the life of the experiment. A new prime focus corrector for theKPNO Mayall telescope delivers light to 5000 fiber optic positioners. The fibers in turn feed ten broad-bandspectrographs. We describe key aspects and lessons learned from the development, delivery and installation ofthe fiber system at the Mayall telescope.

FRAM (F/Photometric Robotic Atmospheric Monitor) is a robotic telescope operated at the Pierre Auger Observatory in Argentina for the purposes of atmospheric monitoring using stellar photometry. As a passive system which does not produce any light that could interfere with the observations of the fluorescence telescopes of the observatory, it complements the active monitoring systems that use lasers. We discuss the applications of stellar photometry for atmospheric monitoring at optical observatories in general and the particular modes of operation employed by the Auger FRAM. We describe in detail the technical aspects of FRAM, the hardware and software requirements for a successful operation of a robotic telescope for such a purpose and their implementation within the FRAM system.

Oscillons are spatially localized structures that appear in scalar field theories and exhibit extremely long life-times. We go beyond single-field analyses and study oscillons comprised of multiple interacting fields, each having an identical potential with quadratic, quartic and sextic terms. We consider quartic interaction terms of either attractive or repulsive nature. In the two-field case, we construct semi-analytical oscillon profiles for different values of the potential parameters and coupling strength using the two-timing small-amplitude formalism. We show that the interaction sign, attractive or repulsive, leads to different oscillon solutions, albeit with similar characteristics, like the emergence of "flat-top" shapes. In the case of attractive interactions, the oscillons can reach higher values of the energy density and smaller values of the width. For repulsive interactions we identify a threshold for the coupling strength, above which oscillons do not exist within the two-timing small-amplitude framework. We extend the Vakhitov-Kolokolov (V-K) stability criterion, which has been used to study single-field oscillons, and show that the symmetry of the potential leads to similar equations as in the single-field case, albeit with modified terms. We explore the basin of attraction of stable oscillon solutions numerically to test the validity of the V-K criterion and show that, depending on the initial perturbation size, unstable oscillons can either completely disperse or relax to the closest stable configuration. Similarly to the V-K criterion, the decay rate and lifetime of two-field oscillons are found to be qualitatively and quantitatively similar to their single-field counterparts. Finally, we generalize our analysis to multi-field oscillons and show that the governing equations for their shape and stability can be mapped to the ones arising in the two-field case.

In traditional models only an order one fraction of energy is transferred from the inflaton to radiation through nonperturbative resonance production in preheating immediately after inflation, due to backreaction effects. We propose a particle production mechanism that could improve the depletion of the inflaton energy density by up to four orders of magnitude. The improvement comes from the fast perturbative decays of resonantly produced daughter particles. They act as a "spillway" to drain these daughter particles, reducing their backreaction on the inflaton and keeping the resonant production effective for a longer period. Thus we dub the scenario "spillway preheating". We also show that the fraction of energy density remaining in the inflaton has a simple inverse power-law scaling in the scenario. In general, spillway preheating is a much more efficient energy dissipation mechanism, which may have other applications in model building for particle physics.

The radiation of linear momentum imparts a recoil (or "kick") to the center of mass of a merging black hole binary system. Recent numerical relativity calculations have shown that eccentricity can lead to an approximate 25% increase in recoil velocities for equal-mass, spinning binaries with spins lying in the orbital plane ("superkick" configurations) [PRD 10 (2020) 024044 (arXiv:1910.01598)]. Here we investigate the impact of nonzero eccentricity on the kick magnitude and gravitational-wave emission of nonspinning, unequal-mass black hole binaries. We confirm that nonzero eccentricities at merger can lead to kicks which are larger by up to ~25% relative to the quasicircular case. We also find that the kick velocity $v$ has an oscillatory dependence on eccentricity, that we interpret as a consequence of changes in the angle between the infall direction at merger and the apoapsis (or periapsis) direction.

In this work, we explore the connection between the critical curves ("shadows") and the quasinormal mode frequencies (in the eikonal limit) of Kerr black holes. This mapping has been previously established for non-rotating black holes. We show that, the shadow seen by an distant observer at a given inclination angle, can be mapped to a family of quasinormal modes with $m/(\ell+1/2)$ bounded within certain range, where $m$ is is azimuthal node number and $\ell$ is the angular node number. We discuss the possibility of testing such relation with space-borne gravitational wave detectors and the next-generation Event Horizon Telescope.

We perform a new test of general relativity (GR) with signals from GWTC-2, the LIGO and Virgo catalog of gravitational wave detections. We search for the presence of amplitude birefringence, in which left versus right circularly polarized modes of gravitational waves are exponentially enhanced and suppressed during propagation. Such an effect is present in various beyond-GR theories but is absent in GR. We constrain the amount of amplitude birefringence consistent with the data through an opacity parameter $\kappa$, which we bound to be $\kappa \lesssim 0.74 \textrm{ Gpc}^{-1}$. We then use these theory-agnostic results to constrain Chern-Simons gravity, a beyond-GR theory with motivations in quantum gravity. We bound the canonical Chern-Simons lengthscale to be $\ell_0 \lesssim 1.0 \times 10^3$ km, improving on previous long-distance measurement results by a factor of two.

We present a scenario where an axion-like field drives inflation until a potential barrier, which keeps a waterfall field at the origin, disappears and a waterfall transition occurs. Such a barrier separates the scale of inflation from that of the waterfall transition. We find the observed spectrum of the cosmic microwave background indicates that the decay constant of the inflaton is well below the Planck scale, with the inflationary Hubble parameter spanning a wide range. Further, our model involves dark matter candidates including the inflaton itself. Also, for a complex waterfall field, we can determine cosmologically the Peccei-Quinn scale associated with the strong CP problem.

Being generated, the relic neutrino background contained equal fractions of electron $\nu_e$, muon $\nu_\mu$, and taon $\nu_\tau$ neutrinos. We show that the gravitational field of our Galaxy and other nearby cosmic objects changes this composition near the Solar System, enriching it with the heaviest neutrino $nu_3$. This mass state is almost free of the electron component (only $\sim 2\%$ of $\nu_e$) and contains more muon component than the tau one. As a result, the relic background becomes enriched with taon and particularly muon neutrinos. The electron relic neutrinos are the rarest for a terrestrial observer: instead of $1/3$, the relic background may contain only $\gtrsim 20\%$ of them.

The KArlsruhe TRItium Neutrino (KATRIN) experiment aims to measure the neutrino mass with a sensitivity of $0.2\,eV$ ($90\,\%$ CL). This will be achieved by a precision measurement of the endpoint region of the $\beta$-electron spectrum of tritium decay. The electrons from tritium $\beta$-decay are produced in the Windowless Gaseous Tritium Source (WGTS) and guided magnetically through the beamline. In order to accurately extract the neutrino mass the source properties, in particular the activity, are required to be stable and known to a high precision. The WGTS therefore undergoes constant extensive monitoring from several measurement systems. The Forward Beam Monitor (FBM) is one such monitoring system. The FBM system comprises a complex mechanical setup capable of inserting a detector board into the KATRIN beamline inside the Cryogenic Pumping Section with a positioning precision of better than $0.3\,mm$. The electron flux density at that position is on the order of $10^{6}\,s^{-1}mm^{-2}$. The detector board contains a hall sensor, a temperature gauge, and two silicon detector chips of $\textit{p}$-$\textit{i}$-$\textit{n}$ diode type which can measure the $\beta$-electron flux from the source with a precision of $0.1\,\%$ in less than a minute with an energy resolution of FWHM = $2\,keV$.

We present an introduction to some of the state of the art in reduced order and surrogate modeling in gravitational wave (GW) science. Approaches that we cover include Principal Component Analysis, Proper Orthogonal Decomposition, the Reduced Basis approach, the Empirical Interpolation Method, Reduced Order Quadratures, and Compressed Likelihood evaluations. We divide the review into three parts: representation/compression of known data, predictive models, and data analysis. The targeted audience is that one of practitioners in GW science, a field in which building predictive models and data analysis tools that are both accurate and fast to evaluate, especially when dealing with large amounts of data and intensive computations, are necessary yet can be challenging. As such, practical presentations and, sometimes, heuristic approaches are here preferred over rigor when the latter is not available. This review aims to be self-contained, within reasonable page limits, with little previous knowledge (at the undergraduate level) requirements in mathematics, scientific computing, and other disciplines. Emphasis is placed on optimality, as well as the curse of dimensionality and approaches that might have the promise of beating it. We also review most of the state of the art of GW surrogates. Some numerical algorithms, conditioning details, scalability, parallelization and other practical points are discussed. The approaches presented are to large extent non-intrusive and data-driven and can therefore be applicable to other disciplines. We close with open challenges in high dimension surrogates, which are not unique to GW science.