Papers for Monday, Jun 28 2021

Jovian planet formation has been shown to be strongly correlated with host star metallicity, which is thought to be a proxy for disk solids. Observationally, previous works have indicated that jovian planets preferentially form around stars with solar and super solar metallicities. Given these findings, it is challenging to form planets within metal-poor environments, particularly for hot Jupiters that are thought to form via metallicity-dependent core accretion. Although previous studies have conducted planet searches for hot Jupiters around metal-poor stars, they have been limited due to small sample sizes, which are a result of a lack of high-quality data making hot Jupiter occurrence within the metal-poor regime difficult to constrain until now. We use a large sample of halo stars observed by TESS to constrain the upper limit of hot Jupiter occurrence within the metal-poor regime (-2.0 $\leq$ [Fe/H] $\leq$ -0.6). Placing the most stringent upper limit on hot Jupiter occurrence, we find the mean 1-$\sigma$ upper limit to be 0.18 $\%$ for radii 0.8 -2 R$_{\rm{Jupiter}}$ and periods $0.5- 10$ days. This result is consistent with previous predictions indicating that there exists a certain metallicity below which no planets can form.

For the first time, we present the simultaneous detection and characterization of three distinct phases at $>10^5$ K in $z=0$ absorption, using deep $\it{Chandra}$ observations toward Mrk 421. The extraordinarily high signal-to-noise ratio ($\geqslant60$) of the spectra has allowed us to detect a $\it{hot}$ phase of the Milky Way circumgalactic medium (CGM) at 3.2$^{+1.5}_{-0.5}\times$ 10$^7$ K, coexisting with a $\textit{warm-hot}$ phase at 1.5$\pm$0.1$\times$10$^6$ K and a $\textit{warm}$ phase at 3.0$\pm$0.4$\times$10$^5$ K. The $\textit{warm-hot}$ phase is at the virial temperature of the Galaxy, and the $\textit{warm}$ phase may have cooled from the $\textit{warm-hot}$ phase, but the super-virial $\textit{hot}$ phase remains a mystery. We find that [C/O] in the $\textit{warm}$ and $\textit{warm-hot}$ phases, [Mg/O] in the $\textit{warm-hot}$ phase and [Ne/O] in the $\textit{hot}$ phase are super-solar, and the $\textit{hot}$ and the $\textit{warm-hot}$ phases are $\alpha-$enhanced. Non-thermal line broadening is evident in the $\textit{warm-hot}$ and the $\textit{hot}$ phases and it dominates the total line broadening. Our results indicate that the $>10^5$ K CGM is a complex ecosystem. It provides insights on the thermal and chemical history of the Milky Way CGM, and theories of galaxy evolution.

Scalar Field Dark Matter (SFDM) comprised of ultralight ($\gtrsim 10^{-22}$ eV) bosons is an alternative to standard, collisionless Cold Dark Matter (CDM) that is CDM-like on large scales but inhibits small-scale structure formation. As a Bose-Einstein condensate, its free-field ("fuzzy") limit (FDM) suppresses structure below the de Broglie wavelength, $\lambda_\text{deB}$, creating virialized haloes with central cores of radius $\sim\lambda_\text{deB}$, surrounded by CDM-like envelopes, and a halo mass function (HMF) with a sharp cut-off on small scales. With a strong enough repulsive self-interaction (SI), structure is inhibited, instead, below the Thomas-Fermi (TF) radius, $R_\text{TF}$ (the size of an SI-pressure-supported ($n=1$)-polytrope), when $R_\text{TF} > \lambda_\text{deB}$. Previously, we developed tools to describe SFDM dynamics on scales above $\lambda_\text{deB}$ and showed that SFDM-TF haloes formed by Jeans-unstable collapse from non-cosmological initial conditions have $R_\text{TF}$-sized cores, surrounded by CDM-like envelopes. Revisiting SFDM-TF in the cosmological context, we simulate halo formation by cosmological infall and collapse, and derive its transfer function from linear perturbation theory to produce cosmological initial conditions and predict statistical measures of structure formation, such as the HMF. Since FDM and SFDM-TF transfer functions both have small-scale cut-offs, we can align them to let observational constraints on FDM proxy for SFDM-TF, finding FDM with particle masses $1 \lesssim m/(10^{-22} \text{ eV}/c^2) \lesssim 30$ corresponds to SFDM-TF with $10 \gtrsim R_\text{TF}/(1 \text{ pc}) \gtrsim 1$, favoring sub-galactic (sub-kpc) core-size. The SFDM-TF HMF cuts off gradually, however, leaving more small-mass haloes: its Jeans mass shrinks so fast, scales filtered early can still recover and grow!

Molecular cooling is essential for studying the formation of sub-structure of dissipative dark-matter halos that may host compact objects such as black holes. Here, we analyze the reaction rates relevant for the formation, dissociation, and transition of hydrogenic molecules while allowing for different values of the physical parameters: the coupling constant, the proton mass, and the electron mass. For all cases, we re-scale the reaction rates for the standard molecular hydrogen, so our results are valid as long as the dark matter is weakly coupled and one of the fermions is much heavier than the other. These results will allow a robust numerical treatment of cosmic structure, in particular for mini-halos for which molecular cooling is important, in a dissipative dark matter scenario.

Star-forming and starburst galaxies, which are well-known cosmic-rays reservoirs, are expected to emit gamma-rays and neutrinos predominantly via hadronic collisions. In this Letter, we analyze the 10-year Fermi-LAT spectral energy distributions of 13 nearby galaxies by means of a physical model which accounts for high-energy proton transport in starburst nuclei and includes the contribution of primary and secondary electrons. In particular, we test the hypothesis that the observed gamma-ray fluxes are mostly due to star-forming activity, in agreement with the available star formation rates coming from IR and UV observations. Through this observation-based approach, we determine the most-likely neutrino counterpart from star-forming and starburst galaxies and quantitatively assess the ability of current and upcoming neutrino telescopes to detect them as point-like sources. Remarkably, we find that the cores of the Small Magellanic Cloud and the Circinus galaxy are potentially observable by KM3NeT/ARCA with 6 years of observation. Moreover, most of the nearby galaxies are likely to be just a factor of a few below the KM3NeT and IceCube-Gen2 point-like sensitivities. After investigating the prospects for detection of gamma-rays above TeV energies from these sources, we conclude that the joint observations of high-energy neutrinos and gamma-rays with upcoming telescopes will be an objective test for our emission model and may provide compelling evidence of star-forming activity as a tracer of neutrino production.

