Papers for Wednesday, Sep 21 2022

Current prescriptions for supernova natal kicks in rapid binary population synthesis simulations are based on fits of simple functions to single pulsar velocity data. We explore a new parameterization of natal kicks received by neutron stars in isolated and binary systems developed by Mandel & M\"uller, which is based on 1D and 3D supernova simulations and accounts for the physical correlations between progenitor properties, remnant mass, and the kick velocity. We constrain two free parameters in this model using very long baseline interferometry velocity measurements of Galactic single pulsars. We find that the inferred values of natal kick parameters do not differ significantly between single and binary evolution scenarios. The best-fit values of these parameters are $v_{\rm ns} = 520$ km s$^{-1}$ for the scaling pre-factor for neutron star kicks, and $\sigma_{\rm ns}=0.3$ for the fractional stochastic scatter in the kick velocities.

Long Gamma-Ray Bursts (LGRBs), the most powerful events in the Universe, are generated by jets that emerge from dying massive stars. Highly beamed geometry and immense energy make jets promising gravitational wave (GW) sources. However, their sub-Hertz GW emission is outside of ground based GW detector (LIGO) frequency band. Using a 3D general-relativistic magnetohydrodynamic simulation of a dying star, we show that jets inflate a turbulent, energetic bubble-cocoon that emits strong quasi-spherical GW emission within the LIGO band, $ 0.1-0.6 $ kHz, over the characteristic jet activity timescale, $ \approx 10-100 $ s. This is the first non-inspiral GW source detectable by LIGO out to hundreds of Mpc, with $ \approx 0.1 - 10 $ detectable events expected during LIGO observing run O4. These GWs are likely accompanied by detectable energetic core-collapse supernova and cocoon electromagnetic emission, making jetted stellar explosions promising multi-messenger sources.

A critical early stage for the origin of life on Earth may have involved the production of hydrogen cyanide (HCN) in a reducing, predominantly H$_2$ atmosphere. HCN is crucial for the origin of life as it is a possible precursor to several biomolecules that make up RNA and proteins including nucleobases, nucleotides, amino acids, and ribose. In this work, we perform an in depth experimental and theoretical investigation of HCN production in reducing atmospheric conditions (89-95% H$_2$) possibly representing the earliest stages of the Hadean eon, ~4.5-4.3 billion years ago. We make use of cold plasma discharges - a laboratory analog to shortwave UV radiation - to simulate HCN production in the upper layers of the atmosphere for CH$_4$ abundances ranging from 0.1-6.5%. We then combine experimental mass spectrum measurements with our theoretical plasma models to estimate the HCN concentrations produced in our experiments. We find that upper atmospheric HCN production scales linearly with CH$_4$ abundance with the relation [HCN] = 0.13 $\pm$ 0.01[CH$_4$]. Concentrations of HCN near the surface of the Hadean Earth are expected to be about 2-3 orders of magnitude lower. The addition of 1% water to our experiments results in a ~50% reduction in HCN production. We find that four reactions are primarily responsible for HCN production in our experiments: (i) $^4$N + CH$_3$ -> H$_2$CN + H -> HCN + H$_2$, (ii) $^4$N + CH -> CN + H followed by CN + CH$_4$ -> HCN + CH$_3$, (iii) C$_2$H$_4$ + $^4$N -> HCN + CH$_3$, and (iv) $^4$N + $^3$CH$_2$ -> HCN + H. The most prebiotically favorable Hadean atmosphere would have been very rich in CH$_4$ (> 5%), and as a result of greenhouse effects the surface would be likely very hot. In such a prebiotic scenario, it may have been important to incorporate HCN into organic hazes that could later release biomolecules and precursors into the first ponds.

Early-type dwarfs (ETDs) are the end points of the evolution of low-mass galaxies whose gas supply has been extinguished. The cessation of star-formation lays bare the ancient stellar populations. A wealth of information is stored in the colours, magnitudes, metallicities and abundances of resolved stars of the dwarf spheroidal and ultra-faint galaxies around the Milky Way, allowing their chemistry and stellar populations to be studied in great detail. Here, we summarize our current understanding, which has advanced rapidly over the last decade thanks to the flourishing of large-scale astrometric, photometric and spectroscopic surveys. We emphasise that the primeval stellar populations in the ETDs provide a unique laboratory to study the physical conditions on small scales at epochs beyond z=2. We highlight the observed diversity of star-formation and chemical enrichment histories in nearby dwarfs. These data can not yet be fully deciphered to reveal the key processes in the dwarf evolution but the first successful attempts have been made to pin down the sites of heavy element production.

Well-resolved galaxy clusters often show a large-scale quasi-spiral structure in deprojected density $\rho$ and temperature $T$ fields, delineated by a tangential discontinuity known as a cold front, superimposed on a universal radial entropy profile with a linear $K(r)\propto T\rho^{-2/3}\propto r$ adiabat. We show that a spiral structure provides a natural quasi-stationary solution for the mixed intracluster medium (ICM), introducing a modest pressure spiral that confines the locally buoyant or heavy plasma phases. The solution persists in the presence of uniform or differential rotation, and can accommodate both an inflow and an outflow. Hydrodynamic adiabatic simulations with perturbations that deposit angular momentum and mix the plasma thus asymptote to a self-similar spiral structure. We find similar spirals in Eulerian and Lagrangian simulations of 2D and 3D, merger and offset, clusters. The discontinuity surface is given in spherical coordinates $\{r,\theta,\phi\}$ by $\phi\propto \Phi(r)$, where $\Phi$ is the gravitational potential, combining a trailing spiral in the equatorial ($\theta=\pi/2$) plane and semicircles perpendicular to the plane, in resemblance of a snail shell. A local convective instability can develop between spiral windings, driving a modified global instability in sublinear $K(r)$ regions; evolved spirals thus imprint the observed $K\propto r$ onto the ICM even after they dissipate. The spiral structure brings hot and cold phases to close proximity, suggesting that the observed fast outflows could sustain the structure even in the presence of radiative cooling.

We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC~253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central $\sim$250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and non-spherical outflow geometry. Our derived temperatures and densities yield few-Myr cooling times in the nuclear region, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates on the mass outflow rates which are $2.8\:M_{\odot}/\rm{yr}$ in the northern outflow and $3.2\:M_{\odot}/\rm{yr}$ in the southern outflow. Additionally, we detect emission from charge exchange and find it has a significant contribution ($20-42$\%) to the total broad-band ($0.5-7$~keV) X-ray emission in the central and southern regions of the outflow.

While dwarf galaxies observed in the field are overwhelmingly star-forming, dwarf galaxies in environments as dense or denser than the Milky Way are overwhelmingly quenched. In this paper, we explore quenching in the lower density environment of the Small-Magellanic-Cloud-mass galaxy NGC 3109 ($\text{M}_* \sim 10^8 \, \text{M}_\odot$), which hosts two known dwarf satellite galaxies (Antlia and Antlia B), both of which are HI deficient compared to similar galaxies in the field and have recently stopped forming stars. Using a new semi-analytic model in concert with the measured star formation histories and gas masses of the two dwarf satellite galaxies, we show that they could not have been quenched solely by direct ram pressure stripping of their interstellar media, as is common in denser environments. Instead, we find that separation of the satellites from pristine gas inflows, coupled with stellar-feedback-driven outflows from the satellites (jointly referred to as the starvation quenching model), can quench the satellites on timescales consistent with their likely infall times into NGC 3109's halo. It is currently believed that starvation is caused by "weak" ram pressure that prevents low-density, weakly-bound gas from being accreted onto the dwarf satellite, but cannot directly remove the denser interstellar medium. This suggests that star-formation-driven outflows serve two purposes in quenching satellites in low-mass environments: outflows from the host form a low-density circumgalactic medium that cannot directly strip the interstellar media from its satellites, but is sufficient to remove loosely-bound gaseous outflows from the dwarf satellites driven by their own star formation.

One of the major science goals of the Square Kilometre Array (SKA) is to understand the role played by atomic hydrogen (HI) gas in the evolution of galaxies throughout cosmic time. The hyperfine transition line of the hydrogen atom at 21-cm is one of the best tools to detect and study the properties of HI gas associated with galaxies. In this article, we review our current understanding of HI gas and its relationship with galaxies through observations of the 21-cm line both in emission and absorption. In addition, we provide an overview of the HI science that will be possible with SKA and its pre-cursors and pathfinders, i.e. HI 21-cm emission and absorption studies of galaxies from nearby to high redshifts that will trace various processes governing galaxy evolution.

