Papers for Friday, Apr 01 2022

Minor mergers are common in galaxy clusters. They have the potential to create sloshing cold fronts (SCFs) in the intracluster medium (ICM) of the cluster. However, the resilience of SCFs to subsequent minor mergers is unknown. Here we investigate the extent to which SCFs established by an off-axis minor merger are disrupted by a subsequent minor merger. We perform a suite of 13 hydrodynamic + N-body simulations of idealised triple cluster mergers in which we vary the approach direction and impact parameter of the tertiary cluster. Except for ~1 Gyr after the first core passage of the tertiary cluster, clear SCFs are present in all merger configurations. Subsequent head-on minor mergers reduce the number of SCFs significantly, while subsequent off-axis minor mergers only moderately reduce the number of SCFs. In particular, outer (>~500 kpc) SCFs are resilient. The results of this work indicate that SCFs are easily formed in the course of a minor merger and are long-lived even if a further minor merger takes place. SCFs thus should be ubiquitous, but deriving the merger history of a given cluster based on its observed SCFs might be more complex than previously thought.

The majority of massive stars are found in close binaries which: (i) are prone to merge and (ii) are accompanied by another distant tertiary star (triples). Here, we study the evolution of the stellar post-merger binaries composed of the merger product and the tertiary companion. We find that post-merger binaries originating from compact stellar triples with outer semi-major axes $a_{\rm out,init}\lesssim10^1\,-\,10^2\,\rm AU$ provide a new way to form binary black hole mergers in the galactic field. By means of a population synthesis, we estimate their contribution to the total black hole merger rate to be $\mathcal{R}(z=0)=0.3\,-\,25.2\,\rm Gpc^{-3}\,yr^{-1}$. Merging binary black holes that form from stellar post-merger binaries have exceptionally low mass ratios. We identify a critical mass ratio $q\simeq0.5$ below which they dominate the total black hole merger rate in the field. We show that after including their additional contribution, the mass ratio distribution of binary black hole mergers in the galactic field scenario is in better agreement with that inferred from gravitational wave detections.

