Papers for Tuesday, Nov 30 2021

The Nuclear Stellar Disc (NSD) is a flattened high-density stellar structure that dominates the gravitational field of the Milky Way at Galactocentric radius $30\lesssim R\lesssim 300$ pc. We construct axisymmetric self-consistent equilibrium dynamical models of the NSD in which the distribution function is an analytic function of the action variables. We fit the models to the normalised kinematic distributions (line-of-sight velocities + VIRAC2 proper motions) of stars in the NSD survey of Fritz et al., taking the foreground contamination due to the Galactic Bar explicitly into account using an $N$-body model. The posterior marginalised probability distributions give a total mass of $M_{\rm NSD} = 10.5^{+0.8}_{-0.7} \times 10^8 \,{\rm M_\odot}$, roughly exponential radial and vertical scale-lengths of $R_{\rm disc} = 88.6^{+9.2}_{-6.9}$ pc and $H_{\rm disc}=28.4^{+2.2}_{-2.1}$ pc respectively, and a velocity dispersion $\sigma \simeq 70$ km/s that decreases with radius. We find that the assumption that the NSD is axisymmetric provides a good representation of the data. We quantify contamination from the Galactic Bar in the sample, which is substantial in most observed fields. Our models provide the full 6D (position+velocity) distribution function of the NSD, which can be used to generate predictions for future surveys. We make the models publicly available as part of the software package \textsc{Agama}.

Hot and Warm Jupiters (HJs\&WJs) are gas-giant planets orbiting their host stars at short orbital periods, posing a challenge to an efficient in-situ formation of such planets. Most of the HJs/WJs populations are thought to have migrated from an initially farther-out birth locations. Both disc-migration (driven by gas-dissipation) and eccentric-migration (driven by tidal evolution) fail to produce the occurrence rate and orbital properties of HJs/WJs. Here we study the role of the thermal evolution and its coupling to tidal evolution. We use the AMUSE numerical environment and MESA planetary evolution modeling to model in detail the coupled internal and orbital evolution of gas-giants during their eccentric migration. In this paper we use a numerical full planetary evolution thermal-dynamical model. In a companion paper, we use a simple semi-analytic model, validated by our numerical model, and run a population synthesis study. We consider the initially inflated radii of gas-giants (expected following their formation), as well study the effects of the potential slowed contraction and even re-inflation of gas-giants (due to tidal and radiative heating) on the eccentric migration. Tidal forces that drive eccentric-migration are highly sensitive to the planetary structure and radius. Consequently, we find that this form of inflated eccentric-migration operates on significantly (up to an order of magnitude) shorter timescales than previously studied eccentric migration models. Thereby, inflated eccentric migration gives rise to more rapid formation of HJs/WJs and higher occurrence rates of WJs, as well as higher rates of tidal disruptions, compared with previous eccentric migration models which consider constant ~Jupiter radii for HJ/WJ progenitors. Coupled thermal-dynamical evolution of eccentric gas-giants can therefore play a key-role in their evolution.

Hot and warm Jupiters (HJs and WJs correspondingly) are gas-giants orbiting their host stars at very short orbital periods ($P_{HJ}<10$ days; $10<P_{WJ}<200$ days). HJs and a significant fraction of WJs are thought to have migrated from an initially farther-out birth locations. While such migration processes have been extensively studied, the thermal evolution of gas-giants and its coupling to the the migration processes are usually overlooked. In particular, gas-giants end their core-accretion phase with large radii and then contract slowly to their final radii. Moreover, intensive heating can slow the contraction at various evolutionary stages. The initial inflated large radii lead to faster tidal migration due to the strong dependence of tides on the radius. Here we explore this accelerated migration channel, which we term inflated eccentric migration, using a semi-analytical self-consistent modeling of the thermal-dynamical evolution of the migrating gas-giants, later validated by our numerical model (see a companion paper, paper II). We demonstrate our model for specific examples and carry a population synthesis study. Our results provide a general picture of the properties of the formed HJs\&WJs via inflated migration, and the dependence on the initial parameters/distributions. We show that tidal migration of gas-giants could occur far more rapidly then previously thought and lead to accelerated destruction and formation of HJs and enhanced formation rate of WJs. Accounting for the coupled thermal-dynamical evolution is therefore critical to the understanding of HJs/WJs formation, evolution and final properties of the population and play a key role in their migration process.

NGC 6388 is one of the most massive Galactic globular clusters (GC) and it is an old, metal-rich, Galactic bulge cluster. By exploiting previous spectroscopic observations we were able to bypass the uncertainties in membership related to the strong field stars contamination. We present the abundance analysis of 12 new giant stars with UVES spectra and 150 giants with GIRAFFE spectra acquired at the ESO-VLT. We derived radial velocities, atmospheric parameters and iron abundances for all stars. When combined to previous data, we obtain a grand total of 185 stars homogeneously analysed in NGC 6388 from high-resolution spectroscopy. The average radial velocity of the 185 stars is 81.2+/-0.7, rms = 9.4 km/s. We obtain an average metallicity [Fe/H]=-0.480 dex, rms = 0.045 dex (35 stars) and [Fe/H]=-0.488 dex, rms = 0.040 dex (150 stars) from the UVES and GIRAFFE samples, respectively. Comparing these values to internal errors in abundance, we exclude the presence of a significant intrinsic metallicity spread within the cluster. Since about a third of giants in NGC 6388 is claimed to belong to the "anomalous red giants" in the HST pseudo-colour map defining the so-called type-II GCs, we conclude that either enhanced metallicity is not a necessary requisite to explain this classification (as also suggested by the null iron spread for NGC 362) or NGC 6388 is not a type-II globular cluster.

Over the course of their lifetimes, the rotation of solar-type stars goes through different phases. Once they reach the zero-age main sequence (ZAMS), their global rotation rate decreases during the main sequence until at least the solar age, approximately following the empirical Skumanich law and enabling gyrochronology. Older solar-type stars might then reach a point of transition when they stop braking, according to recent results of asteroseismology. Additionally, recent 3D numerical simulations of solar-type stars show that different regimes of differential rotation can be characterized with the Rossby number. In particular, anti-solar differential rotation (fast poles, slow equator) may exist for high Rossby number (slow rotators). If this regime occurs during the main sequence and, in general, for slow rotators, we may consider how magnetic generation through the dynamo process might be impacted. In particular, we consider whether slowly rotating stars are indeed subject to magnetic cycles. We find that kinematic $\alpha$ $\Omega$ dynamos allow for the presence of magnetic cycles and global polarity reversals for both rotation regimes, but only if the $\alpha$-effect is saddled on the tachocline. If it is distributed in the convection zone, solar-type cases still possess a cycle and anti-solar cases do not. Conversely, we have not found any possibility for sustaining a magnetic cycle with the traditional Babcock-Leighton flux-transport dynamos in the anti-solar differential rotation regime due to flux addition. Graphic interpretations are proposed in order to illustrate these cases. However, we find that hybrid models containing both prescriptions can still sustain local polarity reversals at some latitudes. We conclude that stars in the anti-solar differential rotation regime can sustain magnetic cycles only for very specific dynamo processes.

We examine whether stellar flyby simulations can initiate FU Orionis outbursts using 3D hydrodynamics simulations coupled to live Monte Carlo radiative transfer. We find that a flyby where the secondary penetrates the circumprimary disc triggers a 1-2 year rise in the mass accretion rate to $10^{-4}~{\rm M_{\odot}~ yr^{-1}}$ that remains high ($\gtrsim 10^{-5}~{\rm M_{\odot}~yr^{-1}}$) for more than a hundred years, similar to the outburst observed in FU Ori. Importantly, we find that the less massive star becomes the dominant accretor, as observed.

Cosmic Microwave Background (CMB) observations have been used extensively to constrain key properties of neutrinos, such as their mass. However, these inferences are typically dependent on assumptions about the cosmological model, and in particular upon the distribution function of neutrinos in the early Universe. In this paper, we aim to assess the full extent to which CMB experiments are sensitive to the shape of the neutrino distribution. We demonstrate that Planck and CMB-S4-like experiments have no prospects for detecting particular features in the distribution function. Consequently, we take a general approach and marginalise completely over the form of the neutrino distribution to derive constraints on the relativistic and non-relativistic neutrino energy densities, characterised by $N_\mathrm{eff} = 3.0 \pm 0.4$ and $\rho_{\nu,0}^{\rm NR} < 14 \, \mathrm{eV}\,\mathrm{cm}^{-3}$ at 95% CL, respectively. The fact that these are the only neutrino properties that CMB data can constrain has important implications for neutrino mass limits from cosmology. Specifically, in contrast to the $\Lambda$CDM case where CMB and BAO data tightly constrain the sum of neutrinos masses to be $\sum m_\nu < 0.12 \, \mathrm{eV}$, we explicitly show that neutrino masses as large as $\sum m_\nu \sim 3 \, \mathrm{eV}$ are perfectly consistent with this data. Importantly, for this to be the case, the neutrino number density should be suitably small such that the bound on $\rho_{\nu,0}^\mathrm{NR} = \sum m_\nu n_{\nu,0}$ is still satisfied. We conclude by giving an outlook on the opportunities that may arise from other complementary experimental probes, such as galaxy surveys, neutrino mass experiments and facilities designed to directly detect the cosmic neutrino background.

