Papers for Friday, Jan 29 2021

The delay-time distribution (DTD) is the occurrence rate of a class of objects as a function of time after a hypothetical burst of star formation. DTDs are mainly used as a statistical test of stellar evolution scenarios for supernova progenitors, but they can be applied to many other classes of astronomical objects. We calculate the first DTD for RR Lyrae variables using 29,810 RR Lyrae from the OGLE-IV survey and a map of the stellar-age distribution (SAD) in the Large Magellanic Cloud (LMC). We find that $\sim 46\%$ of the OGLE-IV RR Lyrae are associated with delay-times older than 8 Gyr (main-sequence progenitor masses less than 1 M$_{\odot}$), and consistent with existing constraints on their ages, but surprisingly about $51\%$ of RR Lyrae appear have delay times $1.2-8$ Gyr (main-sequence masses between $1 - 2$ M$_{\odot}$ at LMC metallicity). This intermediate-age signal also persists outside the Bar-region where crowding is less of a concern, and we verified that without this signal, the spatial distribution of the OGLE-IV RR Lyrae is inconsistent with the SAD map of the LMC. Since an intermediate-age RR Lyrae channel is in tension with the lack of RR Lyrae in intermediate-age clusters (noting issues with small-number statistics), and the age-metallicity constraints of LMC stars, our DTD result possibly indicates that systematic uncertainties may still exist in SAD measurements of old-stellar populations, perhaps stemming from the construction methodology or the stellar evolution models used. We described tests to further investigate this issue.

We study the phenomenon of gravitational particle production as applied to a scalar spectator field in the context of $\alpha$-attractor inflation. Assuming that the scalar has a minimal coupling to gravity, we calculate the abundance of gravitationally-produced particles as a function of the spectator's mass $m_\chi$ and the inflaton's $\alpha$ parameter. If the spectator is stable and sufficiently weakly coupled, such that it does not thermalize after reheating, then a population of spin-0 particles is predicted to survive in the universe today, providing a candidate for dark matter. Inhomogeneities in the spatial distribution of dark matter correspond to an isocurvature component, which can be probed by measurements of the cosmic microwave background anisotropies. We calculate the dark matter-photon isocurvature power spectrum and by comparing with upper limits from Planck, we infer constraints on $m_\chi$ and $\alpha$. If the scalar spectator makes up all of the dark matter today, then for $\alpha = 10$ and $T_\text{RH} = 10^4 \ \mathrm{GeV}$ we obtain $m_\chi > 1.8 \times 10^{13} \ \mathrm{GeV} \approx 1.2 \, m_\phi$, where $m_\phi$ is the inflaton's mass.

The environments of explosive transients link their progenitors to the underlying stellar population, providing critical clues for their origins. However, some Ca-rich supernovae (SNe) and short gamma ray burst (sGRBs) appear to be located at remote locations, far from the stellar population of their host galaxy, challenging our understanding of their origin and/or physical evolution. These findings instigated models suggesting that either large velocity kicks were imparted to their progenitors, allowing them to propagate to large distances and attain their remote locations; or that they formed in dense globular clusters residing in the halos. Here we show that instead, the large spatial-offsets of these transients are naturally explained by the observations of highly extended underlying stellar populations in (mostly early type) galaxy halos, typically missed since they can only be identified through ultra-deep/stacked images. Consequently, no large velocity kicks, nor halo globular cluster environments are required in order to explain the origin of these transients. These findings support thermonuclear explosions on white-dwarfs, for the origins of Ca-rich SNe progenitors, and no/small-natal-kicks given to sGRB progenitors. Since early-type galaxies contain older stellar populations, transient arising from older stellar populations would have larger fractions of early-type galaxy hosts, and consequently larger fractions of transients at large offsets. This is verified by our results for sGRBs and Ca-rich SNe showing different offset distributions in early vs. late-type galaxies. Furthermore, once divided between early and late type galaxies, the offsets' distributions of the different transients are consistent with each other. Finally, we point out that studies of other transients' offset distribution (e.g. Ia-SNe or FRBs) should similarly consider the host galaxy-type.