Transiting planet systems offer the best opportunity to measure the masses and radii of a large sample of planets and their host stars. However, relative photometry and radial velocity measurements alone only constrain the density of the host star. Thus, there is a one-parameter degeneracy in the mass and radius of the host star, and by extension the planet. Several theoretical, semi-empirical, and nearly empirical methods have been used to break this degeneracy and independently measure the mass and radius of the host star and planets(s). As we approach an era of few percent precisions on some of these properties, it is critical to assess whether these different methods are providing accuracies that are of the same order, or better than, the stated statistical precisions. We investigate the differences in the planet parameter estimates inferred when using the Torres empirical relations, YY isochrones, MIST isochrones, and a nearly-direct empirical measurement of the radius of the host star using its spectral energy distribution, effective temperature, and \textit{Gaia} parallax. We focus our analysis on modelling KELT-15b, a fairly typical hot Jupiter, using each of these methods. We globally model TESS photometry, optical-to-NIR flux densities of the host star, and \textit{Gaia} parallaxes, in conjunction with extant KELT ground-based follow-up photometric and radial velocity measurements. We find systematic differences in several of the inferred parameters of the KELT-15 system when using different methods, including a $\sim 6\%$ ($\sim 2\sigma$) difference in the inferred stellar and planetary radii between the MIST isochrones and SED fitting.

We report the detection of the ground state rotational emission of ammonia, ortho-NH$_3$ $(J_K=1_0\rightarrow0_0)$ in a gravitationally lensed, intrinsically hyperluminous, star-bursting galaxy at $z=2.6$. The integrated line profile is consistent with other molecular and atomic emission lines which have resolved kinematics well-modelled by a 5 kpc-diametre rotating disc. This implies that the gas responsible for NH$_3$ emission is broadly tracing the global molecular reservoir, but likely distributed in pockets of high density ($n\gtrsim5\times10^4$ cm$^{-3}$). With a luminosity of $2.8\times10^{6}$ $L_\odot$, the NH$_3$ emission represents $2.5\times10^{-7}$ of the total infrared luminosity of the galaxy, comparable to the ratio observed in the Kleinmann-Low nebula in Orion and consistent with sites of massive star formation in the Milky Way. If $L_{\rm NH_3}/L_{\rm IR}$ serves as a proxy for the 'mode' of star formation, this hints that the nature of star formation in extreme starbursts in the early Universe is similar to that of Galactic star-forming regions, with a large fraction of the cold interstellar medium in this state, plausibly driven by a storm of violent disc instabilities in the gas-dominated disc. This supports the 'full of Orions' picture of star formation in the most extreme galaxies seen close to the peak epoch of stellar mass assembly.

Asymmetry in the spatially integrated, 1D HI global profiles of galaxies can inform us on both internal (e.g. outflows) and external (e.g. mergers, tidal interactions, ram pressure stripping) processes that shape galaxy evolution. Understanding which of these primarily drive HI profile asymmetry is of particular interest. In the lead-up to SKA pathfinder and SKA HI emission surveys, hydrodynamical simulations have proved to be a useful resource for such studies. Here we present the methodology behind, as well as first results, of ASymba: Asymmetries in HI of Simba galaxies, the first time this simulation suite has been used for this type of study. We generate mock observations of the HI content of these galaxies and calculate the profile asymmetries using three different methods. We find that $M_{\rm HI}$ has the strongest correlation with all asymmetry measures, with weaker correlations also found with the number of mergers a galaxy has undergone, and gas and galaxy rotation. We also find good agreement with the xGASS sample, in that galaxies with highly asymmetric profiles tend to have lower HI gas fractions than galaxies with symmetric profiles, and additionally find the same holds in sSFR parameter space. For low HI mass galaxies, it is difficult to distinguish between asymmetric and symmetric galaxies, but this becomes achievable in the high HI mass population. These results showcase the potential of ASymba and provide the groundwork for further studies, including comparison to upcoming large HI emission surveys.

Narrow-line Seyfert 1 (NLS1) galaxies are unevolved active galactic nuclei (AGN) that exist predominantly in spiral galaxies. However, mostly due to the small number of sources studied, it has been under debate whether also the hosts of jetted NLS1 galaxies, a particular subclass of these sources hosting a relativistic jet, are disk-like, or elliptical, as the hosts of more powerful jetted AGN. We studied the host morphologies of 14 NLS1 galaxies, 11 of which have been detected at 37 GHz indicating that these sources harbour relativistic jets. The J- and Ks-band data used in this study were obtained with the Nordic Optical Telescope (NOT). We performed the photometric decomposition of the host galaxy using the band that gave a better fit, and additionally, created colour maps of all sources that had both a J- and a Ks-band observation. We were able to successfully model 12 sources, nine of which most likely have disk-like morphology. Of the remaining sources, one source could possibly be hosted either in a disk-like or a dwarf galaxy, and in two cases the results are inconclusive. Only one of our sources shows clear signs of interaction, but the colour maps of most of our sources hint at ample dust in the nuclei, possibly indicating earlier minor mergers, that can go unnoticed due to the limited resolution of these observations. Our results further support disk-like galaxies as the predominant host type of jetted NLS1 galaxies. Most importantly, with the number of modelled hosts of jetted NLS1s now exceeding 50, with only a few elliptical hosts, it seems to be safe to conclude that also disk-like galaxies are able to launch and maintain relativistic jets, and that the traditional jet paradigm stating that only massive elliptical galaxies are capable of hosting relativistic jets is severely outdated.

Massive binary stars play a crucial role in many astrophysical fields. Investigating the statistical properties of massive binary stars is essential to trace the formation of massive stars and constrain the evolution of stellar populations. However, no consensus has been achieved on the statistical properties of massive binary stars, mainly due to the lack of a large and homogeneous sample of spectroscopic observations. We study the intrinsic binary fraction $f_{\rm b}^{\rm in}$ and distributions of mass ratio $f(q)$ and orbital period $f(P)$ of early-type stars (comprised of O-, B-, and A-type stars) and investigate their dependences on effective temperature $T_{\rm eff}$, stellar metallicity [M/H], and the projection velocity $v\sin{i}$, based on the homogeneous spectroscopic sample from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data Release Eight (DR8). We found that $f_{\rm b}^{\rm in}$ increases with increasing $T_\mathrm{eff}$. The binary fraction is positively correlated with metallicity for spectra in the sample. Over all the $v\sin{i}$ values we considered, the $f_{\rm b}^{\rm in}$ have constant values of $\sim$50\%. It seems that the binary population is relatively evenly distributed over a wide range of $v\sin{i}$ values, while the whole sample shows that most of the stars are concentrated at low values of $v\sin{i}$ (probably from strong wind and magnetic braking of single massive stars) and at high values of $v\sin{i}$ (likely from the merging of binary stars). Stellar evolution and binary interaction may be partly responsible for this.There are no correlations found between $\pi$($\gamma$) and $T_{\rm eff}$, nor for $\pi$($\gamma$) and [M/H]. The uncertainties of the distribution decrease toward a larger sample size with higher observational cadence.

We present a robust representation of the chemical and thermal structure in the galaxy group NGC 5813 using archival, deep X-ray observations, and employing a multi-temperature spectral model based on up to date atomic line emission databases. The selection of our target is motivated by the fact that NGC 5813 has a very relaxed morphology, making it a promising candidate for the study of the AGN feedback's influence in the intra-group medium (IGrM). Our results showcase a prominent, extended distribution of cool gas along the group's NE-SW direction, correlating with the direction along which the supermassive black hole in the group's central galaxy is known to interact with the IGrM. Our analysis indicates gas being uplifted from the group's centre as the probable origin of the cool gas, although alternative scenarios, such as in-situ cooling can not be explicitly ruled out. Regarding the chemical structure of the IGrM, and unlike previous findings in massive clusters, we find no evidence for recent metal transport by jets/lobes from the central AGN. Instead, elemental abundances remain near Solar on average across the group. The distribution of elements appears to be independent of galactocentric radius, azimuth and the thermodynamics of the gas, suggesting that the IGrM has been efficiently mixed. The large scale uniformity of the abundance distribution implies the presence of complex dynamical processes in NGC 5813, despite its overall relaxed morphology. Past events of extreme AGN feedback or sloshing could be the primary mechanisms behind this.

