Papers for Friday, Oct 07 2022

Distinguishing signal of young gas rich circumstellar disks from stellar signal in near infrared light is a difficult task. Current techniques such as Angular Differential Imaging (ADI) and Polarimetric Differential Imaging (PDI) cope with drawbacks such as self-subtraction. To address these drawbacks we explore Iterative Angular Differential Imaging (IADI) techniques to increase signal throughput in total intensity observations. This work aims to explore the effectiveness of IADI to recover the self-subtracted regions of disks by applying ADI techniques iteratively. To determine the effectiveness of IADI a model of a disk image is made and post-processed with IADI. In addition, masking based on polarimetric images and a signal threshold for feeding back signal are explored. Asymmetries are a very important factor in recovering the disk due to less overlap of the disk in the data set. In some cases, a factor 75 more flux could be recovered with IADI compared to ADI. The Procrustes distance is used to quantify the impact of the algorithm on the scattering phase function. Depending on the level of noise and the ratio between the stellar signal and disk signal, the phase function can be recovered a factor 6.4 in Procrustes distance better than standard ADI. The amplification and smearing of noise over the image due to many iterations did occur and by using binary masks and a dynamic threshold this feedback was mitigated, but it still is a problem in the final pipeline. Lastly observations of protoplanetary disks made with VLT/SPHERE were processed with IADI giving rise to very promising results. While IADI has problems with low signal-to-noise observations due to noise amplification and star reconstruction, higher signal-to-noise observations show promising results with respect to standard ADI.

We use the Magneticum suite of state-of-the-art hydrodynamical simulations to identify cosmic voids based on the watershed technique and investigate their most fundamental properties across different resolutions in mass and scale. This encompasses the distributions of void sizes, shapes, and content, as well as their radial density and velocity profiles traced by the distribution of cold dark matter particles and halos. We also study the impact of various tracer properties, such as their sparsity and mass, and the influence of void merging on these summary statistics. Our results reveal that all of the analyzed void properties are physically related to each other and describe universal characteristics that are largely independent of tracer type and resolution. Most notably, we find that the motion of tracers around void centers is perfectly consistent with linear dynamics, both for individual, as well as stacked voids. Despite the large range of scales accessible in our simulations, we are unable to identify the occurrence of nonlinear dynamics even inside voids of only a few Mpc in size. This suggests voids to be among the most pristine probes of cosmology down to scales that are commonly referred to as highly nonlinear in the field of large-scale structure.

In a General Relativistic framework, Gravitational Waves (GW) and Electromagnetic (EM) waves are expected to respond in the same way to the effects of matter perturbations between the emitter and the observer. A different behaviour might be a signature of alternative theories of gravity. In this work we study the cross-correlation of resolved GW events (from compact objects mergers detected by the Einstein Telescope, either assuming or excluding the detection of an EM counterpart) and EM signals (coming both from the Intensity Mapping of the neutral hydrogen distribution and resolved galaxies from the SKA Observatory), considering weak lensing, angular clustering and their cross term ($\mathrm{L \times C}$) as observable probes. Cross-correlations of these effects are expected to provide promising information on the behaviour of these two observables, hopefully shedding light on beyond GR signatures. We perform a Fisher matrix analysis with the aim of constraining the $\{\mu_0,\eta_0,\Sigma_0\}$ parameters, either opening or keeping fixed the background parameters $\{w_0,w_a\}$. We find that, although lensing-only forecasts provide significantly unconstrained results, the combination with angular clustering and the cross-correlation of all three considered tracers (GW, IM, resolved galaxies) leads to interesting and competitive constraints. This offers a novel and alternative path to both multi-tracing opportunities for Cosmology and the Modified Gravity sector.

Observationally, the mysterious fast radio bursts (FRBs) are classified as repeating ones and apparently non-repeating ones. While repeating FRBs cannot be classified into the non-repeating group, it is unknown whether the apparently non-repeating FRBs are actually repeating FRBs whose repetitions are yet to be discovered, or whether they belong to another physically distinct type from the repeating ones. In a series of two papers, we attempt to disentangle this mystery with machine learning methods. In this first paper, we focus on an array of supervised machine learning methods. We train the machine learning algorithms with a fraction of the observed FRBs in the first CHIME/FRB catalog, telling them which ones are apparently non-repeating and which ones are repeating. We then let the trained models predict the repetitiveness of the rest of the FRB data with the observed parameters, and we compare the predictions with the observed repetitiveness. We find that the models can predict most FRBs correctly, hinting towards distinct mechanisms behind repeating and non-repeating FRBs. We also find that the two most important distinguishing factors between non-repeating and repeating FRBs are brightness temperature and rest-frame frequency bandwidth. By applying the trained models back to the entire first CHIME catalog, we further identify some potentially repeating FRBs currently reported as non-repeating. We recommend a list of these bursts as targets for future observing campaigns to search for repeated bursts in a combination with the results presented in Paper II using unsupervised machine learning methods.

Photometric and spectral variability of brown dwarfs probes heterogeneous temperature and cloud distribution and traces the atmospheric circulation patterns. We present a new 42-hr Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectral time series of VHS 1256$-$1257 b, a late L-type planetary-mass companion that has been shown to have one of the highest variability amplitudes among substellar objects. The light curve is rapidly evolving and best-fit by a combination of three sine waves with different periods and a linear trend. The amplitudes of the sine waves and the linear slope vary with wavelength, and the corresponding spectral variability patterns match the predictions by models invoking either heterogeneous clouds or thermal profile anomalies. Combining these observations with previous HST monitoring data, we find that the peak-to-valley flux difference is $33\pm2$% with an even higher amplitude reaching 38% in the $J$ band, the highest amplitude ever observed in a substellar object. The observed light curve can be explained by maps that are composed of zonal waves, spots, or a mixture of the two. Distinguishing the origin of rapid light curve evolution requires additional long-term monitoring. Our findings underscore the essential role of atmospheric dynamics in shaping brown dwarf atmospheres and highlight VHS 1256$-$1257 b as one of the most favorable targets for studying atmospheres, clouds, and atmospheric circulation of planets and brown dwarfs.

The resonantly scattered Lyman-$\alpha$ line illuminates the extended halos of neutral hydrogen in the circumgalactic medium of galaxies. We present integral field Keck Cosmic Web Imager observations of double-peaked, spatially extended Ly$\alpha$ emission in 12 relatively low-mass ($M_{\star} \sim10^9 \, M_{\odot}$) $z\sim2$ galaxies characterized by extreme nebular emission lines. Using individual spaxels and small bins as well as radially binned profiles of larger regions, we find that for most objects in the sample the Ly$\alpha$ blue-to-red peak ratio increases, the peak separation decreases, and the fraction of flux emerging at line center increases with radius. We use new radiative transfer simulations to model each galaxy with a clumpy, multiphase outflow with radially varying outflow velocity, and self-consistently apply the same velocity model to the low ionization interstellar absorption lines. These models reproduce the trends of peak ratio, peak separation and trough depth with radius, and broadly reconcile outflow velocities inferred from Ly$\alpha$ and absorption lines. The galaxies in our sample are well-described by a model in which neutral, outflowing clumps are embedded in a hotter, more highly ionized inter-clump medium (ICM), whose residual neutral content produces absorption at the systemic redshift. The peak ratio, peak separation and trough flux fraction are primarily governed by the line-of-sight component of the outflow velocity, the HI column density, and the residual neutral density in the ICM respectively. Azimuthal asymmetries in the line profile further suggest non-radial gas motions at large radii and variations in the HI column density in the outer halos.

Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$\alpha$ forest provides powerful independent constraints on the impact of AGN feedback. We show that TNG over predicts the number density of absorbers at column densities $N_{\rm HI} < 10^{14}$ cm$^{-2}$ compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducial Simba model including AGN jet feedback provides an excellent fit to the observed CDD of the $z = 0.1$ Lyman-$\alpha$ forest across five orders of magnitude in column density. When removing the jet feedback mode in Simba we recover similar results as TNG, reminiscent of the ionizing "photon underproduction crisis", a problem which arose from simulations lacking efficient heating/ionization of intergalactic medium (IGM) gas on large scales. AGN jets in Simba are high speed, collimated, weakly-interacting with the interstellar medium (via brief hydrodynamic decoupling) and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the IGM when compared to TNG kinetic feedback, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low redshift Lyman-$\alpha$ forest provides plausible evidence for long-range AGN jet feedback.