The temperature of an atmosphere decreases with increasing altitude, unless a shortwave absorber exists that causes a temperature inversion. Ozone plays this role in the Earth`s atmosphere. In the atmospheres of highly irradiated exoplanets, shortwave absorbers are predicted to be titanium oxide (TiO) and vanadium oxide (VO). Detections of TiO and VO have been claimed using both low and high spectral resolution observations, but later observations have failed to confirm these claims or overturned them. Here we report the unambiguous detection of TiO in the ultra-hot Jupiter WASP-189b using high-resolution transmission spectroscopy. This detection is based on applying the cross-correlation technique to many spectral lines of TiO from 460 to 690 nm. Moreover, we report detections of metals, including neutral and singly ionised iron and titanium, as well as chromium, magnesium, vanadium and manganese (Fe, Fe+, Ti, Ti+, Cr, Mg, V, Mn). The line positions of the detected species differ, which we interpret as a consequence of spatial gradients in their chemical abundances, such that they exist in different regions or dynamical regimes. This is direct observational evidence for the three-dimensional thermo-chemical stratification of an exoplanet atmosphere derived from high-resolution ground-based spectroscopy.

High-Velocity Clouds (HVCs) are believed to be an important source of gas accretion for star formation in the Milky Way. Earlier numerical studies have found that the Galactic magnetic field and radiative cooling strongly affects accretion. However, these effects have not previously been included together in the context of clouds falling through the Milky Way's gravitational potential. We explore this by simulating an initially stationary cloud falling through the hot hydrostatic corona towards the disc. This represents a HVC that has condensed out of the corona. We include the magnetic field in the corona to examine its effect on accretion of the HVC and its associated cold gas. Remnants of the original cloud survive in all cases, although a strong magnetic field causes it to split into several fragments. We find that mixing of cold and hot gas leads to cooling of coronal gas and an overall growth with time in cold gas mass, despite the low metallicity of the cloud and corona. The role of the magnetic field is to (moderately to severely) suppress the mixing and subsequent cooling, which in turn leads to less accretion compared to when the field is absent. A stronger field leads to less suppression of condensation because it enhances Rayleigh-Taylor instability. However, magnetic tension in a stronger field substantially decelerates condensed cloudlets. These have velocities typically a factor 3-8 below the velocity of the main cloud remnants by the end of the simulation. Some of these cloudlets likely disperse before reaching the disc.

The distribution of orbital energies imparted into stellar debris following the close encounter of a star with a supermassive black hole is the principal factor in determining the rate of return of debris to the black hole, and thus in determining the properties of the resulting lightcurves from such events. We present simulations of tidal disruption events for a range of $\beta\equiv r_{\rm t}/r_{\rm p}$ where $r_{\rm p}$ is the pericentre distance and $r_{\rm t}$ the tidal radius. We perform these simulations at different spatial resolutions to determine the numerical convergence of our models. We compare simulations in which the heating due to shocks is included or excluded from the dynamics. For $\beta \lesssim 8$ the simulation results are well-converged at sufficiently moderate-to-high spatial resolution, while for $\beta \gtrsim 8$ the breadth of the energy distribution can be grossly exaggerated by insufficient spatial resolution. We find that shock heating plays a non-negligible role only for $\beta \gtrsim 4$, and that typically the effect of shock heating is mild. We show that self-gravity can modify the energy distribution over time after the debris has receded to large distances for all $\beta$. Primarily, our results show that across a range of impact parameters, while the shape of the energy distribution varies with $\beta$, the width of the energy spread imparted to the bulk of the debris is closely matched to the canonical spread, $\Delta E = GM_\bullet R_\star/r_{\rm t}^2$, for the range of $\beta$ we have simulated.

We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star, $r_{\rm p}$, is well within the tidal radius of the black hole, $r_{\rm t}$, i.e., when $\beta \equiv r_{\rm t}/r_{\rm p} \gg 1$. We find that shocks form for $\beta \gtrsim 3$, but they are weak (with Mach numbers $\sim 1$) for all $\beta$, and that they reach the center of the star prior to the time of maximum adiabatic compression for $\beta \gtrsim 10$. The maximum density and temperature reached during the TDE follow much shallower relations with $\beta$ than the previously predicted $\rho_{\rm max} \propto \beta^3$ and $T_{\rm max} \propto \beta^2$ scalings. Below $\beta \simeq 10$, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the freefalling gas, while above $\beta \simeq 10$ we find that shocks prematurely halt the compression and yield the scalings $\rho_{\rm max} \propto \beta^{1.62}$ and $T_{\rm max} \propto \beta^{1.12}$. We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic vs.~non-polytropic stars, and general relativistic effects, which become important when the pericenter of the star nears the direct capture radius, at $\beta \sim 12.5$ (2.7) for a solar-like star disrupted by a $10^6M_{\odot}$ ($10^{7}M_{\odot}$) supermassive black hole.

We present an analysis of the formation and chemical evolution of stellar haloes around (a) Milky Way Analogue (MWA) galaxies and (b) galaxy clusters in the L-Galaxies 2020 semi-analytic model of galaxy evolution. Observed stellar halo properties are better reproduced when assuming a gradual stripping model for the removal of cold gas and stars from satellites, compared to an instantaneous stripping model. The slope of the stellar mass -- metallicity relation for MWA stellar haloes is in good agreement with that observed in the local Universe. This extends the good agreement between L-Galaxies 2020 and metallicity observations from the gas and stars inside galaxies to those outside. Halo stars contribute on average only $\sim{}0.1$ per cent of the total circumgalactic medium (CGM) enrichment by $z=0$ in MWAs, ejecting predominantly carbon produced by AGB stars. Around a quarter of MWAs have a single `significant progenitor' with a mean mass of $\sim{}2.3\times{}10^{9}M_{\odot}$, similar to that measured for Gaia Enceladus. For galaxy clusters, L-Galaxies 2020 shows good correspondence with observations of stellar halo mass fractions, but slightly over-predicts stellar masses. Assuming a Navarro-Frenk-White profile for the stellar halo mass distribution provides the best agreement. On average, the intracluster stellar component (ICS) is responsible for 5.4 per cent of the total intracluster medium (ICM) iron enrichment, exceeding the contribution from the brightest cluster galaxy (BCG) by $z=0$. We find that considering gradual stripping of satellites and realistic radial profiles is crucial for accurately modelling stellar halo formation on all scales in semi-analytic models.

Cosmic birefringence -- the rotation of the polarization of the cosmic microwave background (CMB) photons as they travel to the Earth -- is a smoking gun for axion-like particles (ALPs) that interact with the photon. It has recently been suggested that topological defects in the ALP field called cosmic strings can cause polarization rotation in quantized amounts that are proportional to the electromagnetic chiral anomaly coefficient $\mathcal A$, which holds direct information about physics at very high energies. In this work, we study the detectability of a random network of cosmic strings through estimating rotation using quadratic estimators (QEs). We show that the QE derived from the maximum likelihood estimator is equivalent to the recently proposed global-minimum-variance QE; the classic Hu-Okamoto QE equals the global-minimum-variance QE under special conditions, but is otherwise still nearly globally optimal. We calculate the sensitivity of QEs to cosmic birefringence from string networks, for the Planck satellite mission, as well as for third- and fourth-generation ground-based CMB experiments. Using published binned rotation power spectrum derived from the Planck 2015 polarization data, we obtain a constraint $\mathcal A^2\,\xi_0 < 0.93$ at the 95\% confidence level, where $\xi_0$ is the total length of strings in units of the Hubble scale per Hubble volume, for a phenomenological but reasonable string network model describing a continuous distribution of string sizes. Looking forward, experiments such as the Simons Observatory and CMB-S4 will either discover or falsify the existence of an ALP string network for the theoretically plausible parameter space $\mathcal A^2\,\xi_0 \gtrsim 0.01$.

The Milky Way Galaxy hosts a four million solar mass black hole, Sgr A*, that underwent a major accretion episode approximately 3-6 Myr ago. During the episode, hundreds of young massive stars formed in a disc orbiting Sgr A* in the central half parsec. The recent discovery of a hypervelocity star S5-HVS1, ejected by Sgr A* five Myr ago with a velocity vector consistent with the disc, suggests that this event also produced binary star disruptions. The initial stellar disc has to be rather eccentric for this to occur. Such eccentric disks can form from the tidal disruptions of molecular clouds. Here we perform simulations of such disruptions, focusing on gas clouds on rather radial initial orbits. As a result, stars formed in our simulations are on very eccentric orbits ($\bar{e}\sim 0.6$) with a lopsided configuration. For some clouds counter-rotating stars are formed. As in previous work, we find that such discs undergo a secular gravitational instability that leads to a moderate number of particles obtaining eccentricities of 0.99 or greater, sufficient for stellar binary disruption. We also reproduce the mean eccentricity of the young disk in the Galactic centre, though not the observed surface density profile. We discuss missing physics and observational biases that may explain this discrepancy. We conclude that observed S-stars, hypervelocity stars, and disc stars tightly constrain the initial cloud parameters, indicating a cloud mass between a few$\times 10^4$ and $10^5 M_{\odot}$, and a velocity between $\sim 40$ and 80 km s$^{-1}$ at 10 pc.

LUCI is an general-purpose spectral line-fitting pipeline which natively integrates machine learning algorithms to initialize fit functions. LUCI currently uses point-estimates obtained from a convolutional neural network (CNN) to inform optimization algorithms; this methodology has shown great promise by reducing computation time and reducing the chance of falling into a local minimum using convex optimization methods. In this update to LUCI, we expand upon the CNN developed in Rhea et al. 2020 so that it outputs Gaussian posterior distributions of the fit parameters of interest (the velocity and broadening) rather than simple point-estimates. Moreover, these posteriors are then used to inform the priors in a Bayesian inference scheme, either emcee or dynesty. The code is publicly available at https://github.com/crhea93/LUCI.

PANOSETI (Pulsed All-Sky Near-infrared Optical Search for Extra Terrestrial Intelligence) is a dedicated SETI (Search for Extraterrestrial Intelligence) observatory that is being designed to observe 4,441 sq. deg. to search for nano- to milli-second transient events. The experiment will have a dual observatory system that has a total of 90 identical optical 0.48 m telescopes that each have a 99 square degree field of view. The two observatory sites will be separated by 1 km distance to help eliminate false positives and register a definitive signal. We discuss the overall mechanical design of the telescope modules which includes a Fresnel lens housing, a shutter, three baffles, an 32x32 array of Hamamatsu Multi-Photon Pixel Counting (MPPC) detectors that reside on a linear stage for focusing. Each telescope module will be housed in a triangle of a 3rd tessellation frequency geodesic dome that has the ability to have directional adjustment to correct for manufacturing tolerances and astrometric alignment to the second observatory site. Each observatory will have an enclosure to protect the experiment, and an observatory room for operations and electronics. We will review the overall design of the geodesic domes and mechanical telescope attachments, as well as the overall cabling and observatory infrastructure layout.

New periodic comet P/2019 LD2 (ATLAS) located on unstable quasi-Trojan orbit still an interesting object to study in the last years. We present the results of broadband observations of comet P/2019 LD2 (ATLAS) performed at the Sanglokh observatory of the Institute of Astrophysics, National Academy of Sciences of Tajikistan for 5 nights in August 2020. The dependence of the Afrho parameters of the comet on the aperture radius is measured. We obtained the coma corrected upper limit of the cometary radius Rn not more than 6.1 +- 0.1 km and calculated the absolute magnitude H0 = 11.41 +- 0.03. Finson - Probstein diagram method was used to explain a dust tail appearance. A subgroup of comets with orbital parameters close to P/2019 LD2 (ATLAS) is analyzed.

Joint ranking statistics are used to distinguish real from random coincidences, ideally considering whether shared parameters are consistent with each other as well as whether the individual candidates are distinguishable from noise. We expand on previous works to include additional shared parameters, we use galaxy catalogues aspriors for sky localization and distance, and avoid some approximations previously used. We develop methods to calculate this statistic both in low-latency using HEALPix sky maps, as well as with posterior samples. We show that these changes lead to a factor of one to two orders of magnitude improvement for GW170817-GRB 170817A depending on the method used, placing this significant event further into the foreground. We also examined the more tenuous joint candidate GBM-GW150914, which was largely penalized by these methods. Finally, we performed a simplistic simulation that argues these changes could better help distinguish between real and random coincidences in searches, although more realistic simulations are needed to confirm this.

ATLASGAL is a 870-mircon dust survey of 420 square degrees of the inner Galactic plane and has been used to identify ~10 000 dense molecular clumps. Dedicated follow-up observations and complementary surveys are used to characterise the physical properties of these clumps, map their Galactic distribution and investigate the evolutionary sequence for high-mass star formation. The analysis of the ATLASGAL data is ongoing: we present an up-to-date version of the catalogue. We have classified 5007 clumps into four evolutionary stages (quiescent, protostellar, young stellar objects and HII regions) and find similar numbers of clumps in each stage, suggesting a similar lifetime. The luminosity-to-mass (L/M) ratio curve shows a smooth distribution with no significant kinks or discontinuities when compared to the mean values for evolutionary stages indicating that the star-formation process is continuous and that the observational stages do not represent fundamentally different stages or changes in the physical mechanisms involved. We compare the evolutionary sample with other star-formation tracers (methanol and water masers, extended green objects and molecular outflows) and find that the association rates with these increases as a function of evolutionary stage, confirming that our classification is reliable. This also reveals a high association rate between quiescent sources and molecular outflows, revealing that outflows are the earliest indication that star formation has begun and that star formation is already ongoing in many of the clumps that are dark even at 70 micron.

We present new reflection models specifically tailored to model the X-ray radiation reprocessed in accretion disks around neutron stars, in which the primary continuum is characterized by a single temperature blackbody spectrum, emitted either at the surface of the star, or at the boundary layer. These models differ significantly from those with a standard power-law continuum, typically observed in most accreting black holes. We show comparisons with earlier reflection models, and test their performance in the NuSTAR observation of the neutron star 4U 1705-44. Simulations of upcoming missions such as XRISM-Resolve and Athena X-IFU are shown to highly the diagnostic potential of these models for high-resolution X-ray reflection spectroscopy. These new reflection models xillverNS, and their relativistic counterpart relxillNS, are made publicly available to the community as an additional flavor in the relxill suite of reflection models.

