Papers for Wednesday, Jul 14 2021

In magnetic Cataclysmic Variables (mCVs), X-ray radiation originates from the shock heated multi-temperature plasma in the post-shock region near the white dwarf surface. These X-rays are modified by a complex distribution of absorbers in the pre-shock region. The presence of photo-ionized lines and warm absorber features in the soft X-ray spectra of these mCVs suggests that these absorbers are ionized. We developed the ionized complex absorber model zxipab, which is represented by a power-law distribution of ionized absorbers in the pre-shock flow. Using the ionized absorber model zxipab along with a cooling flow model and a reflection component, we model the broadband Chandra/HETG and NuSTAR spectra of two IPs: NY Lup and V1223 Sgr. We find that this model describes well many of the H and He like emission lines from medium Z elements, which arises from the collisionally excited plasma. However the model fails to account for some of the He like triplets from medium Z elements, which points towards its photo-ionization origin. We do not find a compelling evidence for a blackbody component to model the soft excess seen in the residuals of the Chandra/HETG spectra, which could be due to the uncertainties in estimation of the interstellar absorption of these sources using Chandra/HETG data and/or excess fluxes seen in some photo-ionized emission lines which are not accounted by the cooling flow model. We describe the implications of this model with respect to the geometry of the pre-shock region in these two IPs.

The properties of quasar-host galaxies might be determined by the growth and feedback of their supermassive (SMBH, $10^{8-10}$ M$_{\odot}$) black holes. We investigate such connection with a suite of cosmological simulations of massive (halo mass $\approx 10^{12}$ M$_{\odot}$) galaxies at $z\simeq 6$ which include a detailed sub-grid multiphase gas and accretion model. BH seeds of initial mass $10^5$ M$_{\odot}$ grow mostly by gas accretion, and become SMBH by $z=6$ setting on the observed $M_{\rm BH} - M_{\star}$ relation without the need for a boost factor. Although quasar feedback crucially controls the SMBH growth, its impact on the properties of the host galaxy at $z=6$ is negligible. In our model, quasar activity can both quench (via gas heating) or enhance (by ISM over-pressurization) star formation. However, we find that the star formation history is insensitive to such modulation as it is largely dominated, at least at $z>6$, by cold gas accretion from the environment that cannot be hindered by the quasar energy deposition. Although quasar-driven outflows can achieve velocities $> 1000~\rm km~s^{-1}$, only $\approx 4$% of the outflowing gas mass can actually escape from the host galaxy. These findings are only loosely constrained by available data, but can guide observational campaigns searching for signatures of quasar feedback in early galaxies.

Structure formation in our Universe creates non-Gaussian random fields that will soon be observed over almost the entire sky by the Euclid satellite, the Vera-Rubin observatory, and the Square Kilometre Array. An unsolved problem is how to analyze such non-Gaussian fields best, e.g. to infer the physical laws that created them. This problem could be solved if a parametric non-Gaussian sampling distribution for such fields were known, as this distribution could serve as likelihood during inference. We therefore create a sampling distribution for non-Gaussian random fields. Our approach is capable of handling strong non-Gaussianity, while perturbative approaches such as the Edgeworth expansion cannot. To imitate cosmological structure formation, we enforce our fields to be (i) statistically isotropic, (ii) statistically homogeneous, and (iii) statistically independent at large distances. We generate such fields via a Monte Carlo Markov Chain technique and find that even strong non-Gaussianity is not necessarily visible to the human eye. We also find that sampled marginals for pixel pairs have an almost generic Gauss-like appearance, even if the joint distribution of all pixels is markedly non-Gaussian. This apparent Gaussianity is a consequence of the high dimensionality of random fields. We conclude that vast amounts of non-Gaussian information can be hidden in random fields that appear nearly Gaussian in simple tests, and that it would be short-sighted not to try and extract it.

We present a comparative study of X-ray and IR selected AGNs at $z\approx2$ to highlight the importance of the AGN selection effects on the distributions of star formation (SF) and morphological properties of the host galaxies. We find that while the median SF of X-ray AGN hosts are similar to non-AGN star forming galaxies (SFGs), the incidence of X-ray AGNs, q$_{\rm{AGN}}$, is higher among galaxies with suppressed SF and larger stellar mass surface density within both the half-light radius ($\Sigma_e$) and the central 1 kpc ($\Sigma_{\rm{1kpc}}$), IR AGN hosts are different. They are less massive, have elevated SF and share similar distributions of colors, $\Sigma_e$ and $\Sigma_{\rm{1kpc}}$ with normal SFGs. Given that $\Sigma_e$ and $\Sigma_{\rm{1kpc}}$ strongly correlate with stellar mass (M$_*$), we introduce $\frac{M_{\rm{1kpc}}}{M_*}$, the fractional mass within central 1 kpc, to quantify galaxy compactness, which is independent on M$_*$. Both AGN populations have similar $\frac{M_{\rm{1kpc}}}{M_*}$ distributions to normal SFGs'. We show that while q$_{\rm{AGN}}$ increases with both $\Sigma_e$ and $\Sigma_{\rm{1kpc}}$, it remains constant with $\frac{M_{\rm{1kpc}}}{M_*}$, indicating that the trend of increasing q$_{\rm{AGN}}$ with $\Sigma$ is driven by M$_*$. While our findings are not in conflict with the scenario of AGN quenching, they do not directly imply it either, because the incidence of AGNs being hosted by transitional galaxies depends crucially on AGN selections. The additional evidence that no clear correlation is observed between SF and AGN bolometric luminosity, regardless of the selection, calls into question the notion that AGNs are the direct cause of quenching in $z\approx2$ massive galaxies.

We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) to study the dark matter content of galaxy groups in the Universe using weak gravitational lensing. As our lens sample, we use galaxy groups that have been spectroscopically selected from the Galaxy Mass and Assembly galaxy survey in approximately 100 sq. degrees of the sky that overlap with the HSC survey. We restrict our analysis to the 1587 groups with at least five group members. We divide these galaxy groups into six bins each of galaxy group luminosity and group member velocity dispersion and measure the coherent tangential ellipticity pattern on background HSC galaxies imprinted by weak gravitational lensing. We measure the weak lensing signal with a signal-to-noise ratio of 55 and 51 for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our galaxy groups binned in the two different observable properties and obtain constraints on the power-law scaling relation between mean halo masses and the two group observable properties. We obtain a 5 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\langle M\rangle = (0.81\pm0.04)\times10^{14}h^{-1}M_\odot$ for $L_{\rm grp}=10^{11.5}h^{-2}L_\odot$, and a power-law index of $\alpha=1.01\pm0.07$. We also obtain a 5-percent constraint on the amplitude of the scaling relation between halo mass and velocity dispersion with $\langle M\rangle=(0.93\pm0.05)\times10^{14}h^{-1}M_\odot$ for $\sigma=500{\,\rm kms}^{-1}$ and a power-law index $\alpha=1.52\pm0.10$, although these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature indicate that our results are consistent, but have significantly reduced errors.

Context: In spiral galaxies, star formation tends to trace features of the spiral pattern, including arms, spurs, feathers, and branches. However, in our own Milky Way, it has been challenging to connect individual star-forming regions to their larger Galactic environment owing to our perspective from within the disk. One feature in nearly all modern models of the Milky Way is the Sagittarius Arm, located inward of the Sun with a pitch angle of ~12 deg. Aims: We map the 3D locations and velocities of star-forming regions in a segment of the Sagittarius Arm using young stellar objects (YSOs) from the Spitzer/IRAC Candidate YSO (SPICY) catalog to compare their distribution to models of the arm. Methods: Distances and velocities for these objects are derived from Gaia EDR3 astrometry and molecular line surveys. We infer parallaxes and proper motions for spatially clustered groups of YSOs and estimate their radial velocities from the velocities of spatially associated molecular clouds. Results: We identify 25 star-forming regions in the Galactic longitude range l~4.0-18.5 deg arranged in a narrow, ~1 kpc long linear structure with a high pitch angle of $\psi = 56$ deg and a high aspect ratio of ~7:1. This structure includes massive star-forming regions such as M8, M16, M17, and M20. The motions in the structure are remarkably coherent, with velocities in the direction of Galactic rotation of $240\pm3$ km/s (slightly higher than average) and slight drifts toward the Galactic center (-4.3 km/s) and in the negative Z direction (-2.9 km/s). The rotational shear experienced by the structure is 4.6 km/s/kpc. Conclusions: The observed 56 deg pitch angle is remarkably high for a segment of the Sagittarius Arm. We discuss possible interpretations of this feature as a substructure within the lower pitch angle Sagittarius Arm, as a spur, or as an isolated structure.

