Papers for Friday, Jul 09 2021

Previous studies of the stellar mean metallicity and [Mg/Fe] values of massive elliptical (E)~galaxies suggest that their stars were formed in a very short timescale which cannot be reconciled with estimates from stellar population synthesis (SPS) studies and with hierarchical-assembly. Applying the previously developed chemical evolution code, GalIMF, which allows an environment-dependent stellar initial mass function (IMF) to be applied in the integrated galaxy initial mass function (IGIMF) theory instead of an invariant canonical IMF, the star formation timescales (SFT) of E galaxies are re-evaluated. The code's uniqueness lies in it allowing the galaxy-wide IMF and associated chemical enrichment to evolve as the physical conditions in the galaxy change. The calculated SFTs become consistent with the independent SPS results if the number of type Ia supernovae (SNIa) per unit stellar mass increases for more massive E~galaxies. This is a natural outcome of galaxies with higher star-formation rates producing more massive star clusters, spawning a larger number of SNIa progenitors per star. The calculations show E~galaxies with a stellar mass $\approx 10^{9.5} M_\odot$ to have had the longest mean SFTs of $\approx2\,$Gyr. The bulk of more massive E~galaxies were formed faster (SFT$\,\approx 1\,$Gyr) leading to domination by M~dwarf stars and larger dynamical mass-to-light ratios as observed, while lower-mass galaxies tend to lose their gas supply more easily due to their shallower potential and therefore also have similarly-short mean SFTs. This work achieves, for the first time, consistency of the SFTs for early-type galaxies between chemical-enrichment and SPS modelling and leads to an improved understanding of how the star formation environment may affect the total number of SNIa per unit stellar mass formed.

We discuss the usefulness and theoretical consistency of different entropy variables used in the literature to describe isocurvature perturbations in multifield inflationary models with a generic curved field space. We clarify which is the proper entropy variable to be used to match the evolution of isocurvature modes during inflation to the one after the reheating epoch in order to compare with observational constraints. In particular, we find that commonly used variables, as the relative entropy perturbation or the one associated to the decomposition in tangent and normal perturbations with respect to the inflationary trajectory, even if more useful to perform numerical studies, can lead to results which are wrong by several orders of magnitude, or even to apparent destabilisation effects which are unphysical for cases with light kinetically coupled spectator fields.

Using a novel suite of cosmological simulations zooming in on a Mpc-scale intergalactic sheet or "pancake" at z~3-5, we conduct an in-depth study of the thermal properties and HI content of the warm-hot intergalactic medium (WHIM) at those redshifts. The simulations span nearly three orders of magnitude in gas-cell mass, from ~(7.7x10^6-1.5x10^4)Msun, one of the highest resolution simulations of such a large patch of the inter-galactic medium (IGM) to date. At z~5, a strong accretion shock develops around the main pancake following a collision between two smaller sheets. Gas in the post-shock region proceeds to cool rapidly, triggering thermal instabilities and the formation of a multiphase medium. We find neither the mass, nor the morphology, nor the distribution of HI in the WHIM to be converged at our highest resolution. Interestingly, the lack of convergence is more severe for the less dense, more metal-poor, intra-pancake medium (IPM) in between filaments and far from any star-forming galaxies. As the resolution increases, the IPM develops a shattered structure, with ~kpc scale clouds containing most of the HI. From our lowest to highest resolution, the covering fraction of metal-poor (Z<10^{-3}Zsun) Lyman-limit systems (NHI>10^{17.2}/cm^2) in the IPM at z~4 increases from (3-15)%, while that of Damped Lyman-alpha Absorbers (NHI>10^{20}/cm^2) with similar metallicity increases threefold, from (0.2-0.6)%, with no sign of convergence. We find that a necessary condition for the formation of a multiphase, shattered structure is resolving the cooling length, lcool=cs*tcool, at T~10^5K. If this scale is unresolved, gas "piles up" at these temperatures and cooling to lower temperatures becomes very inefficient. We conclude that state-of-the-art cosmological simulations are still unable to resolve the multi-phase structure of the low-density IGM, with potentially far-reaching implications.

We investigate the evolutionary nexus between the morphology and internal kinematics of the central regions of collisional, rotating, multi-mass stellar systems, with special attention to the spatial characterisation of the process of mass segregation. We report results from idealized, purely $N$-body simulations that show multi-mass, rotating, and spherical systems rapidly form an oblate, spheroidal massive core, unlike single-mass rotating or multi-mass non-rotating configurations with otherwise identical initial properties, indicating that this evolution is a result of the interplay between the presence of a mass spectrum and angular momentum. This feature appears to be long-lasting, preserving itself for several relaxation times. The degree of flattening experienced by the systems is directly proportional to the initial degree of internal rotation. In addition, this morphological effect has a clear characterisation in terms of orbital architecture, as it lowers the inclination of the orbits of massive stars. We offer an idealised dynamical interpretation that could explain the mechanism underpinning this effect and we highlight possible useful implications, from kinematic hysteresis to spatial distribution of dark remnants in dense stellar systems.

We comprehensively study the variability of Miras in the Large Magellanic Cloud (LMC) by simultaneous analysing light curves in 14 bands in the range of 0.5$-$24 microns. We model over 20-years-long, high cadence $I$-band light curves collected by The Optical Gravitational Lensing Experiment (OGLE) and fit them to light curves collected in the remaining optical/near-infrared/mid-infrared bands to derive both the variability amplitude ratio and phase-lag as a function of wavelength. We show that the variability amplitude ratio declines with the increasing wavelength for both oxygen-rich (O-rich) and carbon-rich (C-rich) Miras, while the variability phase-lag increases slightly with the increasing wavelength. In a significant number of Miras, mostly the C-rich ones, the spectral energy distributions (SEDs) require a presence of a cool component (dust) in order to match the mid-IR data. Based on SED fits for a golden sample of 140 Miras, we calculated synthetic period-luminosity relations (PLRs) in 42 bands for the existing and future sky surveys that include OGLE, The VISTA Near-Infrared $YJK_\mathrm{s}$ Survey of the Magellanic Clouds System (VMC), Legacy Survey of Space and Time (LSST), Gaia, Spitzer, The Wide-field Infrared Survey Explorer (WISE), The James Webb Space Telescope (JWST), The Nancy Grace Roman Space Telescope (formerly WFIRST), and The Hubble Space Telescope (HST). We show that the synthetic PLR slope decreases with increasing wavelength for both the O-rich and C-rich Miras in the range of 0.1$-$40 microns. Finally, we show the location and motions of Miras on the color-magnitude (CMD) and color-color (CCD) diagrams.

Characterizing a planet detected by microlensing is hard if the planetary signal is weak or the lens-source relative trajectory is far from caustics. However, statistical analyses of planet demography must include those planets to accurately determine occurrence rates. As part of a systematic modeling effort in the context of a $>10$-year retrospective analysis of MOA's survey observations to build an extended MOA statistical sample, we analyze the light curve of the planetary microlensing event MOA-2014-BLG-472. This event provides weak constraints on the physical parameters of the lens, as a result of a planetary anomaly occurring at low magnification in the light curve. We use a Bayesian analysis to estimate the properties of the planet, based on a refined Galactic model and the assumption that all Milky Way's stars have an equal planet-hosting probability. We find that a lens consisting of a $1.9^{+2.2}_{-1.2}\,\mathrm{M}_\mathrm{J}$ giant planet orbiting a $0.31^{+0.36}_{-0.19}\,\mathrm{M}_\odot$ host at a projected separation of $0.75\pm0.24\,\mathrm{au}$ is consistent with the observations and is most likely, based on the Galactic priors. The lens most probably lies in the Galactic bulge, at $7.2^{+0.6}_{-1.7}\mathrm{kpc}$ from Earth. The accurate measurement of the measured planet-to-host star mass ratio will be included in the next statistical analysis of cold planet demography detected by microlensing.