Many short Gamma-Ray Bursts (sGRBs) have a prolonged plateau in the X-ray afterglow lasting up to tens of thousands of seconds. A central engine injecting energy into the remnant may fuel the plateau. A simple analytic model describing the interaction of the magnetized relativistic wind from a rapidly-rotating magnetar with the surrounding environment can reproduce X-ray plateaux and instantaneous spectra. The model is analogous to classic, well-established models of young supernova remnants and applies the underlying physics to sGRB remnants. The light curve and spectra produced by the model are compared to observations of GRB 130603B. The spectra are also used to estimate parameters of the magnetar including its poloidal field strength and angular frequency. If combined with a gravitational wave signal, this model could provide insight into multimessenger astronomy and neutron star physics.

We predict the three-dimensional intensity power spectrum (PS) of the [CII] 158$\,\mu$m line throughout the epoch of (and post) reionization at redshifts from $\approx$ 3.5 to 8. We study the detectability of the PS in a line intensity mapping (LIM) survey with the Fred Young Submillimeter Telescope (FYST). We created mock [CII] tomographic scans in redshift bins at $z\approx$ 3.7, 4.3, 5.8, and 7.4 using the Illustris TNG300-1 $\Lambda$CDM simulation and adopting a relation between the star formation activity and the [CII] luminosity ($L_{[CII]}$) of galaxies. A star formation rate (SFR) was assigned to a dark matter halo in the Illustris simulation in two ways: (i) we adopted the SFR computed in the Illustris simulation and, (ii) we matched the abundance of the halos with the SFR traced by the observed dust-corrected ultraviolet luminosity function of high-redshift galaxies. The $L_{[CII]}$ is related to the SFR from a semi-analytic model of galaxy formation, from a hydrodynamical simulation of a high-redshift galaxy, or from a high-redshift [CII] galaxy survey. The [CII] intensity PS was computed from mock tomographic scans to assess its detectability with the anticipated observational capability of the FYST. The amplitude of the predicted [CII] intensity power spectrum varies by more than a factor of 10, depending on the choice of the halo-to-galaxy SFR and the SFR-to-$L_{[CII]}$ relations. In the planned $4^{\circ} \times 4^{\circ}$ FYST LIM survey, we expect a detection of the [CII] PS up to $z \approx$ 5.8, and potentially even up to $z \approx $ 7.4. The design of the envisioned FYST LIM survey enables a PS measurement not only in small (<10 Mpc) shot noise-dominated scales, but also in large (>50 Mpc) clustering-dominated scales making it the first LIM experiment that will place constraints on the SFR-to-$L_{[CII]}$ and the halo-to-galaxy SFR relations simultaneously.

Based on the data release of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey (LAMOST DR5) and the \emph{Gaia} Early Data Release 3 (\emph{Gaia} EDR3), we construct a sample containing 46,109 giant (log~$\emph{g}$ $\leqslant$ 3.5\,dex) stars with heliocentric distance d $\leqslant$ 4 kpc, and the sample is further divided into two groups of the inner ({\RGC} $<$ 8.34 \,kpc) and outer region ({\RGC} $>$ 8.34\,kpc). The {\LZ} distributions of our program stars in the panels with different [Fe/H] and [$\alpha$/Fe] suggest that the thick-disk consists of two distinct components with different chemical compositions and kinematic properties. For the inner region, the metal-weak thick disk (MWTD) contributes significantly when [$\alpha$/Fe] $>+$0.2\,dex and [Fe/H] $<-$0.8\,dex, while the canonical thick-disk (TD) dominates when [Fe/H]$>-$0.8\,dex. However, MWTD clearly appears only when [$\alpha$/Fe] $>+$0.2\,dex and [Fe/H] $<-$1.2\,dex for the outer region, and its proportion is lower than that of the inner region within the same metallicity. Similar result can be obtained from the {\VP} distribution. The high fraction of MWTD in the inner than that in the outer region imply that MWTD may form in the inner disk, and is an observational evidence about the inside-out disk formation scenario.

Optical and infrared polarization mapping and recent Planck observations of the filamentary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to the long axis of the cloud. We use the HAWC+ polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of $36^\circ$, about five times that measured on a larger scale by Planck. Applying the Davis-Chandrasekhar-Fermi (DCF) method with velocity information obtained from IRAM 30m C$^{18}$O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfv\'enic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfv\'enic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments.

The radiative properties of interstellar dust are affected not only by the grain size distribution but also by the grain porosity. We develop a model for the evolution of size-dependent grain porosity and grain size distribution over the entire history of galaxy evolution. We include stellar dust production, supernova dust destruction, shattering, coagulation, and accretion. Coagulation is {assumed to be} the source of grain porosity. We use a one-zone model with a constant dense gas fraction ($\eta_\mathrm{dense}$), which regulates the balance between shattering and coagulation. We find that porosity develops after small grains are sufficiently created by the interplay between shattering and accretion (at age $t\sim 1$ Gyr for star formation time-scale $\tau_\mathrm{SF}=5$ Gyr) and are coagulated. The filling factor drops down to 0.3 at grain radii $\sim 0.03~\mu$m for $\eta_\mathrm{dense}=0.5$. The grains are more porous for smaller $\eta_\mathrm{dense}$ because small grains, from which porous coagulated grains form, are more abundant. We also calculate the extinction curves based on the above results. The porosity steepens the extinction curve significantly for silicate, but not much for amorphous carbon. The porosity also increases the collisional cross-sections and produces slightly more large grains through the enhanced coagulation; however, the extinction curve does not necessarily become flatter because of the steepening effect by porosity. We also discuss the implication of our results for the Milky Way extinction curve.

We show that rotating white dwarfs admixed with dark matter have interesting properties that may be used to reveal the presence of dark matter. Even though such objects follow universal relations among the $I$ (moment of inertia), Love (tidal Love number), and $Q$ (quadrupole moment) that are robust with respect to different normal matter equations of state, these relations are sensitive to the dark matter fraction. Since each white dwarf may have a different dark matter fraction, the $I-\text{Love}-Q$ relations for dark matter admixed white dwarfs span bands above those without dark matter admixture. Furthermore, the limiting mass of dark matter admixed rotating white dwarfs can be increased beyond those without any dark matter for some rotational rules. Ultra-massive white dwarfs with a total mass of at least $2.6$ $M_{\odot}$ could be formed. The accretion-induced collapse of such an object may lead to a $2.6$ $M_{\odot}$ compact object, such as the one discovered in the gravitational-wave event GW190814.

The climates of terrestrial planets with a small amount of water on their surface, called land planets, are significantly different from the climates of planets having a large amount of surface water. Land planets have a higher runaway greenhouse threshold than aqua planets, which extends the inner edge of the habitable zone inward. Land planets also have the advantage of avoiding global freezing due to drier tropics, leading to a lower planetary albedo. In this study, we systematically investigate the complete freezing limit for various surface water distribution using a three-dimensional dynamic atmospheric model. As in a previous study, we found that a land planet climate has dry tropics that result in less snow and fewer clouds. The complete freezing limit decreases from that for aqua planets (92% S0, where S0 is Earth's present insolation) to that for land planets (77% S0) with an increasing dry area. Values for the complete freezing limit for zonally uniform surface water distributions are consistently lower than those for meridionally uniform surface water distribution. This is because the surface water distribution in the tropics in the meridionally uniform cases causes ice-albedo feedback until a planet lapses into the complete freezing state. For a surface water distribution using the topographies of the terrestrial planets, the complete freezing limit has values near those for the meridionally uniform cases. Our results indicate that the water distribution is important for the onset of a global ice-covered state for Earth-like exoplanets.

We present the results of 3D hydrodynamic simulations of gamma-ray burst (GRB) jet emanating from a massive star with a particular focus on the formation of high-velocity quasi-spherical ejecta and the jet-induced chemical mixing. Recent early-time optical observations of supernovae associated with GRBs (e.g., GRB 171205A/SN 2017iuk) indicate a considerable amount of heavy metals in the high-velocity outer layers of the ejecta. Using our jet simulations, we show that the density and chemical structure of the outer ejecta implied by observations can be naturally reproduced by a powerful jet penetrating the progenitor star. We consider three representative jet models with a stripped massive star, a standard jet, a weak jet, and a jet choked by an extended circumstellar medium, to clarify the differences in the dynamical evolution and the chemical properties of the ejected materials. The standard jet successfully penetrates the progenitor star and creates a quasi-spherical ejecta component (cocoon). The jet-induced mixing significantly contaminates the cocoon with heavy elements that have been otherwise embedded in the inner layer of the ejecta. The weak and choked jet models fail to produce an ultra-relativistic jet but produce a quasi-spherical cocoon with different chemical properties. We discuss the impact of the different jet-star interactions on the expected early-time electromagnetic signatures of long GRBs and how to probe the jet dynamics from observations.

Power spectra of spatial fluctuations of X-ray emission may impose constraints on the origins of the emission independent of that from the energy spectra. We generated spatial power spectrum densities (PSD) of blank X-ray skies observed with Suzaku X-ray observatory utilizing the modified $\Delta$-variance method. Using the total measured count rate as the diagnostic tool, we found that a model consisting of the sum of two components, one for the unresolved faint point sources and one for the uniform flat-field emission, can well represent the observed PSD in three different energy bands (0.2-0.5 keV, 0.5-2 keV, and 2-10 keV); only an upper limit is obtained for the latter component in 2-10 keV. X-ray counting rates corresponding to the best-fit PSD model functions and diffuse emission fractions were estimated, and we confirmed that the sum of the counting rates of two model components is consistent with those actually observed with the detector for all energy bands. The ratio of the flat-field counting rate to the total in 0.5-2 keV, however, is significantly larger than the diffuse emission fraction estimated from the model fits of energy spectra. We discussed that this discrepancy can be reconciled by systematic effects in the PSD and energy spectrum analyses. The present study demonstrates that the spatial power spectrum is powerful in constraining the origins of the X-ray emission.

We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.

We propose that fast radio bursts (FRBs) can be used as the probes to constrain the possible anisotropic distribution of baryon matter in the Universe. Monte Carlo simulations show that, 400 (800) FRBs are enough to detect the anisotropy at 95\% (99\%) confidence level, if the dipole amplitude is at the order of magnitude 0.01. However, much more FRBs are required to tightly constrain the dipole direction. Even 1000 FRBs are far from enough to constrain the dipole direction within angular uncertainty $\Delta\theta<40^{\circ}$ at 95\% confidence level. The uncertainty on the dispersion measure of host galaxy does not significantly affect the results. If the dipole amplitude is in the level of 0.001, however, 1000 FRBs are not enough to correctly detect the anisotropic signal.

As a space X-ray imaging mission dedicated to time-domain astrophysics, the Einstein Probe (EP) carries two kinds of scientific payloads, the wide-field X-ray telescope (WXT) and the follow-up X-ray telescope (FXT). FXT utilizes Wolter-I type mirrors and the pn-CCD detectors. In this work, we investigate the in-orbit background of FXT based on Geant4 simulation. The impact of various space components present in the EP orbital environment are considered, such as the cosmic photon background, cosmic ray primary and secondary particles (e.g. protons, electrons and positrons), albedo gamma rays, and the low-energy protons near the geomagnetic equator. The obtained instrumental background at 0.5-10 keV, which is mainly induced by cosmic ray protons and cosmic photon background, corresponds to a level of $\sim$3.1$\times$10$^{-2}$ counts s$^{-1}$ keV$^{-1}$ in the imaging area of the focal plane detector (FPD), i.e. 3.7$\times$10$^{-3}$ counts s$^{-1}$ keV$^{-1}$ cm$^{-2}$ after normalization. Compared with the instrumental background, the field of view (FOV) background, which is induced by cosmic photons reflected by the optical mirror, dominates below 2 keV. Based on the simulated background level within the focal spot (a 30$^{\prime\prime}$-radius circle), the sensitivity of FXT is calculated, which could theoretically achieve several $\mu$crab (in the order of 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$) in 0.5-2 keV and several tens of $\mu$crab (in the order of 10$^{-13}$ erg cm$^{-2}$ s$^{-1}$) in 2-10 keV for a pointed observation with an exposure of 25 minutes. This sensitivity becomes worse by a factor of $\sim2$ if additional 10% systematic uncertainty of the background subtraction is included.

Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting $r$-process nucleosynthesis and kilonova emission. In this work, we analyze neutrino emission from a large sample of merger radiation hydrodynamics simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $10^{53}{\rm erg~s^{-1}}$) decreases for increasing reduced tidal deformability, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterization of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