Background: Most of the stars in the Universe will end their evolution by losing their envelope during the thermally pulsing asymptotic giant branch (TP-AGB) phase, enriching the interstellar medium of galaxies with heavy elements, partially condensed into dust grains formed in their extended circumstellar envelopes. Among these stars, carbon-rich TP-AGB stars (C-stars) are particularly relevant for the chemical enrichment of galaxies. We here investigated the role of the metallicity in the dust formation process from a theoretical viewpoint. Methods: We coupled an up-to-date description of dust growth and dust-driven wind, which included the time-averaged effect of shocks, with FRUITY stellar evolutionary tracks. We compared our predictions with observations of C-stars in our Galaxy, in the Magellanic Clouds (LMC and SMC) and in the Galactic Halo, characterised by metallicity between solar and 1/10 of solar. Results: Our models explained the variation of the gas and dust content around C-stars derived from the IRS Spitzer spectra. The wind speed of the C-stars at varying metallicity was well reproduced by our description. We predicted the wind speed at metallicity down to 1/10 of solar in a wide range of mass-loss rates.

We present redshift space two-point ($\xi$), three-point ($\zeta$) and reduced three-point (Q) correlation of Ly$\alpha$ absorbers (i.e Voigt profile components having HI column density $N_{HI}>10^{13.5}$cm$^{-2}$) over three redshift bins spanning $1.7<z<3.5$ using high resolution spectra of 292 quasars. We detect positive $\xi$ up to 8$h^{-1}$cMpc in all three redshift bins. The strongest detection of $\zeta$ is seen in $z=1.7-2.3$ redshift bin at $1-2h^{-1}$cMpc with an amplitude of $1.81\pm0.59$ ($\sim3.1\sigma$ level). The corresponding Q is found to be $0.68\pm0.23$. The measured $\xi$ and $\zeta$ values show an increasing trend with $N_{HI}$, while Q remains relatively independent of $N_{HI}$. We find $\xi$ and $\zeta$ to evolve strongly with $z$ over the redshift range studied. Using hydrodynamical simulations, we find that the $\xi$ and especially $\zeta$ seen in real space may be strongly amplified by peculiar velocities in redshift space. Simulations also suggest that while feedback, thermal and pressure smoothing effects affect the clustering of Ly$\alpha$ absorbers at small scales, i.e $<0.5h^{-1}$cMpc, the HI photo-ionization rates ($\Gamma_{HI}$) strongly influence the correlation amplitudes at all scales. We find that the strong redshift evolution shown by $\xi$ and $\zeta$ is primarily sourced by the redshift evolution of the relationship between $N_{HI}$ and baryon overdensity ($\Delta$). Our simulations that uses available best fitted $\Gamma_{HI}(z)$ measurements produce consistent clustering signals with observations at $z\sim2$ but under-predict the clustering at higher redshifts. One possible remedy is to have higher values of $\Gamma_{HI}$ at higher redshifts compared to the existing measurements. Alternatively the discrepancy could be related to non-equilibrium and inhomogeneous conditions prevailing during HeII reionization not captured by our simulations.

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $\eta$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($\eta$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature \Tx, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged)

We investigate the origin of the abundance ratios and scatter of $\alpha$- and neutron-capture elements of old, metal-poor stars, using cosmological, hydrodynamical simulations of galaxy formation. For this, we implement a novel treatment for the production and distribution of chemical products of Type II supernovae, which considers the effects of the rotation of massive stars on the chemical yields and the effects of the different life-times of stars that are progenitors of this type of supernovae. We focus on the stellar halo of a Milky Way-mass galaxy, studying the abundances and scatter of [O/Fe], [Mg/Fe], [Si/Fe], [Sr/Fe], [Eu/Fe] and [Ba/Fe]. Our model is able, for the first time in a cosmological simulation, to describe at the same time the low scatter in the abundances of $\alpha$-elements and the higher scatter associated to neutron-capture elements in the halo stars, as suggested by observations of the Milky Way. We also reproduce the scatter observed in the [Sr/Ba] ratio, which results from the treatment of the fast-rotating stars and the dependence of the chemical yields on the metallicity, mass and rotational velocities. Our simulations show that such scatter patterns appear naturally if the different ejection times associated to stars of different mass are properly described, without the need to invoke for additional mixing mechanisms or a distinct treatment of the alpha- and neutron-capture elements. Simulations of this type will help characterizing and identifying the past accretion debris as well as the pristine in-situ populations in the Galaxy unveiled by Gaia and spectroscopic data.