We use three campaigns of K2 observations to complete the census of rotation in low-mass members of the benchmark, $\approx$670-Myr-old open cluster Praesepe. We measure new rotation periods (\prot) for 220 $\lesssim$1.3~\Msun\ Praesepe members and recover periods for $97\%$ (793/812) of the stars with a \prot\ in the literature. Of the 19 stars for which we do not recover a \prot, 17 were not observed by K2. As K2's three Praesepe campaigns took place over the course of three years, we test the stability of our measured \prot\ for stars observed in more than one campaign. We measure \prot\ consistent to within $10\%$ for $>95\%$ of the 331 likely single stars with $\geq$2 high-quality observations; the median difference in \prot\ is $0.3\%$, with a standard deviation of $2\%$. Nearly all of the exceptions are stars with discrepant \prot\ measurements in Campaign 18, K2's last, which was significantly shorter than the earlier two ($\approx$50~d rather than $\approx$75~d). This suggests that, despite the evident morphological evolution we observe in the light curves of $38\%$ of the stars, \prot\ measurements for low-mass stars in Praesepe are stable on timescales of several years. A \prot\ can therefore be taken to be representative even if measured only once.

By measuring the elemental abundances of a star, we can gain insight into the composition of its initial gas cloud -- the formation site of the star and its planets. Planet formation requires metals, the availability of which is determined by the elemental abundance. In the case where metals are extremely deficient, planet formation can be stifled. To investigate such a scenario requires a large sample of metal-poor stars and a search for planets therein. This paper focuses on the selection and validation of a halo star sample. We select ~17,000 metal-poor halo stars based on their Galactic kinematics, and confirm their low metallicities ([Fe/H] < -0.5), using spectroscopy from the literature. Furthermore, we perform high-resolution spectroscopic observations using LBT/PEPSI and conduct detailed metallicity ([Fe/H]) analyses on a sample of 13 previously known halo stars that also have hot kinematics. We can use the halo star sample presented here to measure the frequency of planets and to test planet formation in extremely metal-poor environments. The result of the planet search and its implications will be presented and discussed in a companion paper by Boley et al.

Understanding the interaction of massive black hole binaries with their gaseous environment is crucial since at sub-parsec scales the binary is too wide for gravitational wave emission to take over and to drive the two black holes to merge. We here investigate the interaction between a massive black hole binary and a self-gravitating circumbinary disc using 3D smoothed particle hydrodynamics simulations. We find that, when the disc self-gravity regulates the angular momentum transport, the binary semi-major axis decreases regardless the choice of disc masses and temperatures, within the range we explored. In particular, we find that the disc initial temperature (hence the disc aspect ratio) has little effect on the evolution of the binary since discs with the same mass self-regulate towards the same temperature. Initially warmer discs cause the binary to shrink on a slightly shorter timescale until the disc has reached the self-regulated equilibrium temperature. More massive discs drive the binary semi-major axis to decrease at a faster pace compared to less massive discs and result in faster binary eccentricity growth even after the initial-condition-dependent transient evolution. Finally we investigate the effect that the initial cavity size has on the binary-disc interaction and we find that, in the self-gravitating regime, an initially smaller cavity leads to a much faster binary shrinking, as expected. Our results are especially important for very massive black hole binaries such as those in the PTA band, for which gas self gravity cannot be neglected.

We review the formation and evolution of fossil groups and clusters from both the theoretical and the observational points of view. In the optical band, these systems are dominated by the light of the central galaxy. They were interpreted as old systems that had enough time to merge all the M* galaxies within the central one. During the last two decades many observational studies were performed to prove the old and relaxed state of fossil systems. The majority of these studies, that spans a wide range of topics including halos global scaling relations, dynamical substructures, stellar populations, and galaxy luminosity functions, seem to challenge this scenario. The general picture that can be obtained by reviewing all the observational works is that the fossil state could be transitional. Indeed, the formation of the large magnitude gap observed in fossil systems could be related to internal processes rather than an old formation.

Galaxy groups and poor clusters are more common than rich clusters, and host the largest fraction of matter content in the Universe. Hence, their studies are key to understand the gravitational and thermal evolution of the bulk of the cosmic matter. Moreover, because of their shallower gravitational potential, galaxy groups are systems where non-gravitational processes (e.g., cooling, AGN feedback, star formation) are expected to have a higher impact on the distribution of baryons, and on the general physical properties, than in more massive objects, inducing systematic departures from the expected scaling relations. Despite their paramount importance from the astrophysical and cosmological point of view, the challenges in their detection have limited the studies of galaxy groups. Upcoming large surveys will change this picture, reassigning to galaxy groups their central role in studying the structure formation and evolution in the Universe, and in measuring the cosmic baryonic content. Here, we review the recent literature on various scaling relations between X-ray and optical properties of these systems, focusing on the observational measurements, and the progress in our understanding of the deviations from the self-similar expectations on groups' scales. We discuss some of the sources of these deviations, and how feedback from supernovae and/or AGNs impacts the general properties and the reconstructed scaling laws. Finally, we discuss future prospects in the study of galaxy groups.

Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories. Widely used cosmological simulations with ~(100 Mpc)^3 volumes contain statistical samples of groups that provide important tests of galaxy evolution influenced by environmental processes. Larger volumes capable of reproducing LSS while following the redistribution of baryons by cooling and feedback are essential tools necessary to constrain cosmological parameters. Higher resolution simulations can currently model satellite interactions, the processing of cool (T~10^4 K) multi-phase gas, and non-thermal physics including turbulence, magnetic fields, and cosmic ray transport. We review simulation results regarding the gas and stellar contents of groups, cooling flows and the relation to the central galaxy, the formation and processing of multi-phase gas, satellite interactions with the intragroup medium, and the impact of groups for cosmological parameter estimation. Cosmological simulations provide evolutionarily consistent predictions of these observationally difficult-to-define objects, and have untapped potential to accurately model their gaseous, stellar, and dark matter distributions.

Galaxy groups host the majority of matter and more than half of all the galaxies in the Universe. Their hot ($10^7$ K), X-ray emitting intra-group medium (IGrM) reveals emission lines typical of many elements synthesized by stars and supernovae. Because their gravitational potentials are shallower than those of rich galaxy clusters, groups are ideal targets for studying, through X-ray observations, feedback effects, which leave important marks on their gas and metal contents. Here, we review the history and present status of the chemical abundances in the IGrM probed by X-ray spectroscopy. We discuss the limitations of our current knowledge, in particular due to uncertainties in the modeling of the Fe-L shell by plasma codes, and coverage of the volume beyond the central region. We further summarize the constraints on the abundance pattern at the group mass scale and the insight it provides to the history of chemical enrichment. Parallel to the observational efforts, we review the progress made by both cosmological hydrodynamical simulations and controlled high-resolution 3D simulations to reproduce the radial distribution of metals in the IGrM, the dependence on system mass from group to cluster scales, and the role of AGN and SN feedback in producing the observed phenomenology. Finally, we highlight future prospects in this field, where progress will be driven both by a much richer sample of X-ray emitting groups identified with eROSITA, and by a revolution in the study of X-ray spectra expected from micro-calorimeters onboard XRISM and ATHENA.

