Papers for Thursday, Jun 17 2021

We propose a modification of the Natural Inflation (NI) potential in such a way that the spontaneous symmetry breaking scale $f$ can take values less than one (in Planck units). The proposed potential seems simple enough, however, its consequences are difficult to calculate analytically. Therefore, we illustrate the feasibility of the model by considering some numerical examples that easily satisfy the conditions imposed on the observables $n_s$ and $r$ by the most recent observations, while at the same time maintaining the number of e-folds during the inflationary epoch within the expected range.

We present time-series spectroscopy and photometry of Gaia DR2 6097540197980557440, a new deeply-eclipsing hot subdwarf B (sdB) + M dwarf (dM) binary. We discovered this object during the course of the Eclipsing Reflection Effect Binaries from Optical Surveys (EREBOS) project, which aims to find new eclipsing sdB+dM binaries (HW Vir systems) and increase the small sample of studied systems. In addition to the primary eclipse, which is in excess of $\sim$5 magnitudes in the optical, the light curve also shows features typical for other HW Vir binaries such as a secondary eclipse and strong reflection effect from the irradiated, cool companion. The orbital period is 0.127037 d ($\sim$3 hr), falling right at the peak of the orbital period distribution of known HW Vir systems. Analysis of our time-series spectroscopy yields a radial velocity semi-amplitude of $K_{\rm sdB}=100.0\pm2.0\,{\rm km\,s}^{-1}$, which is amongst the fastest line-of-sight velocities found to date for an HW Vir binary. State-of-the-art atmospheric models that account for deviations from local thermodynamic equilibrium are used to determine the atmospheric parameters of the sdB. Although we cannot claim a unique light curve modeling solution, the best-fitting model has an sdB mass of $M_{\rm sdB} = 0.47\pm0.03\,M_{\odot}$ and a companion mass of $M_{\rm dM} = 0.18\pm0.01\,M_{\odot}$. The radius of the companion appears to be inflated relative to theoretical mass-radius relationships, consistent with other known HW Vir binaries. Additionally, the M dwarf is one of the most massive found to date amongst this type of binary.

Observations of hyper-luminous quasars at $z>6$ reveal the rapid growth of supermassive black holes (SMBHs $>10^9 \rm M_{\odot}$) whose origin is still difficult to explain. Their progenitors may have formed as remnants of massive, metal free stars (light seeds), via stellar collisions (medium-weight seeds) and/or massive gas clouds direct collapse (heavy seeds). In this work we investigate for the first time the relative role of these three seed populations in the formation of $z>6$ SMBHs within an Eddington-limited gas accretion scenario. To this aim, we implement in our semi-analytical data-constrained model a statistical description of the spatial fluctuations of Lyman-Werner (LW) photo-dissociating radiation and of metal/dust enrichment. This allows us to set the physical conditions for BH seeds formation, exploring their relative birth rate in a highly biased region of the Universe at $z>6$. We find that the inclusion of medium-weight seeds does not qualitatively change the growth history of the first SMBHs: although less massive seeds ($<10^3 \rm M_\odot$) form at a higher rate, the mass growth of a $\sim 10^9 \rm M_\odot$ SMBH at $z<15$ is driven by efficient gas accretion (at a sub-Eddington rate) onto its heavy progenitors ($10^5 \rm M_\odot$). This conclusion holds independently of the critical level of LW radiation and even when medium-weight seeds are allowed to form in higher metallicity galaxies, via the so-called super-competitive accretion scenario. Our study suggests that the genealogy of $z \sim 6$ SMBHs is characterized by a rich variety of BH progenitors, which represent only a small fraction ($< 10 - 20\%$) of all the BHs that seed galaxies at $z > 15$.

The nuclear obscurer of Active Galactic Nuclei (AGN) is poorly understood in terms of its origin, geometry and dynamics. We investigate whether physically motivated geometries emerging from hydro-radiative simulations can be differentiated with X-ray reflection spectroscopy. For two new geometries, the radiative fountain model of Wada (2012) and a warped disk, we release spectral models produced with the ray tracing code XARS. We contrast these models with spectra of three nearby AGN taken by NuSTAR and Swift/BAT. Along heavily obscured sight-lines, the models present different 4-20keV continuum spectra. These can be differentiated by current observations. Spectral fits of the Circinus Galaxy favor the warped disk model over the radiative fountain, and clumpy or smooth torus models. The necessary reflector (NH>10^25/cm^2) suggests a hidden population of heavily Compton-thick AGN amongst local galaxies. X-ray reflection spectroscopy is a promising pathway to understand the nuclear obscurer in AGN.