We comprehensively study the effects of bubble wall thickness and speed on the gravitational wave emission spectrum of collisions of two vacuum bubbles. We numerically simulate a large dynamical range, making use of symmetry to reduce the dimensionality. The high-frequency slope of the gravitational wave spectrum is shown to depend on the thickness of the bubble wall, becoming steeper for thick-wall bubbles, in agreement with recent fully 3+1 dimensional lattice simulations of many-bubble collisions. This dependence is present, even for highly relativistic bubble wall collisions. We use the reduced dimensionality as an opportunity to investigate dynamical phenomena which may underlie the observed differences in the gravitational wave spectra. These phenomena include `trapping', which occurs most for thin-wall bubbles, and oscillations behind the bubble wall, which occur for thick-wall bubbles.

We present an analysis of a densely repeating sample of bursts from the first repeating fast radio burst, FRB 121102. We detected a total of 133 bursts in 3 hours of data at a center frequency of 1.4 GHz using the Arecibo Telescope, and develop robust modeling strategies to constrain the spectro-temporal properties of all the bursts in the sample. Most of the burst profiles show a scattering tail, and burst spectra are well modeled by a Gaussian with a median width of 230 MHz. We find a lack of emission below 1300 MHz, consistent with previous studies of FRB 121102. We also find that the peak of the log-normal distribution of wait times decreases from 207 s to 75 s using our larger sample of bursts. Our observations do not favor either Poissonian or Weibull distributions for the burst rate distribution. We searched for periodicity in the bursts using multiple techniques, but did not detect any significant period. The cumulative burst energy distribution exhibits a broken power-law shape, with the lower and higher-energy slopes of $-0.4\pm0.1$ and $-1.8\pm0.2$, with the break at $(2.3\pm0.2)\times 10^{37}$ ergs. We provide our burst fitting routines as a python package \textsc{burstfit}. All the other analysis scripts and results are publicly available.

The CHIME/FRB experiment has detected thousands of Fast Radio Bursts (FRBs) due to its sensitivity and wide field of view; however, its low angular resolution prevents it from localizing events to their host galaxies. Very Long Baseline Interferometry (VLBI), triggered by FRB detections from CHIME/FRB will solve the challenge of localization for non-repeating events. Using a refurbished 10-m radio dish at the Algonquin Radio Observatory located in Ontario Canada, we developed a testbed for a VLBI experiment with a theoretical ~<30 masec precision. We provide an overview of the 10-m system and describe its refurbishment, the data acquisition, and a procedure for fringe fitting that simultaneously estimates the geometric delay used for localization and the dispersive delay from the ionosphere. Using single pulses from the Crab pulsar, we validate the system and localization procedure, and analyze the clock stability between sites, which is critical for phase-referencing an FRB event. We find a localization of 50 masec is possible with the performance of the current system. Furthermore, for sources with insufficient signal or restricted wideband to simultaneously measure both geometric and ionospheric delays, we show that the differential ionospheric contribution between the two sites must be measured to a precision of 1e-8 pc/cc to provide a reasonable localization from a detection in the 400--800 MHz band. Finally we show detection of an FRB observed simultaneously in the CHIME and the Algonquin 10-m telescope, the first FRB cross-correlated in this very long baseline. This project serves as a testbed for the forthcoming CHIME/FRB Outriggers project.

During a period of strong $\gamma$-ray flaring activity from BL Lacertae, we organized Swift, NICER, and NuSTAR follow-up observations. The source has been monitored by Swift-XRT between 2020 August 11 and October 16, showing a variability amplitude of 65, with a flux varying between 1.0 $\times$ 10$^{-11}$ and 65.3 $\times$ 10$^{-11}$ erg cm$^{-2}$ s$^{-1}$. On 2020 October 6, Swift-XRT has observed the source during its historical maximum X-ray flux. A softer-when-brighter behaviour has been observed by XRT, suggesting an increasing importance of the synchrotron emission in the X-ray part of the spectrum covered by XRT during bright states. Rapid variability in soft X-rays has been observed with both the Swift-XRT and NICER observations with a minimum variability time-scale of 60 s and 240 s, respectively, suggesting very compact emitting regions (2.4 $\times$ 10$^{13}$ cm and 9.5 $\times$ 10$^{13}$ cm). At hard X-rays, a minimum variability time-scale of $\sim$ 5.5 ks has been observed by NuSTAR. We report the first simultaneous NICER and NuSTAR observations of BL Lacertae during 2020 October 11-12. The joint NICER and NuSTAR spectra are well fitted by a broken power-law with a significant difference of the photon index below (2.10) and above (1.60) an energy break at $\sim$ 2.7 keV, indicating the presence of two different emission components (i.e, synchrotron and inverse Compton) in the broad band X-ray spectrum. Leaving the total hydrogen column density toward BL Lacertae free to vary, a value of N$_{H,tot}$ = (2.58 $\pm$ 0.09) $\times$ 10$^{21}$ cm$^{-2}$ has been estimated.

(Abridged) Dynamical friction can be used to distinguish Newtonian gravity and modified Newtonian dynamics (MOND) because it works differently in these frameworks. This concept, however, has yet to be explored very much with MOND. Previous simulations showed weaker dynamical friction during major mergers for MOND than for Newtonian gravity with dark matter. Analytic arguments suggest the opposite for minor mergers. In this work, we verify the analytic predictions for MOND by high-resolution $N$-body simulations of globular clusters (GCs) moving in isolated ultra-diffuse galaxies (UDGs). We test the MOND analog of the Chandrasekhar formula for the dynamical friction proposed by S\'anchez-Salcedo on a single GC. We also explore whether MOND allows GC systems of isolated UDGs to survive without sinking into nuclear star clusters. The simulations are run using the adaptive-mesh-refinement code Phantom of Ramses. The mass resolution is $20\,M_\odot$ and the spatial resolution $50\,$pc. The GCs are modeled as point masses. Simulations including a single GC reveal that, as long as the apocenter of the GC is over about 0.5 effective radii, the S\'anchez-Salcedo formula works excellently, with an effective Coulomb logarithm increasing with orbital circularity. Once the GC reaches the central kiloparsec, its sinking virtually stops, likely because of the core stalling mechanism. In simulations with multiple GCs, many of them sink toward the center, but the core stalling effect seems to prevent them from forming a nuclear star cluster. The GC system ends up with a lower velocity dispersion than the stars of the galaxy. By scaling the simulations, we extend these results to most UDG parameters, as long as these UDGs are not external-field dominated.

(Abridged) NGC 985 was observed by XMM-Newton twice in 2015, revealing that the source was coming out from a soft X-ray obscuration event that took place in 2013. These kinds of events are possibly recurrent since a previous XMM-Newton archival observation in 2003 also showed signatures of partial obscuration. We have analyzed the high-resolution X-ray spectra of NGC 985 obtained by the RGS in 2003, 2013, and 2015 in order to characterize the ionized absorbers superimposed to the continuum and to study their response as the ionizing flux varies. We found that up to four warm absorber (WA) components were present in the grating spectra of NGC 985, plus a mildy ionized (log xi ranging between 0.2 and 0.5) obscuring (log N(H) of about 22.3) wind outflowing at about 6000 km/s. The absorbers have a log N(H) ranging from 21 to about 22.5, and ionization parameters ranging from 1.6 to 2.9. The most ionized component is also the fastest, moving away at 5100 km/s, while the others outflow in two kinematic regimes, at about 600 and 350 km/s. These components showed variability at different time scales in response to changes in the ionizing continuum. Assuming that these changes are due to photoionization we have obtained upper and lower limits on the density of the gas and therefore on its distance, finding that the closest two components are at pc-scale distances, while the rest may extend up to tens of pc from the central source. The fastest, most ionized WA component accounts for the bulk of the kinetic luminosity injected back into the ISM of the host galaxy, which is on the order of 0.8% of the bolometric luminosity of NGC 985. According to the models, this amount of kinetic energy per unit time would be sufficient to account for cosmic feedback.

Young massive stars warm up the large amount of gas and dust which condenses in their vicinity, exciting a forest of lines from different molecular species. Their line brightness is a diagnostic tool of the gas physical conditions locally, which we use to set constraints on the environment where massive stars form. We made use of the Atacama Large Millimeter/submillimeter Array at frequencies near 349 GHz, with an angular resolution of $0.1''$, to observe the methyl cyanide (CH$_3$CN) emission which arises from the accretion disk of a young massive star. We sample the disk midplane with twelve distinct beams, where we get an independent measure of the gas (and dust) physical conditions. The accretion disk extends above the midplane showing a double-armed spiral morphology projected onto the plane of the sky, which we sample with ten additional beams: along these apparent spiral features, gas undergoes velocity gradients of about $\rm 1 km s^{-1}$ per 2000 au. The gas temperature (T) rises symmetrically along each side of the disk, from about 98 K at 3000 au to 289 K at 250 au, following a power law with radius, R$^{-0.43}$. The CH$_3$CN column density (N) increases from $\rm 9.2\times10^{15} cm^{-2}$ to $\rm 8.7\times10^{17} cm^{-2}$ at the same radii, following a power law with radius, R$^{-1.8}$. In the framework of a circular gaseous disk observed approximately edge-on, we infer an H$_2$ volume density in excess of $\rm 4.8\times10^9 cm^{-3}$ at a distance of 250 au from the star. We study the disk stability against fragmentation following the methodology by Kratter et al. (2010), appropriate under rapid accretion, and we show that the disk is marginally prone to fragmentation along its whole extent.

We report the detection of three RR Lyrae (RRL) stars (two RRc and one RRab) in the ultra-faint dwarf (UFD) galaxy Centaurus I (Cen I) and two Milky Way (MW) $\delta$ Scuti/SX Phoenicis stars based on multi-epoch $giz$ DECam observations. The two RRc stars are located within 2 times the half-light radius (r$_h$) of Cen I, while the RRab star (CenI-V3) is at $\sim6$ r$_h$. The presence of three distant RRL stars clustered this tightly in space represents a 4.7$\sigma$ excess relative to the smooth distribution of RRL in the Galactic halo. Using the newly detected RRL stars, we obtain a distance modulus to Cen I of $\mu_0 = 20.35 \pm 0.05$ mag (D$_{\odot} = 118 \pm 3$ kpc). The location of the RRL stars in the Bailey diagram allows us to classify Cen I as an Oosterhoff II system, in agreement with other UFD galaxies. Finally, we study the relative rate of RRc+RRd (RRcd) stars ($f_{cd}$) in UFD and classical dwarf galaxies. The full sample of MW dwarf galaxies gives a mean of $f_{cd} = 0.28$. While several UFD galaxies, such as Cen I, present higher RRcd ratios, if we combine the RRL populations of all UFD galaxies, the RRcd ratio is similar to the one obtained for the classical dwarfs ($f_{cd}$ $\sim$ 0.3). Therefore, there is no evidence for a different fraction of RRcd stars in UFD and classical dwarf galaxies.

Both simulation and observational data have shown that the spin and shape of dark matter halos are correlated with their nearby large-scale environment. As structure formation on different scales is strongly coupled, it is trick to disentangle the formation of halo with the large-scale environment, making it difficult to infer which is the driving force for the correlation between halo spin/shape with the large-scale structure. In this paper, we use N-body simulation to produce twin Universes that share the same initial conditions on small scales but different on large scales. This is achieved by changing the random seeds for the phase of those k modes smaller than a given scale in the initial conditions. In this way, we are able to disentangle the formation of halo and large-scale structure, making it possible to investigate how halo spin and shape correspond to the change of environment on large scales. We identify matching halo pairs in the twin simulations as those sharing the maximum number of identical particles within each other. Using these matched halo pairs, we study the cross match of halo spin and their correlation with the large-scale structure. It is found that when the large-scale environment changes (eigenvector) between the twin simulations, the halo spin has to rotate accordingly, although not significantly, to maintain the universal correlation seen in each simulation. Our results suggest that the large-scale structure is the main factor to drive the correlation between halo properties and their environment.

Rapid binary population synthesis codes are often used to investigate the evolution of compact-object binaries. They typically rely on analytical fits of single-star evolutionary tracks and parameterized models for interactive phases of evolution (e.g., mass-transfer on thermal timescale, determination of dynamical instability, and common envelope) that are crucial to predict the fate of binaries. These processes can be more carefully implemented in stellar structure and evolution codes such as MESA. To assess the impact of such improvements, we compare binary black hole mergers as predicted in models with the rapid binary population synthesis code COSMIC to models ran with MESA simulations through mass transfer and common-envelope treatment. We find that results significantly differ in terms of formation paths, the orbital periods and mass ratios of merging binary black holes, and consequently merger rates. While common-envelope evolution is the dominant formation channel in COSMIC, stable mass transfer dominates in our MESA models. Depending upon the black hole donor mass, and mass-transfer and common-envelope physics, at sub-solar metallicity COSMIC overproduces the number of binary black hole merges by factors of 2--35 with a significant fraction of them having merger times orders of magnitude shorter than the binary black holes formed when using detailed MESA models. Therefore we find that some binary black hole merger rate predictions from rapid population syntheses of isolated binaries may be overestimated by factors of ~5--500. We conclude that the interpretation of gravitational-wave observations requires the use of detailed treatment of these interactive binary phases.

The PolyOculus technology produces large-area-equivalent telescopes by using fiber optics to link modules of multiple semi-autonomous, small, inexpensive, commercial-off-the-shelf telescopes. Crucially, this scalable design has construction costs which are > 10x lower than equivalent traditional large-area telescopes. We have developed a novel photonic lantern approach for the PolyOculus fiber optic linkages which potentially offers substantial advantages over previously considered free-space optical linkages, including much higher coupling efficiencies. We have carried out the first laboratory tests of a photonic lantern prototype developed for PolyOculus, and demonstrate broadband efficiencies of ~91%, confirming the outstanding performance of this technology.

We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through the sedimentation factor $f_\mathrm{sed}$ plays a larger role in shaping cloud thickness and its radiative feedback on the local gas temperatures -- e.g. hot spot shift and day-to-night side temperature gradient -- than the switch in mixing treatment. We demonstrate using our new mixing treatments that simulations with cloud scales which are a fraction of the pressure scale height improve agreement with the observed transmission spectra, the emission spectra, and the Spitzer 4.5 $\mathrm{\mu m}$ phase curve, although our models are still unable to reproduce the optical and UV transmission spectra. We also find that the inclusion of cloud increases the transit asymmetry in the optical between the east and west limbs, although the difference remains small ($\lesssim 1\%$).

We report the results of ongoing monitoring of the 6.7 GHz CH$_3$OH masers associated with G188.95+0.89. In these observations five features are periodically varying and at least two exhibit evidence of velocity drifts. It is not clear the cause of these velocity drifts. The spectra have varied significantly since detection in 1991. The 11.45 km s$^1$ feature has decreased exponentially from 2003. Complementary ALMA 1.3 mm continuum and line observational results are also presented. Eight continuum cores (MM1 - MM8) were detected in G188.95+0.89. We derived the masses of the detected cores. G188.95+0.89 MM2 was resolved into 2 continuum cores (separated by 0.1 arcsec) in ALMA band 7 observations. Also CH$_3$OH (4$_{(2,2)}$-3$_{(1,2)})$ thermal emission associated with MM2 is double peaked. We propose the presence of multiple (at least binary) young stellar objects in MM2. SiO emission exhibit a bow-shock morphology in MM2 while strong emission of $^{12}$CO at the east and west of MM2 suggest the presence of an east-west bipolar outflow.

We analyze surveys of molecular cloud structures defined by tracers ranging from CO $J = 1-0$ through $^{13}$CO $J = 1-0$ to dust emission together with NH$_3$ data. The mean value of the virial parameter and the fraction of mass in bound structures depends on the method used to identify structures. Generally, the virial parameter decreases and the fraction of mass in bound structures increases with the effective density of the tracer, the surface density and mass of the structures, and the distance from the center of a galaxy. For the most complete surveys of structures in the Galaxy defined by CO $J = 1-0$, the fraction of mass that is in bound structures is 0.19. For catalogs of other galaxies based on CO $J = 2-1$, the fraction is 0.35. These results offer substantial alleviation of the fundamental problem of slow star formation. If only clouds found to be bound are counted and they are assumed to collapse in a free-fall time at their mean cloud density, the sum over all clouds in a complete survey of the Galaxy yields a predicted star formation rate of 46 solar masses per year, a factor of 6.5 less than if all clouds are bound.

We report high-cadence photometry of the ultra-fast ($t_2\sim1.2$ d) nova V1674 Her during its rise to maximum light ($V\sim6.3$) and the beginning of its subsequent decline. These observations from Evryscope and the Mount Laguna Observatory All-Sky Camera reveal a plateau in the pre-maximum light curve at $g\sim14$ ($\sim$8 mag below peak) that lasted for at least three hours. Similar features (so-called pre-maximum halts) have been observed in some novae near maximum light, but to our knowledge the detection of a plateau in the light curve $\sim$8 mag below peak is unprecedented.

Supernovae Ia (SN) are among the brightest objects we can observe and can provide a unique window on the large scale structure of the Universe at redshifts where other observations are not available. The photons emitted by SNe are in fact affected by the density field between the source and the observer, and from the observed luminosity distance it is possible to solve the inversion problem (IP), i.e. to reconstruct the density field which produced those effects. So far the IP was only solved assuming some restrictions about the geometry of the problem, such as spherical symmetry for example, and the approach was based on solving complicated systems of differential equations which required smooth function as inputs, while observational data is not smooth, due to its discrete nature. In order to overcome these limitations we develop for the first time an inversion method which is not assuming any symmetry, and can be applied directly to observational data, without the need of any data smoothing procedure. The method is based on the use of convolutional neural networks (CNN) trained on simulated data, and it shows quite accurate results. The training data set is obtained by first generating random density and velocity profiles, and then computing their effects on the luminosity distance. The CNN is then trained to reconstruct the density field from the luminosity distance. The CNN is a modified version of U-Net to account for the tridimensionality of the data, and can reconstruct the density and velocity fields with a good level of accuracy. The use of neural networks to analyze observational data from future SNe catalogues will allow to reconstruct the large scale structure of the Universe to an unprecedented level of accuracy, at a redshift at which few other observations are available.

Magnetohydrodynamical evolution of the Richtmyer-Meshkov instability (RMI) is investigated by two-dimensional MHD simulations. The RMI is suppressed by a strong magnetic field, whereas the RMI amplifies an ambient magnetic field by many orders of magnitude if the seed field is weak. We have found that the suppression and amplification processes can be evaluated continuously along with the amplitude of the Alfv\'en number $R_A$, which is defined as the ratio of the linear growth velocity of the RMI to the Alfv\'en speed at the interface. When the Alfv\'en number is less than unity, the Lorentz force acting on the fluid mitigates the unstable motion of the RMI significantly, and the interface oscillates stably in this limit. If $R_A \gtrsim 1$, on the other hand, the surface modulation increases due to the growth of the RMI. The maximum strength of the magnetic field is enhanced up to by a factor of $R_A$. This critical feature is universal and independent of the initial Mach number of the incident shock, the Atwood number, corrugation amplitude, and even the direction of the initial magnetic field.

At present, 19 double neutron star (DNS) systems are detected by radio timing and 2 merging DNS systems are detected by kilo-hertz gravitational waves. Because of selection effects, none of them has an orbital period $P_b$ in the range of a few tens of minutes. In this paper we consider a multimessenger strategy proposed by Kyutoku et al. (2019), jointly using the Laser Interferometer Space Antenna (LISA) and the Square Kilometre Array (SKA) to detect and study Galactic pulsar-neutron star (PSR-NS) systems with $P_b \sim$ 10-100 min. We assume that we will detect PSR-NS systems by this strategy. We use standard pulsar timing software to simulate times of arrival of pulse signals from these binary pulsars. We obtain the precision of timing parameters of short-orbital-period PSR-NS systems whose orbital period $P_b \in (8,120)\,$min. We use the simulated uncertainty of the orbital decay, $\dot{P}_{b}$, to predict future tests for a variety of alternative theories of gravity. We show quantitatively that highly relativistic PSR-NS systems will significantly improve the constraint on parameters of specific gravity theories in the strong field regime. We also investigate the orbital periastron advance caused by the Lense-Thirring effect in a PSR-NS system with $P_b = 8\,$min, and show that the Lense-Thirring effect will be detectable to a good precision.

During the last decade in star formation research, many studies have targeted low- and high-mass star formation regions located at different distances, with different telescopes having specific angular resolution capabilities. We present a systematic investigation of the angular resolution effects, with special attention being paid to the derived masses of sources as well as the shape of the resulting source mass functions (SMFs). We tested the impact of angular resolution, from 0.6 down to 0.02 pc, in two star-forming regions observed with Herschel (NGC6334 and Aquila), and three (magneto)-hydrodynamical simulations. We detected and measured sources at each resolution using getsf and we analysed the derived masses and sizes of the sources. We find that the number of sources does not converge from 0.6 to 0.05 pc. It increases by about two when the angular resolution increases with a similar factor. Below 0.05 pc, the number of source still increases by about 1.3 when the angular resolution increases by two, suggesting that we are close to, but not yet at, convergence. We find that the measured sizes and masses of sources linearly depend on the angular resolution with no sign of convergence to a resolution-independent value. The corresponding SMF peak also shifts with angular resolution, while the slope of the high-mass tail of the SMFs remains almost invariant. If prestellar cores, physically distinct from their background, exist in cluster-forming molecular clouds, we conclude that their mass must be lower than reported so far in the literature. We discuss various implications for the studies of star formation: the problem of determining the mass reservoirs involved in the star-formation process; the inapplicability of the Gaussian beam deconvolution to infer source sizes; and the impossibility to determine the efficiency of the mass conversion from the cores to the stars.

We investigate the escape process from a perpendicular shock region of a spherical shock in the interstellar medium (ISM). The diffusive shock acceleration in the perpendicular shock of supernova remnants (SNRs) has been expected to accelerate cosmic rays (CRs) to the PeV scale without an upstream magnetic field amplification. We estimate the maximum energy of CRs limited by the escape from the perpendicular shock region. By performing test particle simulations, we confirm the theoretical estimation, showing that the escape-limited maximum energy in the perpendicular shock is several 10 TeV for the typical type Ia SNRs. Therefore, in order for SNRs in the ISM to accelerate CRs to the PeV scale, an upstream magnetic field amplification is needed. The characteristic energy scale of several 10 TeV could be the origin of the spectral break around 10 TeV, which was reported by recent direct CR observations. In addition, we show that, in the free expansion phase, the rapid perpendicular shock acceleration works on about 20% area of the whole shock surface, which is larger than the size of the superluminal shock region. We also discuss the escape of CR electrons from the perpendicular shock.

The recently launched TESS mission is for the first time giving us the potential to perform inference asteroseismology across the whole sky. TESS observed the Kepler field entirely in its Sector 14 and partly in Sector 15. Here, we seek to detect oscillations in the red giants observed by TESS in the Kepler field of view. Using the full 4-yr Kepler results as the ground truth, we aim to characterise how well the seismic signal can be detected using TESS data. Because our data are based on one and two sectors of observation, our results will be representative of what one can expect for the vast majority of the TESS data. We detect clear oscillations in $\sim$3000 stars with another $\sim$1000 borderline (low S/N) cases, all of which yield a measurement of the frequency of maximum acoustic power, numax. In comparison, a simple calculation predicts $\sim$4500 stars would show detectable oscillations. Of the clear detections we reliably measure the frequency separation between overtone radial modes, dnu, in 570 stars, meaning an overall dnu yield of 20%, which splits into a one-sector yield of 14% and a two-sector yield of 26%. These yields imply that typical (1-2 sector) TESS data will result in significant detection biases. Hence, to boost the number of stars, one might need to use only numax as the seismic input for stellar property estimation. On the up side, we find little or no bias in the seismic measurements and typical scatter relative to the Kepler `truth' is about 5-6% in numax and 2-3% in dnu. These values, coupled with typical uncertainties in parallax, Teff, and Fe/H in a grid-based approach, would provide internal uncertainties of 3% in inferred stellar radius, 6% in mass and 20% in age. Finally, despite relatively large pixels of TESS, we find red giant seismology is not expected to be significantly affected by blending for stars with Tmag < 12.5.