The co-evolution between supermassive black holes and their environment is most directly traced by the hot atmospheres of dark matter halos. Cooling of the hot atmosphere supplies the central regions with fresh gas, igniting active galactic nuclei (AGN) with long duty cycles. Outflows from the central engine tightly couple with the surrounding gaseous medium and provide the dominant heating source preventing runaway cooling by carving cavities and driving shocks across the medium. The AGN feedback loop is a key feature of all modern galaxy evolution models. Here we review our knowledge of the AGN feedback process in the specific context of galaxy groups. Galaxy groups are uniquely suited to constrain the mechanisms governing the cooling-heating balance. Unlike in more massive halos, the energy supplied by the central AGN to the hot intragroup medium can exceed the gravitational binding energy of halo gas particles. We report on the state-of-the-art in observations of the feedback phenomenon and in theoretical models of the heating-cooling balance in galaxy groups. We also describe how our knowledge of the AGN feedback process impacts on galaxy evolution models and on large-scale baryon distributions. Finally, we discuss how new instrumentation will answer key open questions on the topic.

The growing trove of precision astrometric observations from the Gaia space telescope and other surveys is revealing the structure and dynamics of the Milky Way in ever more exquisite detail. We summarize the current status of our understanding of the structure and the characteristics of the Milky Way, and we review the emerging picture: the Milky Way is evolving through interactions with the massive satellite galaxies that stud its volume, with evidence pointing to a cataclysmic past. It is also woven with stellar streams, and observations of streams, satellites, and field stars offer new constraints on its dark matter, both on its spatial distribution and its fundamental nature. The recent years have brought much focus to the study of dwarf galaxies found within our Galaxy's halo and their internal matter distributions. In this review, we focus on the predictions of the cold dark matter paradigm at small mass scales through precision astrometric measurements, and we summarize the modern consensus on the extent to which small-scale probes are consistent with this paradigm. We note the discovery prospects of these studies, and also how they intertwine with probes of the dynamics and evolution of the Milky Way in various and distinct ways.