We present results on the measured shapes of 832 galaxies in 11 galaxy clusters at 1.0 < z <1.4 from the GOGREEN survey. We measure the axis ratio ($q$), the ratio of the minor to the major axis, of the cluster galaxies from near-infrared Hubble Space Telescope imaging using S\'ersic profile fitting and compare them with a field sample. We find that the median $q$ of both star-forming and quiescent galaxies in clusters increases with stellar mass, similar to the field. Comparing the axis ratio distributions between clusters and the field in four mass bins, the distributions for star-forming galaxies in clusters are consistent with those in the field. Conversely, the distributions for quiescent galaxies in the two environments are distinct, most remarkably in $10.1\leq\log(M/{\rm M}_{\odot})<10.5$ where clusters show a flatter distribution, with an excess at low $q$. Modelling the distribution with oblate and triaxial components, we find that the cluster and field sample difference is consistent with an excess of flattened oblate quiescent galaxies in clusters. The oblate population contribution drops at high masses, resulting in a narrower $q$ distribution in the massive population than at lower masses. Using a simple accretion model, we show that the observed $q$ distributions and quenched fractions are consistent with a scenario where no morphological transformation occurs for the environmentally quenched population in the two intermediate mass bins. Our results suggest that environmental quenching mechanism(s) likely produce a population that has a different morphological mix than those resulting from the dominant quenching mechanism in the field.

Normal-mode oscillation frequencies computed from stellar models differ from those which would be measured from stars with identical interior structures, because of modelling errors in the near-surface layers. These frequency differences are referred to as the asteroseismic "surface term". The vast majority of solar-like oscillators which have been observed, and which are expected to be observed in the near future, are evolved stars which exhibit mixed modes. For these evolved stars, the inference of stellar properties from these mode frequencies has been shown to depend on how this surface term is corrected for. We show that existing parametrisations of the surface term account for mode mixing only to first order in perturbation theory, if at all, and therefore may not be adequate for evolved stars. Moreover, existing nonparametric treatments of the surface term do not account for mode mixing. We derive both a first-order construction, and a more general approach, for one particular class of nonparametric methods. We illustrate the limits of first-order approximations from both analytic considerations and using numerical injection-recovery tests on stellar models. First-order corrections for the surface term are strictly only applicable where the size of the surface term is much smaller than both the coupling strength between the mixed p- and g-modes, as well as the local g-mode spacing. Our more general matrix construction may be applied to evolved stars, where perturbation theory cannot be relied upon.

Understanding the impact of environment on the formation and evolution of dark matter halos and galaxies is a crucial open problem. Studying statistical correlations in large simulated populations sheds some light on these impacts, but the causal effect of an environment on individual objects is harder to pinpoint. Addressing this, we present a new method for resimulating a single dark matter halo in multiple large-scale environments. In the initial conditions, we 'splice' (i.e. insert) the Lagrangian region of a halo into different Gaussian random fields, while enforcing consistency with the statistical properties of $\Lambda$CDM. Applying this technique, we demonstrate that the mass of halos is primarily determined by the density structure inside their Lagrangian patches, while the halos' concentration is more strongly affected by environment. The splicing approach will also allow us to study, for example, the impact of the cosmic web on accretion processes and galaxy quenching.

Low-mass ($M_{\rm{500}}<5\times10^{14}{\rm{M_\odot}}$) galaxy clusters have been largely unexplored in radio observations, due to the inadequate sensitivity of existing telescopes. However, the upgraded GMRT (uGMRT) and the Low Frequency ARray (LoFAR), with unprecedented sensitivity at low frequencies, have paved the way to closely study less massive clusters than before. We have started the first large-scale programme to systematically search for diffuse radio emission from low-mass galaxy clusters, chosen from the Planck Sunyaev-Zel'dovich cluster catalogue. We report here the detection of diffuse radio emission from four of the 12 objects in our sample, shortlisted from the inspection of the LoFAR Two Meter Sky Survey (LoTSS-I), followed up by uGMRT Band 3 deep observations. The clusters PSZ2~G089 (Abell~1904) and PSZ2~G111 (Abell~1697) are detected with relic-like emission, while PSZ2~G106 is found to have an intermediate radio halo and PSZ2~G080 (Abell~2018) seems to be a halo-relic system. PSZ2~G089 and PSZ2~G080 are among the lowest-mass clusters discovered with a radio relic and a halo-relic system, respectively. A high ($\sim30\%$) detection rate, with powerful radio emission ($P_{1.4\ {\rm GHz}}\sim10^{23}~{\rm{W~Hz^{-1}}}$) found in most of these objects, opens up prospects of studying radio emission in galaxy clusters over a wider mass range, to much lower-mass systems.

We present a test platform for the Athena X-IFU detection chain, which will serve as the first demonstration of the representative end-to-end detection and readout chain for the X-IFU, using prototypes of the future flight electronics and currently available subsystems. This test bench, housed in a commercial two-stage ADR cryostat, includes a focal plane array placed at the 50 mK cold stage of the ADR with a kilopixel array of transition-edge sensor microcalorimeter spectrometers and associated cold readout electronics. Prototype room temperature electronics for the X-IFU provide the readout, and will evolve over time to become more representative of the X-IFU mission baseline. The test bench yields critical feedback on subsystem designs and interfaces, in particular the warm readout electronics, and will provide an in-house detection system for continued testing and development of the warm readout electronics and for the validation of X-ray calibration sources. In this paper, we describe the test bench subsystems and design, characterization of the cryostat, and current status of the project.

Ultra-short-period (USP) planets reside inside the expected truncation radius for typical T Tauri disks. As a result, their current orbital locations require an explanation beyond standard disk migration or in situ formation. Modern theories of planet-disk interactions indicate that once a planet migrates close to the disk's inner truncation radius, Type I torques vanish or switch direction, depending on the stellar and disk conditions, so that the planet is expected to stop its orbital decay and become trapped. In this work, we show that that magnetically-driven sub-Keplerian gas flow in the inner disk can naturally counteract these effects and produce systems with USP planets at their observed orbital radii. The sub-Keplerian gas flow provides a headwind to small planets, and the resulting torque can overcome the effects of outward Type I migration near the co-rotation radius. For suitable disk and planet parameters, the torques due to the sub-Keplerian gas flow lead to inward migration on a rapid timescale. Over the time span of an FU Ori outburst, which moves the disk truncation radius inward, the rapid headwind migration can place planets in USP orbits. The combination of headwind migration and FU Ori outbursts thus provides a plausible mechanism to move small planets from $a=0.05-0.1$ AU down to $a=0.01-0.02$ AU. This effect is amplified for low-mass planets, consistent with existing observations.

We present a detailed study of the shapes and alignments of different galaxy cluster components using hydrodynamical simulations. We compute shape parameters from the Dark Matter (DM) distribution, the galaxy members and the intra-cluster light (ICL). We assess how well the DM cluster shape can be constrained by means of the identified galaxy member positions and the ICL. Further, we address the dilution factor introduced when estimating the cluster elongation using weak-lensing stacking techniques, which arises due to the misalignment between the total surface mass distribution and the distribution of luminous tracers. The dilution is computed considering the alignment between the DM and the Brightest Cluster Galaxy, the galaxy members and the ICL. Our study shows that distributions of galaxy members and ICL are less spherical than the DM component, although both are well aligned with the semi-major axis of the later. We find that the distribution of galaxy members hosted in more concentrated subhalos is more elongated than the distribution of the DM. Moreover, these galaxies are better aligned with the dark matter component compared to the distribution of galaxies hosted in less concentrated subhalos. We conclude that the positions of galaxy members can be used as suitable tracers to estimate the cluster surface density orientation, even when a low number of members is considered. Our results provide useful information for interpreting the constraints on the shapes of galaxy clusters in observational studies.