In the present work, recent characterization results of the 4K$\times$4K CCD Imager (a first light instrument of the 3.6m Devasthal Optical Telescope; DOT) and photometric calibrations are discussed along with measurements of the extinction coefficients and sky brightness values at the location of the 3.6m DOT site based on the imaging data taken between 2016 to 2021. For the 4K$\times$4K CCD Imager, all given combinations of gains (1, 2, 3, 5, and 10 e$^-$/ADU) and readout noise values for the three readout speeds (100 kHz, 500 kHz, and 1 MHz) are verified using the sky flats and bias frames taken during early 2021; measured values resemble well with the theoretical ones. Using colour-colour and colour-magnitude transformation equations, colour coefficients ($\alpha$) and zero-points ($\beta$) are determined to constrain and examine their long-term consistencies and any possible evolution based on $UBVRI$ observations of several Landolt standard fields observed during 2016-2021. Our present analysis exhibits consistency among estimated $\alpha$ values within the 1$\sigma$ and does not show any noticeable trend with time. We also found that the photometric errors and limiting magnitudes computed using the data taken using the CCD Imager follow the simulated ones published earlier. The average extinction coefficients, their seasonal variations, and zenith night-sky brightness values for the moon-less nights for all ten Bessell and SDSS filters are also estimated and found comparable to those reported for other good astronomical sites.

In this work, we present a catalog of cataclysmic variables (CVs) identified from the Sixth Data Release (DR6) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). To single out the CV spectra, we introduce a novel machine-learning algorithm called UMAP to screen out a total of 169,509 H$\alpha$-emission spectra, and obtain a classification accuracy of the algorithm of over 99.6$\%$ from the cross-validation set. We then apply the template matching program PyHammer v2.0 to the LAMOST spectra to obtain the optimal spectral type with metallicity, which helps us identify the chromospherically active stars and potential binary stars from the 169,509 spectra. After visually inspecting all the spectra, we identify 323 CV candidates from the LAMOST database, among them 52 objects are new. We further discuss the new CV candidates in subtypes based on their spectral features, including five DN subtype during outbursts, five NL subtype and four magnetic CVs (three AM Her type and one IP type). We also find two CVs that have been previously identified by photometry, and confirm their previous classification by the LAMOST spectra.

Ultra-hot jupiters (UHJs) are giant planets on short orbital periods with high equilibrium temperature (T_eq) values. Their hot, extended atmospheres are perfect laboratories for transmission spectroscopy studies based on high-resolution spectrographs. In recent years, a variety of atoms and molecules were found in their atmospheres, using different methods such as cross-correlation or transmission and emission spectroscopy. Here, we present the studies of six ultra-hot Jupiters: HAT-P-57b, KELT-7b, KELT-17b, KELT-21b, MASCARA-1b, and WASP-189b, based on high-resolution observations obtained with HARPS-N and HARPS spectrographs. By applying line and cross-correlation transmission spectroscopy methods, we searched for the absorption features of a broad range of atomic and molecular species. We did not detect any absorption features in our sample of UHJs, with the exception of WASP-189b, for which we detected FeI, FeII, and TiI using cross-correlation. The transmission spectroscopy of single lines for WASP-189b revealed several absorption features (including Halpha, Hbeta and Ca H&K), but they remain tentative pending a better modelling of the gravity darkening deformation of the Rossiter-McLaughlin effect. The non-detections with regard to the rest of the planets can be explained via a combination of stellar pulsations and the Rossiter-McLaughlin effect, which mask possible planetary signals for most of the planets, and by the low signal-to-noise ratios (S/N) of the observations for KELT-21b. Here, we compare our results with the known population of planets for which atmospheric detections have been reported in the literature. We find that the empirical frontier between hot and ultra-hot planets, based on the detection of atomic and ionized species in their atmospheres, can be established as Teq = 2150 K.

The second catalogue of Planck Sunyaev-Zeldovich (SZ) sources, hereafter PSZ2, represents the largest galaxy cluster sample selected by means of their SZ signature in a full-sky survey. Using telescopes at the Canary Island observatories, we conducted the long-term observational program 128- MULTIPLE-16/15B (hereafter LP15), a large and complete optical follow-up campaign of all the unidentified PSZ2 sources in the northern sky, with declinations above $-15^\circ$ and no correspondence in the first Planck catalogue PSZ1. This paper is the third and last in the series of LP15 results, after Streblyanska et al. (2019) and Aguado-Barahona et al. (2019), and presents all the spectroscopic observations of the full program. We complement these LP15 spectroscopic results with Sloan Digital Sky Survey (SDSS) archival data and other observations from a previous program (ITP13-08), and present a catalog of 388 clusters and groups of galaxies including estimates of their velocity dispersion. The majority of them (356) are the optical counterpart of a PSZ2 source. A subset of 297 of those clusters is used to construct the $M_{\rm SZ}-M_{\rm dyn}$ scaling relation, based on the estimated SZ mass from Planck measurements and our dynamical mass estimates. We discuss and correct for different statistical and physical biases in the estimation of the masses, such as the Eddington bias when estimating $M_{SZ}$ and the aperture and the number of galaxies used to calculate $M_{dyn}$. The SZ-to-dynamical mass ratio for those 297 PSZ2 clusters is $(1-B) = 0.80\pm0.04$ (stat) $\pm 0.05$ (sys), with only marginal evidence for a possible mass dependence of this factor. Our value is consistent with previous results in the literature, but presents a significantly smaller uncertainty due to the use of the largest sample size for this type of studies.

Recently, a covariant formulation of non-equilibrium phenomena in the context of General Relativity was proposed in order to explain from first principles the observed accelerated expansion of the Universe, without the need for a cosmological constant, leading to the GREA theory. Here, we confront the GREA theory against the latest cosmological data, including type Ia supernovae, baryon acoustic oscillations, the cosmic microwave background (CMB) radiation, Hubble rate data from the cosmic chronometers and the recent $H_0$ measurements. We perform Markov Chain Monte Carlo analyses and a Bayesian model comparison, by estimating the evidence via thermodynamic integration, and find that when all the aforementioned data are included, but no prior on $H_0$, the difference in the log-evidence is $\sim -9$ in favor of GREA, thus resulting in overwhelming support for the latter over the cosmological constant and cold dark matter model ($\Lambda$CDM). When we also include priors on $H_0$, either from Cepheids or the Tip of the Red Giant Branch measurements, then due to the tensions with CMB data the GREA theory is found to be statistically equivalent with $\Lambda$CDM.

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

We introduce in this work the Solar System Object Search Service (SSOSS), a service aimed at providing the scientific community with a search service for all potential detections of SSOs among the ESA astronomy archival imaging data. We illustrate its functionalities using the case of asteroid (16) Psyche, for which no information in the far-IR (70-500 {\mu}m) has previously been reported, to derive its thermal properties in preparation for the upcoming NASA Psyche mission. This service performs a geometrical cross-match of the orbital path of each object with respect to the public high-level imaging footprints stored in the ESA archives. For this first release, three missions were chosen: XMM-Newton, the Hubble Space Telescope (HST), and Herschel Observatory. We present a catalog listing all potential detections of asteroids within estimated limiting magnitude or flux limit in Herschel, XMM-Newton, and HST archival imaging data, including 909 serendipitous detections in Herschel images, 985 in XMM-Newton Optical Monitor camera images, and over 32,000 potential serendipitous detections in HST images. We also present a case study: the analysis of the thermal properties of Psyche from four serendipitous Herschel detections, combined with previously published thermal IR measurements. We see strong evidence for an unusual drop in (hemispherical spectral) emissivity, from 0.9 at 100 {\mu}m down to about 0.6 at 350 {\mu}m, followed by a possible but not well-constrained increase towards 500 {\mu}m, comparable to what was found for Vesta. The combined thermal data set puts a strong constraint on Psyche's thermal inertia (between 20 to 80$J m^{-2} s^{-1/2} K^{-1}$) and favours an intermediate to low level surface roughness (below 0.4 for the rms of surface slopes).

Neutrino telescopes are unrivaled tools to explore the Universe at its most extreme. The current generation of telescopes has shown that very high energy neutrinos are produced in the cosmos, even with hints of their possible origin, and that these neutrinos can be used to probe our understanding of particle physics at otherwise inaccessible regimes. The fluxes, however, are low, which means newer, larger telescopes are needed. Here we present the Pacific Ocean Neutrino Experiment, a proposal to build a multi-cubic-kilometer neutrino telescope off the coast of Canada. The idea builds on the experience accumulated by previous sea-water missions, and the technical expertise of Ocean Networks Canada that would facilitate deploying such a large infrastructure. The design and physics potential of the first stage and a full-scale P-ONE are discussed.

Warm dark matter (WDM) can potentially explain small-scale observations that currently challenge the cold dark matter (CDM) model, as warm particles suppress structure formation due to free-streaming effects. Observing small-scale matter distribution provides a valuable way to distinguish between CDM and WDM. In this work, we use observations from the Dark Energy Survey and PanSTARRS1, which observe 270 Milky-Way satellites after completeness corrections. We test WDM models by comparing the number of satellites in the Milky Way with predictions derived from the Semi-Analytical SubHalo Inference ModelIng (SASHIMI) code, which we develop based on the extended Press-Schechter formalism and subhalos' tidal evolution prescription. We robustly rule out WDM with masses lighter than 4.4 keV at 95% confidence level for the Milky-Way halo mass of $10^{12} M_\odot$. The limits are a weak function of the (yet uncertain) Milky-Way halo mass, and vary as $m_{\rm WDM}>3.6$-$5.1$ keV for $(0.6$-$2.0) \times 10^{12} M_\odot$. For the sterile neutrinos that form a subclass of WDM, we obtain the constraints of $m_{\nu_s}>11.6$ keV for the Milky-Way halo mass of $10^{12} M_{\odot}$. These results based on SASHIMI do not rely on any assumptions of galaxy formation physics or are not limited by numerical resolution. The models, therefore, offer a robust and fast way to constrain the WDM models. By applying a satellite forming condition, however, we can rule out the WDM mass lighter than 9.0 keV for the Milky-Way halo mass of $10^{12} M_\odot$.

In binary systems, the helium accretion onto carbon-oxygen (CO) white dwarfs (WDs) plays a vital role in many astrophysical scenarios, especially in supernovae type Ia. Moreover, ignition density for accretion rate $\dot{M} \lesssim 10^{-9} M_{\odot}$ $yr^{-1}$ in helium accreting CO white dwarfs decides the triggering mechanism of a supernova explosion which could be either off-centre helium flash or central carbon flash. We aim to study the accretion of helium with a slow accretion rate $5 \times 10^{-10} M_{\odot}$ $yr^{-1}$ onto relatively cool and hot white dwarfs of different abundances of carbon and oxygen. The simulation code Modules for Experiments in Stellar Astrophysics (MESA) has been used for our study. We analyze the variation in several properties like surface gravity (g), helium luminosity ($L_{He}$), and effective temperature ($T_{eff}$) during the accretion phase of the white dwarfs. We also calculate the ignition density (${\rho}_{He}$) and ignition temperature ($T_{He}$) of helium burning. As expected, the size of WD decreases and g increases during the accretion. However, a red-giant-like expansion is observed after the rapid ignition towards the end. The dependence of helium accreting WD evolution on its composition has also been explored in this study. We find that white dwarfs of the lower abundance of carbon accrete slightly longer before the onset of helium ignition.