Fast radio bursts (FRBs) are one of the most mysterious astronomical transients. Observationally, they can be classified into repeaters and apparently non-repeaters. However, due to the lack of continuous observations, some apparently repeaters may have been incorrectly recognized as non-repeaters. In a series of two papers, we intend to solve such problem with machine learning. In this second paper of the series, we focus on an array of unsupervised machine learning methods. We apply multiple unsupervised machine learning algorithms to the first CHIME/FRB catalog to learn their features and classify FRBs into different clusters without any premise about the FRBs being repeaters or non-repeaters. These clusters reveal the differences between repeaters and non-repeaters. Then, by comparing with the identities of the FRBs in the observed classes, we evaluate the performance of various algorithms and analyze the physical meaning behind the results. Finally, we recommend a list of most credible repeater candidates as targets for future observing campaigns to search for repeated bursts in combination of the results presented in Paper I using supervised machine learning methods.

We present BIFROST, an extended version of the GPU-accelerated hierarchical fourth-order forward symplectic integrator code FROST. BIFROST (BInaries in FROST) can efficiently evolve collisional stellar systems with arbitrary binary fractions up to $f_\mathrm{bin}=100\%$ by using secular and regularised integration for binaries, triples, multiple systems or small clusters around black holes within the fourth-order forward integrator framework. Post-Newtonian (PN) terms up to order PN3.5 are included in the equations of motion of compact subsystems with optional three-body and spin-dependent terms. PN1.0 terms for interactions with black holes are computed everywhere in the simulation domain. The code has several merger criteria (gravitational-wave inspirals, tidal disruption events and stellar and compact object collisions) with the addition of relativistic recoil kicks for compact object mergers. We show that for systems with $N$ particles the scaling of the code remains good up to $N_\mathrm{GPU} \sim 40\times N / 10^6$ GPUs and that the increasing binary fractions up to 100 per cent hardly increase the code running time (less than a factor $\sim 1.5$). We also validate the numerical accuracy of BIFROST by presenting a number of star clusters simulations the most extreme ones including a core collapse and a merger of two intermediate mass black holes with a relativistic recoil kick.

Context: Early-type galaxies (ETGs) are divided into slow and fast rotators (FRs and SRs) according to the degree of ordered rotation of their stellar populations. Cosmological hydrodynamical simulations indicate that galaxies are formed as FRs before their rotational support decreases, usually because of mergers. Aims: We aimed to investigate this process observationally for galaxies outside of clusters. Methods: We made use of the fact that different merger types leave different traces that have different lifetimes. We statistically analyzed multiple characteristics of galaxies that are expected to be influenced by mergers: tidal features, kinematically distinct cores, stellar age, etc. They were taken from the MATLAS and \atlas databases. We identified through multilinear regression the quantities that, at a fixed mass and environmental density of the galaxy, significantly correlate with a measure of the ordered rotation of the galaxy, \rotsup. Results: We found a negative correlation of the rotational support with the occurrence of tidal disturbances and kinematic substructures and a positive correlation with metallicity and metallicity gradients. For massive galaxies, the rotational support correlates negatively with the abundance of alpha elements, and for the galaxies in low-density environments, it correlates negatively with the central photometric cuspiness. These and additional literature observational constraints are explained the easiest if the mergers that decreased the rotational support of ETGs were typically minor, wet and happening at $z\approx 2$. They did not form the currently observed tidal features. The observed frequency of tidal features implies a merging rate of 0.07-0.2 per Gyr. This is insufficient for explaining the observed growth of radii of ETGs with redshift by mergers.

The launch of NuSTAR and the increasing number of binary black hole (BBH) mergers detected through gravitational wave (GW) observations have exponentially advanced our understanding of black holes. Despite the simplicity owed to being fully described by their mass and angular momentum, black holes have remained mysterious laboratories that probe the most extreme environments in the Universe. While significant progress has been made in the recent decade, the distribution of spin in black holes has not yet been understood. In this work, we provide a systematic analysis of all known black holes in X-ray binary systems (XB) that have previously been observed by NuSTAR, but have not yet had a spin measurement using the "relativistic reflection" method obtained from that data. By looking at all the available archival NuSTAR data of these sources, we measure ten new black hole spins: IGR J17454-2919 -- $a=0.97^{+0.03}_{-0.17}$; GRS 1758-258 -- $a=0.991^{+0.007}_{-0.019}$; MAXI J1727-203 -- $a=0.986^{+0.012}_{-0.159}$; MAXI J0637-430 -- $a=0.97\pm0.02$; Swift J1753.5-0127 -- $a=0.997^{+0.001}_{-0.003}$; V4641 Sgr -- $a=0.86^{+0.04}_{-0.06}$; 4U 1543-47 -- $a=0.98^{+0.01}_{-0.02}$; 4U 1957+11 -- $a=0.95^{+0.02}_{-0.04}$; H 1743-322 -- $a=0.98^{+0.01}_{-0.02}$; MAXI J1820+070 -- $a=0.988^{+0.006}_{-0.028}$ (all uncertainties are at the $1\sigma$ confidence level). We discuss the implications of our measurements on the entire distribution of stellar mass black hole spins in XB, and we compare that with the spin distribution in BBH, finding that the two distributions are clearly in disagreement. Additionally, we discuss the implications of this work on our understanding of how the "relativistic reflection" spin measurement technique works, and discuss possible sources of systematic uncertainty that can bias our measurements.

The search for habitable planets has revealed many planets that can vary greatly from an Earth analog environment. These include highly eccentric orbits, giant planets, different bulk densities, relatively active stars, and evolved stars. This work catalogs all planets found to reside in the HZ and provides HZ boundaries, orbit characterization, and the potential for spectroscopic follow-up observations. Demographics of the HZ planets are compared with a full catalog of exoplanets. Extreme planets within the HZ are highlighted, and how their unique properties may affect their potential habitability. Kepler-296 f is the most eccentric <2 $R_\oplus$ planet that spends 100% of its orbit in the HZ. HD 106270 b and HD 38529 c are the most massive planets (<13 $M_J$) that orbit within the HZ, and are ideal targets for determining the properties of potential hosts of HZ exomoons. These planets, along with the others highlighted, will serve as special edge-cases to the Earth-based scenario and observations of these targets will help test the resilience of habitability outside the standard model. The most promising observational targets are HD 102365 b and 55 Cnc f, and the best candidates that are <2 $R_\oplus$ are GJ 667 C c, Wolf 1061 c, Teegarden's Star b, and Proxima Cen b.

The high-ionization mid-IR lines, excited in the Narrow Line Regions (NLR) of Active Galactic Nuclei (AGN), barely affected by stellar excitation and dust extinction, trace the AGN bolometric power. We used the complete 12 micron sample of Seyfert galaxies, for which 100/116 objects have reliable 2-10keV observations. The [NeV] and [OIV] mid-IR lines linearly correlate with several AGN bolometric indicators (intrinsic 2-10keV and observed 14-195keV X-ray emission, compact nuclear 12 micron emission, [OIII] 5007A line emission), both in terms of flux and luminosity. No evidence of systematic differences in these correlations is found among the Seyfert populations, including type 1 and type 2, and Compton thick and thin AGN. Nevertheless, we find that a sequence of high-to-low Eddington ratio together with strong-to-weak line excitation (traced by the [OIV]/[Ne II] line ratio) encompasses from type 1 through type 2 AGN to Low Ionization Nuclear Emission Line Region (LINER) galaxies, showing intrinsic differences in these three AGN populations. A positive correlation between the Black Hole Accretion Rate (BHAR) and the Star Formation Rate (SFR) is found, but no correlation between the specific SFR (sSFR) and the ratio BHAR/M_(BH), simply reflecting the fact that the more massive is a galaxy, the more it is forming stars and feeding its central black hole. The JWST telescope, just beginning operations, will allow large samples of AGN to be observed in these lines in the nearby Universe (z<0.9).

We describe our first attempt to systematically simulate the solar wind during different phases of the last solar cycle with the Alfv\'en Wave Solar atmosphere Model (AWSoM) developed at the University of Michigan. Key to this study is the determination of the optimal values of one of the most important input parameter of the model, the Poynting flux, which prescribes the energy flux passing through the chromospheric boundary of the model in form of Alfv\'en wave turbulence. It is found that the optimal value of the Poynting flux parameter is correlated with: 1) the open magnetic flux with the linear correlation coefficient of 0.913; 2) the area of the open magnetic field regions with the linear correlation coefficient of 0.946. These highly linear correlations could shed light on understanding how Alfv\'en wave turbulence accelerates the solar wind during different phases of the solar cycle and estimating the Poynting flux parameter for real-time solar wind predictions with AWSoM.