Early dark energy (EDE) offers a particularly interesting theoretical approach to the Hubble tension, albeit one that introduces its own set of challenges, including a new `why then' problem related to the EDE injection time at matter-radiation equality, and a mild worsening of the large-scale structure (LSS) tension. Both these challenges center on the properties of dark matter, which becomes the dominant component of the Universe at EDE injection and is also responsible for seeding LSS. Motivated by this, we explore the potential of couplings between EDE and dark matter to address these challenges, focusing on a mechanism similar to chameleon dark energy theories, deeming this chameleon early dark energy (CEDE). We study the cosmological implications of such theories by fitting to the CMB, BAO, supernovae and the local value of $H_0$. We find that the Hubble tension is resolved by CEDE with $H_0 = 71.19(71.85)\pm 0.99$ km/s/Mpc. Further, the model provides an excellent fit to all the data, with no change to the CMB $\chi^2$ relative to a $\Lambda$CDM fit to just the CMB, BAO and SNe (i.e. excluding the $H_0$ tension for $\Lambda$CDM). We find a mild preference $(\sim 2\sigma)$ for the chameleon coupling constant $\beta >0$.

Binary systems of a hot subdwarf B (sdB) star + a white dwarf (WD) with orbital periods less than 2-3 hours can come into contact due to gravitational waves and transfer mass from the sdB star to the WD before the sdB star ceases nuclear burning and contracts to become a WD. Motivated by the growing class of observed systems in this category, we study the phases of mass transfer in these systems. We find that because the residual outer hydrogen envelope accounts for a large fraction of an sdB star's radius, sdB stars can spend a significant amount of time ($\sim$10s of Myr) transferring this small amount of material at low rates ($\sim 10^{-10}$-$10^{-9}\ M_\odot\,\rm yr^{-1}$) before transitioning to a phase where the bulk of their He transfers at much faster rates ($\gtrsim 10^{-8}\ M_\odot\,\rm yr^{-1}$). These systems therefore spend a surprising amount of time with Roche-filling sdB donors at orbital periods longer than the range associated with He star models without an envelope. We predict that the envelope transfer phase should be detectable by searching for ellipsoidal modulation of Roche-filling objects with $P_{\rm orb}=30$-$100$ min and $T_{\rm eff}=20{,}000$-$30{,}000$ K, and that many ($\geq$10) such systems may be found in the Galactic plane after accounting for reddening. We also argue that many of these systems may go through a phase of He transfer that matches the signatures of AM CVn systems, and that some AM CVn systems associated with young stellar populations likely descend from this channel.

We present analysis of two type-I X-ray bursts observed by NuSTAR originating from the very faint transient neutron star low-mass X-ray binary GRS 1741.9-2853 during a period of outburst in May 2020. We show that the persistent emission can be modeled as an absorbed, Comptonized blackbody in addition to Fe K$\alpha$ emission which can be attributed to relativistic disk reflection. We measure a persistent bolometric, unabsorbed luminosity of $L_{\mathrm{bol}}=7.03^{+0.04}_{-0.05}\times10^{36}\,\mathrm{erg\,s^{-1}}$, assuming a distance of 7 kpc, corresponding to an Eddington ratio of $4.5\%$. This persistent luminosity combined with light curve analysis leads us to infer that the bursts were the result of pure He burning rather than mixed H/He burning. Time-resolved spectroscopy reveals that the bolometric flux of the first burst exhibits a double-peaked structure, placing the source within a small population of accreting neutron stars which exhibit multiple-peaked type-I X-ray bursts. We find that the second, brighter burst shows evidence for photospheric radius expansion (PRE) and that at its peak, this PRE event had an unabsorbed bolometric flux of $F_{\mathrm{peak}}=2.94^{+0.28}_{-0.26}\times10^{-8}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$. This yields a new distance estimate of $d=9.0\pm0.5$ kpc, assuming that this corresponds to the Eddington limit for pure He burning on the surface of a canonical neutron star. Additionally, we performed a detailed timing analysis which failed to find evidence for quasiperiodic oscillations or burst oscillations, and we place an upper limit of $16\%$ on the rms variability around 589 Hz, the frequency at which oscillations have previously been reported.

We employ MUSE/VLT data to study the ionised and highly ionised gas phases of the feedback in Circinus, the closest Seyfert 2 galaxy to us. The analysis of the nebular emission allowed us to detect a remarkable high-ionisation gas outflow beyond the galaxy plane traced by the coronal lines [Fe VII] $\lambda$6089 and [Fe X] $\lambda$6374, extending up to 700 pc and 350 pc NW from the nucleus, respectively. This is the first time that the [Fe X] emission is observed at such distances from the central engine in an AGN. The gas kinematics reveals expanding gas shells with velocities of a few hundred km s$^{-1}$, spatially coincident with prominent hard X-ray emission detected by Chandra. Density and temperature sensitive line ratios show that the extended high-ionisation gas is characterized by a temperature reaching 25000 K and an electron density > 10$^2$ cm$^{-3}$. We found that local gas excitation by shocks produced by the passage of a radio jet leads to the spectacular high-ionisation emission in this object. This hypothesis is fully supported by photoionisation models that accounts for the combined effects of the central engine and shocks. They reproduce the observed emission line spectrum at different locations inside and outside of the NW ionisation cone. The energetic outflow produced by the radio jet is spatially located close to an extended molecular outflow recently reported using ALMA which suggests that they both represent different phases of the same feedback process acting on the AGN.

A zero point calibration of the Red Giant Branch Tip (TRGB) in the $I$-band is determined from OGLE photometry of the Magellanic Clouds (MCs). It is shown that TRGB measurements made in star-forming regions, with concomitantly high quantities of gas and dust, are less precise and biased to fainter magnitudes, as compared to the same measurements made in quiescent regions. Once these low accuracy fields are excluded from consideration, the TRGB can be used for the first time to constrain the three-dimensional plane geometry of the LMC. Composite CMDs are constructed for the SMC and LMC from only those fields with well-defined TRGB features, and the highest accuracy TRGB zero point calibration to date is presented. The $I$-band TRGB magnitude is measured to be flat over the color range $ 1.45 < (V-I)_0 < 1.95$ mag, with a modest slope introduced when including metal-rich (up to $(V-I)_0 = 2.2$ mag) Tip stars into the fit. Both the flat, blue zero point and the shallow slope calibration are consistent with the canonical value of $-4.05$ mag for the old, metal-poor TRGB, and would appear to resolve a recent debate in the literature over the method's absolute calibration.

The investigation of the wind in the solar corona initiated with the observations of the resonantly scattered UV emission of the coronal plasma obtained with UVCS-SOHO, designed to measure the wind outflow speed by applying the Doppler dimming diagnostics. Metis on Solar Orbiter complements the UVCS spectroscopic observations, performed during solar activity cycle 23, by simultaneously imaging the polarized visible light and the HI Ly-alpha corona in order to obtain high-spatial and temporal resolution maps of the outward velocity of the continuously expanding solar atmosphere. The Metis observations, on May 15, 2020, provide the first HI Ly-alpha images of the extended corona and the first instantaneous map of the speed of the coronal plasma outflows during the minimum of solar activity and allow us to identify the layer where the slow wind flow is observed. The polarized visible light (580-640 nm), and the UV HI Ly-alpha (121.6 nm) coronal emissions, obtained with the two Metis channels, are combined in order to measure the dimming of the UV emission relative to a static corona. This effect is caused by the outward motion of the coronal plasma along the direction of incidence of the chromospheric photons on the coronal neutral hydrogen. The plasma outflow velocity is then derived as a function of the measured Doppler dimming. The static corona UV emission is simulated on the basis of the plasma electron density inferred from the polarized visible light. This study leads to the identification, in the velocity maps of the solar corona, of the high-density layer about +/-10 deg wide, centered on the extension of a quiet equatorial streamer present at the East limb where the slowest wind flows at about (160 +/- 18) km/s from 4 Rs to 6 Rs. Beyond the boundaries of the high-density layer, the wind velocity rapidly increases, marking the transition between slow and fast wind in the corona.

We present a model where some proportion of extraterrestrial civilizations expand uniformly over time to reach a cosmological scale. We then ask what humanity could infer if a sky survey were to find zero, one, or more such civilizations. We show how the results of this survey, in combination with an approach to anthropics called the Self Indication Assumption (SIA), would shift any prior estimates of two quantities: 1) The chance that a technological civilization like ours survives to embark on such expansion, and 2) the maximum feasible speed at which it could expand. The SIA gives pessimistic estimates for both, but survey results (even null results) can reverse some of the effect.

We present high angular resolution imaging observations of 517 host stars of TESS exoplanet candidates using the `Alopeke and Zorro speckle cameras at Gemini North and South. The sample consists mainly of bright F, G, K stars at distances of less than 500 pc. Our speckle observations span angular resolutions of ~20 mas out to 1.2 arcsec, yielding spatial resolutions of <10 to 500 AU for most stars, and our contrast limits can detect companion stars 5-9 magnitudes fainter than the primary at optical wavelengths. We detect 102 close stellar companions and determine the separation, magnitude difference, mass ratio, and estimated orbital period for each system. Our observations of exoplanet host star binaries reveal that they have wider separations than field binaries, with a mean orbital semi-major axis near 100 AU. Other imaging studies have suggested this dearth of very closely separated binaries in systems which host exoplanets, but incompleteness at small separations makes it difficult to disentangle unobserved companions from a true lack of companions. With our improved angular resolution and sensitivity, we confirm that this lack of close exoplanet host binaries is indeed real. We also search for a correlation between planetary orbital radii vs. binary star separation, but given the very short orbital periods of the TESS planets, we do not find any clear trend. We do note that in exoplanet systems containing binary host stars, there is an observational bias against detecting Earth-size planet transits due to transit depth dilution caused by the companion star.

We present an innovative and widely applicable approach for the detection and classification of stellar clusters, developed for the PHANGS-HST Treasury Program, an $NUV$-to-$I$ band imaging campaign of 38 spiral galaxies. Our pipeline first generates a unified master source list for stars and candidate clusters, to enable a self-consistent inventory of all star formation products. To distinguish cluster candidates from stars, we introduce the Multiple Concentration Index (MCI) parameter, and measure inner and outer MCIs to probe morphology in more detail than with a single, standard concentration index (CI). We improve upon cluster candidate selection, jointly basing our criteria on expectations for MCI derived from synthetic cluster populations and published cluster catalogues, yielding model and empirical selection regions (respectively). Selection purity (confirmed clusters versus candidates, assessed via human-based classification) is high (up to 70\%) for moderately luminous sources in the empirical selection region, and somewhat lower overall (outside the region or fainter). The number of candidates rises steeply with decreasing luminosity, but pipeline-integrated Machine Learning (ML) classification prevents this from being problematic. We quantify the performance of our PHANGS-HST methods in comparison to LEGUS for a sample of four galaxies in common to both surveys, finding overall agreement with 50-75% of human verified star clusters appearing in both catalogues, but also subtle differences attributable to specific choices adopted by each project. The PHANGS-HST ML-classified Class 1 or 2 catalogues reach $\sim1$ magnitude fainter, $\sim2\times$ lower stellar mass, and are $2{-}5\times$ larger in number, than attained in the human classified samples.