We use photometric redshifts and statistical background subtraction to measure stellar mass functions in galaxy group-mass ($4.5-8\times10^{13}~\mathrm{M}_\odot$) haloes at $1<z<1.5$. Groups are selected from COSMOS and SXDF, based on X-ray imaging and sparse spectroscopy. Stellar mass ($M_{\mathrm{stellar}}$) functions are computed for quiescent and star-forming galaxies separately, based on their rest-frame $UVJ$ colours. From these we compute the quiescent fraction and quiescent fraction excess (QFE) relative to the field as a function of $M_{\mathrm{stellar}}$. QFE increases with $M_{\mathrm{stellar}}$, similar to more massive clusters at $1<z<1.5$. This contrasts with the apparent separability of $M_{\mathrm{stellar}}$ and environmental factors on galaxy quiescent fractions at $z\sim 0$. We then compare our results with higher mass clusters at $1<z<1.5$ and lower redshifts. We find a strong QFE dependence on halo mass at fixed $M_{\mathrm{stellar}}$; well fit by a logarithmic slope of $\mathrm{d}(\mathrm{QFE})/\mathrm{d}\log (M_{\mathrm{halo}}) \sim 0.24 \pm 0.04$ for all $M_{\mathrm{stellar}}$ and redshift bins. This dependence is in remarkably good qualitative agreement with the hydrodynamic simulation BAHAMAS, but contradicts the observed dependence of QFE on $M_{\mathrm{stellar}}$. We interpret the results using two toy models: one where a time delay until rapid (instantaneous) quenching begins upon accretion to the main progenitor ("no pre-processing") and one where it starts upon first becoming a satellite ("pre-processing"). Delay times appear to be halo mass dependent, with a significantly stronger dependence required without pre-processing. We conclude that our results support models in which environmental quenching begins in low-mass ($<10^{14}M_\odot$) haloes at $z>1$.

Radio astronomy and Space Infrastructures in the Azores have a great scientific and industrial interest because they benefit from a unique geographical location in the middle of the North Atlantic allowing a vast improvement in the sky coverage. This fact obviously has a very high added value for: i) the establishment of space tracking and communications networks for the emergent global small satellite fleets ii) it is invaluable to connect the radio astronomy infrastructure networks in Africa, Europe and America continents using Very Large Baseline Interferometry (VLBI) techniques, iii) it allows excellent potential for monitoring space debris and Near Earth Objects (NEOs). There is in S. Miguel island a 32-metre SATCOM antenna that could be integrated in advanced VLBI networks and be capable of additional Deep Space Network ground support. This paper explores the space science opportunities offered by the upgrade of the S. Miguel 32-metre SATCOM antenna into a world-class infrastructure for radio astronomy and space exploration: it would enable a Deep Space Network mode and would constitute a key space facility for data production, promoting local digital infrastructure investments and the testing of cutting-edge information technologies. Its Atlantic location also enables improvements in angular resolution, provides many baseline in East-West and North-South directions connecting the emergent VLBI stations in America to Europe and Africa VLBI arrays therefore contributing for greater array imaging capabilities especially for sources or well studied fields close to or below the celestial equator, where ESO facilities, ALMA, SKA and its precursors do or will operate and observe in the coming decades.

A dynamical model for large near-Earth asteroids (NEAs) is developed here to understand the occurrence rate and nature of Cretaceous-Paleogene (K/Pg) scale impacts on the Earth. We find that 16--32 (2--4) impacts of diameter $D>5$ km ($D>10$ km) NEAs are expected on the Earth in 1 Gyr, with about a half of impactors being dark primitive asteroids (most of which start with semimajor axis $a>2.5$ au). These results explain why the Chicxulub crater, the third largest impact structure found on the Earth (diameter $\simeq180$ km), was produced by impact of a carbonaceous chondrite. They suggest, when combined with previously published results for small ($D \lesssim 1$ km) NEAs, a size-dependent sampling of the main belt. We conclude that the impactor that triggered the K/Pg mass extinction $\simeq 66$ Myr ago was a main belt asteroid that quite likely ($\simeq 60$\% probability) originated beyond 2.5 au.

We analysed the open clusters Czernik 2 and NGC 7654 using CCD UBV photometric and Gaia Early Data Release 3 (EDR3) photometric and astrometric data. Structural parameters of the two clusters were derived, including the physical sizes of Czernik 2 being r=5 and NGC 7654 as 8 min. We calculated membership probabilities of stars based on their proper motion components as released in the Gaia EDR3. To identify member stars of the clusters, we used these membership probabilities taking into account location and the impact of binarity on main-sequence stars. We used membership probabilities higher than $P=0.5$ to identify 28 member stars for Czernik 2 and 369 for NGC 7654. We estimated colour-excesses and metallicities separately using two-colour diagrams to derive homogeneously determined parameters. The derived $E(B-V)$ colour excess is 0.46(0.02) mag for Czernik 2 and 0.57(0.04) mag for NGC 7654. Metallicities were obtained for the first time for both clusters, -0.08(0.02) dex for Czernik 2 and -0.05(0.01) dex for NGC 7654. Keeping the reddening and metallicity as constant quantities, we fitted PARSEC models using colour-magnitude diagrams, resulting in estimated distance moduli and ages of the two clusters. We obtained the distance modulus for Czernik 2 as 12.80(0.07) mag and for NGC 7654 as 13.20(0.16) mag, which coincide with ages of 1.2(0.2) Gyr and 120(20) Myr, respectively. The distances to the clusters were calculated using the Gaia EDR3 trigonometric parallaxes and compared with the literature. We found good agreement between the distances obtained in this study and the literature. Present day mass function slopes for both clusters are comparable with the value of Salpeter (1955), being X=-1.37(0.24) for Czernik 2 and X=-1.39(0.19) for NGC 7654.

In an effort to improve our understanding of the spiral arm structure of the Milky Way, we use Classical Cepheids (CCs) to increase the number of young tracers on the far side of the Galactic disk with accurately determined distances. We use a sample of 30 CCs, discovered using near-infrared photometry from the VISTA Variables in the V\'ia L\'actea survey (VVV) and classified based on their radial velocities and metallicities. We combine them with another 20 CCs from the literature for which VVV photometry is available. The compiled sample of CCs with homogeneously computed distances based on VVV infrared photometry was employed as a proof of concept to trace the spiral structure in the poorly explored far side of the disk. Although the use of CCs has some caveats, these variables are currently the only available young tracers in the far side disk for which a numerous sample with accurate distances can be obtained. Therefore, a larger sample could allow us to make a significant step forward in our understanding of the Milky Way disk as a whole. We present preliminary evidence that CCs favor: a spiral arm model with two main arms (Perseus and Scutum-Centaurus) branching out into four arms at galactocentric distances, $R_\mathrm {GC}\gtrsim5-6\,\mathrm{kpc}$; the extension of the Scutum-Centaurus arm behind the Galactic center; a possible connection between the Perseus arm and the Norma tangency direction. The current sample of CCs in the far side of the Galaxy are in the mid-plane, arguing against the presence of a severely warped disk at small Galactocentric distances ($R_\mathrm {GC}\lesssim12\,\mathrm{kpc}$) in the studied area. The discovery and characterization of CCs at near-IR wavelengths appears to be a promising tool to complement studies based on other spiral arm tracers and extend them to the far side of our Galaxy.

Cold dark matter (CDM) has faced a number of challenges mainly at small scales, such as the too-big-to-fail problem, and core-cusp density profile of dwarf galaxies. Such problems were argued to have a solution either in the baryonic physics sector or in modifying the nature of dark matter to be self-interacting, or self-annihilating, or ultra-light. Here we present a new challenge for CDM by showing that some of Milky Way's satellites are too dense, requiring the formation masses and redshifts of halos in CDM not compatible with being a satellite. These too-dense-to-be-satellite systems are dominated by dark matter and exhibit a surface density above mean dark matter cosmic surface density $\sim\Omega_{dm} \rho_c c/H_0\approx 200~\rm M_{\odot}/pc^2$. This value corresponds to dark matter pressure of $\approx 10^{-10}{\rm erg/cm^3}$. This problem, unlike other issues facing CDM, has no solution in the baryonic sector and none of the current alternatives of dark matter can account for it. The too-dense-to-be-satellite problem presented in this work provides a new clue for the nature of dark matter, never accounted for before.

Various nucleosynthesis studies have pointed out that the r-process elements in very metal-poor (VMP) halo stars might have different origins. By means of familiar concepts from statistics (correlations, cluster analysis, rank tests of elemental abundances), we look for causally correlated elemental abundance patterns and attempt to link them to astrophysical events. Some of these events produce the r-process elements jointly with iron, while others do not have any significant iron contribution. In the early stage of our Galaxy, at least three r-process nucleosynthesis sites have been active. The first two produce and eject iron and the majority of the lighter r-process elements. We assign them to two different types of core-collapse events, not identical to regular core-collapse supernovae (CCSNe), which produce only light trans-Fe elements. The third category is characterized by a strong r-process and responsible for the major fraction of the heavy main r-process elements without a significant co-production of Fe. It does not appear to be connected to CCSNe, in fact the Fe found in the related r-process enriched stars must come from previously occurring CCSNe. The existence of actinide boost stars indicates a further division among strong r-process sites. We assign these two strong r-process sites to neutron star mergers without fast black hole formation and to events where the ejecta are dominated by black hole accretion disk outflows. Indications from the lowest-metallicity stars hint at a connection to massive single stars (collapsars) forming black holes in the early Galaxy.