The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$CO, C$^{18}$O and optically thin C$^{17}$O $J$=2-1 lines, and use additional high angular resolution Atacama Large Millimeter Array millimeter continuum and CO data to construct a representative model of LkCa 15. The transitions that constrain the gas mass and carbon abundance most are C$^{17}$O 2-1, N${_2}$H$^+$ 3-2 and HD 1-0. Using these three molecules we find that the gas mass in the LkCa 15 disk is $M_\mathrm{g}=0.01 ^{+0.01}_{-0.004} M_{\odot}$, a factor of six lower than estimated before. The carbon abundance is C/H = ($3 \pm 1.5) \times10^{-5}$, implying a moderate depletion of elemental carbon by a factor of 3-9. All other analyzed transitions also agree with these numbers, within a modeling uncertainty of a factor of two. Using the resolved \ce{C2H} image we find a C/O ratio of $\sim$1, which is consistent with literature values of H$_2$O depletion in this disk. The lack of severe carbon depletion in the LkCa 15 disk is consistent with the young age of the disk, but contrasts with the higher depletions seen in older cold transition disks. Combining optically thin CO isotopologue lines with N$_2$H$^+$ is promising to break the degeneracy between gas mass and CO abundance. The moderate level of depletion for this source with a cold, but young disk, suggests that long carbon transformation timescales contribute to the evolutionary trend seen in the level of carbon depletion among disk populations, rather than evolving temperature effects and presence of dust traps alone.

We perform non-hydrodynamical 2.5D simulations to study the dynamics of material above accretion disk based on the disk radiation pressure acting on dust. We assume a super-accreting underlying disk with the accretion rate of 10 times the Eddington rate with central black hole mass ranging from $10^7$ up to $10^9 M_{\odot}$. Such high accretion rates are characteristic for extreme sources. We show that for high accretors radiatively dust-driving mechanism based on FRADO model always leads to a massive outflow from the disk surface, and the failed wind develops only at larger radii. The outflow rate strongly depends on the black hole mass, and in optically-thick energy-driven solution can exceed the accretion rate for masses larger than $10^ 8 M_{\odot}$ but momentum-driven outflow does not exceed the accretion rate even for super-Eddington accretion, therefore not violating the adopted stationarity of the disk. However, even in this case the outflow from the disk implies a strong mechanical feedback.

Ground-level particle detection is now a well-established approach to TeV {\gamma}-ray astronomy. Detection of Cherenkov light produced in water-filled detection units is a proven and cost-effective approach. Here we discuss the optimization of the units towards the future Southern Wide-field Gamma-ray Observatory (SWGO). We investigate a configuration in which each Water Cherenkov Detector (WCD) unit in the array comprises two chambers with black or reflective walls and a single PMT in each chamber. A shallow lower chamber with a PMT facing downwards enables muon tagging and the identification of hadron-induced air showers, which are the primary source of background in {\gamma}-ray astronomy. We investigate how {\gamma}/hadron separation power and achievable angular resolution depend on the geometry and wall reflectivity of the detector units in this configuration. We find that excellent angular resolution, background rejection power and low-energy response are achievable in this double-layer configuration, with some reflective surfaces in both chambers.

We describe the first results from the All-sky BRIght, Complete Quasar Survey (AllBRICQS), which aims to discover the last remaining optically bright quasars. We present 116 quasars (105 newly identified) having |b| > 10deg and Gaia magnitudes brighter than B_P =16.5 or R_P =16 mag (plus another four at slightly fainter magnitudes), which span a redshift range of z = 0.07 - 2.45. The AllBRICQS sources have been selected by combining data from the Gaia and WISE all-sky satellite missions, and we successfully identify quasars not flagged as candidates by Gaia Data Release 3. We expect the completeness to be approximately 96% within our magnitude and latitude limits, while the preliminary results indicate a selection purity of approximately 97%. The optical spectroscopy used for source classification will also enable detailed quasar characterisation, including black hole mass measurements and identification of foreground absorption systems. The AllBRICQS sources will greatly enhance the number of quasars available for high-signal-to-noise follow-up with present and future facilities.

One of the key processes driving galaxy evolution during the Cosmic Dawn is supernova feedback. This likely helps regulate star formation inside of galaxies, but it can also drive winds that influence the large-scale intergalactic medium. Here, we present a simple semi-analytic model of supernova-driven galactic winds and explore the contributions of different phases of galaxy evolution to metal enrichment in the high-redshift (z > 6) Universe. We show that models calibrated to the observed galaxy luminosity function at z~6-8 have filling factors ~1% at z~6 and ~0.1% at z~12, with different star formation prescriptions providing about an order of magnitude uncertainty. Despite the small fraction of space filled by winds, these scenarios still provide more than enough enriched volume to explain the observed abundance of metal-line absorbers in quasar spectra at z > 5. We also consider enrichment through winds driven by Pop III star formation in minihaloes. We find that these can dominate the total filling factor at z > 10 and even compete with winds from normal galaxies at z~6, at least in terms of the total enriched volume. But these regions have much lower overall metallicities, because each one is generated by a small burst of star formation. Finally, we show that Compton cooling of these supernova-driven winds at z > 6 has only a small effect on the cosmic microwave background.

In order to understand grain-surface chemistry, one must have a good understanding of the reaction rate parameters. For diffusion-based reactions, these parameters are binding energies of the reacting species. However, attempts to estimate these values from grain-surface abundances using Bayesian inference are inhibited by a lack of enough sufficiently constraining data. In this work, we use the Massive Optimised Parameter Estimation and Data (MOPED) compression algorithm to determine which species should be prioritised for future ice observations to better constrain molecular binding energies. Using the results from this algorithm, we make recommendations for which species future observations should focus on.

The kinematics of the Milky Way (MW) and M31, the dominant galaxies in the Local Group (LG), can be used to estimate the LG total mass. New results on the M31 proper motion have recently been used to improve that estimate. Those results are based on kinematic priors that are sometimes guided and evaluated using cosmological N-body simulations. However, the kinematic properties of simulated LG analogues could be biased due to the effective power spectrum truncation induced by the small size of the parent simulation. Here we explore the dependence of LG kinematics on the simulation box size to argue that cosmological simulations need a box size on the order of 1 Gpc in order to claim convergence on the LG kinematic properties. Using a large enough simulation, we find M31 tangential and radial velocities relative to the MW to be in the range $v_{\mathrm {tan}}=105^{+94}_{-59} $ km/s and $v_{\mathrm {rad}}=-108^{+68}_{-81}$ km/s, respectively. This study highlights that LG kinematics derived from N-body simulations have to be carefully interpreted taking into account the size of the parent simulation.

Delta ($\delta$)-spots are active regions (ARs) in which positive and negative umbrae share a penumbra. They are known to be the source of strong flares. We introduce a new quantity, the degree of $\delta$ (Do$\delta$), to measure the fraction of umbral flux participating in the $\delta$-configuration and to isolate the dynamics of the magnetic knot, i.e. adjacent umbrae in the $\delta$-configuration. Using Helioseismic and Magnetic Imager data, we analyze 19 $\delta$-spots and 11 $\beta$-spots in detail, and 120 $\delta$-spots in less detail. We find that $\delta$-regions are not in a $\delta$-configuration for the entire time but spend 55$\%$ of their observed time as $\delta$-spots with an average, maximum Do$\delta$ of 72$\%$. Compared to $\beta$-spots, $\delta$-spots have 2.6$\times$ the maximum umbral flux, 1.9$\times$ the flux emergence rate, 2.6$\times$ the rotation, and 72$\times$ the flare energy. On average, the magnetic knots rotate 17$^{\circ}$ day$^{-1}$ while the $\beta$-spots rotate 2$^{\circ}$ day$^{-1}$. Approximately 72$\%$ of the magnetic knots present anti-Hale or anti-Joy tilts, contrasting starkly with only 9$\%$ of the $\beta$-spots. A positive correlation exists between $\phi_{Do\delta}$ and the flare energy emitted by that region. The $\delta$-spots obey the hemispheric current helicity rule 64$\%$ of the time. 84$\%$ of the $\delta$-spots are formed by single flux emergence events and 58$\%$ have a quadrupolar magnetic configuration. The $\delta$-spot characteristics are consistent with the formation mechanism signatures as follows: 42$\%$ with the kink instability or Sigma effect, 32$\%$ with multi-segment buoyancy, 16$\%$ with collisions and two active regions that are unclassified but consistent with a rising O-ring.