The Southern Wide-field Gamma-ray Observatory (SWGO) Collaboration is currently engaged in design and prototyping work towards the realisation of this future gamma-ray facility. SWGO will complement CTA and the existing ground-particle based-detectors of the Northern Hemisphere (HAWC and LHAASO) with a very wide field and high duty cycle view of the southern sky. Here I summarise the status of the project and plans for the future, including expectations for sensitivity and science targets as well as the status of the site search and technological developments.

The formation of compact objects--neutron stars, black holes, and supermassive black holes--and its connection to the chemical composition of the galaxies is one of the central questions in astrophysics. We propose a novel data-driven, multi-messenger technique to address this question by exploiting the inevitable correlation between gravitational waves and atomic/molecular emission line signals. We show that the minimum delay time of $0.5$ Gyr and the power-law index $\kappa=1$ of the fiducial scenario of the probability distribution function $p(t_d)\propto t_d^{-\kappa}$ is possible to measure with a standard deviation $0.12$ (and $0.45$) and $0.06$ (and $0.34$) respectively from five years of LIGO-Virgo-KAGRA observation in synergy with SPHEREx line intensity mapping (and DESI emission-line galaxies).

We investigate the failed partial eruption of a filament system in NOAA AR 12104 on 2014 July 5, using multiwavelength EUV, magnetogram, and H$\alpha$ observations, as well as magnetic field modeling. The filament system consists of two almost co-spatial segments with different end points, both resembling a C shape. Following an ejection and a precursor flare related to flux cancellation, only the upper segment rises and then displays a prominent twisted structure, while rolling over toward its footpoints. The lower segment remains undisturbed, indicating that the system possesses a double-decker structure. The erupted segment ends up with a reverse-C shape, with material draining toward its footpoints, while losing its twist. Using the flux rope insertion method, we construct a model of the source region that qualitatively reproduces key elements of the observed evolution. At the eruption onset, the model consists of a flux rope atop a flux bundle with negligible twist, which is consistent with the observational interpretation that the filament possesses a double-decker structure. The flux rope reaches the critical height of the torus instability during its initial relaxation, while the lower flux bundle remains in stable equilibrium. The eruption terminates when the flux rope reaches a dome-shaped quasi-separatrix layer that is reminiscent of a magnetic fan surface, although no magnetic null is found. The flux rope is destroyed by reconnection with the confining overlying flux above the dome, transferring its twist in the process.

This paper is the second part of a series of studies discussing a novel attitude determination method for nano-satellites. Our approach is based on the utilization of thermal imaging sensors to determine the direction of the Sun and the nadir with respect to the satellite with sub-degree accuracy. The proposed method is planned to be applied during the Cubesats Applied for MEasuring and LOcalising Transients (CAMELOT) mission aimed at detecting and localizing gamma-ray bursts with an efficiency and accuracy comparable to large gamma-ray space observatories. In this paper we introduce a simulation model aimed at testing the applicability of our attitude determination approach. Its first part simulates the orbit and rotation of a satellite with arbitrary initial conditions while its second part applies our attitude determination algorithm which is based on a multiplicative extended Kalman filter. The simulated satellite is assumed to be equipped with a GPS system, MEMS gyroscopes and the infrasensors. These instruments provide the required data input for the Kalman filter. We demonstrate the applicability of our attitude determination algorithm by simulating the motion of a nano-satellite on Low Earth Orbit. Our results show that the attitude determination may have a 1$\sigma$ error of $\sim30'$ even with a large gyroscope drift during the orbital periods when the infrasensors provide both the direction of the Sun and the Earth (the nadir). This accuracy is an improvement on the point source detection accuracy of the infrasensors. However, the attitude determination error can get as high as 25$^{\circ}$ during periods when the Sun is occulted by the Earth. We show that following an occultation period the attitude information is immediately recovered by the Kalman filter once the Sun is observed again.

We have constructed a catalog of active galactic nuclei (AGNs) with $z < 0.13$, based on optical spectroscopy, from the parent sample of galaxies in the 6dF galaxy survey (Final Release of 6dFGS), a census of the Southern hemisphere. This work is an extension of our all sky AGN catalog \citet [ZCF, here after]{ZCF19}. The ZCF is based on 43,533 galaxies with $\rm K_s \leq$ 11.75 ($z \leq 0.09$) in the 2MASS Redshift Survey (2MRS). The parent catalog of this work, the 6dF catalog, consists of 136,304 publicly available digital spectra for 125,071 galaxies with $\rm Dec \leq 0^\circ$ and $\rm K_s \leq 12.65$ (median $z = 0.053$). Our AGN catalog consists of 3109 broad line AGNs and 12,156 narrow line AGNs which satisfy the \citet{Kauffmann03} criteria, of which 3865 also satisfy the \citet{Kewley01} criteria. We also provide emission line widths, fluxes, flux errors, and signal-to-noise ratios of all the galaxies in our spectroscopic sample, allowing users to customize the selection criteria. In addition, we provide AGN likelihood for the rest of galaxies based on the availability and quality of their spectra. These likelihood values can be used for rigorous statistical analyses.

We report the detection of radio emission from the known X-ray flaring star EXO 040830$-$7134.7 during MeerKAT observations of the nearby cataclysmic variable VW Hydri. We have three epochs of MeerKAT observations, where the star is not detected in the first epoch, is detected in the second epoch, and is marginally detected in the third epoch. We cannot distinguish whether the detection is quiescent emission or a transient radio burst. If we assume the radio detection is quiescent emission the source lies somewhat to the right of the G\"udel-Benz relation; however, if we assume the upper-limit on the radio non-detection is indicative of the quiescent emission then the source lies directly on the relation. Both cases are broadly consistent with the relation. We use archival spectral energy distribution data and new SALT high-resolution spectroscopy to confirm that EXO 040830$-$7134.7 is a chromospherically active M-dwarf with a temperature of 4000$\pm$200 K of spectral type M0V. We use ASAS, ASAS-SN and TESS optical photometry to derive an improved rotational period of 5.18$\pm$0.04 days. This is the first radio detection of the source, and the first MeerKAT detection of an M-dwarf.

We study an impact of asymmetric fermionic dark matter on neutron star properties, including tidal deformability, mass, radius, etc. We present the conditions at which dark matter particles tend to form a compact structure in a core of the star or create an extended halo around it. We show that compact core of dark matter leads to a decrease of the total gravitational mass and tidal deformability compared to a pure baryonic star, while presence of a dark matter halo increases those observable quantities. By imposing an existing astrophysical and gravitational wave constraints set by LIGO/Virgo Collaboration together with the recent results on the spatial distribution of dark matter in the Milky Way we determine a new upper limit on the mass and fraction of dark matter particles inside compact stars. Furthermore, we show that the formation of an extended halo around a NS is incompatible with the GW170817 tidal deformability constraint.

Single dish observations of AGB star L$_2$ Pup have revealed exceptionally low mass loss rate and expansion velocity, challenging interpretations in terms of standard wind models. Recent VLT and ALMA observations have drawn a detailed picture of the circumstellar envelope within $\sim$20 au from the centre of the star: a nearly edge-on rotating disc of gas and dust, probably hosting a planetary companion near the star. However, these observations provide no direct information on the wind escaping the gravity of the star. The present article uses ALMA observations of the $^{12,13}$CO(3-2), $^{29}$SiO(8-7), $^{12}$CO(2-1) and $^{28}$SiO(5-4) line emissions to shed new light on this issue. It shows the apparent normality of L$_2$ Pup in terms of the formation of the nascent wind, with important line broadening within 4 au from the centre of the star, but no evidence for a wind flowing along the disc axis. At larger distances, up to some 200 au from the centre of the star, the wind morpho-kinematics is dominated by a disc, or equatorial enhancement, expanding isotropically and radially with a velocity not exceeding some 5 km s$^{-1}$ , inclined in the north-west/south-east direction with respect to the plane of the sky. In addition, outflows of lower density are observed on both sides of the disc, covering large solid angles about the disc axis, contributing about half the flux of the disc. Such morphology is at strong variance with the expectation of a pair of back-to-back outflows collimated by the central gas-and-dust disc.

The observed late-type WC Wolf-Rayet stars (WC7-9) with low luminosity below $\rm \log L/L_{\odot} < 5.4$ in the HR diagram cannot be reproduced satisfactorily by the evolutionary track of single stars. The mass transfer due to Roche lobe overflow drastically modifies the internal structure and surface compositions of two components. Therefore, binaries provide a very promising evolutionary channel to produce these WC stars.

The origin of the variable X-ray emission in the $(0.3-30)\,$keV energy range of ultraluminous X-ray sources (ULXs) remains unclear, making it difficult to constrain the mass of the central compact object. X-ray luminosities of bright ULXs can be explained with sub-critical accretion ($L<L_{\rm Edd}$) on to an intermediate-mass BH, with the alternative being super-critical accretion on to a stellar-mass BH. Broadband X-ray emission in the former scenario can be explained using the canonical disk plus Comptonizing corona model, whereas in the latter scenario radiation pressure driven massive winds lead to complex spectra that are inclination angle dependent. Here we fit the broadband (optical/UV to X-ray) spectrum of the persistently bright ULX Holmberg IX X-1 with the disk-corona plus irradiated outer disk model in an effort to constrain the BH mass. We use a one-zone time-dependent numerical code to exactly solve for the steady-state properties of the optically thick coronal photon-electron-positron plasma. Our modelling suggests that Holmberg IX X-1 hosts a stellar mass BH, with mass $4\lesssim(\hat M_{\rm BH}\equiv\alpha M_{\rm BH}/M_\odot)\lesssim10$ where $1/6\leq\alpha<1$ for a spinning (Kerr) BH, undergoing super-critical accretion ($L_{\rm Bol}/L_{\rm Edd}\sim20\alpha$). In our model, the X-ray spectrum below $10\,$keV is explained with an absorbed multi-colour disk spectrum having inner disk temperature $k_BT_{\rm in}\sim(2.2-2.9)\,$keV. An additional cooler thermal spectral component, as found in many works and not included in our modeling, is required. The hard excess above $10\,$keV, as seen by NuSTAR, arises in a photon-rich optically-thick Comptonizing spherical corona with optical depth $\tau_T\sim3.5$ and particle temperature $k_BT_e\sim14\,$keV.

We have observed the 70 $\mu$m dark infrared dark cloud (IRDC) G14.492-00.139 in the N$_2$D$^+$ $J$=3--2, DCO$^+$ $J$=3--2, DCN $J$=3--2, and C$^{18}$O $J$=2--1 lines, using the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). We find that the spatial distribution is different among the observed emission from the deuterated molecular lines. The N$_2$D$^+$ emission traces relatively quiescent regions, while both the DCO$^+$ and DCN emission emanates mainly from regions with signs of active star-formation. In addition, the DCO$^+$/N$_2$D$^+$ ratio is found to be lower in several dense cores than in starless cores embedded in low-mass star-forming regions. By comparing the observational results with chemical model calculations, we discuss the origin of the low DCO$^+$/N$_2$D$^+$ ratio in this IRDC clump. The low DCO$^+$/N$_2$D$^+$ ratio can be explained if the temperature of the dense cores is in the range between the sublimation temperature of N$_2$ ($\sim$20 K) and CO ($\sim$25 K). The results suggest that the dense cores in G14.492-00.139 are warmer and denser than the dense cores in low-mass star-forming regions.