Polarimetric observations of Fast Radio Bursts (FRBs) are a powerful resource for better understanding these mysterious sources by directly probing the emission mechanism of the source and the magneto-ionic properties of its environment. We present a pipeline for analysing the polarized signal of FRBs captured by the triggered baseband recording system operating on the FRB survey of The Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB). Using a combination of simulated and real FRB events, we summarize the main features of the pipeline and highlight the dominant systematics affecting the polarized signal. We compare parametric (QU-fitting) and non-parametric (rotation measure synthesis) methods for determining the Faraday rotation measure (RM) and find the latter method susceptible to systematic errors from known instrumental effects of CHIME/FRB observations. These errors include a leakage artefact that appears as polarized signal near $\rm{RM\sim 0 \; rad \, m^{-2}}$ and an RM sign ambiguity introduced by path length differences in the system's electronics. We apply the pipeline to a bright burst previously reported by \citet[FRB 20191219F;][]{Leung2021}, detecting an $\mathrm{RM}$ of $\rm{+6.074 \pm 0.006 \pm 0.050 \; rad \, m^{-2}}$ with a significant linear polarized fraction ($\gtrsim0.87$) and strong evidence for a non-negligible circularly polarized component. Finally, we introduce an RM search method that employs a phase-coherent de-rotation algorithm to correct for intra-channel depolarization in data that retain electric field phase information, and successfully apply it to an unpublished FRB, FRB 20200917A, measuring an $\mathrm{RM}$ of $\rm{-1294.47 \pm 0.10 \pm 0.05 \; rad \, m^{-2}}$ (the second largest unambiguous RM detection from any FRB source observed to date).

Recent discoveries of anomalously bright radar reflections below the Mars South Polar Layered Deposit (SPLD) have sparked new speculation that liquid water may be present below the ice cap. The reflections, discovered in data acquired by the Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on board the Mars Express orbiter, were interpreted as reflections from damp materials or even subsurface ponds and lakes similar to those found beneath Earth's ice sheets. Recent studies, however, have questioned the feasibility of melting and maintaining liquid water below the SPLD. Herein, we compare radar simulations to MARSIS observations in order to present an alternate hypothesis: that the bright reflections are the result of interference between multiple layer boundaries, with no liquid water present. This new interpretation is more consistent with known conditions on modern Mars.

Data analysis in modern science using extensive experimental and observational facilities, such as a gravitational wave detector, is essential in the search for novel scientific discoveries. Accordingly, various techniques and mathematical principles have been designed and developed to date. A recently proposed approximate correlation method based on the information theory is widely adopted in science and engineering. Although the maximal information coefficient (MIC) method remains in the phase of improving its algorithm, it is particularly beneficial in identifying the correlations of multiple noise sources in gravitational-wave detectors including non-linear effects. This study investigates various prospects for determining MIC parameters to improve the reliability of handling multi-channel time-series data, reduce high computing costs, and propose a novel method of determining optimized parameter sets for identifying noise correlations in gravitational wave data.

Understanding the characteristics of the solar magnetic field is essential for interpreting solar activities and dynamo. In this research, we investigated the asymmetric distribution of the solar photospheric magnetic field values, using synoptic charts constructed from space-borne high-resolution magnetograms. It is demonstrated that the Lorentzian function describes the distribution of magnetic field values in the synoptic charts much better than the Gaussian function, and this should reflect the gradual decay process from strong to weak magnetic fields. The asymmetry values are calculated under several circumstances, and the results generally show two periodicities related to the variation of the solar B$_0$ angle and the solar cycle, respectively. We argue that it is the small-scale magnetic fields, the inclination of the solar axis, the emergence and evolution of magnetic flux, and the polar fields that are responsible for the features of asymmetry values. We further determined the polar field reversal time of solar cycles 23 and 24 with the flip of asymmetry values. Specifically, for cycle 24, we assert that the polar polarities of both hemispheres reversed at the same time - in March 2014; as to cycle 23, the reversal time of the S-hemisphere is March 2001, while the determination of the N-hemisphere is hampered by missing data.