We carry out a test of the radial acceleration relation (RAR) for a sample of 10 dynamically relaxed and cool-core galaxy clusters imaged by the Chandra X-ray telescope, which was studied in Giles et al. For this sample, we observe that the best-fit RAR shows a very tight residual scatter equal to 0.09 dex. We obtain an acceleration scale of $1.59 \times 10^{-9} m/s^2$, which is about an order of magnitude higher than that obtained for galaxies. Furthermore, the best-fit RAR parameters differ from those estimated from some of the previously analyzed cluster samples, which indicates that the acceleration scale found from the RAR could be of an emergent nature, instead of a fundamental universal scale.

Many supermassive black holes (SMBH) of mass $10^{6\sim9}M_{\odot}$ are observed at the center of each galaxy even in the high redshift ($z\approx7$) Universe. To explain the early formation and the common existence of SMBH, we proposed previously the SMBH formation scenario by the gravitational collapse of the coherent dark matter (DM) composed from the Bose-Einstein Condensed (BEC) objects. A difficult problem in this scenario is the inevitable angular momentum which prevents the collapse of BEC. To overcome this difficulty, in this paper, we consider the very early Universe when the BEC-DM acquires its proper angular momentum by the tidal torque mechanism. The balance of the density evolution and the acquisition of the angular momentum determines the mass of the SMBH as well as the mass ratio of BH and the surrounding dark halo (DH). This ratio turns out to be $M_{BH}/M_{DH}\approx10^{-3\sim-5}(M_{tot}/10^{12}\mathrm{M}_{\odot})^{-1/2}$ assuming simple density profiles of the initial DM cloud. This estimate turns out to be consistent with the observations at $z\approx0$ and $z\approx6$, although the data scatter is large. Thus the angular momentum determines the separation of black and dark, \textsl{i.e. }SMBH and DH, in the original DM cloud.