Tracing dust in small dense molecular cores is a powerful tool to study the conditions required for ices to form during the pre-stellar phase. To study these environments, five molecular cores were observed: three with ongoing low-mass star formation (B59, B335, and L483) and two starless collapsing cores (L63 and L694-2). Deep images were taken in the infrared JHK bands with the United Kingdom Infrared Telescope (UKIRT) WFCAM (Wide Field Camera) instrument and IRAC channels 1 and 2 on the Spitzer Space Telescope. These five photometric bands were used to calculate extinction along the line of sight toward background stars. After smoothing the data, we produced high spatial resolution extinction maps ($\sim$13-29") . The maps were then projected into the third dimension using the AVIATOR algorithm implementing the inverse Abel transform. The volume densities of the total hydrogen were measured along lines of sight where ices (H$_2$O, CO, and CH$_3$OH) have previously been detected. We find that lines of sight with pure CH$_3$OH or a mixture of CH$_3$OH with CO have maximum volume densities above 1.0$\times$10$^5$ cm$^{-3}$. These densities are only reached within a small fraction of each of the cores ($\sim$0.3-2.1%). CH$_3$OH presence may indicate the onset of complex organic molecule formation within dense cores and thus we can constrain the region where this onset can begin. The maximum volume densities toward star-forming cores in our sample ($\sim$1.2-1.7$\times$10$^6$ cm$^{-3}$) are higher than those toward starless cores ($\sim$3.5-9.5$\times$10$^5$ cm$^{-3}$).

We compare the sizes and luminosities of faint $z=6$-8 galaxies magnified by the Hubble Frontier Fields (HFF) clusters with star-forming regions, as well as more evolved objects, in the nearby universe. Our high-redshift comparison sample includes 333 z=6-8 galaxies, for which size measurements were made as part of a companion study where lensing magnifications were estimated from various public models. Accurate size measurements for these sources are complicated by the lens model uncertainties, but other results and arguments suggest that faint galaxies are small, as discussed in a companion study. The measured sizes for sources in our comparison sample range from <50 pc to ~500 pc. For many of the lowest luminosity sources, extremely small sizes are inferred, reaching individual sizes as small as 10-30 pc, with several sources in the 10-15 pc range with our conservative magnification limits. The sizes and luminosities are similar to those of single star cluster complexes like 30 Doradus in the lower-redshift universe and -- in a few cases -- super star clusters. The identification of these compact, faint star-forming sources in the z~6-8 universe also allows us to set upper limits on the proto-globular cluster LF at z~6. By comparisons of the counts and sizes with recent models, we rule out (with some caveats) proto-globular cluster formation scenarios favoring substantial (xi=10) post-formation mass loss and set useful upper limits on others. Our size results suggest we may be very close to discovering a bona-fide population of forming globular clusters at high redshift.

We present observations of ionized gas outflows in eleven z$ =1.39-2.59$ radio-loud quasar host galaxies. Data was taken with the integral field spectrograph (IFS) OSIRIS and the adaptive optics system at the W.M. Keck Observatory targeting nebular emission lines (H$\beta$, [OIII], H$\alpha$, [NII] and [SII]) redshifted into the near-infrared (1-2.4 \micron). Outflows with velocities of 500 - 1700 km\,s$^{-1}$ are detected in 10 systems on scales ranging from $<1$ kpc to 10 kpc with outflow rates from 8-2400 M$_\odot$yr$^{-1}$. For five sources, the outflow momentum rates are 4-80 times $L_{AGN}$/c, consistent with outflows being driven by an energy conserving shock. The five other outflows are either driven by radiation pressure or an isothermal shock. The outflows are the dominant source of gas depletion, and we find no evidence for star formation along the outflow paths. For eight objects, the outflow paths are consistent with the orientation of the jets. Yet, given the calculated pressures, we find no evidence of the jets currently doing work on these galactic-scale ionized outflows. We find that galactic-scale feedback occurs well before galaxies establish a substantial fraction of their stellar mass, as expected from local scaling relationships.

This work presents an in-depth analysis of four gravitationally lensed red galaxies at z = 1.6-3.2. The sources are magnified by factors of 2.7-30 by foreground clusters, enabling spectral and morphological measurements that are otherwise challenging. Our sample extends below the characteristic mass of the stellar mass function and is thus more representative of the quiescent galaxy population at z > 1 than previous spectroscopic studies. We analyze deep VLT/X-SHOOTER spectra and multi-band Hubble Space Telescope photometry that cover the rest-frame UV-to-optical regime. The entire sample resembles stellar disks as inferred from lensing-reconstructed images. Through stellar population synthesis analysis we infer that the targets are young (median age = 0.1-1.2 Gyr) and formed 80% of their stellar masses within 0.07-0.47 Gyr. Mg II $\lambda\lambda 2796,2803$ absorption is detected across the sample. Blue-shifted absorption and/or redshifted emission of Mg II is found in the two youngest sources, indicative of a galactic-scale outflow of warm ($T\sim10^{4}$ K) gas. The [O III] $\lambda5007$ luminosity is higher for the two young sources (median age less than 0.4 Gyr) than the two older ones, perhaps suggesting a decline in nuclear activity as quenching proceeds. Despite high-velocity ($v\approx1500$ km s$^{-1}$) galactic-scale outflows seen in the most recently quenched galaxies, warm gas is still present to some extent long after quenching. Altogether our results indicate that star formation quenching at high redshift must have been a rapid process (< 1 Gyr) that does not synchronize with bulge formation or complete gas removal. Substantial bulge growth is required if they are to evolve into the metal-rich cores of present-day slow-rotators.

Over the last two decades, asteroseismology has increasingly proven to be the observational tool of choice for the study of stellar physics, aided by the high quality of data available from space-based missions such as CoRoT, Kepler, K2 and TESS. TESS in particular will produce more than an order of magnitude more such data than has ever been available before. While the standard TESS mission products include light curves from 120-sec observations suitable for both exoplanet and asteroseismic studies, they do not include light curves for the vastly larger number of targets observed by the mission at a longer 1800-sec cadence in Full Frame Images (FFIs). To address this lack, the TESS Data for Asteroseismology (T'DA) group under the TESS Asteroseismic Science Consortium (TASC), has constructed an open-source pipeline focused on producing light curves for all stars observed by TESS at all cadences, currently including stars down to a TESS magnitude of 15. The pipeline includes target identification, background estimation and removal, correction of FFI timestamps, and a range of potential photometric extraction methodologies, though aperture photometry is currently the default approach. For the brightest targets, we transparently apply a halo photometry algorithm to construct a calibrated light curve from unsaturated pixels in the image. In this paper, we describe in detail the algorithms, functionality, and products of this pipeline, and summarize the noise metrics for the light curves. Companion papers will address the removal of systematic noise sources from our light curves, and a stellar variability classification from these.

In the BH-galaxy co-evolution framework, most of the star-formation (SF) and the black hole (BH) accretion is expected to take place in highly obscured conditions. Thus, obscured AGN are difficult to identify in optical or X-ray bands, but shine bright in the IR. Moreover, X-ray background (XRB) synthesis models predict that a large fraction of the yet-unresolved XRB is due to the most obscured (Compton thick, CT) of these AGN. In this work, we investigate the synergies between putative IR missions (using SPICA, proposed for ESA/M5 but withdrawn in October 2020, and Origins Space Telescope, OST, as `templates') and the X-ray mission Athena, which should fly in early 2030s, in detecting and characterizing AGN, with a particular focus on the most obscured ones. Using an XRB synthesis model, we estimated the number of AGN and the number of those which will be detected in the X-rays. For each AGN we associated an optical-to-FIR SED from observed AGN with both X-ray data and SED decomposition, and used these SEDs to check if the AGN will be detected by SPICA-like or OST at IR wavelengths. We expect that, with the deepest Athena and SPICA-like (or OST) surveys, we will be able to detect in the IR more than $90\,\%$ of all the AGN (down to L$_{2-10\text{keV}} \sim 10^{42}\,$erg/s and up to $z \sim 10$) predicted by XRB synthesis modeling, and we will detect at least half of them in the X-rays. Athena will be extremely powerful in detecting and discerning moderate- and high-luminosity AGN. We find that the most obscured and elusive CT-AGN will be exquisitely sampled by SPICA-like mission or OST and that Athena will allow a fine characterization of the most-luminous ones. This will provide a significant step forward in the process of placing stronger constraints on the yet-unresolved XRB and investigating the BH accretion rate evolution up to very high redshift ($z \ge 4$).

A free-floating planet is a planetary-mass object that orbits around a non-stellar massive object (e.g. a brown dwarf) or around the Galactic Center. The presence of exomoons orbiting free-floating planets has been theoretically predicted by several models. Under specific conditions, these moons are able to retain an atmosphere capable of ensuring the long-term thermal stability of liquid water on their surface. We model this environment with a one-dimensional radiative-convective code coupled to a gas-phase chemical network including cosmic rays and ion-neutral reactions. We find that, under specific conditions and assuming stable orbital parameters over time, liquid water can be formed on the surface of the exomoon. The final amount of water for an Earth-mass exomonoon is smaller than the amount of water in Earth oceans, but enough to host the potential development of primordial life. The chemical equilibrium time-scale is controlled by cosmic rays, the main ionization driver in our model of the exomoon atmosphere.

Radio relics are the manifestation of electrons presumably being shock (re-)accelerated to high energies in the outskirts of galaxy clusters. However, estimates of the shocks' strength yield different results when measured with radio or X-ray observations. In general, Mach numbers obtained from radio observations are larger than the corresponding X-ray measurements. In this work, we investigate this Mach number discrepancy. For this purpose, we used the cosmological code ENZO to simulate a sample of galaxy clusters that host bright radio relics. For each relic, we computed the radio Mach number from the integrated radio spectrum and the X-ray Mach number from the X-ray surface brightness and temperature jumps. Our analysis suggests that the differences in the Mach number estimates follow from the way in which different observables are related to different parts of the underlying Mach number distribution: radio observations are more sensistive to the high Mach numbers present only in a small fraction of a shock's surface, while X-ray measurements reflect the average of the Mach number distribution. Moreover, X-ray measurements are very sensitive to the relic's orientation. If the same relic is observed from different sides, the measured X-ray Mach number varies significantly. On the other hand, the radio measurements are more robust, as they are unaffected by the relic's orientation.

As the fundamental physical process with many astrophysical implications, the diffusion of cosmic rays (CRs) is determined by their interaction with magnetohydrodynamic (MHD) turbulence. We consider the magnetic mirroring effect arising from MHD turbulence on the diffusion of CRs. Due to the intrinsic superdiffusion of turbulent magnetic fields, CRs with large pitch angles that undergo mirror reflection, i.e., bouncing CRs, are not trapped between magnetic mirrors, but move diffusively along the magnetic field, leading to a new type of parallel diffusion. This diffusion is in general slower than the diffusion of non-bouncing CRs with small pitch angles that undergo gyroresonant scattering. The critical pitch angle at the balance between magnetic mirroring and pitch-angle scattering is important for determining the diffusion coefficients of both bouncing and non-bouncing CRs and their scalings with the CR energy. We find non-universal energy scalings of diffusion coefficients, depending on the properties of MHD turbulence.

It is still an open question how the Solar system is structured beyond 100 au from the Sun. Our understanding of this vast region remains very limited and only recently we have become aware of the existence there of a group of enigmatic bodies known as the extreme trans-Neptunian objects (ETNOs) that have large orbits with perihelia beyond the orbit of Neptune. Four ETNOs -- Sedna, Leleakuhonua, 2012 VP113, and 2013 SY99 -- have perihelia beyond 50 au. The study of the ETNOs may provide much needed information on how this remote region is organized. Here, we apply machine-learning techniques to the sample of 40 known ETNOs to identify statistically significant clusters that may signal the presence of true dynamical groupings and study the distribution of the mutual nodal distances of the known ETNOs that measure how close two orbits can get to each other. Machine-learning techniques show that the known ETNOs may belong to four different populations. Results from the analysis of the distribution of nodal distances show that 41 per cent of the known ETNOs have at least one mutual nodal distance smaller than 1.45 au (1st percentile of the distribution), perhaps hinting at past interactions. In this context, the peculiar pair of ETNOs made of 505478 (2013 UT15) and 2016 SG58 has a mutual ascending nodal distance of 1.35 au at 339 au from the Sun. In addition, the known ETNOs exhibit a highly statistically significant asymmetry between the distributions of object pairs with small ascending and descending nodal distances that might be indicative of a response to external perturbations.

While it has been relatively easy to determine solar-cycle related changes in solar dynamics, determining changes in structure in the deeper layers of the Sun has proved to be difficult. By using helioseismic data obtained over two solar cycles, and sacrificing resolution in favour of lower uncertainties, we show that there are significant changes in the solar convection zone, and perhaps even below it. Using MDI data, we find a relative squared sound-speed difference of $(2.56\pm 0.71)\times 10^{-5}$ at the convection-zone base between the maximum of solar Cycle~23 and the minimum between Cycles~23 and 24. The squared sound-speed difference for the maximum of Cycle~24 obtained with HMI data is $(1.95\pm 0.69)\times 10^{-5}$. GONG data support these results. We also find that the sound speed in the solar convection zone decreases compared to the sound speed below it as the Sun becomes more active. We find evidence of changes in the radial derivative of the sound-speed difference between the solar minimum and other epochs at the base of the convection zone implying possible small changes in the position of the convection-zone base, however, the results are too noisy to make any definitive estimates of the change.

NASA's Kepler, K2 and TESS missions employ Simple Aperture Photometry (SAP) to derive time-series photometry, where an aperture is estimated for each star, and pixels containing each star are summed to create a single light curve. This method is simple, but in crowded fields the derived time-series can be highly contaminated. The alternate method of fitting a Point Spread Function (PSF) to the data is able to account for crowding, but is computationally expensive. In this paper, we present a new approach to extracting photometry from these time-series missions, which fits the PSF directly, but makes simplifying assumptions in order to greatly reduce the computation expense. Our method fixes the scene of the field in each image, estimates the PSF shape of the instrument with a linear model, and allows only source flux and position to vary. We demonstrate that our method is able to separate the photometry from blended targets in the Kepler dataset that are separated by less than a pixel. Our method is fast to compute, and fully accounts for uncertainties from degeneracies due to crowded fields. We name the method described in this work Linearized Field Deblending (LFD). We demonstrate our method on the false positive Kepler target \koi. We are able to separate the photometry of the two sources in the data, and demonstrate the contaminating transiting signal is consistent with a small, sub-stellar companion with a radius of $2.67R_{jup}$ ($0.27R_{sol}$). Our method is equally applicable to extracting photometry from NASA's TESS mission.

Employing a two-parameter model for representing the radiation field, the theory of cosmic-ray acceleration by cyclotron autoresonance is analytically generalized here to include any state of polarization. The equations are derived rigorously and used to investigate the dynamics of the nuclides $_1$H$^1$, $_2$He$^4$, $_{26}$Fe$^{56}$, and $_{28}$Ni$^{62}$, in severe astrophysical conditions. Single-particle calculations and many-particle simulations show that these nuclides can reach ZeV energies ($1 ~ZeV = 10^{21}$ eV) due to interaction with superintense radiation of wavelengths $\lambda=1~$ and $10~ \mu$m, and $\lambda=50$ pm, and magnetic fields of strengths at the mega- and gigatesla levels. Examples employing radiation intensities in the range $10^{32}-10^{42}$ W/m$^2$ are discussed.

Within the MOJAVE VLBA program (Monitoring of Jets in AGN with VLBA Experiments), we have accumulated observational data at 15 GHz for hundreds of jets in $\gamma$-ray bright active galactic nuclei since the beginning of the Fermi scientific observations in August 2008. We investigated a time delay between the flux density of AGN parsec-scale radio emission at 15 GHz and 0.1$-$300 GeV Fermi LAT photon flux, taken from constructed light curves using weekly and adaptive binning. The correlation analysis shows that radio is lagging $\gamma$-ray radiation by up to 8 months in the observer's frame, while in the source frame, the typical delay is about 2-3 months. If the jet radio emission, excluding the opaque core, is considered, significant correlation is found at greater time lags. We supplement these results with VLBI kinematics and core shift data to conclude that the dominant high-energy production zone is typically located within the 15 GHz VLBA core at a distance of a few parsecs from the central nucleus.

We review numerical methods for simulations of cosmic ray (CR) propagation on galactic and larger scales. We present the development of algorithms designed for phenomenological and self-consistent models of CR propagation in kinetic description based on numerical solutions of the Fokker-Planck equation. The phenomenological models assume a stationary structure of the galactic interstellar medium and incorporate diffusion of particles in physical and momentum space together with advection, spallation, production of secondaries and various radiation mechanisms. The self-consistent propagation models of CRs include the dynamical coupling of the CR population to the thermal plasma. The CR transport equation is discretized and solved numerically together with the set of magneto-hydrodynamic (MHD) equations in various approaches treating the CR population as a separate relativistic fluid within the two-fluid approach or as a spectrally resolved population of particles evolving in physical and momentum space. The relevant processes incorporated in self-consistent models include advection, diffusion and streaming well as adiabatic compression and several radiative loss mechanisms. We discuss applications of the numerical models for the interpretation of CR data collected by various instruments. We present example models of astrophysical processes influencing galactic evolution such as galactic winds, the amplification of large-scale magnetic fields and instabilities of the interstellar medium.

Galaxy-galaxy lensing is a powerful probe of the connection between galaxies and their host dark matter halos, which is important both for galaxy evolution and cosmology. We extend the measurement and modeling of the galaxy-galaxy lensing signal in the recent Dark Energy Survey Year 3 cosmology analysis to the highly nonlinear scales ($\sim 100$ kpc). This extension enables us to study the galaxy-halo connection via a Halo Occupation Distribution (HOD) framework for the two lens samples used in the cosmology analysis: a luminous red galaxy sample (redMaGiC) and a magnitude-limited galaxy sample (MagLim). We find that redMaGiC (MagLim) galaxies typically live in dark matter halos of mass $\log_{10}(M_{h}/M_{\odot}) \approx 13.7$ which is roughly constant over redshift ($13.3-13.5$ depending on redshift). We constrain these masses to $\sim 15\%$, approximately $1.5$ times improvement over previous work. We also constrain the linear galaxy bias more than 5 times better than what is inferred by the cosmological scales only. We find the satellite fraction for redMaGiC (MagLim) to be $\sim 0.1-0.2$ ($0.1-0.3$) with no clear trend in redshift. Our constraints on these halo properties are broadly consistent with other available estimates from previous work, large-scale constraints and simulations. The framework built in this paper will be used for future HOD studies with other galaxy samples and extensions for cosmological analyses.

Are critical points important in the Solar Probe Mission? This is a brief discussion of the nature of critical points in solar wind models, what this means physically in the 'real' solar wind, and what can be expected along a nominal Solar Probe Orbit. The conclusion is that the regions where the wind becomes transonic and trans-Alfvenic, which may be irregular and varying, may reveal interesting physics, but the mathematically defined critical points themselves are of less importance.

We present the Time Lag/Delay Reconstructor (TLDR), an algorithm for reconstructing velocity delay maps in the Maximum A Posteriori framework for reverberation mapping. Reverberation mapping is a tomographical method for studying the kinematics and geometry of the broad-line region of active galactic nuclei at high spatial resolution. Leveraging modern image reconstruction techniques, including Total Variation and Compressed Sensing, TLDR applies multiple regularization schemes to re-construct velocity delay maps using the Alternating Direction Method of Multipliers. Along with the detailed description of the TLDR algorithm we present test reconstructions from TLDR applied to synthetic reverberation mapping spectra as well as a preliminary reconstruction of the H\b{eta}feature of Arp 151 from the 2008 Lick Active Galactic Nuclei Monitoring Project.

The Galactic bar plays a critical role in the evolution of the Milky Way's Central Molecular Zone (CMZ), as its potential drives mass inward toward the Galactic Center via gas flows known as dust lanes. To explore the interaction between the CMZ and the dust lanes, we run hydrodynamic simulations in Arepo, modeling the potential of the Milky Way's bar in the absence of gas self-gravity and star formation physics, and we study the flows of mass using Monte Carlo tracer particles. We estimate the efficiency of the inflow via the dust lanes, finding that only about a third (30 +/- 12%) of the dust lanes' mass initially accretes onto the CMZ, while the rest overshoots and accretes later. Given observational estimates of the amount of gas within the Milky Way's dust lanes, this suggests that the true total inflow rate onto the CMZ is 0.8 +/- 0.6 Msun yr$^{-1}$. Clouds in this simulated CMZ have sudden peaks in their average density near apocenter, where they undergo violent collisions with inflowing material. While these clouds tend to counter-rotate due to shear, co-rotating clouds occasionally occur (~52% are strongly counter-rotating, and ~7% are strongly co-rotating of the 44 cloud sample), likely due to the injection of momentum from collisions with inflowing material. We investigate the formation and evolution of these clouds, finding that they are fed by many discrete inflow events, providing a consistent source of gas to CMZ clouds even as they collapse and form stars.

We consider transport properties of the hypernuclear matter in neutron star cores. In particular, we calculate the thermal conductivity, the shear viscosity, and the momentum transfer rates for np$\Sigma^{-}\Lambda e\mu$ composition of dense matter in $\beta$--equilibrium for baryon number densities in the range $0.1-1$~fm$^{-3}$. The calculations are based on baryon interactions treated within the framework of the non-relativistic Brueckner-Hartree-Fock theory. Bare nucleon-nucleon (NN) interactions are described by the Argonne v18 phenomenological potential supplemented with the Urbana IX three-nucleon force. Nucleon-hyperon (NY) and hyperon-hyperon (YY) interactions are based on the \new{NSC97e and NSC97a models} of the Nijmegen group. We find that the baryon contribution to transport coefficients is dominated by the neutron one as in the case of neutron star cores containing only nucleons. In particular, we find that neutrons dominate the total thermal conductivity over the whole range of densities explored and that, due to the onset of $\Sigma^-$ which leads to the deleptonization of the neutron star core, they dominate also the shear viscosity in the high density region, in contrast with the pure nucleonic case where the lepton contribution is always the dominant one.

We present variability analyses of twenty pointed XMM-Newton observations of the high energy peaked TeV blazar PG 1553+113 taken during 2010 to 2018. We found intraday variability in the total X-ray energy range (0.3 -- 10 keV) in 16 out of 19 light curves or a duty cycle of ~84%. A discrete correlation function analysis of the intraday light curves in the soft and hard X-ray bands peaks on zero lag, showing that the emission in hard and soft bands are co-spatial and emitted from the same population of leptons. Red-noise dominates the power spectral density (PSD) of all the LCs although the PSDs have a range of spectral slopes from -2.36 to -0.14. On longer timescales, the optical and UV variability patterns look almost identical and well correlated, as are the soft and hard X-ray bands, but the optical/UV variations are not correlated to those in the X-ray band, indicating that the optical/UV and X-ray emissions are emitted by two different populations of leptons. We briefly discuss physical mechanisms which may be capable of explaining the observed flux and spectral variability of PG 1553+113 on these diverse timescales.

KIC 9406652, one of the recently identified IW And-type dwarf novae, is the best target for studying the tilted disk in cataclysmic variable stars. In a previous paper by Kimura, Osaki, and Kato (2020), we analyzed its Kepler light curves and found that its orbital light curves during the brightening stage were dominated by the reflection effect of the secondary star and varied with the orientation of the tilted disk; the amplitude was maximized at the minimum of the super-orbital signal and the phase of the light maximum shifted to an earlier one with the advance of the super-orbital phase. We argued there that this was the direct evidence of the retrogradely precessing tilted disk as the secondary star acts like a reflecting object. In order to confirm this interpretation, we have performed numerical modeling of orbital light curves in this paper. We have succeeded in reproducing the main characteristics of the observed orbital light curves by a simple model in which the secondary star is irradiated by the tilted disk. We have also constrained the inclination angle, $i$, of the binary system and the tilt angle, $\theta$, of the disk purely from photometric considerations. The best-fitting parameter set is found to be $i \sim$45~deg and $\theta \sim$2.0~deg, respectively. The orbital inclination thus estimated is consistent with that obtained from the spectroscopic considerations within the uncertainty limit. On the other hand, the tilt angle of the disk could be underestimated by using only the semi-amplitude of super-orbital signals.

The intracluster medium (ICM) is expected to experience on average about three passages of weak shocks with low sonic Mach numbers, $M\lesssim 3$, during the formation of galaxy clusters. Both protons and electrons could be accelerated to become high energy cosmic rays (CRs) at such ICM shocks via diffusive shock acceleration (DSA). We examine the effects of DSA by multiple shocks on the spectrum of accelerated CRs by including {\it in situ} injection/acceleration at each shock, followed by repeated re-acceleration at successive shocks in the test-particle regime. For simplicity, the accelerated particles are assumed to undergo adiabatic decompression without energy loss and escape from the system, before they encounter subsequent shocks. We show that in general the CR spectrum is flattened by multiple shock passages, compared to a single episode of DSA, and that the acceleration efficiency increases with successive shock passages. However, the decompression due to the expansion of shocks into the cluster outskirts may reduce the amplification and flattening of the CR spectrum by multiple shock passages. The final CR spectrum behind the last shock is determined by the accumulated effects of repeated re-acceleration by all previous shocks, but it is relatively insensitive to the ordering of the shock Mach numbers. Thus multiple passages of shocks may cause the slope of the CR spectrum to deviate from the canonical DSA power-law slope of the current shock.

Terrestrial planets with large water inventories are likely ubiquitous and will be among the first Earth-sized planets to be characterized with upcoming telescopes. It has previously been argued that waterworlds-particularly those possessing more than 1% H$_2$O-experience limited melt production and outgassing due to the immense pressure overburden of their overlying oceans, unless subject to high internal heating. But an additional, underappreciated obstacle to outgassing on waterworlds is the high solubility of volatiles in high-pressure melts. Here, we investigate this phenomenon and show that volatile solubilities in melts probably prevent almost all magmatic outgassing from waterworlds. Specifically, for Earth-like gravity and oceanic crust composition, oceans or water ice exceeding 10-100 km in depth (0.1-1 GPa) preclude the exsolution of volatiles from partial melt of silicates. This solubility limit compounds the pressure overburden effect as large surface oceans limit both melt production and degassing from any partial melt that is produced. We apply these calculations to Trappist-1 planets to show that, given current mass and radius constraints and implied surface water inventories, Trappist-1f and -1g are unlikely to experience volcanic degassing. While other mechanisms for interior-surface volatile exchange are not completely excluded, the suppression of magmatic outgassing simplifies the range of possible atmospheric evolution trajectories and has implications for interpretation of ostensible biosignature gases, which we illustrate with a coupled model of planetary interior-climate-atmosphere evolution.

EZ Cnc, or EPIC 212182292, is a non-Blazhko RRab variable star located in the field of K2 Campaign 16. Its atmospheric parameters ($T_{\rm eff}$, $\log{g}$, [M/H]) and radial velocities are measured from the 55 high-quality LAMOST medium-resolution spectra. The fundamental frequency of pulsation is derived as $f=1.8323(17)$ d$^{-1}$ from the K2 light curves. The amplitude ratios $R_{21} = 0.5115(15), 0.490(8)$, $R_{31} = 0.3249(20), 0.279(7)$ and Fourier phase differences $\varphi_{21}=2.7550(20), 2.764(16)$, $\varphi_{31}=5.7194(25), 5.719(31)$ are determined from the Fourier decomposition of K2 light curve and LAMOST radial velocity curve, respectively. Through the constraints of the parameters, four optimal models are obtained in a time-dependent turbulent convection model survey for EPIC 212182292. The parameters of EPIC 212182292 are derived as $M=0.48\pm0.03$ M$_{\odot}$, $L = 42\pm2$ L$_{\odot}$, $T_{\rm eff}=6846\pm50$ K, $\log{g}=2.79\pm0.01$ dex, and $Z = 0.006\pm0.002$, respectively. The precisely determined parameters for RRab variable stars like EPIC 212180092 might help to better understand the period-luminosity relationship of RR Lyrae stars.

We present Atacama Large Millimetre Array (ALMA) 2mm continuum observations of a complete and unbiased sample of 99 870micron-selected sub-millimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 microJy/beam. We measure the flux densities for 70 sources, for which we obtain a typical 870micron-to-2mm flux ratio of 14 +/- 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2mm observations, the dust emissivity index is well-constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is $\beta\simeq1.9\pm0.4$, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust grain model predictions. We also find a negative correlation between the dust temperature and $\beta$, similar to low-redshift observational and theoretical studies. Our results indicate that $\beta\simeq2$ in high-redshift dusty star-forming galaxies, implying little evolution in dust grain properties between our SMGs and local dusty galaxy samples, and suggesting these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their ISM.

Interstellar objects (ISOs), the parent population of 1I/Oumuamua and 2I/Borisov, are abundant in the interstellar medium of the Milky Way. This means that the interstellar medium, including molecular cloud regions, has three components: gas, dust, and ISOs. From the observational constraints for the field density of ISOs drifting in the solar neighbourhood, we infer a typical molecular cloud of 10 pc diameter contains some 10$^{18}$ ISOs. At typical sizes ranging from hundreds of metres to tens of km, ISOs are entirely decoupled from the gas dynamics in these molecular clouds. Here we address the question of whether ISOs can follow the collapse of molecular clouds. We perform low-resolution simulations of the collapse of molecular clouds containing initially static ISO populations toward the point where stars form. In this proof-of-principle study, we find that the interstellar objects definitely follow the collapse of the gas --- and many become bound to the new-forming numerical approximations to future stars (sinks). At minimum, 40\% of all sinks have one or more ISO test particles gravitationally bound to them for the initial ISO distributions tested here. This value corresponds to at least $10^{10}$ actual interstellar objects being bound after three initial free-fall times. Thus, ISOs are a relevant component of star formation. We find that more massive sinks bind disproportionately large fractions of the initial ISO population, implying competitive capture of ISOs. Sinks can also be solitary, as their ISOs can become unbound again --- particularly if sinks are ejected from the system. Emerging planetary systems will thus develop in remarkably varied environments, ranging from solitary to richly populated with bound ISOs.

Modern hydrodynamical simulations reproduce many properties of the real universe. These simulations model various physical processes, but many of these are included using `subgrid models' due to resolution limits. Although different subgrid models have been successful in modelling the effects of supernovae (SNe) feedback on galactic properties, it remains unclear if, and by how much, these differing implementations affect observable halo gas properties. In this work, we use `zoom-in' cosmological initial conditions of two volumes selected to resemble the Local Group (LG) evolved with both the Auriga and EAGLE galaxy formation models. While the subgrid physics models in both simulations reproduce realistic stellar components of $L^\star$ galaxies, they exhibit different gas properties. Namely, Auriga predicts that the Milky Way (MW) is almost baryonically closed, whereas EAGLE suggests that only half of the expected baryons reside within the halo. Furthermore, EAGLE predicts that this baryon deficiency extends to the LG, ($r \leq 1 \mathrm{~Mpc}$). The baryon deficiency in EAGLE is likely due to SNe feedback at high redshift, which generates halo-wide outflows, with high covering fractions and radial velocities, which both eject baryons and significantly impede cosmic gas accretion. Conversely, in Auriga, gas accretion is almost unaffected by feedback. These differences appear to be the result of the different energy injection methods from SNe to gas. Our results suggest that both quasar absorption lines and fast radio burst dispersion measures could constrain these two regimes with future observations.

Under cold conditions in dense cores, gas-phase molecules and atoms are depleted from the gas-phase to the surface of interstellar grains. Considering the time scales and physical conditions within these cores, a portion of these molecules has to be brought back into the gas-phase to explain their observation by milimeter telescopes. We tested the respective efficiencies of the different mechanisms commonly included in the models. We also tested the addition of sputtering of ice grain mantles via a collision with cosmic rays in the electronic stopping power regime. The ice sputtering induced by cosmic rays has been added to the Nautilus gas-grain model while the other processes were already present. Each of these processes were tested on a 1D physical structure determined by observations in TMC1 cold cores. The resulting 1D chemical structure was also compared to methanol gas-phase abundances observed in these cores. We found that all species are not sensitive in the same way to the non-thermal desorption mechanisms, and the sensitivity also depends on the physical conditions. Thus, it is mandatory to include all of them. Chemical desorption seems to be essential in reproducing the observations for H densities smaller than $4\times 10^4$~cm$^{-3}$, whereas sputtering is essential above this density. The models are, however, systematically below the observed methanol abundances. A more efficient chemical desorption and a more efficient sputtering could better reproduce the observations. In conclusion, the sputtering of ices by cosmic-rays collisions may be the most efficient desorption mechanism at high density (a few $10^4$~cm$^{-3}$ under the conditions studied here) in cold cores, whereas chemical desorption is still required at smaller densities. Additional works are needed on both mechanisms to assess their efficiency with respect to the main ice composition.

Measuring time lags between time-series or lighcurves at different wavelengths from a variable or transient source in astronomy is an essential probe of physical mechanisms causing multiwavelength variability. Time-lags are typically quantified using discrete correlation functions (DCF) which are appropriate for linear relationships. However, in variable sources like X-ray binaries, active galactic nuclei (AGN) and other accreting systems, the radiative processes and the resulting multiwavelength lightcurves often have non-linear relationships. For such systems it is more appropriate to use non-linear information-theoretic measures of causation like mutual information, routinely used in other disciplines. We demonstrate with toy models loopholes of using the standard DCF & show improvements when using the mutual information correlation function (MICF). For non-linear correlations, the latter accurately & sharply identifies the lag components as opposed to the DCF which can be erroneous. Following that we apply the MICF to the multiwavelength lightcurves of AGN NGC 4593. We find that X-ray fluxes lead UVW2 fluxes by ~0.2 days, closer to model predictions from reprocessing by the accretion disk than the DCF estimate. The uncertainties with the current lightcurves are too large though to rule out -ve lags. Additionally, we find another delay component at ~-1 day i.e. UVW2 leading X-rays consistent with inward propagating fluctuations in the accretion disk scenario. This is not detected by the DCF. Keeping in mind the non-linear relation between X-ray & UVW2, this is worthy of further theoretical investigation. From both toy models & real observations, it is clear that the mutual information based estimator is highly sensitive to complex non-linear correlations. With sufficiently high temporal resolution, we will precisely detect each of the lag features corresponding to these correlations.

Forthcoming radio surveys will include full polarisation information, which can be potentially useful for weak lensing observations. We propose a new method to measure the (integrated) gravitational field between a source and the observer, by looking at the angle between the morphology of a radio galaxy and the orientation of the polarisation. For this we use the fact that, while the polarisation of a photon is parallel transported along the photon geodesic, the infinitesimal shape of the source, e.g. its principal axis in the case of an ellipse, is Lie transported. As an example, we calculate the rotation of the shape vector with respect to the polarisation direction which is generated by lensing by a distribution of foreground Schwarzschild lenses. For radio galaxies, the intrinsic morphological orientation of a source and its polarised emission are correlated. It follows that observing both the polarisation and the morphological orientation provides information on both the unlensed source orientation and on the gravitational potential along the line of sight.

We report here the discovery of a hot Jupiter at an orbital period of $3.208666\pm0.000016$ days around TOI-1789 (TYC 1962-00303-1, $TESS_{mag}$ = 9.1) based on the TESS photometry, ground-based photometry, and high-precision radial velocity observations. The high-precision radial velocity observations were obtained from the high-resolution spectrographs, PARAS at Physical Research Laboratory (PRL), India, and TCES at Th\"uringer Landessternwarte Tautenburg (TLS), Germany, and the ground-based transit observations were obtained using the 0.43~m telescope at PRL with the Bessel-$R$ filter. The host star is a slightly evolved ($\log{g_*}$ = $3.939^{+0.024}_{-0.046}$), late F-type ($T_{eff}$ = $5984^{+55}_{-57}$ K), metal-rich star ([Fe/H] = $0.370^{+0.073}_{-0.089}$ dex) with a radius of {\ensuremath{$R_{*}$}} = $2.172^{+0.037}_{-0.035}$ \(R_\odot\) located at a distance of $223.56^{+0.91}_{-0.90}$ pc. The simultaneous fitting of the multiple light curves and the radial velocity data of TOI-1789 reveals that TOI-1789b has a mass of $M_{P}$ = $0.70\pm0.16 $ $M_{J}$, a radius of $R_{P}$ = $1.40^{+0.22}_{-0.13}$ $R_{J}$, and a bulk density of $\rho_P$ = $0.31^{+0.15}_{-0.13}$ g cm$^{-3}$ with an orbital separation of a = $0.04873^{+0.00065}_{-0.0016}$ AU. This puts TOI-1789b in the category of inflated hot Jupiters. It is one of the few nearby evolved stars with a close-in planet. The detection of such systems will contribute to our understanding of mechanisms responsible for inflation in hot Jupiters and also provide an opportunity to understand the evolution of planets around stars leaving the main sequence branch.

We analyze the physical conditions, chemical composition and other properties of the photoionized Herbig-Haro object HH~204 through Very Large Telescope (VLT) echelle spectroscopy and Hubble Space Telescope (\textit{HST}) imaging. We kinematically isolate the high-velocity emission of HH~204 from the emission of the background nebula and study the sub-arcsecond distribution of physical conditions and ionic abundances across the HH object. We find that low and intermediate-ionization emission arises exclusively from gas at photoionization equilibrium temperatures, whereas the weak high-ionization emission from HH~204 shows a significant contribution from higher temperature shock-excited gas. We derive separately the ionic abundances of HH~204, the emission of the Orion Nebula and the fainter Diffuse Blue Layer.In HH~204, the O$^{+}$ abundance determined from Collisional Excited Lines (CELs) matches the one based on Recombination Lines (RLs), while the O$^{2+}$ abundance is very low, so that the oxygen abundance discrepancy is zero. The ionic abundances of Ni and Fe in HH~204 have similar ionization and depletion patterns, with total abundances that are a factor of 3.5 higher than in the rest of the Orion Nebula due to dust destruction in the bowshock. We show that a failure to resolve the kinematic components in our spectra would lead to significant error in the determination of chemical abundances (for instance, 40\% underestimate of O), mainly due to incorrect estimation of the electron density.

Gaussian processes offers a convenient way to perform nonparametric reconstructions of observational data assuming only a kernel which describes the covariance between neighbouring points in a data set. We approach the ambiguity in the choice of kernel in Gaussian processes with two methods -- (a) approximate Bayesian computation with sequential Monte Carlo sampling and (b) genetic algorithm -- in order to address the often ad hoc choice of the kernel and use the overall resulting method to reconstruct the cosmic chronometers and supernovae type Ia data sets. The results have shown that the Mat\'{e}rn$\left( \nu = 5/2 \right)$ kernel emerges on top of the two-hyperparameter family of kernels for both cosmological data sets. On the other hand, we use the genetic algorithm in order to select a most naturally-fit kernel among a competitive pool made up of a ten-hyperparameters class of kernels. Imposing a Bayesian information criterion-inspired measure of the fitness, the results have shown that a hybrid of the radial basis function and the Mat\'{e}rn$\left( \nu = 5/2 \right)$ kernel best represented both data sets.

It is well known that low mass young stellar objects (LMYSOs) gain a significant portion of their final mass through episodes of very rapid accretion, with mass accretion rates up to $\dot M_* \sim 10^{-4} M_{\odot}$~yr$^{-1}$. Recent observations of high mass young stellar objects (HMYSO) with masses $M_* \gtrsim 10 M_{\odot}$ uncovered outbursts with accretion rates exceeding $\dot M_*\sim 10^{-3}M_{\odot}$~yr$^{-1}$. Here we examine which scenarios proposed in the literature so far to explain accretion bursts of LMYSOs can apply to the episodic accretion in HMYSOs. We utilise a 1D time dependent models of protoplanetary discs around HMYSOs to study burst properties. We find that discs around HMYSOs are much hotter than those around their low mass cousins. As a result, much more extended regions of the disc are prone to the thermal hydrogen ionisation and MRI activation instabilities. The former in particular is found to be ubiquitous in a very wide range of accretion rates and disc viscosity parameters. The outbursts triggered by these instabilities, however, always have too low $\dot M_*$, and are one to several orders of magnitude too long compared to those observed from HMYSOs so far. On the other hand, bursts generated by tidal disruptions of gaseous giant planets formed by the gravitational instability of the protoplanetary discs yield properties commensurate with observations, provided that the clumps are in the post-collapse configuration with planet radius $R_{\rm p} \gtrsim 10 $ Jupiter radii. Furthermore, if observed bursts are caused by disc ionisation instabilities then they should be periodic phenomena with the duration of the quiescent phase comparable to that of the bursts. This may yield potentially observable burst periodicity signatures in the jets, the outer disc, or the surrounding diffuse material of massive HMYSOs. (abridged)

The Airborne Infrared Spectrometer (AIR-Spec) was commissioned during the 2017 total solar eclipse, when it observed five infrared coronal emission lines from the Gulfstream V High-performance Instrumented Airborne Platform for Environmental Research (GV HIAPER), a research jet owned by the National Science Foundation (NSF) and operated by the National Center for Atmospheric Research (NCAR). The second AIR-Spec research flight took place during the July 2, 2019 total solar eclipse across the south Pacific. The 2019 eclipse flight resulted in seven minutes of observations, during which the instrument measured all four of its target emission lines: S XI 1.393 $\mu$m, Si X 1.431 $\mu$m, S XI 1.921 $\mu$m, and Fe IX 2.853 $\mu$m. The 1.393 $\mu$m line, half of a density-sensitive S XI line pair, was detected for the first time. The 2017 AIR-Spec detection of Fe IX was confirmed and the first observations were made of the Fe IX intensity as a function of solar radius. Observations of S XI and Si X were used to estimate the temperature and density above the east and west limbs, the subject of a future paper. Atmospheric absorption was significant in the 2019 data, and atmospheric modeling was required to extract accurate line intensities. Telluric absorption features were used to calibrate the wavelength mapping, instrumental broadening, and throughput of the instrument. AIR-Spec underwent significant upgrades in preparation for the 2019 eclipse flight. The thermal background was reduced by a factor of 30, providing a 5.5x improvement in signal-to-noise ratio, and the pointing stability was improved by a factor of five to $<$10 arcsec RMS after image co-alignment. In addition, two imaging artifacts were identified and resolved, making the 2019 data easier to interpret and improving the spectral resolution by up to 50%.

We investigate the coronal magnetic energy and helicity budgets of ten solar ARs, around the times of large flares. In particular, we are interested in a possible relation of the derived quantities to the particular type of the flares that the AR produces, i.e., whether they are associated with a CME or they are confined. Using an optimization approach, we employ time series of 3D nonlinear force-free magnetic field models of ten ARs, covering a time span of several hours around the time of occurrence of large solar flares (GOES class M1.0 and larger). We subsequently compute the 3D magnetic vector potentials associated to the model 3D coronal magnetic field using a finite-volume method. This allows us to correspondingly compute the coronal magnetic energy and helicity budgets, as well as related (intensive) quantities such as the relative contribution of free magnetic energy, $E_{\mathrm{F}}/{E}$ (energy ratio), the fraction of non-potential (current-carrying) helicity, $|H_{\mathrm{J}}|/|{H_{V}}|$ (helicity ratio), and the normalized current-carrying helicity, $|H_{\mathrm{J}}|/{\phi^{\prime}}^{2}$. The total energy and helicity budgets of flare-productive ARs (extensive parameters) cover a broad range of magnitudes, with no obvious relation to the eruptive potential of the individual ARs, i.e., whether or not a CME is produced in association with the flare. The intensive eruptivity proxies, $E_{\mathrm{F}}/{E}$ and $|H_{\mathrm{J}}|/|{H_{V}}|$, and $|H_{\mathrm{J}}|/{\phi^{\prime}}^{2}$, however, seem to be distinctly different for ARs that produced CME-associated large flares compared to those which produced confined flares. For the majority of ARs in our sample, we are able to identify characteristic pre-flare magnitudes of the intensive quantities, clearly associated to subsequent CME-productivity.

We studied a sample of 274 radio and X-ray selected quasars (XQSOs) detected in the COSMOS and XXL-S radio surveys at 3 GHz and 2.1 GHz, respectively. This sample was identified by adopting a conservative threshold in X-ray luminosity, Lx [2-10\ keV] >= 10^44 erg/s, selecting only the most powerful quasars. Using available multiwavelength data, we examined various criteria for the selection of radio-loud (RL) and radio-quiet (RQ) XQSOs, finding that the number of RL/RQ XQSOs changes significantly depending on the chosen criterion. This discrepancy arises due to the different criteria tracing different physical processes and due to our sample being selected from flux-limited radio and X-ray surveys. Another approach to study the origin of radio emission in XQSOs is via their radio luminosity function (RLF). We constructed the XQSO 1.4 GHz radio luminosity functions (RLFs) in six redshift bins at 0.5 <= z <= 3.7. The lower-1.4 GHz luminosity end shows a higher normalization than expected only from AGN contribution in all studied redshift bins. The found "bump" is mostly dominated by emission due to star-forming (SF) processes within the XQSO host galaxies. As expected, AGN-related radio emission dominates at the higher-luminosity end of RLF. The evolution of XQSO RLF was studied via combination of analytic forms from the literature to constrain the lower-luminosity "bump" and the higher-luminosity AGN part of the RLF. We defined two 1.4 GHz luminosity thresholds, L_th,SF and L_th,AGN, below and above which more than 80% of sources contributing to the RLF are dominated by SF and AGN-related activity, respectively. These thresholds evolve with redshift, most likely due to the strong evolution of SFRs of the XQSO host galaxies.

A late (t $\sim$ 1,500 days) multi-wavelength (UV, optical, IR, and X-ray) flare was found in PS1-10adi, a tidal disruption event (TDE) candidate that took place in an active galactic nucleus (AGN). TDEs usually involve super-Eddington accretion, which drives fast mass outflow (disk wind). So here we explore a possible scenario that such a flare might be produced by the interaction of the disk wind with a dusty torus for TDEs in AGN. Due to the high velocity of the disk wind, strong shocks will emerge and convert the bulk of the kinetic energy of the disk wind to radiation. We calculate the dynamics and then predict the associated radiation signatures, taking into account the widths of the wind and torus. We compare our model with the bolometric light curve of the late flare in PS1-10adi constructed from observations. We find from our modeling that the disk wind has a total kinetic energy of about $10^{51}$ erg and a velocity of 0.1 c (i.e., a mass of 0.3 $M_{\odot}$); the gas number density of the clouds in the torus is $3\times 10^{7}$ $\rm cm^{-3}$. Observation of such a late flare can be an evidence of the disk wind in TDEs and can be used as a tool to explore the nuclear environment of the host.

Aims. We report in this paper the test of a plane holographic optical element to be used as an aberration-corrected grating for a slitless spectrograph, inserted in a convergent telescope beam. Our long term objective is the optimisation of a specific hologram to switch the auxiliary telescope imager of the Vera Rubin Observatory into an accurate slitless spectrograph, dedicated to the atmospheric transmission measurement. We present and discuss here the promising results of tests performed with prototype holograms at the CTIO $0.9\,$m telescope during a run of 17 nights in May-June 2017. Methods. After their on-sky geometrical characterisation, the performances of the holograms as aberration-balanced dispersive optical elements have been established by analysing spectra obtained from spectrophotometric standard stars and narrow-band emitter planetary nebulae. Results. Thanks to their additional optical function, our holographic disperser prototypes allow to produce significantly better focused spectra within the full visible wavelength domain $[370,1050]\,$nm than a regular grating, which suffers from strong defocusing and aberrations when used in similar conditions. We show that the resolution of our slitless on-axis spectrograph equipped with the hologram approaches its theoretical performance. Conclusions. While estimating the benefits of an hologram for the spectrum resolution, the roadmap to produce a competitive holographic element for the Vera Rubin Observatory auxiliary telescope has been established.

New MMT/Hectospec spectroscopy centered on the galaxy cluster A2626 and covering a ${\sim} 1.8\,\text{deg}^2$ area out to $z \sim 0.46$ more than doubles the number of galaxy redshifts in this region. The spectra confirm four clusters previously identified photometrically. A2625, which was previously thought to be a close neighbor of A2626, is in fact much more distant. The new data show six substructures associated with A2626 and five more associated with A2637. There is also a highly collimated collection of galaxies and galaxy groups between A2626 and A2637 having at least three and probably four substructures. At larger scales, the A2626--A2637 complex is not connected to the Pegasus--Perseus filament.

The high cosmological precision offered by the next generation of galaxy surveys hinges on improved corrections for Galactic dust extinction. We explore the possibility of estimating both the dust extinction and large-scale structure from a single photometric galaxy survey, making use of the predictable manner in which Milky Way dust affects the measured brightness and colors of galaxies in a given sky location in several redshift bins. To test our method, we use a synthetic catalog from a cosmological simulation designed to model the Vera C. Rubin Observatory Legacy Survey of Space and Time. At high Galactic latitude ($|b|\gtrsim20^\circ$) and a resolution of $1^\circ$ ($7'$), we predict the uncertainty in the measurement of dust extinction, $E(B-V)$, to be $0.005\ \mathrm{mag}$ ($0.015\ \mathrm{mag}$). This is similar to the uncertainty of existing dust maps, illustrating the feasibility of our method. Simultaneous estimation of large-scale structure is predicted to recover the galaxy overdensity $\delta$ with a precision of $\sim0.01$ ($\sim0.05$) at $1^\circ$ ($7'$) resolution. We also introduce a Bayesian formalism that combines prior information from existing dust maps with the likelihood of Galactic dust extinction determined from the excursion of observed galaxy properties.

We perform a comprehensive determination of the systemic proper motions of 74 dwarf galaxies and dwarf galaxy candidates in the Local Group based on Gaia early data release 3. The outputs of the analysis for each galaxy, including probabilities of membership, will be made publicly available. The analysis is augmented by a determination of the orbital properties of galaxies within 500 kpc. We adopt the flexible Bayesian methodology presented by McConnachie \& Venn (2020), which takes into account the location of the stars on the sky, on the colour-magnitude diagram and on the proper motion plane. We apply some modifications, in particular to the way the colour-magnitude diagram and spectroscopic information are factored in, e.g. by including stars in several evolution phases. The bulk motions are integrated in three gravitational potentials: two where the Milky Way is treated in isolation and has a mass 0.9 \& 1.6 $\times 10^{12}$M$_{\odot}$ and the time-varying potential by Vasiliev et al. (2021), which includes the infall of a massive Large Magellanic Cloud (LMC). We are able to determine bulk proper motions for 73 systems, and we consider reliable 66 of these measurements. For the first time, systemic motions are presented for galaxies out to a distance of 1.4 Mpc, in the NGC~3109 association. The inclusion of the infall of a massive LMC significantly modifies the orbital trajectories of the objects, with respect to orbit integration in static Milky Way-only potentials, and leads to 6 galaxies being likely associated to the LMC, 3 possibly associated and 1 recently captured object. We discuss the results of the orbit integration in the context of the relation of the galaxies to the system of Milky Way satellites, implications for the too-big-to-fail problem, impact on star formation histories, and tidal disruption.

Tidal disruption events (TDEs) can uncover the quiescent super-massive black holes (SMBHs) at the center of galaxies and also offer a promising method to study them. After the disruption of a star by a SMBH, the highly elliptical orbit of the debris stream will be gradually circularized due to the self-crossing, and then the circularized debris will form an accretion disk. The recent TDE candidate AT 2019avd has double peaks in its optical light curve, and the X-ray emerges near the second peak. The durations of the peaks are about 400 and 600 days, respectively, and the separation between them is ~ 700 days. We fit and analyse its spectral energy distribution (SED) in optical/UV, mid-infrared, and X-ray bands. We find that this source can be interpreted as the circularization process in the first phase plus the delayed accretion process in the second phase. Under this two-phase scenario, we use the succession of self-crossing circularization model to fit the first peak, and the delayed accretion model to fit the second peak. The fitting results are consistent with the partial disruption of a 0.9 M_sun star by a 7 * 10^6 M_sun SMBH with the penetration factor \beta ~ 0.6. Furthermore, we find the large-amplitude (by factors up to ~ 5) X-ray variability in AT 2019avd can be interpreted as the rigid-body precession of the misaligned disk due to the Lense-Thirring effect of a spinning SMBH, with the disk precession period of 10 - 25 days.

We investigate the impact on convective numerical simulations of thermo-compositional diabatic processes. We focus our study on simulations with a stabilizing temperature gradient and a destabilizing mean-molecular weight gradient. We aim to establish the possibility for a reduced temperature-gradient in such setups. A suite of 3D simulations were conducted using a numerical hydrodynamic code. We used as a simplified test case, a sample region of the secondary atmosphere of a hot rocky exoplanet within which the chemical transition CO+O $\leftrightarrow$ CO$_{2}$ could occur. Newtonian cooling and a chemical source term was used to maintain a negative mean molecular weight gradient. Our results demonstrate that this setup can reduce the temperature gradient, a result which does not converge away with resolution or over time. We also show that the presence of the reduced temperature gradient is a function of the forcing timescales. The above transition leads to a bifurcation of the temperature profile when the chemical forcing is fast, reminiscent of the bifurcation seen in the boiling crisis for steam/liquid convection. With the reduced temperature gradient in these idealized setups, there exists the possibility for an analogy of the reddening (currently observed in the spectra of brown dwarfs) in the spectra of rocky exoplanet atmospheres. Detailed 1D modelling is needed, in order to characterize the equilibrium thermal and compositional gradients, the timescales, and the impact of a realistic equation of state, in order to assess if the regime identified here will develop in realistic situations. This possibility cannot, however, be excluded a priori. This prediction is new for terrestrial atmospheres and represents strong motivation for the use of diabatic models when analysing atmospheric spectra of rocky exoplanets that will be observed with e.g. the James Webb Space Telescope.

Measurements of the clustering of galaxies in Fourier space, and at low wavenumbers, offer a window into the early Universe via the possible presence of scale dependent bias generated by Primordial Non Gaussianites. On such large scales a Newtonian treatment of density perturbations might not be sufficient to describe the measurements, and a fully relativistic calculation should be employed. The interpretation of the data is thus further complicated by the fact that relativistic effects break statistical homogeneity and isotropy and are potentially divergent in the Infra-Red (IR). In this work we compute for the first time the ensemble average of the most used Fourier space estimator in spectroscopic surveys, including all general relativistic (GR) effects, and allowing for an arbitrary choice of angular and radial selection functions. We show that any observable is free of IR sensitivity once all the GR terms, individually divergent, are taken into account, and that this cancellation is a consequence of the presence of the Weinberg adiabatic mode as a solution to Einstein's equations. We then study the importance of GR effects, including lensing magnification, in the interpretation of the galaxy power spectrum multipoles, finding that they are in general a small, less than ten percent level, correction to the leading redshift space distortions term. This work represents the baseline for future investigations of the interplay between Primordial Non Gaussianities and GR effects on large scales and in Fourier space.

[ABRIDGED] AIMS: The MONOS project is collecting information and studying O-type spectroscopic binaries with delta > -20 deg. In this 2nd paper, we tackle the study of the 35 single-line spectroscopic binary (SB1) systems identified in the previous paper of the series (arXiv:1904.11385) by analyzing our data and reviewing the literature orbits of such systems. METHODS: We have measured the radial velocities for the ~700 spectra in our database using two different methods: Gaussian fitting for several diagnostic lines per object and cross-correlation using synthetic spectra. We also explored the TESS database and analyzed the light curves for 31 of the systems. RESULTS: We have confirmed 21 SB1 systems, discarded the binary nature of 6 stars (9 Sge, HD 192 281, HDE 229 232 AB, 68 Cyg, HD 108 and \alpha Cam), and left 6 stars as inconclusive due to lack of data. The remaining two stars are 15 Mon Aa which has been classified as SB2, and Cyg OB2-22 C, for which we find evidence that it is most likely a triple system where the O star is orbiting an eclipsing SB1. We have also recalculated 20 new orbital solutions, including the first spectroscopic orbital solution for V747 Cep. For Cyg OB2-22 C we have obtained new ephemerides but no new orbit.

In the last years, a new generation of large-scale imaging surveys have probed wide field regions around some nearby galaxies at unprecedented low surface brightness regime (~28.0-29.0 mag arcsec^-2). This offers a chance of discovering very faint dwarf satellites by means of visual inspection of these public deep images. We report the first results of a systematic survey of faint dwarf spheroidal galaxies in the vicinity of the bright late-type spiral NGC 253 galaxy by means of a visual inspection of the images taken by the Dark Energy Survey. Three new dwarf galaxies have been discovered in the vicinity of the brightest member of the Sculptor filament, the late-type spiral NGC 253. Assuming they are companions of NGC 253, their total absolute V-magnitudes fall in the -7 to -9 mag range, which is typical for dwarf satellites in the local Universe. The central surface brightness tend to be extremely low for all the discovered dwarfs and fall roughly in the range of 25-26 mag arcsec^-2 in g-band. Using known data on distances and velocities of galaxies, we estimate the total virial mass of the NGC 253 group to be 8 x 10^11 Mo, which gives the virial radius R_200 = 186 kpc and the turn-around radius of 706 kpc. We also discuss the possible existence of a spatially flattened and velocity-correlated satellite system around NGC 253. This large-scale structure is orientated almost edge-on to line of sight. The possible plane of satellites is only 31 kpc thick with the minor-to-major axis ratio of 0.14. Four out of five galaxies with measured velocities follow a common velocity trend similar to those observed in the planes of satellites around the Andromeda and Centaurus A galaxies. However, the small number of galaxies with known velocities prevents to reach a definitive conclusion about the formation scenario of the structure and its possible relation to the surrounding cosmic web.

The X-ray and extreme-ultra-violet (EUV) emissions from the low-mass stars significantly affect the evolution of the planetary atmosphere. It is, however, observationally difficult to constrain the stellar high-energy emission because of the interstellar extinction. In this work, we simulate the XUV (X-ray+EUV) emission from the Sun-like stars by extending the solar coronal heating model that self-consistently solves the surface-to-corona energy transport, turbulent energy dissipation, and coronal thermal response by conduction and radiation with sufficiently high resolution. The simulations are performed for a range of loop lengths and magnetic filling factors at the stellar surface. When applied to the solar corona, our model is found to reproduce the observed solar XUV spectrum below the Lyman edge, which validates the capability of our model in predicting the XUV spectra of other Sun-like stars. The nearly-linear relation between the unsigned magnetic flux and X-ray luminosity is also reproduced self-consistently. From the simulation runs with various loop lengths and filling factors, the following scaling relations are found. $\log L_{\rm EUV} = 9.93 + 0.67 \log L_{\rm X}$, $\log \Phi^{\rm EUV}_{\rm photon} = 20.40 + 0.66 \log L_{\rm X}$, where $L_{\rm EUV}$ and $L_{\rm X}$ are the cgs-unit luminosity in the EUV and X-ray range, respectively, and $\Phi_{\rm photon}^{\rm EUV}$ is the total number of EUV photons emitted per second. This study demonstrates a refined picture of solar and stellar coronal heating and provides the above observable relations that will be useful for estimating the luminosity of the hidden stellar EUV from X-ray observations.

The precision of the parallax measurements by Gaia is unprecedented. As of Gaia Data Release 2, the number of known nearby open clusters has increased. Some of the clusters appear to be relatively close to each other and form aggregates, which makes them interesting objects to study. We study the aggregates of clusters which share several of the assigned member stars in relatively narrow volumes of the phase space. Using the most recent list of open clusters, we compare the cited central parallaxes with the histograms of parallax distributions of cluster aggregates. The aggregates were chosen based on the member stars which are shared by multiple clusters. Many of the clusters in the aggregates have been assigned parallaxes which coincide with the histograms. However, clusters that share a large number of members in a small volume of the phase space display parallax distributions which do not coincide with the values from the literature. This is the result of ignoring a possibility of assigning multiple probabilities to a single star. We propose that this small number of clusters should be analysed anew.

Adopting three physically-motivated scales (micro - meso - macro, which refer to mpc - kpc - Mpc, respectively) is paramount for achieving a unified theory of multiphase active galactic nuclei feeding and feedback, and it represents a keystone for astrophysical simulations and observations in the upcoming years. In order to promote this multi-scale idea, we have decided to adopt an interdisciplinary approach, exploring the possible conceptual similarities between supermassive black hole feeding and feedback cycles and the dynamics occurring in human cancer microenvironment.

When the light from a distant stellar explosion passes very near to a foreground galaxy or cluster, gravitational lensing can cause it to appear as multiple images on the sky. Such strongly-lensed supernovae can be used to constrain the cosmic expansion rate and dark energy models. Achieving these cosmological goals will require many lensed supernovae with precise time delay measurements. Here we report the discovery of a multiply-imaged supernova that will enable a time delay measurement with an uncertainty of $<1\%$. It appeared in an evolved galaxy at z = 1.95, gravitationally lensed by a massive foreground galaxy cluster. It is likely a Type Ia supernova---the explosion of a low-mass stellar remnant, whose light curve can be used to measure cosmic distances. In archival Hubble Space Telescope imaging, three lensed images of the supernova are detected with relative time delays of $<$200 days. We predict a fourth image will appear close to the cluster core in the year 2037$\pm$2. The SN classification and the predicted reappearance time could be improved with further lens modeling and a comprehensive analysis of systematic uncertainties. Observation of the fourth image could provide a time delay precision of $\approx$7 days over an extraordinary 20 year baseline.

Feedback on the ASTRONET Science Vision and Infrastructure Roadmap from the CTA Consortium.

We present a search for gravitational waves from the coalescence of sub-solar mass black hole binaries using data from the first half of Advanced LIGO and Virgo's third observing run. The observation of a sub-solar mass black hole merger is a clear indication of primordial origin; primordial black holes may contribute to the dark matter distribution. We search for black hole mergers where the primary mass is $0.1-7 M_{\odot}$ and the secondary mass is $0.1-1 M_{\odot}$. A variety of models predict the production and coalescence of binaries containing primordial black holes; some involve dynamical assembly which may allow for residual eccentricity to be observed. For component masses $>0.5 M_{\odot}$, we also search for sources in eccentric orbits, measured at a reference gravitational-wave frequency of 10 Hz, up to $e_{10}\sim 0.3$. We find no convincing candidates and place new upper limits on the rate of primordial black hole mergers. The merger rate of 0.5-0.5 (1.0-1.0)~$M_{\odot}$ sources is $<7100~(1200)$ Gpc$^{-3}$yr$^{-1}$. Our limits are $\sim3-4\times$ more constraining than prior analyses.

There are hints suggesting that properties of galaxy populations in dark matter haloes may depend on their large-scale environment. Recent works point out that very low-density environments influence halo occupation distribution (HOD), however there is not a similar analysis focused on high-density environments. Here we use a simulated set of future virialized superstructures (FVS) to analyse the occupation of galaxies in haloes within these high globally dense regions. We use a publicly available simulated galaxy set constructed with a semi-analytical model to identify FVS in the simulation. Then, we computed the HOD within these superstructures for different absolute magnitude thresholds and make several analysis including the comparison to the global HOD results. We study the dependence on the results on properties of the FVS such as density and volume as well as consider the morphology of galaxies. We also analysed the properties of the stellar content of galaxies and the formation time of the haloes inside FVS. We find a significant increase in the HOD inside FVS. This result is present for all absolute magnitude thresholds explored. The effect is larger in the densest regions of FVS, but does not depend on the volume of the superstructure. We also find that the stellar-mass content of galaxies considerably differs inside the superstructures. Low mass haloes have their central and satellite galaxies with a higher stellar mass content (50%), and exhibit mean star ages (20%) older than average. For massive haloes in FVS we find that only the stellar mass of satellite galaxies varies considerably corresponding to a decrease of 50%. We find a significant statistical difference between the formation times of haloes in FVS and the average population. Haloes residing in superstructures formed earlier, a fact that leads to several changes in the HOD and their member galaxy properties.

Intermediate redshifts between galaxy surveys and the cosmic microwave background (CMB) remain unexplored territory. Line intensity mapping (LIM) offers a way to probe the $z\gtrsim 1$ Universe, including the epoch of reionization and the dark ages. Via exact nulling of the lensing kernel, we show that LIM lensing, in combination with galaxy (resp., CMB) lensing, can uniquely probe the $z\gtrsim 1$ (resp., pre-reionization) Universe. However, LIM foregrounds are a key hurdle to this futuristic technique. While continuum foregrounds can be controlled by discarding modes perpendicular to the line of sight (low $k_\parallel$ modes), interloper foregrounds haven't been addressed in the context of LIM lensing. In this paper, we quantify the interloper bias to LIM lensing for the first time, and derive a ''LIM-pair'' estimator which avoids it exactly after cross-correlating with CMB lensing. This new quadratic lensing estimator works by combining two intensity maps in different lines, from the same redshift, whose interlopers are uncorrelated. As a result, this foreground avoidance method is robust to even large changes in the amplitude of the interloper power and non-Gaussianity. The cross-spectrum of the LIM-pair estimator with CMB lensing is thus robust to the currently large theoretical uncertainties in LIM modeling at high redshift.

The Eddington ratio is a key parameter that governs the diversity of quasar properties. It can be scaled with a strong anti-correlation between optical Fe II and [O III] emission. In search of such indicators in the far-UV band, the HST far-UV spectra of 150 low-redshift quasars are analyzed in combination with their optical SDSS counterparts. The strength of Fe II+Fe III 1123 emission is significantly correlated with that of optical Fe II. A moderate correlation may also exist between Fe II 1071 and optical Fe II. The finding opens the possibility that far-UV Fe II emission may serve as a new gauge of the Eddington ratios. The high- and low-ionization lines in the far-UV band display different patterns: for the quasars with higher Eddington ratios, the low-ionization UV lines are stronger, and the high-ionization lines are broader and weaker.

We point out a novel role for the Standard Model neutrino in dark matter phenomenology where the exchange of neutrinos generates a long-range potential between dark matter particles. The resulting dark matter self interaction could be sufficiently strong to impact small-scale structure formation, without the need of any dark force carrier. This is a generic feature of theories where dark matter couples to the visible sector through the neutrino portal. It is highly testable with improved decay rate measurements at future $Z$, Higgs, and $\tau$ factories, as well as precision cosmology.

We use a science-grade Skipper Charge Coupled Device (Skipper-CCD) operating in a low-radiation background environment to develop a semi-empirical model that characterizes the origin of single-electron events in CCDs. We identify, separate, and quantify three independent contributions to the single-electron events, which were previously bundled together and classified as ``dark counts'': dark current, amplifier light, and spurious charge. We measure a dark current, which depends on exposure, of (5.89+-0.77)x10^-4 e-/pix/day, and an unprecedentedly low spurious charge contribution of (1.52+-0.07)x10^-4 e-/pix, which is exposure-independent. In addition, we provide a technique to study events produced by light emitted from the amplifier, which allows the detector's operation to be optimized to minimize this effect to a level below the dark-current contribution. Our accurate characterization of the single-electron events allows one to greatly extend the sensitivity of experiments searching for dark matter or coherent neutrino scattering. Moreover, an accurate understanding of the origin of single-electron events is critical to further progress in ongoing R\&D efforts of Skipper and conventional CCDs.

We introduce Bilby-MCMC, a Markov-Chain Monte-Carlo sampling algorithm tuned for the analysis of gravitational waves from merging compact objects. Bilby-MCMC provides a parallel-tempered ensemble Metropolis-Hastings sampler with access to a block-updating proposal library including problem-specific and machine learning proposals. We demonstrate that learning proposals can produce over a 10-fold improvement in efficiency by reducing the autocorrelation time. Using a variety of standard and problem-specific tests, we validate the ability of the Bilby-MCMC sampler to produce independent posterior samples and estimate the Bayesian evidence. Compared to the widely-used dynesty nested sampling algorithm, Bilby-MCMC is less efficient in producing independent posterior samples and less accurate in its estimation of the evidence. However, we find that posterior samples drawn from the Bilby-MCMC sampler are more robust: never failing to pass our validation tests. Meanwhile, the dynesty sampler fails the difficult-to-sample Rosenbrock likelihood test, over constraining the posterior. For CBC problems, this highlights the importance of cross-sampler comparisons to ensure results are robust to sampling error. Finally, Bilby-MCMC can be embarrassingly and asynchronously parallelised making it highly suitable for reducing the analysis wall-time using a High Throughput Computing environment. Bilby-MCMC may be a useful tool for the rapid and robust analysis of gravitational-wave signals during the advanced detector era and we expect it to have utility throughout astrophysics.

Negative ions have been detected in abundance in recent years by spacecraft across the solar system. These detections were, however, made by instruments not designed for this purpose and, as such, significant uncertainties remain regarding the prevalence of these unexpected plasma components. In this article, the phenomenon of photodetachment is examined and experimentally and theoretically derived cross-sections are used to calculate photodetachment rates for a range of atomic and molecular negative ions subjected to the solar photon spectrum. These rates are applied to negative ions outflowing from Europa, Enceladus, Titan, Dione and Rhea and their trajectories are traced to constrain source production rates and the extent to which negative ions are able to pervade the surrounding space environments. Predictions are also made for further negative ion populations in the outer solar system with Triton used as an illustrative example. This study demonstrates how, at increased heliocentric distances, negative ions can form stable ambient plasma populations and can be exploited by future missions to the outer solar system.

Starting from the evidence that dark matter indeed exists and permeates the entire cosmos, various bounds on its properties can be estimated. Beginning with the cosmic microwave background and large scale structure, we summarize bounds on the ultralight bosonic dark matter (UBDM) mass and cosmic density. These bounds are extended to larger masses by considering galaxy formation and evolution, and the phenomenon of black hole superradiance. We then discuss the formation of different classes of UBDM compact objects including solitons/axion stars and miniclusters. Next, we consider astrophysical constraints on the couplings of UBDM to Standard Model particles, from stellar cooling (production of UBDM) and indirect searches (decays or conversion of UBDM). Throughout, there are short discussions of "hints and opportunities" in searching for UBDM in each area.

Gravitational wave observation has the potential of probing ultralight bosonic fields such as axion. Axion forms a cloud around a rotating black hole (BH) by superradiant instability and should affect the gravitational waveform from binary BHs. On the other hand, considering the cloud associated with a BH in a binary system, tidal interaction depletes the cloud in some cases during the inspiral phase. We made the exhaustive study of cloud depletion numerically in a wide parameter range for equal mass binaries, assuming only the quadrupolar tidal perturbation is at work. We found that clouds can avoid disappearing due to the tidal effect only when $l=1$ mode is the fastest growing mode and when the binary orbit is counter-rotating in the non-relativistic parameter region.

Detection of many compact binary coalescences (CBCs) is one of the primary goals of the present and future ground-based gravitational-wave (GW) detectors. While increasing the detectors' sensitivities will be crucial in achieving this, efficient data analysis strategies can play a vital role. With given computational power in hand, efficient data analysis techniques can expand the size and dimensionality of the parameter space to search for a variety of GW sources. Matched filtering based analyses that depend on modeled signals to produce adequate signal-to-noise ratios for signal detection may miss them if the parameter space is too restrained. Specifically, the CBC search is currently limited to non-precessing binaries only, where the spins of the components are either aligned or anti-aligned to the orbital angular momentum. A hierarchical search for CBCs is thus well motivated. The first stage of this search is performed by matched filtering coarsely sampled data with a coarse template bank to look for candidate events. These candidates are then followed up for a finer search around the vicinity of an event's parameter space. Performing such a search leads to enormous savings in computational cost. Here we report the first successful implementation of the hierarchical search as a PyCBC-based production pipeline to perform a complete analysis of LIGO observing runs. With this, we analyze Advanced LIGO's first and second observing run data. We recover all the events detected by the PyCBC (flat) search in the first GW catalog, GWTC-1, published by the LIGO-Virgo collaboration, with nearly the same significance using a scaled background. In the analysis, we get an impressive factor of 20 speed-up in computation compared to the flat search. With a standard injection study, we show that the sensitivity of the hierarchical search remains comparable to the flat search within the error bars.

In a brief article published in 1931 and expanded in 1935, the Indian astrophysicist Subrahmanyan Chandrasekhar shared an important astronomical discovery where he introduced what is now known as Chandrasekhar limit. This limit establishes the maximum mass that a white dwarf can reach, which is the stellar remnant that is generated when a low mass star has used up its nuclear fuel. The present work has a double purpose. The first is to present a heuristic derivation of the Chandrasekhar limit. The second is to clarify the genesis of the discovery of Chandrasekhar, as well as the conceptual aspects of the subject. The exhibition only uses high school algebra, as well as some general notions of classical physics and quantum theory.

The underlying geometri of spacetime algebra allows one to derive a force by contracting the relativistic generalization of angular momentum, M, with the mass-current, mw, where w is a proper 4-vector velocity. By applying this force to a cosmological object, a repulsive inverse distance-square law is found, which is proportional to the velocity dispersion squared of that structure. It is speculated if this finding may be relevant to the recent suggestion, that such a force may accelerate the expanding universe with no need for a cosmological constant.