We use the hydrodynamical EAGLE simulation to predict the numbers, masses and radial distributions of tidally stripped galaxy nuclei in massive galaxy clusters, and compare these results to observations of ultra-compact dwarf galaxies (UCDs) in the Virgo cluster. We trace the merger trees of galaxies in massive galaxy clusters back in time and determine the numbers and masses of stripped nuclei from galaxies disrupted in mergers. The spatial distribution of stripped nuclei in the simulations is consistent with those of UCDs surrounding massive galaxies in the Virgo cluster. Additionally, the numbers of stripped nuclei are consistent with the numbers of M > $10^{7}~M_{\odot}$ UCDs around individual galaxies and in the Virgo cluster as a whole. The mass distributions in this mass range are also consistent. We find that the numbers of stripped nuclei surrounding individual galaxies correlates better with the stellar or halo mass of individual galaxies than the total cluster mass. We conclude that most high mass (M > $10^{7}~M_{\odot}$ UCDs are likely stripped nuclei. It is difficult to draw reliable conclusions about low mass (M < $10^{7}~M_{\odot}$ UCDs because of observational selection effects. We additionally predict that a few hundred stripped nuclei below a mass of $2~\times~10^{6}~M_{\odot}$ should exist in massive galaxies that will overlap in mass with the globular cluster population. Approximately 1-3 stripped nuclei in the process of forming also exist per massive galaxy.

Signatures of vertical disequilibrium have been observed across the Milky Way's disk. These signatures manifest locally as unmixed phase-spirals in $z$--$v_z$ space ("snails-in-phase") and globally as nonzero mean $z$ and $v_z$ which wraps around as a physical spiral across the $x$--$y$ plane ("snails-in-space"). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle Milky Way (MW)-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using a $z$--$v_z$ asymmetry metric we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal action $J_\phi$, excited by a single event (in our case, a satellite disk-crossing); (b) these distinct phase-spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase-spirals offer a global view of the MW's perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)-like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars within $R \lesssim 10$ kpc. Hence such perturbations have widespread influence which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW-Sgr interaction. While Sgr has not yet crossed the disk (currently, $z_{Sgr} \approx -6$ kpc, $v_{z,Sgr} \approx 210$ km/s), we demonstrate that the peak of the impact has already passed. Sgr's pull over the past 150 Myr creates a global $v_z$ signature with amplitude $\propto M_{Sgr}$, which might be detectable in future spectroscopic surveys.

We perform a systematic study of merging black hole (BH) binaries with compact star (CS) companions, including black hole--white dwarf (BH--WD), black hole--neutron star (BH--NS) and black hole--black hole (BH--BH) systems. Previous studies have shown that mass transfer stability and common envelope evolution can significantly affect the formation of merging BH--CS binaries through isolated binary evolution. With detailed binary evolution simulations, we obtain easy-to-use criteria for the occurrence of the common envelope phase in mass-transferring BH binaries with a nondegenerate donor, and incorporate into population synthesis calculations. To explore the impact of possible mass gap between NSs and BHs on the properties of merging BH--CS binary population, we adopt different supernova mechanisms involving the \textit{rapid}, \textit{delayed} and \textit{stochastic} prescriptions to deal with the compact remnant masses and the natal kicks. Our calculations show that there are $ \sim 10^{5} -10^{6}$ BH--CS binaries in the Milky Way, among which dozens are observable by future space-based gravitational wave detectors. We estimate that the local merger rate density of all BH--CS systems is $ \sim 60-200 \,\rm Gpc^{-3}yr^{-1}$. While there are no low-mass BHs formed via \textit{rapid} supernovae, both \textit{delayed} and \textit{stochastic} prescriptions predict that $ \sim 100\% $/$ \sim 70\% $/$ \sim 30\% $ of merging BH--WD/BH--NS/BH--BH binaries are likely to have BH components within the mass gap.

The first solids that form as a white dwarf (WD) starts to crystallize are expected to be greatly enriched in actinides. Previously [PRL 126, 1311010] we found that these solids might support a nuclear fission chain reaction that could ignite carbon burning and provide a new Type Ia supernova (SN Ia) mechanism involving an {\it isolated} WD. Here we explore this fission mechanism in more detail and calculate the final temperature and density after the chain reaction and discuss a number of open physics questions.

Primordial black hole (PBH) is a fascinating candidate for the origin of binary merger events observed by LIGO-Virgo collaboration. The spatial distribution of PBHs at formation is an important feature to estimate the merger rate. We investigate the clustering of PBHs formed by Affleck-Dine (AD) baryogenesis, where dense baryon bubbles collapse to form PBHs. We found that formed PBHs show a strong clustering due to the stochastic dynamics of the AD field. Including the clustering, we evaluate the merger rate and isocurvature perturbations of PBHs, which show significant deviations from those without clustering.

In this work, we studied the impact of galaxy morphology on photometric redshift (photo-$z$) probability density functions (PDFs). By including galaxy morphological parameters like the radius, axis-ratio, surface brightness and the S\'ersic index in addition to the $ugriz$ broadbands as input parameters, we used the machine learning photo-$z$ algorithm ANNz2 to train and test on galaxies from the Canada-France-Hawaii Telescope Stripe-82 (CS82) Survey. Metrics like the continuous ranked probability score (CRPS), probability integral transform (PIT), Bayesian odds parameter, and even the width and height of the PDFs were evaluated, and the results were compared when different number of input parameters were used during the training process. We find improvements in the CRPS and width of the PDFs when galaxy morphology has been added to the training, and the improvement is larger especially when the number of broadband magnitudes are lacking.

Gravitational waves are ripples in the fabric of space-time that travel at the speed of light. The detection of gravitational waves by LIGO is a major breakthrough in the field of astronomy. Deep Learning has revolutionized many industries including health care, finance and education. Deep Learning techniques have also been explored for detection of gravitational waves to overcome the drawbacks of traditional matched filtering method. However, in several researches, the training phase of neural network is very time consuming and hardware devices with large memory are required for the task. In order to reduce the extensive amount of hardware resources and time required in training a neural network for detecting gravitational waves, we made SpecGrav. We use 2D Convolutional Neural Network and spectrograms of gravitational waves embedded in noise to detect gravitational waves from binary black hole merger and binary neutron star merger. The training phase of our neural network was of about just 19 minutes on a 2GB GPU.

We present timing and spectral results of 2018 outburst of Cepheus X-4, observed twice by AstroSat at luminosity of $2.04 \times 10^{37}$ erg s$^{-1}$ and $1.02 \times 10^{37}$ erg s$^{-1}$ respectively. The light curves showed strong pulsation and co-related X-ray intensity variation in SXT (0.5--8.0 keV) and LAXPC (3--60 keV) energy bands. Spin-period and spin-down rate of the pulsar were determined from two observations as $65.35080\pm0.00014$ s , $(-2.10\pm0.8)\times10^{-12}$ Hz s$^{-1}$ at an epoch MJD 58301.61850 and $65.35290\pm0.00017$ s, $(-1.6\pm0.8)\times10^{-12}$ Hz s$^{-1}$ for an epoch MJD 58307.40211. Pulse-shape studies with AstroSat showed energy and intensity dependent variations. The pulsar showed an overall continuous spin-down, over 30 years at an average-rate of $(-2.455\pm0.004)\times10^{-14}$ Hz s$^{-1}$, attributed to propeller-effect in the subsonic-regime of the pulsar, in addition to variations during its outburst activities. Spectra between 0.7--55 keV energy band were well fitted by two continuum models, an absorbed compTT-model and an absorbed power-law with a Fermi-Dirac cutoff (FD-cutoff) model with a black-body. These were combined with an iron-emission line and a cyclotron absorption line. The prominent cyclotron resonance scattering features with a peak absorption energy of $30.48^{+0.33}_{-0.34}$ keV and $30.68^{+0.45}_{-0.44}$ keV for FD-cutoff-model and $30.46^{+0.32}_{-0.28}$ keV and $30.30^{+0.36}_{-0.34}$ keV for compTT-model were detected during two AstroSat observations. These when compared with earlier results, showed long term stability of its average value of $30.23 \pm 0.22$ keV. The pulsar showed pulse-phase as well as luminosity dependent variations in cyclotron-line energy and width and in plasma optical-depth of its spectral continuum.

We propose a model to explain the time delay between the peak of the optical and X-ray luminosity, \dt hereafter, in UV/optically-selected tidal disruption events (TDEs). The following picture explains the observed \dt in several TDEs as a consequence of the circularization and disk accretion processes as long as the sub-Eddington accretion. At the beginning of the circularization, the fallback debris is thermalized by the self-crossing shock caused by relativistic precession, providing the peak optical emission. During the circularization process, the mass fallback rate decreases with time to form a ring around the supermassive black hole (SMBH). The formation timescale corresponds to the circularization timescale of the most tightly bound debris, which is less than a year to several decades, depending mostly on the penetration factor, the circularization efficiency, and the black hole mass. The ring will subsequently evolve viscously over the viscous diffusion time. We find that it accretes onto the SMBH on a fraction of the viscous timescale, which is $2$ years for given typical parameters, leading to X-ray emission at late times. The resultant \dt\,is given by the sum of the circularization timescale and the accretion timescale and significantly decreases with increasing penetration factor to several to $\sim10$ years typically. Since the X-ray luminosity substantially decreases as the viewing angle between the normal to the disk plane and line-of-sight increases from $0^\circ$ to $90^\circ$, a low late-time X-ray luminosity can be explained by an edge-on view. We also discuss the super-Eddington accretion scenario, where \dt\,is dominated by the circularization timescale.

We critically discuss recent articles by S. Hoffmann and N. Vogt on historical novae and supernovae (SNe) as well as their list of `24 most promising events' `with rather high probability to be a nova' (Hoffmann et al. 2020). Their alleged positional accuracy of previously suggested historical nova/SN records is based on inhomogeneous datasets (Vogt et al. 2019), but then used for the nova search in Hoffmann et al. (2020). Their claim that previously only `point coordinates' for nova/SN candidates were published, is fabricated. Their estimate of expected nova detection rates is off by a factor of 10 due to miscalculation. They accept counterparts down to 4 to 7 mag at peak, which is against the consensus for the typical limit of naked-eye discovery. When they discuss previously suggested identifications of historical novae, which they all doubt, they do not present new facts (Hoffmann 2019). Their catalog of `24 most promising events' for novae (Hoffmann et al. 2020) neglects important recent literature (e.g. Pankenier et al. 2008 and Stephenson and Green 2009), the claimed methods are not followed, etc. At least half of their short-list candidates were and are to be considered comets. For many of the others, duration of more than one night and/or a precise position is missing and/or the sources were treated mistakenly. Two highlights, a fabricated SN AD 667/668 and a presumable recurrent nova in AD 891, are already rejected in detail in Neuhaeuser et al. (2021); in both cases, all evidence speaks in favor of comets. There remains only one reliable case, where close to one (possible) historically reported position, a nova shell was already found (AD 1437, Shara et al. 2017). Follow-up observations cannot be recommended.

How a black hole accretes matter and how this process is regulated are fundamental but unsolved questions in astrophysics. In transient black-hole binaries, a lot of mass stored in an accretion disk is suddenly drained to the central black hole because of thermal-viscous instability. This phenomenon is called an outburst and is observable at various wavelengths (Frank et al., 2002). During the outburst, the accretion structure in the vicinity of a black hole shows dramatical transitions from a geometrically-thick hot accretion flow to a geometrically-thin disk, and the transition is observed at X-ray wavelengths (Remillard, McClintock, 2006; Done et al., 2007). However, how that X-ray transition occurs remains a major unsolved problem (Dunn et al., 2008). Here we report extensive optical photometry during the 2018 outburst of ASASSN-18ey (MAXI J1820$+$070), a black-hole binary at a distance of 3.06 kpc (Tucker et al., 2018; Torres et al., 2019) containing a black hole and a donor star of less than one solar mass. We found optical large-amplitude periodic variations similar to superhumps which are well observed in a subclass of white-dwarf binaries (Kato et al., 2009). In addition, the start of the stage transition of the optical variations was observed 5 days earlier than the X-ray transition. This is naturally explained on the basis of our knowledge regarding white dwarf binaries as follows: propagation of the eccentricity inward in the disk makes an increase of the accretion rate in the outer disk, resulting in huge mass accretion to the black hole. Moreover, we provide the dynamical estimate of the binary mass ratio by using the optical periodic variations for the first time in transient black-hole binaries. This paper opens a new window to measure black-hole masses accurately by systematic optical time-series observations which can be performed even by amateur observers.

It is pointed out that MOND defines a fiducial specific angular momentum (SAM) for a galaxy of total (baryonic) mass $\mathcal{M}$: $j_M(\mathcal{M})\equiv\mathcal{M}^{3/4}(G^3/a_0)^{1/4}\approx 383(\mathcal{M}/10^{10}M_\odot)^{3/4}{\rm kpc~km/s}$. It plays important roles in disc-galaxy dynamics and evolution: It underlies scaling relations in virialized galaxies that involve their angular-momentum. I show that the disc SAM should be $j_D\approx[\langle r\rangle/r_M(\mathcal{M})]j_M(\mathcal{M})=[\Sigma_M/\langle \Sigma\rangle]^{1/2}j_M(\mathcal{M})$, with $\langle r\rangle$ the mean radius of the disc, $\langle \Sigma\rangle=\mathcal{M}/2\pi\langle r\rangle^2$ some mean surface density of the galaxy, $r_M=(\mathcal{M} G/a_0)^{1/2}$ is the MOND radius of the galaxy, and $\Sigma_M=a_0/2\pi G$ is the (universal) MOND surface density. So, e.g., for a fixed $\langle \Sigma\rangle$, $j_D\propto \mathcal{M}^{3/4}$, while for a fixed $\langle r\rangle$, $j_D\propto \mathcal{M}^{1/4}$. Furthermore, $j_M(\mathcal{M})$ is a reference predictor of the type of galaxy a protogalaxy will settle into, if it evolves in isolation: A protogalaxy of mass $\mathcal{M}$, and SAM $j\gg j_M(\mathcal{M})$ should settle into a low-surface-density disc -- with mean acceleration $\langle a\rangle/a_0\approx j_M/j\ll 1$. While a protogalaxy with $j\lesssim j_M(\mathcal{M})$ should end up with a disc of mass $\mathcal{M}_D\approx j\mathcal{M}/j_M(\mathcal{M})$, having a SAM $j_D\approx j_M(\mathcal{M})$, which is tantamount to $\langle a\rangle\approx a_0$ (i.e., at the `Freeman limit'); it should also develop a low-SAM bulge, taking up the rest of the mass $\mathcal{M}_B\approx\mathcal{M}-\mathcal{M}_D$.

We present infrared spectral indices (1.0-2.3 um) of Galactic late-type giants and red supergiants (RSGs). We used existing and new spectra obtained at resolution power R=2000 with SpeX on the IRTF telescope. While a large CO equivalent width (EW), at 2.29 um ([CO, 2.29]>45 AA) is a typical signature of RSGs later than spectral type M0, [CO] of K-type RSGs and giants are similar. In the [CO, 2.29] versus [Mg I, 1.71] diagram, RSGs of all spectral types can be distinguished from red giants, because the Mg I line weakens with increasing temperature and decreasing gravity. We find several lines that vary with luminosity, but not temperature: Si I (1.59 um), Sr (1.033 um), Fe+Cr+Si+CN (1.16 um), Fe+Ti (1.185 um), Fe+Ti (1.196 um), Ti+Ca (1.28 um), and Mn (1.29 um). Good markers of CN enhancement are the Fe+Si+CN line at 1.087 um and CN line at 1.093 um. Using these lines, at the resolution of SpeX, it is possible to separate RSGs and giants. Contaminant O-rich Mira and S-type AGBs are recognized by strong molecular features due to water vapor features, TiO band heads, and/or ZrO absorption. Among the 42 candidate RSGs that we observed, all but one were found to be late-types. 21 have EWs consistent with those of RSGs, 16 with those of O-rich Mira AGBs, and one with an S-type AGB. These infrared results open new, unexplored, potential for searches at low-resolution of RSGs in the highly obscured innermost regions of the Milky Way.

We present an analysis of the spatial clustering of 695 Ly$\alpha$-emitting galaxies (LAE) in the MUSE-Wide survey. All objects have spectroscopically confirmed redshifts in the range $3.3<z<6$. We employ the K-estimator of Adelberger et al. (2005), adapted and optimized for our sample. We also explore the standard two-point correlation function approach, which is however less suited for a pencil-beam survey such as ours. The results from both approaches are consistent. We parametrize the clustering properties by, (i) modelling the clustering signal with a power law (PL), and (ii) adopting a Halo Occupation Distribution (HOD) model. Applying HOD modeling, we infer a large-scale bias of $b_{\rm{HOD}}=2.80^{+0.38}_{-0.38}$ at a median redshift of the number of galaxy pairs $\langle z_{\rm pair}\rangle\simeq3.82$, while the PL analysis results in $b_{\rm{PL}}=3.03^{+1.51}_{-0.52}$ ($r_0=3.60^{+3.10}_{-0.90}\;h^{-1}$Mpc and $\gamma=1.30^{+0.36}_{-0.45}$). The implied typical dark matter halo (DMH) mass is $\log(M_{\rm{DMH}}/[h^{-1}\rm{M}_\odot])=11.34^{+0.23}_{-0.27}$. We study possible dependencies of the clustering signal on object properties by bisecting the sample into disjoint subsets, considering Ly$\alpha$ luminosity, UV absolute magnitude, Ly$\alpha$ equivalent width, and redshift as variables. We find a suggestive trend of more luminous Ly$\alpha$ emitters residing in more massive DMHs than their lower Ly$\alpha$ luminosity counterparts. We also compare our results to mock LAE catalogs based on a semi-analytic model of galaxy formation and find a stronger clustering signal than in our observed sample. By adopting a galaxy-conserving model we estimate that the LAEs in the MUSE-Wide survey will typically evolve into galaxies hosted by halos of $\log(M_{\rm{DMH}}/[h^{-1}\rm{M}_\odot])\approx13.5$ at redshift zero, suggesting that we observe the ancestors of present-day galaxy groups.

DA-type white dwarfs account for 80% of all white dwarfs and represent, for most of them, the ultimate outcome of the typical evolution of low-to-intermediate mass stars. Their internal chemical stratification is strongly marked by passed, often uncertain, stellar evolution processes that occurred during the helium (core and shell) burning phases, i.e., from the horizontal branch through AGB and post-AGB stages. Pulsating white dwarfs, in particular the "cool" DA-type ZZ Ceti variables, offer an outstanding opportunity to dig into these stars by fully exploiting their asteroseismic potential. With our most recent tools dedicated to that purpose, we show that a complete cartography of the stratification of the main constituents of a white dwarf can be inferred, leading in particular to strong constraints on the C/O core structure produced by the processes mentioned above. This opens up the way toward a systematic exploration of white-dwarf internal properties.

One of the most promising tracers of the Galactic potential in the halo region are stellar streams. However, individual stream fits can be limited by systematic biases. To study these individual stream systematics, we fit streams in Milky Way-like galaxies from FIRE cosmological galaxy formation simulations with an analytic gravitational potential by maximizing the clustering of stream stars in action space. We show that for coherent streams the quality of the constraints depends on the orbital phase of the observed stream stars, despite the fact that the phase information is discarded in action-clustering methods. Streams on intermediate phases give the most accurate results, whereas pericentre streams can be highly biased. This behaviour is tied to the amount of correlation present between positions and momenta in each stream's data: weak correlation in pericentre streams prohibits efficient differentiation between potentials, while strong correlation in intermediate streams promotes it. Although simultaneous fitting of multiple streams is generally prescribed as the remedy to combat individual stream biases, we find that combining multiple pericentric streams is not enough to yield a bias-free result. We finally show that adopting the two-component St\"ackel model does not fundamentally induce a biased mass estimate. With our full data set of two multi-wrap streams, we recovered the true rotation curve of the simulated galaxy within $12\%$ over the entire range of radii covered by our set of stars (10 - 176 kpc) and within $6.5\%$ between the 5 and 95-percentile distance range (23 - 109 kpc).

We performed a systematic search for X-ray bursts of the SGR J1935+2154 using the Fermi Gamma-ray Burst Monitor continuous data dated from Jan 2013 to July 2021. Seven major bursting phases, which consist of a total of 253 individual bursts, are identified. We further analyze the periodic properties of our sample using the Lomb-Scargle spectrum and a novel model (named "Simple Period Model") developed by ourselves. Two methods yield the same results in that those bursts exhibit a period of ~237 days with a ~58.6% duty cycle. Based on our analysis, we further predict two upcoming active windows of the X-ray bursts. As of July 7th, 2021, the beginning date of our first perdition has been confirmed by the ongoing X-ray activities of the SGR J1935+2154.

The CARMENES exoplanet survey of M dwarfs has obtained more than 18 000 spectra of 329 nearby M dwarfs over the past five years as part of its guaranteed time observations (GTO) program. We determine planet occurrence rates with the 71 stars from the GTO program for which we have more than 50 observations. We use injection-and-retrieval experiments on the radial-velocity (RV) time series to measure detection probabilities. We include 27 planets in 21 planetary systems in our analysis. We find 0.06+0.04-0.03 giant planets (100 M_Earth < M_pl sin i < 1000 M_Earth) per star in periods of up to 1000 d, but due to a selection bias this number could be up to a factor of five lower in the whole 329-star sample. The upper limit for hot Jupiters (orbital period of less than 10 d) is 0.03 planets per star, while the occurrence rate of planets with intermediate masses (10 M_Earth < M_pl sin i < 100 M_Earth) is 0.18+0.07-0.05 planets per star. Less massive planets with 1 M_Earth < M_pl sin i < 10 M_Earth are very abundant, with an estimated rate of 1.32+0.33-0.31 planets per star for periods of up to 100 d. When considering only late M dwarfs with masses M_star < 0.34 M_sol, planets more massive than 10 M_Earth become rare. Instead, low-mass planets with periods shorter than 10 d are significantly overabundant. For orbital periods shorter than 100 d, our results confirm the known stellar mass dependences from the Kepler survey: M dwarfs host fewer giant planets and at least two times more planets with M_pl sin i < 10 M_Earth than G-type stars. In contrast to previous results, planets around our sample of very low-mass stars have a higher occurrence rate in short-period orbits of less than 10 d. Our results demonstrate the need to take into account host star masses in planet formation models.

We perform a study on the impacts on Ceres and Vesta. We aim to determine the size-frequency distribution (SFD) of impactors and to identify and quantify the contribution of each source region, as well as the craters produced and fragments ejected in these impact events. We used a multipart collisional evolution model of the MB called ACDC (Asteroid Collisions and Dynamic Computation) that simulates the collisional evolution of the MB, which is split into six regions (namely Inner, Middle, Pristine, Outer, Cybele, and High-Inclination belts). Furthermore, it includes the Yarkovsky effect as a dynamical remotion mechanism. The six regions of the MB provide, to a greater or lesser extent, impactors on Ceres and Vesta. The Outer belt is the main source of impactors smaller than 10 km on Ceres, providing more than half of the impacts, while the Middle belt is the secondary source. On Vesta, the relative impactor contribution of the Inner, Middle, and Outer belts is almost even. We are able to reproduce the craters larger than 100 km in Vesta and identify two large depressions identified in Ceres as impact craters: one called Vendimia Planitia of 900 km and a second one of 570 km. As an outcome of these impacts, Ceres and Vesta eject fragments into the MB. We obtain fragmentation rates of tens of fragments larger than 1 m per year for both bodies, to tens of fragments larger than 100 m per million years for Vesta and a factor of 4 greater for Ceres. We find that hundreds of bodies larger than 10 km should have been ejected from Ceres and Vesta during their history

We consider the oblique fragmentation of the current layer as a result of the thermal instability described in Ledentsov (Sol. Phys. 296, 74, 2021a). It is shown that the fragmentation transverse to the current is a natural feature of the model. The fragmentation tilt does not exceed a few degrees for realistic preflare parameters of the coronal plasma. As a consequence, oblique fragmentation generally does not have a strong impact on the simulation results, however, extreme changes can reach an order of magnitude. Thus, oblique fragmentation can lead to a decrease in the estimate of the spatial period of the location of elementary energy release in solar flares to 0.1-1 Mm instead of 1-10 Mm obtained earlier.

We have done photometric calibration of the 60 cm Nedeljkovi\'{c} telescope equipped with FLI PL 230 CCD camera, mounted at the Astronomical Station Vidojevica (Serbia), using standard stars from the Landolt's catalog. We have imaged 31 fields of standard stars using Johnson's $BVRI$ filters during three nights in August 2019. We have measured both extinction and color correction. Relating our calibrated magnitudes to the magnitudes of the standard stars from the Landolt's catalog, we have achieved accuracy of 2\%-5\% for the $BVRI$ magnitudes.

Solar Energetic Particle events (SEPs) are among the most dangerous transient phenomena of solar activity. As hazardous radiation, SEPs may affect the health of astronauts in outer space and adversely impact current and future space exploration. In this paper, we consider the problem of daily prediction of Solar Proton Events (SPEs) based on the characteristics of the magnetic fields in solar Active Regions (ARs), preceding soft X-ray and proton fluxes, and statistics of solar radio bursts. The machine learning (ML) algorithm uses an artificial neural network of custom architecture designed for whole-Sun input. The predictions of the ML model are compared with the SWPC NOAA operational forecasts of SPEs. Our preliminary results indicate that 1) for the AR-based predictions, it is necessary to take into account ARs at the western limb and on the far side of the Sun; 2) characteristics of the preceding proton flux represent the most valuable input for prediction; 3) daily median characteristics of ARs and the counts of type II, III, and IV radio bursts may be excluded from the forecast without performance loss; and 4) ML-based forecasts outperform SWPC NOAA forecasts in situations in which missing SPE events is very undesirable. The introduced approach indicates the possibility of developing robust "all-clear" SPE forecasts by employing machine learning methods.