Volatile elements and compounds found in extra-terrestrial environments are often the target of In Situ Resource Utilization (ISRU) studies. Although water and hydroxide are most commonly the focus of these studies as they can be used for propellant and human consumption; we instead focus on the possible exploitation of sulfur and how it could be utilized to produce building materials on the Moon, Mars and Asteroids. We describe the physical and chemical pathways for extracting sulfur from native sulfide minerals, manufacturing sulfuric acid in situ, and using the produced acid to dissolve native silicate minerals. The final products of this process, which we call the Silicate-Sulfuric Acid Process (SSAP), include iron metal, silica, oxygen and metal oxides, all of which are crucial in the scope of a sustainable, space-based economy. Although our proposed methodology requires an initial investment of water, oxygen, and carbon dioxide, we show that all of these volatiles are recovered and reused in order to repeat the process. We calculate the product yield from this process if it were enacted on the lunar highlands, lunar mare, Mars, as well as an array of asteroid types.

We present a first catalog of sources detected by the Mikhail Pavlinsky ART-XC telescope aboard the SRG observatory in the 4-12 keV energy band during its on-going all-sky survey. The catalog comprises 867 sources detected on the combined map of the first two 6-month scans of the sky (Dec. 2019 - Dec. 2020) - ART-XC sky surveys 1 and 2, or ARTSS12. The achieved sensitivity to point sources varies between ~5x10-12 erg/s/cm2 near the Ecliptic plane and better than 10-12 erg/s/cm2 (4-12 keV) near the Ecliptic poles, and the typical localization accuracy is ~15 arcsec. Among the 750 sources of known or suspected origin in the catalog, 56% are extragalactic (mostly active galactic nuclei (AGN) and clusters of galaxies) and the rest are Galactic (mostly cataclysmic variables (CVs) and low- and high-mass X-ray binaries). For 116 sources ART-XC has detected X-rays for the first time. Although the majority of these (~80) are expected to be spurious (for the adopted detection threshold), there can be a significant number of newly discovered astrophysical objects. We have started a program of optical follow-up observations of the new and previously unidentified X-ray sources, which has already led to the identification of several AGN and CVs. With the SRG all-sky survey planned to continue for a total of 4 years, we can expect the ART-XC survey in the 4-12 keV band to significantly surpass the previous surveys carried out in similar (medium X-ray) energy bands in terms of the combination of angular resolution, sensitivity, and sky coverage.

We report the first unambiguous observational evidence in the radio range of the reflection of a coronal shock wave at the boundary of a coronal hole. The event occurred above an active region located at the northwest limb of the Sun and was characterized by an eruptive prominence and an extreme-ultraviolet (EUV) wave steepening into a shock. The EUV observations acquired by the Atmospheric Imaging Assembly (AIA) instrument on board the Solar Dynamics Observatory(SDO) and the Extreme Ultraviolet Imager (EUVI) instrument on board the Solar TErrestrial RElations Observatory(STEREO-A) were used to track the development of the EUV front in the inner corona. Metric type II radio emission, a distinguishing feature of shock waves propagating in the inner corona, was simultaneously recorded by ground-based radio spectrometers. The radio dynamic spectra displayed an unusual reversal of the type II emission lanes, together with type III-like herringbone emission, indicating shock-accelerated electron beams. Combined analysis of imaging data from the two space-based EUV instruments and the Nancay Radioheliograph (NRH) evidences that the reverse-drifting typeiiemission was produced at the intersection of the shock front, reflected at a coronal hole boundary, with an intervening low-Alfv\'en-speed region characterized by an open field configuration. We also provide an outstanding data-driven reconstruction of the spatiotemporal evolution in the inner corona of the shock-accelerated electron beams produced by the reflected shock.

FeI and NiI emission lines have recently been found in the spectra of 17 Solar System comets observed at heliocentric distances between 0.68 and 3.25 au and in the interstellar comet 2I/Borisov. The blackbody equilibrium temperature at the nucleus surface is too low to vaporize the refractory dust grains that contain metals, making the presence of iron and nickel atoms in cometary atmospheres a puzzling observation. Moreover, the measured NiI/FeI abundance ratio is on average one order of magnitude larger than the solar photosphere value. We report new measurements of FeI and NiI production rates and abundance ratios for the Jupiter-family comet (JFC) 46P/Wirtanen in its 2018 apparition and from archival data of the Oort-cloud comet (OCC) C/1996 B2 (Hyakutake). The comets were at geocentric distances of 0.09 au and 0.11 au, respectively. The emission line surface brightness was found to be inversely proportional to the projected distance to the nucleus, confirming that FeI and NiI atoms are ejected from the surface of the nucleus or originate from a short-lived parent. Considering the full sample of 20 comets, we find that the range of NiI/FeI abundance ratios is significantly larger in JFCs than in OCCs. We also unveil significant correlations between NiI/FeI and C$_2$/CN, C$_2$H$_6$/H$_2$O, and NH/CN. Carbon-chain- and NH-depleted comets show the highest NiI/FeI ratios. The existence of such relations suggests that the diversity of NiI/FeI abundance ratios in comets could be related to the cometary formation rather than to subsequent processes~in~the~coma.