We study the demographics of z ~ 6 broad-line quasars in the black hole (BH) mass-luminosity plane using a sample of more than 100 quasars at 5.7 < z < 6.5. These quasars have well quantified selection functions and nearly one third of them also have virial BH masses estimated from near-IR spectroscopy. We use forward modeling of parameterized intrinsic distributions of BH masses and Eddington ratios, and account for the sample flux limits and measurement uncertainties of the BH masses and luminosities. We find significant differences between the intrinsic and observed distributions of the quantities due to measurement uncertainties and sample flux limits. There is also marginal evidence that the virial BH masses are susceptible to a positive luminosity-dependent bias (BH mass is overestimated when luminosity is above the average), and that the mean Eddington ratio increases with BH mass. Our models provide reliable constraints on the z ~ 6 black hole mass function at M_BH > 10^8.5 M_Sun, with a median 1-sigma uncertainty of ~0.5 dex in abundance. The intrinsic Eddington ratio distribution of M_BH > 10^8.5 M_Sun quasars can be approximated by a mass-dependent Schechter model, with a broad peak around log(L_bol/L_Edd}) ~ -0.9. We also find that, at 4.5 < z < 6, the number densities of more massive BHs tend to decline more rapidly with increasing redshift, contrary to the trend at 2.5 < z < 4.5 reported previously.

In the following short review we will outline some of the possible interaction processes of lower solar atmospheric plasma with the embedded small-scale solar magnetic fields. After introducing the topic, important types of small-scale solar magnetic field elements are outlined to then focus on their creation and evolution, and finally end up describing foremost processes these magnetic fields are involved in, such as the reconnection of magnetic field lines and the creation of magneto-hydrodynamic waves. The occurrence and global coverage in the solar atmosphere of such small-scale phenomena surpass on average those of the more explosive and intense events, mainly related to solar active regions, and therefore their key role as building blocks of solar activity even during the weaker phases of the 11-year solar cycle. In particular, understanding the finest ingredients of solar activity from the lower to the upper solar atmosphere could be determinant to fully understand the heating of the solar corona, which stands out as one of the most intriguing problems in astrophysics nowadays.

We present a third epoch of Chandra observations of the Type Ia Large Magellanic Cloud Supernova remnant (SNR) 0509-67.5. With these new observations from 2020, the baseline for proper motion measurements of the expansion has grown to 20 years (from the earliest Chandra observations in 2000). We report here the results of these new expansion measurements. The lack of nearby bright point sources renders absolute image alignment difficult. However, we are able to measure the average expansion of the diameter of the remnant along several projection directions. We find that the remnant is expanding with an average velocity of 6120 (4900 -- 7360) km s$^{-1}$. This high shock velocity is consistent with previous works, and also consistent with the inference that 0509-67.5 is expanding into a very low density surrounding medium. At the distance of the LMC, this velocity corresponds to an undecelerated age of 600 yrs, with the real age somewhat smaller.

We study the broadband emission of the TeV blazar Mrk501 using multi-wavelength (MWL) observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. During this period, Mrk501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Despite the low activity, significant flux variations are detected at all wavebands, with the highest variations occurring at X-rays and VHE $\gamma$-rays. A significant correlation (>3$\sigma$) between X-rays and VHE $\gamma$-rays is measured, supporting leptonic scenarios to explain the variable parts of the spectral energy distribution (SED), also during low activity states. Extending our data set to 12-years (from 2008 to 2020), we find significant correlations between X-rays and HE $\gamma$-rays, indicating, for the first time, a common physical origin driving the variability between these two bands. We additionally find a correlation between HE $\gamma$-rays and radio, with the radio emission lagging the HE $\gamma$-ray emission by more than 100 days. This is consistent with the $\gamma$-ray emission zone being located upstream of the radio-bright regions of the Mrk501 jet. Furthermore, Mrk501 showed a historically low activity in both X-rays and VHE $\gamma$-rays from mid-2017 to mid-2019 with a stable VHE flux (>2TeV) of 5% the emission of the Crab Nebula. The broadband SED of this 2-year long low-state, the potential baseline emission of Mrk501, can be adequately characterized with a one-zone leptonic model, and with (lepto)-hadronic models that fulfill the neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the historically low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock.

We present the discovery and analysis of SN 2022oqm, a Type Ic supernova (SN) detected <1 day after explosion. The SN rises to a blue and short-lived (2 days) initial peak. Early spectral observations of SN 2022oqm show a hot (40,000 K) continuum with high-ionization C and O absorption features at velocities of 4,000 km s$^{-1}$, while its photospheric radius expands at 20,000 km s$^{-1}$, indicating a pre-existing distribution of expanding C/O material, likely ejected around 2 weeks before the explosion. After around 2.5 days, both the spectrum and light curves evolve into those of a typical SN Ic, with line velocities of 10,000 km s$^{-1}$, in agreement with the photospheric radius evolution. The optical light curves reach a second peak around t ~15 days. By t=60 days, the spectrum of SN 2022oqm becomes nearly nebular, displaying strong C II and [Ca II] emission with no detectable [O I] and marking this event as Ca-rich. The early behavior can be explained by $10^{-3}$ solar mass of optically thin circumstellar material (CSM) surrounding either (1) a massive compact progenitor such as a Wolf-Rayet star, (2) a massive stripped progenitor with an extended envelope, or (3) a binary system with a white dwarf. We propose that the early-time light curve is powered by a combination of interaction of the ejecta with the optically thin CSM and shock cooling (in the massive-star scenario), until the radioactive decay of $^{56}$Ni becomes the dominant power source. The observations can be explained by CSM that is optically thick to X-ray photons which are down converted, is optically thick in the lines as seen in the spectra, and is optically thin to visible-light continuum photons that come either from down-converted X-rays or from the shock-heated ejecta. Calculations show that this scenario is self-consistent.

This paper introduces the emerging field of astrobotics, that is, a recently-established branch of robotics to be of service to astrophysics and observational astronomy. We first describe a modern requirement of dark matter studies, i.e., the generation of the map of the observable universe, using astrobots. Astrobots differ from conventional two-degree-of-freedom robotic manipulators in two respects. First, the dense formation of astrobots give rise to the extremely overlapping dynamics of neighboring astrobots which make them severely subject to collisions. Second, the structure of astrobots and their mechanical specifications are specialized due to the embedded optical fibers passed through them. We focus on the coordination problem of astrobots whose solutions shall be collision-free, fast execution, and complete in terms of the astrobots' convergence rates. We also illustrate the significant impact of astrobots assignments to observational targets on the quality of coordination solutions To present the current state of the field, we elaborate the open problems including next-generation astrophysical projects including 20,000 astrobots, and other fields, such as space debris tracking, in which astrobots may be potentially used

Astrobot swarms are used to capture astronomical signals to generate the map of the observable universe for the purpose of dark energy studies. The convergence of each swarm in the course of its coordination has to surpass a particular threshold to yield a satisfactory map. The current coordination methods do not always reach desired convergence rates. Moreover, these methods are so complicated that one cannot formally verify their results without resource-demanding simulations. Thus, we use support vector machines to train a model which can predict the convergence of a swarm based on the data of previous coordination of that swarm. Given a fixed parity, i.e., the rotation direction of the outer arm of an astrobot, corresponding to a swarm, our algorithm reaches a better predictive performance compared to the state of the art. Additionally, we revise our algorithm to solve a more generalized convergence prediction problem according to which the parities of astrobots may differ. We present the prediction results of a generalized scenario, associated with a 487-astrobot swarm, which are interestingly efficient and collision-free given the excessive complexity of this scenario compared to the constrained one.

Moons orbiting rocky exoplanets in compact orbits about other stars experience an accelerated tidal evolution, and can either merge with their parent planet or reach the limit of dynamical instability within a Hubble time. We review the parameter space over which moons become unbound, including the effects of atmospheric tides on the planetary spin. We find that such tides can change the final outcome from merger to escape, albeit over a limited parameter space. We also follow the further evolution of unbound moons, and demonstrate that the overwhelmingly most likely long-term outcome is that the unbound moon returns to collide with its original parent planet. The dust released by such a collision is estimated to reach optical depths approximately 0.001, exhibit characteristic temperatures of a few hundred degrees Kelvin, and last for a few thousand years. These properties make such events an attractive model for the emerging class of middle-aged main sequence stars that are observed to show transient clouds of warm dust. Furthermore, a late collision between a planet and a returning moon on a hyperbolic orbit may sterilise an otherwise habitable planet.