Highlights are presented about the science to be done with SKA. as well as state of the art science already done today with its precursors (MeerKAT, ASKAP) and pathfinders (LOFAR, NenuFAR), with accent on the expected breakthroughs.

Extremely low-mass white dwarfs (ELM WDs) are the result of binary evolution in which a low-mass donor star is stripped by its companion leaving behind a helium-core white dwarf. We explore the formation of ELM WDs in binary systems considering the Convection And Rotation Boosted magnetic braking treatment. Our evolutionary sequences were calculated using the MESA code, with initial masses of 1.0 and 1.2 Msun (donor), and 1.4 (accretor), compatible with low mass X-ray binaries (LMXB) systems. We obtain ELM models in the range 0.15 to 0.27 Msun from a broad range of initial orbital periods, 1 to 25 d. The bifurcation period, where the initial period is equal to the final period, ranges from 20 to 25 days. In addition to LMXBs, we show that ultra-compact X-ray binaries (UCXB) and wide-orbit binary millisecond pulsars can also be formed. The relation between mass and orbital period obtained is compatible with the observational data from He white dwarf companions to pulsars.

Zahn's widely-used model for turbulent mixing induced by rotational shear has recently been validated (with some caveats) in non-rotating shear flows. It is not clear, however, whether his model remains valid in the presence of rotation, even though this was its original purpose. Furthermore, new instabilities arise in rotating fluids, such as the Goldreich-Schubert-Fricke (GSF) instability. Which instability dominates when more than one can be excited, and how they influence each other, were open questions that this paper answers. To do so, we use direct numerical simulations of diffusive stratified shear flows in a rotating triply-periodic Cartesian domain located at the equator of a star. We find that either the GSF instability or the shear instability tends to take over the other in controlling the system, suggesting that stellar evolution models only need to have a mixing prescription for each individual instability, together with a criterion to determine which one dominates. However, we also find that it is not always easy to predict which instability "wins" for given input parameters, because the diffusive shear instability is subcritical, and only takes place if there is a finite-amplitude turbulence ``primer'' to seed it. Interestingly, we find that the GSF instability can in some cases play the role of this primer, thereby providing a pathway to excite the subcritical shear instability. This can also drive relaxation oscillations, that may be observable. We conclude by proposing a new model for mixing in the equatorial regions of stellar radiative zones due to differential rotation.