Highly r-process enhanced metal-poor stars (MP r-II, $\rm [Eu/Fe]>1$ and $\rm [Fe/H]\lesssim-1.5$) have been observed in ultra-faint dwarf (UFD) galaxy, specifically in Reticulum~II (Ret~II). The fact that only a few UFDs contain such stars implies that the r-process site may reflect very rare, but individually prolific events, such as neutron star mergers (NSMs). Considering the relatively short star formation history (SFH) of UFDs, it is puzzling how they could experience such a rare phenomenon. In this work, we show the results of cosmological hydrodynamic zoom-in simulations of isolated UFDs ($M_{vir}\approx10^7-10^8$ solar mass and $M_{\ast}\approx10^3-10^4$ solar mass at $z=0$) to explain the formation of MP r-II stars in UFDs. We employ a simple toy model for NSM events, adopting parameters consistent with observations, such as the NSM rate (1 per $M_{\ast}\approx10^5$ solar mass) and europium (Eu) mass ($M_{Eu}\approx10^{-5}$ solar mass). We identify only one simulated galaxy ($ M_{vir}\approx4.6\times10^7$, $M_{\ast}\approx 3.4\times 10^3$ solar mass at $z=0$) with abundances similar to Ret~II in a simulation volume that hosts $\sim30$ UFD analogs, indicating that such abundances are possible but rare. By exploring a range of key parameters, we demonstrate that the most important factor in determining the formation of MP r-II stars in UFDs is how quickly subsequent stars can be formed out of r-process enriched gas. We find that it takes between 10 to 100~Myr to form the first and second burst of MP r-II stars. Over this period, Eu-polluted gas maintains the required high abundance ratios of $\rm [Eu/Fe]>1$.

Continued follow-up of WISEA J153429.75-104303.3, announced in Meisner et al (2020), has proven it to have an unusual set of properties. New imaging data from Keck/MOSFIRE and HST/WFC3 show that this object is one of the few faint proper motion sources known with J-ch2 > 8 mag, indicating a very cold temperature consistent with the latest known Y dwarfs. Despite this, it has W1-W2 and ch1-ch2 colors ~1.6 mag bluer than a typical Y dwarf. A new trigonometric parallax measurement from a combination of WISE, Spitzer, and HST astrometry confirms a nearby distance of $16.3^{+1.4}_{-1.2}$ pc and a large transverse velocity of $207.4{\pm}15.9$ km/s. The absolute J, W2, and ch2 magnitudes are in line with the coldest known Y dwarfs, despite the highly discrepant W1-W2 and ch1-ch2 colors. We explore possible reasons for the unique traits of this object and conclude that it is most likely an old, metal-poor brown dwarf and possibly the first Y subdwarf. Given that the object has an HST F110W magnitude of 24.7 mag, broad-band spectroscopy and photometry from JWST are the best options for testing this hypothesis.

We report on detailed analysis of the hard X-ray and GeV gamma-ray spectra of LS 5039, one of the brightest gamma-ray binary system in the Galaxy. The NuSTAR observation covering its entire orbit in 2016 allowed us for the first time to study the orbital variability of the spectrum above 10 keV. The hard X-ray spectrum is well described with a single power-law component up to 78 keV. The X-ray flux showed a slight deviation from those observed previously with Suzaku in 2007. The fast X-ray brightening observed with Suzaku, around the inferior conjunction, was not observed in this observation. We also analyzed 11 years of Fermi Large Area Telescope data of LS 5039. The GeV spectrum around the inferior conjunction was well described with two non-thermal components; a power law with a photon index of $\sim 3$ and a cut-off power law with a cutoff energy of $\sim 2$ GeV. The orbital flux variability also changed gradually around a few GeV. These results indicate that there are two emission components in the GeV band, and the dominant component above $\sim 1$ GeV does not depend on the orbital phase. By combining these results, we update the spectral energy distribution of LS 5039 with the highest available statistics. Theoretical models proposed so far cannot explain the obtained multi-wavelength spectrum, especially the emission from $\sim$ 1 MeV to $\sim$ 400 MeV, and we discuss a possibility that particle acceleration in LS 5039 is different from the shock acceleration.

Recently Kipping (2021) identified the so-called "exomoon corridor", a potentially powerful new tool for identifying possible exomoon hosts, enabled by the observation that fully half of all planets hosting an exomoon will exhibit transit timing variation (TTV) periodicities of 2-4 epochs. One key outstanding problem in the search for exomoons, however, is the question of how well the methods we have developed under the single moon assumption extend to systems with multiple moons. In this work we use $N$-body simulations to examine the exomoon corridor effect in the more general case of $N \geq 1$ moons, generating realistic TTVs produced by satellite systems more akin to those seen in the outer Solar System. We find that indeed the relationship does hold for systems with up to 5 moons in both resonant and non-resonant chain configurations. Our results suggest an observational bias against finding systems with large numbers of massive moons; as the number of moons increases, total satellite mass ratios are generally required to be significantly lower in order to maintain stability, or architectures must be more finely tuned to survive. Moons produced in impact or capture scenarios may therefore dominate early detections. Finally, we examine the distribution of TTV periods measured for a large number of Kepler Objects of Interest (KOIs) and find the same characteristic exomoon corridor distribution in several cases. This could be dynamical evidence for an abundance of moons in the field, though we caution against strong inferences based on this result.

Deuterated molecules are good tracers of the evolutionary stage of star-forming cores. During the star formation process, deuterated molecules are expected to be enhanced in cold, dense pre-stellar cores and to deplete after protostellar birth. In this paper we study the deuteration fraction of formaldehyde in high-mass star-forming cores at different evolutionary stages to investigate whether the deuteration fraction of formaldehyde can be used as an evolutionary tracer. Using the APEX SEPIA Band 5 receiver, we extended our pilot study of the $J$=3$\rightarrow$2 rotational lines of HDCO and D$_2$CO to eleven high-mass star-forming regions that host objects at different evolutionary stages. High-resolution follow-up observations of eight objects in ALMA Band 6 were performed to reveal the size of the H$_2$CO emission and to give an estimate of the deuteration fractions HDCO/H$_2$CO and D$_2$CO/HDCO at scales of $\sim$6" (0.04-0.15 pc at the distance of our targets). Our observations show that singly- and doubly deuterated H$_2$CO are detected toward high-mass protostellar objects (HMPOs) and ultracompact HII regions (UCHII regions), the deuteration fraction of H$_2$CO is also found to decrease by an order of magnitude from the earlier HMPO phases to the latest evolutionary stage (UCHII), from $\sim$0.13 to $\sim$0.01. We have not detected HDCO and D$_2$CO emission from the youngest sources (high-mass starless cores, HMSCs). Our extended study supports the results of the previous pilot study: the deuteration fraction of formaldehyde decreases with evolutionary stage, but higher sensitivity observations are needed to provide more stringent constraints on the D/H ratio during the HMSC phase. The calculated upper limits for the HMSC sources are high, so the trend between HMSC and HMPO phases cannot be constrained.