Symbiotic stars belong to a group of interacting binaries that display a wide variety of phenomena, including prominent outbursts connected with mass transfer, as well as stellar winds, jets, eclipses, or intrinsic variability of the components. Dozens of new symbiotic stars and candidates have been discovered in recent years. However, there are many objects which are still poorly studied. Some symbiotic candidates suspected in the literature have never been studied spectroscopically. In this contribution, we present the first results of the ongoing campaign focused on symbiotic candidates. In the first paper in the series, we study the nature of ten candidate classical symbiotic stars suspected based on their photometric behaviour, colours or abundance pattern. To confirm or reject the symbiotic nature of the studied candidates, we obtained new spectra and analysed them in detail together with available multi-frequency photometric and spectroscopic observations of the objects. Hen 3-860 and V2204 Oph are genuine symbiotic systems showing typical spectral features of burning symbiotic stars and outbursts in the last 100 years. The first object belongs to the uncommon group of eclipsing symbiotic stars. V1988 Sgr cannot be classified as a genuine burning symbiotic star, but the scenario of an accreting-only symbiotic system cannot be ruled out. Hen 4-204 might be a bona-fide symbiotic star due to its similarity with the known symbiotic binary BD Cam. Six other symbiotic candidates (V562 Lyr, IRAS 19050+0001, EC 19249-7343, V1017 Cyg, PN K1-6, V379 Peg) are either single dwarf or giant stars or non-symbiotic binaries.

We study the density structures of dark matter subhalos for both cold dark matter and self-interacting dark matter models using high-resolution cosmological $N$-body simulations. We quantify subhalo's central density at 150 pc from the center of each subhalo at the classical dwarf spheroidal and ultrafaint dwarf scales. By comparing them with observations, we find that the self-interacting scattering cross-section of $\sigma/m<3\ \rm{cm^{2}g^{-1}}$ is favored. Due to the combination of hosts' tide and self-interactions, the central density of subhalos with small pericenter shows a noticeable difference between the cold and the self-interacting models, indicating that dwarf satellites with small pericenter are ideal sites to further constrain the nature of dark matter by future large spectroscopic surveys.

We have revisited the problem of metallicity prediction of RR Lyrae stars from their near-infrared light curves in the Cousins I waveband. Our study is based on high-quality time-series photometry and state-of-the-art high-resolution spectroscopic abundance measurements of 80 fundamental-mode (RRab) and 24 first-overtone (RRc) stars, spanning $\sim$[$-2.7$,$+0.18$] dex and $\sim$[$-3$,$-0.5$] dex ranges, respectively. Employing machine-learning methods, we investigated various light-curve representations and regression models to identify their optimal form for our objective. Accurate new empirical relations between the [Fe/H] iron abundance and the light-curve parameters have been obtained using Bayesian regression for both RRab and RRc stars with mean absolute prediction errors of 0.16 and 0.18 dex, respectively. We found that earlier $I$-band [Fe/H] estimates had a systematic positive bias of up to $\sim 0.4$ dex in the metal-poor regime. Our new predictive models were deployed on large ensembles of RR Lyrae stars to obtain photometric metallicity distribution functions (MDFs) for various old stellar populations in and around the Milky Way. We find that the mode of the old bulge component's MDF is approximately $-1.4$ dex, in remarkable agreement with the latest spectroscopic result. Furthermore, we derive MDF modes of $-1.83$, $-2.13$, and $-1.77$ dex for the Large and Small Magellanic Clouds, and the Sagittarius dwarf galaxy, respectively.

Collisional threats posed by Near-Earth Objects (NEOs) are increasingly being confirmed by National Aeronautics and Space Administration (NASA) and European Space Agency (ESA) sky surveys. Efforts to develop tools to perform modelling, analysis and prediction of possible future impact events are ongoing. The aim of this research is to design a model for Large Earth impact events, describing atmospheric entry, ground impact and post-impact dynamics for impactors large enough to cause harm to the Earth. We present results of Large Earth Impact (LEI) events involving a number of objects with densities ranging between $1000 - 8000 kg m^{-3}$. We use Python code to solve differential equations and perform calculations on analytical models built into the collision simulation algorithms while feeding relevant physical input in terms of object velocity and diameter. We find that denser impactors have slower rates of decrease in momentum while hitting the atmosphere. The denser impactors which exhibit greater resistance to atmospheric disruption result in more energetic ground impacts events in terms of crater formation, thermal radiation generation, seismic effects, ejecta displacement and air-blast effects. We conclude that impacting objects with high masses maintain greater kinetic energies compared to objects with less mass thus posing greater collisional threats to the Earth.

A model of the brightness of the visorsat Starlink spacecraft is presented based on published information on the engineering design and from analysis of 131 observations of individual visorsats in late 2020. Comments are offered on the implications of this model on the visibility of visorsat spacecraft across the sky. This is an updated and expanded version of analysis published in Research Notes of the AAS (2020). An additional section has been added in this version to consider observations made in June 2021 which indicate brighter visorsat magnitudes.

Observing simultaneously pulsed radio and gamma-ray emission from these stars helps to constrain the geometry and radiation mechanisms within their magnetosphere and to localize the multiple photon production sites. In this paper, we fit the time-aligned gamma-ray light-curves of young radio-loud gamma-ray pulsars. We assume a dipole force-free magnetosphere where radio photons emanate from high altitude above the polar caps and gamma-rays originate from outside the light-cylinder, within the striped wind current sheet. We compute a full atlas of radio and gamma-ray pulse profiles depending on the magnetic axis obliquity and line of sight inclination with respect to the neutron star rotation axis. By applying a $\rchi^2$ fitting technique, we are able to pin down accurately the magnetosphere geometry. Further constrains are obtained from radio polarization measurement following the rotating vector model, including aberration and retardation effects. We found a good agreement between our model and the time-aligned single or double peaked gamma-ray pulsar observations. We deduce the magnetic inclination angle and the observer line of sight with respect to the rotation axis within a small error bar. The distinction between radio-loud or radio-quiet gamma-ray pulsars or only radio pulsars can entirely be related to the geometry of the associated emitting regions. The high altitude polar cap model combined with the striped wind represents a minimalistic approach able to reproduce a wealth of gamma-ray pulse profiles for young radio pulsars. Based on self-consistent force-free simulations, it gives a full geometrical picture of the emission properties without resorting to detailed knowledge of the individual particle dynamics and energetics.

Some catastrophic stellar explosions, such as thermonuclear supernovae, core collapse supernovae, compact binary coalescences, and micro-tidal disruption events, are believed to be embedded in the accretion disks of active galactic nuclei (AGN). We show high-energy neutrinos can be produced efficiently through $pp$-interactions between shock-accelerated cosmic rays and AGN disk materials shortly after the explosion ejecta shock breaks out of the disk. AGN stellar explosions are ideal targets for joint neutrino and electromagnetic (EM) multimessenger observations. Future EM follow-up observations of neutrino bursts can help us search for yet-discovered AGN stellar explosions. We suggest that AGN stellar explosions could potentially be important astrophysical neutrino sources. The contribution from AGN stellar explosions to the observed diffuse neutrino background depends on the uncertain local event rate densities of these events in AGN disks. By considering binary white dwarf mergers and neutron star mergers in AGN disks with known theoretical local event rate densities, we show that these events can contribute to $\lesssim10\%$ of the observed diffuse neutrino background.

The 'cold' main classical Kuiper Belt consists of those small solar system bodies with low orbital inclinations and orbital semi-major axes between 42.4 and 47.7~au. Various arguments suggest that these objects formed \textit{in situ} and the original population has experienced minimal collisional modification since their formation. Using the Outer Solar System Origins Survey (OSSOS) ensemble sample and characterization, combined with constraints on the number of small cold classical objects from deeper surveys and supported by evidence from the Minor Planet Center catalog, we determine the absolute magnitude $H_r$ distribution of the cold classical belt from $H_r\simeq5$ to 12 (roughly diameters of 400 km to 20 km). We conclude that the cold population's size distribution exhibits an exponential cutoff at large sizes. Exponential cutoffs at large sizes are not a natural outcome of pair-wise particle accretion but exponentially tapered power-law size distributions are a feature of numerical simulations of planetesimal formation via a streaming instability. Our observation of an exponential cutoff agrees with previous observational inferences that no large objects ($D \gtrsim 400$~km) exist in the cold population. Studies of the transneptunian region are providing the parameters that will enable future streaming-instability studies to determine the initial conditions of planetesimal formation in the $\approx 45$~au region of the Sun's protoplanetary disk.