The majority of planetary nebulae (PNe) show axisymmetric morphologies, whose causes are not well understood. In this work, we present spatially resolved kinematic observations of 14 Galactic PNe surrounding Wolf-Rayet ([WR]) and weak emission-line stars ($wels$) based on the H$\alpha$ and [N II] emission taken with the Wide Field Spectrograph on the ANU 2.3-m telescope. Velocity-resolved channel maps and position--velocity diagrams, together with archival Hubble Space Telescope ($HST$) and ground-based images, are employed to construct three-dimensional morpho-kinematic models of 12 objects using the program SHAPE. Our results indicate that these 12 PNe have elliptical morphologies with either open or closed outer ends. Kinematic maps also illustrate on-sky orientations of elliptically symmetric morphologies of the interior shells in NGC 6578 and NGC 6629, and the compact ($\leq 6$ arcsec) PNe Pe1-1, M3-15, M1-25, Hen2-142, and NGC 6567, in agreement with the high-resolution $HST$ images containing morphological details. Point-symmetric knots in Hb4 exhibit deceleration with distance from the nebular center that could be due to shock collisions with the ambient medium. Velocity dispersion maps of Pe1-1 disclose point-symmetric knots similar to those in Hb4. Collimated outflows are also visible in the position--velocity diagrams of M3-30, M1-32, M3-15, and K2-16, which are reconstructed by tenuous prolate ellipsoids extending upwardly from thick toroidal shells in our models.