Geometric modelling of Coronal Mass Ejections (CMEs) is a widely used tool for assessing their kinematic evolution. Furthermore, techniques based on geometric modelling, such as ELEvoHI, are being developed into forecast tools for space weather prediction. These models assume that solar wind structure does not affect the evolution of the CME, which is an unquantified source of uncertainty. We use a large number of Cone CME simulations with the HUXt solar wind model to quantify the scale of uncertainty introduced into geometric modelling and the ELEvoHI CME arrival times by solar wind structure. We produce a database of simulations, representing an average, a fast, and an extreme CME scenario, each independently propagating through 100 different ambient solar wind environments. Synthetic heliospheric imager observations of these simulations are then used with a range of geometric models to estimate the CME kinematics. The errors of geometric modelling depend on the location of the observer, but do not seem to depend on the CME scenario. In general, geometric models are biased towards predicting CME apex distances that are larger than the true value. For these CME scenarios, geometric modelling errors are minimised for an observer in the L5 region. Furthermore, geometric modelling errors increase with the level of solar wind structure in the path of the CME. The ELEvoHI arrival time errors are minimised for an observer in the L5 region, with mean absolute arrival time errors of $8.2\pm1.2$~h, $8.3\pm1.0$~h, and $5.8\pm0.9$~h for the average, fast, and extreme CME scenarios.

The Trappist-1 system contains seven roughly Earth-sized planets locked in a multi-resonant orbital configuration, which has enabled precise measurements of the planets' masses and constrained their compositions. Here we use the system's fragile orbital structure to place robust upper limits on the planets' bombardment histories. We use N-body simulations to show how perturbations from additional objects can break the multi-resonant configuration by either triggering dynamical instability or simply removing the planets from resonance. The planets cannot have interacted with more than ${\sim 5\%}$ of an Earth mass (${M_\oplus}$) in planetesimals -- or a single rogue planet more massive than Earth's Moon -- without disrupting their resonant orbital structure. This implies an upper limit of ${10^{-4}}$ to ${10^{-2} M_\oplus}$ of late accretion on each planet since the dispersal of the system's gaseous disk. This is comparable to or less than the late accretion on Earth after the Moon-forming impact, and demonstrates that the Trappist-1 planets' growth was complete in just a few million years, roughly an order of magnitude faster than Earth's. Our results imply that any large water reservoirs on the Trappist-1 planets must have been incorporated during their formation in the gaseous disk.

Type Ibn supernovae (SNe) are a mysterious class of transients whose spectra exhibit persistently narrow HeI lines, and whose bolometric light curves are typically fast evolving and overluminous at peak relative to standard Type Ibc SNe. We explore the interaction scenario of such Type Ibn SNe by performing radiation-hydrodynamics and radiative-transfer calculations. We find that standard-energy helium-star explosions within dense wind-like circumstellar material (CSM) can reach on day timescales a peak luminosity of a few 10^44 erg/s, reminiscent of exceptional events like AT2018cow. Similar interactions but with weaker winds can lead to Type Ibc SNe with double-peak light curves and peak luminosities in the range ~10^42.2 to ~10^43 erg/s. In contrast, the narrow spectral lines and modest peak luminosities of most Type Ibn SNe are suggestive of a low-energy explosion in an initially <~5 Msun helium star, most likely arising from interacting binaries, and colliding with a massive helium-rich, probably ejecta-like, CSM at ~10^15 cm. Nonlocal thermodynamic equilibrium radiative-transfer simulations of a slow-moving dense shell born out and powered by the interaction compare favorably to Type Ibn SNe like 2006jc, 2011hw, or 2018bcc at late times and suggest a composition made of about 50% helium, a solar metallicity, and a total ejecta/CSM mass of 1-2Msun. A lower fractional helium abundance leads to weak or absent HeI lines and thus excludes more massive configurations for observed Type Ibn SNe. Further, the dominance of FeII emission below 5500A seen in Type Ibn SNe at late times is not predicted at low metallicity. Hence, despite their promising properties, Type Ibn SNe from pulsational-pair instability in very massive stars, which requires low metallicity, have probably not yet been observed.

Our paper presents the results of the second census of pulsars in decametre wave range at UTR-2 radio telescope. Over the past ten years, the number of discovered nearby pulsars in the world has doubled, which has made it urgent to search for a low-frequency radio emission from newly discovered sources. To increase this census sensitivity, the integration time was doubled compared with the first census of 2010-2013. As a result, the decametre radio emission of 20 pulsars was detected, their flux densities and the shape of pulses were obtained. The dispersion measure for 10 pulsars and the rotation period for 8 pulsars were refined. For several pulsars the scattering time constant and FWHM were estimated in decametre wave range. Upper limits of flux densities of 102 not yet detected pulsars were also estimated.

The association of a short gamma-ray burst with a core-collapse supernova seems to challenge current scenarios for the origin of these extreme events. But how much can we rely on observed duration for pinpointing their progenitors?

The observed durations of prompt gamma-ray emission from Gamma-Ray Bursts (GRBs) are often used to infer the progenitors and energetics of the sources. Inaccurate duration measurements will have a significant impact on constraining the processes powering the bursts. The "tip-of-the-iceberg" effect describes how the observed signal is lost into background noise; lower instrument sensitivity leads to higher measurement bias. In this study, we investigate how observing conditions, such as the number of enabled detectors, background level, and incident angle of the source relative to the detector plane affect the measured duration of GRB prompt emission observed with the Burst Alert Telescope on board the Neil Gehrels Swift Observatory (Swift/BAT). We generate "simple-pulse" light curve from an analytical Fast Rise Exponential Decay (FRED) function and from a sample of eight real GRB light curves, we fold these through the Swift/BAT instrument response function to simulate light curves Swift/BAT would have observed for specific observing conditions. We find duration measurements are highly sensitive to observing conditions and the incident angle of the source has the highest impact on measurement bias. In most cases duration measurements of synthetic light curves are significantly shorter than the true burst duration. For the majority of our sample, the percentage of duration measurements consistent with the true duration is as low as $\sim10\%-30\%$. In this article, we provide quantification of the tip-of-the-iceberg effect on GRB light curves due to Swift/BAT instrumental effects for several unique light curves.

Dust emission mechanisms as one aspect of wind-driven particle motion on planetary surfaces are still poorly understood. The microphysics is important though as it determines dust sizes and morphologies which set sedimentation speeds and optical properties. We consider the effects of tribocharging in this context as grains in wind driven granular matter charge significantly. This leads to large electric fields above the granular bed. Airborne dielectric grains are polarized in these electric fields, which leads to attractive forces between grains. To simulate aggregation under these conditions we carried out drop tower experiments using tracer particles, mimicking the gas coupling behavior of small dust grains in terms of high surface to mass ratios and efficient gas drag. Under microgravity, the particles are released into an observation chamber in which an alternating electric field up to 80 kV/m is applied. Without electric field no aggregation can be observed on timescales of seconds. However, polarization instantly leads to aggregation of particles when the field is switched on and long chains aligned to the electric field form. Scaled to dust entrained into planetary atmospheres, fine and coarse grain fractions might readily form aggregates after being liberated. Under certain natural conditions, aggregates might therefore start chain-like or at least a chain-like appearance is favored. If atmospheric influences on their stability are small, aerodynamic and optical properties might depend on this.

The Galactic transient black hole candidate GX 339-4 is a very interesting object to study as it showed both complete and failed types of outbursts. We studied both spectral and temporal properties of the 2017-18 outburst of the source using archival data of NICER and AstroSat instruments. This 2017-18 outburst is found to be failed in nature, as during the entire period of the outburst, the source was only in the hard spectral state. Source spectra were highly dominated with the non-thermal fluxes. When we tried to fit spectra with phenomenological models, most of the spectra were fitted with only the powerlaw model, and only six spectra required disk black body plus powerlaw models. While fitting spectra with the physical two-component advective flow (TCAF) model, we observed that the flow was highly dominated by the sub-Keplerian halo rate. The presence of stronger shock at a larger radius from the black hole was also observed during the rising and declining phases of the outburst. A prominent signature of $0.31$~Hz QPO is observed with its four harmonics. Mass of the black hole was also estimated from our spectral analysis with the TCAF model as $11.24^{+0.59}_{-1.25}~M_\odot$.

Complete asymptotic expansions are developed for slow, Alfv\'en and fast magnetohydrodynamic waves at the base of an isothermal three-dimensional (3D) plane stratified atmosphere. Together with existing convergent Frobenius series solutions about $z=\infty$, matchings are numerically calculated that illuminate the fates of slow and Alfv\'en waves injected from below. An Alfv\'en wave in a two-dimensional model is 2.5D in the sense that the wave propagates in the plane of the magnetic field but its polarization is normal to it in an ignorable horizontal direction, and the wave remains an Alfv\'en wave throughout. The rotation of the plane of wave propagation away from the vertical plane of the magnetic field pushes the plasma displacement vector away from horizontal, thereby coupling it to stratification. It is shown that potent slow-Alfv\'en coupling occurs in such 3D models. It is found that about 50% of direction-averaged Alfv\'en wave flux generated in the low atmosphere at frequencies comparable to or greater than the acoustic cutoff can reach the top as Alfv\'en flux for small magnetic field inclinations $\theta$, and this increases to 80% or more with increasing $\theta$. On the other hand, direction-averaged slow waves can be 40% effective in converting to Alfv\'en waves at small inclination, but this reduces sharply with increasing $\theta$ and wave frequency. Together with previously explored fast-slow and fast-Alfv\'en couplings, this provides valuable insights into which injected transverse waves can reach the upper atmosphere as Alfv\'en waves, with implications for solar and stellar coronal heating and solar/stellar wind acceleration.

The PBHs clusters can be source of gravitational waves, and the merger rate depends on the spatial distribution of PBHs in the cluster which changes over time. It is well known that gravitational collisional systems experience the core collapse, that leads to significantly increase of the central density and shrinking of the core. After core collapse, the cluster expands almost self-similarly (i.e. density profile extends in size without changing its shape). These dynamic processes affect at the merger rate of PBHs. In this paper the dynamics of the PBH cluster is considered using the Fokker-Planck equation. We calculate the merger rate of PBHs on cosmic time scales and show that its time dependence has a unique signature. Namely, it grows by about an order of magnitude at the moment of core collapse which depends on the characteristic of the cluster, and then decreases according to the dependence $\mathcal{R} \propto t^{-1.48}$. It was obtained for monochromatic and power-law PBH mass distributions with some fixed parameters. Obtained results can be used to test the model of the PBH clusters via observation of gravitational waves at high redshift.