Using $\it{Herschel}$-SPIRE imaging and the Canada-France Imaging Survey (CFIS) Low Surface Brightness data products from the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS), we present a cross-correlation between the cosmic infrared background and cosmic optical background fluctuations. With a combined sky area of $91\,{\rm deg}^2$ and a minimum resolved scale of $18\,$arcsec, the cross-spectrum is measured for two cases: all galaxies are kept in the images; or all individually-detected galaxies are masked to produce `background' maps. We report the detection of the cross-correlation signal at $\gtrsim 22\,\sigma$ ($\gtrsim 18\,\sigma$ for the background map). The part of the optical brightness variations that are correlated with the submm emission translates to an rms brightness of $\simeq 32.5\,{\rm mag}\,{\rm arcsec}^{-2}$ in the $r$ band, a level normally unreachable for individual sources. A critical issue is determining what fraction of the cross-power spectrum might be caused by emission from Galactic cirrus. For one of the fields, the Galactic contamination is approximately a factor of 10 higher than the extragalactic signal, with the contamination being estimated using a linear regression from several external survey maps; however, for the other fields, the contamination is typically around 20 per cent. An additional discriminant is that the cross-power spectrum is of the approximate form $P(k)\propto 1/k$, much shallower than that of Galactic cirrus. We interpret the results in a halo-model framework, which shows good agreement with independent measurements for the scalings of star-formation rates in galaxies. The approach presented in this study holds great promise for future surveys such as FYST/CCAT-prime combined with ${\it Euclid}$ or the Vera Rubin Observatory (LSST), which will enable a detailed exploration of the evolution of star formation in galaxies.

The current generation of galaxy simulations can resolve individual giant molecular clouds, the progenitors of dense star clusters. But the evolutionary fate of these young massive clusters (YMCs), and whether they can become the old globular clusters (GCs) observed in many galaxies, is determined by a complex interplay of internal dynamical processes and external galactic effects. We present the first star-by-star $N$-body models of massive ($N\sim10^5-10^7$) star clusters formed in a FIRE-2 MHD simulation of a Milky Way-mass galaxy, with all of the relevant initial conditions and galactic tidal effects extracted directly from the cosmological simulation. We randomly select 895 ($\sim 30\%$) of the YMCs with $ > 6\times10^4M_{\odot}$ from Grudi\'c et al. 2022 and integrate them to the present day using the Cluster Monte Carlo Code, CMC. This procedure predicts a MW-like system with 148 GCs, most of which were formed during the early, bursty mode of star formation in the galaxy. Our GCs are younger, less massive, and more core collapsed than clusters in the Milky Way or M31. This is a direct result of the assembly history and age-metallicity relationship of the GCs' host galaxy: younger clusters are preferentially born in stronger galactic tidal fields and initially retain fewer stellar-mass black holes, causing them to lose mass faster and reach core collapse sooner than their older counterparts. Our results suggest that the masses and core/half-light radii of GCs are shaped not only by internal dynamical processes, but by the specific evolutionary history of their host galaxies as well. These results emphasize that $N$-body studies with realistic stellar physics are crucial to understanding the evolution and present-day properties of galactic GC systems.

Using sedimentary and eclipse-based measurements of the lunar recession velocity, we derive a new local-Universe measurement of the Hubble constant ($H_0$) from the recession rate of Earth's Moon. Taking into account the effects of tides, we find a value of $H_{0}$ = 63.01 $\pm$ 1.79 km s$^{-1}$ Mpc$^{-1}$, which is in approximate agreement with the Planck space mission's measurement using the cosmic microwave background (CMB) and base $\Lambda$CDM. Our new measurement represents the first ever model-independent, single-step measurement of the Universe's current expansion rate. This is also the first major local Universe measurement of $H_0$ which is below the measurement from Planck. Importantly, it is robust to the systematic errors that may be present in other $H_0$ measurements using other cosmological probes such as type Ia supernovae, baryon acoustic oscillations, or lensed quasars. Our work provides key evidence towards the notion that the existing Hubble tension may indeed be a result of systematic uncertainties in the local distance ladder.

We present a collection of 65,981 Mira-type variable stars found in the Optical Gravitational Lensing Experiment (OGLE) project database. Two-thirds of our sample (40,356 objects) are located in the Galactic bulge fields, whereas 25,625 stars are in the Galactic disk. The vast majority of the collection (47,532 objects) are new discoveries. We provide basic observational parameters of the Mira variables: equatorial coordinates, pulsation periods, $I$-band and $V$-band mean magnitudes, $I$-band brightness amplitudes, and identifications in other catalogs of variable stars. We also provide the $I$-band and $V$-band time-series photometry collected since 1997 during the OGLE-II, OGLE-III, and OGLE-IV phases. The old-fashioned selection process, i.e., mostly based on the visual inspection of light curves by experienced astronomers, led us to the high purity of the catalog. As a result, this collection can be used as a training set in the machine learning classification algorithms. Using overlapping parts of adjacent OGLE fields, we estimate the completeness of the catalog to be about 96%. We compare and discuss the statistical features of Miras located in different regions of the Milky Way. We show examples of stars that change their type over time, from a semi-regular variable to Mira and vice versa. This dataset is perfectly suited to study the three-dimensional structure of the Milky Way, and it may help to explain the puzzle of the X-shaped bulge.

Some protostellar systems show little or no millimetre line emission of complex organics. This can be interpreted as a low abundance of these molecules, alternatively they could be present in the system but are not seen in the gas. The goal is to investigate the latter hypothesis for methanol. We will attempt to answer the question: Does the presence of a disk and optically thick dust reduce methanol emission even if methanol is abundant in the ices and gas? Using the radiative transfer code RADMC-3D, methanol emission lines from an envelope-only model and an envelope-plus-disk model are calculated and compared with each other and the observations. Methanol emission from the envelope-only model is always stronger than from the envelope-plus-disk model by at least a factor 2 as long as the disk radius is larger than 30 au (for L=8 L$_{\odot}$). In most cases, this is due to lower temperatures and, hence, the smaller amount of warm methanol inside the snow surface of the envelope-plus-disk model. The intensities drop by more than an order of magnitude for models including high mm opacity dust grains and disk radii of at least 50 au (for L=8 L$_{\odot}$) due to continuum over-subtraction. The line intensities from the envelope-only models overproduce the observations of protostars with lower methanol emission even with large dust optical depth effects. The envelope-plus-disk models can explain the bulk of the observations. However, they can only reproduce the observations of sources with high luminosities and low methanol emission when dust optical depth effects especially continuum over-subtraction in the disk becomes effective. Therefore, both the effects of disk and dust optical depth should be considered to explain the observations. In conclusion: Absence of methanol emission does not imply absence of methanol molecules in either gas or ice.

The outflows from neutrino-cooled black-hole (BH) accretion disks formed in neutron-star mergers or cores of collapsing stars are expected to be neutron-rich enough to explain a large fraction of elements created by the rapid neutron-capture (r-) process, but their precise chemical composition remains elusive. Here, we investigate the role of fast neutrino flavor conversion, motivated by the findings of our post-processing analysis that shows evidence of electron-neutrino lepton-number (ELN) crossings deep inside the disk, hence suggesting possibly non-trivial effects due to neutrino flavor mixing. We implement a parametric, dynamically self-consistent treatment of fast conversion in time-dependent simulations and examine the impact on the disk and its outflows. By activating the, otherwise inefficient, emission of heavy-lepton neutrinos, fast conversions enhance the disk cooling rates and reduce the absorption rates of electron-type neutrinos, causing a reduction of the electron fraction in the disk by 0.03-0.06 and in the ejected material by 0.01-0.03. The r-process yields are enhanced by typically no more than a factor of two, rendering the overall impact of fast conversions modest. The kilonova is prolonged as a net result of increased lanthanide opacities and enhanced radioactive heating rates. We observe only mild sensitivity to the disk mass, the condition for the onset of flavor conversion, and to the considered cases of flavor mixing. Remarkably, parametric models of flavor mixing that conserve the lepton numbers per family result in an overall smaller impact than models invoking three-flavor equipartition, often assumed in previous works.

We explore the effects of rapid rotation on the properties of neutrino-heated winds from proto-neutron stars (PNS) formed in core-collapse supernovae or neutron-star mergers by means of three-dimensional general-relativistic hydrodynamical simulations with M0 neutrino transport. We focus on conditions characteristic of a few seconds into the PNS cooling evolution when the neutrino luminosities obey $L_{\nu_e} + L_{\bar{\nu}_e} \approx 7\times 10^{51}$ erg s$^{-1}$, and over which most of the wind mass-loss will occur. After an initial transient phase, all of our models reach approximately steady-state outflow solutions with positive energies and sonic surfaces captured on the computational grid. Our non-rotating and slower-rotating models (angular velocity relative to Keplerian $\Omega/\Omega_{\rm K} \lesssim 0.4$; spin period $P \gtrsim 2$ ms) generate approximately spherically symmetric outflows with properties in good agreement with previous PNS wind studies. By contrast, our most rapidly spinning PNS solutions ($\Omega/\Omega_{\rm K} \gtrsim 0.75$; $P \approx 1$ ms) generate outflows focused in the rotational equatorial plane with much higher mass-loss rates (by over an order of magnitude), lower velocities, lower entropy, and lower asymptotic electron fractions, than otherwise similar non-rotating wind solutions. Although such rapidly spinning PNS are likely rare in nature, their atypical nucleosynthetic composition and outsized mass yields could render them important contributors of light neutron-rich nuclei compared to more common slowly rotating PNS birth. Our calculations pave the way to including the combined effects of rotation and a dynamically-important large-scale magnetic field on the wind properties within a 3D GRMHD framework.

The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a method to determine the collisional excitation conditions from [Fe II] line ratios independently from the ionisation state and find that it cannot be applied to highly-ionised models due to the influence of recombination cascades on Fe$^{+}$ level populations. When the ionisation state is artificially lowered, the line ratios (and excitation conditions) are too similar to distinguish between explosion scenarios. We investigate changes to the treatment of non-thermal energy deposition as a way to reconcile over-ionised theoretical models with observations and find that a simple work function approximation provides closer agreement with the data for sub-Mch models than a detailed Spencer-Fano treatment with widely-used cross section data. To quantify the magnitude of additional heating processes that would be required to sufficiently reduce ionisation from fast leptons, we artificially boost the rate of energy loss to free electrons. We find that the equivalent of as much as an eight times increase to the plasma loss rate would be needed to reconcile the sub-Mch model with observed spectra. Future studies could distinguish between reductions in the non-thermal ionisation rates and increased recombination rates, such as by clumping.

Dark matter models can be classified according to their impact on the properties of galaxies, including cold dark matter (CDM), warm dark matter (WDM), self-interacting dark matter (SIDM) and fuzzy dark matter (FDM). In celebration of April Fool's Day, and also of the 1-year anniversary of the start of the 2022 volcanic eruption at Fagradalsfjall here in Iceland, we explore fresh lava as a candidate for WDM specifically. We verify first hand that lava is indeed warm (exhibits free-streaming and retains temperature for several months after the eruption ends, is 1000K, sets fire to grass, can feel one's eyebrows singe at a distance of 4m) and dark once sufficiently decoupled from its source of production.

We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optical/near-infrared filters (g, r, i, z). DELVE DR2 provides point-source and automatic aperture photometry for ~2.5 billion astronomical sources with a median 5{\sigma} point-source depth of g=24.3, r=23.9, i=23.5, and z=22.8 mag. A region of ~17,000 deg^2 has been imaged in all four filters, providing four-band photometric measurements for ~618 million astronomical sources. DELVE DR2 covers more than four times the area of the previous DELVE data release and contains roughly five times as many astronomical objects. DELVE DR2 is publicly available via the NOIRLab Astro Data Lab science platform.

High-energy emission associated with star formation has been proposed as a significant source of interstellar medium (ISM) ionization in low-metallicity starbursts and an important contributor to the heating of the intergalactic medium (IGM) in the high-redshift ($z > 8$) Universe. Using Chandra observations of a sample of 30 galaxies at $D \approx$~200--450 Mpc that have high specific star-formation rates of 3--9 Gyr$^{-1}$ and metallicities near $Z \approx 0.3 Z_\odot$, we provide new measurements of the average 0.5--8 keV spectral shape and normalization per unit star-formation rate (SFR). We model the sample-combined X-ray spectrum as a combination of hot gas and high-mass X-ray binary (HMXB) populations and constrain their relative contributions. We derive scaling relations of $\log L_{\rm 0.5-8 keV}^{\rm HMXB}$/SFR $= 40.19 \pm 0.06$ and $\log L_{\rm 0.5-2 keV}^{\rm gas}$/SFR $= 39.58^{+0.17}_{-0.28}$; significantly elevated compared to local relations. The HMXB scaling is also somewhat higher than $L_{\rm 0.5-8 keV}^{\rm HMXB}$-SFR-$Z$ relations presented in the literature, potentially due to our galaxies having relatively low HMXB obscuration and young and X-ray luminous stellar populations. The elevation of the hot gas scaling relation is at the level expected for diminished attenuation due to a reduction of metals; however, we cannot conclude that an $L_{\rm 0.5-2 keV}^{\rm gas}$-SFR-$Z$ relation is driven solely by changes in ISM metal content. Finally, we present SFR-scaled spectral models (both emergent and intrinsic) that span the X-ray--to--IR band, providing new benchmarks for studies of the impact of ISM ionization and IGM heating in the early Universe.

Ultralight axion-like particles (ALPs) are compelling dark matter candidates because of their potential to resolve small-scale discrepancies between $\Lambda$CDM predictions and cosmological observations. Axion-photon coupling induces a polarization rotation in linearly polarized photons traveling through an ALP field; thus, as the local ALP dark matter field oscillates in time, distant static polarized sources will appear to oscillate with a frequency proportional to the ALP mass. We use observations of the cosmic microwave background from SPT-3G, the current receiver on the South Pole Telescope, to set upper limits on the value of the axion-photon coupling constant $g_{\phi\gamma}$ over the approximate mass range $10^{-22} - 10^{-19}$ eV, corresponding to oscillation periods from 12 hours to 100 days. For periods between 1 and 100 days ($4.7 \times 10^{-22} \text{ eV} \leq m_\phi \leq 4.7 \times 10^{-20} \text{ eV}$), where the limit is approximately constant, we set a median 95% C.L. upper limit on the amplitude of on-sky polarization rotation of 0.071 deg. Assuming that dark matter comprises a single ALP species with a local dark matter density of $0.3 \text{ GeV/cm}^3$, this corresponds to $g_{\phi\gamma} < 1.18 \times 10^{-12}\text{ GeV}^{-1} \times \left( \frac{m_{\phi}}{1.0 \times 10^{-21} \text{ eV}} \right)$. These new limits represent an improvement over the previous strongest limits set using the same effect by a factor of ~3.8.

We report the broadband spectro-temporal study of the poorly studied accreting neutron star (NS) low mass X-ray binary (LMXB) 4U 1724-30 using data from Soft X-ray Telescope (SXT) and Large Area X-ray Proportional Counters (LAXPC) instruments on board AstroSat. The dim persistent LMXB source was observed with AstroSat over 4 epochs in 2017, all of which corresponded to a low-luminosity non-thermal emission dominated (hard/island) emission state with modest spectral evolution. All the X-ray broadband spectra can be modelled by a combination of thermal emission from the NS boundary layer (BL) or NS surface and a non-thermal emission component possibly originating from the inverse Comptonization of the disc seed photons. We investigate the presence of frequency and energy-dependent variabilities to probe the origin of the disc/coronal fluctuations. We also report the detection of a Type-I X-ray burst displaying a photospheric radius expansion (PRE). During the burst, a hard X-ray shortage in the 30-80 keV energy band and the enhancement of the persistent emission reveal the burst feedback on the overall accretion process. Using the touch-down burst flux $\sim$ $4.25 \times 10^{-8}$ erg s$^{-1}$ cm$^{-2}$, the distance of the source is estimated as $\sim$ 8.4 kpc.

Interferometric observations of HI in galaxies played a pivotal role in studies of nearby galaxies. Compared to single-dish observation, it provides resolved distribution of gas in galaxies with unprecedented resolution. Several extensive HI surveys of nearby galaxies have been performed in the past; however, most of them consist of less than 100 galaxies due to individual efforts. On the other hand, present-day archives of the radio telescopes include data for at least several hundred galaxies. To utilize these data sets to their full potential, we construct a sample including all galaxies observed by the Giant Meter wave Radio Telescope (GMRT) in HI. This results in a total of 515 galaxies, the largest sample to date. We intend to analyze all the data uniformly and carry out different exciting science. As a pilot project, we analyze data from 11 galaxies and present the data products in this paper. We further investigate the neutral ISM in these galaxies and extract cold and warm phases using a multi-Gaussian decomposition method. This pilot project assures the quality of the data, which will enable us to perform critical science investigations using the full sample.

We present a sample of 32 stars of spectral types G and K, and luminosity classes I to V, with moderate activity levels, covering four orders of magnitude of surface gravity and a representative range of effective temperature. For each star we obtained high S/N TIGRE-HEROS spectra with a spectral resolving power of $R\approx20,000$ and have measured the Ca II K line-widths of interest, $W_0$ and $W_1$. The main physical parameters are determined by means of iSpec synthesis and Gaia EDR3 parallaxes. Mass estimates are based on matching to evolution models. Using this stellar sample, that is highly uniform in terms of spectral quality and assessment, we derive the best-fit relation between emission line width and gravity $g$, including a notable dependence on effective temperature $T_{\rm eff}$, of the form $W_1 \propto g^{-0.229} T_{\rm eff}^{+2.41}$. This result confirms the physical interpretation of the Wilson-Bappu effect as a line saturation and photon redistribution effect in the chromospheric Ca II column density, under the assumption of hydrostatic equilibrium at the bottom of the chromosphere. While the column density (and so $W_1$) increases towards lower gravities, the observed temperature dependence is then understood as a simple ionization effect -- in cooler stars, Ca II densities decrease in favor of Ca I.

Caffeination can open tired eyes and enhance focus. Over-caffeination, furthermore, can lead to errors but also to unexpected discoveries that might not have happened without 30 hours of sleep deprivation and 500mg of caffeine in our bodies. This paper presents exactly such a discovery. Upon much staring into our coffee cups, empty anew, the thought struck us: coffee in space. Caffeine may not be the only key. HL Tau, Taurus, bull... Taurine! We grinded some red bourbon for a new pour-over, and developed the new, coffee-groundsbreaking Large Astrocomical Taurine Tester Experiment (LATTE) in just 1/4 of a day. We felt bull-ish about our chances of making a great discovery! We installed LATTE, aimed it at the well-known young star HL Tau, and there it was: an abundance of taurine gas beautifully outlining a cup of cosmic flat white, with the ring structure of HL Tau turning out to be latte art performed by a skillful cosmic barista. The first Robusta discovery of coffee in space. Speaking of coffee, we hope you have a nice hot cup with you, and we encourage you to pun-tinue all the way to the end of this bean-grinding paper.

We present a novel method for the detection and removal of Radio Frequency Interference (RFI) from the raw digitized signal in the signal processing chain of a typical radio astronomy experiment. The main advantage of our method is that it does not require a training set. Instead, our method relies on the fact that the true signal of interest coming from astronomical sources is thermal and therefore described as a Gaussian random process, which cannot be compressed. We employ a variational encoder/decoder network to find the compressible information in the datastream that can explain the most variance with the fewest degrees of freedom. We demonstrate it on a set of toy problems and stored ringbuffers from the Baryon Mapping eXperiment (BMX) prototype. We discuss advantages and limitations of this method and possible implementation in the front-end of future radio experiments.

More often than not a lunch time conversation will veer off into bizarre and uncharted territories. In rare instances these frontiers of conversation can lead to deep insights about the Universe we inhabit. This paper details the fruits of one such conversation. In this paper we will answer the question: How many cows do you need to form a planetoid entirely comprised of cows, which will support a methane atmoosphere produced by the planetary herd? We will not only present the necessary assumptions and theory underpinning the cow-culations, but also present a thorough (and rather robust) discussion of the viability of, and implications for accomplishing, such a feat.

The ice shell of Enceladus exhibits strong asymmetry between its hemispheres, with all known geysers concentrated over the south pole, even though its orbital configuration is almost perfectly symmetric. By exploring ocean circulation across a range of ocean salinities and core/shell heating partitions, we study the role of ice-ocean interaction in hemispheric symmetry breaking. We find that: (i) asymmetry is enhanced by cross-equatorial ocean heat transport when the ice shell is the major heat source and vice versa, (ii) the magnitude of ocean heat transport is comparable to the global heat production, significantly affecting the ice shell evolution and equilibrium state and (iii) more than one equilibrium state can exist due to a positive feedback between melting and ocean circulation.

The ice shell and subsurface ocean on icy worlds are strongly coupled together -- heat and salinity flux from the ice shell induced by the ice thickness gradient drives circulation in the ocean, and in turn, the heat transport by ocean circulation shapes the ice shell. Since measurements in the near future are likely to remain constrained to above the ice shell, understanding this ocean-ice interaction is crucial. Using an ocean box model and a series of experiments simulating the 2D ocean circulation, we find that large icy moons with strong gravity tend to have stronger ocean heat transport under the same ice-shell topography. As a result, the equilibrium ice shell geometry is expected to be flatter on moons with larger size, and vice versa. This finding is broadly consistent with the observed ice shell geometry for Enceladus and Europa.

Steadily accreting white dwarfs (WDs) are efficient sources of ionization and thus, are able to create extended ionized nebulae in their vicinity. These nebulae represent ideal tools for the detection of accreting WDs, given that in most cases the source itself is faint. In this work, we combine radiation transfer simulations with known H and He accreting WD models, providing for the first time the ionization state and the emission line spectra of the formed nebulae as a function of the WD mass, the accretion rate and the chemical composition of the accreted material. We find that the nebular optical line fluxes and radial extent vary strongly with the WD's accretion properties, peaking in systems with WD masses of 0.8 - 1.2 $\rm~M_{\odot}$. Projecting our results on the 'BPT' diagnostic diagrams, we show that accreting WDs nebulae possess characteristics distinct from those of H II-like regions, while they share similar line ratios with the galactic low-ionization emission-line regions. Finally, we compare our results to the relevant constraints imposed by the lack of ionized nebulae in the vicinity of supersoft X-ray sources (SSSs) and Type Ia supernova remnants - sources which are related to steadily accreting WDs. The large discrepancies uncovered by our comparison rule out any steadily accreting WD as a potential progenitor of the studied remnants and additionally require the ambient medium around the SSSs to be less dense than 0.2 $\rm~cm^{-3}$. We discuss possible alternatives that could bridge the incompatibility between the theoretical expectations and the relevant observations.

We use $\Lambda$CDM cosmological hydrodynamical simulations to explore the kinematics of gaseous discs in late-type dwarf galaxies. We create high-resolution 21-cm 'observations' of simulated dwarfs produced in two variations of the EAGLE galaxy formation model: one where supernova-driven gas flows redistribute dark matter through gravitational interaction and form constant-density central 'cores', and another where the central 'cusps' survive intact. We 'observe' each galaxy along multiple sight lines and derive a rotation curve for each observation using a conventional tilted-ring approach to model the gas kinematics. We find that the modelling process introduces systematic discrepancies between the recovered rotation curve and the actual circular velocity curve driven primarily by (i) non-circular gas orbits within the discs; (ii) the finite thickness of gaseous discs, which leads to overlap of different radii in projection; and (iii) departures from dynamical equilibrium. Dwarfs with dark matter cusps often appear to have a core, whilst the inverse error is less common. These effects naturally reproduce an observed trend which other models struggle to explain: late-type dwarfs with more steeply-rising rotation curves appear to be dark matter-dominated in the inner regions, whereas the opposite seems to hold in galaxies with core-like rotation curves. We conclude that if similar effects affect the rotation curves of observed dwarfs, a late-type dwarf population in which all galaxies have sizeable dark matter cores is most likely incompatible with current measurements.

Some recent literature has claimed there to be an evolution in galaxies' dust temperatures towards warmer (or colder) spectral energy distributions (SEDs) between low and high redshift. These conclusions are driven by both theoretical models and empirical measurement. Such claims sometimes contradict one another and are prone to biases in samples or SED fitting techniques. What has made direct comparisons difficult is that there is no uniform approach to fitting galaxies' infrared/millimeter SEDs. Here we aim to standardize the measurement of galaxies' dust temperatures with a python-based SED fitting procedure, MCIRSED. We draw on reference datasets observed by IRAS, Herschel, and Scuba-2 to test for redshift evolution out to $z\sim2$. We anchor our work to the L$_{IR}$-$\lambda_{peak}$ plane, where there is an empirically observed anti-correlation between IR luminosity and rest-frame peak wavelength (an observational proxy for luminosity-weighted dust temperature) such that $\lambda_{peak}=\lambda_{t}({L_{IR}}/{L_{t}})^{\eta}$ where $\eta=-0.09\pm0.01$, L$_{t}=10^{12}$ L$_{\odot}$, and $\lambda_{t}=92\pm2\mu$m. We find no evidence for redshift evolution of galaxies' temperatures, or $\lambda_{peak}$, at fixed L$_{IR}$ from $0<z<2$ with >99.99% confidence. Our finding does not preclude evolution in dust temperatures at fixed stellar mass, which is expected from a non-evolving L$_{IR}$-$\lambda_{peak}$ relation and a strongly evolving SFR-M$_\star$ relation. The breadth of dust temperatures at a given L$_{IR}$ is likely driven by variation in galaxies' dust geometries and sizes and does not evolve. Testing for L$_{IR}$-$\lambda_{peak}$ evolution toward higher redshift ($z\sim5-6$) requires better sampling of galaxies' dust SEDs near their peaks (observed $\sim$200-600$\mu$m) with $<$1 mJy sensitivity. This poses a significant challenge to current instrumentation.

Current estimates of the normalized accretion rates of quasars (L/L_Edd), rely on measuring the velocity widths of broad optical-UV emission lines (e.g., H$\beta$ and Mg II $\lambda2800$). However, such lines tend to be weak or inaccessible in the most distant quasars, leading to increasing uncertainty in L/L_Edd estimates at $z > 6$. Utilizing a carefully selected sample of 53 radio-quiet quasars that have H$\beta$ and C IV $\lambda1549$ spectroscopy as well as {\sl Chandra} coverage, we searched for a robust accretion-rate indicator for quasars, particularly at the highest-accessible redshifts ($z \sim 6-7$). Our analysis explored relationships between the H$\beta$-based L/L_Edd, the equivalent width (EW) of C IV, and the optical-to-X-ray spectral slope (a_ox). Our results show that EW(C IV) is the strongest indicator of the H$\beta$-based L/L_Edd parameter, consistent with previous studies, although significant scatter persists particularly for sources with weak C IV lines. We do not find evidence for the a_ox parameter improving this relation, and we do not find a significant correlation between a_ox and H$\beta$-based L/L_Edd. This absence of an improved relationship may reveal a limitation in our sample. X-ray observations of additional luminous sources, found at $z \gtrsim 1$, may allow us to mitigate the biases inherent in our archival sample and test whether X-ray data could improve L/L_Edd estimates. Furthermore, deeper X-ray observations of our sources may provide accurate measurements of the hard-X-ray power-law photon index ($\Gamma$), which is considered an unbiased L/L_Edd indicator. Correlations between EW(C IV) and a_ox with $\Gamma$-based L/L_Edd may yield a more robust prediction of a quasar normalized accretion rate.

High mass stars play an important role in the Interstellar Medium, but much remains to be known about their formation. Class I methanol masers may be unique tracers of an early stage of high mass star formation, and a better understanding of such masers will allow them to be used as more effective probes of the high mass star forming process. We present an investigation of the long-term variability of Class I methanol masers at 44 GHz toward the high mass star forming region DR21(OH). We compare observations taken in 2017 to 2012 and also to 2001 data from the literature. A total of 57 maser spots were found in the 2017 data, with center velocities ranging between -8.65 km/s to +2.56 km/s. The masers are arranged in a western and an eastern lobe with two arcs in each lobe that look like bowshocks, consistent with previous observations. The general trend is an increase in intensity from 2001 to 2012, and a decrease from 2012 to 2017. Variability appears to be more prevalent in the inner arc of the western lobe than in the outer arc. We speculate that this may be a consequence of episodic accretion, in which a later accretion event has resulted in ejection of material whose shock reached the inner arc at some point in time after 2001. We conclude that class I methanol masers are variable on long timescales (of the order of 5-10 years).

Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observed properties. Binary population synthesis with full stellar-structure and evolution simulations are computationally expensive requiring a large number of mass-transfer sequences. The recently developed binary population synthesis code POSYDON incorporates grids of MESA binary star simulations which are then interpolated to model large-scale populations of massive binaries. The traditional method of computing a high-density rectilinear grid of simulations is not scalable for higher-dimension grids, accounting for a range of metallicities, rotation, and eccentricity. We present a new active learning algorithm, psy-cris, which uses machine learning in the data-gathering process to adaptively and iteratively select targeted simulations to run, resulting in a custom, high-performance training set. We test psy-cris on a toy problem and find the resulting training sets require fewer simulations for accurate classification and regression than either regular or randomly sampled grids. We further apply psy-cris to the target problem of building a dynamic grid of MESA simulations, and we demonstrate that, even without fine tuning, a simulation set of only $\sim 1/4$ the size of a rectilinear grid is sufficient to achieve the same classification accuracy. We anticipate further gains when algorithmic parameters are optimized for the targeted application. We find that optimizing for classification only may lead to performance losses in regression, and vice versa. Lowering the computational cost of producing grids will enable future versions of POSYDON to cover more input parameters while preserving interpolation accuracies.

Radio wave scattering can cause severe reductions in detection sensitivity for surveys of Galactic and extragalactic fast ($\sim$ms duration) transients. While Galactic sources like pulsars are subject to scattering in the Milky Way interstellar medium (ISM), extragalactic fast radio bursts (FRBs) can also experience scattering in their host galaxies and other galaxies intervening their lines-of-sight. We assess Galactic and extragalactic scattering horizons for fast radio transients using a combination of NE2001 to model the dispersion measure (DM) and scattering time ($\tau$) contributed by the Milky Way, and independently constructed electron density models for other galaxies' ISMs and halos that account for different galaxy morphologies, masses, densities, and strengths of turbulence. For FRB source redshifts $z_{\rm s} \lesssim 1$, an all-sky, isotropic FRB population has values of $\tau$ ranging between $\sim 1\ \mu$s and $\sim 2$ ms at 1 GHz (observer frame) that are dominated by host galaxies. For a hypothetical, high-redshift ($z_{\rm s}\sim5$) FRB population, $\tau$ ranges from $\sim 0.01 - 100$s of ms at 1 GHz, and is largely dominated by intervening galaxies. About $20\%$ of these high-redshift FRBs are predicted to have $\tau > 5$ ms at 1 GHz (observer frame), and $\gtrsim 40\%$ of FRBs between $z_{\rm s} \sim 0.5 - 5$ are predicted to have $\tau \gtrsim 1$ ms for $\nu\leq 800$ MHz. The percentage of FRBs selected against from scattering may be substantially larger because our scattering predictions are conservative compared to localized FRBs, and if circumgalactic turbulence causes density fluctuations larger than those observed from nearby halos.

Mid-Infrared (MIR) images of the Galactic center show extended gas and dust features along with bright IRS sources. Some of these dust features are a part of ionized clumpy streamers orbiting Sgr~A*, known as the mini-spiral. We present their proper motions over 12 year time period and report their flux densities in $N$-band filters {and derive their spectral indices}. The observations were carried out by VISIR at ESO VLT. High-pass filtering led to the detection of several resolved filaments and clumps along the mini-spiral. Each source was fit by a 2-D Gaussian profile to determine the offsets and aperture sizes. We perform aperture photometry to extract fluxes in two different bands. We present the proper motions of the largest consistent set of resolved and reliably determined sources. In addition to stellar orbital motions, we identify a stream-like motion of extended clumps along the mini-spiral. We also detect MIR counterparts of the radio tail components of the IRS7 source. They show a clear kinematical deviation with respect to the star. They likely represent Kelvin-Helmholtz instabilities formed downstream in the shocked stellar wind. We also analyze the shape and the orientation of the extended late-type IRS3 star that is consistent with the ALMA sub-mm detection of the source. Its puffed-up envelope with the radius of $\sim 2\times 10^6\,R_{\odot}$ could be the result of the red-giant collision with a nuclear jet, which was followed by the tidal prolongation along the orbit.

The high-magnification microlensing event KMT-2021-BLG-1077 exhibits a subtle and complex anomaly pattern in the region around the peak. We analyze the lensing light curve of the event with the aim of revealing the nature of the anomaly. We test various models in combination with several interpretations. We find that the anomaly cannot be explained by the usual three-body (2L1S and 1L2S) models. The 2L2S model improves the fit compared to the three-body models, but it still leaves noticeable residuals. On the other hand, the 3L1S interpretation yields a model explaining all the major anomalous features in the lensing light curve. According to the 3L1S interpretation, the estimated mass ratios of the lens companions to the primary are $\sim 1.56 \times 10^{-3}$ and $\sim 1.75 \times 10^{-3}$, which correspond to $\sim 1.6$ and $\sim 1.8$ times the Jupiter/Sun mass ratio, respectively, and therefore the lens is a multiplanetary system containing two giant planets. With the constraints of the event time-scale and angular Einstein radius, it is found that the host of the lens system is a low-mass star of mid-to-late M spectral type with a mass of $M_{\rm h} = 0.14^{+0.19}_{-0.07}~M_\odot$, and it hosts two gas giant planets with masses of $M_{\rm p_1}=0.22^{+0.31}_{-0.12}~M_{\rm J}$ and $M_{\rm p_2}=0.25^{+0.35}_{-0.13}~M_{\rm J}$. The planets lie beyond the snow line of the host with projected separations of $a_{\perp, {\rm p}_1}=1.26^{+1.41}_{-1.08}~{\rm AU}$ and $a_{\perp, {\rm p}_2}=0.93^{+1.05}_{-0.80}~{\rm AU}$. The planetary system resides in the Galactic bulge at a distance of $D_{\rm L}=8.24^{+1.02}_{-1.16}~{\rm kpc}$. The lens of the event is the fifth confirmed multiplanetary system detected by microlensing following OGLE-2006-BLG-109L, OGLE-2012-BLG-0026L, OGLE-2018-BLG-1011L, and OGLE-2019-BLG-0468L.

The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emission span from 50 to 1000\,au. This range could be explained by viscous evolution models with different $\alpha$ values (mostly of $10^{-4}-10^{-3}$) and/or a spread of initial conditions. The CO sizes for most disks are also consistent with MHD wind models that directly remove disk angular momentum, but very large initial disk sizes would be required to account for the very extended CO disks in the sample. As no CO size evolution is observed across stellar ages of 0.5--20\,Myr in this sample, determining the dominant mechanism of disk evolution will require a more complete sample for both younger and more evolved systems. We find that the CO emission is universally more extended than the continuum emission by an average factor of $2.9\pm1.2$. The ratio of the CO to continuum sizes does not show any trend with stellar mass, mm continuum luminosity, or the properties of substructures. The GO Tau disk has the most extended CO emission in this sample, with an extreme CO to continuum size ratio of 7.6. Seven additional disks in the sample show high size ratios ($\gtrsim4$) that we interpret as clear signs of substantial radial drift.

Pinpointing the neutrino sources is crucial to unveil the mystery of high-energy cosmic rays. The search for neutrino-source candidates from coincident neutrino-photon signatures and electromagnetic objects with peculiar flaring behaviors have the potential to increase our chances of finding neutrino emitters. In this paper, we first study the temporal correlations of astrophysical flares with neutrinos, considering a few hundreds of multi-frequency sources from ALMA, WISE, Swift, and Fermi in the containment regions of IceCube high-energy alerts. Furthermore, the spatial correlations between blazars and neutrinos are investigated using the subset of 10-year IceCube track-like neutrinos with around 250 thousand events. In the second test, we account for 2700 blazars with different types of flaring phases in addition to sole position. No significant neutrino emissions were found from our analyses. Our results indicate an interesting trend showing the infrared flaring stages of WISE blazars might be correlated with arrival times of the neutrino alerts. Possible overflow of neutrinos associated with two of our blazar sub-samples are also illustrated. One is characterized by a significant flaring lag in infrared with respect to gamma-rays, like seen for TXS0506+056, and the other is characterized by highly simultaneous infrared and gamma-ray flares. These phenomena suggest the need to improve current multi-frequency light-curve catalogs to pair with the advent of more sensitive neutrino observatories.

Presently, Pluto and Charon have been in tidal locking and the ratio of Pluto and Charon's spin velocity to the orbital mean motion is $1$. We test the effect of tidal evolution at Pluto-Charon and register the possibilities of the initial states with $\Delta t$ model and $Q$ model. In $\Delta t$ model, the tidal evolution equation will have an unnecessary overshoot, so we find a more exact tidal evolution equation for the $Q$ model. Besides we consider the influence of inclination and discuss the model with a high semi-major axis and high eccentricity.

Low surface brightness (LSB) galaxies are galaxies with central surface brightness fainter than the night sky. Due to the faint nature of LSB galaxies and the comparable sky background, it is difficult to search LSB galaxies automatically and efficiently from large sky survey. In this study, we established the Low Surface Brightness Galaxies Auto Detect model (LSBG-AD), which is a data-driven model for end-to-end detection of LSB galaxies from Sloan Digital Sky Survey (SDSS) images. Object detection techniques based on deep learning are applied to the SDSS field images to identify LSB galaxies and estimate their coordinates at the same time. Applying LSBG-AD to 1120 SDSS images, we detected 1197 LSB galaxy candidates, of which 1081 samples are already known and 116 samples are newly found candidates. The B-band central surface brightness of the candidates searched by the model ranges from 22 mag arcsec $^ {- 2} $ to 24 mag arcsec $^ {- 2} $, quite consistent with the surface brightness distribution of the standard sample. 96.46\% of LSB galaxy candidates have an axis ratio ($b/a$) greater than 0.3, and 92.04\% of them have $fracDev\_r$\textless 0.4, which is also consistent with the standard sample. The results show that the LSBG-AD model learns the features of LSB galaxies of the training samples well, and can be used to search LSB galaxies without using photometric parameters. Next, this method will be used to develop efficient algorithms to detect LSB galaxies from massive images of the next generation observatories.

To study the demographics of interstellar ices in the Central Molecular Zone (CMZ) of the Milky Way, we obtain near-infrared spectra of $109$ red point sources using NASA IRTF/SpeX at Maunakea. We select the sample from near- and mid-infrared photometry, including $12$ objects in the previous paper of this series, to ensure that these sources trace a large amount of absorption through clouds in each line of sight. We find that most of the sample ($100$ objects) show CO band-head absorption at $2.3\ \mu$m, tagging them as red (super-) giants. Despite the photospheric signature, however, a fraction of the sample with $L$-band spectra ($9/82=0.11$) exhibit large H$_2$O ice column densities ($N > 2\times10^{18}\ {\rm cm}^{-2}$), and six of them also reveal CH$_3$OH ice absorption. As one of such objects is identified as a young stellar object (YSO) in our previous work, these ice-rich sight lines are likely associated with background stars in projection to an extended envelope of a YSO or a dense cloud core. The low frequency of such objects in the early stage of stellar evolution implies a low star-formation rate ($<0.02\ M_\odot$ yr$^{-1}$), reinforcing the previous claim on the suppressed star-formation activity in the CMZ. Our data also indicate that the strong "shoulder" CO$_2$ ice absorption at $15.4\ \mu$m observed in YSO candidates in the previous paper arises from CH$_3$OH-rich ice grains having a large CO$_2$ concentration [$N {\rm (CO_2)} / N {\rm (CH_3OH)} \approx 1/3$].

Most types of massive stars display X-ray emission that is affected by the properties of their stellar winds. Single non-magnetic OB stars have an X-ray luminosity that scales with their bolometric luminosity and their emission is thought to arise from a distribution of wind-embedded shocks. The lack of significant short-term stochastic variability indicates that the winds consist of a large number of independent fragments. Detailed variability studies unveiled a connection between the photosphere and the wind: well-studied O-type stars exhibit a ~ 10% modulation of their emission on timescales consistent with the rotation period, and a few early B-type pulsators display ~ 10% modulations of their X-ray flux with the pulsation period. Unlike OB stars, their evolved descendants (WR and LBV stars) lack a well-defined relation between their X-ray and bolometric luminosities, and several subcategories of objects remain undetected. These properties most likely stem from the combined effects of wind optical depth and wind velocity. Magnetic OB stars display an enhanced X-ray emission frequently modulated by the rotation of the star. These properties are well explained by the magnetically confined wind shock model and an oblique magnetic rotator configuration. Some massive binaries display phase-dependent excess emission arising from the collision between the winds of the binary components. Yet, the majority of the massive binaries do not show such an emission, probably as a consequence of radiative cooling of the shock-heated plasma. Finally, a growing subset of the Be stars, the so-called gamma Cas stars, feature an unusually hard and strong thermal X-ray emission that varies over a wide range of timescales. Several scenarios have been proposed to explain these properties, but the origin of the phenomenon remains currently one of the major unsolved puzzles in stellar X-ray astrophysics.

Active Galactic Nuclei (AGN) are highly energetic astrophysical sources powered by accretion onto supermassive black holes in galaxies, which present unique observational signatures covering the full electromagnetic spectrum (and more) over about twenty orders of magnitude in frequency. We first review the main AGN properties and diversities and show that they can be explained by a small number of parameters. We then discuss the so-called Unification Models for non-jetted AGN, according to which these sources are believed to have the same nuclear engine and circumnuclear matter, with the same geometry for the obscuring structure. This simplified scenario, however, cannot explain all the observed complexities, such as the presence of multiple absorbers on different physical scales, including recent X-ray observations of circumnuclear matter. Finally, we touch upon AGN evolution in the X-ray and $\gamma$-ray bands.

The study of exoplanets and especially their atmospheres can reveal key insights on their evolution by identifying specific atmospheric species. For such atmospheric investigations, high-resolution transmission spectroscopy has shown great success, especially for Jupiter-type planets. Towards the atmospheric characterization of smaller planets, the super-Earth exoplanet 55 Cnc e is one of the most promising terrestrial exoplanets studied to date. Here, we present a high-resolution spectroscopic transit observation of this planet, acquired with the PEPSI instrument at the Large Binocular Telescope. Assuming the presence of Earth-like crust species on the surface of 55 Cnc e, from which a possible silicate-vapor atmosphere could have originated, we search in its transmission spectrum for absorption of various atomic and ionized species such as Fe , Fe+, Ca , Ca+, Mg and K , among others. Not finding absorption for any of the investigated species, we are able to set absorption limits with a median value of 1.9 x RP. In conclusion, we do not find evidence of a widely extended silicate envelope on this super-Earth reaching several planetary radii.

In the binary-driven hypernova (BdHN) scenario, long gamma-ray bursts (GRBs) originate in a cataclysmic event that occurs in a binary system composed of a carbon-oxygen (CO) star and a neutron star (NS) companion in close orbit. The collapse of the CO star generates at its center a newborn NS ($\nu$NS), and a supernova (SN) explosion. Matter from the ejecta is accreted both onto the $\nu$NS because of fallback and onto the NS companion, leading to the collapse of the latter into a black hole (BH). Each of the ingredients of the above system leads to observable emission episodes in a GRB. In particular, the $\nu$NS is expected to show up (hereafter $\nu$NS-rise) in the early GRB emission, nearly contemporary or superimposed to the ultrarelativistic prompt emission (UPE) phase, but with a different spectral signature. Following the $\nu$NS-rise, the $\nu$NS powers the afterglow emission by injecting energy into the expanding ejecta leading to synchrotron radiation. We here show that the $\nu$NS-rise and the subsequent afterglow emission in both systems, GRB 180720B and GRB 190114C, are powered by the release of rotational energy of a Maclaurin spheroid, starting from the bifurcation point to the Jacobi ellipsoid sequence. This implies that the $\nu$NS evolves from a triaxial Jacobi configuration, prior to the $\nu$NS-rise, into the axially symmetric Maclaurin configuration observed in the GRB. The triaxial $\nu$NS configuration is short-lived (less than a second) due to a copious emission of gravitational waves, before the GRB emission, and it could be in principle detected for sources located at distances closer than $100$ Mpc. This appears to be the sole process of emission of gravitational waves in long GRBs.

Many transient and variable sources detected at multiple wavelengths are also observed to vary at radio frequencies. However, these samples are typically biased towards sources that are initially detected in wide-field optical, X-ray or gamma-ray surveys. Many sources that are insufficiently bright at higher frequencies are therefore missed, leading to potential gaps in our knowledge of these sources and missing populations that are not detectable in optical, X-rays or gamma-rays. Taking advantage of new state-of-the-art radio facilities that provide high quality wide-field images with fast survey speeds, we can now conduct unbiased surveys for transient and variable sources at radio frequencies. In this paper, we present an unbiased survey using observations obtained by MeerKAT, a mid-frequency ($\sim$1.4 GHz) radio array in South Africa's Karoo Desert. The observations used were obtained as part of a weekly monitoring campaign for X-ray binaries (XRBs) and we focus on the field of MAXI J1820+070. We develop methods to optimally filter transient and variable candidates that can be directly applied to other datasets. In addition to MAXI J1820+070, we identify four likely active galactic nuclei, one source that could be a Galactic source (pulsar or quiescent X-ray binary) or an AGN, and one variable pulsar. No transient sources, defined as being undetected in deep images, were identified leading to a transient surface density of $<3.7\times10^{-2}$ deg$^{-2}$ at a sensitivity of 1 mJy on timescales of one week at 1.4 GHz.

We analyse five early science datasets from the APERture Tile in Focus (Apertif) phased array feed system to verify the polarisation capabilities of Apertif in view of future larger data releases. We aim to characterise the source population of the polarised sky in the L-Band using polarised source information in combination with IR and optical data. We use automatic routines to generate full field-of-view Q- and U-cubes and perform RM-Synthesis, source finding, and cross-matching with published radio, optical, and IR data to generate polarised source catalogues. SED-fitting routines were used to determine photometric redshifts, star-formation rates, and galaxy masses. IR colour information was used to classify sources as AGN or star-forming-dominated and early- or late-type. We surveyed an area of 56deg$^2$ and detected 1357 polarised source components in 1170 sources. The fraction of polarised sources is 10.57% with a median fractional polarisation of 4.70$\pm$0.14%. We confirmed the reliability of the Apertif measurements by comparing them with polarised cross-identified NVSS sources. Average RMs of the individual fields lie within the error of the best Milky Way foreground measurements. All of our polarised sources were found to be dominated by AGN activity in the radio regime with most of them being radio-loud (79%) and of the FRII class (87%). The host galaxies of our polarised source sample are dominated by intermediate disc and star-forming disc galaxies. The contribution of star formation to the radio emission is on the order of a few percent for $\approx$10% of the polarised sources while for $\approx$90% it is completely dominated by the AGN. We do not see any change in fractional polarisation for different star-formation rates of the AGN host galaxies.

Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes in the 21-cm and 11-cm bands. Additionally, we aim to identify future possible lines of study using this long-term scintillation dataset. We present the long-term time series of $\nu_{\rm d}$ and $\tau_{\rm d}$ for 13 pulsars. Sanity-checks and comparisons indicate that the scintillation parameters of our work and previously published works are mostly consistent. For two pulsars, PSRs~J1857+0943 and J1939+2134, we were able to obtain measurements of the $\nu_{\rm d}$ at both bands, which allows us to derive the time series of frequency scaling indices with a mean and a standard deviation of 2.82$\pm$1.95 and 3.18$\pm$0.60, respectively. We found some interesting features which will be studied in more detail in subsequent papers in this series: (i) in the time series of PSR~J1939+2134, where the scintillation bandwidth sharply increases or decreases associated with a sharp change of dispersion measure; (ii) PSR~J0613$-$0200 and PSR~J0636+5126 show a strong annual variation in the time series of the $\tau_{\rm d}$; (iii) PSR~J1939+2134 shows a weak anti-correlation between scintillation timescale and dispersion in WSRT data.

The novel coronavirus, dubbed COVID-19, upended our lives may be in irreversible ways during its initial spread throughout the world in March 2020. It forced us all, willingly or unwillingly, to keep social distance from each other to slow down the spread of COVID-19. As scientists, we started speculating what kind of separation is between the constitutes of different objects in the Universe. In this work, we study the "social" distance between elements inside various objects, no matter their size, mass, and nature. We consider things ranging from diamond, baseball to Saturn, asteroid belt or M87 Black Hole, to name a few. We show our results in the form of a fascinating mass/"social" distance plot, where a cool cartoon figure represents each object.

We use a sample of 536 \hii\ regions located in nearby spirals, with an homogeneous determination of their $T_e$-based abundances, to obtain new empirical calibrations of the N2O2, N2S2, O3N2, and N2 strong-line indices to estimate the nitrogen-to-oxygen abundance ratio when auroral lines are not detected. All indices are strongly correlated with the $T_e$-based $\log$(N/O) for our \hii\ region sample, even more strongly than with $12+\log$(O/H). N2O2 is the most strongly correlated index, and the best fit to the $\log$(N/O)-N2O2 relation is obtained with a second-order polynomial. The derived relation has a low dispersion ({\em rms}$<$0.09~dex), being valid in the range $-1.74 < $ N2O2 $< 0.62$ (or $-1.81 < $ $\log$(N/O) $< -0.13$). We have compared our calibration with previous ones and have discussed the differences between them in terms of the nature of the objects used as calibrators.

We present, for the first time, the relationship between local stellar mass surface density, $\mathrm{\Sigma_{*}}$, and N/O derived from SDSS-IV MaNGA data, using a sample of $792765$ high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modelling and partial correlation analysis, we find that $\mathrm{\Sigma_{*}}$ alone is insufficient to predict the N/O in MaNGA spaxels, and that there is an additional dependence on the local star formation rate surface density, $\mathrm{\Sigma_{SFR}}$. This effect is a factor of $3$ stronger than the dependence of 12+log(O/H) on $\mathrm{\Sigma_{SFR}}$. Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed $\Sigma_{*}$ as well as the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when $\mathrm{\Sigma_{*}}$ and $\mathrm{\Sigma_{SFR}}$ are controlled for. We show that $\sim50\%$ of the variance of N/O is explained by the total stellar mass and size. Thus, the evolution of nitrogen in galaxies is set by more than just local effects and does not simply track the build up of oxygen in galaxies. The precise form of the N/O-O/H relation is therefore sensitive to the sample of galaxies from which it is derived. This result casts doubt on the universal applicability of nitrogen-based strong-line metallicity indicators derived in the local universe.

Decayless kink oscillations of solar coronal loops (or decayless oscillations for short) have attracted great attention since their discovery. Coronal bright points (CBPs) are mini-active regions and consist of loops with a small size. However, decayless oscillations in CBPs have not been widely reported. In this study, we identified this kind of oscillations in some CBPs using 171 \AA\, images taken by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). After using the motion magnification algorithm to increase oscillation amplitudes, we made time-distance maps to identify the oscillatory signals. We also estimated the loop lengths and velocity amplitudes. We analysed 23 CBPs, and found 31 oscillation events in 16 of them. The oscillation periods range from 1 to 8 minutes (on average about 5 minutes), and the displacement amplitudes have an average value of 0.07 Mm. The average loop length and velocity amplitude are 23 Mm and 1.57 \kms, respectively. Relationships between different oscillation paraments are also examined. Additionally, we performed a simple forward model to illustrate how these sub-pixel oscillation amplitudes (less than 0.4 Mm) could be detected. Results of the model confirm the reliability of our data processing methods. Our study shows for the first time that decayless oscillations are common in small-scale loops of CBPs. These oscillations allow for seismological diagnostics of the Alfv\'{e}n speed and magnetic field strength in the corona.

Photometric redshift (photo-z) is a fundamental parameter for multi-wavelength photometric surveys, while galaxy clusters are important cosmological probers and ideal objects for exploring the dense environmental impact on galaxy evolution. We extend our previous work on estimating photo-z and detecting galaxy clusters to the latest data releases of the Dark Energy Spectroscopic Instrument (DESI) imaging surveys, Dark Energy Survey (DES), and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) imaging surveys and make corresponding catalogs publicly available for more extensive scientific applications. The photo-z catalogs include accurate measurements of photo-z and stellar mass for about 320, 293, and 134 million galaxies with $r<23$, $i<24$, and $i<25$ in DESI DR9, DES DR2, and HSC-SSP PDR3 data, respectively. The photo-z accuracy is about 0.017, 0.024, and 0.029 and the general redshift coverage is $z<1$, $z<1.2$, and $z<1.6$, respectively for those three surveys. The uncertainties of the logarithmic stellar mass that is inferred from stellar population synthesis fitting is about 0.2 dex. With the above photo-z catalogs, galaxy clusters are detected using a fast cluster-finding algorithm. A total of 532,810, 86,963, and 36,566 galaxy clusters with the number of members larger than 10 are discovered for DESI, DES, and HSC-SSP, respectively. Their photo-z accuracy is at the level of 0.01. The total mass of our clusters are also estimated by using the calibration relations between the optical richness and the mass measurement from X-ray and radio observations. The photo-z and cluster catalogs are available at ScienceDB (https://www.scidb.cn/s/2AjaEb) and PaperData Repository (https://nadc.china-vo.org/article/20200722160959?id=101089).

21cm Intensity Mapping (IM) has been proposed about 15 years ago as a cost effective method to carry out cosmological surveys and to map the 3D distribution of matter in the universe, over a large range of post EoR redshifts, from z=0 to z=6. Since then a number of pathfinder instruments have been built, such as CHIME or Tianlai. Several other ones will be commissioned in the next few years (HIRAX, CHORD, BINGO), while even larger arrays, with several thousand antennae are being considered for the next generation experiments. We will briefly review the 21cm cosmology of the Epoch of Reionisation (EoR), and we will then focus on IM for late time cosmology. After presenting some of the promises of this technique to constrain the cosmological model, dark energy and inflation, we will review some of the instrumental and scientific challenges of IM surveys. The second part of the paper presents an overview of the ongoing and future experiments, as well as recent results by GBT, CHIME and Tianlai.

Using smoothed particle hydrodynamics we model giant impacts of Super-Earth mass rocky planets between an atmosphere-less projectile and an atmosphere-rich target. In this work we present results from head-on to grazing collisions. The results of the simulations fall into two broad categories: 1) one main post-collision remnant containing material from target and projectile; 2) two main post-collision remnants resulting from `erosive hit-and-run' collisions. All collisions removed at least some of the target atmosphere, in contrast to the idealised hit-and-run definition in which the target mass is unchanged. We find that the boundary between `hit-and-run' collisions and collisions that result in the projectile and target accreting/merging to be strongly correlated with the mutual escape velocity at the predicted point of closest approach. Our work shows that it is very unlikely for a single giant impact to remove all of the atmosphere. For all the atmosphere to be removed, head-on impacts require roughly the energy of catastrophic disruption (i.e. permanent ejection of half the total system mass) and result in significant erosion of the mantle. We show that higher impact angle collisions, which are more common, are less efficient at atmosphere removal than head-on collisions. Therefore, single collisions that remove all the atmosphere without substantially disrupting the planet are not expected during planet formation.