Since the 1950s, quasi-periodic oscillations have been studied in the terrestrial equatorial stratosphere. Other planets of the solar system present (or are expected to present) such oscillations, like the Jupiter Equatorial Oscillation(JEO) and the Saturn Semi-Annual Oscillation (SSAO). In Jupiter's stratosphere, the equatorial oscillation of its relative temperature structure about the equator, is characterized by a quasi-period of 4.4 years. The stratospheric wind field in Jupiter's equatorial zone has never been directly observed. In this paper, we aim at mapping the absolute wind speeds in Jupiter's equatorial stratosphere to quantify vertical and horizontal wind and temperature shear. Assuming geostrophic equilibrium, we apply the thermal wind balance using nearly simultaneous stratospheric temperature measurements between 0.1 and 30 mbar performed with Gemini/TEXES and direct zonal wind measurements derived at 1 mbar from ALMA observations, all carried out between March 14th and 22nd, 2017. We are thus able to calculate self-consistently the zonal wind field in Jupiter's stratosphere where the JEO occurs. We obtain stratospheric map of the zonal wind speeds as a function of latitude and pressure about Jupiter's equator for the first time. The winds are vertically layered with successive eastward and westward jets. We find a 200 m/s westward jet at 4 mbar at the equator, with a typical longitudinal variability on the order of ~50 m/s. By extending our wind calculations to the upper troposphere, we find a wind structure qualitatively close to the wind observed using cloud-tracking techniques. Nearly simultaneous temperature and wind measurements, both in the stratosphere, are a powerful tool for future investigations of the JEO (and other planetary equatorial oscillations) and its temporal evolution.

We aim to study the structure and kinematics of the broad line region (BLR) of a sample of 27 gravitationally lensed quasars with up to five different epochs of observation. This sample is composed of ~100 spectra from the literature plus 22 unpublished spectra of 11 systems. We measure the magnitude differences in the broad emission line (BEL) wings and statistically model the distribution of microlensing magnifications to determine a maximum likelihood estimate for the sizes of the C IV, C III], and Mg II emitting regions. The BELs in lensed quasars are expected to be magnified differently owing to the different sizes of the regions from which they originate. Focusing on the most common BELs in our spectra (C IV, C III], and Mg II), we find that the low-ionization line Mg II is only weakly affected by microlensing. In contrast, the high-ionization line C IV shows strong microlensing in some cases, indicating that its emission region is more compact. Thus, the BEL profiles are deformed differently depending on the geometry and kinematics of the corresponding emitting region. We detect microlensing in either the blue or the red wing (or in both wings with different amplitudes) of C IV in more than 50% of the systems and find outstanding asymmetries in the wings of QSO 0957+561, SDSS J1004+4112, SDSS J1206+4332, and SDSS J1339+1310. This observation indicates that the BLR is, in general, not spherically symmetric and supports the existence of two regions in the BLR, one insensitive to microlensing and another that only shows up when it is magnified by microlensing.

We use a smoothed particle hydrodynamics (SPH) code to examine the effects of a binary companion on a Be star disk for a range of disk viscosities and misalignment angles, i.e. the angle between the orbital plane and the primary's spin axis. The density structures in the disk due to the tidal interaction with the binary companion are investigated. Expanding on our previous work, the shape and density structure of density enhancements due to the binary companion are analyzed and the changes in observed interferometric features due to these orbiting enhancements are also predicted. We find that larger misalignment angles and viscosity values result in more tightly wound spiral arms with densities that fall-off more slowly with radial distance from the central star. We show that the orbital phase has very little effect on the structure of the spiral density enhancements. We demonstrate that these spiral features can be detected with an interferometer in H$\alpha$ and K-band emission. We also show that the spiral features affect the axis ratios determined by interferometry depending on the orientation of these features and the observer. For example, our simulations show that the axis ratios can vary by 20% for our co-planar binary disk system depending on the location of the disk density enhancements.

We consider galaxy halos formed by dark matter bosons with mass in the range of about a few tens or hundreds eV. A major part of the particles is in a noncondensed state and described under the Thomas-Fermi approach. Derived equations are solved numerically to find the halo density profile. The noncondensed state is supported in the entire halo except compact gravitationally bounded Bose-Einstein condensates. Although the size of these compact objects, also known as Bose stars, depends on interactions between the particles, its upper limit is only about 100 astronomical units. The Bose stars collect the condensed bosons providing a density cusp avoidance in the halo as well as a natural mechanism to prevent overproduction of small halos. Clusters of the Bose stars can also contribute to the halo density profile. The model is analyzed by confronting its predictions with observations of galaxy rotation curves. We employ 22 low surface brightness galaxies and obtain that the model is consistent with the observational data when the particle mass is in the range above about 50 eV and the best fit corresponds to the mass $m=86$ eV. This mass is appropriate for relic dark matter bosons, which decouple just after QCD phase transition.

Many globular clusters (GCs) are known to host multiple populations distinguishable by their light-element content. Less common are GCs displaying iron abundance spreads which are seen as evidence for enrichment through core collapse supernovae (SNe). We present a simple analytical method to estimate the number of SNe required to have occurred in a GC from its metallicity and iron abundance spread. We then use this result to estimate how long star formation (SF) lasted to build the GC. We apply our method to up-to-date measurements and find that, assuming the correctness of these measurements, multiple SNe (up to $10^5$) are required in most GCs with iron abundance spreads. The number of SNe events which contributed to the enrichment of the GCs studied here is typically a factor of 10 less than the expected number of SNe in a canonical initial mass function (IMF). This indicates that gas expulsion from the forming GC occurred after the first 10 per cent of SNe exploded. We compute that for the GCs typically SF ends after only a few Myr (extending up to $\approx 30 ~\rm Myr$ in a few cases). We also discuss possible improvements of this method and especially its sensitivity to the error of iron abundance measurements of individual stars of a GC. The method presented here can quickly give an estimate for the number of SNe required to explain the iron abundance spread in a GC without the requirement of any hydrodynamical simulations.

We present radio observations (1--40 GHz) for 36 classical novae, representing data from over five decades compiled from the literature, telescope archives, and our own programs. Our targets display a striking diversity in their optical parameters (e.g., spanning optical fading timescales, t_2 = 1--263 days), and we find a similar diversity in the radio light curves. Using a brightness temperature analysis, we find that radio emission from novae is a mixture of thermal and synchrotron emission, with non-thermal emission observed at earlier times. We identify high brightness temperature emission (T_B > 5x10^4 K) as an indication of synchrotron emission in at least 9 (25%) of the novae. We find a class of synchrotron-dominated novae with mildly evolved companions, exemplified by V5589 Sgr and V392 Per, that appear to be a bridge between classical novae with dwarf companions and symbiotic binaries with giant companions. Four of the novae in our sample have two distinct radio maxima (the first dominated by synchrotron and the later by thermal emission), and in four cases the early synchrotron peak is temporally coincident with a dramatic dip in the optical light curve, hinting at a common site for particle acceleration and dust formation. We publish the light curves as tables and encourage use of these data by the broader community in multi-wavelength studies and modeling efforts.

We conducted interferometric observations with the CHARA Array of transiting super-Earth host HD 97658 and measured its limb-darkened angular diameter to be $\theta_{\text{LD}}=0.314\pm0.004$ mas. The combination of the angular diameter with the Gaia EDR3 parallax value with zero-point correction ($\pi=46.412\pm0.022$ mas, $d=21.546\pm0.011$ pc) yields a physical radius of $R_\star=0.728\pm0.008$ $R_\odot$. We also measured the bolometric flux of the star to be $F_\text{bol}=2.42\pm 0.05\times 10^{-8}$erg s$^{-1}$ cm$^{-2}$, which, together with angular size, allows a measurement of the effective temperature $T_{\text{eff}}=5212\pm43$ K. Our directly determined physical stellar properties are in good agreement with previous estimates derived from spectroscopy. We used our measurements in combination with stellar evolutionary models and properties of the transit of HD 97658 b to determine the mass and age of HD 97658 as well as constrain the properties of the planet. Our results and our analysis of the TESS lightcurve on the planet (TOI-1821) corroborate previous studies of this system with tighter uncertainties.

Asteroid surfaces are subjected to mechanical weathering processes that result in the development and evolution of regolith. Two proposed mechanisms--impact bombardment and thermal fatigue--have been proposed as viable and dominant weathering processes. Previously, we compiled and estimated thermal inertias of several hundred asteroids (mostly in the main-belt) for which we determined dependencies on temperature, diameter, and rotation period. In this work, we estimate grain sizes of asteroid regoliths from this large thermal inertia dataset using thermal conductivity models. Following our previous work we perform multi-variate linear model fits to the grain size dataset and quantify its dependency on diameter and rotation period. We find that the best-fit model fit indicates that asteroid grain sizes are inversely dependent on object size for <10 km asteroids and exhibits no relationship above this size cutoff. Rotation period and grain size show a positive relationship when the rotation period is greater than ~5 hr, and an inverse relationship below this rotation period. We conclude that both impact weathering and thermal fatigue are significant regolith evolution mechanisms. Furthermore, we run post-hoc t-tests between spectral groups to identify compositional differences among our asteroid set. Notably, suspected metal-rich, M-type and E-type asteroids have larger than expected grain sizes, and P-types have distinctly smaller grains than other groups