We present the results of NH3 (1,1), (2,2) and (3,3) and H2O maser simultaneous mapping observations toward the high-mass star-forming region W33 with the Nobeyama 45-m radio telescope. W33 has six dust clumps and one of which, W33 Main, is associated with a compact HII region. To investigate star-forming feedback activity on its surroundings, the spatial distribution of the physical parameters was established. The distribution of the rotational temperature shows a systematic change from west to east in our observed region. The high-temperature region obtained in the region near W33 Main is consistent with interaction between the compact HII region and the periphery molecular gas. The size of the interaction area is estimated to be approximately 1.25 pc. NH3 absorption features are detected toward the centre of the HII region. Interestingly, the absorption feature was detected only in the NH3 (1,1) and (2,2) transitions, with no absorption feature seen in the (3,3) transition. These complex profiles in NH3 are difficult to explain by a simple model and may suggest that the gas distribution around the HII region is highly complicated.

Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived from total intensity pulse profiles of pulsars via cross-correlation with an idealised `1D' template of that profile. While a number of results have been presented in the literature relying on the ever increasing sensitivity of such pulsar timing experiments, there is no consensus on the most reliable methods for TOA creation and, more importantly, the associated TOA uncertainties for each scheme. In this article, we present a comprehensive comparison of TOA determination practices, focusing on the creation of timing templates, TOA determination methods and the most useful TOA bandwidth. The aim is both to present a possible approach towards TOA optimisation as well as the (partial) identification of an optimal TOA-creation scheme and the demonstration of optimisation differences between pulsars and data sets. We compare the values of data-derived template profiles as compared to analytic profiles and evaluate the three most commonly used template-matching methods. Finally, we study the relation between timing precision and TOA bandwidth to identify any potential breaks in that relationship. As a practical demonstration, we apply our selected methods to European Pulsar Timing Array data on the three test pulsars PSRs\ J0218+4232, J1713+0747 and J2145$-$0750.

We present the results obtained from timing and spectral studies of the highly obscured low luminosity active galactic nucleus NGC~4941 using data obtained from the {\it Nuclear Spectroscopic Telescope Array} and the {\it Neil Gehrels Swift} Observatories. We find similar variability in $3-10$~keV and $10-60$~keV energy ranges with fractional rms variability of $\sim$14\%. We investigate broad-band spectral properties of the source in 3-150 keV range, using data from {\it NuSTAR} and {\it Swift}/BAT, with phenomenological slab model and physically motivated {\tt mytorus} model. From the spectral analysis, we find heavy obscuration with the global average column density of the obscured material as $3.09^{+1.68}_{-1.01} \times 10^{24}$ cm$^{-2}$. Evidence of a strong reflection component is observed in the spectrum. We detect a strong iron line with an equivalent width of $\sim$1~keV. From the slab model, we obtain the exponential cutoff energy as $177^{+92}_{-16}$~keV. From this, we estimate the Compton cloud properties with the hot electron temperature $kT_{\rm e} = 59^{+31}_{-5}$~keV and the optical depth $\tau = 2.7^{+0.2}_{-1.6}$.

We present results of new Chandra High-Energy Transmission Grating (HETG) observations (2019 November - December) of the massive Wolf-Rayet (WR) binary WR 48a. Analysis of these high-quality data showed that the spectral lines in this massive binary are broadened (FWHM = 1400 km/s) and marginally blushifted (~ -100 km/s). A direct modelling of these high-resolution spectra in the framework of the standard colliding stellar wind (CSW) picture provided a very good correspondence between the shape of the theoretical and observed spectra. Also, the theoretical line profiles are in most cases an acceptable representation of the observed ones. We applied the CSW model to the X-ray spectra of WR 48a from previous observations: Chandra-HETG (2012 October) and XMM-Newton (2008 January). From this expanded analysis, we find that the observed X-ray emission from WR48a is variable on the long timescale (years) and the same is valid for its intrinsic X-ray emission. This requires variable mass-loss rates over the binary orbital period. The X-ray absorption (in excess of that from the stellar winds in the binary) is variable as well. We note that lower intrinsic X-ray emission is accompanied by higher X-ray absorption. A qualitative explanation could be that the presence of clumpy and non-spherically symmetric stellar winds may play a role.

Europa has been spotted to have water outgassing activities by the space and ground-based telescopes as well as reanalysis of the Galileo data (Roth et al. 2014; Sparks et al. 2016, 2017; Paganini et al. 2020; Jia et al. 2018; Arnold et al. 2019). However, these observations only provided limited information about plume dynamics, which is critical in understanding the eruption mechanism and preparation of future exploration. We adopt a 3D DSMC model to investigate the plume characteristics of Europa assuming supersonic expansion originated from the undersurface vent. The main goal is to understand the physical processes and structures of Europa water vapor plumes, which can play a key role on probing its undersurface vent condition and outgassing mechanism. With a parametric study of the total gas production rate and initial gas bulk velocity, the gas number density, temperature and velocity information of the outgassing plumes from the various case studies are derived. Our results show that the plume gases experience acceleration through mutual collisions and adiabatic cooling when exiting and expanding from the surface. The central part of the plume with the relatively large gas production rates (of 1029 and 1030 H2O s-1) is found to sustain thermal equilibrium and nearly continuum condition. Column density maps integrated along two different viewing angles are presented to demonstrate the importance of the projection effect on remote sensing diagnostics. Finally, the density profiles at different altitudes are provided to prepare for observations of Europa plumes including the upcoming spacecraft missions such as JUICE and Europa Clipper.

We present the characterisation of the massive cluster ClG-J$104803.7+313843$ at $z=0.76$ performed using a serendipitous XMM-Newton observation. High redshift and massive objects represent an ideal laboratory to benchmark our understanding of how cluster form and assembly formation driven mainly by gravity.Leveraging the high throughput of XMM-Newton we were firstly able to determine the redshift of the object, shedding light on ambiguous photometric redshift associations. We investigated the morphology of this cluster which shows signs of merging activities in the outskirts and a flat core. We also measured the radial density profile up to $R_{500}$. With these quantities in hand, we were able to determine the mass, $M_{500}=5.64^{+0.79}_{-0.62} \times 10^{14}M_{\odot}$, using the YX proxy. This quantity improves previous measurement of the mass of this object by a factor of $\sim 3.5$. The characterisation of one cluster at such mass and redshift regime is fundamental as these objects are intrinsically rare, the number of objects discovered so far being less than $\sim 25$. Our study highlights the importance of using X-ray observations in combination with ancillary multi-wavelength data to improve our understanding of high-z and massive clusters

Hot, compact, hydrogen-deficient pre-white dwarfs (pre-WDs) with effective temperatures of Teff > 70,000 K and a surface gravity of 5.0 < log g < 7.0 are rather rare objects despite recent and ongoing surveys. It is believed that they are the outcome of either single star evolution (late helium-shell flash or late helium-core flash) or binary star evolution (double WD merger). Their study is interesting because the surface elemental abundances reflect the physics of thermonuclear flashes and merger events. Spectroscopically they are divided in three different classes, namely PG1159, O(He), or He-sdO. We present a spectroscopic analysis of five such stars that turned out to have atmospheric parameters in the range Teff = 70,000-80,000 K and log g = 5.2-6.3. The three investigated He-sdOs have a relatively high hydrogen mass fraction (10%) that is unexplained by both single (He core flash) and binary evolution (He-WD merger) scenarios. The O(He) star JL9 is probably a binary helium-WD merger, but its hydrogen content (6%) is also at odds with merger models. We found that RL 104 is the 'coolest' (Teff = 80,000 K) member of the PG1159 class in a pre-WD stage. Its optical spectrum is remarkable because it exhibits C IV lines involving Rydberg states with principal quantum numbers up to n = 22. Its rather low mass (0.48 +0.03/-0.02 Msun) is difficult to reconcile with the common evolutionary scenario for PG1159 stars due to it being the outcome of a (very) late He-shell flash. The same mass-problem faces a merger model of a close He-sdO plus CO WD binary that predicts PG1159-like abundances. Perhaps RL 104 originates from a very late He-shell flash in a CO/He WD formed by a merger of two low-mass He-WDs.

The High-Energy Lightcurve Generator (HILIGT) is a new web-based tool which allows the user to generate long-term lightcurves of X-ray sources. It provides historical data and calculates upper limits from image data in real-time. HILIGT utilizes data from twelve satellites, both modern missions such as XMM-Newton and Swift, and earlier facilities such as ROSAT, EXOSAT, Einstein or Ariel V. Together, this enables the user to query 50 years of X-ray data and, for instance, study outburst behavior of transient sources. In this paper we focus on the individual back-end servers for each satellite, detailing the software layout, database design, catalog calls, and image footprints. We compile all relevant calibration information of these missions and provide an in-depth summary of the details of X-ray astronomical instrumentation and data.

We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including discs and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, whereas in the spirals it is only 5-12%. J1010+0612 and J1430+1339 show double-peaked CO morphologies which do not have optical counterparts. Based on our analysis of the ionized and molecular kinematics and mm continuum emission, these CO morphologies could be produced by AGN feedback in the form of outflows and/or jets/shocks. The CO kinematics of the QSO2s are dominated by rotation but also reveal non-circular motions. According to our analysis of the kinematics, these non-circular motions correspond to molecular outflows mostly coplanar with the CO discs in four of the QSO2s, and either to a coplanar inflow or vertical outflow in the case of J1010+0612. These molecular outflows represent 0.2-0.7% of the QSO2s total molecular gas mass, have maximum velocities of 200-350 km/s, radii of 0.4-1.3 kpc and outflow rates of 8-16 Msun/yr. These properties are intermediate between those of the mild molecular outflows measured for Seyferts, and the fast and energetic outflows of ULIRGs. This suggests that it is not only AGN luminosity what drives massive molecular outflows. Other factors such as jet power, coupling between wind/jet/ionized outflows and CO discs, and amount/geometry of dense gas in the nuclear regions might be also relevant. Although we do not find evidence for a significant impact of the QSO2s molecular outflows/AGN feedback on the total molecular gas reservoirs and SFRs, they would be modifying the distribution of cold molecular gas in the central kpc of the galaxies.

Low angular momentum, general relativistic, axially symmetric accretion of hydrodynamic fluid onto Schwarzschild black holes may undergo more than one critical transition. To obtain the stationary integral solutions corresponding to such multi-critical accretion flow, one needs to employ numerical solutions of the corresponding fluid dynamics equations. In the present work, we develop a completely analytical solution scheme which may be used to find several trans-critical flow behaviours of aforementioned accretion, without explicitly solving the flow equations numerically. We study all possible geometric configurations of the flow profile, governed by all possible thermodynamic equations of state. We use Sturm's chain algorithm to find out how many physically acceptable critical points the accretion flow can have, and discuss the transition from the mono to the multi-critical flow profile, and related bifurcation phenomena. We thus illustrate, completely analytically, the application of certain aspects of the dynamical systems theory in the field of large scale astrophysical flow under the influence of strong gravity. Our work may possibly be generalized to calculate the maximal number of equilibrium points certain autonomous dynamical systems can have in general.

Understanding nonlinear structure formation is crucial for the full exploration of the data generated by forthcoming stage IV surveys, which will require accurate modelling of the matter power spectrum. This is particularly challenging for deviations from $\Lambda$CDM. We need to ensure that alternative cosmologies are well tested and avoid false detections. In this work we present an extension of the halo model reaction framework for interacting dark energy. We describe modifications to the halo model to include the additional force acting on dark matter within the Dark Scattering model, and implement them into the code ReACT. The halo model reaction is then combined with a pseudo spectrum from the EuclidEmulator2 and compared to Dark Scattering N-body simulations. Using a standard halo mass function and concentration-mass relation, we find these predictions to be 1% accurate at $z=0$ up to $k=0.8~h/{\rm Mpc}$ for the largest interaction strength tested ($\xi=50$ b/GeV), improving to $k=2~h/{\rm Mpc}$ at $z=1$. For a smaller interaction strength ($\xi=10$ b/GeV), we find 1%-agreement at $z=1$ up to at least $k=3.5~h/{\rm Mpc}$, being close to $k=1~h/{\rm Mpc}$ at $z=0$. Finally, we improve our predictions with the inclusion of baryonic feedback and massive neutrinos and search for degeneracies between the nonlinear effects of these two contributions and those of the dark sector interaction. By limiting the scales to those where our modelling is 1% accurate, we find a degeneracy between the effects of the interaction and those of baryonic feedback, but not with those generated by massive neutrinos. We expect the degeneracy with baryonic feedback to be resolvable when smaller scales are included. This work represents the first analytical tool for calculating the nonlinear matter power spectrum for coupled dark matter - dark energy models.

The binary neutron star (BNS) mass distribution measured with gravitational-wave observations has the potential to reveal information about the dense matter equation of state, supernova physics, the expansion rate of the universe, and tests of General Relativity. As most current gravitational-wave analyses measuring the BNS mass distribution do not simultaneously fit the spin distribution, the implied population-level spin distribution is the same as the spin prior applied when analyzing individual sources. In this work, we demonstrate that introducing a mismatch between the implied and true BNS spin distributions can lead to biases in the inferred mass distribution. This is due to the strong correlations between the measurements of the mass ratio and spin components aligned with the orbital angular momentum for individual sources. We find that applying a low-spin prior which excludes the true spin magnitudes of some sources in the population leads to significantly overestimating the maximum neutron star mass and underestimating the minimum neutron star mass at the population level with as few as six BNS detections. The safest choice of spin prior that does not lead to biases in the inferred mass distribution is one which allows for high spin magnitudes and tilts misaligned with the orbital angular momentum.

Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion time scale, this has not been directly constrained using observations. Using the unprecedented amount and quality of data of the inner regions of the Milky Way that has recently been produced by APOGEE and Gaia, we aim to derive strong constraints on the infall history of the inner (less than 6 kpc) Galaxy (with a focus on stars between 4-6 kpc, which we show is an appropriate proxy for the entire inner disc). We have implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt-Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time-[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius between 4-6 kpc from APOGEE. Our results point to a fast early gas accretion, with an upper limit of accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2-3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.

We present high-speed optical photometry from SAAO and SALT on the black hole LMXB \MAXI (ASSASN-18ey), some of it simultaneous with NICER, Swift and Insight-HXMT X-ray coverage. We detect optical Quasi-Periodic Oscillations (QPOs) that move to higher frequencies as the outburst progresses, tracking both the frequency and evolution of similar X-ray QPOs previously reported. Correlated X-ray/optical data reveal a complex pattern of lags, including an anti-correlation and a sub-second lag that evolve over the first few weeks of outburst. They also show correlated components separated by a lag equal to the QPO period roughly centered on zero lag, implying that the inter-band variability is strongly and consistently affected by these QPOs at a constant phase lag of roughly +/- pi. The synchronisation of X-ray and optical QPOs indicates that they must be produced in regions physically very close to each other; we thus propose that they can be explained by a precessing jet model, based on analogies with V404 Cyg and MAXI J1348-630.

We explore the possibility that relativistic protons in the extremely powerful jets of blazars may boost via elastic collisions the dark matter particles in the surroundings of the source to high energies. We concentrate on two sample blazars, TXS 0506+056 - towards which IceCube recently reported evidence for a high-energy neutrino flux - and BL Lacertae, a representative nearby blazar. We find that the dark matter flux at Earth induced by these sources may be sizeable, larger than the flux associated with the analogous process of DM boosted by galactic cosmic rays, and relevant to access direct detection for dark matter particle masses lighter than 1 GeV. From the null detection of a signal by XENON1T, MiniBooNE, and Borexino, we derive limits on dark matter-nucleus spin-independent and spin-dependent cross sections which, depending on the modelization of the source, improve on other currently available bounds for light DM candidates of one up to five orders of magnitude.

Dark sector particles at the GeV scale carrying baryon number provide an attractive framework for understanding the origin of dark matter and the matter-antimatter asymmetry of the universe. We demonstrate that dark decays of hadronic states containing strange quarks -- hyperons -- offer excellent prospects for discovering such dark baryons. Building up on novel calculations of the matrix elements relevant for hyperon dark decays, and in view of various collider, flavor, and astrophysical constraints, we determine the expected rates at hyperon factories like BESIII and LHCb. We also highlight the interesting theoretical connections of hyperon dark decays to the neutron lifetime anomaly and Mesogenesis.

We resurrect the infamous harmonic mean estimator for computing the marginal likelihood (Bayesian evidence) and solve its problematic large variance. The marginal likelihood is a key component of Bayesian model selection since it is required to evaluate model posterior probabilities; however, its computation is challenging. The original harmonic mean estimator, first proposed in 1994 by Newton and Raftery, involves computing the harmonic mean of the likelihood given samples from the posterior. It was immediately realised that the original estimator can fail catastrophically since its variance can become very large and may not be finite. A number of variants of the harmonic mean estimator have been proposed to address this issue although none have proven fully satisfactory. We present the learnt harmonic mean estimator, a variant of the original estimator that solves its large variance problem. This is achieved by interpreting the harmonic mean estimator as importance sampling and introducing a new target distribution. The new target distribution is learned to approximate the optimal but inaccessible target, while minimising the variance of the resulting estimator. Since the estimator requires samples of the posterior only it is agnostic to the strategy used to generate posterior samples. We validate the estimator on a variety of numerical experiments, including a number of pathological examples where the original harmonic mean estimator fails catastrophically. In all cases our learnt harmonic mean estimator is shown to be highly accurate. The estimator is computationally scalable and can be applied to problems of dimension $\mathcal{O}(10^3)$ and beyond. Code implementing the learnt harmonic mean estimator is made publicly available.

Gravitational waves provide a novel way to probe axions or axion-like particles coupled to a dark photon field, even in the absence of couplings to Standard Model particles. In the conventional misalignment mechanism, the generation of an observable stochastic gravitational wave background, however, requires large axion decay constants. We here investigate the gravitational wave signal generated within the kinetic misalignment scenario, where the axion is assumed to have a large initial velocity. Its kinetic energy then provides a sufficiently high energy budget to generate a detectable gravitational wave signal also at lower values of the decay constant. We obtain an analytic estimate as well as perform numerical simulations of the corresponding gravitational wave signal, and evaluate its detectability at current and future gravitational wave observatories. We further present the corresponding projected constraints on the parameter space of the model, along with the parameter regions in which the dark photon or axion constitute dark matter, or in which the baryon asymmetry of the Universe is generated via the axiogenesis mechanism. Finally we compute the GW production from the fragmentation of rotating axions, which is however difficult to observe experimentally.

We report precision mass measurements of neutron-deficient gallium isotopes approaching the proton drip line. The measurements of $^{60-63}$Ga performed with TITAN's multiple-reflection time-of-flight mass spectrometer provide a more than threefold improvement over the current literature mass uncertainty of $^{61}$Ga and mark the first direct mass measurement of $^{60}$Ga. The improved precision of the $^{61}$Ga mass has important implications for the astrophysical rp-process, as it constrains essential reaction Q-values near the $^{60}$Zn waiting point. Based on calculations with a one-zone model, we demonstrate the impact of the improved mass data on prediction uncertainties of X-ray burst models. The first-time measurement of the $^{60}$Ga ground-state mass establishes the proton-bound nature of this nuclide; thus, constraining the location of the proton drip line along this isotopic chain. Including the measured mass of $^{60}$Ga enables us to extend the evaluated $T=1$ isobaric multiplet mass equation up to $A=60$.

We present an inequality between two types of distance measures to a single source in general relativity. It states that for a given emitter and observer the distance between them measured by the trigonometric parallax is never shorter than the angular diameter distance provided that the null energy condition holds and that there are no focal points in between. This result is independent of the details of the spacetime geometry or the motions of the observer and the source. The proof is based on the geodesic bilocal operator formalism together with well known properties of infinitesimal light ray bundles. Observation of the violation of the distance inequality would mean that on large scales either the null energy condition does not hold or that light does not travel along null geodesics.

The various materials of test masses, and the difference of arm lengths of global ground-based gravitational-wave interferometer detectors offer a unique approach to test Newton's second law, weak equivalence principle, and Einstein equivalence principle with dynamical space-time effects in terms of the interaction of gravitational waves with detectors. We proposed a novel test strategy for the interaction between gravitational waves and detectors, which is independent of particular gravitation theory. A new population level of the Fisher-Matrix approach for multiple sources and multiple detectors case is formalized to evaluate the prospects for a binary neutron star and binary black hole coalescences. Through a generalized detector response, we found more sources could break the parameter degeneracy and one could constrain the interaction and gravitational-inertial mass ratio parameters with the standard deviation $\ls 1\%$ with about 10 compact binary coalescence sources with future third-generation detectors network.

Primordial black holes can be produced by density fluctuations generated from delayed vacuum decays of first-order phase transition. The primordial black holes generated at the electroweak phase transition have masses of about $10^{-5}$ solar mass. Such primordial black holes in the mass range can be tested by current and future microlensing observations, such as Subaru HSC, OGLE, PRIME and Roman telescope. Therefore, we may be able to explore new physics models with strongly first-order electroweak phase transition via primordial black holes. We examine this possibility by using models with first-order electroweak phase transition in the standard model effective field theory with dimension 6 and 8 operators. We find that depending on parameters of the phase transition a sufficient number of primordial black holes can be produced to be observed by above mentioned experiments. Our results would suggest that primordial black holes can be used as a new probe of models with strongly first-order electroweak phase transition, which has complementarity with measurements of the triple Higgs boson coupling at future collider experiments and observations of gravitational waves at future space-based interferometers.

Simulation-based inference with conditional neural density estimators is a powerful approach to solving inverse problems in science. However, these methods typically treat the underlying forward model as a black box, with no way to exploit geometric properties such as equivariances. Equivariances are common in scientific models, however integrating them directly into expressive inference networks (such as normalizing flows) is not straightforward. We here describe an alternative method to incorporate equivariances under joint transformations of parameters and data. Our method -- called group equivariant neural posterior estimation (GNPE) -- is based on self-consistently standardizing the "pose" of the data while estimating the posterior over parameters. It is architecture-independent, and applies both to exact and approximate equivariances. As a real-world application, we use GNPE for amortized inference of astrophysical binary black hole systems from gravitational-wave observations. We show that GNPE achieves state-of-the-art accuracy while reducing inference times by three orders of magnitude.

An effect of the Lorentz symmetry breaking is pointed out in the cosmological context. Using a Bianchi I geometry coupled to the Kalb-Ramond field, a consequence of the Lorentz symmetry violation is indicated by a different rate of expansion in a given spatial direction. This article focuses on the coupling constant $\xi_1$, which generates, from the Kalb-Ramond field, all three coefficients that give rise to the Lorentz violation in the gravity sector of the minimal Standard Model Extension. The coupling constant $\xi_1$ increases the rate of expansion of the universe in a given direction during a dark energy era. As a consequence, a range of validity of that coupling constant is also obtained.

The high-frequency gravitons can be absorbed by the first and second viscosities of the post-inflationary plasma as the corresponding wavelengths reenter the Hubble radius prior to big-bang nucleosynthesis. When the total sound speed of the medium is stiffer than radiation the rate of expansion still exceeds the shear rate but the bulk viscosity is not negligible. Depending on the value of the entropy density at the end of inflation the spectral energy density of the relic gravitons gets modified in comparison with the inviscid result when the frequency ranges between the kHz band and the GHz region. In the nHz domain the spectrum inherits a known suppression due to neutrino free-streaming but also a marginal spike potentially caused by a sudden outbreak of the bulk viscosity around the quark-hadron phase transition, as suggested by the hadron spectra produced in the collisions of heavy ions.

Throughout recorded history, humans have crossed national borders to seek safety in nearby countries. The reasons for displacement have been generated by phenomena of terrestrial origin, but exposure to unexpected extra-terrestrial threats poses a different scenario. An asteroid impact warning implies a change of paradigm which would represent a historic precedent. In this regard, the analogies with natural disasters must be considered, along with multiple possible scenarios, and legal aspects related to a) the legal framework to regulate this situation; b) the action and responsibility of the states; and c) the definition of impact refugee and the reconfiguration of traditional concepts such as deterritorialized states. In addition, the decision-making process and the actors involved must be led by a cooperative effort to improve international law. These new circumstances should be established with a consideration of inequalities between the states, and an aim of protecting humanity through democratic solutions using the safest, most effective techniques.

We perform several consistency tests between different phases of binary black hole dynamics; the inspiral, the merger, and the ringdown on the gravitational wave events GW150914 and GW170814. These tests are performed explicitly in the time domain, without any spectral leakage between the different phases. We compute posterior distributions on the mass and spin of the initial black holes and the final black hole. We also compute the initial areas of the two individual black holes and the final area from the parameters describing the remnant black hole. This facilitates a test of Hawking's black hole area theorem. We use different waveform models to quantify systematic waveform uncertainties for the area increase law with the two events. We find that these errors may lead to overstating the confidence with which the area theorem is confirmed. For example, we find $>99\%$ agreement with the area theorem for GW150914 if a damped sinusoid consisting of a single-mode is used at merger to estimate the final area. This is because this model overestimates the final mass. Including an overtone of the dominant mode decreases the confidence to $\sim94\%$; using a full merger-ringdown model further decreases the confidence to $\sim 85-90\%$. We find that comparing the measured change in the area to the expected change in area yields a more robust test, as it also captures over estimates in the change of area. We find good agreement with GR when applying this test to GW150914 and GW170814.