Mass-radius relations of homogeneous cold spheres are obtained for six solid materials commonly found in terrestrial planets. An additional degeneracy in the (exo-)planets' profiles is discussed together with their properties concluded from our findings in the framework of Palatini $f(\mathcal R)$ gravity. Moreover, a new test of gravity has been proposed: The results presented here will allow to test and to constrain models of gravity by the use of seismic data acquired from earthquakes and marsquakes.

Intermediate/Extreme mass ratio inspiral (IMRI/EMRI) system provides a good tool to test the nature of gravity in strong field. We construct the self-force and use the self-force method to generate accurate waveform templates for IMRIS/EMRIs on quasi-elliptical orbits in Brans-Dicke theory. The extra monopole and dipole emissions in Brans-Dicke theory accelerate the orbital decay, so the observations of gravitational waves may place stronger constraint on Brans-Dicke theory. With a two-year observations of gravitational waves emitted from IMRIs/EMRIs with LISA, we can get the most stringent constraint on the Brans-Dicke coupling parameter $\omega_0>10^5$.

A pseudo Nambu-Goldstone boson (such as an axion-like particle) is a theoretically well-motivated inflaton as it features a naturally flat potential (natural inflation). This is because Goldstone's theorem protects its potential from sizable quantum corrections. Such corrections, however, generically generates an $R^2$ term in the action, which leads to another inflaton candidate because of the equivalence between the $R^2$ term and a scalar field, the scalaron, with a quasi flat potential (Starobinsky inflation). Here it is investigated a new multifield scenario in which both the scalaron and a pseudo Nambu-Goldstone boson are active (natural-scalaron inflation). For generality, also a non-minimal coupling is included, which is shown to emerge from microscopic theories. It is demonstrated that a robust inflationary attractor is present even when the masses of the two inflatons are comparable. Moreover, the presence of the scalaron allows to satisfy all observational bounds in a large region of the parameter space, unlike what happens in pure-natural inflation.

It is often said that asymmetric dark matter is light compared to typical weakly interacting massive particles. Here we point out a simple scheme with a neutrino portal and $\mathcal{O}(60 \text{ GeV})$ asymmetric dark matter which may be ''added'' to any standard baryogenesis scenario. The dark sector contains a copy of the Standard Model gauge group, as well as (at least) one matter family, Higgs, and right-handed neutrino. After baryogenesis, some lepton asymmetry is transferred to the dark sector through the neutrino portal where dark sphalerons convert it into a dark baryon asymmetry. Dark hadrons form asymmetric dark matter and may be directly detected due to the vector portal. Surprisingly, even dark anti-neutrons may be directly detected if they have a sizeable electric dipole moment. The dark photons visibly decay at current and future experiments which probe complementary parameter space to dark matter direct detection searches. Exotic Higgs decays are excellent signals at future $e^+ e^-$ Higgs factories.

For a stationary, axisymmetric, asymptotically flat, ultra-compact [$i.e.$ containing light-rings (LRs)] object, with a $\mathbb{Z}_2$ north-south symmetry fixing an equatorial plane, we establish that the structure of timelike circular orbits (TCOs) in the vicinity of the equatorial LRs, for either rotation direction, depends exclusively on the stability of the LRs. Thus, an unstable LR delimits a region of unstable TCOs (no TCOs) radially above (below) it; a stable LR delimits a region of stable TCOs (no TCOs) radially below (above) it. Corollaries are discussed for both horizonless ultra-compact objects and black holes. We illustrate these results with a variety of exotic stars examples and non-Kerr black holes, for which we also compute the efficiency associated with converting gravitational energy into radiation by a material particle falling under an adiabatic sequence of TCOs. For most objects studied, it is possible to obtain efficiencies larger than the maximal efficiency of Kerr black holes, $i.e.$ larger than $42\%$.

We consider a cosmological scenario in which the very early Universe experienced a transient epoch of matter domination due to the formation of a large population of primordial black holes (PBHs) with masses $M \lesssim 10^{9}\,\textrm{g}$, that evaporate before Big Bang nucleosynthesis. In this context, Hawking radiation would be a non-thermal mechanism to produce a cosmic background of axion-like particles (ALPs). We assume the minimal scenario in which these ALPs couple only with photons. In the case of ultralight ALPs ($m_a \lesssim 10^{-9}\,\textrm{eV}$) the cosmic magnetic fields might trigger ALP-photon conversions, while for masses $m_a \gtrsim 10\,\textrm{eV}$ spontaneous ALP decay in photon pairs would be effective. We investigate the impact of these mechanisms on the cosmic X-ray background, on the excess in X-ray luminosity in Galaxy Clusters, and on the process of cosmic reionization.

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

We perform a complete next-to-leading order calculation of the non-Abelian electric field correlator in a SU($N_c$) plasma, which encodes properties of the plasma relevant for heavy particle bound state formation and dissociation. The calculation is carried out in the real-time formalism of thermal field theory and includes both vacuum and finite temperature contributions. By working in the $R_\xi$ gauge, we explicitly show the results are gauge independent, infrared and collinear safe. The previous results on the renormalization of the electric field correlator are also confirmed. Our next-to-leading order calculation can be directly applied to any dipole singlet-adjoint transition of heavy particle pairs. For example, it can be used to describe dissociation and (re)generation of heavy quarkonia inside the quark-gluon plasma well below the melting temperature, as well as heavy dark matter pairs (or charged co-annihilating partners) in the early universe.

We explore the effect of forcing on the linear shear flow or plane Couette flow, which is also the background flow in the very small region of the Keplerian accretion disk. We show that depending on the strength of forcing and boundary conditions suitable for the systems under consideration, the background plane shear flow and, hence, accretion disk velocity profile modifies to parabolic flow, which is plane Poiseuille flow or Couette-Poiseuille flow, depending on the frame of reference. In the presence of rotation, plane Poiseuille flow becomes unstable at a smaller Reynolds number under pure vertical as well as threedimensional perturbations. Hence, while rotation stabilizes plane Couette flow, the same destabilizes plane Poiseuille flow faster and forced local accretion disk. Depending on the various factors, when local linear shear flow becomes Poiseuille flow in the shearing box due to the presence of extra force, the flow becomes unstable even for the Keplerian rotation and hence turbulence will pop in there. This helps in resolving a long standing problem of sub-critical transition to turbulence in hydrodynamic accretion disks and laboratory plane Couette flow.