We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization and the growth of PMFs. We run a series of $(5{\rm Mpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field and its spectral index. We find that PMFs with a spectral index ($n_B$) and a comoving amplitude ($B_0$) that have $n_B>-0.562\log_{10}(B_0/1{\rm n}G) - 3.35$ produce electron optical depths ($\tau_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Ly$\alpha$ forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a particularly sensitive indicator of the properties of the PMFs, even for those that satisfy our $\tau_e$ constraint. Furthermore, strong PMFs can increase the escape of LyC photons by up to 25% and shrink the effective radii of galaxies by 44% which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500{\rm A}}=-13}$ can probe the sources that provide the 50% of photons for reionization out to $z=12$.

We present the spatial analysis of five Compton thick (CT) active galactic nuclei (AGNs), including MKN 573, NGC 1386, NGC 3393, NGC 5643, and NGC 7212, for which high resolution Chandra observations are available. For each source, we find hard X-ray emission (>3 keV) extending to ~kpc scales along the ionization cone, and for some sources, in the cross-cone region. This collection represents the first, high-signal sample of CT AGN with extended hard X-ray emission for which we can begin to build a more complete picture of this new population of AGN. We investigate the energy dependence of the extended X-ray emission, including possible dependencies on host galaxy and AGN properties, and find a correlation between the excess emission and obscuration, suggesting a connection between the nuclear obscuring material and the galactic molecular clouds. Furthermore, we find that the soft X-ray emission extends farther than the hard X-rays along the ionization cone, which may be explained by a galactocentric radial dependence on the density of molecular clouds due to the orientation of the ionization cone with respect to the galactic disk. These results are consistent with other CT AGN with observed extended hard X-ray emission (e.g., ESO 428-G014 and the Ma et al. 2020 CT AGN sample), further demonstrating the ubiquity of extended hard X-ray emission in CT AGN.

The matter density field at $z\sim 6$ is very challenging to probe. One of the traditional probes of the low density IGM that works successfully at lower redshifts is the Lyman-alpha forest in quasar spectra. However, at the end of reionization, the residual neutral hydrogen usually creates saturated absorption, thus much of the information about the gas density is lost. Luckily, in a quasar proximity zone, the ionizing radiation is exceptionally intense, thus creating a large region with non-zero transmitted flux. In this study we use the synthetic spectra from simulations to investigate how to recover the density fluctuations inside the quasar proximity zones. We show that, under ideal conditions, the density can be recovered accurately with a small scatter. We also discuss how systematics such as the quasar continuum fitting and reionization models affect the results. This study shows that by analyzing the absorption features inside quasar proximity zones we can potentially constrain quasar properties and the environments they reside in.

Starlight subtraction algorithms based on the method of Karhunen-Lo\`eve eigenimages have proved invaluable to exoplanet direct imaging. However, they scale poorly in runtime when paired with differential imaging techniques. In such observations, reference frames and frames to be starlight-subtracted are drawn from the same set of data, requiring a new subset of references (and eigenimages) for each frame processed to avoid self-subtraction of the signal of interest. The data rates of extreme adaptive optics instruments are such that the only way to make this computationally feasible has been to downsample the data. We develop a technique that updates a pre-computed singular value decomposition of the full dataset to remove frames (i.e. a "downdate") without a full recomputation, yielding the modified eigenimages. This not only enables analysis of much larger data volumes in the same amount of time, but also exhibits near-linear scaling in runtime as the number of observations increases. We apply this technique to archival data and investigate its scaling behavior for very large numbers of frames $N$. The resulting algorithm provides speed improvements of $2.6\times$ (for 200 eigenimages at $N = 300$) to $140 \times$ (at $N = 10^4$) with the advantage only increasing as $N$ grows. This algorithm has allowed us to substantially accelerate KLIP even for modest $N$, and will let us quickly explore how KLIP parameters affect exoplanet characterization in large $N$ datasets.

We report the temporal and spectral analysis of three thermonuclear X-ray bursts from 4U 1608-52, observed by the Neutron Star Interior Composition Explorer (NICER) during and just after the outburst observed from the source in 2020. In two of the X-ray bursts, we detect secondary peaks, 30 and 18 seconds after the initial peaks. The secondary peaks show a fast rise exponential decay-like shape resembling a thermonuclear X-ray burst. Time-resolved X-ray spectral analysis reveals that the peak flux, blackbody temperature, and apparent emitting radius values of the initial peaks are in agreement with X-ray bursts previously observed from 4U 1608-52, while the same values for the secondary peaks tend toward the lower end of the distribution of bursts observed from this source. The third X-ray burst, which happened during much lower accretion rates did not show any evidence for a deviation from an exponential decay and was significantly brighter than the previous bursts. We present the properties of the secondary peaks and discuss the events within the framework of short recurrence time bursts or bursts with secondary peaks. We find that the current observations do not fit in standard scenarios and challenge our understanding of flame spreading.

The Gaia mission is delivering exquisite astrometric data for 1.47 billion sources, which are revolutionizing many fields in astronomy. For a small fraction of these sources the astrometric solutions are poor, and the reported values and uncertainties may not apply. For many analyses it is important to recognize and excise these spurious results, commonly done by means of quality flags in the Gaia catalog. Here we devise and apply a path to separating 'good' from 'bad' astrometric solutions that is an order-of-magnitude cleaner than any single flag: we achieve a purity of 99.7% and a completeness of 97.6% as validated on our test data. We devise an extensive sample of manifestly bad astrometric solutions: sources whose inferred parallax is 'negative' at >= 4.5 sigma; and a corresponding sample of presumably good solutions: the sources in HEALPix patches of the sky that do not contain extremely negative parallaxes. We then train a neural net that uses 14 pertinent Gaia catalog entries to discriminate these two samples, captured in a single 'astrometric fidelity' parameter. An extensive and diverse set of verification tests show that our approach to assessing astrometric fidelity works very cleanly also in the regime where no negative parallaxes are involved; its main limitations are in the very low S/N regime. Our astrometric fidelities for all EDR3 can be queried via the Virtual Observatory. In the spirit of open science, we make our code and training/validation data public, so that our results can be easily reproduced.

Hydrodynamical (or 3D) and non-local thermodynamic equilibrium (NLTE) effects are known to affect abundance analyses. However, there are very few observational abundance testsof 3D and NLTE models. We developed a new way of testing the abundance predictions of 3D and NLTE models, taking advantage of large spectroscopic survey data. We use a line-by-line analysis of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) spectra (H band) with the Brussels Automatic Code for Characterizing High accUracy Spectra (BACCHUS). We compute line-by-line abundances of Mg, Si, Ca, and Fe for a large number of globular cluster K giants in the APOGEE survey. We compare this line-by-line analysis against NLTE and 3D predictions. While the 1D-NLTE models provide corrections in the right direction, there are quantitative discrepancies between different models. We observe a better agreement with the data for the models including reliable collisional cross-sections. The agreement between data and models is not always satisfactory when the 3D spectra are computed in LTE. However, we note that for a fair comparison, 3D corrections should be computed with self-consistently derived stellar parameters, and not on 1D models with identical stellar parameters. Finally, we focus on 3D and NLTE effects on Fe lines in the H band, where we observe a systematic difference in abundance relative to the value from the optical. Our results suggest that the metallicities obtained from the H band are more accurate in metal-poor giants. More atomic data and extended self-consistent 3D-NLTE computations need to be made. The method we have developed for testing 3D and NLTE models could be extended to other lines and elements, and is particularly suited for large spectroscopic surveys.

The application of Bayesian techniques to astronomical data is generally non-trivial because the fitting parameters can be strongly degenerated and the formal uncertainties are themselves uncertain. An example is provided by the contradictory claims over the presence or absence of a universal acceleration scale (g$_\dagger$) in galaxies based on Bayesian fits to rotation curves. To illustrate the situation, we present an analysis in which the Newtonian gravitational constant $G_N$ is allowed to vary from galaxy to galaxy when fitting rotation curves from the SPARC database, in analogy to $g_{\dagger}$ in the recently debated Bayesian analyses. When imposing flat priors on $G_N$, we obtain a wide distribution of $G_N$ which, taken at face value, would rule out $G_N$ as a universal constant with high statistical confidence. However, imposing an empirically motivated log-normal prior returns a virtually constant $G_N$ with no sacrifice in fit quality. This implies that the inference of a variable $G_N$ (or g$_{\dagger}$) is the result of the combined effect of parameter degeneracies and unavoidable uncertainties in the error model. When these effects are taken into account, the SPARC data are consistent with a constant $G_{\rm N}$ (and constant $g_\dagger$).

The actual mechanism(s) powering galactic outflows in active galactic nuclei (AGN) is still a matter of debate. At least two physical models have been considered in the literature: wind shocks and radiation pressure on dust. Here we provide a first quantitative comparison of the AGN radiative feedback scenario with observations of galactic outflows. We directly compare our radiation pressure-driven shell models with the observational data from the most recent compilation of molecular outflows on galactic scales. We show that the observed dynamics and energetics of galactic outflows can be reproduced by AGN radiative feedback, with the inclusion of radiation trapping and/or luminosity evolution. The predicted scalings of the outflow energetics with AGN luminosity can also quantitatively account for the observational scaling relations. Furthermore, sources with both ultra-fast and molecular outflow detections are found to be located in the `forbidden' region of the $N_\mathrm{H} - \lambda$ plane. Overall, an encouraging agreement is obtained over a wide range of AGN and host galaxy parameters. We discuss our results in the context of recent observational findings and numerical simulations. In conclusion, AGN radiative feedback is a promising mechanism for driving galactic outflows that should be considered, alongside wind feedback, in the interpretation of future observational data.

We present ALMA observations of CO isotopologues and high-density molecular tracers (HCN, HCO+, CN, etc.) in NGC 7465, an unusually gas-rich early-type galaxy that acquired its cold gas recently. In the inner 300 pc, the molecular gas kinematics are misaligned with respect to all other galaxy components; as the gas works its way inward it is torqued into polar orbits about the stellar kinematically-decoupled core (KDC), indicating that the stellar KDC is not related to the current gas accretion event. The galaxy also exhibits unusually high 12CO/13CO line ratios in its nucleus but typical 13CO/C18O ratios. Our calculations show that this result does not necessarily indicate an unusual [12CO/13CO] abundance ratio but rather that 12CO (1-0) is optically thin due to high temperatures and/or large linewidths associated with the inner decoupled, misaligned molecular structure. Line ratios of the higher-density tracers suggest that the densest phase of molecular gas in NGC 7465 has a lower density than is typical for nearby galaxies, possibly as a result of the recent gas accretion. All of the observed molecular properties of NGC 7465 are consistent with it having acquired its molecular (and atomic) gas from a spiral galaxy. Further detailed studies of the CO isotopologues in other early-type galaxies would be valuable for investigating the histories of those that may have acquired their gas from dwarfs. Finally, these ALMA data also show an unidentified line source that is probably a background galaxy similar to those found at z=1-3 in blind CO surveys.

We estimate the black hole spin parameter in GRS 1915+105 using the continuum-fitting method with revised mass and inclination constraints based on the very long baseline interferometric parallax measurement of the distance to this source. We fit Rossi X-ray Timing Explorer observations selected to be accretion disk-dominated spectral states as described in McClinotck et al. (2006) and Middleton et al. (2006), which previously gave discrepant spin estimates with this method. We find that, using the new system parameters, the spin in both datasets increased, providing a best-fit spin of $a_*=0.86$ for the Middleton et al. data and a poor fit for the McClintock et al. dataset, which becomes pegged at the BHSPEC model limit of $a_*=0.99$. We explore the impact of the uncertainties in the system parameters, showing that the best-fit spin ranges from $a_*= 0.4$ to 0.99 for the Middleton et al. dataset and allows reasonable fits to the McClintock et al. dataset with near maximal spin for system distances greater than $\sim 10$ kpc. We discuss the uncertainties and implications of these estimates.

In search for a connection between the formation of stars and the formation of planets, a new semi-automatic spectral analysis method using \textsf{iSpec} was developed for the TIGRE telescope installed in Guanajuato, Mexico. TIGRE is a 1.2m robotic telescope, equipped with an Echelle spectrograph (HEROS), with a resolution R $\simeq 20000$. \textsf{iSpec} is a synthetic spectral fitting program for stars that allows to determine in an homogeneous way their fundamental parameters: effective temperature, $T_{\rm eff}$, surface gravity, $\log g$, metallicities, $[M/H]$ and $[Fe/H]$, and rotational velocity, $V \sin i$. In this first article we test our method by analysing the spectra of 46 stars, host of exoplanets, obtained with the TIGRE.

12 stable single sunspots were investigated in two aspects. The outer boundaries of the sunspots umbra were determined by an independent mathematical method and the vertical component of the magnetic field at these boundaries was found. We analyzed data from the SDO station using the segments of the continuum and magnetic field. It is shown that in a wide range of sunspot fields, the field average along the contour changes weakly, while the vertical component of the magnetic field itself, in a particular spot, changes significantly along the defined contour which reflects the fibrous structure of the penumbral field.

A sample of 46 stars, host of exoplanets, is used to search for a connection between their formation process and the formation of the planets rotating around them. Separating our sample in two, stars hosting high-mass exoplanets (HMEs) and low-mass exoplanets (LMEs), we found the former to be more massive and to rotate faster than the latter. We also found the HMEs to have higher orbital angular momentum than the LMEs and to have lost more angular momentum through migration. These results are consistent with the view that the more massive the star and higher its rotation, the more massive was its protoplanetarys disk and rotation, and the more efficient the extraction of angular momentum from the planets.

$\tau^{9}$ Eri is a Bp star that was previously reported to be a single-lined spectroscopic binary. Using 17 ESPaDOnS spectropolarimetric (Stokes $V$) observations we identified the weak spectral lines of the secondary component and detected a strong magnetic field in the primary. We performed orbital analysis of the radial velocities of both components to find a slightly eccentric orbit ($e= 0.129$) with a period of $5.95382(2)$ days. The longitudinal magnetic field ($B_\ell$) of the primary was measured from each of the Stokes $V$ profiles, with typical error bars smaller than 10 G. Equivalent widths (EWs) of LSD profiles corresponding to only the Fe lines were also measured. We performed frequency analysis of both the $B_\ell$ and EW measurements, as well as of the Hipparcos, SMEI, and TESS photometric data. All sets of photometric observations produce two clear, strong candidates for the rotation period of the Bp star: 1.21 days and 3.82 days. The $B_\ell$ and EW measurements are consistent with only the 3.82-day period. We conclude that HD 25267 consists of a late-type Bp star (M= $3.6_{-0.2}^{+0.1} M_\odot$, T= $12580_{-120}^{+150}$ K) with a rotation period of 3.82262(4) days orbiting with a period of 5.95382(2) days with a late-A/early-F type secondary companion (M= $1.6\pm 0.1 M_\odot$, T= $7530_{-510}^{+580}$ K). The Bp star's magnetic field is approximately dipolar with $i= 41\pm 2^{\circ}$, $\beta= 158\pm 5^{\circ}$ and $B_{\rm d}= 1040\pm 50$ G. All evidence points to the strong $1.209912(3)$ day period detected in photometry, along with several other weaker photometric signals, as arising from $g$-mode pulsations in the primary.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne far-infrared telescope that will survey star formation history over cosmological time scales to improve our understanding of why the star formation rate declined at redshift z < 2, despite continued clustering of dark matter. Specifically,EXCLAIM will map the emission of redshifted carbon monoxide and singly-ionized carbon lines in windows over a redshift range 0 < z < 3.5, following an innovative approach known as intensity mapping. Intensity mapping measures the statistics of brightness fluctuations of cumulative line emissions instead of detecting individual galaxies, thus enabling a blind, complete census of the emitting gas. To detect this emission unambiguously, EXCLAIM will cross-correlate with a spectroscopic galaxy catalog. The EXCLAIM mission uses a cryogenic design to cool the telescope optics to approximately 1.7 K. The telescope features a 90-cm primary mirror to probe spatial scales on the sky from the linear regime up to shot noise-dominated scales. The telescope optical elements couple to six {\mu}-Spec spectrometer modules, operating over a 420-540 GHz frequency band with a spectral resolution of 512 and featuring microwave kinetic inductance detectors. A Radio Frequency System-on-Chip (RFSoC) reads out the detectors in the baseline design. The cryogenic telescope and the sensitive detectors allow EXCLAIM to reach high sensitivity in spectral windows of low emission in the upper atmosphere. Here, an overview of the mission design and development status since the start of the EXCLAIM project in early 2019 is presented.

We review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average inter-particle separation in a galaxy like the Milky Way, and are well described as classical waves. We outline the particle physics motivations for them, including the QCD axion and ultra-light axion-like-particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: (1) Wave interference leads to order unity density fluctuations on de Broglie scale. A manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. (2) For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and random walks, another manifestation of wave interference. The halo/subhalo abundance is suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. (3) For ultra-light ~$10^{-22}$ eV dark matter, the wave interference substructures can be probed by tidal streams/gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius/image separation. (4) Axion detection experiments are sensitive to interference substructures for moderately light masses. The stochastic nature of the waves affects the interpretation of experiments and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological/galactic/black-hole observations, are discussed.

The Dark Energy Spectroscopic Instrument (DESI) is an ongoing spectroscopic survey to measure the dark energy equation of state to unprecedented precision. We describe the DESI Sky Continuum Monitor System, which tracks the night sky brightness as part of a system that dynamically adjusts the spectroscopic exposure time to produce more uniform data quality and to maximize observing efficiency. The DESI dynamic exposure time calculator (ETC) will combine sky brightness measurements from the Sky Monitor with data from the guider system to calculate the exposure time to achieve uniform signal-to-noise ratio (SNR) in the spectra under various observing conditions. The DESI design includes 20 sky fibers, and these are split between two identical Sky Monitor units to provide redundancy. Each Sky Monitor unit uses an SBIG STXL-6303e CCD camera and supports an eight-position filter wheel. Both units have been completed and delivered to the Mayall Telescope at the Kitt Peak National Observatory. Commissioning results show that the Sky Monitor delivers the required performance necessary for the ETC.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will map carbon monoxide and singly-ionized carbon emission lines across redshifts from 0 to 3.5, using an intensity mapping approach. EXCLAIM will broaden our understanding of these elemental and molecular gases and the role they play in star formation processes across cosmic time scales. The focal plane of EXCLAIM's cryogenic telescope features six {\mu}-Spec spectrometers. {\mu}-Spec is a compact, integrated grating-analog spectrometer, which uses meandered superconducting niobium microstrip transmission lines on a single-crystal silicon dielectric to synthesize the grating. It features superconducting aluminum microwave kinetic inductance detectors (MKIDs), also in a microstrip architecture. The spectrometers for EXCLAIM couple to the telescope optics via a hybrid planar antenna coupled to a silicon lenslet. The spectrometers operate from 420 to 540 GHz with a resolving power R={\lambda}/{\Delta}{\lambda}=512 and employ an array of 355 MKIDs on each spectrometer. The spectrometer design targets a noise equivalent power (NEP) of 2x10-18W/\sqrt{Hz} (defined at the input to the main lobe of the spectrometer lenslet beam, within a 9-degree half width), enabled by the cryogenic telescope environment, the sensitive MKID detectors, and the low dielectric loss of single-crystal silicon. We report on these spectrometers under development for EXCLAIM, providing an overview of the spectrometer and component designs, the spectrometer fabrication process, fabrication developments since previous prototype demonstrations, and the current status of their development for the EXCLAIM mission.

Complex A is a high-velocity cloud that is traversing through the Galactic halo toward the Milky Way's disk. We combine both new and archival Green Bank Telescope observations to construct a spectroscopically resolved HI~21-cm map of this entire complex at a $17.1\lesssim\log{\left({N_{\rm HI},\,1\sigma}/{\rm cm}^{-2}\right)}\lesssim17.9$ sensitivity for a ${\rm FWHM}=20~{\rm km}\,{\rm s}^{-1}$ line and $\Delta\theta=9.1\,{\rm arcmins}$ or $17\lesssim\Delta d_{\theta}\lesssim30~\rm pc$ spatial resolution. We find that that Complex A is has a Galactic standard of rest frame velocity gradient of $\Delta\rm v_{GSR}/\Delta L=25~{\rm km}\,{\rm s}^{-1}/{\rm kpc}$ along its length, that it is decelerating at a rate of $\langle a\rangle_{\rm GSR}=55~{\rm km}/{\rm yr}^2$, and that it will reach the Galactic plane in $\Delta t\lesssim70~{\rm Myrs}$ if it can survive the journey. We have identify numerous signatures of gas disruption. The elongated and multi-core structure of Complex A indicates that either thermodynamic instabilities or shock-cascade processes have fragmented this stream. We find Rayleigh-Taylor fingers on the low-latitude edge of this HVC; many have been pushed backward by ram-pressure stripping. On the high-latitude side of the complex, Kelvin-Helmholtz instabilities have generated two large wings that extend tangentially off Complex A. The tips of these wings curve slightly forward in the direction of motion and have an elevated \hi\ column density, indicating that these wings are forming Rayleigh-Taylor globules at their tips and that this gas is becoming entangled with unseen vortices in the surrounding coronal gas. These observations provide new insights on the survivability of low-metallicity gas streams that are accreting onto $L_\star$ galaxies.

The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km$^2$ large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down to about 10$^{17}$ eV. At the depth of 2.3 m the electromagnetic component of cosmic ray showers is almost entirely absorbed so that the buried scintillators provide an independent and direct measurement of the air showers muon content. This work describes the design and implementation of the AMIGA embedded system, which provides centralized control, data acquisition and environment monitoring to its detectors. The presented system was firstly tested in the engineering array phase ended in 2017, and lately selected as the final design to be installed in all new detectors of the production phase. The system was proven to be robust and reliable and has worked in a stable manner since its first deployment.

Accurate estimation of cratering asymmetry on the Moon is crucial for understanding Moon evolution history. Early studies of cratering asymmetry have omitted the contributions of high lunar obliquity and inclination. Here, we include lunar obliquity and inclination as new controlling variables to derive the cratering rate spatial variation as a function of longitude and latitude. With examining the influence of lunar obliquity and inclination on the asteroids population encountered by the Moon, we then have derived general formulas of the cratering rate spatial variation based on the crater scaling law. Our formulas with addition of lunar obliquity and inclination can reproduce the lunar cratering rate asymmetry at the current Earth-Moon distance and predict the apex/ant-apex ratio and the pole/equator ratio of this lunar cratering rate to be 1.36 and 0.87, respectively. The apex/ant-apex ratio is decreasing as the obliquity and inclination increasing. Combining with the evolution of lunar obliquity and inclination, our model shows that the apex/ant-apex ratio does not monotonically decrease with Earth-Moon distance and hence the influences of obliquity and inclination are not negligible on evolution of apex/ant-apex ratio. This model is generalizable to other planets and moons, especially for different spin-orbit resonances.

Is there oceanic superrotation on exoplanets? Atmospheric superrotation, characterized by west-to-east winds over the equator, is a common phenomenon in the atmospheres of Venus, Titan, Saturn, Jupiter, and tidally locked exoplanets. The stratospheric atmosphere of Earth is also superrotating during the westerly phase of the quasi-biennial oscillation (QBO). However, whether the same phenomenon can occur in ocean is poorly known. Through numerical simulations, here we show that oceanic superrotation does occur on tidally locked terrestrial planets around low-mass stars. Its formation (spun-up from rest) is associated with surface winds, the equatorward momentum convergence by Rossby waves, and the eastward propagation of Kelvin waves in the ocean. Its maintenance is driven by equatorward momentum transports of coupled Rossby-Kelvin waves in the ocean excited from the uneven stellar radiation distribution. The width of the superrotation is mainly constrained by the Rossby deformation radius in the ocean, while its strength is more complex. Many factors can influence the strength, including planetary rotation rate, stellar flux, greenhouse gas concentration, seawater salinity, bottom drag, and a scaling theory is lack. This work confirms that superrotation can occur on tidally locked terrestrial planets with seawater oceans and suggests that it may also occur on tidally locked hot planets with magma oceans that will possibly be observed in the near future.

The Dark Energy Spectroscopic Instrument (DESI) is a Stage IV ground-based dark energy experiment that will measure the expansion history of the Universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasars over 14000 square degrees will be measured during the life of the experiment. We describe the installation of the major elements of the instrument at the Mayall 4m telescope, completed in late 2019. The previous prime focus corrector, spider vanes, and upper rings were removed from the Mayall's Serrurier truss and replaced with the newly-constructed DESI ring, vanes, cage, hexapod, and optical corrector. The new corrector was optically aligned with the primary mirror using a laser tracker system. The DESI focal plane system was integrated to the corrector, with each of its ten 500-fiber-positioner petal segments installed using custom installation hardware and the laser tracker. Ten DESI spectrographs with 30 cryostats were installed in a newly assembled clean room in the Large Coude Room. The ten cables carrying 5000 optical fibers from the positioners in the focal plane were routed down the telescope through cable wraps at the declination and hour angle axes, and their integral slitheads were integrated with the ten spectrographs. The fiber view camera assembly was installed to the Mayall's primary mirror cell. Servers for the instrument control system replaced existing computer equipment. The fully integrated instrument has been commissioned and is ready to start its operations phase.

Carbon stars, enhanced in carbon and neutron-capture elements, provide wealth of information about the nucleosynthesis history of the Galaxy. In this work, we present the first ever detailed abundance analysis of carbon star LAMOSTJ091608.81+230734.6 and a detailed abundance analysis of neutron-capture elements for the object LAMOSTJ151003.74+305407.3. Updates on the abundances of elements C, O, Mg, Ca, Cr, Mn and Ni for LAMOSTJ151003.74+305407.3 are also presented. Our analysis is based on high resolution spectra obtained using Hanle Echelle Spectrograph (HESP) attached to the Himalayan Chandra Telescope (HCT), IAO, Hanle. The stellar atmospheric parameters (T$_{eff}$, logg, micro-turbulance ${\zeta}$, metallicity [Fe/H]) are found to be (4820, 1.43, 1.62, $-$0.89) and (4500, 1.55, 1.24, $-$1.57) for these two objects respectively. The abundance estimates of several elements, C, N, O, Na, $\alpha$-elements, Fe-peak elements and neutron-capture elements Rb, Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm and Eu are presented. Our analysis shows the star LAMOSTJ151003.74+305407.3 to be a CEMP-r/s star, and LAMOSTJ091608.81+230734.6 a CH giant. We have examined if the i-process model yields ([X/Fe]) of heavy elements could explain the observed abundances of the CEMP-r/s star based on a parametric model based analysis. The negative values obtained for the neutron density dependent [Rb/Zr] ratio confirm former low-mass AGB companions for both the stars. Kinematic analysis shows that LAMOSTJ151003.74+305407.3 belongs to the Galactic halo population and LAMOSTJ091608.81+230734.6 to the disc population.

All the elements heavier than Fe are produced either by slow (-s) or rapid (-r) neutron-capture process. Neutron density prevailing in the stellar sites is one of the major factors that determines the type of neutron-capture processes. We present the results based on the estimates of corrected value of absolute carbon abundance, [C/N] ratio, carbon isotopic ratio and [hs/ls] ratio obtained from the high resolution spectral analysis of six stars that include both CH stars and CEMP stars. All the stars show enhancement of neutron-capture elements. Location of these objects in the A(C) vs. [Fe/H] diagram shows that they are Group I objects, with external origin of carbon and neutron-capture elements. Low values of carbon isotopic ratios estimated for these objects may also be attributed to some external sources. As the carbon isotopic ratio is a good indicator of mixing, we have used the estimates of 12C/13C ratios to examine the occurance of mixing in the stars. While the object HD 30443 might have experienced an extra mixing process that usually occurs after red giant branch (RGB) bump for stars with log(L/L0) > 2.0, the remaining objects do not show any evidence of having undergone any such mixing process. The higher values of [C/N] ratios obtained for these objects also indicate that none of these objects have experienced any strong internal mixing processes. Based on the estimated abundances of carbon and the neutron-capture elements, and the abundance ratios, we have classified the objects into different groups. While the objects HE 0110-0406, HD 30443 and CD-38 2151 are found to be CEMP-s stars, HE 0308-1612 and HD 176021 show characteristic properties of CH stars with moderate enhancement of carbon. The object CD-28 1082 with enhancement of both r- and s-process elements is found to belong to the CEMP-r/s group.

Elliptical galaxies have dynamically hot ($\sigma_{\rm 1D}$ ~ 100--300 km s$^{-1}$) populations of stars, and presumably, smaller objects like comets. Because interstellar minor bodies are moving much faster, they hit planets harder and more often than in the local Galaxy. I estimate the rates for Chicxulub-scale impacts on an Earth-size planet in elliptical galaxies as a potential habitability constraint on intelligent life. Around most stars in a normal elliptical galaxy, these planets receive only ~ 0.01--0.1 Gyr$^{-1}$, although hazardous rates may be common in certain compact early-type galaxies. About ~10% of the stellar mass is in a region where the rate is >10 Gyr$^{-1}$, large enough to dominate the mass extinction rate. This suggests that elliptical galaxies have an exclusion zone several hundred parsecs in radius around their centers for the evolution of intelligent life.

The cohesion of particles has a significant effect on the properties of small bodies. In this study, we measured in open air, the cohesive forces of tens of micron-sized irregularly shaped meteorite, silica sand, glass powder, and spherical glass particles, using a centrifugal method. In addition, we estimated the amount of water vapor adsorbed on the particles under the measurement conditions. The measured cohesive forces of the meteorite particles are tens of times smaller than those of an ideally spherical silica particle and correspond to the submicron-scale effective (or equivalent) curvature radius of the particle surface. Moreover, based on the estimated amount of water vapor adsorbed on the particles, we expect the cohesive forces of the particles in airless bodies to be approximately 10 times larger than those measured in open air. Based on the measurement results, we estimate that the cohesive forces of the particles on asteroids are typically in the sub-micro-Newton range, and that the particles on fast-rotating asteroids are tens of microns in size.

We report on the detection of source noise in the time domain at 162MHz with the Murchison Widefield Array. During the observation the flux of our target source Virgo A (M87) contributes only $\sim$1\% to the total power detected by any single antenna, thus this source noise detection is made in an intermediate regime, where the source flux detected by the entire array is comparable with the noise from a single antenna. The magnitude of source noise detected is precisely in line with predictions. We consider the implications of source noise in this moderately strong regime on observations with current and future instruments.

The Andromeda (M31) galaxy displays several substructures in its inner halo whose origin as remnants of accreted satellites or perturbations of the pre-existing disc are encoded in the properties of their stellar populations (SPs), leaving traces on their deep [OIII] 5007 \AA planetary nebulae luminosity functions (PNLFs). By characterizing the morphology of the PNLFs, we constrain their origin. From our 54 sq. deg. deep narrow-band [OIII] survey of M31, we identify planetary nebulae (PNe) in the M31 disc and six major inner-halo substructures - the Giant Stream, North East Shelf, G1-Clump, Northern Clump, Western Shelf and Stream-D. We measure PNLF parameters from cumulative fits and statistically compare the PNLFs in each substructure and the disc. We link the PNLF parameters and those for the Large Magellanic Cloud to published metallicities and age measurements for their parent SPs. The absolute magnitudes of the PNLF bright cut-off ($M^{*}$) for these sub-populations span a significant magnitude range, despite having similar distance and line-of-sight extinction. $M^{*}$ for the Giant Stream, W-shelf and Stream-D PNLFs are fainter than those predicted by PN evolution models for the metallicity of the parent SPs. The faint-end slope of the PNLF increases linearly with decreasing fraction of stellar mass younger than 5 Gyr across the M31 regions and the LMC. From their PNLFs, the Giant Stream and NE-shelf are consistent with being stellar debris from an infalling satellite, while the G1 Clump appears to be linked with the pre-merger disc. The SPs of the substructures are consistent with those predicted by simulations of a single massive merger event that took place 2--3 Gyr ago in M31. Stream-D has an unrelated, distinct, origin. Furthermore, this study provides independent evidence that the faint-end of the PNLF is preferentially populated by PNe evolved from older stars.

We investigate the cosmological observational test of the extended quintessence model, i.e. a scalar-tensor gravity model with a scalar field potential serving as dark energy, by using the Planck 2018 cosmic microwave background (CMB) data, together with the baryon acoustic oscillations (BAO) and redshift-space distortion (RSD) data. As an example, we consider the model with a Brans-Dicke kinetic term and a quadratic scalar potential, which reduces to general relativity (GR) in the limit $\omega(\phi) \to \infty$, and the cosmological constant in the limit $B=C=0$. In such a model the scalar field typically rolls down the potential and oscillates around the minimum of $V (\phi)$. We find that the model parameter estimate for the CMB+BAO+RSD data set is given by $\lg \alpha = -3.6 _{-0.54}^{+0.66}~ (68\%)$, corresponding to $ 3.8 \times 10^5 < \omega_0 < 9.5 \times 10^7~(68\%)$, and $\lg C = 4.9 \pm 1.4~ (68\%) $. However, the GR $\Lambda$CDM model can fit the data almost as good as this extended quintessence model, and is favored by the Akaike information criterion (AIC). The variation of the gravitational constant since the epoch of Recombination is constrained to be $0.97 < G_{\rm rec}/G_0 < 1.03~ (1 \sigma)$. In light of recent report that the CMB data favors a closed universe, we consider the case with non-flat geometry in our fit, and find that the mean value of $\Omega_k$ shifts a little bit from $-0.049$ to $-0.036$, and the parameters in our model are not degenerate with $\Omega_k$.

Gamma ray bursts (GRBs) are astronomical phenomena detected at highest energies. The gamma ray photons carry energies on the order of mega-electronovolts and arrive to us from the point-like sources that are uniformly distributed on the sky. A typical burst has a form of a pulse that lasts for about a minute. As the Earth atmosphere is not transparent to the very high energy radiation, the bursts are detected by means of the telescopes onboard satellites that are placed on the orbit. The total energetics of GRB events, which is given by the integrated energy flux by the detector unit area, implies that we are witnessing very powerful explosions, where an enormously great power is released within a short time. There is only one way to obtain such huge energies in cosmos: the disruption of a star.

Following an eruptive accretion event in NGC6334I-MM1, flares in the various maser species, including water masers, were triggered. We report the observed relative proper motion of the highly variable water masers associated with the massive star-forming region, NGC6334I. High velocity H$_2$O maser proper motions were detected in 5 maser clusters, CM2-W2 (bow-shock structure), MM1-W1, MM1-W3, UCHII-W1 and UCHII-W3. The overall average of the derived relative proper motion is 85 km s$^{-1}$. This mean proper motion is in agreement with the previous results from VLA multi-epoch observations. Our position and velocity variance and co-variance matrix analyses of the maser proper motions show its major axis to have a position angle of $-$79.4$^\circ$, cutting through the dust cavity around MM1B and aligned in the northwest-southeast direction. We interpret this as the axis of the jet driving the CM2 shock and the maser motion. The complicated proper motions in MM1-W1 can be explained by the combined influence of the MM1 northeast-southwest bipolar outflow, CS(6-5) north-south collimated bipolar outflow, and the radio jet. The relative proper motions of the H$_2$O masers in UCHII-W1 are likely not driven by the jets of MM1B protostar but by MM3-UCHII. Overall, the post-accretion burst relative proper motions of the H$_2$O masers trace shocks of jet motion.

By fitting stellar populations to SDSS-IV MaNGA survey observations of ~7000 suitably-weighted individual galaxies, we reconstruct the star-formation history of the Universe, which we find to be in reasonable agreement with previous studies. Dividing the galaxies by their present-day stellar mass, we demonstrate the downsizing phenomenon, whereby the more massive galaxies hosted the most star-formation at earlier times. Further dividing the galaxy sample by colour and morphology, we find that a galaxy's present-day colour tells us more about its historical contribution to the cosmic star formation history than its current morphology. We show that downsizing effects are greatest among galaxies currently in the blue cloud, but that the level of downsizing in galaxies of different morphologies depends quite sensitively on the morphological classification used, due largely to the difficulty in classifying the smaller low-mass galaxies from their ground-based images. Nevertheless, we find agreement that among galaxies with stellar masses $M_{\star}>6\times10^{9}\,M_{\odot}$, downsizing is most significant in spirals. However, there are complicating factors. For example, for more massive galaxies, we find that colour and morphology are predictors of the past star formation over a longer timescale than in less massive systems. Presumably this effect is reflecting the longer period of evolution required to alter these larger galaxies' physical properties, but shows that conclusions based on any single property don't tell the full story.

The past year has seen numerous publications underlining the importance of a space mission to the ice giants in the upcoming decade. Proposed mission plans involve a Jupiter swing-by followed by a $\sim$10 year cruise time to the ice giants. The cruise time can be utilized to search for low-frequency gravitational waves by observing the Doppler shift caused by them in the Earth-spacecraft radio link. We calculate the sensitivity of prospective ice giant missions to gravitational waves. Then, adopting a steady-state black hole binary population, we derive a conservative estimate for the detection rate of extreme mass ratio inspirals (EMRIs), supermassive (SMBH) and stellar mass binary black hole (sBBH) mergers. We link the SMBHs population to the fraction of quasars $f_\rm{bin}$ resulting from galaxy mergers that pair SMBHs to a binary. For a total of ten 40-day observations during the cruise of a single spacecraft, $\mathcal{O}(f_\rm{bin})\sim0.5$ detections of SMBH mergers are likely, if Allan deviation of Cassini-era noise is improved by $\sim 10^2$. For EMRIs we estimate this number to lie between $\sim\mathcal{O}(0.1) - \mathcal{O}(100)$, and for sBBHs around $\sim \mathcal{O}(0.01)$.

We gave an extensive study for the quasi-periodic perturbations on the time profiles of the line of sight (LOS) magnetic field in 10x10 sub-areas in a solar plage region (corresponds to a facula on the photosphere). The perturbations are found to be associated with enhancement of He I 10830 A absorption in a moss region, which is connected to loops with million-degree plasma. FFT analysis to the perturbations gives a kind of spectrum similar to that of Doppler velocity: a number of discrete periods around 5 minutes. The amplitudes of the magnetic perturbations are found to be proportional to magnetic field strength over these sub-areas. In addition, magnetic perturbations lag behind a quarter of cycle in phase with respect to the p-mode Doppler velocity. We show that the relationships can be well explained with an MHD solution for the magneto-acoustic oscillations in high-\b{eta} plasma. Observational analysis also shows that, for the two regions with the stronger and weaker magnetic field, the perturbations are always anti-phased. All findings show that the magnetic perturbations are actually magneto-acoustic oscillations on the solar surface, the photosphere, powered by p-mode oscillations. The findings may provide a new diagnostic tool for exploring the relationship between magneto-acoustic oscillations and the heating of solar upper atmosphere, as well as their role in helioseismology.

MUSE-based emission-line maps of the spiral galaxy NGC 4030 reveal the existence of unresolved sources with forbidden line emission enhanced with respect to those seen in its own HII regions. This study reports our efforts to detect and isolate these objects and identify their nature. Candidates are first detected as unresolved sources on an image of the second principal component of the Hb, [OIII]5007, Ha, [NII]6584, [SII]6716, 6731 emission-line data cube, where they stand out clearly against both the dominant HII region population and the widespread diffuse emission. The intrinsic emission is then extracted accounting for the highly inhomogeneous emission-line "background" throughout the field of view. Collisional to recombination line ratios like [SII]/Ha, [NII]/Ha, and [OI]/Ha tend to increase when the background emission is corrected for. We find that many (but not all) sources detected with the principal component analysis have properties compatible with supernova remnants (SNRs). Applying a combined [SII]/Ha and [NII]/Ha classification criterion leads to a list of 59 sources with SNR-like emission lines. Many of them exhibit conspicuous spectral signatures of SNRs around 7300 Angs, and a stacking analysis shows that these features are also present, except weaker, in other cases. At nearly 30 Mpc, these are the most distant SNRs detected by optical means to date. We further report the serendipitous discovery of a luminous (M_V ~ -12.5), blue, and variable source, possibly associated with a supernova impostor.

Light scattering at near-infrared wavelengths has been used to study the optical properties of the interstellar dust grains, but these studies are limited by the assumptions on the strength of the radiation field. On the other hand, thermal dust emission can be used to constrain the properties of the radiation field, although this is hampered by uncertainty about the dust emissivity. We test if current dust models allow us to model a molecular cloud simultaneously in the near infrared (NIR) and far infrared (FIR) wavelengths and compare the results with observations. Our aim is to place constraints on the properties of the dust grains and the strength of the radiation field. We present computations of dust emission and scattered light of a quiescent molecular cloud LDN1512. We construct radiative transfer models for LDN1512 that include an anisotropic radiation field and a three-dimensional cloud model. We are able to reproduce the observed FIR observations, with a radiation field derived from the DIRBE observations, with all of the tested dust models. However, with the same density distribution and the assumed radiation field, the models fail to reproduce the observed NIR scattering in all cases except for models that take into account dust evolution via coagulation and mantle formation. We find that the column densities derived from our radiative transfer modelling can differ by a factor of up to two, compared to the column densities derived from the observations with modified blackbody fits. The discrepancy in the column densities is likely caused because of temperature difference between a modified blackbody fit and the real spectra. We show that the observed dust emission can be reproduced with several different assumptions about the properties of the dust grains. However, in order to reproduce the observed scattered surface brightness dust evolution must be taken into account.

Only 20\% of old field stars have detectable debris discs, leaving open the question of what disc, if any, is present around the remaining 80\%. Young moving groups allow to probe this population, since discs are expected to have been brighter early on. This paper considers the population of F~stars in the 23~Myr-old BPMG where we find that 9/12 targets possess discs. We also analyse archival ALMA data to derive radii for 4 of the discs, presenting the first image of the 63au radius disc of HD~164249. Comparing the BPMG results to disc samples from $\sim45$~Myr and $\sim150$~Myr-old moving~groups, and to discs found around field stars, we find the disc incidence rate in young moving~groups is comparable to that of the BPMG and significantly higher than that of field~stars. The BPMG discs tend to be smaller than those around field~stars. However, this difference is not statistically significant due to the small number of targets. Yet, by analysing the fractional luminosity vs disc radius parameter space we find that the fractional luminosities in the populations considered drop by two orders of magnitude within the first 100~Myr. This is much faster than expected by collisional evolution, implying a decay equivalent to $1/\text{age}^2$. We attribute this depletion to embedded planets which would be around 170~$M_\text{earth}$ to cause a depletion on the appropriate timescale. However, we cannot rule out that different birth environments of nearby young clusters result in brighter debris discs than the progenitors of field~stars which likely formed in a more dense environment.

We present the results of our multifrequency observations for 10 ultraluminous infrared galaxies (ULIRGs) made by the Karl G. Jansky Very Large Array at 1.4, 5.5, 9.0, and 14.0 GHz. Our sample was selected from ULIRGs whose active galactic nuclei (AGNs) are not found at optical (~70% of the entire ULIRGs), but whose presence is suggested at mid-infrared or submillimetre (>50% of the non-AGN ULIRGs at optical). The statistical properties of the targets are similar to those of the entire ULIRG sample, which implies that ULIRGs have common radiative processes regardless of the presence of optical AGNs, and thus AGNs might equally contribute to the radio emissions of every ULIRG. Although their spectra are mainly explained by starburst, some individual sources suggest contributions from AGNs. IRAS 00091-0738, IRAS 00188-0856, and IRAS 01298-0744 show spectral breaks at high frequencies, which can be explained by synchrotron aging of nonthermal plasma emitted from AGNs. In addition, we found 100-kpc scale extended emission associated with IRAS 01004-2237. The two-sided morphology and absence of extended X-ray emission suggest that this system is not induced by a merger in a cluster but originates from AGN activity.

We investigate a recent claim that observed galaxy clusters produce an order of magnitude more galaxy-galaxy strong lensing (GGSL) than simulated clusters in a LCDM cosmology. We take galaxy clusters from the C-EAGLE hydrodynamical simulations and calculate the expected amount of GGSL for sources placed behind the clusters at different redshifts. The probability of a source lensed by one of the most massive C-EAGLE clusters being multiply imaged by an individual cluster member is in good agreement with that inferred for observed clusters. We show that numerically converged results for the GGSL probability require higher resolution simulations than had been used previously. On top of this, different galaxy formation models predict cluster substructures with different central densities, such that the GGSL probabilities in LCDM cannot yet be robustly predicted. Overall, we find that galaxy-galaxy strong lensing within clusters is not currently in tension with the LCDM cosmological model.

The GRAVITY instrument on the ESO VLTI pioneers the field of high-precision near-infrared interferometry by providing astrometry at the $10 - 100\,\mu$as level. Measurements at such high precision crucially depend on the control of systematic effects. Here, we investigate how aberrations introduced by small optical imperfections along the path from the telescope to the detector affect the astrometry. We develop an analytical model that describes the impact of such aberrations on the measurement of complex visibilities. Our formalism accounts for pupil-plane and focal-plane aberrations, as well as for the interplay between static and turbulent aberrations, and successfully reproduces calibration measurements of a binary star. The Galactic Center observations with GRAVITY in 2017 and 2018, when both Sgr A* and the star S2 were targeted in a single fiber pointing, are affected by these aberrations at a level of less than 0.5 mas. Removal of these effects brings the measurement in harmony with the dual beam observations of 2019 and 2020, which are not affected by these aberrations. This also resolves the small systematic discrepancies between the derived distance $R_0$ to the Galactic Center reported previously.

Strong gravitational lensing observations can test structure formation models by constraining the masses and concentrations of subhaloes in massive galaxy clusters. Recent work has concluded that cluster subhaloes are more abundant and/or concentrated than predicted by $\Lambda$CDM simulations; this finding has been interpreted as a arising from unidentified issues with simulations or an incorrect understanding of the nature of dark matter. We test these hypotheses by comparing observed subhalo masses and maximum circular velocities $v_\mathrm{max}$ to predictions from the high resolution Hydrangea galaxy cluster simulation suite, which is based on the successful EAGLE galaxy formation model. The simulated subhalo mass distribution and mass-$v_\mathrm{max}$ relation agrees well with observations, due to the presence of baryons during tidal stripping. Similar agreement is found for the lower-resolution Illustris-TNG300 simulation. In combination, our results suggest that the abundance and concentration of cluster substructures are not in tension with $\Lambda$CDM, but may provide useful constraints for the refinement of baryon physics models in simulations.

We present the first spectroscopically resolved \ha\ emission map of the Large Magellanic Cloud's (LMC) galactic wind. By combining new Wisconsin H-alpha Mapper (WHAM) observations ($I_{\rm H\alpha}\gtrsim10~{\rm mR}$) with existing \hicm\ emission observations, we have (1) mapped the LMC's near-side galactic wind over a local standard of rest (LSR) velocity range of $+50\le\rm v_{LSR}\le+250~{\rm km}~{\rm s}^{-1}$, (2) determined its morphology and extent, and (3) estimated its mass, outflow rate, and mass-loading factor. We observe \ha\ emission from this wind to typically 1-degree off the LMC's \hi\ disk. Kinematically, we find that the diffuse gas in the warm-ionized phase of this wind persists at both low ($\lesssim100~{\rm km}~{\rm s}^{-1}$) and high ($\gtrsim100~{\rm km}~{\rm s}^{-1}$) velocities, relative to the LMC's \hi\ disk. Furthermore, we find that the high-velocity component spatially aligns with the most intense star-forming region, 30~Doradus. We, therefore, conclude that this high-velocity material traces an active outflow. We estimate the mass of the warm ($T_e\approx10^4~\rm K$) ionized phase of the near-side LMC outflow to be $\log{\left(M_{\rm ionized}/M_\odot\right)=7.51\pm0.15}$ for the combined low and high velocity components. Assuming an ionization fraction of 75\% and that the wind is symmetrical about the LMC disk, we estimate that its total (neutral and ionized) mass is $\log{\left(M_{\rm total}/M_\odot\right)=7.93}$, its mass-flow rate is $\dot{M}_{\rm outflow}\approx1.43~M_\odot~\rm yr^{-1}$, and its mass-loading factor is $\eta\approx4.54$. Our average mass-loading factor results are roughly a factor of 2.5 larger than previous \ha\ imaging and UV~absorption line studies, suggesting that those studies are missing nearly half the gas in the outflows.

As exoplanetary science matures into its third decade, we are increasingly offered the possibility of pre existing, archival observations for newly detected candidates. This is particularly poignant for the TESS mission, whose survey spans bright, nearby dwarf stars in both hemispheres, which are precisely the types of sources targeted by previous radial velocity (RV) surveys. On this basis, we investigated whether any of the TESS Objects of Interest (TOIs) coincided with such observations, from which we find 18 single planet candidate systems. Of these, one exhibits an RV signature that has the correct period and phase matching the transiting planetary candidate with a false alarm probability of less than 1 percent. After further checks, we exploit this fact to validate HD 183579b (TOI-1055b). This planet is less than 4 Earth Radii and has better than 33 percent planetary mass measurements, thus advancing the TESS primary objective of finding 50 such worlds. We find that this planet is amongst the most accessible small transiting planets for atmospheric characterization. Our work highlights that the efforts to confirm and even precisely measure the masses of new transiting planet candidates need not always depend on acquiring new observations - that in some instances these tasks can be completed with existing data.

Aims. We aim to determine whether dissociative excitation of cometary neutrals by electron impact is the major source of far-ultraviolet (FUV) emissions at comet 67P/Churyumov-Gerasimenko in the southern hemisphere at large heliocentric distances, both during quiet conditions and impacts of corotating interaction regions observed in the summer of 2016. Methods. We combined multiple datasets from the Rosetta mission through a multi-instrument analysis to complete the first forward modelling of FUV emissions in the southern hemisphere of comet 67P and compared modelled brightnesses to observations with the Alice FUV imaging spectrograph. We modelled the brightness of OI1356, OI1304, Lyman-$\beta$, CI1657, and CII1335 emissions, which are associated with the dissociation products of the four major neutral species in the coma: CO$_2$, H$_2$O, CO, and O$_2$. The suprathermal electron population was probed by RPC/IES and the neutral column density was constrained by several instruments: ROSINA, MIRO and VIRTIS. Results. The modelled and observed brightnesses of the FUV emission lines agree closely when viewing nadir and dissociative excitation by electron impact is shown to be the dominant source of emissions away from perihelion. The CII1335 emissions are shown to be consistent with the volume mixing ratio of CO derived from ROSINA. When viewing the limb during the impacts of corotating interaction regions, the model reproduces brightnesses of OI1356 and CI1657 well, but resonance scattering in the extended coma may contribute significantly to the observed Lyman-$\beta$ and OI1304 emissions. The correlation between variations in the suprathermal electron flux and the observed FUV line brightnesses when viewing the comet's limb suggests electrons are accelerated on large scales and that they originate in the solar wind. This means that the FUV emissions are auroral in nature.

ACO2163 is one of the hottest (mean $kT=12-15.5$ keV) and extremely X-ray overluminous merging galaxy clusters which is located at $z=0.203$. The cluster hosts one of the largest giant radio halos which are observed in most of the merging clusters, and a candidate radio relic. Recently, three merger shock fronts were detected in this cluster, explaining its extreme temperature and complex structure. Furthermore, previous XMM-Newton and Chandra observations hinted at the presence of a shock front that is associated with the gas `bullet' crossing the main cluster in the west-ward direction, and which heated the intra-cluster medium, leading to adiabatic compression of the gas behind the 'bullet'. The goal of this paper is to report on the detection of this shock front as revealed by the temperature discontinuity in the X-ray XMM-Newton image, and the edge in the Very Large Array (VLA) radio image. We also report on the detection of a relic source in the north-eastern region of the radio halo in the KAT-7 data, confirming the presence of an extended relic in this cluster. The brightness edge in the X-rays corresponds to a shock front with a Mach number $M= 2.2\pm0.3$, at a distance of 0.2 Mpc from the cluster centre. An estimate from the luminosity jump gives $M=1.9\pm0.4$. We consider a simple explanation for the electrons at the shock front, and for the observed discrepancy between the average spectral index of the radio halo emission and that predicted by the $M=2.2$ shock which precedes the 'bullet'.

The recently commissioned Dark Energy Spectroscopic Instrument (DESI) will measure the expansion history of the Universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasars over 14000 sqdeg will be measured during the life of the experiment. A new prime focus corrector for the KPNO Mayall telescope delivers light to 5000 fiber optic positioners. The fibers in turn feed ten broad-band spectrographs. A consortium of Aix-Marseille University (AMU) and CNRS laboratories (LAM, OHP and CPPM) together with LPNHE (CNRS, IN2P3, Sorbonne Universit\'e and Universit\'e de Paris) and the WINLIGHT Systems company based in Pertuis (France), were in charge of integrating and validating the performance requirements of the ten full spectrographs, equipped with their cryostats, shutters and other mechanisms. We present a summary of our activity which allowed an efficient validation of the systems in a short-time schedule. We detail the main results. We emphasize the benefits of our approach and also its limitations.

We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advected to Eulerian space using the displacement field of numerical simulations. The cosmology-dependence is captured thanks to a cosmology-rescaling algorithm. We show that our emulator is capable of describing the power spectrum of galaxy formation simulations for a sample mimicking that of an typical Emission-Line survey at $z \sim 1$ with an accuracy of $1-2\%$ up to nonlinear scales $k \sim 0.7 h \mathrm{Mpc}^{-1}$.

In recent years, weak lensing of the cosmic microwave background (CMB) has emerged as a powerful tool to probe fundamental physics, such as neutrino masses, primordial non-Gaussianity, dark energy, and modified gravity. The prime target of CMB lensing surveys is the lensing potential, which is reconstructed from observed CMB temperature $T$ and polarization $E$ and $B$ fields. Until very recently, this reconstruction has been performed with quadratic estimators (QEs), which, although known to be suboptimal for high-sensitivity experiments, are numerically efficient, and useful to make forecasts and cross-check the results of more sophisticated likelihood-based methods. It is expected that ongoing and near-future CMB experiments such as AdvACT, SPT-3G and the Simons Observatory (SO), will also rely on QEs. Here, we review different QEs, and clarify their differences. In particular, we show that the Hu-Okamoto (HO02) estimator is not the absolute optimal lensing estimator that can be constructed out of quadratic combinations of $T, E$ and $B$ fields. Instead, we derive the global-minimum-variance (GMV) lensing quadratic estimator. Although this estimator can be found elsewhere in the literature, it was erroneously described as equivalent to the HO02 estimator, and has never been used in real data analyses. Here, we show explicitly that the HO02 estimator is suboptimal to the GMV estimator, with a reconstruction noise larger by up to $\sim 9\%$ for a SO-like experiment. We further show that the QE used in the Planck, and recent SPT lensing analysis is suboptimal to both the HO02 and GMV estimator, and would have a reconstruction noise up to $\sim 11\%$ larger than that of the GMV estimator for a SO-like experiment. In addition to clarifying differences between different QEs, this work should thus provide motivation to implement the GMV estimator in future lensing analyses relying on QEs.

Context. Young and embedded stellar populations are important probes of the star formation process. Paradoxically, we have a better census of nearby embedded young populations than the slightly more evolved optically visible young populations. The high accuracy measurements and all-sky coverage of Gaia data are about to change this situation. Aims. This work aims to construct the most complete sample to date of YSOs in the $\rho$ Oph region. Methods. We compile a catalog of 761 Ophiuchus YSOs from the literature and crossmatch it with the Gaia EDR3, Gaia-ESO and APOGEE-2 surveys. We apply a multivariate classification algorithm to this catalog to identify new, co-moving population candidates. Results. We find 173 new YSO candidates in the Gaia EDR3 catalog belonging to the $\rho$ Oph region. The new sources appear to be mainly Class III M-stars and substellar objects and are less extincted than the known members. We find 19 previously unknown sources with disks. The analysis of the proper motion distribution of the entire sample reveals a well-defined bimodality, implying two distinct populations sharing a similar 3D volume. The first population comprises young stars' clusters around the $\rho$ Ophiuchi star and the main Ophiuchus clouds (L1688, L1689, L1709). In contrast, the second population is older ($\sim$ 10 Myr), dispersed, has a distinct proper motion, and is likely part of the Upper-Sco group. The two populations are moving away from each other at about 3.8 km/s, and will no longer overlap in about 4 Myr. Finally, we flag 47 sources in the literature as impostors, which are sources that exhibit large deviations from the average distance and proper motion properties of the $\rho$ Oph population. Our results show the importance of accurate 3D space and motion information for improved stellar population analysis. (Abridged)

Stellar tidal streams are sensitive tracers of the properties of the gravitational potential in which they orbit and detailed observations of their density structure can be used to place stringent constraints on fluctuations in the potential caused by, e.g., the expected populations of dark matter subhalos in the standard cold dark matter paradigm (CDM). Simulations of the evolution of stellar streams in live $N$-body halos without low-mass dark-matter subhalos, however, indicate that streams exhibit significant perturbations on small scales even in the absence of substructure. Here we demonstrate, using high-resolution $N$-body simulations combined with sophisticated semi-analytic and simple analytic models, that the mass resolutions of $10^4$--$10^5\,\rm{M}_{\odot}$ commonly used to perform such simulations cause spurious stream density variations with a similar magnitude on large scales as those expected from a CDM-like subhalo population and an order of magnitude larger on small, yet observable, scales. We estimate that mass resolutions of $\approx100\,\rm{M}_{\odot}$ ($\approx1\,\rm{M}_{\odot}$) are necessary for spurious, numerical density variations to be well below the CDM subhalo expectation on large (small) scales. That streams are sensitive to a simulation's particle mass down to such small masses indicates that streams are sensitive to dark matter clustering down to these low masses if a significant fraction of the dark matter is clustered or concentrated in this way, for example, in MACHO models with masses of $10$--$100\,\rm{M}_{\odot}$.

If the Peccei-Quinn symmetry associated to an axion has ever been restored after inflation, axion strings inevitably produce a contribution to the stochastic gravitational wave background. Combining effective field theory analysis with numerical simulations, we show that the resulting gravitational wave spectrum has logarithmic deviations from a scale invariant form with an amplitude that is significantly enhanced at low frequencies. As a result, a single ultralight axion-like particle with a decay constant larger than $10^{14}~{\rm GeV}$ and any mass between $10^{-18}~{\rm eV}$ and $10^{-28}~{\rm eV}$ leads to an observable gravitational wave spectrum and is compatible with constraints on the post-inflationary scenario from dark matter overproduction, isocurvature and dark radiation. Since the spectrum extends over a wide range of frequencies, the resulting signal could be detected by multiple experiments. We describe straightforward ways in which the Peccei-Quinn symmetry can be restored after inflation for such decay constants. We also comment on the recent possible NANOgrav signal in light of our results.

The scale of small-field inflation cannot be constrained via primordial gravitational waves through measurement of tensor-to-scalar ratio $r$. In this study, I show that if cosmic strings are produced after symmetry breaking at the end of hilltop supernatural inflation, this small-field inflation model can be tested through the production of gravitational waves from cosmic strings. Future experiments of gravitational wave detectors will determine or further constrain the parameter space in the model.

We demonstrate that global observations of high-energy cosmic rays contribute to understanding unique characteristics of a large-scale magnetic flux rope causing a magnetic storm in August 2018. Following a weak interplanetary shock on 25 August 2018, a magnetic flux rope caused an unexpectedly large geomagnetic storm. It is likely that this event became geoeffective because the flux rope was accompanied by a corotating interaction region and compressed by high-speed solar wind following the flux rope. In fact, a Forbush decrease was observed in cosmic-ray data inside the flux rope as expected, and a significant cosmic-ray density increase exceeding the unmodulated level before the shock was also observed near the trailing edge of the flux rope. The cosmic-ray density increase can be interpreted in terms of the adiabatic heating of cosmic rays near the trailing edge of the flux rope, as the corotating interaction region prevents free expansion of the flux rope and results in the compression near the trailing edge. A northeast-directed spatial gradient in the cosmic-ray density was also derived during the cosmic-ray density increase, suggesting that the center of the heating near the trailing edge is located northeast of Earth. This is one of the best examples demonstrating that the observation of high-energy cosmic rays provides us with information that can only be derived from the cosmic ray measurements to observationally constrain the three-dimensional macroscopic picture of the interaction between coronal mass ejections and the ambient solar wind, which is essential for prediction of large magnetic storms.

We study some aspects of dynamical compactification scenario where stabilisation of extra dimensions occurs due to presence the Gauss-Bonnet term and non-zero spatial curvature. In the framework of the model under consideration there exists two-stages scenario of evolution of a Universe: on the first stage, the space evolves from a totally anisotropic state to the state with 3-dimensional (corresponding to our "real", world) expanding and $D$-dimensional contracting isotropic subspaces; on the second stage, constant curvature of extra dimensions begins to play role and provide compactification of extra dimensions. It is already known that such a scenario is realizable when constant curvature of extra dimensions is negative. Here we show that a range of coupling constants for which exponential solutions with 3-dimensional expanding and $D$-dimensional contracting isotropic subspaces are stable is located in a zone where compactification solutions with positively curved extra space are unstable, so that two-stage scenario analogous to the one described above is \emph{not} realizable. Also we study "nearly-Friedmann", regime for the case of arbitrary constant curvature of extra dimensions and describe new parametrization of the general solution for the model under consideration which provide elegant way of describing areas of existence over parameters space.

We report results of a search for an isotropic gravitational-wave background (GWB) using data from Advanced LIGO's and Advanced Virgo's third observing run (O3) combined with upper limits from the earlier O1 and O2 runs. Unlike in previous observing runs in the advanced detector era, we include Virgo in the search for the GWB. The results are consistent with uncorrelated noise, and therefore we place upper limits on the strength of the GWB. We find that the dimensionless energy density $\Omega_{\rm GW}\leq 5.8\times 10^{-9}$ at the 95% credible level for a flat (frequency-independent) GWB, using a prior which is uniform in the log of the strength of the GWB, with 99% of the sensitivity coming from the band 20-76.6 Hz; $\leq 3.4 \times 10^{-9}$ at 25 Hz for a power-law GWB with a spectral index of 2/3 (consistent with expectations for compact binary coalescences), in the band 20-90.6 Hz; and $\leq 3.9 \times 10^{-10}$ at 25 Hz for a spectral index of 3, in the band 20-291.6 Hz. These upper limits improve over our previous results by a factor of 6.0 for a flat GWB. We also search for a GWB arising from scalar and vector modes, which are predicted by alternative theories of gravity; we place upper limits on the strength of GWBs with these polarizations. We demonstrate that there is no evidence of correlated noise of magnetic origin by performing a Bayesian analysis that allows for the presence of both a GWB and an effective magnetic background arising from geophysical Schumann resonances. We compare our upper limits to a fiducial model for the GWB from the merger of compact binaries. Finally, we combine our results with observations of individual mergers andshow that, at design sensitivity, this joint approach may yield stronger constraints on the merger rate of binary black holes at $z \lesssim 2$ than can be achieved with individually resolved mergers alone. [abridged]

This letter presents a combined measurement of the energy spectra of atmospheric $\nu_e$ and $\nu_\mu$ in the energy range between $\sim$100 GeV and $\sim$50 TeV with the ANTARES neutrino telescope. The analysis uses 3012 days of detector livetime in the period 2007--2017, and selects 1016 neutrinos interacting in (or close to) the instrumented volume of the detector, yielding shower-like events (mainly from $\nu_e+\overline \nu_e$ charged current plus all neutrino neutral current interactions) and starting track events (mainly from $\nu_\mu + \overline \nu_\mu$ charged current interactions). The contamination by atmospheric muons in the final sample is suppressed at the level of a few per mill by different steps in the selection analysis, including a Boosted Decision Tree classifier. The distribution of reconstructed events is unfolded in terms of electron and muon neutrino fluxes. The derived energy spectra are compared with previous measurements that, above 100 GeV, are limited to experiments in polar ice and, for $\nu_\mu$, to Super-Kamiokande.