We study to what extend LISA can observe features of gravitational wave spectra originating from cosmological first-order phase transitions. We focus on spectra which are of the form of double-broken power laws. These spectra are predicted by hydrodynamic simulations and also analytical models such as the sound shell model. We argue that the ratio of the two break frequencies is an interesting observable since it can be related to the wall velocity while overall amplitude and frequency range are often degenerate for the numerous characteristics of the phase transition. Our analysis uses mock data obtained from the power spectra predicted by the simplified simulations and the sound shell model and analyzes the detection prospects using $\chi^2$-minimization and likelihood sampling. We point out that the prospects of observing two break frequencies from the electroweak phase transition is hindered by a shift of the spectrum to smaller frequencies for strong phase transitions. Finally, we also highlight some differences between signals from the sound shell model compared to simulations.

In this work, we consider the scalar-tensor theory that contains a scalar field with its kinetic and potential terms minimally coupled to gravity, while the scalar field is assumed to have a coulombic form. In the context of this theory, we obtain analytic, asymptotically flat and regular (ultra-compact) black-hole solutions with non-trivial scalar hair of secondary type. At first, we examine the properties of the static and spherically symmetric black-hole solution and we find that in the causal region of the spacetime the stress-energy tensor, needed to support our solution, satisfies the strong energy conditions. Then, by using the slow-rotating approximation, we generalize the static solution into a slowly rotating one, and we determine explicitly its angular velocity $\omega(r)$. We also find that the angular velocity of our ultra-compact solution is always larger compared to the angular velocity of the corresponding equally massive slow-rotating Schwarzschild black hole. In addition, we investigate the axial perturbations of the derived solutions by determining the Schr\"{o}ndiger-like equation and the effective potential. We show that there is a region in the parameter space of the free parameters of our theory, which allows for the existence of stable ultra-compact black hole solutions. Specifically, we calculate that the most compact and stable black hole solution is 0.551 times smaller than the Schwarzschild one, while it rotates 2.491 times faster compared to the slow-rotating Schwarzschild black hole. Finally, we present without going into details the generalization of the derived asymptotically flat solutions to asymptotically (A)dS solutions.

Constraints on the indirect detection of dark matter are usually obtained from observations of astrophysical objects -- the Galactic Center, dwarf galaxies, M31, etc. Here we propose instead to look for the annihilation or decay of dark matter particles taking place inside detectors searching \emph{directly} for dark matter or in large neutrino experiments. We show that the data from XENON1T and Borexino set limits on the annihilation and decay rate of dark matter particles with masses in the keV to few MeV range. All relevant final states are considered: annihilation into $\gamma\gamma$ and $e^-e^+$, and decays into $\gamma\gamma$, $\gamma\nu$ and $e^-e^+$. The expected sensitivities in XENONnT, DARWIN, JUNO and THEIA are also computed. Though weaker than current astrophysical bounds, the laboratory limits (and projections) obtained are free from the usual astrophysical uncertainties associated with $J$-factors and unknown backgrounds, and may thus offer a complementary probe of the dark matter properties. We point out that current and future (astro)particle physics detectors might also be used to set analogous limits for different decays and dark matter masses above a few MeV.

Photons emitted by light sources in the neighbourhood of a black hole can wind several times around it before fleeing towards the observer. For spherically symmetric black holes, two infinite sequences of images are created for any given source, asymptotically approaching the shadow border with decreasing magnitude. These sequences are reflected by a characteristic staircase structure in the complex visibility function that may be used to decode the properties of the black hole metric. Recalling the formalism of gravitational lensing in the strong deflection limit, we derive analytical formulae for the height, the width and the periodicities of the steps in the visibility as functions of the black hole parameters for the case of a single compact source. With respect to diffuse emission by the whole accretion flow, this ideal framework provides clean insight and model-independent information on the metric. These basic formulae can then be used to build visibilities for more complicated sources and track the changes induced by alternative metrics and ultimately test General Relativity. As simple examples, we include visibilities for Reissner-Nordstr\"om and Janis-Newman-Winicour metrics.

We study some consequences of the loop quantization of the outermost dust shell in the Lema\^itre-Tolman-Bondi spacetime with a homogeneous dust density using different quantization strategies motivated by loop quantum gravity. Prior work has dealt with loop quantizing this model by employing holonomies and the triads, following the procedure in standard loop quantum cosmology. In this work we compare this quantization with the one in which holonomies and gauge-covariant fluxes are used. While both of the quantization schemes resolve the central singularity, they lead to different mass gaps at which a trapped surface forms. This trapped surface which is matched to an exterior generalized Vaidya spacetime disappears when the density of the dust shell is in the Planck regime. We find that the quantization based on holonomies and gauge-covariant fluxes generically results in an asymmetric evolution of the dust shell in which the effective mass associated with the white hole as seen by an external observer is $2/\pi$ of the one for the black hole. This effective difference in masses results from difference in the classical limits in pre- and post-bounce regimes in the two quantizations. This distinctive feature rules out formation of any black hole-white hole twin in presence of gauge-covariant flux modifications which is in contrast to the quantization using holonomies and triads where the gravitational collapse always leads to a black hole-white hole twins. In another striking difference, for the quantization based on holonomies and gauge-covariant fluxes there can be situations in which during a non-singular collapse only a black hole forms without a white hole.

We study dynamical friction in the Newtonian regime of nonlocal gravity (NLG), which is a classical nonlocal generalization of Einstein's theory of gravitation. The nonlocal aspect of NLG simulates dark matter. The attributes of the resulting effective dark matter are described and the main physical predictions of nonlocal gravity, which has a characteristic lengthscale of order 1 kpc, for galactic dynamics are presented. Within the framework of NLG, we derive the analogue of Chandrasekhar's formula for dynamical friction. The astrophysical implications of the results for the apparent rotation of a central bar subject to dynamical friction in a barred spiral galaxy are briefly discussed.

We propose a novel $k$-Gauss-Bonnet model, in which a kinetic term of scalar field is allowed to non-minimally couple to the Gauss-Bonnet topological invariant in the absence of a potential of scalar field. As a result, this model is shown to admit an isotropic power-law inflation provided that the scalar field is phantom. Furthermore, stability analysis based on the dynamical system method is performed to indicate that this inflation solution is indeed stable and attractive. More interestingly, a gradient instability in tensor perturbations is shown to disappear in this model.

An analytical solution representing traversable asymptotically flat and symmetric wormholes was obtained without adding exotic matter in two different theories independently: in the Einstein-Maxwell-Dirac theory and in the second Randall-Sundrum brane-world model. Further, a smooth normalizable asymmetric wormhole solution has been recently obtained numerically in the Einstein-Maxwell-Dirac theory. Using the time-domain integration method we study quasinormal ringing of all these wormholes with emphasis to the regime of mimicking the near extremal Reissner-Nordstr\"om black holes, which is characterised by echoes. In addition, we calculate radius of shadows cast by these wormholes.

Renormalization group (RG) applications to cosmological problems often encounter difficulties in the interpretation of the field independent term in the effective potential. While this term is constant with respect to field variations, it generally depends on the RG scale k. Since the RG running could be associated with the temporal evolution of the Universe according to the identification $k \sim 1/t$, one can treat the field independent constant, i.e., the $\Lambda$ term in Einstein's equations as a running (scale-dependent) parameter. Its scale dependence can be described by nonperturbative RG, but it has a serious drawback, namely $k^4$ and $k^2$ terms appear in the RG flow in its high-energy (UV) limit which results in a rampant divergent behaviour. Here, we propose a subtraction method to handle this unphysical UV scaling and provides us a framework to build up a reliable solution to the cosmological constant problem.

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not only probed in astrophysical observations, but also in terrestrial heavy-ion collision experiments. In this work, we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars and from heavy-ion collisions of gold nuclei at relativistic energies with microscopic nuclear theory calculations to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent NICER observations. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment, and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

We find that sub-GeV neutrino portal bosons that carry lepton number can condense inside a proto-neutron star (newly born neutron star). These bosons are produced copiously and form a Bose-Einstein condensate for a range of as yet unconstrained coupling strengths to neutrinos. The condensate is a lepton number superfluid with transport properties that differ dramatically from those encountered in ordinary dense baryonic matter. We discuss how this phase could alter the evolution of proto-neutron stars and comment on the implications for neutrino signals and nucleosynthesis.