It has previously been suggested that the rate of exoplanet discovery, doubling roughly every 39 months, is indicative of a runaway increase in the number of exoplanets in our galaxy. In this paper, we posit that - due to the finite nature of space in the milky way - it will become increasingly likely that one of these exoplanets will be found within our solar system. We calculate the odds of this occurring pass 50\% on Friday 9th December 2146. We go on to suggest novel methods for influencing where this exoplanet may be discovered, note possible drawbacks of the discovery, and finally explore how the previously-hypothesized `exoplanet singularity' (both figurative and literal) could be averted.

The reconstruction of the differential emission measure (DEM) from observations of spectral line intensities provides information on the temperature distribution of the emission measure in the region observed. The inversion process is known to be highly unstable and it has been necessary to provide additional constraints, such as requiring that the DEM should be smooth. However, this is a non-physical constraint. The goal of this analysis is to make an empirical determination of the ability of a set of emission line intensities to constrain the reconstruction. Here, a simple model is used, by means of a Monte-Carlo-Markov-Chain process, to arrive at solutions that reproduces the observed intensities in a region of the quiet Sun and a solar active region. These solutions are compared by means of the reduced chi-squared. The conclusion from this analysis is that the observations are only capable of constraining model consisting of 4 temperature-emission measure pairs plus a determination of the standard deviation of the model from the observed line intensities. A more complex model with 5 temperature-emission measure pairs does not improve the fit and leads to parameters that are irrelevant. A more general conclusion is that the information content of a set of observed emission lines is limited with respect to the determination of the emission measure distribution.

We have carried out detailed time-resolved timing analyses of three cataclysmic variables (CVs) namely LS Cam, V902 Mon, and SWIFT J0746.3-1608, using the long-baseline, high-cadence optical photometric data from the Transiting Exoplanet Survey Satellite (TESS). Our analysis of LS Cam observations hints the presence of a superorbital period of $\sim$ 4.025$\pm$0.007 d along with negative and positive superhump periods of $\sim$ 3.30 h and 3.70 h, respectively. These results can be explained as an interaction of nodal and apsidal precession of the accretion disc with orbital motion. For the other two sources, V902 Mon and SWIFT J0746.3-1608, we have found evidence of a beat period of 2387.0$\pm$0.6 s and 2409.5$\pm$0.7 s, respectively, which were not found in earlier studies. Our results presented in this study indicate the change in the accretion mode during the entire observing period for both sources. For V902 Mon, an apparent orbital period derivative of (6.09 $\pm$ 0.60) $\times 10^{-10}$ was also found. Moreover, the second harmonic of orbital frequency dominates the power spectrum of SWIFT J0746.3-1608, suggestive of ellipsoidal modulation of the secondary star. Present analyses suggest that LS Cam could be a superhumping CV whereas V902 Mon and SWIFT J0746.3-1608 are likely to be variable disc-overflow accreting intermediate polars.

I present for your appraisal three independent cases of the manuscript referee process conducted by a venerable peer-reviewed scientific journal. Each case involves a little pig, who submitted for consideration a theoretical plan for a house to be constructed presently, in a faraway land. An anonymous big bad wolf was assigned by the journal to assess the merit of these manuscripts. The pigs proposed three distinct construction frameworks, which varied in physical and mathematical sophistication. The first little pig submitted a model of straw, based on the numerical method of toe-counting. His design included odd features, such as spilled millet and cloven-hoofprints on the window sill -- possibly a ploy to distract the wolf from the manuscript's facile mathematical foundation. The second little pig used a more advanced approach, employing Newton's classical laws of motion, to propose a house of sticks. This pig included in her manuscript copious citations by a specific wolf, possibly aiming to ensure acceptance by flattering the wolf whom she anticipated would be the referee. The third little pig described an ostentatious house of bricks based on an elaborate dynamical systems and stability analysis, possibly scheming to dazzle and impress. The big bad wolf did not appear moved by any of the pigs' tactics. His recommendations were, for straw: the minor revision of water-proofing; for sticks: the major revision of fire-proofing, given concerns surrounding climate change; for bricks: unequivocal rejection, accompanied by multiple derogatory comments regarding "high-and-mighty theorists." I describe each case in detail, and suggest that the wolf's reports might be driven as much by self interest as the manuscripts themselves -- namely, that at the time the wolf wrote the reviews, he was rather hungry. Finally, I examine morals learned, if any.

Despite all its well-known flaws and calls for its dismissal, the notorious $h$-index is still used in many instances when awarding grants, or promoting and hiring scientists. To address this, I set out to devise a better index, with the twofold aim of taking into account the authors' respective contributions and considerably reducing the pollution of the scientific literature. Finally, I present a strategy that is guaranteed to be best for all researchers.

The Zwicky Transient Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT2019dsg and likely TDE AT2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that -21 can contribute no more than 87% of the total, while transients brighter than -22 can contribute no more than 58% of the total, neglecting the effect of extinction. These are the the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, suggesting that most astrophysical neutrinos are not produced during such optical flares. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory.

We consider an anisotropic search for the stochastic gravitational-wave (GW) background by decomposing the gravitational-wave sky into its spherical harmonics components. Previous analyses have used the diffraction limit to define the highest-order spherical harmonics components used in this search. We investigate whether the angular resolution of this search is indeed diffraction-limited by testing our ability to detect and localize simulated GW signals. We show that while using low-order spherical harmonics modes is optimal for initially detecting GW sources, the detected sources can be better localized with higher-order spherical harmonics than expected based on the diffraction limit argument. Additionally, we discuss how the ability to recover simulated GW sources is affected by the number of detectors in the network, the frequency range over which the search is performed, and the method by which the covariance matrix of the GW skymap is regularized. While we primarily consider point-source signals in this study, we briefly apply our methodology to spatially-extended sources and discuss potential future modifications of our analysis for such signals.

We present ALMA observations at a spatial resolution of 0.2 arcsec (60 pc) of CO emission from the Taffy galaxies (UGC 12914/5). The observations are compared with narrow-band Pa$\alpha$, mid-IR, radio continuum and X-ray imaging, plus optical spectroscopy. The galaxies have undergone a recent head-on collision, creating a massive gaseous bridge which is known to be highly turbulent. The bridge contains a complex web of narrow molecular filaments and clumps. The majority of the filaments are devoid of star formation, and fall significantly below the Kennicutt-Schmidt relationship for normal galaxies, especially for the numerous regions undetected in Pa$\alpha$ emission. Within the loosely connected filaments and clumps of gas we find regions of high velocity dispersion which appear gravitationally unbound for a wide range of likely values of $X_{\rm CO}$. Like the "Firecracker" region in the Antennae system, they would require extremely high external dynamical or thermal pressure to stop them dissipating rapidly on short crossing timescales of 2-5~Myrs. We suggest that the clouds may be transient structures within a highly turbulent multi-phase medium which is strongly suppressing star formation. Despite the overall turbulence in the system, stars seem to have formed in compact hotspots within a kpc-sized extragalactic HII region, where the molecular gas has a lower velocity dispersion than elsewhere, and shows evidence for a collision with an ionized gas cloud. Like the shocked gas in the Stephan's Quintet group, the conditions in the Taffy bridge shows how difficult it is to form stars within a turbulent, multi-phase, gas.

Solar wind charge exchange (SWCX) is the primary contamination to soft X-ray emission lines from the Milky Way (MW) hot gas. We report a solar-cycle ($\approx 10$ yr) temporal variation of observed \ion{O}{7} and \ion{O}{8} emission line measurements in the {\it XMM-Newton} archive, which is tightly correlated with the solar cycle traced by the sunspot number (SSN). This temporal variation is expected to be associated with the heliospheric SWCX. Another observed correlation is that higher solar wind (SW) fluxes lead to higher O VII or O VIII fluxes, which is due to the magnetospheric SWCX. We construct an empirical model to reproduce the observed correlation between the line measurements and the solar activity (i.e., the SW flux and the SSN). With this model we discovered a lag of $0.91_{-0.22}^{+0.20}$ yr between the O VII flux and the SSN. This time lag is a combination of the SW transit time within the heliosphere, the lag of the neutral gas distribution responding to solar activity, and the intrinsic lag between the SSN and the launch of a high-energy SW (i.e., $\rm O^{7+}$ and $\rm O^{8+}$). MW O VII and O VIII fluxes have mean values of 5.4 L.U. and 1.7 L.U., which are reduced by $50\%$ and $30\%$, compared to studies where the SWCX contamination is not removed. This correction also changes the determination of the density distribution and the temperature profile of the MW hot gas.

The census of known exoplanets exhibits a variety of physical parameters, including densities that are measured to span the range from less dense than styrofoam to more dense than iron. These densities represent a large diversity of interior structures. Despite this staggering diversity, recent analyses have shown that the densities of planets that orbit a common star exhibit remarkable uniformity. A fascinating exception to this is the system HIP 41378 (also known as K2-93), which contains a super-puff planet, HIP 41378 f, as well as several planets with more typical bulk densities. The range of densities in this system begs the question of what physical processes are responsible for the disparate planetary structures in this system. In this paper, we consider how the densities of the planets in the HIP 41378 system would have changed over time as the host star evolved and the planets' atmospheres were subsequently affected by the evolving insolation level. We also present a range of allowable core masses for HIP 41378 f based on the measured planet parameters, and comment on the feasibility of the proposed existence of planetary rings around HIP 41378 f as an explanation for its current low density.

With the proliferation of online Zoom meetings as a means of doing science in the 2020s, astronomers have made new and unexpected Target of Opportunity (ToO) observations. Chief among these ToOs are observations of exop(lan)ets, or "exopets." Building on the work of Mayorga et al. (2021) - whose work characterized the rotational variations of "floofy" objects - we model exopets using methods similar to those used for exoplanetary transits. We present data collected for such exopet Zoom transits through a citizen science program in the month of February 2022. The dataset includes parameters like exopet color, floofiness, transit duration, and percentage of Zoom screen covered during the event. For some targets, we also present microlensing and direct imaging data. Using results from our modelling of 62 exopet observations as transits, microlensing, and direct imaging events, we discuss our inferences of exopet characteristics like their masses, sizes, orbits, colors, and floofiness.

Primordial non-Gaussianities of the scalar(tensor)-tensor-tensor type supporting a non-trivial squeezed component are known to induce anisotropies in the stochastic gravitational wave background. We derive the explicit form of such anisotropies by making use, for the first time in this context, of the in-in formalism for cosmological correlation functions. After illustrating the general method and using it for the minimal single-field slow-roll case, we apply it to multi-field models, providing both a tree-level and a one-loop example. First, we make contact with previous results on anisotropies due to the presence of an extra spin-2 field during inflation. Secondly, we calculate the 1-loop scalar-tensor-tensor three-point function in the context of so-called supersolid inflation. The corresponding gravitational wave anisotropy is induced atop a gravitational signal that may be sufficiently large for detection.

We present optical photometric and spectroscopic observations of two Cataclysmic Variables (CVs), namely RBS 0490 and SDSS J075939.79+191417.3. The optical variations of RBS 0490 have been found to occur at the period of 1.689$\pm$0.001 hr which appears to be a probable orbital period of the system. Present photometric observations of SDSS J075939.79+191417.3 confirm and refine the previously determined orbital period as 3.14240928$\pm$0.00000096 hr. The presence of long-duration eclipse features in the light curves of SDSS J075939.79+191417.3 indicates eclipses might be due to an accretion disc and bright spot. The orbital inclination of SDSS J075939.79+191417.3 is estimated to be $\sim$ 78 $^\circ$ using the eclipse morphology. The phased-light curve variations during the orbital cycle of RBS 0490 provide evidence of the emission from an independent second accretion region or a second fainter pole. Optical spectra of RBS 0490 and SDSS J075939.79+191417.3 show the presence of strong Balmer, weak He II ($\lambda$4686) emission lines, along with the detection of strong $H\beta$ emission lines with a large value of equivalent width. The characteristic features of RBS 0490 seem to favour low-field polars, while SDSS J075939.79+191417.3 appears to be similar to the non-magnetic systems.

We investigate using three-point statistics in constraining the galaxy-halo connection. We show that for some galaxy samples, the constraints on the halo occupation distribution parameters are dominated by the three-point function signal (over its two-point counterpart). We demonstrate this on mock catalogs corresponding to the Luminous Red Galaxies (LRGs), Emission-Line Galaxies (ELG), and quasars (QSOs) targeted by the Dark Energy Spectroscopic Instrument (DESI) Survey. The projected three-point function for triangle sides less up to 20$h^{-1}$ Mpc measured from a cubic Gpc of data can constrain the characteristic minimum mass of the LRGs with a precision of $0.46$ %. For comparison, similar constraints from the projected two-point function are $1.55$ %. The improvements for the ELGs and QSOs targets are more modest. In the case of the QSOs it is caused by the high shot-noise of the sample, and in the case of the ELGs, this is caused by the range of halo masses of the host halos. The most time-consuming part of our pipeline is the measurement of the three-point functions. We adopt a tabulation method, proposed in earlier works for the two-point function, to reduce significantly the required compute time for the three-point analysis.

Advanced LIGO and Virgo have reported ninety confident gravitational-wave (GW) observations from compact-binary coalescences from their three observation runs. In addition, numerous subthreshold gravitational-wave candidates have been identified. Binary neutron star (BNS) mergers can produce gravitational waves and short-gamma ray bursts, as confirmed by GW170817/GRB 170817A. There may be electromagnetic counterparts recorded in archival observations associated with subthreshold gravitational-wave candidates. The CHIME/FRB collaboration has reported the first large sample of fast radio bursts (FRBs), millisecond radio transients detected up to cosmological distances; a fraction of these may be associated with BNS mergers. This work searches for coincident gravitational waves and FRBs from BNS mergers using candidates from the 4th-Open Gravitational-wave Catalog (4-OGC) and the first CHIME/FRB catalog. We use a ranking statistic for GW/FRB association which combines the gravitational-wave detection statistic with the odds of temporal and spatial association. We analyze gravitational-wave candidates and non-repeating FRBs from 2019 April 1 to 2019 July 1, when both the Advanced LIGO/Virgo gravitational-wave detectors and the CHIME radio telescope were observing. The most significant coincident candidate has a false alarm rate of 0.29 per observation time, which is consistent with a null observation. The null results imply at most $\mathcal{O}(0.01)\%$ - $\mathcal{O}(1)\%$ of FRBs are produced from the BNS mergers.

Within the context of hot big-bang cosmology, a cosmic background of presently low energy neutrinos is predicted to exist in concert with the photons of the cosmic background radiation. The number density of the cosmological neutrinos is of the same order as that of the photons of the cosmic background radiation. That makes neutrinos the second most abundant particle species in the universe. In the early universe, when these neutrinos were highly relativistic, their effects in determining the ultimate structure and evolution of the universe were significant.

Here we propose a simple explanation why echoes from wormholes mimicking black holes may be so small that they cannot be observed. The essence of the effect is in the redistribution of the initial energy of gravitational wave among multiple universes, connected by a wormhole.

$\textit{The Bachelor}$ is a reality TV dating show in which a single bachelor selects his wife from a pool of approximately 30 female contestants over eight weeks of filming (American Broadcasting Company 2002). We collected the following data on all 422 contestants that participated in seasons 11 through 25: their Age, Hometown, Career, Race, Week they got their first 1-on-1 date, whether they got the first impression rose, and what "place" they ended up getting. We then trained three machine learning models to predict the ideal characteristics of a successful contestant on $\textit{The Bachelor}$. The three algorithms that we tested were: random forest classification, neural networks, and linear regression. We found consistency across all three models, although the neural network performed the best overall. Our models found that a woman has the highest probability of progressing far on $\textit{The Bachelor}$ if she is: 26 years old, white, from the Northwest, works as an dancer, received a 1-on-1 in week 6, and did not receive the First Impression Rose. Our methodology is broadly applicable to all romantic reality television, and our results will inform future $\textit{The Bachelor}$ production and contestant strategies. While our models were relatively successful, we still encountered high misclassification rates. This may be because: (1) Our training dataset had fewer than 400 points or (2) Our models were too simple to parameterize the complex romantic connections contestants forge over the course of a season.

Quantum field theory, which is generally used to describe the origin of large-scale gravitational perturbations during cosmic inflation, has been shown to omit an important physical effect in curved space-time, the nonlocal entanglement among quantized modes from their gravitational effect on causal structure. It is argued here that in a different model of quantum gravity that coherently preserves nonlocal directional and causal relationships, primordial perturbations originate instead from coherent quantum distortions of emergent inflationary horizons; and moreover, that causal constraints account for approximate symmetries of cosmic microwave background correlations measured at large angular separations, which are highly anomalous in the standard picture. Thus, symmetries already apparent in the large-angle CMB pattern may be unique signatures of the emergence of locality and causal structure from quantum gravity.

In this work I reason that in expanding space only those quantum modes contribute to the measured vacuum energy that do not transcend the observable volume. Since all quantised field modes have various observable consequences, when a gravitational horizon causally confines an observer to a finite volume quantised modes should be restricted to the observable patch to remain consistent with gravity. Within the observable patch of Friedmann-Lemaitre-Robertson-Walker (FLRW) space the vacuum expectation value of the energy-momentum tensor can be expressed as a sum over discrete field modes. Friedmann's first equation provides a straightforward ultraviolet cut-off allowing only a finite number of modes in the sum. The finite volume acts as an infrared regulator and the calculation of the vacuum energy density is tractable without regularisation and renormalisation. To test the validity of this idea I quantise a scalar field on an FLRW background and calculate its vacuum energy density in the vacuum dominated, conformal, holographic limit. In this limit I show that the quantum vacuum energy density scales with the square of the Hubble parameter, consistently with gravity. In this example quantum vacuum expands space while the horizon of the expanding space limits the energy density of the vacuum to the observed value.

We give conditions for an exoplanetary system to function as an ideal amusement park/vacation resort (with its separate parking lot, of course); in case of massive human interplanetary colonization. Our considerations stem from the fact that an amusement park needs a parking lot of roughly the same surface area, thus the best option for its construction would be a system with at least 2 planets close to each other for easy tourist transportation. We also discuss the likelihood of finding such a system out there to cut down on construction costs.

We study the properties of pasta structures and their influence on the neutron star observables employing the effective relativistic mean-field theory (E-RMF). The compressible liquid drop model is used to incorporate the finite size effects, considering the possibility of non-spherical structures in the inner crust. The unified equation of state are constructed for several E-RMF parameters to study various properties such as pasta mass and thickness in the neutron star's crust. The majority of the pasta properties are sensitive to the symmetry energy in the subsaturation density region. Using the results from Monte Carlo simulations, we estimate the shear modulus of the crust in context of quasiperiodic oscillations from soft gamma-ray repeaters and calculate the frequency of fundamental torsional oscillation mode in the inner crust. Global properties of the neutron star such as mass-radius profile, the moment of inertia, crustal mass, crustal thickness and fractional crustal moment of inertia are worked out. The results are consistent with various observational and theoretical constraints.

We investigate the well-known phenomenon of the beam-plasma instability in the gravitational sector, when a fast population of particles interacts with the massive scalar mode of an Horndeski theory of gravity, resulting into the linear growth of the latter amplitude. Following the approach used in the standard electromagnetic case, we start from the dielectric representation of the gravita- tional plasma, as introduced in a previous analysis of the Landau damping for the scalar Horndeski mode. Then, we set up the modified Vlasov-Einstein equation, using at first a Dirac delta func- tion to describe the fast beam distribution. This way, we provide an analytical expression for the dispersion relation and we demonstrate the existence of non-zero growth rate for the linear evolu- tion of the Horndeski scalar mode. A numerical investigation is then performed with a trapezoidal beam distribution function, which confirms the analytical results and allows to demonstrate how the growth rate decreases as the beam spread increases.

Though it may be a behavior that has been observed and documented for millennia, and despite the connection between it and the full moon, the astronomical community has afforded very little attention to lycanthropy. We hope to address this deficiency by using the population of known exoplanets as a natural experiment to better characterize what properties of the moon are necessary to trigger a transformation into a werewolf. We additionally investigate which exoplanets are most likely to have exomoons which may cause werewolves, with a particular focus on LHS 1140 b. We also propose a new mission called the Werewolves From Infrared Radiation and Spectral-typing Telescope, or WFIRST, in order to better characterize exoplanetary systems. This will allow us to explore the impact of stellar type on lycanthropy more than it has traditionally been considered. We believe this represents a major step forward in our understanding and recognition of the burgeoning field of exocryptozoology.

We outline a theoretical approach supporting strong phase transitions from normal nuclear matter to the deconfined quark-gluon plasma, in the equation of state (EOS) for compact star matter. Implications of this hypothesis are discussed for astrophysical applications. Special emphasis is devoted to potentially detectable signatures, which can be directly related with the occurrence of a sufficiently strong phase transition. Therefore, simulations of core-collapse supernovae and binary compact star mergers are considered, including the subsequent emission of gravitational waves and, in the case of supernova, in addition the neutrinos play the role of messengers.

We exploit a new theory of gravity proposed by Finzi, which gives stronger interaction at large scales, to study the thermodynamic description of galaxy clusters. We employ a statistical model to deduce various thermodynamics equations of state. In addition, we analyze the behavior of clustering parameter in comparison to its standard (Newtonian) counterpart. The general distribution function and its behavior with varying strength of clustering parameter are also studied. The possibility of phase transition is also investigated and observed that a phase transition is possible though hierarchically.

We study the effect of nonfactorizable background geometry on the thermodynamics of the clustering of galaxies. A canonical partition function is derived for the gravitating system of galaxies treated as point particles contained in cells of appropriate dimensions. Various thermodynamic equations of state, like Helmholtz free energy and entropy, among others, are also obtained. We also estimate the effect of the corrected Newton law on the distribution function of galaxies. Remarkably, the effect of the modified Newton law is seen only in the clustering parameter while the standard structure of the equations is preserved. A comparison of the modified clustering parameter with that of the original clustering parameter is made to visualize the effect of the correction on the time scale of clustering. The possibility of system symmetry breaking is also analyzed by investigating the behavior of the specific heat with increasing system temperature.

The phrase "more is different" is often used to refer to the new, unexpected collective phenomena that can arise when the number of states in a given system is large. In this contribution to the Snowmass 2021 Study, we describe 13 unexpected collective phenomena that can arise when the dark sector contains a large number of states, contrary to the usual assumptions. These 13 take-away lessons stretch across all of the domains of relevance for dark-matter physics, including collider signatures, direct-detection signatures, indirect-detection signatures, new perspectives on dark-matter complementarity, and even unexpected astrophysical and cosmological phenomena that transcend those normally associated with single-component dark-matter scenarios. These lessons -- and the phenomena on which they are based -- thereby illustrate the need to maintain a broad perspective when contemplating the possible signatures and theoretical possibilities associated with non-minimal dark sectors.