Context. Photospheric bright points (BPs), as the smallest magnetic element of the photosphere and the footpoint tracer of the magnetic flux tube, are of great significance to the study of BPs. Compared with the study of the characteristics and evolution of a few specific BPs, the study of BPs groups can provide us with a better understanding of the characteristics and overall activities of BPs groups. Aims. We aim to find out the evolution characteristics of the brightness and number of BPs groups at different brightness levels, and how these characteristics differ between quiet and active regions. Methods. We propose a hybrid BPs detection model (HBD Model) combining traditional technology and neural network. The Model is used to detect and calculate the BPs brightness characteristics of each frame of continuous high resolution image sequences of active and quiet regions in TiO-band of a pair of BBSO. Using machine learning clustering method, the PBs of each frame was divided into four levels groups (level1-level4) according to the brightness from low to high. Finally, Fourier transform and inverse Fourier transform are used to analyze the evolution of BPs brightness and quantity in these four levels groups. Results. The activities of BPs groups are not random and disorderly. In different levels of brightness, their quantity and brightness evolution show complex changes. Among the four levels of brightness, BPs in the active region were more active and intense than those in the quiet region. However, the quantity and brightness evolution of BPs groups in the quiet region showed the characteristics of large periodic changes and small periodic changes in the medium and high brightness levels (level3 and level4). The brightness evolution of PBs group in the quiet region has obvious periodic changes, but the active region is in a completely random and violent fluctuation state.

Local early-type galaxies with directly-measured black hole masses, $M_{\rm bh}$, have been reported to represent a biased sample relative to the population at large. Such galaxies with Spitzer Space Telescope imaging have been purported to possess velocity dispersions, $\sigma$, at least $\sim$0.1 dex larger for a given galaxy stellar mass, $M_{\rm *,gal}$, than is typically observed among thousands of early-type galaxies imaged by the Sloan Digital Sky Survey. This apparent offset led Shankar et al. to reduce the normalisation of the observed $M_{\rm bh} \propto \sigma^5$ relation by at least $\sim$0.5 dex to give their "intrinsic relations", including $\sigma$-based modifications to the $M_{\rm bh}$-$M_{\rm *,gal}$ relation. These modifications were based on the untested assumption that the stellar masses had been derived consistently between the two samples. Here, we provide the necessary check using galaxies common to the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) and the Sloan Digital Sky Survey (SDSS). We find that the stellar masses of galaxies with and without directly measured black holes had appeared offset from each other due to the use of inconsistent stellar mass-to-light ratios, $\Upsilon_*$, for the optical and infrared data. We briefly discuss the "intrinsic relations" and why some of these will at times appear to have had partial success when applied to data based on similarly inconsistent values of $\Upsilon_*$. Finally, we reiterate the importance of the $\upsilon$ (lower-case $\Upsilon$) term, which we previously introduced into the $M_{\rm bh}$-$M_{\rm *}$ relations to help avoid $\Upsilon_*$-related mismatches.

We show that chirping gravitational waves in the LIGO frequency band $f=1 - 5000$ Hz can be gravitationally diffracted by the Sun, due to the coincidence of its Fresnel length $r_F \propto \sqrt{1\, {\rm AU}/f}$ and the solar radius $r_\odot$. This solar diffraction is detectable through its frequency-dependent amplification of the wave, albeit with low event rates. Furthermore, we find that solar diffraction allows probing the inner solar profile with the chirping evolution of frequencies. A similar phenomenon can also help discover non-relativistic wave dark matter, as studied in a sequel. This work not only presents an interesting opportunity with ongoing and future LIGO-band missions but also develops the diffractive lensing of long-wavelength waves in the universe.

The origin of spectral curvature at energies $E\simeq 10$ keV in ultraluminous X-ray sources is not well understood. In this paper, we propose a novel mechanism based on synchrotron radiation to explain this cutoff. We show that relativistic plasma can give rise to observed spectral curvature for neutron star magnetic fields due to the variation in the latitude of synchrotron radiation. We analyze the NuSTAR data of two bright pulsar ULXs, NGC 5907 ULX1 and NGC 7793 P13, and provide estimates of the physical parameters of these sources. We fit the data for synchrotron emission at various latitudes and show that the spectral cutoff in these cases can be explained for a large range of acceptable physical parameters, e.g., a semi-relativistic plasma with $\gamma \simeq 10$ and $B\simeq 10^{12} \, \rm G$. We also discuss how such an emission mechanism can be distinguished from other proposed models. A corollary to our study is that most ULXs might be neutron stars as they display such a spectral cutoff.

The acceleration of a jet to relativistic velocities produces a unique memory type gravitational waves (GW) signal: {\it Jet-GW}. I discuss here resent result concerning properties of these GWs and consider their detectability in current and proposed detectors. Classical sources are long and short Gamma-ray bursts as well as hidden jets in core-collapse supernovae. Detection of jet-GWs from these sources will require detectors, such as the proposed BBO, DECIGO and lunar based detectors, that will operate in the deciHz band. The current LVK detectors could detect jet-GWs from a Galactic SGR flare if it is sufficiently asymmetric. Once detected these signals could reveal information concerning jet acceleration and collimation that cannot be explored otherwise.

Thanks to Rosetta orbiter's and Philae lander's data our knowledge of cometary nuclei composition has experienced a great advancement. The properties of 67P/CG nucleus are discussed and compared with other comets explored in the past by space missions. Cometary nuclei are made by a collection of ices, minerals, organic matter, and salts resulting in very dark and red-colored surfaces. When far from the Sun, exposed water and carbon dioxide ices are found only in few locations of 67P/CG where the exposure of pristine subsurface layers or the recondensation of volatile species driven by the solar heating and local terrain morphology can sustain their temporary presence on the surface. The nucleus surface appears covered by a dust layer of variable thickness. Dust grains appear mostly dehydrated and are made by an assemblage of minerals, organic matter, and salts. Spectral analysis shows that the mineral phase is dominated by silicates, fine-grained opaques and ammoniated salts. Aliphatic and aromatic groups, with the presence of the strong hydroxyl group, are identified within the organic matter. The surface composition of cometary nuclei evolves with heliocentric distance and seasonal cycling: approaching perihelion the increase of the solar flux boosts the activity through the sublimation of volatiles which in turn causes the erosion of surface layers, the exposure of ices, the activity in cliffs and pits, the collapse of overhangs and walls, and the mobilization and redistribution of dust. The evolution of color, composition, and texture changes occurring across different morphological regions of the nucleus are correlated with these processes. In this chapter we discuss 67P/CG nucleus composition and evolutionary processes as observed by Rosetta mission in the context of other comets previously explored by space missions or observed from Earth.

In this paper, we use a newly compiled sample of ultra-compact structure in radio quasars and strong gravitational lensing systems with quasars acting as background sources to constrain six spatially flat and non-flat cosmological models ($\Lambda$CDM, PEDE and DGP). These two sets of quasar data (the time-delay measurements of six strong lensing systems and 120 intermediate-luminosity quasars calibrated as standard rulers) could break the degeneracy between cosmological parameters ($H_0$, $\Omega_m$ and $\Omega_k$) and therefore provide more stringent cosmological constraints for the six cosmological models we study. A joint analysis of the quasar sample provides model-independent estimations of the Hubble constant $H_0$, which is strongly consistent with that derived from the local distance ladder by SH0ES collaboration in the $\Lambda$CDM and PEDE model. However, in the framework of a DGP cosmology (especially for the flat universe), the measured Hubble constant is in good agreement with that derived from the the recent Planck 2018 results. In addition, our results show that zero spatial curvature is supported by the current lensed and unlensed quasar observations and there is no significant deviation from a flat universe. For most of cosmological model we study (the flat $\Lambda$CDM, non-flat $\Lambda$CDM, flat PEDE, and non-flat PEDE models), the derived matter density parameter is completely consistent with $\Omega_m\sim 0.30$ in all the data sets, as expected by the latest cosmological observations. Finally, according to the the statistical criteria DIC, although the joint constraints provide substantial observational support to the flat PEDE model, they do not rule out dark energy being a cosmological constant and non-flat spatial hypersurfaces.