We report the new detection of $^7$Be II in the ultraviolet spectra of V5669 Sgr during its early decline phase ($+24$ and $+28$ d). We identified three blue-shifted absorption systems in our spectra. The first two, referred to as low- and high-velocity components, were noticeably identified among H I Balmer, Na I D, and Fe II whose lower energies of transients are low ($<4$ eV). The third absorption component was identified among N II, He I, and C II lines whose lower energy levels are relatively high (9--21 eV). The absorption lines of $^7$Be II at $3130.583$ {\AA}, and $3132.228$ {\AA} were identified as the first and second components in our observations. No evidence suggested the existence of Li I at 6708 {\AA} in any velocity components. The estimated number density ratio of lithium relative to hydrogen, which was finally produced by this object using the equivalent widths of $^7$Be and Ca II K, $N({\rm ^{7}Li})/N({\rm H})_{\rm final}$ is $4.0\pm0.7\times10^{-6}$. This value is an order of magnitude lower than the average observed values for classical novae wherein $^7$Be has been detected, and is comparable to the most optimistic value of theoretical predictions.

We searched for X-ray supernova remnants (SNRs) in the starburst region of M82, using archival data from the Chandra X-ray Observatory with a total effective exposure time of 620 ks with an X-ray spectroscopic selection. Strong line-emission from Fe xxv at 6.7 keV is a characteristic spectral feature of hot, shocked gas of young SNRs and distinctive among the discrete sources in the region populated by X-ray binaries. We selected candidates using narrow-band imaging aimed at the line excess and identified six (and possibly a seventh) X-ray SNRs. Two previously known examples were recovered by our selection. Five of them have radio counterparts, including the radio supernova SN2008iz, which was discovered as a radio transient in 2008. It shows a hard X-ray spectrum with a blueshifted Fe K feature with v ~ -2700 km/s, both of which suggest its youth. The 4-8 keV luminosities of the selected SNRs are in the range of (0.3-3)e38 erg/s. We made a crude estimate of the supernova rate, assuming that more luminous SNRs are younger, and found 0.06 (0.03-0.13) /yr, in agreement with the supernova rates estimated by radio observations and the generally believed star formation rate of M82, although the validity of the assumption is questionable. A sum of the Fe xxv luminosity originating from the selected X-ray SNRs consists of half of the total Fe xxv luminosity observed in the central region of M82. We briefly discuss its implications for starburst winds and the Fe xxv emission in more luminous starburst galaxies.

A rotating black hole can be clouded by light bosons via superradiance, and thus acquire an atom-like structure. If such a gravitational atom system is companioned with a pulsar, the pulsar can trigger transitions between energy levels of the gravitational atom, and these transitions can be detected by pulsar timing. We show that in such pulsar-black hole systems, fine and hyperfine structure transitions are more likely to be probed than the Bohr transition. Also, the calculation of these fine and hyperfine structure transitions are under better analytic control. Thus, these fine and hyperfine structure transitions are more ideal probes in the search for gravitational collider signals in pulsar-black hole systems.

Electromagnetic waves due to electron-positron clouds (bunches), created by cascading processes in pulsar magnetospheres, have been proposed to explain the pulsar radio emission. In order to verify this hypothesis, we utilized for the first time Particle-in-Cell (PIC-) code simulations to study the nonlinear evolution of electron-positron bunches in dependence on the relative drift speeds of electrons and positrons, on the initial plasma temperature, and on the initial distance between the bunches. For this sake, we utilized the PIC-code ACRONYM with a high-order field solver and particle weighting factor, appropriate to describe relativistic pair plasmas. We found that the bunch expansion is mainly determined by the relative electron-positron drift speed. Finite drift speeds were found to cause the generation of strong electric fields that reach up to $E \sim 7.5 \times 10^{5}$ V/cm ($E / (m_\mathrm{e} c \omega_\mathrm{p} e^{-1}) \sim 4.4$) and strong plasma heating. As a result, up to 15~\% of the initial kinetic energy is transformed into the electric field energy. Assuming the same electron- and positron-distributions we found that the fastest (in the bunch reference frame) particles of consecutively emitted bunches eventually overlap in the momentum (velocity) space. This overlap causes two-stream instabilities that generate (electrostatic) subluminal L-mode waves with electric field amplitudes reaching up to $E \sim 1.9\times 10^{4}$ V/cm ($E / (m_\mathrm{e} c \omega_\mathrm{p} e^{-1}) \sim 0.11$). We found that the interaction of electron-position bunches leads to plasma heating, to the generation of strong electric fields and of intense superluminal L-mode waves which, in principle, can be behind the observed electromagnetic emissions of pulsars in the radio wave range.

Solar activity is an important driver of long-term climate trends and must be accounted for in climate models. Unfortunately, direct measurements of this quantity over long periods do not exist. The only observation related to solar activity whose records reach back to the seventeenth century are sunspots. Surprisingly, determining the number of sunspots consistently over time has remained until today a challenging statistical problem. It arises from the need of consolidating data from multiple observing stations around the world in a context of low signal-to-noise ratios, non-stationarity, missing data, non-standard distributions and many kinds of errors. The data from some stations experience therefore severe and various deviations over time. In this paper, we propose the first systematic and thorough statistical approach for monitoring these complex and important series. It consists of three steps essential for successful treatment of the data: smoothing on multiple timescales, monitoring using block bootstrap calibrated CUSUM charts and classifying of out-of-control situations by support vector techniques. This approach allows us to detect a wide range of anomalies (such as sudden jumps or more progressive drifts), unseen in previous analyses. It helps us to identify the causes of major deviations, which are often observer or equipment related. Their detection and identification will contribute to improve future observations. Their elimination or correction in past data will lead to a more precise reconstruction of the world reference index for solar activity: the International Sunspot Number.

We outline a new method suggested by Conway (2016) for solving the two-body problem for solid bodies of spheroidal or ellipsoidal shape. The method is based on integrating the gravitational potential of one body over the surface of the other body. When the gravitational potential can be analytically expressed (as for spheroids or ellipsoids), the gravitational force and mutual gravitational potential can be formulated as a surface integral instead of a volume integral, and solved numerically. If the two bodies are infinitely thin disks, the surface integral has an analytical solution. The method is exact as the force and mutual potential appear in closed-form expressions, and does not involve series expansions with subsequent truncation errors. In order to test the method, we solve the equations of motion in an inertial frame, and run simulations with two spheroids and two infinitely thin disks, restricted to torque-free planar motion. The resulting trajectories display precession patterns typical for non-Keplerian potentials. We follow the conservation of energy and orbital angular momentum, and also investigate how the spheroid model approaches the two cases where the surface integral can be solved analytically, i.e. for point masses and infinitely thin disks.