The orbital distribution of transneptunian objects (TNOs) in the distant Kuiper Belt (with semimajor axes beyond the 2:1 resonance, roughly $a$=50-100 au) provides constraints on the dynamical history of the outer solar system. Recent studies show two striking features of this region: 1) a very large population of objects in distant mean-motion resonances with Neptune, and 2) the existence of a substantial detached population (non-resonant objects largely decoupled from Neptune). Neptune migration models are able to implant some resonant and detached objects during the planet migration era, but many fail to match a variety of aspects of the orbital distribution. In this work, we report simulations carried out using an improved version of the GPU-based code GLISSE, following 100,000 test particles per simulation in parallel while handling their planetary close encounters. We demonstrate for the first time that a 2 Earth-mass rogue planet temporarily present during planet formation can abundantly populate both the distant resonances and the detached populations, surprisingly even without planetary migration. We show how weak encounters with the rogue greatly increase the efficiency of filling the resonances, while also dislodging TNOs out of resonance once they reach high perihelia. The rogue's secular gravitational influence simultaneously generates numerous detached objects observed at all semimajor axes. These results suggest that the early presence of additional planet(s) reproduces the observed TNO orbital structure in the distant Kuiper Belt.

Asteroseismology has opened a window on the internal physics of thousands of stars, by relating oscillation spectra properties to the internal physics of stars. Mode identification, namely the process of associating a measured oscillation frequency to the corresponding mode geometry and properties, is the cornerstone of this analysis of seismic spectra. In rapidly rotating stars this identification is a challenging task that remains incomplete, as modes assume complex geometries and regular patterns in frequencies get scrambled under the influence of the Coriolis force and centrifugal flattening. In this article, I will first discuss the various classes of mode geometries that emerge in rapidly-rotating stars and the related frequency and period patterns, as predicted by ray dynamics, complete (non-)adiabatic calculations, or using the traditional approximation of rotation. These patterns scale with structural quantities and help us derive crucial constraints on the structure and evolution of these stars. I will summarize the amazing progress accomplished over the last few years for the deciphering of gravity-mode pulsator oscillation spectra, and recent developments based on machine-learning classification techniques to distinguish oscillation modes and pattern analysis strategies that let us access the underlying physics of pressure-mode pulsators. These approaches pave the way to ensemble asteroseismology of classical pulsators. Finally, I will highlight how these recent progress can be combined to improve forward seismic modelling. I will focus on the example of Rasalhague, a well-known rapid rotator, to illustrate the process and the needed advances to obtain \`a-la-carte modelling of such stars.

The surface of airless bodies like asteroids in the Solar System are known to be affected by space weathering. Experiments simulating space weathering are essential for studying the effects of this process on meteorite samples, but the problem is that the time spent to reproduce space weathering in these experiments is billions of times shorter than the actual phenomenon. In December 2010, the T-type asteroid 596 Scheila underwent a collision with a few-tens-of-meters impactor. A decade later, there is an opportunity to study how the surface layer of this asteroid is being altered by space weathering after the impact. To do so, we performed visible spectrophotometric and near-infrared spectroscopic observations of 596 Scheila. The acquired spectrum is consistent with those observed shortly after the 2010 impact event within the observational uncertainty range. This indicates that the surface color of dark asteroids is not noticeably changed by space weathering over a 10-year period. This study is the first to investigate color changes due to space weathering on an actual asteroid surface in the Solar System. Considering that fresh layers are regularly created on asteroid surfaces by collisions, we suggest a genetic link between D/T-type and dark (low albedo) X-complex asteroids and very red objects such as 269 Justitia, 732 Tjilaki (and 203 Pompeja). New observations show that 203 Pompeja has a X-type-like surface, with some local surface areas exhibiting a very red spectrum.

The aim of this study is to construct a simple stellar model with non-uniform polytropic index. We find that the Emden equation cannot deal with the polytrope gas sphere with non-uniform polytropic index in a real star, and then we construct a realistic stellar model. The key point is that we should solve the two independent equations for density and pressure due to the essence of polytropic relation, but not the Emden equation which combines the hydrostatic balance and polytropic relation. We take the Sun for a computational example to find that this simple model yields quite a good result compared to the MESA code. The advantage of this simple model lies in its much simpler equation of state than that in the standard stellar model.

We report the discovery of a triple-giant-planet system around an evolved star HD 184010 (HR 7421, HIP 96016). This discovery is based on observations from Okayama Planet Search Program, a precise radial velocity survey, undertaken at Okayama Astrophysical Observatory between 2004 April and 2021 June. The star is K0 type and located at beginning of the red-giant branch. It has a mass of $1.35_{-0.21}^{+0.19} M_{\odot}$, a radius of $4.86_{-0.49}^{+0.55} R_{\odot}$, and a surface gravity $\log g$ of $3.18_{-0.07}^{+0.08}$. The planetary system is composed of three giant planets in a compact configuration: The planets have minimum masses of $M_{\rm{b}}\sin i = 0.31_{-0.04}^{+0.03} M_{\rm{J}}$, $M_{\rm{c}}\sin i = 0.30_{-0.05}^{+0.04} M_{\rm{J}}$, and $M_{\rm{d}}\sin i = 0.45_{-0.06}^{+0.04} M_{\rm{J}}$, and orbital periods of $P_{\rm{b}}=286.6_{-0.7}^{+2.4}\ \rm{d}$, $P_{\rm{c}}=484.3_{-3.5}^{+5.5}\ \rm{d}$, and $P_{\rm{d}}=836.4_{-8.4}^{+8.4}\ \rm{d}$, respectively, which are derived from a triple Keplerian orbital fit to three sets of radial velocity data. The ratio of orbital periods are close to $P_{\rm{d}}:P_{\rm{c}}:P_{\rm{b}} \sim 21:12:7$, which means the period ratios between neighboring planets are both lower than $2:1$. The dynamical stability analysis reveals that the planets should have near-circular orbits. The system could remain stable over 1 Gyr, initialized from co-planar orbits, low eccentricities ($e=0.05$), and planet masses equal to the minimum mass derived from the best-fit circular orbit fitting. Besides, the planets are not likely in mean motion resonance. HD 184010 system is unique: it is the first system discovered to have a highly evolved star ($\log g < 3.5$ cgs) and more than two giant planets all with intermediate orbital periods ($10^2\ \rm{d} < P < 10^3\ \rm{d}$).

Laser frequency combs are fast becoming critical to reaching the highest radial velocity precisions. One shortcoming is the highly variable brightness of the comb lines across the spectrum (up to 4-5 orders of magnitude). This can result in some lines saturating while others are at low signal and lost in the noise. Losing lines to either of these effects reduces the precision and hence effectiveness of the comb. In addition, the brightness of the comb lines can vary with time which could drive comb lines with initially reasonable SNR's into the two regimes described above. To mitigate these two effects, laser frequency combs use optical flattener's. Flattener's are typically bulk optic setups that disperse the comb light with a grating, and then use a spatial light modulator to control the amplitude across the spectrum before recombining the light into another single mode fiber and sending it to the spectrograph. These setups can be large (small bench top), expensive (several hundred thousand dollars) and have limited stability. To address these issues, we have developed an all-photonic spectrum flattener on a chip. The device is constructed from optical waveguides on a SiN chip. The light from the laser frequency comb's output optical fiber can be directly connected to the chip, where the light is first dispersed using an arrayed waveguide grating. To control the brightness of each channel, the light is passed through a Mach-Zehnder interferometer before being recombined with a second arrayed waveguide grating. Thermo-optic phase modulators are used in each channel before recombination to path length match the channels as needed. Here we present the results from our first generation prototype. The device operates from 1400-1800 nm (covering the H band), with 20, 20 nm wide channels.

Previous numerical simulations have shown that strong winds can be produced in the hot accretion flows around black holes. Most of those studies focus only on the region close to the central black hole, therefore it is unclear whether the wind production stops at large radii around Bondi radius. Bu et al. 2016 studied the hot accretion flow around the Bondi radius in the presence of nuclear star gravity. They find that when the nuclear stars gravity is important/comparable to the black hole gravity, winds can not be produced around the Bondi radius. However, for some galaxies, the nuclear stars gravity around Bondi radius may not be strong. In this case, whether winds can be produced around Bondi radius is not clear. We study the hot accretion flow around Bondi radius with and without thermal conduction by performing hydrodynamical simulations. We use the virtual particles trajectory method to study whether winds exist based on the simulation data. Our numerical results show that in the absence of nuclear stars gravity, winds can be produced around Bondi radius, which causes the mass inflow rate decreasing inwards. We confirm the results of Yuan et al. which indicates this is due to the mass loss of gas via wind rather convectional motions.

Samuil Kaplan (1921-1978) was a productive and famous astrophysicist. He was affiliated with a number of scientific centers in different cities of former Soviet Union. The earliest 13 years of his career, namely in the 1948-1961 years, he worked in Lviv University in Ukraine (then it was called the Ukrainian Soviet Socialist Republic). In the present paper, the Lviv period of his life and scientific activity is described on the basis of archival materials and his published studies. Kaplan arrived in Lviv in June 1948, at the same month when he obtained the degree of Candidate of science. He was a head of the astrophysics sector at the Astronomical Observatory of the University, was a professor of department for theoretical physics as well as the founder and head of a station for optical observations of artificial satellites of Earth. He was active in the organization of the astronomical observational site outside of the city. During the years in Lviv, Kaplan wrote more than 80 articles and 3 monographs in 9 areas. The focus of his interests at that time was on stability of circular orbits in the Schwarzschild field, on white dwarf theory, on space gas dynamics, and cosmic plasma physics, and turbulence, on acceleration of cosmic rays, on physics of interstellar medium, on physics and evolution of stars, on cosmology and gravitation, and on optical observations of Earth artificial satellites. Some of his results are fundamental for development of theory in these fields as well as of observational techniques. The complete bibliography of his works published during the Lviv period is presented. Respective scientific achievements of Samuil Kaplan are reviewed in the light of the current state of research in these areas.

We present high angular-resolution measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward two galaxy clusters, RCS J2319+0038 at z=0.9 and HSC J0947-0119 at z=1.1, by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. They are supplemented with available Chandra X-ray data, optical data taken by Hyper Suprime-Cam on Subaru, and millimeter-wave SZE data from the Atacama Cosmology Telescope. Taking into account departures from spherical symmetry, we have reconstructed non-parametrically the inner pressure profile of two clusters as well as electron temperature and density profiles for RCS J2319+0038. This is one of the first such measurements for an individual cluster at $z \gtrsim 0.9$. We find that the inner pressure profile of both clusters is much shallower than that of local cool-core clusters. Our results consistently suggest that RCS J2319+0038 hosts a weak cool core, where radiative cooling is less significant than in local cool cores. On the other hand, HSC J0947-0119 exhibits an even shallower pressure profile than RCS J2319+0038 and is more likely a non-cool-core cluster. The SZE centroid position is offset by more than 140 $h_{70}^{-1}$kpc from the peaks of galaxy distribution in HSC J0947-0119, suggesting a stronger influence of mergers in this cluster. We conclude that these distant clusters are at a very early stage of developing the cool cores typically found in clusters at lower redshifts.

Electromagnetic radiation from human activities, known as man-made Radio Frequency Interference (RFI), adversely affects radio astronomy observations. In the vicinity of the Upgraded Giant Metrewave Radio Telescope (uGMRT) array, the sparking on power lines is the major cause of interference at observing frequencies less than 800 MHz. A real-time broadband RFI detection and filtering system is implemented as part of the uGMRT wideband signal processing backend to mitigate the effect of broadband RFI. Performance analysis techniques used for testing and commissioning the system for observations in the beamformer and correlator modes of the uGMRT are presented. The concept and implementation of recording simultaneous unfiltered and filtered data along with data analysis and interpretation is illustrated using an example. For the beamformer mode, spectrogram, single spectral channel, and its Fourier transform is used for performance analysis whereas, in the correlator mode, the cross-correlation function, closure phase, and visibilities from the simultaneously recorded unfiltered and filtered is carried out. These techniques are used for testing the performance of the broadband RFI filter and releasing it for uGMRT users.

We consider strong gravitational lensing by nearby stars. Using our wave-optical treatment of lensing by a compact mass, we study Einstein rings that may form around such stellar lenses. These large and bright rings are resolvable by existing instruments. Such lensing events take place in hours or days, with peak light amplification lasting several minutes. Many such events may be predicted using the Gaia astrometric catalogue. Serendipitous discoveries are also possible. Fortuitous alignments can be used to confirm or discover and study exoplanets. For lenses that have dense stellar regions in their background, these events may occur annually or more often, warranting their continuous or recurrent monitoring. Resolved imaging and spectroscopy of the evolving morphology of an Einstein ring offers knowledge about both the lens and the source. The angular size of the Einstein ring amounts to a direct measurement of the lens mass. The changing orientation of the major and minor images of the source offers astrometric information. The event duration helps determine the source's size. The sky position of planetary lensing events constrains the planet's orbit. Spectroscopy of the ring allows for direct investigations of the source. The frequency and predictability of these events and the wealth of information that can be obtained by imaging motivate observational campaigns using existing facilities or new instruments dedicated to the search and study of Einstein rings that form around nearby stars. As a specific example, we consider a predicted 2028 lensing of a red giant by $\alpha$ Centauri A.

In this brief report we discuss how continuous changes on the physical parameters that determine the weather conditions may lead to long term climate variability. This variability of the weather patterns are a response to continuous random short period weather excitations that are imprinted in the ocean-atmosphere-cryosphere-land system, the Earth System. Given that Earth System is, in the Anthropocene, dominated by the human action, it responds to the intensity and the rate of change of the humankind activities. Thus, we argue, in the context of a specific model of the Earth System, that this rate of change may admit a chaotic-type behaviour.

Kinetic inductance detectors (KID) have great potential in astronomical observation, such as searching for exoplanets, because of their low noise, fast response and photon counting characteristics. In this paper, we present the design process and simulation results of a microstrip line coupled KIDs array for near-infrared astronomical observation. Compared with coplanar waveguide (CPW) feedlines, microstrip feedlines do not require air bridges, which simplify fabrication process. In the design part, we mainly focus on the impedance transforming networks, the KID structure, and the frequency crosstalk simulations. The test array has a total of 104 resonators with 8 rows and 13 columns, which ranges from 4.899~GHz to 6.194~GHz. The pitch size is about 200~$\mu$m and the frequency crosstalk is less than 50~kHz in simulation.

A pilot project has been proceeded to map 1 deg$^2$ on the Galactic plane for radio recombination lines (RRLs) using the Five hundred meter Aperture Spherical Telescope (FAST). The motivation is to verify the techniques and reliabilities for a large-scale Galactic plane RRL survey with FAST aiming to investigate the ionized environment in the Galaxy. The data shows that the bandpass of the FAST 19 beam L-band is severely affected by radio frequency interferences (RFIs) and standing wave ripples, which can hardly be corrected by traditional low order polynomials. In this paper, we investigate a series of penalized least square (PLS) based baseline correction methods for radio astronomical spectra that usually contain weak signals with high level of noise. Three promising penalized least squares based methods, AsLS, arPLS, and asPLS are evaluated. Adopting their advantages, a modified method named rrlPLS is developed to optimize the baseline fitting to our RRL spectra. To check their effectiveness, the four methods are tested by simulations and further verified using observed data sets. It turns out that the rrlPLS method, with optimized parameter $\lambda = 2 \times 10^8$ , reveals the most sensitive and reliable emission features in the RRL map. By injecting artificial line profiles into the real data cube, a further evaluation of profile distortion is conducted for rrlPLS. Comparing to simulated signals, the processed lines with low signal-to-noise ratio are less affected, of which the uncertainties are mainly caused by the rms noise. The rrlPLS method will be applied for baseline correction in future data processing pipeline of FAST RRL survey. Configured with proper parameters, the rrlPLS technique verified in this work may also be used for other spectroscopy projects.

Context: Rapid rotation is a common feature of early-type stars but which remains a challenge for the models. The understanding of its effect on stellar evolution is however imperative to interpret the observed properties of numerous stars. Aims: We wish to bring more observational constraints on the properties of fast rotating stars, especially on their oscillation modes. Methods: We focus on the nearby star Altair which is known as a very rapidly rotating star with an equatorial velocity estimated recently at 313 km/s. We observed this star with the high-resolution spectropolarimeter Neo-Narval during six nights, with one night of interruption, in September 2020. Results: We detect significant line profile variations on the mean line profile of the spectra. Their time-frequency analysis shows that these variations are induced by gravito-inertial waves propagating at Altair's surface with azimuthal wavenumbers of order $m=10-15$. With a preliminary computation of the eigenspectrum using the most recent concordance model of Altair we can give a first modelling of the observed waves. Conclusions: Altair was known as the brightest $\delta$ Scuti star. We now see that it is the brightest hybrid oscillating star with excited gravito-inertial waves and acoustic waves. Clearly, more observations and more advanced models are needed to explain the observations in greater details

The unidentified TeV source HESS J1702-420 has recently been proposed as a new hadronic PeVatron candidate, based on the discovery of a small-scale emission sub-region with extremely hard gamma-ray spectrum up to 100 TeV (named HESS J1702-420A). Given the difficulty to discriminate between a hadronic or leptonic origin of the TeV emission, based on the H.E.S.S. measurement alone, we opted for a multi-wavelength approach. A deep X-ray observation was carried out using the XMM-Newton satellite, with the goal of probing a possible association with a hidden leptonic accelerator. No evidence of a clear counterpart for HESS J1702-420A was found in the X-ray data. After excluding an association with all nearby X-ray point sources, we derived strict upper limits on the diffuse X-ray emission and average magnetic field in the HESS J1702-420A region. We additionally report the serendipitous discovery of a new extended X-ray source, whose association with HESS J1702-420A is not obvious but cannot be ruled out either. A set of scripts dedicated to the multi-wavelength modeling of X-ray and gamma-ray data, based on Gammapy, Naima and Xspec, was developed in the context of this work and is made publicly available along with this paper.

Pressure profiles are sensitive probes of the thermodynamic conditions and the internal structure of galaxy clusters. The intra-cluster gas resides in hydrostatic equilibrium within the Dark Matter gravitational potential. However, this equilibrium may be perturbed, e.g. as a consequence of thermal energy losses, feedback and non-thermal pressure supports. Accurate measures of the gas pressure over the cosmic times are crucial to constrain the cluster evolution as well as the contribution of astrophysical processes. In this work we presented a novel algorithm to derive the pressure profiles of galaxy clusters from the Sunyaev-Zeldovich (SZ) signal measured on a combination of Planck and South Pole Telescope (SPT) observations. The synergy of the two instruments made it possible to track the profiles on a wide range of spatial scales. We exploited the sensitivity to the larger scales of the Planck High-Frequency Instrument to observe the faint peripheries, and the higher spatial resolution of SPT to solve the innermost regions. We developed a two-step pipeline to take advantage of the specifications of each instrument. We first performed a component separation on the two data-sets separately to remove the background (CMB) and foreground (galactic emission) contaminants. Then we jointly fitted a parametric pressure profile model on a combination of Planck and SPT data. We validated our technique on a sample of 6 CHEX-MATE clusters detected by SPT. We compare the results of the SZ analysis with profiles derived from X-ray observations with XMM-Newton. We find an excellent agreement between these two independent probes of the gas pressure structure.

Polarized dust emission is a key tracer in the study of interstellar medium and of star formation. The observed polarization, however, is a product of magnetic field structure, dust grain properties and grain alignment efficiency, as well as their variations in the line of sight, making it difficult to interpret polarization unambiguously. The comparison of polarimetry at multiple wavelengths is a possible way of mitigating this problem. We use data from HAWC+/SOFIA and from SCUBA-2/POL-2 (from the BISTRO survey) to analyse the NGC 2071 molecular cloud at 154, 214 and 850 $\mu$m. The polarization angle changes significantly with wavelength over part of NGC 2071, suggesting a change in magnetic field morphology on the line of sight as each wavelength best traces different dust populations. Other possible explanations are the existence of more than one polarization mechanism in the cloud or scattering from very large grains. The observed change of polarization fraction with wavelength, and the 214-to-154 $\mu$m polarization ratio in particular, are difficult to reproduce with current dust models under the assumption of uniform alignment efficiency. We also show that the standard procedure of using monochromatic intensity as a proxy for column density may produce spurious results at HAWC+ wavelengths. Using both long-wavelength (POL-2, 850 $\mu$m) and short-wavelength (HAWC+, $\lesssim 200\, \mu$m) polarimetry is key in obtaining these results. This study clearly shows the importance of multi-wavelength polarimetry at submillimeter bands to understand the dust properties of molecular clouds and the relationship between magnetic field and star formation.

We present the star formation activity around the emission nebula Sh2-112. At a distance of $\sim2.1$~kpc, this \ion{H}{2} complex, itself 3~pc in radius, is illuminated by the massive star (O8\,V) BD$+$45\,3216. The associated molecular cloud extends in angular scales of $2\fdg0\times0\fdg83$, corresponding to linear sizes of 73~pc by 30~pc, along the Galactic longitude. The high-resolution ($30\arcsec$) extinction map reveals a chain of dust clumps aligned with the filament-like structure with an average extinction of $A_{V} \sim 2.78$~mag, varying up to a maximum of $\sim17$~mag. Our analysis led to identification of a rich population ($\sim 500$) of young (average age of $\sim 1$~Myr) stars, plus a numerous number ($\sim 350$) of H$\alpha$ emitters, spatially correlated with the filamentary clouds. Located near the edge of the cloud, the luminous star BD$+$45\,3216 has created an arc-like pattern as the ionizing radiation encounters the dense gas, forming a blister-shaped morphology. We found three distinct young stellar groups, all coincident with relatively dense parts of the cloud complex, signifying ongoing star formation. Moreover, the cloud filament (excitation temperature $\sim 10$~K) traced by the CO isotopologues and extending nearly $\sim 80$~pc is devoid of ionized gas except at the dense cores (excitation temperature $\sim$ 28--32~K) wherein significant ionized emission excited by OB stars (dynamical age $\sim$ 0.18--1.0~Myr) pertains. The radial velocity is dynamic (median $\sim -3.65$~km~s$^{-1}$) along the main filament, increasing from Galactic east to west, features mass flow to form the massive stars/clusters at the central hubs.

Radio emission from 23 tidal disruption event host galaxies were searched in the 144 MHz LOFAR-LoTSS2 images. Three host galaxies are detected with diffuse radio emission, which can be interpreted as either galactic synchrotron radio emission or diffuse radio halo in dense galactic environments. Non-detection of (transient) radio emission in majority of galaxies could be due to self-absorption of radio emission associated with the tidal disruption event. The detected objects are AT2018iih, RBS 1032, and NGC 5905. Non-detection for Sw 1644+57 is also reported.

In the multi-scale view of the star formation process the material flows from large molecular clouds down to clumps and cores. In this paradigm it is still unclear if it is gravity or turbulence that drives the observed supersonic non-thermal motions during the collapse, in particular in high-mass regions, and at which scales gravity becomes eventually dominant over the turbulence of the interstellar medium. To investigate this problem we have combined the dynamics of a sample of 70 micron-quiet clumps, selected to cover a wide range of masses and surface densities, with the dynamics of the parent filaments in which they are embedded. We observe a continuous interplay between turbulence and gravity, where the former creates structures at all scales and the latter takes the lead when a critical value of the surface density is reached, Sigma_th = 0.1 g cm^-2. In the densest filaments this transition can occur at the parsec, or even larger scales, leading to a global collapse of the whole region and most likely to the formation of the massive objects.

The ExoClock project has been created with the aim of increasing the efficiency of the Ariel mission. It will achieve this by continuously monitoring and updating the ephemerides of Ariel candidates over an extended period, in order to produce a consistent catalogue of reliable and precise ephemerides. This work presents a homogenous catalogue of updated ephemerides for 450 planets, generated by the integration of $\sim$18000 data points from multiple sources. These sources include observations from ground-based telescopes (ExoClock network and ETD), mid-time values from the literature and light-curves from space telescopes (Kepler/K2 and TESS). With all the above, we manage to collect observations for half of the post-discovery years (median), with data that have a median uncertainty less than one minute. In comparison with literature, the ephemerides generated by the project are more precise and less biased. More than 40\% of the initial literature ephemerides had to be updated to reach the goals of the project, as they were either of low precision or drifting. Moreover, the integrated approach of the project enables both the monitoring of the majority of the Ariel candidates (95\%), and also the identification of missing data. The dedicated ExoClock network effectively supports this task by contributing additional observations when a gap in the data is identified. These results highlight the need for continuous monitoring to increase the observing coverage of the candidate planets. Finally, the extended observing coverage of planets allows us to detect trends (TTVs - Transit Timing Variations) for a sample of 19 planets. All products, data, and codes used in this work are open and accessible to the wider scientific community.

In pre-recombination early dark energy (EDE) resolutions of the Hubble tension, the rise of Hubble constant value $H_0$ is usually accompanied with the exacerbation of so-called $S_8$ tension. Inspired by the swampland conjecture, we investigate what if a fraction $f_*$ of dark matter is coupled to EDE, $m_{cdm}\sim \exp{(-c{|\Delta\phi_{ede}|\over M_{pl}})}$ with $c\sim {\cal O}(1)$. We perform the MCMC analysis for the relevant EDE models with PlanckCMB, BAO, Pantheon and SH0ES dataset, as well as DES-Y1 data, and find that such a fraction helps to alleviate the $S_8$ tension. However, though $c\gtrsim 0.1$ is allowed for a very small $f_*$, which suggests that a small fraction of dark matter has ever faded with EDE, $c\sim0$ is also consistent.

Axions provide a natural and well-motivated dark matter candidate, with the capability to convert directly to photons in the presence of an electromagnetic field. A particularly compelling observational target is the conversion of dark matter axions into photons in the magnetospheres of highly magnetised neutron stars, which is expected to produce a narrow spectral peak centred at the frequency of the axion mass. We point the MeerKAT radio telescope towards the isolated neutron star J0806.4$-$4123 for $10$-hours of observation and obtain the radio spectra in the frequency range $769$-$1051$ MHz. By modelling the conversion process of infalling axion dark matter (DM) we then compare these spectra to theoretical expectations for a given choice of axion parameters. Whilst finding no signal above $5\sigma$ in the data, we provide a unique constraint on the Primakoff coupling of axion DM, $g_{{\rm a}\gamma\gamma}\lesssim 9.3 \times 10^{-12}\,{\rm GeV}^{-1}$ at the $95\%$ confidence level, in the mass range $3.18$-$4.35\,\mu$eV. This result serves the strongest constraint in the axion mass range $4.20$-$4.35\,\mu$eV.

We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of sky that is dominated by complex and bright Galactic cirrus and it has been poorly-studied since its discovery more than 20 years ago. The tidal feature was revealed in a deep map of resolved red giant branch stars that we constructed using data from our Subaru Hyper Suprime-Cam survey of the M81 Group. It has an average surface brightness of $\mu_g \sim 32$ mag arcsec$^{-2}$ and with our current imagery can be traced for over a degree ($\gtrsim$ 60 kpc at the distance of F8D1) to the North-East of the galaxy, pointing in the direction of the dwarf spiral NGC 2976, as well as M81. We revisit the main body properties of F8D1 using deep multiband imagery acquired with MegaCam on CFHT and measure effective radii of 1.7-1.9 kpc and central surface brightnesses of $\mu_{0} \sim 24.7-25.7$ mag. Assuming a symmetric feature on the other side of the galaxy, we calculate that $30-36$% of F8D1's present-day luminosity is contained in the stream. We argue that the most likely origin of F8D1's disruption is a recent close passage to M81, which would have stripped its gas and quenched its star formation. As the only UDG that is close enough to allow studies at extremely low surface brightness, and the first to be unambiguously linked to tidal stripping and heating, F8D1 is of particular importance. Many or most UDGs could be the result of similar processes, with the most obvious signatures of this lurking below current detection limits.

The Einstein Probe (EP) is a mission designed to monitor the sky in the soft X-ray band. It will perform systematic surveys and characterisation of high-energy transients and monitoring of variable objects at unprecedented sensitivity and monitoring cadences. It has a large instantaneous field-of-view (3,600 sq. deg.), that is realised via the lobster-eye micro-pore X-ray focusing optics. EP also carries a conventional X-ray focusing telescope with a larger effective area to perform followup observations and precise positioning of newly-discovered transients. Alerts for transient objects will be issued publicly and timely. The scientific goals of EP are concerned with discovering faint, distant or rare types of high-energy transients and variable sources. Within the confines of a modest-sized mission, EP will cover a wide range of scientific topics, from the nearby to high-redshift Universe. The Einstein Probe is a mission of the Chinese Academy of Sciences, and also an international collaborative project. This paper presents the background, scientific objectives, and the mission design including the micro-pore optics and CMOS technologies adopted, the instruments and their expected performance, and the mission profile. The development status of the project is also presented.

We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy-galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering auto-correlation and galaxy-galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in $\Lambda$CDM and $w$CDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlations functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the auto-correlation of the lens bins, indicating systematics other than magnification may be the cause of the discrepancy. We show adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns.

MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA.

One of the major discoveries of asteroseismology is the signature of a strong extraction of angular momentum (AM) in the radiative zones of stars across the entire Hertzsprung-Russell diagram, resulting in weak core-to-surface rotation contrasts. Despite all efforts, a consistent AM transport theory, which reproduces both the internal rotation and mixing probed thanks to the seismology of stars, remains one of the major open problems in modern stellar astrophysics. A possible key ingredient to figure out this puzzle is magnetic field with its various possible topologies. Among them, strong axisymmetric toroidal fields, which are subject to the so-called Tayler MHD instability, could play a major role. They could trigger a dynamo action in radiative layers while the resulting magnetic torque allows an efficient transport of AM. But is it possible to detect signatures of these deep toroidal magnetic fields? The only way to answer this question is asteroseismology and the best laboratories of study are intermediate-mass and massive stars because of their external radiative envelope. Since most of these are rapid rotators during their main-sequence, we have to study stellar pulsations propagating in stably stratified, rotating, and potentially strongly magnetised radiative zones. For that, we generalise the traditional approximation of rotation, which provides in its classic version a flexible treatment of the adiabatic propagation of gravito-inertial modes, by taking simultaneously general axisymmetric differential rotation and toroidal magnetic fields into account. Using this new non-perturbative formalism, we derive the asymptotic properties of magneto-gravito-inertial modes and we explore the different possible field configurations. We found that the magnetic effects should be detectable for equatorial fields using high-precision asteroseismic data.

Recent observations of substructures such as dust gaps and dust rings in protoplanetary discs have highlighted the importance of including dust into purely gaseous disc models. At the same time, computational difficulties arise with the standard models of simulating the dust and gas separately. These include the cost of accurately simulating the interactions between well coupled dust and gas and issues of dust concentration in areas below resolution of the gas phase. We test a single fluid approach, that incorporates the terminal velocity approximation valid for small particles, which can overcome these difficulties, through modification of FARGO3D. We compare this single fluid model with a multifluid model for a variety of planet masses. We find differences in the dust density distribution in all cases. For high-mass, gap-opening planets we find differences in the amplitude of the resulting dust rings, which we attribute to the failure of the terminal velocity approximation around shocks. For low-mass planets, both models agree everywhere except in the corotation region, where the terminal velocity approximation shows overdense dust lobes. We tentatively interpret these as dusty equivalents of thermal lobes seen in non-isothermal simulations with thermal diffusion, but more work is necessary to confirm this. At the same resolution, the computational time for the terminal velocity approximation model is significantly less than a two fluid model. We conclude that the terminal velocity approximation is a valuable tool for modelling a protoplanetary disc but that care should be taken when shocks are involved.

We present a pipeline for searching for trans-Neptunian objects (TNOs) using data from the TESS mission, that includes a novel optimization-based framework for subtracting the effects of scattered light and pointing jitter. The background subtraction procedure we adopt, when combined with a moving average, allows one to see TNOs such as 90366 SEDNA, 2015 BP519 with the "naked eye." Moreover, this procedure also enabled us to identify two TNO candidates via direct visual observation (subsequently identified to be 2003 UZ413 and 2005 RR43). To automate the extraction of candidate TNOs, we apply a matched filter that can be tuned to objects at different distances and orbital inclinations. We also demonstrate the performance of the algorithm by recovering signals of three trans-Neptunian objects automatically with a high level of confidence. We further validate the approach via synthetic experiments that test recovery rate as a function of magnitude and distance. We find that there is a trade-off between distance and magnitude; controlling for magnitude, it is easier to detect faster moving objects. Our method can detect objects at distances of 250 AU for magnitudes of +21, and closer objects at fainter magnitudes. On a single contemporary GPU (NVIDIA A100) the method can search for 100 trajectories on 1000 2048 x 2048 frames in 5 minutes, dramatically faster than previous approaches. This method can be used to perform large scale fully automated surveys for TNOs and potential Planet Nine candidates.

LiteBIRD is a JAXA-led strategic large-class satellite mission designed to measure the polarization of the cosmic microwave background and Galactic foregrounds from 34 to 448 GHz across the entire sky from L2 in the late 2020s. The scientific payload includes three telescopes which are called the low-, mid-, and high-frequency telescopes each with their own receiver that covers a portion of the mission's frequency range. The low frequency telescope will map synchrotron radiation from the Galactic foreground and the cosmic microwave background. We discuss the design, fabrication, and characterization of the low-frequency focal plane modules for low-frequency telescope, which has a total bandwidth ranging from 34 to 161 GHz. There will be a total of 4 different pixel types with 8 overlapping bands to cover the full frequency range. These modules are housed in a single low-frequency focal plane unit which provides thermal isolation, mechanical support, and radiative baffling for the detectors. The module design implements multi-chroic lenslet-coupled sinuous antenna arrays coupled to transition edge sensor bolometers read out with frequency-domain mulitplexing. While this technology has strong heritage in ground-based cosmic microwave background experiments, the broad frequency coverage, low optical loading conditions, and the high cosmic ray background of the space environment require further development of this technology to be suitable for LiteBIRD. In these proceedings, we discuss the optical and bolometeric characterization of a triplexing prototype pixel with bands centered on 78, 100, and 140 GHz.

BL Lacertae (BL Lac) objects are one type of blazar, distinguished by their featureless optical spectrum. This presents a challenge in measuring the redshift of the BL Lacs. This paper uses the photometric dropout technique to measure the redshifts of BL Lac objects. Space-based telescope \emph{Swift} and ground-based SARA telescopes are employed to provide magnitudes in the $uvw2,\ uvm2,\ uvw1,\ u,\ b,\ v,\ g',\ r',\ i',\ z'$ filters. We observe 60 BL Lacs and report reliable redshift upper limits for 41 sources. We discover four high-$z$ BL Lacs ($z>1.3$), bringing the number of high-$z$ BL Lacs found by this method up to 20. We also discuss the blazar sequence, the \emph{Fermi} blazar divide, and the gamma-ray horizon using the 4LAC catalog and all high-$z$ BL Lacs discovered with the photo-z technique.

Observations of the Lyman-$\alpha$ (Ly$\alpha$) forest from spectroscopic surveys such as BOSS/eBOSS, or the ongoing DESI, offer a unique window to study the growth of structure on megaparsec scales. Interpretation of these measurements is a complicated task, requiring hydrodynamical simulations to model and marginalise over the thermal and ionisation state of the intergalactic medium. This complexity has limited the use of Ly$\alpha$ clustering measurements in joint cosmological analyses. In this work we show that the cosmological information content of the 1D power spectrum ($P_\mathrm{1D}$) of the Ly$\alpha$ forest can be compressed into a simple two-parameter likelihood without any significant loss of constraining power. We simulate $P_\mathrm{1D}$ measurements from DESI using hydrodynamical simulations and show that the compressed likelihood is model independent and lossless, recovering unbiased results even in the presence of massive neutrinos or running of the primordial power spectrum.

We present the first numerical analysis of causal, stable first-order relativistic hydrodynamics with ideal gas microphysics, based in the formalism developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK theory). The BDNK approach provides definitions for the conserved stress-energy tensor and baryon current, and rigorously proves causality, local well-posedness, strong hyperbolicity, and linear stability (about equilibrium) for the equations of motion, subject to a set of coupled nonlinear inequalities involving the undetermined model coefficients (the choice for which defines the "hydrodynamic frame"). We present a class of hydrodynamic frames derived from the relativistic ideal gas "gamma-law" equation of state which satisfy the BDNK constraints, and explore the properties of the resulting model for a series of (0+1)D and (1+1)D tests in 4D Minkowski spacetime. These tests include a comparison of the dissipation mechanisms in Eckart, BDNK, and Muller-Israel-Stewart theories, as well as investigations of the impact of hydrodynamic frame on the causality and stability properties of Bjorken flow, planar shockwave, and heat flow solutions.

The standard Bayesian technique for searching pulsar timing data for gravitational wave (GW) bursts with memory (BWMs) using Markov Chain Monte Carlo (MCMC) sampling is very computationally expensive to perform. In this paper, we explain the implementation of an efficient Bayesian technique for searching for BWMs. This technique makes use of the fact that the signal model for Earth-term BWMs (BWMs passing over the Earth) is fully factorizable. We estimate that this implementation reduces the computational complexity by a factor of 100. We also demonstrate that this technique gives upper limits consistent with published results using the standard Bayesian technique, and may be used to perform all of the same analyses that standard MCMC techniques can perform.

Supernova spectral time series can be used to reconstruct a spatially resolved explosion model known as supernova tomography. In addition to an observed spectral time series, a supernova tomography requires a radiative transfer model to perform the inverse problem with uncertainty quantification for a reconstruction. The smallest parametrizations of supernova tomography models are roughly a dozen parameters with a realistic one requiring more than 100. Realistic radiative transfer models require tens of CPU minutes for a single evaluation making the problem computationally intractable with traditional means requiring millions of MCMC samples for such a problem. A new method for accelerating simulations known as surrogate models or emulators using machine learning techniques offers a solution for such problems and a way to understand progenitors/explosions from spectral time series. There exist emulators for the TARDIS supernova radiative transfer code but they only perform well on simplistic low-dimensional models (roughly a dozen parameters) with a small number of applications for knowledge gain in the supernova field. In this work, we present a new emulator for the radiative transfer code TARDIS that not only outperforms existing emulators but also provides uncertainties in its prediction. It offers the foundation for a future active-learning-based machinery that will be able to emulate very high dimensional spaces of hundreds of parameters crucial for unraveling urgent questions in supernovae and related fields.

We briefly discuss explicit compact object solutions in higher-order scalar-tensor theories. We start by so-called stealth solutions, whose metric are General Relativity (GR) solutions, but accompanied by a non-trivial scalar field, in both spherically-symmetric and rotating cases. The latter then enables to construct an analytic stationary solution of scalar tensor theory which is called disformed Kerr metric. This solution constitutes a measurable departure from the usual Kerr geometry of GR. We finally consider a scalar-tensor theory stemming from a Kaluza-Klein reduction of a higher-dimensional Lovelock theory, and which enables to obtain non-stealth black holes, highly compact neutron stars and finally wormhole solutions.

We consider the impact of a kinetic pole of order one or two on the non-supersymmetric model of hybrid inflation. These poles arise due to logarithmic Kaehler potentials which control the kinetic mixing of the inflaton field and parameterize hyperbolic manifolds with scalar curvature related to the coefficient (-N)<0 of the logarithm. Inflation is associated with the breaking of a local SU(2)xU(1) symmetry, which does not produce any cosmological defects after it, and remains largely immune from the minimal possible radiative corrections to the inflationary potential. For N=1 and equal values of the relevant coupling constants, lambda and kappa, the achievement of the observationally central value of ns requires the mass parameter, m, and the symmetry breaking scale, M, to be of the order of 10^12 GeV and 10^17 GeV respectively. Increasing N above unity the tensor-to-scalar ratio r increases above 0.002 and reaches its maximal allowed value for N~10-20.

We review the composition and the equation of state of the hyperonic core of neutron stars at finite temperature within a relativistic mean-field approach. We make use of the new FSU2H* model, which is built upon the FSU2H scheme by improving on the Xi potential according to the recent analysis on the Xi atoms, and we extend it to include finite temperature corrections. The calculations are done for a wide range of densities, temperatures and charge fractions, thus exploring the different conditions that can be found in protoneutron stars, binary mergers remnants and supernovae explosions. The inclusion of hyperons has a strong effect on the composition and the equation of state at finite temperature, which consequently would lead to significant changes in the properties and evolution of hot neutron stars.

In this work we consider a combined theoretical framework comprised by $F(R)$ gravity and a kinetic scalar field. The kinetic energy of the scalar field dominates over its potential for all cosmic times, and the kinetic scalar potential is chosen to be small and non-trivial. In this case, we show that the primordial gravitational wave energy spectrum of vacuum $F(R)$ gravity is significantly enhanced and can be detectable in future interferometers. The kinetic scalar thus affects significantly the inflationary era, since it extends its duration, but also has an overall amplifying effect on the energy spectrum of pure $F(R)$ gravity primordial gravitational waves. The form of the signal is characteristic for all these theories, since it is basically flat and should be detectable from all future gravitational wave experiments for a wide range of frequencies, unless some unknown damping factor occurs due to some unknown physical process.

Dark matter (DM) can be captured in celestial bodies after scattering and losing sufficient energy to become gravitationally bound. We derive a general framework that describes the current DM distribution inside celestial objects, which self-consistently includes the effects of concentration diffusion, thermal diffusion, gravity, and capture accumulation. For DM with sufficient interactions, we show that a significant DM population can thermalize and sit towards the celestial-body surface. This floating distribution allows for new phenomenology for DM searches in a wide range of celestial bodies, including the Sun, Earth, Jupiter, Brown Dwarfs, and Exoplanets.