Dense gas is the key for understanding star formation in galaxies. We present high resolution ($\sim3''$) observations of CN 2-1 and CS 5-4 as dense gas tracers toward Arp 299, a mid stage major merger of galaxies, with the Submillimeter Array (SMA). The spatial distribution of CN 2-1 and CS 5-4 are generally consistent with each other, as well as HCN 1-0 in literature. However, different line ratios of CS 5-4 and CN 2-1 are found in A, B, and C regions, with highest value in B. Dense gas fraction decreases from IC 694 (A), to NGC 3690 (B) and the overlap starburst region (C and C$'$), which indicates that circum-nuclear upcoming starburst in A and B will be more efficient than that in the overlap region of Arp 299.

Coronal Mass Ejections (CMEs), as they can inject a large amounts of mass and magnetic flux into the interplanetary space, are the primary source of space weather phenomena on the Earth. The present review first briefly introduces the solar surface signatures of the origins of CMEs and then focuses on the attempts to understand the kinematic evolution of CMEs from the Sun to the Earth. CMEs have been observed in the solar corona in white-light from a series of space missions over the last five decades. In particular, LASCO/SOHO has provided almost continuous coverage of CMEs for more than two solar cycles until today. However, the observations from LASCO suffered from projection effects and limited field of view (within 30 Rs from the Sun). The launch in 2006 of the twin STEREO spacecraft made possible multiple viewpoints imaging observations, which enabled us to assess the projection effects on CMEs. Moreover, heliospheric imagers (HIs) onboard STEREO continuously observed the large and unexplored distance gap between the Sun and Earth. Finally, the Earth-directed CMEs that before have been routinely identified only near the Earth at 1 AU in in situ observations from ACE and WIND, could also be identified at longitudes away from the Sun-Earth line using the in situ instruments onboard STEREO. Our review presents the frequently used methods for estimation of the kinematics of CMEs and their arrival time at 1 AU using primarily SOHO and STEREO observations. We emphasize the need of deriving the three-dimensional (3D) properties of Earth-directed CMEs from the locations away from the Sun-Earth line. The results improving the CME arrival time prediction at Earth and the open issues holding back progress are also discussed. Finally, we summarize the importance of heliospheric imaging and discuss the path forward to achieve improved space weather forecasting.

Stellar oscillations can be of topological origin. We reveal this deep and so-far hidden property of stars by establishing a novel parallel between stars and topological insulators. We construct an hermitian problem to derive the expression of the stellar $\mathrm{\textit{acoustic-buoyant}}$ frequency $S$ of non-radial adiabatic pulsations. A topological analysis then connects the changes of sign of the acoustic-buoyant frequency to the existence of Lamb-like waves within the star. These topological modes cross the frequency gap and behave as gravity modes at low harmonic degree $\ell$ and as pressure modes at high $\ell$. $S$ is found to change sign at least once in the bulk of most stellar objects, making topological modes ubiquitous across the Hertzsprung-Russel diagram. Some topological modes are also expected to be trapped in regions where the internal structure varies strongly locally.

The role of electron captures by nuclei in the shallow heating of magnetars is further investigated using both nuclear measurements and the theoretical atomic mass table HFB-27. Starting from the composition of the outer crust in full equilibrium, we have calculated the onset of electron captures and the heat released due to the slow decay of the magnetic field. Numerical results are found to be similar to those previously obtained with the HFB-24 atomic mass model and are consistent with neutron-star cooling data.

We report VLBI monitoring observations of the 22 GHz water (H$_{2}$O) masers around the Mira variable BX Cam, which were carried out as a part of the EAVN Synthesis of Stellar Maser Animations (ESTEMA) project. Data of 37 epochs in total were obtained from 2018 May to 2021 June with a time interval of 3-4 weeks, spanning approximately three stellar pulsation periods ($P= \sim$440 d). In particular, the dual-beam system equipped on the VERA stations was used to measure the kinematics and parallaxes of the H$_{2}$O maser features. The measured parallax, $\pi=1.79\pm 0.08$ mas, is consistent with $Gaia$ EDR3 and previously measured VLBI parallaxes within a 1-$\sigma$ error level. The position of the central star was estimated, based on both the $Gaia$ EDR3 data and the center position of the ring-like 43 GHz silicon-monoxide (SiO) maser distribution imaged with the KVN. The three-dimensional H$_{2}$O maser kinematics indicates that the circumstellar envelope is expanding at a velocity of $13\pm4$ km s$^{-1}$, while there are asymmetries in both the spatial and velocity distributions of the maser features. Furthermore, the H$_{2}$O maser animation achieved by our dense monitoring program manifests the propagation of shock waves in the circumstellar envelope of BX Cam.

The standard $\Lambda$CDM model is recently reported to deviate from the high-redshift Hubble diagram of type Ia supernovae (SNe) and quasars (QSOs) at $\sim4\sigma$ confidence level. In this work, we combine the PAge approximation (a nearly model-independent parameterization) and a high-quality QSO sample to search for the origins of the deviation. By visualizing the $\Lambda$CDM model and the marginalized $3\sigma$ constraints of SNe+QSOs into PAge space, we confirm that the SNe+QSOs constraints in both flat and non-flat PAge cases are in remarkable tension with the standard $\Lambda$CDM cosmology. Next, we investigate the tension from the perspective of redshift-evolution effects. We find that the QSO correlation coefficient $\gamma$ calibrated by SNe+low-z QSOs and SNe+high-z QSOs shows $\sim2.7\sigma$ and $\sim4\sigma$ tensions in flat and non-flat universes, respectively. The tensions for intrinsic dispersion $\delta$ between different data sets are found to be $>4\sigma$ in both flat and non-flat cases. These results indicate that the QSO luminosity correlation suffers from significant redshift evolution and non-universal intrinsic dispersion. Using a redshift-dependence correlation to build QSO Hubble diagram could lead to biases. Thus, the $\sim4\sigma$ deviation from the standard $\Lambda$CDM probably originates from the redshift-evolution effects and non-universal dispersion of the QSO luminosity correlation rather than new physics.

We use three-dimensional radiation hydrodynamic (RHD) simulations to study the formation of massive star clusters under the combined effects of direct ultraviolet (UV) and dust-reprocessed infrared (IR) radiation pressure. We explore a broad range of mass surface density $\Sigma \sim 10^2$-$10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, spanning values typical of weakly star-forming galaxies to extreme systems such as clouds forming super-star clusters, where radiation pressure is expected to be the dominant feedback mechanism. We find that star formation can only be regulated by radiation pressure for $\Sigma \lesssim 10^3 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, but that clouds with $\Sigma \lesssim 10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$ become super-Eddington once high star formation efficiencies ($\sim 80 \%$) are reached, and therefore launch the remaining gas in a steady outflow. These outflows achieve mass-weighted radial velocities of $\sim 15$ - $30 \,\mathrm{km} \, \mathrm{s}^{-1}$, which is $\sim 0.5$ - $2.0$ times the cloud escape speed. This suggests that radiation pressure is a strong candidate to explain recently observed molecular outflows found in young super-star clusters in nearby starburst galaxies. We quantify the relative importance of UV and IR radiation pressure in different regimes, and deduce that both are equally important for $\Sigma \sim 10^3 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, whereas clouds with higher (lower) density are increasingly dominated by the IR (UV) component. Comparison with control runs without either the UV or IR bands suggests that the outflows are primarily driven by the impulse provided by the UV component, while IR radiation has the effect of rendering a larger fraction of gas super-Eddington, and thereby increasing the outflow mass flux by a factor of $\sim 2$.

Sun-as-a-star analyses, in which observational data is spatially integrated, are useful for interpreting stellar data. For future applications to stellar observations, we performed Sun-as-a-star analyses of H$\alpha$ spectra for various active events on the Sun, not only for flares and filament eruptions/surges on the solar disk, but also for eruptions of off limb prominences using H$\alpha$ spectral images taken by the Solar Magnetic Activity Research Telescope / Solar Dynamics Doppler Imager (SMART/SDDI) at Hida Observatory, Kyoto University. All the analyzed events show emission relative to the pre-event state and the changes in their H$\alpha$ equivalent widths are all on the orders of 10$^{-4}$ \r{A}. Sun-as-a-star H$\alpha$ spectra exhibit different features depending on the causes of the emission: (i) Flares show emission at the H$\alpha$ line center, together with red asymmetry and line broadening, as reported in a previous study. (ii) Filament eruptions with and without flares show emission near the H$\alpha$ line center, accompanied by blue-/red-shifted absorption. Notably, disappearance of dark filaments leads to the apparent enhancement of the H$\alpha$ line center emission. (iii) Eruptions of off limb prominences show blue-/red-shifted emission. These spectral features enable us to identify the active phenomena on Sun-like stars. We have also found that even the filament eruptions showing red-shifted absorptions in Sun-as-a-star H$\alpha$ spectra lead to coronal mass ejections (CMEs). This result suggests that even if the falling components of stellar filament eruptions are detected as red-shifted absorptions in H$\alpha$ spectra, such stellar filament eruptions may also develop into CMEs.

We present an updated sample of blue horizontal-branch (BHB) stars selected from the photometric and spectroscopic data from Sloan Digital Sky Survey and its associated project Sloan Extension for Galactic Understanding and Exploration (SEGUE). With this data, we selected candidates for A-type stars in the color-color space and then a mixture modeling technique was implemented in order to distinguish between BHB and main-sequence/blue-straggler stars based on their surface gravity values ($\log \rm{g}$) estimated by the SEGUE Stellar Parameter Pipeline. Our robust approach allows us to attribute individual probabilities of each star truly being in the BHB stage. Hence, our method is advantageous in comparison to previous SEGUE BHB selections that adopted simple $\log \rm{g}$ cuts. We also revisit the color-magnitude relation for these stars and propose two calibrations, based on updated distances for Galactic globular clusters, to estimate absolute magnitudes with $(g-r)_0$ and $(u-r)_0$ colors.

We investigate the ionized gas kinematics relationship with X-ray, radio and accreting properties using a sample of 348 nearby ($z < 0.4$) SDSS-FIRST-X-ray detected AGNs. X-ray properties of our sample are obtained from XMM-$Newton$, $Swift$ and $Chandra$ observations. We unveil the ionized gas outflows in our sample manifested by the non-gravitational broad component in [OIII] $\lambda$5007$\overset{\circ}{A}$. emission line profiles. From the comparison of the correlation of non-parametric outflow velocities (i.e., the velocity width, the maximal velocity of outflow and line dispersion) with X-ray luminosity and radio luminosity, we find that outflow velocities have similarly positive {correlations} with both X-ray and radio luminosity. After correcting for the gravitational component, we find that the [OIII] velocity dispersion normalized by stellar mass also increases with both X-ray luminosity and radio luminosity. We also find that for a given X-ray (radio) luminosity, radio (X-ray) luminous AGNs have higher outflow velocities than non-radio (non-X-ray) luminous AGNs. Therefore, we find no clear preference between X-ray luminosity and radio luminosity in driving high-velocity ionized outflows and conclude that both AGN activity and small-scale jets contribute comparably. Moreover, there is no evidence that our obscured AGNs are preferentially associated with higher velocity outflows. Finally, we find a turning point around log$(\lambda_{Edd}) \simeq -1.3$ when we explore the dependency of outflow velocity on Eddington ratio. It can be interpreted considering the role of high radiation pressure (log$(\lambda_{Edd}) \gtrsim -1.3$) in drastic reduction in the covering factor of the circumnuclear materials.

Context. The physical state of most of the baryonic matter in the local universe is unknown, which is commonly referred to as the ``missing baryon problem". It is theorized that at least half of these missing baryons are in a warm-hot, low-density phase outside of the virialized dark-matter halos. Aims. We make an attempt to find the signature of this warm-hot intergalactic medium (WHIM) phase in the filaments of the nearby Virgo cluster by using optical and Sunyaev-Zeldovich effect data. Methods. Specifically, we use a filament-galaxy catalog created from the HyperLeda database and an all-sky Compton-y map extracted from the Planck satellite data for 2-dimensional cross-correlation analysis by applying spherical harmonics transform. Significance test is based on the null-test simulations which exploits advanced cut-sky analysis tools for a proper map reconstruction. To place upper limits on the WHIM density in the Virgo filaments, realistic baryonic density modelling within the cosmic filaments is done based on state-of-the-art hydro-simulations, and it's done within the signal-boosting routine. Results. The cross-correlation signal is found to be too dim compared to the noise level in the Planck y-map. At 3$\sigma$ confidence level, the upper limit on volume-average WHIM density turns out to be $\left\langle n_e \right\rangle \lt 4\times10^{-4} cm^{-3}$, which is indeed consistent with the WHIM parameter space as predicted from simulations.

An important parameter in the theory of hot accretion flows around black holes is $\delta$, which describes the fraction of ``viscously'' dissipated energy in the accretion flow that directly heats the electrons. The radiative efficiency of a hot accretion flow is determined by its value if other parameters, such as the accretion rate, are determined. Unfortunately, the value of $\delta$ is hard to determine from first principles. The recent Event Horizon Telescope Collaboration (EHTC) results on M87* and Sgr A* provide us with a different way of constraining $\delta$. By combining the mass accretion rates in M87* and Sgr A* estimated by the EHTC with the measured bolometric luminosities of the two sources, we derive good constraints on the radiative efficiencies of the respective accretion flows. In parallel, we use a theoretical model of hot magnetically arrested disks (MAD) to calculate the expected radiative efficiency as a function of $\delta$ (and accretion rate). By comparing the EHTC-derived radiative efficiencies with the theoretical MAD model, we find that Sgr A* requires $\delta > 0.3$, with the most likely value being $\delta \sim 0.5$. A similar comparison in the case of M87* gives inconclusive results, because there is still a large uncertainty in the accretion rate in this source.

Selfsimilarity relations for torsional oscillation frequencies of neutron star crust are discussed. For any neutron star model, the frequencies of fundamental torsional oscillations (with no nodes of radial wave function, i.e. at n=0, and at all possible angular wave numbers l >= 2) is determined by a single constant. Frequencies of ordinary torsional oscillations (at any n>0 with l >= 2) are determined by two constants. These constants are easily calculated through radial integrals over the neutron star crust, giving the simplest method to determine full oscillation spectrum. All constants for a star of fixed mass can be accurately interpolated for stars of various masses (but the same equation of state). In addition, the torsional oscillations can be accurately studied in the flat space-time approximation within the crust. The results can be useful for investigating magneto-elastic oscillations of magnetars which are thought to be observed as quasi-periodic oscillations after flares of soft-gamma repeaters.

Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrounding, geometrically thin disc on this precession, as would be the case in the truncated disc picture of the hard state. To address this, we perform a set of GRMHD simulations of truncated discs misaligned with the spin axes of their central black holes. Our results confirm that the inner-hot flow still undergoes precession, though at a rate that is only 5 percent of what is predicted for an isolated, precessing torus. We find that the exchange of angular momentum between the outer, thin and the inner, thick disc causes this slow-down in the precession rate and discuss its relevance to type-C QPOs.

There is a growing interest in the science uniquely enabled by observations in the MeV range, particularly in light of multi-messenger astrophysics. The Compton Pair (ComPair) telescope, a prototype of the AMEGO Probe-class concept, consists of four subsystems that together detect and characterize gamma rays in the MeV regime. A double-sided strip silicon Tracker gives a precise measure of the first Compton scatter interaction and tracks pair-conversion products. A novel cadmium zinc telluride (CZT) detector with excellent position and energy resolution beneath the Tracker detects the Compton-scattered photons. A thick cesium iodide (CsI) calorimeter contains the high-energy Compton and pair events. The instrument is surrounded by a plastic anti-coincidence (ACD) detector to veto the cosmic-ray background. In this work, we will give an overview of the science motivation and a description of the prototype development and performance.

We present the first clustering measurement of Strong Blended Lyman $\alpha$ (SBLA) absorption systems by measuring their cross-correlation with the Lyman $\alpha$ forest. SBLAs are a new population of absorbers detected within the Lyman $\alpha$ forest. We find a bias of $2.329\pm0.057$, consistent with that of Damped Lyman $\alpha$ absorbers (DLAs). For DLAs, we recover a bias of $2.331\pm0.057$ larger than previously reported (P\'erez-R\`afols et al. 2018b). We also find a redshift space distortion parameter $\beta=0.417\pm0.010$, also consistent with the recovered value for DLAs ($\beta=0.416\pm0.010$). This is consistent with SBLA and DLA systems tracing different portions of the circumgalactic medium of a broadly common population of galaxies. Given these common clustering properties, we combined them to perform a cross-correlation of galaxies in absorption with the Ly$\alpha$ forest. We find that the BAO scale uncertainty of this new measurement is $1.75\times$ that of Ly$\alpha$ auto-correlation and $1.6\times$ that of the quasar cross-correlation with the Ly$\alpha$ forest. We note that the current preferred metal contamination model for fitting the correlation functions with respect to the Ly$\alpha$ forest is not realistic enough for SBLA systems, likely due to their status as high redshift precision sites of high metal enrichment. Mock spectra including SBLA systems and their associated metal absorption are required to understand this sample fully. We conclude that SBLAs have the potential to complement the standard Ly$\alpha$ cosmological analyses in future surveys.

The environment surrounding supermassive black holes (SMBHs) in galactic nuclei (GNs) is expected to harbour stellar-mass binary black hole (BBH) populations. These binaries were suggested to form a hierarchical triple system with the SMBH, and gravitational perturbations from the SMBH can enhance the mergers of BBHs through Lidov-Kozai (LK) oscillations. Previous studies determined the expected binary parameter distribution for this merger channel in single GNs. Here we account for the different spatial distribution and mass distribution models of BBHs around SMBHs and perform direct high-precision regularized N-body simulations, including Post-Newtonian (PN) terms up to order PN2.5, to model merging BBH populations in single GNs. We use a full inspiral-merger-ringdown waveform model of BBHs with nonzero eccentricities and take into account the observational selection effect to determine the parameter distributions of LK-induced BBHs detected with single advanced GW detectors from all GNs in the Universe. We find that the detected mergers' total binary mass distribution is tilted towards lower masses, and the mass ratio distribution is roughly uniform. The redshift distribution peaks between ~0.15-0.55, and the vast majority of binaries merge within redshift ~1.1. The fraction of binaries entering the LIGO/Virgo/KAGRA band with residual eccentricities >0.1 ranges between ~3-12%. We identify a negative correlation between residual eccentricity and mass parameters and a negative correlation between residual eccentricity and source distance. Our results for the parameter distributions and correlations among binary parameters may make it possible to disentangle this merger channel from other BBH merger channels statistically.

We present a three-dimensional (3-D) visualization of jet geometry using numerical methods based on a Markov Chain Monte Carlo (MCMC) and limited memory Broyden-Fletcher-Goldfarb-Shanno (BFGS) optimized algorithm. Our aim is to visualize the 3-D geometry of an active galactic nucleus (AGN) jet using observations, which are inherently two-dimensional (2-D) images. Many AGN jets display complex structures that include hotspots and bends. The structure of these bends in the jet's frame may appear quite different than what we see in the sky frame, where it is transformed by our particular viewing geometry. The knowledge of the intrinsic structure will be helpful in understanding the appearance of the magnetic field and hence emission and particle acceleration processes over the length of the jet. We present the $JetCurry$ algorithm to visualize the jet's 3-D geometry from its 2-D image. We discuss the underlying geometrical framework and outline the method used to decompose the 2-D image. We report the results of our 3-D visualization of the jet of M87, using the test case of the knot D region. Our 3-D visualization is broadly consistent with the expected double helical magnetic field structure of the knot D region of the jet. We also discuss the next steps in the development of the $JetCurry$ algorithm.

We conducted photometric and spectroscopic characterization of near-Earth asteroid (52768) 1998 OR2 during a close approach to the Earth in April of 2020. Our photometric measurements confirm the rotation period of the asteroid to be 4.126 +/- 0.179 hours, consistent with the previously published value of 4.112 +/- 0.001 hours. By combining our visible spectroscopic measurements (0.45 - 0.93 microns) with archival MITHNEOS near infrared spectra (0.78 - 2.49 microns), we classify the asteroid as an Xn-type in the Bus-DeMeo taxonomy. The combined spectrum shows two weak absorption bands: Band I at 0.926 +/- 0.003 microns and Band II at 2.07 +/- 0.02 microns with band depths of 4.5 +/- 0.15% and 4.0 +/- 0.21%, respectively. The band area ratio is 1.13 +/- 0.05. These spectral band parameters plot at the tip of the S(IV) region of the Gaffey S-asteroid subtypes plot suggesting an affinity to ordinary chondrite meteorites. We calculated the chemistry of the olivine and pyroxene using the Band I center to be 20.1 +/- 2.3 mol% fayalite and 18.2 +/- 1.5 mol% ferrosilite, consistent with H chondrites. Principal component analysis of 1998 OR2's combined visible-NIR spectrum fall on the C/X-complex side of the alpha-line, near the end of the shock darkening trend, consistent with its weak absorption bands (band depth < 5%). We use an aerial mixing model with lab measurements of the shock darkened H5 chondrite, Chergach, to constrain the amount of shock darkened material on the asteroid's surface at ~63% dark lithology and ~37% light lithology.

The recent first measurements of the reflection of the surface of a lava world provides an unprecedented opportunity to investigate different stages of rocky planet evolution. The spectral features of the surfaces of rocky lava world exoplanets give insights into their evolution, mantle composition and inner workings. However, no database exists yet that contains spectral reflectivity and emission of a wide range of potential exoplanet surface materials. Here we first synthesized 16 potential exoplanet surfaces, spanning a wide range of chemical compositions based on potential mantle material guided by the metallicity of different host stars. Then we measured their infrared reflection spectrum (2.5 - 28 {\mu}m, 350 - 4000 cm^{-1}), from which we can obtain their emission spectra and establish the link between the composition and a strong spectral feature at 8 {\mu}m, the Christiansen feature (CF). Our analysis suggests a new multi-component composition relationship with the CF, as well as a correlation with the silica content of the exoplanet mantle. We also report the mineralogies of our materials as possibilities for that of lava worlds. This database is a tool to aid in the interpretation of future spectra of lava worlds that will be collected by the James Webb Space Telescope and future missions

The origin of the high-energy astrophysical neutrinos discovered by IceCube remains largely unknown. Multi-messenger studies have indicated that the majority of these neutrinos come from gamma-ray-dark sources. Choked-jet supernovae (cjSNe), which are supernovae powered by relativistic jets stalled in stellar materials, may lead to neutrino emission via photohadronic interactions while the coproduced gamma rays are absorbed. In this paper, we perform an unbinned-maximum-likelihood analysis to search for correlations between IceCube's 10-year muon-track events and our SN Ib/c sample, collected from publicly available catalogs. In addition to the conventional power-law models, we also consider the impacts of more realistic neutrino emission models for the first time, and study the effects of the jet beaming factor in the analyses. Our results show no significant correlation. Even so, the conservative upper limits we set to the contribution of cjSNe to the diffuse astrophysical neutrino flux still allow SNe Ib/c to be the dominant source of astrophysical neutrinos observed by IceCube. We discuss implications to the cjSNe scenario from our results and the power of future neutrino and supernova observations.

Observing Earth-like exoplanets orbiting within the habitable zone of Sun-like stars and studying their atmospheres in reflected starlight requires contrasts of $\sim1\mathrm{e}{-10}$ in the visible. At such high contrast, starlight reflected by exozodiacal dust is expected to be a significant source of contamination. Here, we present high-fidelity simulations of coronagraphic observations of a synthetic Solar System located at a distance of 10 pc and observed with a 12 m and an 8 m circumscribed aperture diameter space telescope operating at 500 nm wavelength. We explore different techniques to subtract the exozodi and stellar speckles from the simulated images in the face-on, the 30 deg inclined, and the 60 deg inclined case and quantify the remaining systematic noise as a function of the exozodiacal dust level of the system. We find that in the face-on case, the exozodi can be subtracted down to the photon noise limit for exozodi levels up to $\sim1000$ zodi using a simple toy model for the exozodiacal disk, whereas in the 60 deg inclined case this only works up to $\sim50$ zodi. We also investigate the impact of larger wavefront errors and larger system distance, finding that while the former have no significant impact, the latter has a strong (negative) impact. Ultimately, we derive a penalty factor as a function of the exozodi level and system inclination that should be considered in exoplanet yield studies as a realistic estimate for the excess systematic noise from the exozodi.

We compare ionised gas and stellar kinematics of 16 star-forming galaxies ($\log(M_\star/M_\odot)=9.7-11.2$, SFR=6-86 $M_\odot/yr$) at $z\sim1$ using near-infrared integral field spectroscopy (IFS) of H$\alpha$ emission from the KMOS$^{\rm 3D}$ survey and optical slit spectroscopy of stellar absorption and gas emission from the LEGA-C survey. H$\alpha$ is dynamically colder than stars, with higher disc rotation velocities (by ~45 per cent) and lower disc velocity dispersions (by a factor ~2). This is similar to trends observed in the local Universe. We find higher rotational support for H$\alpha$ relative to [OII], potentially explaining systematic offsets in kinematic scaling relations found in the literature. Regarding dynamical mass measurements, for six galaxies with cumulative mass profiles from Jeans Anisotropic Multi-Gaussian Expansion (JAM) models the H$\alpha$ dynamical mass models agree remarkably well out to ~10 kpc for all but one galaxy (average $\Delta M_{\rm dyn}(R_{e,\rm F814W})<0.1$ dex). Simpler dynamical mass estimates based on integrated stellar velocity dispersion are less accurate (standard deviation 0.24 dex). Differences in dynamical mass estimates are larger, for example, for galaxies with stronger misalignments of the H$\alpha$ kinematic major axis and the photometric position angle, highlighting the added value of IFS observations for dynamics studies. The good agreement between the JAM models and the dynamical models based on H$\alpha$ kinematics at $z\sim1$ corroborates the validity of dynamical mass measurements from H$\alpha$ IFS observations also for higher redshift rotating disc galaxies.

We present constraints on the flat $\Lambda$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matter field. Our magnitude-limited bright galaxy sample combines 9-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter $S_8 = \sigma_8 \sqrt{\Omega_{\mathrm{m}}/0.3} = 0.773^{+0.028}_{-0.030}$, and the matter density parameter $\Omega_{\mathrm{m}} = 0.290^{+0.021}_{-0.017}$. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint '$3\times2{\mathrm{pt}}$' clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often-excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias and the impact of baryon feedback, this work demonstrates the potential and a way forward to include non-linear scales in cosmological analyses.

We study the possible correspondence between 5-dimensional primordial black holes and massive 5-dimensional KK gravitons as dark matter candidate within the recently proposed dark dimension scenario that addresses the cosmological hierarchy problem. First, we show that in the local universe a population of 5-dimensional black holes with $M_{\rm BH}\sim 7\times 10^{13}~{\rm g}$ would be practically indistinguishable from a KK tower of dark gravitons with $m_{\rm DM} \sim 50~{\rm keV}$. Second, we connect the mass increase of 5-dimensional black holes and the related temperature decrease with the cooling of the tower of massive spin-2 KK excitations of the graviton. The dark gravitons are produced at a mass $\sim 1 - 50~{\rm GeV}$ and the bulk of their mass shifts down to roughly $1 - 100~{\rm keV}$ today. The cooling of the system proceeds via decay to lighter gravitons without losing much total mass density, resembling the intra-tower decays that characterize the cosmological evolution of the dynamical dark matter framework. We associate the intra-tower decays of the graviton gas with the black hole growth through accretion. We also discuss that the primordial black hole $\leftrightharpoons$ dark graviton gas connection can be nicely explained by the bound state picture of black holes in terms of gravitons.

The scheduled launch of the LISA Mission in the next decade has called attention to the gravitational self-force problem. Despite an extensive body of theoretical work, long-time numerical computations of gravitational waves from extreme-mass-ratio-inspirals remain challenging. This work proposes a class of numerical evolution schemes suitable to this problem based on Hermite integration. Their most important feature is time-reversal symmetry and unconditional stability, which enables these methods to preserve symplectic structure, energy, momentum and other Noether charges over long time periods. We apply Noether's theorem to the master fields of black hole perturbation theory on a hyperboloidal slice of Schwarzschild spacetime to show that there exist constants of evolution that numerical simulations must preserve. We demonstrate that time-symmetric integration schemes based on a 2-point Taylor expansion (such as Hermite integration) numerically conserve these quantities, unlike schemes based on a 1-point Taylor expansion (such as Runge-Kutta). This makes time-symmetric schemes ideal for long-time EMRI simulations.

I present a model of inflation and dark energy in which the inflaton potential is constructed by imposing that a scalar field representing the classical energy of the spacetime foam inside the Hubble horizon is an exact solution to the cosmological equations. The resulting potential has the right properties to describe both the early and late expansion epochs of the universe in a unified picture.

The charge conjugation operator in $4k+2$-dimensional space-time posses the peculiar property of preserving chiralities of the spinor fields its acting upon. This nature of $4k+2-$D fermions and its interactions, on compactifying the extra-dimension, leads to 4-dimensional operators with Kaluza Klein (KK) modes at the lowest order, specifically, in the context of baryon and lepton number violation. With the KK mass degeneracy lifted at 1-loop, these operators generate decays/annihilations of KK-1 fermion modes that violate baryon and lepton number. Moreover, including the fermion bilinear interaction with the broken polarization of the hypercharge gauge boson, namely the ``spineless adjoint scalar", these exotic operators generate proton decay and double proton annihilation. In the simplest scenario, this scalar field, being stable, become the Dark Matter (DM) candidate. Here, we discuss the $\Delta B=1 =\Delta L$ and $\Delta B=2=\Delta L$ operators, in $4k+2-$dimensions, in model independent effective field theory and derive the scale of New Physics using $p \to e^+ + \pi$, $DM+p\to DM+ e^+ + \pi$, $p+p \to e^+ +e^+$, $DM+p \to DM+ \bar{p}+e^++e^+$ and $(DM + e^-) + p \to (DM + e^+) + \bar{p} $ processes. Though the terrestrial experiments are not sensitive enough to Dark Matter influenced $\Delta B=2$, their rate can be enhanced substantially at super dense DM clumps and center of galaxy clusters where the Dark Matter distribution is much denser.

We present an ab initio study of the rovibronic spectra of sulfur monoxide ($^{32}$S$^{16}$O) using internally contracted multireference confoguration interaction (ic-MRCI) method and aug-cc-pV5Z basis sets. It covers 13 electronic states $X^{3}\Sigma^{-}$, $a^{1}\Delta$, $b^{1}\Sigma^{+}$, $c^{1}\Sigma^{-}$, $A^{\prime\prime 3}\Sigma^{+}$, $A^{\prime 3}\Delta$, $A^{3}\Pi$, $B^{3}\Sigma^{-}$, $C^{3}\Pi$, $d^{1}\Pi$, $e^{1}\Pi$, $C^{\prime 3}\Pi$, and $(3)^{1}\Pi$ ranging up to 66800 cm$^{-1}$. The ab initio spectroscopic model includes 13 potential energy curves, 23 dipole and transition dipole moment curves, 23 spin-orbit curves, and 14 electronic angular momentum curves. A diabatic representation is built by removing the avoided crossings between the spatially degenerate pairs $C^{3}\Pi - C^{\prime 3}\Pi$ and $e^{1}\Pi - (3)^{1}\Pi$ through a property-based diabatisation method. We also present non-adiabatic couplings and diabatic couplings for these avoided crossing systems. All phases for our coupling curves are defined, and consistent, providing the first fully reproducible spectroscopic model of SO covering the wavelength range longer than 147 nm. Finally, an ab initio rovibronic spectrum of SO is computed.

In this work we examine the effects of a pre-inflationary de Sitter bounce on the energy spectrum of the primordial gravitational waves. Specifically we assume that the Universe is described by several evolution patches, starting with a de Sitter pre-inflationary bounce which is followed by an quasi-de Sitter slow-roll inflationary era, followed by a constant equation of state parameter abnormal reheating era, which is followed by the radiation and matter domination eras and the late-time acceleration eras. The bounce and the inflationary era can be realized by vacuum $f(R)$ gravity and the abnormal reheating and the late-time acceleration eras by the synergy of $f(R)$ gravity and the prefect matter fluids present. Using well-known reconstruction techniques we find which $f(R)$ gravity can realize each evolution patch, except from the matter and radiation domination eras which are realized by the corresponding matter fluids. Accordingly, we calculate the damping factor of the primordial de Sitter bounce, and as we show, the signal can be detected by only one gravitational wave future experiment, in contrast to the case in which the bounce is absent. We discuss in detail the consequences of our results and the future perspectives.

A variational model for the infra-red spectrum of VO is presented which aims to accurately predict the hyperfine structure within the VO $\mathrm{X}\,^4\Sigma^-$ electronic ground state. To give the correct electron spin splitting of the $\mathrm{X}\,^4\Sigma^-$ state, electron spin dipolar interaction within the ground state and the spin-orbit coupling between $\mathrm{X}\,^4\Sigma^-$ and two excited states, $\mathrm{A}\,^4\Pi$ and $\mathrm{1}\,^2\Sigma^+$, are calculated ab initio alongside hyperfine interaction terms. Four hyperfine coupling terms are explicitly considered: Fermi-contact interaction, electron spin-nuclear spin dipolar interaction, nuclear spin-rotation interaction and nuclear electric quadrupole interaction. These terms are included as part of a full variational solution of the nuclear-motion Schr\"odinger equation performed using program DUO, which is used to generate both hyperfine-resolved energy levels and spectra. To improve the accuracy of the model, ab initio curves are subject to small shifts. The energy levels generated by this model show good agreement with the recently derived empirical term values. This and other comparisons validate both our model and the recently developed hyperfine modules in DUO.

We derive a set of no-go theorems and yes-go examples for the parity-odd primordial trispectrum of curvature perturbations. We work at tree-level in the decoupling limit of the Effective Field Theory of Inflation and assume scale invariance and a Bunch-Davies vacuum. We show that the parity-odd scalar trispectrum vanishes in the presence of any number of scalar fields with arbitrary mass and any parity-odd scalar correlator vanishes in the presence of any number of spinning fields with massless de Sitter mode functions, in agreement with the findings of Liu, Tong, Wang and Xianyu [1]. The same is true for correlators with an odd number of conformally-coupled external fields. We derive these results using both the (boostless) cosmological bootstrap, in particular the Cosmological Optical Theorem, and explicit perturbative calculations. We then discuss a series of yes-go examples by relaxing the above assumptions one at the time. In particular, we provide explicit results for the parity-odd trispectrum for (i) violations of scale invariance in single-clock inflation, (ii) the modified dispersion relation of the ghost condensate (non-Bunch-Davies vacuum), and (iii) interactions with massive spinning fields. Our results establish the parity-odd trispectrum as an exceptionally sensitive probe of new physics beyond vanilla inflation.