Context. As the importance of Gravitational Wave (GW) Astrophysics increases rapidly, astronomers in different fields and with different backgrounds can have the need to get a quick idea of which GW source populations can be detected by which detectors and with what measurement uncertainties. Aims. The GW-Toolbox is an easy-to-use, flexible tool to simulate observations on the GW universe with different detectors, including ground-based interferometers (advanced LIGO, advanced VIRGO, KAGRA, Einstein Telescope, and also customised designs), space-borne interferometers (LISA and a customised design), pulsar timing arrays mimicking the current working ones (EPTA, PPTA, NANOGrav, IPTA) and future ones. We include a broad range of sources such as mergers of stellar mass compact objects, namely black holes, neutron stars and black hole-neutron stars; and supermassive black hole binaries mergers and inspirals, Galactic double white dwarfs in ultra-compact orbit, extreme mass ratio inspirals and Stochastic GW backgrounds. Methods. We collect methods to simulate source populations and determine their detectability with the various detectors. The paper aims at giving a comprehensive description on the algorithm and functionality of the GW-Toolbox. Results. The GW-Toolbox produces results that are consistent with more detailed calculations of the different source classes and can be accessed with a website interface (gw-universe.org) or as a python package (https://bitbucket.org/radboudradiolab/gwtoolbox). In the future, it will be upgraded with more functionality.

Recent large-scale spectropolarimetric surveys have established that a small but significant percentage of massive stars host stable, surface dipolar magnetic fields with strengths on the order of kG. These fields channel the dense, radiatively driven stellar wind into circumstellar magnetospheres, whose density and velocity structure can be probed using ultraviolet (UV) spectroscopy of wind-sensitive resonance lines. Coupled with appropriate magnetosphere models, UV spectroscopy provides a valuable way to investigate the wind-field interaction, and can yield quantitative estimates of the wind parameters of magnetic massive stars. We report a systematic investigation of the formation of UV resonance lines in slowly rotating magnetic massive stars with dynamical magnetospheres. We pair the Analytic Dynamical Magnetosphere (ADM) formalism with a simplified radiative transfer technique to produce synthetic UV line profiles. Using a grid of models, we examine the effect of magnetosphere size, the line strength parameter, and the cooling parameter on the structure and modulation of the line profile. We find that magnetic massive stars uniquely exhibit redshifted absorption at most viewing angles and magnetosphere sizes, and that significant changes to the shape and variation of the line profile with varying line strengths can be explained by examining the individual wind components described in the ADM formalism. Finally, we show that the cooling parameter has a negligible effect on the line profiles.

We analysed the chemodynamical evolution of the Galactic disc using precise [Mg/Fe] abundances from a previous study and accurate Gaia data. For this purpose, we estimated ages and dynamical properties for 366 MSTO solar neighbourhood stars from the AMBRE Project using PARSEC isochrones together with astrometric and photometric values from Gaia DR2. We find a radial gradient of -0.099 ${\pm}$ 0.031 dex kpc$^{-1}$ for [M/H] and +0.023 ${\pm}$ 0.009 dex kpc for the [Mg/Fe] abundance. The steeper [Mg/Fe] gradient than that found in the literature is a result of the improvement of the AMBRE [Mg/Fe] estimates in the metal-rich regime. In addition, we find a significant spread of stellar age at any given [Mg/Fe] value, and observe a clear correlated dispersion of the [Mg/Fe] abundance with metallicity at a given age. While for [M/H] < -0.2, a clear age-[Mg/Fe] trend is observed, more metal-rich stars display ages from 3 up to 12 Gyr, describing an almost flat trend in the [Mg/Fe]-age relation. Moreover, we report the presence of radially migrated stars for a wide range of stellar ages, although we note the large uncertainties of the amplitude of the inferred change in orbital guiding radii. Finally, we observe the appearance of a second chemical sequence in the outer disc, 10-12 Gyr ago, populating the metal-poor, low-[Mg/Fe] tail. These stars are more metal-poor than the coexisting stellar population in the inner parts of the disc, and show lower [Mg/Fe] abundances than prior disc stars of the same metallicity, leading to a chemical discontinuity. Our data favour the rapid formation of an early disc that settled in the inner regions, followed by the accretion of external metal-poor gas -- probably related to a major accretion event such as the Gaia-Enceladus/Sausage one -- that may have triggered the formation of the thin disc population and steepened the abundance gradient in the early disc.

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected (-23.0<M_{UV,AB}<-21.3) galaxies identified over a 7-deg**2 area (including the wide-area COSMOS/UltraVISTA, VIDEO/XMM-LSS, and UKIDSS/UDS fields) for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~81% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 20 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization.

We investigate the large-scale clustering of the final spectroscopic sample of quasars from the recently completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). The sample contains $343708$ objects in the redshift range $0.8<z<2.2$ and $72667$ objects with redshifts $2.2<z<3.5$, covering an effective area of $4699~{\rm deg}^{2}$. We develop a neural network-based approach to mitigate spurious fluctuations in the density field caused by spatial variations in the quality of the imaging data used to select targets for follow-up spectroscopy. Simulations are used with the same angular and radial distributions as the real data to estimate covariance matrices, perform error analyses, and assess residual systematic uncertainties. We measure the mean density contrast and cross-correlations of the eBOSS quasars against maps of potential sources of imaging systematics to address algorithm effectiveness, finding that the neural network-based approach outperforms standard linear regression. Stellar density is one of the most important sources of spurious fluctuations, and a new template constructed using data from the Gaia spacecraft provides the best match to the observed quasar clustering. The end-product from this work is a new value-added quasar catalogue with the improved weights to correct for nonlinear imaging systematic effects, which will be made public. Our quasar catalogue is used to measure the local-type primordial non-Gaussianity in our companion paper, Mueller et al. in preparation.

We present measurements of the local primordial non-Gaussianity parameter \fNLloc from the clustering of 343,708 quasars with redshifts 0.8 < z < 2.2 distributed over 4808 square degrees from the final data release (DR16) of the extended Baryon acoustic Oscillation Spectroscopic Survey (eBOSS), the largest volume spectroscopic survey up to date. Our analysis is performed in Fourier space, using the power spectrum monopole at very large scales to constrain the scale dependent halo bias. We carefully assess the impact of systematics on our measurement and test multiple contamination removal methods. We demonstrate the robustness of our analysis pipeline with EZ-mock catalogues that simulate the eBOSS DR16 target selection. We find $f_\mathrm{NL}=-12\pm 21$ (68\% confidence) for the main clustering sample including quasars with redshifts between 0.8 and 2.2, after exploiting a novel neural network scheme for cleaning the DR16 sample and in particular after applying redshift weighting techniques, designed for non-Gaussianity measurement from large scales structure, to optimize our analysis, which improve our results by 37\%.

The next generation of gravitational-wave experiments, such as Einstein Telescope, Cosmic Explorer and LISA, will test the primordial black hole scenario. We provide a forecast for the minimum testable value of the abundance of primordial black holes as a function of their masses for both the unclustered and clustered spatial distributions at formation. In particular, we show that these instruments may test abundances, relative to the dark matter, as low as $10^{-10}$.

We examine whether the Large Magellanic Cloud (LMC) is currently losing its stellar halo to Milky Way (MW) tides. We present a live $N$-body model for the ongoing MW-LMC interaction that predicts a prominent stream of stars tidally stripped from the progenitor LMC. We use this model to define a strategy to search for stripped material in kinematic space. Of the available stellar tracers, we conclude that samples of RR Lyrae stars provide the highest density of kinematic tracers at present. Using a sample of RR Lyrae stars with Gaia EDR3 astrometry we show that the LMC stellar halo in the Southern Galactic hemisphere extends at least out to $\sim 30^\circ$ from the galaxy centre. In addition, several leading arm candidates are found in the Northern hemisphere as far above the disc plane as $b=+34^\circ$ (at 68$^\circ$ from the LMC).

We introduce an Eulerian Perturbation Theory to study the clustering of tracers for cosmologies in the presence of massive neutrinos. Our approach is based on mapping recently-obtained Lagrangian Perturbation Theory results to the Eulerian framework. We add Effective Field Theory counterterms, IR-resummations and a biasing scheme to compute the one-loop redshift-space power spectrum. To assess our predictions, we compare the power spectrum multipoles against synthetic halo catalogues from the Quijote simulations, finding excellent agreement on scales $k\lesssim 0.25 \,h \text{Mpc}^{-1}$. Extending the range of accuracy to higher wave-numbers is possible at the cost of producing an offset in the best-fit linear local bias. We further discuss the implications for the tree-level bispectrum. Finally, calculating loop corrections is computationally costly, hence we derive an accurate approximation wherein we retain only the main features of the kernels, as produced by changes to the growth rate. As a result, we show how FFTLog methods can be used to further accelerate the loop computations with these reduced kernels.

In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated)

Time series analysis is ubiquitous in many fields of science including gravitational-wave astronomy, where strain time series are analyzed to infer the nature of gravitational-wave sources, e.g., black holes and neutron stars. It is common in gravitational-wave transient studies to apply a tapered window function to reduce the effects of spectral artifacts from the sharp edges of data segments. We show that the conventional analysis of tapered data fails to take into account covariance between frequency bins, which arises for all finite time series -- no matter the choice of window function. We discuss the origin of this covariance and show that as the number of gravitational-wave detections grows, and as we gain access to more high signal-to-noise ratio events, this covariance will become a non-negligible source of systematic error. We derive a framework that models the correlation induced by the window function and demonstrate this solution using both data from the first LIGO--Virgo transient catalog and simulated Gaussian noise.

The Einasto profile has been successful in describing the density profiles of dark matter haloes in $\Lambda$CDM N-body simulations. It has also been able to describe multiple components in the surface brightness profiles of galaxies. However, analytically projecting it to calculate quantities under projection is challenging. In this paper, we will see the development of a highly accurate analytical approximation for the mass (or counts) enclosed in an infinitely long cylindrical column for Einasto profiles--also known as the projected mass (or counts)--using a novel methodology. We will then develop a self-consistent high-accuracy model for the surface density from the expression for the projected mass. Both models are quite accurate for a broad family of functions, with a shape parameter $\alpha$ varying by a factor of 100 in the range $0.05 \lesssim \alpha \lesssim 5.0$, with fractional errors $\sim 10^{-6}$ for $\alpha \lesssim 0.4$. Profiles with $\alpha \lesssim 0.4$ have been shown to fit the density profiles of dark matter haloes in N-body simulations as well as the luminosity profiles of the outer components of massive galaxies. Since the projected mass and the surface density are used in gravitational lensing, I will illustrate how these models facilitate (for the first time) analytical computation of several quantities of interest in lensing due to Einasto profiles. The models, however, are not limited to lensing and apply to similar quantities under projection, such as the projected luminosity, the projected (columnar) number counts and the projected density or the surface brightness.

We present a study of the relative orientation between the magnetic field and elongated cloud structures for the $\rho$ Oph A and $\rho$ Oph E regions in L1688 in the Ophiuchus molecular cloud. Combining inferred magnetic field orientation from HAWC+ 154 $\mu$m observations of polarized thermal emission with column density maps created using Herschel submillimeter observations, we find consistent perpendicular relative alignment at scales of $0.02$ pc ($33.6"$ at $d \approx 137$ pc) using the histogram of relative orientations (HRO) technique. This supports the conclusions of previous work using Planck polarimetry and extends the results to higher column densities. Combining this HAWC+ HRO analysis with a new Planck HRO analysis of L1688, the transition from parallel to perpendicular alignment in L1688 is observed to occur at a molecular hydrogen column density of approximately $10^{21.7}$ cm$^{-2}$. This value for the alignment transition column density agrees well with values found for nearby clouds via previous studies using only Planck observations. Using existing turbulent, magnetohydrodynamic simulations of molecular clouds formed by colliding flows as a model for L1688, we conclude that the molecular hydrogen volume density associated with this transition is approximately $\sim10^{4}$ cm$^{-3}$. We discuss the limitations of our analysis, including incomplete sampling of the dense regions in L1688 by HAWC+.

Teleparallel gravity has significantly increased in popularity in recent decades, bringing attention to Einstein's other theory of gravity. In this Review, we relate this form of geometry to the broader metric-affine approach to forming gravitational theories where we describe a systematic way of constructing consistent teleparallel theories that respect certain physical conditions such as local Lorentz invariance. We first use teleparallel gravity to formulate a teleparallel equivalent of general relativity which is dynamically equivalent to general relativity but which may have different behaviors for other scenarios, such as quantum gravity. After setting this foundation, we describe the plethora of modified teleparallel theories of gravity that have been proposed in the literature. In the second part of the Review, we first survey works in teleparallel astrophysics literature where we focus on the open questions in this regime of physics. We then discuss the cosmological consequences for the various formulations of teleparallel gravity. We do this at background level by exploring works using various approaches ranging from dynamical systems to Noether symmetries, and more. Naturally, we then discuss perturbation theory, firstly by giving a concise approach in which this can be applied in teleparallel gravity theories and then apply it to a number of important theories in the literature. Finally, we examine works in observational and precision cosmology across the plethora of proposal theories. This is done using some of the latest observations and is used to tackle cosmological tensions which may be alleviated in teleparallel cosmology. We also introduce a number of recent works in the application of machine learning to gravity, we do this through deep learning and Gaussian processes, together with discussions about other approaches in the literature.

We present an archaeoastronomical study of the orientations of the colonial Christian churches on the island of Fuerteventura, in the Canary Islands, Spain, mostly built from the period of the Norman conquest in the 15th century to the 19th century. Our goal is to analyze the possible astronomical influence on the orientation of these churches. Preliminary results suggest that the vast majority of the island's religious constructions have their axes oriented within the solar range, between the extreme azimuths of the annual movement of the Sun as it crosses the local horizon. This differs from what was found on the islands of Lanzarote and La Gomera (also in the Canaries) previously studied.