Papers for Wednesday, Jan 18 2023

Current observations favour that the massive ultraviolet-bright clumps with a median stellar mass of $\sim 10^7~M_{\odot}$, ubiquitously observed in $z\sim 1-3$ galaxies, are star-forming regions formed in-situ in galaxies. It has been proposed that they result from gas fragmentation due to gravitational instability of gas-rich, turbulent, high-redshift discs. We bring support to this scenario by reporting the new discovery of giant molecular clouds (GMCs) in the strongly lensed, clumpy, main-sequence galaxy, A521-sys1, at $z=1.043$. Its CO(4-3) emission was mapped with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of $0.19''\times 0.16''$, reading down to 30~pc thanks to gravitational lensing. We identified 14 GMCs, most being virialized, with $10^{5.9}- 10^{7.9}~M_{\odot}$ masses and a median $800~M_{\odot}~\rm{pc}^{-2}$ molecular gas mass surface density, that are, respectively, 100 and 10 times higher than for local GMCs. They are also characterized by 10 times higher supersonic turbulence with a median Mach number of 60. They end up to fall above the Larson scaling relations, similarly to the GMCs in another clumpy $z\simeq 1$ galaxy, the Cosmic Snake, although noteworthy differences between the two sets of high-redshift GMCs exist. Altogether they support that GMCs form with properties that adjust to the ambient interstellar medium conditions prevalent in the host galaxy whatever its redshift. The detected A521-sys1 GMCs are massive enough to be the parent gas clouds of stellar clumps, with a relatively high star-formation efficiency per free-fall time of $\sim 11$ per cent.

Globular clusters (GCs) are powerful tracers of the galaxy assembly process, and have already been used to obtain a detailed picture of the progenitors of the Milky Way. Using the E-MOSAICS cosmological simulation of a (34.4 Mpc)$^3$ volume that follows the formation and co-evolution of galaxies and their star cluster populations, we develop a method to link the origin of GCs to their observable properties. We capture this complex link using a supervised deep learning algorithm trained on the simulations, and predict the origin of individual GCs (whether they formed in the main progenitor or were accreted from satellites) based solely on extragalactic observables. An artificial neural network classifier trained on $\sim50,000$ GCs hosted by $\sim 700$ simulated galaxies successfully predicts the origin of GCs in the test set with a mean accuracy of $89$ per cent for the objects with [Fe/H]<-0.5 that have unambiguous classifications. The network relies mostly on the alpha-element abundances, metallicities, projected positions, and projected angular momenta of the clusters to predict their origin. A real-world test using the known progenitor associations of the Milky Way GCs achieves up to $90$ per cent accuracy, and successfully identifies as accreted most of the GCs in the inner Galaxy associated to the Kraken progenitor, as well as all the Gaia-Enceladus GCs. We demonstrate that the model is robust to observational uncertainties, and develop a method to predict the classification accuracy across observed galaxies. The classifier can be optimized for available observables (e.g. to improve the accuracy by including GC ages), making it a valuable tool to reconstruct the assembly histories of galaxies in upcoming wide-field surveys.

We present new 0.3-21 micron photometry of SN 2021aefx in the spiral galaxy NGC 1566 at +357 days after B-band maximum, including the first detection of any SN Ia at >15 micron. These observations follow earlier JWST observations of SN 2021aefx at +255 days after the time of maximum brightness, allowing us to probe the temporal evolution of the emission properties. We measure the fraction of flux emerging at different wavelengths and its temporal evolution. Additionally, the integrated 0.3-14 micron decay rate of $\Delta m_{0.3-14} = 1.35 \pm 0.05$ mag/100 days is higher than the decline rate from the radioactive decay of $^{56}$Co of $\sim 1.2$mag/100 days. The most plausible explanation for this discrepancy is that flux is shifting to >14 micron, and future JWST observations of SNe Ia will be able to directly test this hypothesis. However, models predicting non-radiative energy loss cannot be excluded with the present data.

The formation of planets like Earth is expected to conclude with a series of late-stage giant impacts that generate warm dusty debris, the most anticipated visible signpost of terrestrial planet formation in progress. While there is now evidence that Earth-sized terrestrial planets orbit a significant fraction of solar-type stars, the anticipated dusty debris signature of their formation is rarely detected. Here we discuss several ways in which our current ideas about terrestrial planet formation imply transport mechanisms capable of erasing the anticipated debris signature. A tenuous gas disk may be regenerated via "takeout" (i.e., the liberation of planetary atmospheres in giant impacts) or "delivery" (i.e., by asteroids and comets flung into the terrestrial planet region) at a level sufficient to remove the warm debris. The powerful stellar wind from a young star can also act, its delivered wind momentum producing a drag that removes warm debris. If such processes are efficient, terrestrial planets may assemble inconspicuously, with little publicity and hoopla accompanying their birth. Alternatively, the rarity of warm excesses may imply that terrestrial planets typically form very early, emerging fully formed from the nebular phase without undergoing late-stage giant impacts. In either case, the observable signposts of terrestrial planet formation appear more challenging to detect than previously assumed. We discuss observational tests of these ideas.

We present spatially resolved kinematics of 31 ALMA-identified dust-obscured star-forming galaxies (DSFGs) at $z\sim$1.3-2.6, as traced by H$\alpha$ emission using VLT/KMOS near-infrared integral field spectroscopy from our on-going Large Programme ''KMOS-ALMA Observations of Submillimetre Sources'' (KAOSS). We derive H$\alpha$ rotation curves and velocity dispersion profiles for the DSFGs. Of the 31 sources with bright, spatially extended H$\alpha$ emission, 25 display rotation curves that are well fit by a Freeman disc model, enabling us to measure a median inclination-corrected velocity at 2.2$R_{\rm d}$ of $v_{\rm rot}$ = 190 $\pm$ 30 kms$^{-1}$ and a median intrinsic velocity dispersion of $\sigma_0$ = 87 $\pm$ 6 kms$^{-1}$ for these $\textit{disc-like}$ sources. By comparison with less actively star-forming galaxies, KAOSS DSFGs are both faster rotating and more turbulent, but have similar $v_{\rm rot}/\sigma_0$ ratios, median 2.4 $\pm$ 0.5. We suggest that $v_{\rm rot}/\sigma_0$ alone is insufficient to describe the kinematics of DSFGs, which are not kinematically ''cold'' discs, and that the individual components $v_{\rm rot}$ and $\sigma_0$ indicate that they are in fact turbulent, but rotationally supported systems in $\sim$50 per cent of cases. This turbulence may be driven by star formation or mergers/interactions. We estimate the normalisation of the stellar Tully-Fisher relation (sTFR) for the disc-like DSFGs and compare it with local studies, finding no evolution at fixed slope between $z\sim$2 and $z\sim$0. Finally, we use kinematic estimates of DSFG halo masses to investigate the stellar-to-halo mass relation, finding our sources to be consistent with shock heating and strong feedback which likely drives the declining stellar content in the most massive halos.

When the accreting white dwarf in a magnetic cataclysmic variable star (mCV) has a field strength in excess of 10 MG, it is expected to synchronize its rotational frequency to the binary orbit frequency, particularly at small binary separations, due to the steep radial dependence of the magnetic field. We report the discovery of an mCV (SDSS J134441.83+204408.3; hereafter, J1344) that defies this expectation by displaying asynchronous rotation ($P_{spin}/P_{orb} = 0.893$) in spite of a high surface field strength (B=56 MG) and a short orbital period (114 min). Previously misidentified as a synchronously rotating mCV, J1344 was observed by TESS during sector 50, and the resulting power spectrum shows distinct spin and orbital frequencies, along with various sidebands and harmonics. Although there are several other asynchronous mCVs at short orbital periods, the presence of cyclotron humps in J1344's SDSS spectrum makes it possible to directly measure the field strength in the cyclotron-emitting region; a previously study estimated 65 MG based on its identification of two cyclotron humps, but we revise this to 56$\pm$2 MG based on the detection of a third hump and on our modeling of the cyclotron spectrum. Short-period mCVs with field strengths above 10 MG are normally expected to be synchronous, so the highly asynchronous rotation in J1344 presents an interesting challenge for theoretical studies of spin-period evolution.

We present the mid-IR (MIR) morphologies for 70 star-forming galaxies (SFGs) at $0.2<z<2.5$ with stellar mass $\mathrm{M_*>10^{9}~M_\odot}$ using JWST MIRI observations from the Cosmic Evolution Early Release Science survey (CEERS). The MIRI bands span the MIR (7.7--21 $\mu$m), enabling us to measure the effective radii ($R_{\mathrm{eff}}$) and S\'{e}rsic indexes of these SFGs at rest-frame 6.2 and 7.7 $\mu$m, which contains strong emission from Polycyclic aromatic hydrocarbon (PAH) features, a well-established tracer of star formation in galaxies. We compare the PAH-band morphologies to those in rest-frame Near-UV (NUV) using HST ACS/F435W or ACS/F606W and optical/near-IR using HST WFC3/F160W imaging from UVCANDELS and CANDELS, where the NUV-band and F160W trace the profile of (unobscured) massive stars and the stellar continuum, respectively. The $R_{\mathrm{eff}}$ of galaxies in the PAH-band are slightly smaller ($\sim$10) than those in F160W for more massive galaxies with $\mathrm{M_*\gtrsim10^{9.5}~M_\odot}$ at $z\leq1.2$ and with $\mathrm{M_*\gtrsim10^{10}~M_\odot}$ at $z\geq1.2$, but the PAH-band and F160W have a similar fractions of light within 1 kpc. In contrast, the $R_{\rm{eff}}$ of galaxies in the NUV-band are larger, with lower fractions of light within 1 kpc compared to the F160W for galaxies at $z\leq1.2$. Using the MIRI data to estimate the surface density of the IR-based SFR, we find the correlation between the IR-based SFR surface density and stellar mass has a steeper slope than that of the UV-based surface density and stellar mass, suggesting galaxies with higher stellar masses having increasing amounts of obscured fraction of star formation in their inner regions. This paper demonstrates how the high-angular resolution data from JWST/MIRI can reveal new information about the morphology of obscured-star formation.

In the JWST, Extremely Large Telescopes, and LUVOIR era, we expect to characterize a number of potentially habitable Earth-like exoplanets. However, the characterization of these worlds depends crucially on the accuracy of theoretical models. Validating these models against observations of planets with known properties will be key for the future characterization of terrestrial exoplanets. Due to its sensitivity to the micro- and macro-physical properties of an atmosphere, polarimetry will be an important tool that, in tandem with traditional flux-only observations, will enhance the capabilities of characterizing Earth-like planets. In this paper we benchmark two different polarization-enabled radiative-transfer codes against each other and against unique linear spectropolarimetric observations of the earthshine that cover wavelengths from $\sim$0.4 to $\sim$2.3 $\mu$m. We find that while the results from the two codes generally agree with each other, there is a phase dependency between the compared models. Additionally, with our current assumptions, the models from both codes underestimate the level of polarization of the earthshine. We also report an interesting discrepancy between our models and the observed 1.27 $\mu$m $O_2$ feature in the earthshine, and provide an analysis of potential methods for matching this feature. Our results suggest that only having access to the 1.27 $\mu$m $O_2$ feature coupled with a lack of observations of the $O_2$ A and B bands could result in a mischaracterization of an Earth-like atmosphere. Providing these assessments is vital to aid the community in the search for life beyond the solar system.

Due to its proximity, the stellar populations of the Galactic bulge (GB) can be resolved and can be studied in detail. This allows tracing the bulge metallicity distribution function (MDF) for different spatial regions within the bulge, which may give us clues about the bulge formation and evolution scenarios. In this work, we developed a chemical evolution model (CEM), taking into account the mass distribution in the bulge and disc, to derive the radial dependence of this time-scale in the Galaxy. Since the infall rate depends on that time scale in the CEM, the results of the model were used to test a scenario where the bulge is formed inside-out. The obtained results for the $[\alpha/\mbox{Fe}]$ vs. [Fe/H] relationship, the MDF and the [Fe/H] radial gradient in the bulge have been compared to available data in the literature. The model is able to reproduce most of the observational data: the spread in the relation $[\alpha/\mbox{Fe}]$ vs. [Fe/H], the MDF shape in different regions of the bulge, the [Fe/H] radial gradient inside it and the age-metallicity relation, as well as the [$\alpha$/Fe] evolution with age. The results of the model point to a scenario where the bulk of the bulge stars pre-existed the boxy/peanut X-shape bar formation. As a result, the classical origin of the GB is not ruled out and this scenario may be invoked to explain the chemical properties of the Galactic bulge.

We have re-evaluated recent studies of effects on Earth by cosmic rays (CRs) from nearby supernovae (SNe) at 100 pc and 50 pc, in the diffusive transport CR case, here including an early time suppression at lower CR energies neglected in the previous work. Inclusion of this suppression leads to lower overall CR flux at early times, lower atmospheric ionization, smaller resulting ozone depletion, and lower sea-level muon radiation dose. Differences in atmospheric impacts are most pronounced for the 100 pc case with less significant differences in the 50 pc case. We find a greater discrepancy in the modeled sea-level muon radiation dose, with significantly smaller dose values in the 50 pc case; our results indicate it is unlikely that muon radiation is a significant threat to the biosphere for SNe beyond at least 20 pc. We have also performed new modeling of effects by SNe at 20 pc and 10 pc. Overall, our results indicate that, considering only effects of SN CRs, the "lethal" SN distance may be closer to 20 pc rather than the typically quoted 8-10 pc.

A tremendous amount of radiation is emitted by the Interstellar Medium in the mid- and far-infrared (3-500 {\mu}m) that represents the majority of the light emitted by a galaxy. In this article we motivate ISM studies in the infrared and the construction of large specialized observatories like the Stratospheric Observatory For Infrared Astronomy (SOFIA), which just concluded its mission on a scientific high note, and the newly launched James Webb Space Telescope (JWST) that just begun its exciting scientific mission. We introduce their capabilities, present a few examples of their scientific discoveries and discuss how they complemented each other. We then consider the impact of the conclusion of SOFIA for the field in a historic context and look at new opportunities specifically for far-infrared observatories in space and in the stratosphere.

We investigate the generation and evolution of switchbacks (SBs), the nature of the sub-Alfv\'enic wind observed by Parker Solar Probe (PSP), and the morphology of the Alfv\'enic transition, all of which are key issues in solar wind research. First we highlight a special structure in the pristine solar wind, termed a low Mach-number boundary layer (LMBL). An increased Alfv\'en radius and suppressed SBs are observed within an LMBL. A probable source on the Sun for an LMBL is the peripheral region inside a coronal hole with rapidly diverging open fields. The sub-Alfv\'enic wind detected by PSP is an LMBL flow by nature. The similar origin and similar properties of the sub-Alfv\'enic intervals favor a wrinkled surface for the morphology of the Alfv\'enic transition. We find that a larger deflection angle tends to be associated with a higher Alfv\'en Mach number. The magnetic deflections have an origin well below the Alfv\'en critical point, and deflection angles larger than $90^{\circ}$ seem to occur only when $M_{\rm A} \gtrsim 2$. The velocity enhancement in units of the local Alfv\'en speed generally increases with the deflection angle, which is explained by a simple model. A nonlinearly evolved, saturated state is revealed for SBs, where the local Alfv\'en speed is roughly an upper bound for the velocity enhancement. In the context of these results, the most promising theory on the origin of SBs is the model of expanding waves and turbulence, and the patchy distribution of SBs is attributed to modulation by reductions in the Alfv\'en Mach number. Finally, a picture on the generation and evolution of SBs is created based on the results.

Black hole X-ray binaries are ideal environments to test the accretion phenomena in the presence of strong gravitational potentials. KERRBB held an important place in the X-ray spectral continuum method for measuring the black hole spin modeling the emission from the innermost regions of the accretion disk. In this work, we present the results of X-ray spectral analysis using publicly available RXTE data of GRO J1655-40 obtained during the 2005 outburst with the two relativistic accretion disk models, KERRBB and KYNBB. Our analysis showed that both models provide identical results with black hole spin measurements, disk temperature, and disk luminosity when the inner edge of the accretion disk is set at the innermost stable circular orbit (ISCO) for the same accretion rates. We couldn't obtain reasonable fits for $\sim$ 89\% of the observations with a fixed black hole spin value at $\mathrm{a_{*}=0.7}$ using both models. Allowing the spin parameter to vary improved the fit statistic significantly with reduced $\rm \chi^{2}$ values being reduced from $\sim$ 10-100 to below 2. Both models revealed black hole spin values varying between $\rm 0.52<a_{*}<0.94$, which can be interpreted as a variable inner edge of the disk throughout different accretion states.

During the British Astronomical Association (BAA) 2022 campaign, 27436 photometric observations of the dwarf nova (DN) CG Draconis were made, with 106 eclipses recorded. This work summarizes the new data available and provides updated ephemeris and commentary on the observed eclipse profiles. The orbital period found is P_orb = 4h31m38s +/- 1s. Two types of quasi-periodic outbursts are identified: normal outbursts, of Delta V of approximately 1.25 mag amplitude, and bright, of Delta V of approximately 1.5 mag. The pattern resembles superoutbursts of SU UMa-type DNe, however, no presence of superhumps characterizing these DNe was found. Given CG Dra is located above the period gap, it may represent a new intermediary subtype between SS Cyg and SU UMa-type stars, or provide support to superoutburst models that do not rely on eccentric accretion disks.

The physical properties of star cluster populations offer valuable insights into their birth, evolution, and disruption. However, individual stars in clusters beyond the nearest neighbours of the Milky Way are unresolved, forcing analyses of star cluster demographics to rely on integrated light, a process that is fraught with uncertainty. Here infer the demographics of the cluster population in the benchmark galaxy NGC 628 using data from the Legacy Extra-galactic UV Survey (LEGUS) coupled to a novel Bayesian forward-modelling technique. Our method allows analysis of a total of 1178 clusters in the LEGUS catalogue, roughly a factor of 4 more than previous studies that required severe completeness cuts to the data. Our results indicate that the cluster mass function is truncated at $\sim 10^4$ $ \mathrm{M}_{\odot}$, consistent with proposed relations between truncation mass and star formation surface density. We find that cluster disruption begins early, $\sim 10$ Myr after formation, but that it is relatively mild, with clusters requiring on average $2-3$ times their present age to disrupt; we do not find any evidence for mass dependent disruption. We also do not find convincing evidence for any radial variations in these conclusions, though we find suggestive hints that inner galaxy clusters may be more prone to disruption. Confirming or refuting these hints will require future observations to increase the sample size of outer galaxy clusters.

High-resolution near infrared observations with GRAVITY instrument have revealed rapid orbital motions of a hot spot around Sgr A*, the supermassive black hole in our Galactic center, during its three bright flares. The projected distances of the spot to the black hole are measured and seems to increase with time. The values of distance, combined with the measured orbiting time, imply that the spot is rotating with a super-Keplerian velocity. These results are hard to understand if the spot stay within the accretion flow thus provide strong constraints on theoretical models for flares. Previously we have proposed a "ME" model for the flares by analogy with the coronal mass ejection model in solar physics. In that model, magnetic reconnection occurred at the surface of the accretion flow results in the formation of flux ropes, which are then ejected out. Energetic electrons accelerated in the current sheet flow into the flux rope region and their radiation is responsible for the flares. In this paper, we apply the model to the interpretation of the GRAVITY results by calculating the dynamics of the ejected flux rope, the evolution of the magnetic field and the energy distribution of accelerated electrons, and the radiation of the system. We find that the model can well explain the observed light curve of the flares, the time-dependent distance and the super-Keplerian motion of the hot spot. It also explains why the light curve of some flares have double peaks.

The role and implication of binding energy through the accretion-induced collapse (AIC) of accreting white dwarfs (WDs) for the production of millisecond pulsars (MSPs) are investigated. I examine the binding energy model due to the dynamical process in close binary systems and investigating the possible mass of the companion sufficient to induce their orbital parameters. The deterministic nature of this interaction has a strong sensitivity to the equation of state of the binary systems (where the compactness of a neutron star is proportional to the amount of binding energy) associated with their initial conditions. This behavior will mimic the commonly assumed mass and amount of accreted matter under the instantaneous mass loss ($\Delta M \sim 0.18M_{\odot}$). As a result, this will indicate an increase in the MSP's gravitational mass due to angular momentum losses. The outcome of such a system will then be a circular binary MSP in which the companion is a low-mass WD, thus distinguishing the binary formation scenarios. In addition, the results of this work could provide constraints on the expected mass and binding energy of a neutron star based on the accretion rate

We derive new empirical scaling relations between WISE mid-infrared galaxy photometry and well-determined stellar masses from SED modeling of a suite of optical-infrared photometry provided by the DR4 Catalogue of the GAMA-KiDS-VIKING survey of the southern G23 field. The mid-infrared source extraction and characterization are drawn from the WISE Extended Source Catalogue (WXSC) and the archival ALLWISE catalog, combining both resolved and compact galaxies in the G23 sample to a redshift of 0.15. Three scaling relations are derived: W1 3.4 micron luminosity versus stellar mass, and WISE W1-W2, W1-W3 colors versus mass-to-light ratio (sensitive to a variety of galaxy types from passive to star-forming). For each galaxy in the sample, we then derive the combined stellar mass from these scaling relations, producing Mstellar estimates with better than $\sim$25-30% accuracy for galaxies with $>$10$^{9}$ Msolar and $<$40 - 50% for lower luminosity dwarf galaxies. We also provide simple prescriptions for rest-frame corrections and estimating stellar masses using only the W1 flux and the W1-W2 color, making stellar masses more accessible to users of the WISE data. Given a redshift or distance, these new scaling relations will enable stellar mass estimates for any galaxy in the sky detected by WISE with high fidelity across a range of mass-to-light.

We report Very Large Array (VLA) observations in the Q-band toward 10 ionized jet candidates to search for SiO emission, a well-known shocked gas tracer. We detected 7 mm continuum counterparts toward 90% of the jet candidates. In most cases, the jet candidate is located toward the center of the 7 mm core, and the high masses ($\approx 100\,M_\odot$) and densities ($\approx 10^7\, \text{cm}^{-3}$) of the cores suggest that the central objects are very young high-mass protostars. We detected SiO $J=1-0$ emission associated with 6 target sources. In all cases, the morphology and spectrum of the emission is consistent with what is expected for molecular jets along an outflow axis, thus confirming the jet nature of 60% of our sample. Our data suggest a positive correlation between the SiO luminosity $L_{SiO}$, and both the bolometric luminosity $L_{Bol}$ and the radio luminosity $S_\nu d^2$ of the driving sources.

We present a comprehensive set of physical and geometrical parameters for each of the components of the close visual binary system Hip 11253 (HD14874). We present an analysis for the binary and multiple stellar systems with the aim to obtain a match between the overall observational spectral energy distribution of the system and the spectral synthesis created from model atmospheres. The epoch positions are used to determine the orbital parameters and the total mass.

We report the first $5\sigma$ detection of HI 21 cm emission from a star-forming galaxy at redshift $z\sim1.3$ (nearly 9 billion years ago) using upgraded Giant Metrewave Radio Telescope (uGMRT). This is the highest redshift HI detection in emission from an individual galaxy to date. The emission is strongly boosted by the gravitational lens, an early type elliptical galaxy, at redshift $z \sim 0.13$. The measured HI mass of the galaxy is $\rm M_{HI} = (0.90 \pm 0.14 \pm 0.05) \times 10^{10}M_{\odot}$, which is almost twice the inferred stellar mass of the galaxy, indicating an extended structure of the HI gas inside the galaxy. By fitting two-dimensional Gaussian to the HI signal at the peak of the spectral line, we find the source to be marginally resolved with the position angle consistent with the emission being tangential to the critical curve of the lens mass distribution. This indicates that the solid angle of the approaching HI line flux comes very close to the inner lens caustic and results in very high magnification. These results, for the first time, demonstrate the feasibility of observing high redshift HI in a lensed system with a modest amount of telescope time and open up exciting new possibilities for probing the cosmic evolution of neutral gas with existing and upcoming low-frequency radio telescopes in the near future.

Jet feedback from active galactic nuclei (AGN) is one of the most promising mechanisms for suppressing cooling flows in cool-core clusters. However, the composition of AGN jets and bubbles remains uncertain; they could be thermally dominated, or dominated by cosmic-ray proton (CRp), cosmic-ray electron (CRe), or magnetic energy. In this work, we investigate the evolution and feedback effects of CRp and CRe dominated jets by conducting 3D magnetohydrodynamic simulations of AGN jet-inflated bubbles in the intracluster medium using the FLASH code. We present the evolution of their energies, dynamics and heating, and model their expected cavity-power versus radio-luminosity relation ($P_{\rm cav}-L_R$). We find that bubbles inflated by CRe dominated jets follow a very similar dynamical evolution to CRp dominated bubbles even though CRe within bubbles suffer significantly stronger synchrotron and inverse-Compton cooling. This is because, as CRe lose their energy, the jet-inflated bubbles quickly become thermally dominated within $\sim 30$ Myr. Their total energy stops decreasing with CR energy and evolves similarly to CRp dominated bubbles. The ability of CRe and CRp dominated bubbles to heat the intracluster medium is also comparable; the cold gas formed via local thermal instabilities is well suppressed in both cases. The CRp and CRe bubbles follow different evolutionary trajectories on the $P_{\rm cav}-L_R$ plane, but the values are broadly consistent with observed ranges for FRI sources. We also discuss observational techniques that have potential for constraining the composition of AGN jets and bubbles.

We map the surface of Pluto using an unsupervised machine learning technique using the near-infrared observations of the LEISA/Ralph instrument onboard NASA's New Horizons spacecraft. The principal component reduced Gaussian mixture model was implemented to investigate the geographic distribution of the surface units across the dwarf planet. We also present the likelihood of each surface unit at the image pixel level. Average I/F spectra of each unit were analyzed -- in terms of the position and strengths of absorption bands of abundant volatiles such as N${}_{2}$, CH${}_{4}$, and CO and nonvolatile H${}_{2}$O -- to connect the unit to surface composition, geology, and geographic location. The distribution of surface units shows a latitudinal pattern with distinct surface compositions of volatiles -- consistent with the existing literature. However, previous mapping efforts were based primarily on compositional analysis using spectral indices (indicators) or implementation of complex radiative transfer models, which need (prior) expert knowledge, label data, or optical constants of representative endmembers. We prove that an application of unsupervised learning in this instance renders a satisfactory result in mapping the spatial distribution of ice compositions without any prior information or label data. Thus, such an application is specifically advantageous for a planetary surface mapping when label data are poorly constrained or completely unknown, because an understanding of surface material distribution is vital for volatile transport modeling at the planetary scale. We emphasize that the unsupervised learning used in this study has wide applicability and can be expanded to other planetary bodies of the Solar System for mapping surface material distribution.

Small-scale fluctuating magnetic fields of order $n$G to $\mu$G are observed in supernova shocks and galaxy clusters, where amplifications of the field are likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel three-dimensional magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use $384 \times 512 \times 384$ tracer trajectories to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman power-law kinetic energy decay ($E_{\textrm{kin}} \propto t^{-1.15})$ in the absence of turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution, because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short, with only $\sim0.3$ of a turbulent turnover time. Yet, the growth rates are still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings via the Eulerian-Lagrangian bridge.

Current and upcoming radio-interferometers are expected to produce volumes of data of increasing size that need to be processed in order to generate the corresponding sky brightness distributions through imaging. This represents an outstanding computational challenge, especially when large fields of view and/or high resolution observations are processed. We have investigated the adoption of modern High Performance Computing systems specifically addressing the gridding, FFT-transform and w-correction of imaging, combining parallel and accelerated solutions. We have demonstrated that the code we have developed can support dataset and images of any size compatible with the available hardware, efficiently scaling up to thousands of cores or hundreds of GPUs, keeping the time to solution below one hour even when images of the size of the order of billion or tens of billion of pixels are generated. In addition, portability has been targeted as a primary objective, both in terms of usability on different computing platforms and in terms of performance. The presented results have been obtained on two different state-of-the-art High Performance Computing architectures.

We present two-dimensional hydrodynamic simulations of the accretion-induced collapse (AIC) of rotating white dwarfs admixed with an extended component of dark matter (DM) comprising of sub-GeV degenerate fermionic DM particles. We find that the DM component would follow the collapse of the normal matter (NM) component to become a bound DM core. Thus, we demonstrate how a DM-admixed neutron star could form through DM-admixed AIC (DMAIC) for the first time, with the dynamics of DM taken into account. The gravitational-wave (GW) signature from the DMAIC shows distinctive features. In the diffusive DM limit, the DM admixture indirectly suppresses the post-bounce spectral peak of the NM GWs. In the compact DM limit, the collapse dynamics of the DM in a Milky Way event generate GWs that are strong enough to be detectable by Advanced LIGO as continuous low-frequency ($< 1000$ Hz) signals after the NM core bounce. Our study not only is the first-ever computation of GW from a collapsing DM object but also provides the key features to identify DM in AIC events through future GW detections.

We explore observational constraints on a cosmological model with an interaction between dark energy (DE) and dark matter (DM), using a compilation of 15 measurements of the 2D BAO (i.e., transversal) scale in combination with Planck-CMB data, to explore the parametric space of a class of interacting DE models. We find that 2D BAO measurements can generate different observational constraints compared to the traditional approach of studying the matter clustering in the 3D BAO measurements. The 2D BAO sample provides strong evidence in favor of the IDE model at more than 3$\sigma$. Also, contrary to the observations for the $\Lambda$CDM and IDE models when analyzed with Planck-CMB + 3D BAO data, we note that Planck-CMB + 2D BAO data favor high values of the Hubble constant $H_0$. From the joint analysis with Planck-CMB + 2D BAO + Gaussian prior on $H_0$, we find $H_0 = 73.4 \pm 0.88$ km/s/Mpc. Our results show that Planck-CMB + 2D BAO measurements form a minimal data set that solves the $H_0$ tension, and at the same time, it provides statistical evidence for the IDE cosmologies.

The strong lensing gravitational wave (SLGW) is a promising transient phenomenon containing rich physics. However, the poor sky localization due to the long-wave nature of gravitational waves makes the identification of such events very challenging. We propose a new method based on the wave optics effect of the microlensing field embedded in SLGW data. The microlensing diffraction/interference fringes can produce frequency-dependent random fluctuations in the waveform. To pin down the microlensing induced stochastic features in the waveform, we utilize both the template-independent method, \texttt{cWB}, and the template-dependent method, \texttt{Bilby}, to reconstruct the waveform with and without microlensing imprints. The mismatching degree of these two waveforms can be treated as an indicator of SLGW events. We forecast the identification rate of this method with the third-generation gravitational wave observatory, such as Cosmic Explorer. Our result shows that this method can successfully identify about 2 (out of 180) SLGW events with strong enough microlensing effect per year. This method is entirely data-driven, which is immune to model priors, and can greatly avoid the false positive errors contaminated by the coincident unlensed events.

We report the results of long-term reverberation mapping (RM) campaigns of the nearby active galactic nuclei (AGN) NGC 4151, spanning from 1994 to 2022, based on archived observations of the FAST Spectrograph Publicly Archived Programs and our new observations with the 2.3m telescope at the Wyoming Infrared Observatory. We reduce and calibrate all the spectra in a consistent way, and derive light curves of the broad H$\beta$ line and 5100\,{\AA} continuum. Continuum light curves are also constructed using public archival photometric data to increase sampling cadences. We subtract the host galaxy contamination using {\it HST} imaging to correct fluxes of the calibrated light curves. Utilizing the long-term archival photometric data, we complete the absolute flux-calibration of the AGN continuum. We find that the H$\beta$ time delays are correlated with the 5100\,{\AA} luminosities as $\tau_{\rm H\beta}\propto L_{5100}^{0.46\pm0.16}$. This is remarkably consistent with Bentz et al. (2013)'s global size-luminosity relationship of AGNs. Moreover, the data sets for five of the seasons allow us to obtain the velocity-resolved delays of the H$\beta$ line, showing diverse structures (outflows, inflows and disks). Combining our results with previous independent measurements, we find the measured dynamics of the H$\beta$ broad-line region (BLR) are possibly related to the long-term trend of the luminosity. There is also a possible additional $\sim$1.86 years time lag between the variation in BLR radius and luminosity. These results suggest that dynamical changes in the BLR may be driven by the effects of radiation pressure.

The K2 eclipsing binary candidates EPIC 211982753 (hereinafter called EPIC2753) and EPIC 211915147 (hereinafter called EPIC5147) are characterized with the help of photometric and high-resolution spectroscopic data. The light curve analysis uses the R-band photometric data from the 1.3-m Devasthal Fast Optical Telescope (DFOT, India), ASAS-3 and K2 observations. High-resolution echelle spectra are collected using the HERMES spectrograph at the 1.2-m MERCATOR telescope (La Palma, Spain). The synthetic light and radial velocity curves are generated with the help of the modeling package PHOEBE 1.0. The orbital period analysis based on the ~3.2 years of K2 observations does not show any change in the orbital period of both targets. The component masses M1,2 are estimated as 1.69(0.02) and 1.59(0.02) solar mass for EPIC2753, and 1.48(0.01) and 1.27(0.01) solar mass for EPIC5147. Both systems are high mass-ratio eclipsing binaries with q>0.85. The component radii R1,2 are found to be 1.66(0.02) and 1.53(0.02) solar radius for EPIC2753, and 1.80(0.05) and 1.42(0.05) solar radius for EPIC5147. The distances of EPIC2753 and EPIC5147 are determined as 238(4) and 199(5) pc, respectively. MESA Isochrones and Stellar Tracks are used to understand the evolutionary status of both systems.

The Hubble constant ($H_0$), which represents the expansion rate of the Universe, is one of the most important cosmological parameters. The recent measurements of $H_0$ using the distance ladder methods such as Type Ia Supernovae (SNe Ia) are significantly greater than the CMB measurements by Planck. The difference points to a crisis in the standard model of cosmology termed as Hubble tension. In this work we compare different cosmological models, determine the Hubble constant and comment on the Hubble tension using the data from differential ages of galaxies. The data we use is free from the systematic effects as the absolute age estimation of the galaxies is not needed. We have used the Bayesian approach along with the commonly used maximum likelihood method to estimate $H_0$ and have calculated the AIC scores to compare the different cosmological models.The non-flat cosmological model provides a higher value for matter density as well as the Hubble constant compared to the flat $\Lambda$CDM model. The AIC score is smaller for the flat $\Lambda$CDM cosmology compared to the non-flat model indicating the flat model a better choice. The best-fit value of $H_0$ for both these models are $68.7\pm3.1$ km/s/Mpc and $72.2\pm4$ km/s/Mpc, respectively. Our results are consistent with the CCHP measurements. However, flat model result does not agree with the SH0ES result, while the non-flat result is inconsistent with the Planck value.

We have used the axisymmetric chemo-hydrodynamical code WALKIMYA-2D to numerically model and reproduce the physical and CO emission properties of the jet-driven outflow from the intermediate-mass protostar Cep E, which was observed at $\sim 800$au resolution in the CO $J=2\to 1$ line with the IRAM interferometer. Our simulations take into account the observational constraints available on the physical structure of the protostellar envelope to provide constraints on the dynamics of the inner protostellar environment from the study of the outflow/jet propagation away from the launch region. WALKIMYA-2D successfully reproduces the main qualitative and quantitative features of the Cep E outflow and the jet kinematics, naturally accounting for their time variability. Signatures of internal shocks are detected as knots along the jet. In the early times of the ejection process, the young emitted knots interact with the dense circumstellar envelope through high-velocity, dissociative shocks, which strongly decrease the CO gas abundance in the jet. As time proceeds, the knots propagate more smoothly through the envelope and dissociative shocks disappear after $\sim 10^3$ yr. The distribution of CO abundance along the jet shows that the latter bears memory of the early dissociative phase in the course of its propagation. Analysis of the velocity field shows that the jet material mainly consists of gas entrained from the circumstellar envelope and accelerated away from the protostar at 700 au scale. As a result, the overall jet mass loss rate appears higher than the actual mass ejection rate by a factor $\sim 3$. Numerical modelling of the Cep E jet-driven outflow and comparison with the CO observations have allowed us to peer into the outflow formation mechanism with unprecedented detail and to retrieve the history of the mass-loss events that have shaped the outflow.

We examine the long-term stability (on decade-like time scales) of optical `high polarization' (HP) state with $p_{opt}$ $> 3\%$, which commonly occurs in flat-spectrum (i.e., beamed) radio quasars (FSRQs) and is a prominent marker of blazar state. Using this clue, roughly a quarter of the FSRQ population has been reported to undergo HP $\leftrightarrow$ non-HP state transition on year-like time scales. This work examines the extent to which HP (i.e., blazar) state can endure in a FSRQ, despite these `frequent' state transitions. This is the first attempt to verify, using purely opto-polarimetric data for a much enlarged sample of blazars, the recent curious finding that blazar state in individual quasars persists for {\it at least} a few decades, despite its changing/swinging observed fairly commonly on year-like time scales. The present analysis is based on a well-defined sample of 83 radio quasars, extracted from the opto-polarimetric survey RoboPol (2013-17), for which old opto-polarimetric data taken prior to 1990 could be found in the literature. By a source-wise comparison of these two datasets of the same observable ($p_{opt}$), we find that $\sim$ 90% of the 63 quasars found in blazar state in our RoboPol sample, were also observed to be in that state about 3 decades before. On the other hand, within the RoboPol survey itself, we find that roughly a quarter of the blazars in our sample migrated to the other polarization state on year-like time scales, by crossing the customary $p_{opt}$ = 3% threshold. Evidently, these relatively frequent transitions (in either direction) do not curtail the propensity of a radio quasar to retain its blazar (i.e., HP) state for at least a few decades. The observed transitions/swings of polarization state are probably manifestation of transient processes, like ejections of synchrotron plasma blobs (VLBI radio knots) from the active nucleus.

We investigate the role of galaxy mergers in triggering AGN in the nearby Universe. Our analysis is based on a sample of 79 post-merger remnant galaxies with deep X-ray observations from Chandra/XMM-Newton capable of detecting a low-luminosity AGN of > 10^40.5 erg s^-1. This sample is derived from a visually classified, volume-limited sample of 807 post-mergers identified in the Sloan Digital Sky Survey Data Release 14 with log M*/M_sun > 10.5 and 0.02 < z < 0.06. We find the X-ray AGN fraction in this sample is 55.7% +\- 5.6% compared to 23.6% +\- 2.8% for a mass and redshift matched non-interacting control sample. The multi-wavelength AGN fraction (identified as an AGN in one of X-ray, IR, radio or optical diagnostics) for post-mergers is 76.6% +\- 4.8% compared to 39.1% +\- 3.2% for controls. Thus post-mergers exhibit a high overall AGN fraction with an excess between 2 - 4 depending on the AGN diagnostics used. In addition, we find most optical, IR, and radio AGN are also identified as X-ray AGN while a large fraction of X-ray AGN are not identified in any other diagnostic. This highlights the importance of deep X-ray imaging to identify AGN. We find the X-ray AGN fraction of post-mergers is independent of stellar mass above log M*/M_sun > 10.5 unlike the trend seen in control galaxies. Overall, our results show that post-merger galaxies are a good tracer of the merger-AGN connection and strongly support the theoretical expectations that mergers trigger AGN.

We present stellar parameters and abundances of 13 elements for 18 very metal-poor (VMP; [Fe/H] $<$ -2.0) stars, selected as extremely metal-poor (EMP; [Fe/H] $<$ -3.0) candidates from SDSS and LAMOST survey. High-resolution spectroscopic observations were performed using GEMINI-N/GRACES. We find ten EMP stars among our candidates, and we newly identify three carbon-enhanced metal poor (CEMP) stars with [Ba/Fe] $<$ 0. Although chemical abundances of our VMP/EMP stars generally follow the overall trend of other Galactic halo stars, there are a few exceptions. One Na-rich star ([Na/Fe] = +1.14) with low [Mg/Fe] suggests a possible chemical connection with second-generation stars in a globular cluster. The progenitor of an extremely Na-poor star ([Na/Fe] = -1.02) with an enhancement of K- and Ni-abundance ratios may have undergone a distinct nucleosynthesis episode, associated with core-collapse supernovae (CCSNe) having a high explosion energy. We have also found a Mg-rich star ([Mg/Fe] = +0.73) with slightly enhanced Na and extremely low [Ba/Fe], indicating that its origin is not associated with neutron-capture events. On the other hand, the origin of the lowest Mg abundance ([Mg/Fe] = -0.61) star could be explained by accretion from a dwarf galaxy, or formation in a gas cloud largely polluted by SNe Ia. We have also explored the progenitor masses of our EMP stars by comparing their chemical-abundance patterns with those predicted by Population III SNe models, and find a mass range of 10 - 26 $M_\odot$, suggesting that such stars were primarily responsible for the chemical enrichment of the early Milky Way.

In the Galaxy, close binaries with compact objects are important low-frequency gravitational wave (GW) sources. As potential low-frequency GW sources, neutron star/white dwarf (WD) ultra-compact X-ray binaries (UCXBs) have been investigated extensively. Using the MESA code, we systematically explored the evolution of black hole (BH)-main sequence star (MS) binaries to diagnose whether their descendants can be detected by space-borne GW detectors. Our simulations show that BH-MS binaries with an initial orbital period less than the bifurcation period can evolve into BH UCXBs that can be detected by LISA. Such an evolutionary channel would form compact mass-transferring BH-WD systems rather than detached BH-WD systems. The calculated X-ray luminosities of BH UCXBs detected by LISA at a distance $d=1$ kpc are $\sim10^{33}-10^{35}~\rm erg\,s^{-1}$ ($\sim10^{34}-10^{35}~\rm erg\,s^{-1}$ for $d=10$ kpc), hence it is possible to detect their electromagnetic counterparts. It is worth emphasizing only some BH-MS systems with an initial orbital period very close to the bifurcation period can evolve toward low-frequency GW sources whose chirp masses can be measured. The maximum GW frequency of BH UCXBs forming by BH-MS pathway is about 3 mHz, which is smaller than the minimum GW frequency (6.4 mHz) of mass-transferring BH-WD originating from a dynamic process. Furthermore, we obtain an initial parameter space (donor-star masses and orbital periods) of progenitors of BH UCXB-GW sources, which can be applied to future population synthesis simulations. By a rough estimation, we predict that LISA could detect only a few BH UCXB-GW sources forming by the BH-MS channel.

In a previous paper, we demonstrated a single-rung method for measuring cosmological distances in active galactic nuclei (AGN) that can be used from low redshift (z < 0.1) to high redshift (z > 3). This method relies on the assumption that the variability seen in AGN is constrained by the speed of light during a flare event and can therefore be used to estimate the size of an emitting region. A limitation of this method is that previously, the Doppler factor was required to be known. In this paper, we derive an extension of the `standard speed-gun' method for measuring cosmological distances that depends on the maximum intrinsic brightness temperature that a source can reach, rather than the Doppler factor. If the precise value of the intrinsic brightness temperature does not evolve with redshift and flares are statistically independent, we can in principle improve the errors in measurements of the matter content of the universe (in a flat LambdaCDM model) statistically. We then explored how well a future observing program would constrain cosmological parameters. We found that recovering the input cosmology depends critically on the uncertainty of the intrinsic brightness temperature and the number of flares observed.

We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader. Upgrading an ordinary non-reversible switching integrator to a reversible one is straightforward and introduces no appreciable computational burden in our tests. Our method checks if the integrator during the time step violates a time-symmetric selection condition and redoes the step if necessary. In our experiments a few percent of steps would have violated the condition without our corrections. By eliminating them the algorithm avoids long-term error accumulation, of several orders magnitude in some cases.

We describe a simulation for the distribution of galaxies focusing on the atomic Hydrogen content. We aim to make predictions for surveys of galaxies using the redshifted 21 cm line emission. We take the expected distribution of HI masses, circular velocities, sizes of galaxies and orientations into account for this simulation. We use the sensitivity of ASKAP and MeeKAT radio telescopes to estimate the number of detections of HI galaxies in upcoming surveys. We validate our simulation with earlier estimates carried out by using some of these considerations. We show that unlike earlier simulations that take some of the factors into account, the predicted number of galaxies and their distribution across masses changes significantly when all of these are accounted for. We describe our predictions for the MIGHTEE-HI and WALLABY surveys for blind detection of galaxies using the redshifted 21 cm radiation. We study the dependence of the predicted number of detections on the HI mass function. We also describe our future plans for improving the simulation.

We report the first constraints on the growth rate of the universe, $f(z)\sigma_8(z)$, with intrinsic alignments (IA) of galaxies. We measure the galaxy density-intrinsic ellipticity cross-correlation and intrinsic ellipticity auto-correlation functions over $0.16 < z < 0.7$ from luminous red galaxies (LRG) and LOWZ and CMASS galaxy samples in the Sloan Digital Sky Survey (SDSS) and SDSS-III BOSS survey. We detect clear anisotropic signals of IA due to redshift-space distortions. By combining measured IA statistics with the conventional galaxy clustering statistics, we obtain tighter constraints on the growth rate. The improvement is particularly prominent for the LRG, which is the brightest galaxy sample and known to be strongly aligned with underlying dark matter distribution; we obtain $f\sigma_8 = 0.5297^{ + 0.0310}_{ - 0.0316}$ (68\% C.L.) from the clustering-only analysis and $f\sigma_8 = 0.4871^{ + 0.0218}_{ - 0.0222}$ with clustering and IA, meaning $24\%$ improvement. The constraint is in good agreement with the prediction of general relativity, $f\sigma_8 = 0.4937 $ at $z=0.34$. For LOWZ and CMASS samples, the improvement by adding the IA statistics is found to be $9\%$ and $2\%$, respectively. Our results indicate that the contribution from IA statistics for cosmological constraints can be further enhanced by carefully selecting galaxies for a shape sample.

The "Condor Array Telescope" or "Condor" is a high-performance "array telescope" comprised of six apochromatic refracting telescopes of objective diameter 180 mm, each equipped with a large-format, very low-read-noise ($\approx 1.2$ e$^-$), very rapid-read-time ($< 1$ s) CMOS camera. Condor is located at a very dark astronomical site in the southwest corner of New Mexico, at the Dark Sky New Mexico observatory near Animas, roughly midway between (and more than 150 km from either) Tucson and El Paso. Condor enjoys a wide field of view ($2.29 \times 1.53$ deg$^2$ or 3.50 deg$^2$), is optimized for measuring both point sources and extended, very low-surface-brightness features, and for broad-band images can operate at a cadence of 60 s (or even less) while remaining sky-noise limited with a duty cycle near 100\%. In its normal mode of operation, Condor obtains broad-band exposures of exposure time 60 s over dwell times spanning dozens or hundreds of hours. In this way, Condor builds up deep, sensitive images while simultaneously monitoring tens or hundreds of thousands of point sources per field at a cadence of 60 s. Condor is also equipped with diffraction gratings and with a set of He II 468.6 nm, [O III] 500.7 nm, He I 587.5 nm, H$\alpha$ 656.3 nm, [N II] 658.4 nm, and [S II] 671.6 nm narrow-band filters, allowing it to address a variety of broad- and narrow-band science issues. Given its unique capabilities, Condor can access regions of "astronomical discovery space" that have never before been studied. Here we introduce Condor and describe various aspects of its performance.

Using optical photometric observations from 1.3m Devasthal Fast Optical Telescope and deep near-infrared (NIR) photometric observations from TANSPEC mounted on 3.6m Devasthal Optical Telescope, along with the multi-wavelength archival data, we present our study of open cluster Kronberger 55 to understand the star formation scenario in the region. The distance, extinction and age of the cluster Kronberger 55 are estimated as ~3.5 kpc, E(B-V)~1.0 mag and $\lesssim$55 Myr, respectively. We identified Young Stellar Objects (YSOs) based on their excess infrared (IR) emission using the two-color diagrams (TCDs). The mid-infrared (MIR) images reveal the presence of extended structure of dust and gas emission along with the outflow activities in the region with two peaks, one at the location of cluster Kronberger 55 and another at 5'.35 southwards to it. The association of radio continuum emission with the southern peak, hints towards the formation of massive star/s. The Herschel sub-millimeter maps reveal the presence of two clumps connected with a filamentary strcuture in this region, and such configuration is also evident in the 12CO(1-0) emission map. Our study suggests that this region might be a hub-filament system undergoing star formation due to the 'end-dominated collapse scenario'.

The abundances of deuterated molecules with respect to their main isotopologue counterparts have been determined to be orders of magnitude higher than expected from the cosmic abundance of deuterium relative to hydrogen. The increasing number of singly and multi-deuterated species detections helps us to constrain the interplay between gas-phase and solid-state chemistry and to understand better deuterium fractionation in the early stages of star formation. Acetaldehyde is one of the most abundant complex organic molecules (COMs) in star-forming regions and its singly deuterated isotopologues have already been observed towards protostars. A spectroscopic catalogue for astrophysical purposes is built for doubly deuterated acetaldehyde CHD2CHO from measurements in the laboratory. Submillimetre wave transitions were measured for the non-rigid doubly deuterated acetaldehyde CHD2CHO displaying hindered internal rotation of its asymmetrical CHD2 methyl group. A line position analysis is carried out allowing us to reproduce 853 transition frequencies with a weighted root mean square standard deviation of 1.7, varying 40 spectroscopic constants. A spectroscopic catalogue for astrophysical purposes is built from the analysis results. Using this catalogue we were able to detect for the first time CHD2CHO towards the low-mass protostellar system IRAS 16293-2422 utilizing data from the ALMA Protostellar Interferometric Line Survey. The first detection of the CHD2CHO species allows for the derivation of its column density with a value of 1.3x10^15 cm^-2 and an uncertainty of 10-20%. The resulting D2/D ratio of ~20% is found to be coincident with D2/D ratios derived for other complex organic molecules towards IRAS~16293-2422, pointing at a common formation environment with enhanced deuterium fractionation.

Multi-messenger high-energy astrophysics has currently achieved the potential to unravel the origin of cosmic rays and how sources accelerate them, their relation to the diffuse radiation in the extra-galactic space, and their role to forge their galaxies of origin while they wander in their magnetic fields for millions of years. Neutrino astronomy produced its major scientific milestone with the discovery by IceCube of a diffuse flux at energies above 60~TeV with intensity comparable to a predicted upper limit to the flux from extra-galactic sources of ultra-high energy cosmic rays. More recent results provide the first strong evidence of a standalone neutrino source and a highly probable coincidence of a neutrino alert with gamma rays. These results of IceCube indicate that neutrino astronomy can complement photon astronomy also providing insights into opaque sources of high-energy radiation. Starburst galaxies and jetted black holes in active galaxies are favored candidates to explain the diffuse cosmic neutrino background. Additionally, gamma-ray bursts remain an intriguing mystery now enriched by joint observations of gamma-rays and gravitational waves. The galactic diffuse flux, produced by cosmic ray interactions on the interstellar matter of our galaxy and peaking at lower energies, is within the reach of neutrino detectors. Together with the measured galactic gamma-ray flux up to PeV energies, they will shed light on the knee region of cosmic rays and the possible existence of dark matter in the Galactic plane. In the future, more work will be done in IceCube and deep sea and lake neutrino telescopes to use further low-energy cascades for cosmic source searches thanks to improved descriptions of detection media and deep learning methods.

The rising concern in the Hubble constant tension ($H_0$ tension) of the cosmological models motivates the scientific community to search for alternative cosmological scenarios that could resolve the $H_0$ tension. In this regard, we aim to work on a torsion-based modified theory of gravity which is an alternative description to the coherence model. We find an analytic solution of the Hubble parameter using a linear Lagrangian function of torsion $T$ and trace of energy-momentum tensor $\mathcal{T}$ for the dust case. Further, we constrain the cosmological and model parameters; to do that, we use Hubble and Pantheon samples and Markov Chain Monte Carlo (MCMC) simulation through Bayesian statistics. We obtain the values of Hubble constant as $H_0= 69.9\pm 6.8$ km$s^{-1}$ Mp$c^{-1}$, $H_0= 70.3\pm 6.3$ km$s^{-1}$ Mp$c^{-1}$, and $H_0= 71.4\pm 6.3$ km$s^{-1}$ Mp$c^{-1}$ at confidence level (CL), for Hubble, Pantheon, and their combine analysis, respectively. These outputs of $H_0$ for our model align with the recent observational measurements of $H_0$. In addition, we test the $Om$ diagnostic to check our model's dark energy profile.

We observed a full-orbit phase curve of the hot Jupiter WASP-43b with MIRI/LRS as part of the Transiting Exoplanet Community Early Release Science Program. Here we report preliminary findings for the instrument performance from the team's MIRI Working Group. Overall we find that MIRI's performance for phase curve observations is excellent, with a few minor caveats. The key takeaways for Cycle 2 planning with MIRI/LRS are: (1) long-duration observations (> 24 hours) have now been successfully executed; (2) for phase curves, we recommend including a one-hour burn-in period prior to taking science data to mitigate the effects of the ramp systematic; and (3) we do not yet recommend partial phase curve observations. In addition, we also find that: the position of the spectrum on the detector is stable to within 0.03 pixels over the full 26.5-hour observation; the light curves typically show a systematic downward ramp that is strongest for the first 30 minutes, but continues to decay for hours; from 10.6-11.8 microns, the ramp effect has remarkably different behavior, possibly due to a different illumination history for the affected region of the detector; after trimming the integrations most affected by the initial ramps and correcting the remaining systematics with analytic models, we obtain residuals to the light-curve fits that are typically within 25% of the photon noise limit for 0.5-micron spectroscopic bins; non-linearity correction is not a significant source of additional noise for WASP-43, though it may be an issue for brighter targets; the gain value of 5.5 electrons/DN currently on CRDS and JDox is known to be incorrect, and the current best estimate for the gain is approximately 3.1 electrons/DN; new reference files for the JWST calibration pipeline reflecting these findings are under development at STScI.

We present a rapid timing analysis of optical (HiPERCAM and ULTRACAM) and X-ray (NICER) observations of the X-ray transient Swift J1858.6-0814 during 2018 and 2019. The optical light curves show relatively slow, large amplitude (~1 mags in g$_s$) `blue' flares (i.e. stronger at shorter wavelengths) on time-scales of ~minutes as well as fast, small amplitude (~0.1 mag in g$_s$) `red' flares (i.e. stronger at longer wavelengths) on time-scales of ~seconds. The `blue' and `red' flares are consistent with X-ray reprocessing and optically thin synchrotron emission, respectively, similar to what is observed in other X-ray binaries. The simultaneous optical versus soft- and hard-band X-ray light curves show time- and energy dependent correlations. The 2019 March 4 and parts of the June data show a nearly symmetric positive cross correlations (CCFs) at positive lags consistent with simple X-ray disc reprocessing. The soft- and hard-band CCFs are similar and can be reproduced if disc reprocessing dominates in the optical and one component (disc or synchrotron Comptonization) dominates both the soft and hard X-rays. A part of the 2019 June data shows a very different CCFs. The observed positive correlation at negative lag in the soft-band can be reproduced if the optical synchrotron emission is correlated with the hot flow X-ray emission. The observed timing properties are in qualitative agreement with the hybrid inner hot accretion flow model, where the relative role of the different X-ray and optical components that vary during the course of the outburst, as well as on shorter time-scales, govern the shape of the optical/X-ray CCFs.

In this paper we discuss the implementation of a discrete, piecewise power-law grain-size distribution method into a numerical multifluid MHD code as described in Sumpter (2020). Such a description allows to capture the full size range of dust grains and their dynamical effects. The only assumptions are that grains within a single discrete bin have the same velocity and charge. We test the implementation by modelling plane-parallel C-type shocks and compare the results with shock models of multispecies grain models. We find that both the discrete and multispecies grain models converge to the same shock profile. However, the convergence for the discrete models is faster than for the multispecies grain models. For the pure advection models a single discrete bin is sufficient, while the multispecies grain models need a minimum of 8 grain species. When including grain sputtering the necessary number of discrete bins increases to 4, as the grain distribution cannot be described by a single power-law as in the advection models. The multispecies grain models still need more grain species to model the distribution, but the number does not increase compared to the pure advection models. Our results show that modelling the grain distribution function using a discrete distribution reduces the computational cost needed to capture the grain physics significantly.

In recent years there appeared a need for astronomical observations timed with sub-millisecond accuracy. These include e.g. timing stellar occultations by small, sub-km or fast Near Earth Asteroids, but also tracking artificial satellites at Low Earth Orbit using optical sensors. Precise astrometry of fast-moving satellites, and accurate timing of stellar occultations have parallel needs, requiring reliable time source and good knowledge of camera delays. Thus a need for an external device that would enable equipment and camera testing, to check if they reach the required accuracy in time. We designed, constructed and thoroughly tested a New EXposure Timing Analyser (NEXTA): a GNSS-based precise timer (Global Navigation Satellite System), allowing to reach the accuracy of 0.1 millisecond, which is an order of magnitude better than in previously available tools. The device is a simple strip of blinking diodes, to be imaged with a camera under test and compare imaged time with internal camera time stamp. Our tests spanned a range of scientific cameras widely used for stellar occultations and ground-based satellite tracking. The results revealed high reliability of both NEXTA and most of the tested cameras, but also pointed that practically all cameras had internal time bias of various level. NEXTA can serve the community, being easily reproducible with inexpensive components. We provide all the necessary schemes and usage instructions.

Astrophysical jets are relativistic outflows that remain collimated for remarkably many orders of magnitude. Despite decades of research, the origin of cosmic rays (CRs) remains unclear, but jets launched by both supermassive black holes in the centre of galaxies and stellar-mass black holes harboured in X-ray binaries (BHXBs) are among the candidate sources for CR acceleration. When CRs accelerate in astrophysical jets, they initiate particle cascades that form {\gamma}-rays and neutrinos. In the so-called hadronic scenario, the population of accelerated CRs requires a significant amount of energy to properly explain the spectral constraints similarly to a purely leptonic scenario. The amount of energy required often exceeds the Eddington limit, or even the total energy available within the jets. The exact energy source for the accelerated protons is unclear, but due to energy conservation along the jets, it is believed to come from the jet itself via transfer of energy from the magnetic fields, or kinetic energy from the outflow. To address this hadronic energy issue and to self-consistently evolve the energy flux along the flows, we explore a novel treatment for including hadronic content, in which instabilities along the jet/wind border play a critical role. We discuss the impact of the different jet composition on the jet dynamics for a pair dominated and an electron-proton jet, and consequently the emitted spectrum, accounting for both leptonic and hadronic processes. Finally, we discuss the implications of this mass-loading scenario to address the proton energy issue.

This study is done in the context of the project titled Interplanetary Laser Tri-lateration Network (ILTN) proposed by \cite{2018P&SS..153..127S} and investigated more in details by \cite{2022P&SS..21405415B} and \cite{2022P&SS..21505423B}. The original idea was to propose interplanetary measurements (in this case between Venus, Mars and the earth) as a way to measure the solar system expansion. But some recent interests on the measurement of asteroid masses and more generally the study of the mass distribution in the outer solar system appear with the ILTN. In this work, we are investigating how different possible configurations of interplanetary measurements of distances can be introduced in planetary ephemeris construction and how they improve our knowledge of planet orbits and other related parameters.

The relationship between galaxies and haloes is central to the description of galaxy formation, and a fundamental step towards extracting precise cosmological information from galaxy maps. However, this connection involves several complex processes that are interconnected. Machine Learning methods are flexible tools that can learn complex correlations between a large number of features, but are traditionally designed as deterministic estimators. In this work, we use the IllustrisTNG300-1 simulation and apply neural networks in a binning classification scheme to predict probability distributions of central galaxy properties, namely stellar mass, colour, specific star formation rate, and radius, using as input features the halo mass, concentration, spin, age, and the overdensity on a scale of 3 $h^{-1}$ Mpc. The model captures the intrinsic scatter in the relation between halo and galaxy properties, and can thus be used to quantify the uncertainties related to the stochasticity of the galaxy properties with respect to the halo properties. In particular, with our proposed method, one can define and accurately reproduce the properties of the different galaxy populations in great detail. We demonstrate the power of this tool by directly comparing traditional single-point estimators and the predicted joint probability distributions, and also by computing the power spectrum of a large number of tracers defined on the basis of the predicted colour-stellar mass diagram. We show that the neural networks reproduce clustering statistics of the individual galaxy populations with excellent precision and accuracy.

Extreme stripped-envelope supernovae (SESNe), including Type Ic superluminous supernovae (SLSNe), broad-line Type Ic SNe (SNe Ic-BL), and fast blue optical transients (FBOTs), are widely believed to harbor a newborn fast-spinning highly-magnetized neutron star (``magnetar''), which can lose its rotational energy via spin-down processes to accelerate and heat the ejecta. The progenitor(s) of these magnetar-driven SESNe, and the origin of considerable angular momentum (AM) in the cores of massive stars to finally produce such fast-spinning magnetars upon core-collapse are still under debate. Popular proposed scenarios in the literature cannot simultaneously explain their event rate density, SN and magnetar parameters, and the observed metallicity. Here, we perform a detailed binary evolution simulation that demonstrates that tidal spin-up helium stars with efficient AM transport mechanism in close binaries can form fast-spinning magnetars at the end of stars' life to naturally reproduce the universal energy-mass correlation of these magnetar-driven SESNe. Our models are consistent with the event rate densities, host environments, ejecta masses, and energetics of these different kinds of magnetar-driven SESNe, supporting that the isolated common-envelope formation channel could be a major common origin of magnetar-driven SESNe. The remnant compact binary systems of magnetar-driven SESNe are progenitors of some galactic systems and gravitational-wave transients.

Despite the fact that the progenitor of fast blue optical transients (FBOTs) is still up for debate, FBOTs are sometimes suggested to originate from the core-collapse of ultra-stripped stars and be powered by a spinning-down neutron star. Following this consideration, it is expected that the late-time evolution of the progenitor stars can leave important imprints in the circumstellar material (CSM) of the FBOTs, due to the strong mass loss of the stars. The interaction of the FBOT ejecta with the CSM can drive a long-lasting shock to generate radio emission, which thus enables us to probe the CSM properties through radio observation although such observations are still rare. Within the framework of the magnetar-powered model, Liu et al.(2022) fitted the multi-band optical light curves of 40 FBOTs and, hence, the statistical distributions of the FBOT magnetar and ejecta parameters were obtained. Based on these FBOT population results, we investigate the dependence of the radio emission on the mass-loss rate of the progenitors and evaluate the detectability of radio emission from FBOTs with current and future telescopes. It is found that the distribution of the peak time and peak luminosity of the emission at 8.4 GHz are primarily in the regions of $t_{\rm{peak},\nu}=10^{2.12\pm0.63}$ days and $L_{\rm{peak},\nu}=10^{28.73\pm0.83}$ erg s$^{-1}$ Hz$^{-1}$, respectively. A joint detection of the Zwicky Transient Facility and Very Large Array could achieve success in about 8.7% FBOTs of $z\leq1$. Furthermore, if considering a joint observation of the Chinese Space Station Telescope and Square Kilometer Array, this fraction can be increased to about 23.9%.

Fast radio bursts (FRBs) are one of the most exciting new mysteries of astrophysics. Their origin is still unknown, but recent observations seem to link them to soft gamma repeaters and, in particular, to magnetar giant flares (MGFs). The recent detection of a MGF at GeV energies by the Fermi Large Area Telescope (LAT) motivated the search for GeV counterparts to the >100 currently known FRBs. To date, none of these has a known gamma-ray counterpart. Taking advantage of more than 12 years of Fermi-LAT data, we perform a search for gamma-ray emission from almost all the reported repeating and non-repeating FRBs. We analyze on different time scales the Fermi-LAT data for each individual source separately and perform a cumulative analysis on the repeating ones. In addition, we perform the first stacking analysis at GeV energies of this class of sources in order to constrain the gamma-ray properties of the FRBs. The stacking analysis is a powerful method that allows for a possible detection from below-threshold FRBs providing important information on these objects. In this proceeding we present preliminary results of our study and we discuss their implications for the predictions of gamma-ray emission from this class of sources.

We use archival ALMA observations of the HCN and CO $J=1-0$ transitions, in addition to the radio continuum at 93 GHz, to assess the relationship between dense gas, star formation, and gas dynamics in ten, nearby (U)LIRGs and late-type galaxy centers. We frame our results in the context of turbulent and gravoturbulent models of star formation to assess if the HCN/CO ratio tracks the gravitationally-bound, star-forming gas in molecular clouds ($f_\mathrm{grav}$) at sub-kpc scales in nearby galaxies. We confirm that the HCN/CO ratio is a tracer of gas above $n_\mathrm{SF}\approx10^{4.5}$ cm$^{-3}$, but the sub-kpc variations in HCN/CO do not universally track $f_\mathrm{grav}$. We find strong evidence for the use of varying star formation density threshold models, which are able to reproduce trends observed in $t_\mathrm{dep}$ and $\epsilon_\mathrm{ff}$ that fixed threshold models cannot. Composite lognormal and powerlaw models outperform pure lognormal models in reproducing the observed trends, even when using a fixed powerlaw slope. The ability of the composite models to better reproduce star formation properties of the gas provides additional indirect evidence that the star formation efficiency per free-fall time is proportional to the fraction of gravitationally-bound gas.

Massive stars blow powerful winds and eventually explode as supernovae. By doing so, they inject energy and momentum in the circumstellar medium, which is pushed away from the star and piles up to form a dense and expanding shell of gas. The effect is larger when many massive stars are grouped together in bound clusters or associations. Large cavities form around clusters as a result of the stellar feedback on the ambient medium. They are called superbubbles and are characterised by the presence of turbulent and supersonic gas motions. This makes star clusters ideal environments for particle acceleration, and potential contributors to the observed Galactic cosmic ray intensity. The acceleration of particles at star clusters and in their surroundings may provide a major contribution to the observed CR flux. Moreover, it may explain the fine structures observed in the chemical composition of these particles, and possibly provide a solution to the puzzle of the origin of cosmic rays of energies in the PeV range and beyond.

Using relativistic supernova simulations of massive progenitor stars with a quark-hadron equation of state (EoS) and a purely hadronic EoS, we identify a distinctive feature in the gravitational-wave signal that originates from a buoyancy-driven mode (g-mode) below the proto-neutron star convection zone. The mode frequency lies in the range $200\lesssim f\lesssim 800\,\text{Hz}$ and decreases with time. As the mode lives in the core of the proto-neutron star, its frequency and power are highly sensitive to the EoS, in particular the sound speed around twice saturation density.

The Very Large Telescope Interferometer is one of the most proficient observatories in the world for high angular resolution. Since its first observations, it has hosted several interferometric instruments operating in various bandwidths in the infrared. As a result, the VLTI has yielded countless discoveries and technological breakthroughs. Here, we introduce a new concept for the VLTI, Asgard: an instrumental suite comprised of four natively collaborating instruments: BIFROST, a combiner whose main science case is studying the formation processes and properties of stellar and planetary systems; NOTT, a nulling interferometer dedicated to imaging young nearby planetary systems in the L band; HEIMDALLR, an all-in-one instrument performing both fringe tracking and stellar interferometry with the same optics; Baldr, a Strehl optimiser. These instruments share common goals and technologies. The goals are diverse astrophysical cases such as the study of the formation and evolution processes of binary systems, exoplanetary systems and protoplanetary disks, the characterization of orbital parameters and spin-orbit alignment of multiple systems, the characterization of the exoplanets, and the study of exozodiacal disks. Thus, the idea of this suite is to make the instruments interoperable and complementary to deliver unprecedented sensitivity and accuracy from the J to M bands to meet these goals. The interoperability of the Asgard instruments and their integration in the VLTI are major challenges for this project.

Measurements of the plasma parameters of coronal mass ejections (CMEs), particularly the magnetic field and non-thermal electron population entrained in the CME plasma, are crucial to understand their propagation, evolution, and geo-effectiveness. Spectral modeling of gyrosynchrotron (GS) emission from CME plasma has been regarded as one of the most promising remote sensing technique for estimating spatially resolved CME plasma parameters. Imaging the very low flux density CME GS emission in close proximity to the Sun with orders of magnitude higher flux density, however, has proven to be rather challenging. This challenge has only recently been met using the high dynamic range imaging capability of the Murchison Widefield Array (MWA). Although routine detection of GS is now within reach, the challenge has shifted to constraining the large number of free parameters in GS models, a few of which are degenerate, using the limited number of spectral points at which the observations are typically available. These degeneracies can be broken using polarimetric imaging. For the first time, we demonstrate this using our recently developed capability of high fidelity polarimetric imaging on the data from the MWA. We show that, in addition to breaking the degeneracies, spectro-polarimetric spectroscopic imaging also yields tighter constraints on the plasma parameters of key interest than possible with total intensity spectroscopic imaging alone.

In order to provide spectral ground truth data for remote sensing applications, we have measured midinfrared spectra (2 to 18 micron) of three typical, well defined lithologies from the Chelyabinsk meteorite. These lithologies are classified as (a) moderately shocked, light lithology, (b) shock darkened lithology, and (c) impact melt lithology. Analyses were made from bulk material in four size fractions (0 to 25 micron, 25 to 63 micron, 63 to 125 micron, and 125 to 250 micron), and from additional thin sections. Characteristic infrared features in the powdered bulk material of the moderately shocked, light lithology, dominated by olivine, pyroxene and feldspathic glass, are a Christiansen feature (CF) between 8.5 and 8.8 micron; a transparency feature (TF) in the finest size fraction at about 13 micron, and strong reststrahlen bands (RB) at about 9.1 micron, 9.5 micron, 10.3 micron, 10.8 micron, 11.2 to 11.3 micron, 12 micron, and between 16 and 17 micron. The ranges of spectral features for the micro FTIR spots show a wider range than those obtained in diffuse reflectance, but are generally similar. With increasing influence of impact shock from pristine LL5 (or LL6) material (which have a low or moderate degree of shock) to the shock-darkened lithology and the impact melt lithology as endmembers, we observe the fading or disappearing of spectral features. Most prominent is the loss of a twin peak feature between 10.8 and 11.3 micron, which turns into a single peak. In addition, in the pure impact melt endmember lithology features at about 9.6 micron and about 9.1 micron are also lost. These losses are most likely correlated with decreasing amounts of crystal structure as the degree of shock melting increases.

56 UMa is a wide binary system containing a chemically peculiar red giant and a faint companion. Due to its surface chemical abundances, the red giant was classified as a Barium (Ba) star. This implies that the companion has to be a white dwarf, since Ba stars form when mass is transferred to them from an s-process rich Asymptotic Giant Branch (AGB) star. However, in the case of 56 UMa, the companion might be too massive to be the progeny of an AGB star that efficiently produced s-process elements like barium. In this letter, we revisit the orbital parameters of the system and perform a full spectral analysis with the goal of investigating the Ba-star classification of the giant and unravelling the nature of its faint companion. We combine radial-velocity and astrometric data to refine the orbital parameters of the system, including the orbital inclination and the companion mass. Then, we redetermine the stellar parameters of the giant and its chemical abundances using high-resolution HERMES spectra. Finally, we investigate the morphology of the interstellar gas in the vicinity of the system. The faint component in 56 UMa has a mass of $1.31 \pm 0.12$ M$_{\odot}$, which, together with the mixed s+r abundance profile of the red giant, confirms that the giant is not a standard Barium star. Additionally, the clear identification of a cavity surrounding 56 UMa could indicate that a supernova explosion occurred about 10$^5$ years ago in the system, suggesting that the faint companion might be a neutron star. However, finding an evolutionary scenario that explains all the observables is not trivial, so we discuss different possible configurations of the system and their respective merits.

Mrk 1239 is a highly polarized NLS1 in the optical band, whose $0.3-3$ keV spectrum has remained remarkably consistent over more than two decades of observation. Previous analysis of this object suggested that the soft X-ray band was dominated by emission lines (collisionally and/or photoionized) from the distant host galaxy as the X-ray emission from the central engine was highly obscured. New XMM-Newton data of Mrk 1239 are presented here to investigate the soft X-ray band of this galaxy with high resolution. The first RGS spectra of this source reveal a plethora of ionized emission lines originating from two distinct plasmas, one collisionally ionized and the other photoionized at approximately equal brightness. The best fit model uses {\sc apec} and {\sc xstar} grids to account for the collisionally ionized and photoionized components, respectively. The fit improves significantly if the photoionized material is allowed to outflow at $\approx 500$ km s$^{-1}$, matching the outflow velocity of the forbidden O{\sc vii} emission line. From constraints on the ionization and density of the photoionized material we can estimate the location of it to be no further than a few pc from the central source, around the outer radius of the torus, which is consistent with the O{\sc vii}$(f)$ emission line. Properties of the collisionally ionized plasma are consistent with star formation rate (SFR) of $\approx 3 M_{\odot} \textrm{yr}^{-1}$, which is comparable with several previous measurements of the SFR in this galaxy.

We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-{\gamma} line-emission and Lp (3.8 {\mu}m) imaging data, we provide the first estimate of its orbital parameters and quantitative characterization of the evolution of its morphology and mass. We find that the leading edge of X7 appears to be on a mildly eccentric (e~0.3), relatively short-period (170 years) orbit and is headed towards periapse passage, estimated to occur in ~2036. Furthermore, our kinematic measurements rule out the earlier suggestion that X7 is associated with the stellar source S0-73 or with any other point source that has overlapped with X7 during our monitoring period. Over the course of our observations, X7 has (1) become more elongated, with a current length-to-width ratio of 9, (2) maintained a very consistent long-axis orientation (position angle of 50 deg), (3) inverted its radial velocity differential from tip to tail from -50 to +80 km/sec, and (4) sustained its total brightness (12.8 Lp magnitudes at the leading edge) and color temperature (425 K), which suggest a constant mass of ~50 MEarth. We present a simple model showing that these results are compatible with the expected effect of tidal forces exerted on it by the central black hole and we propose that X7 is the gas and dust recently ejected from a grazing collision in a binary system.

We present a numerical study of relativistic Bondi-Hoyle-Lyttleton (BHL) accretion onto an asymptotically flat black hole with synchronized hair. The hair is sourced by an ultralight, complex scalar field, minimally coupled to Einstein's gravity. Our simulations consider a supersonic flow parametrized by the asymptotic values of the fluid quantities and a sample of hairy black holes with different masses, angular momenta, and amount of scalar hair. For all models, steady-state BHL accretion solutions are attained that are characterized by the presence of a shock-cone and a stagnation point downstream. For the models of the sample with the largest component of scalar field, the shock-cone envelops fully the black hole, transitioning into a bow-shock, and the stagnation points move further away downstream. Analytical expressions for the mass accretion rates are derived, which can be used to analyze black-hole formation scenarios in the presence of ultralight scalar fields. The formation of a shock-cone leads to regions where sound waves can be trapped and resonant oscillations excited. We measure the frequencies of such quasi-periodic oscillations and point out a possible association with quasi-periodic oscillations in the X-ray light curve of Sgr~A* and microquasars.

The origin of the diffuse flux of TeV-PeV astrophysical neutrinos is still unknown. The $\gamma$-ray blazar PKS 0735+178, located within a 90% error region expanded by a factor of 1.3 of the neutrino event IC-211208A, was found to be flaring in multi-frequency wavebands at the time of detection by IceCube Observatory. In addition to leptonic synchrotron (SYN) and synchrotron self-Compton (SSC) emission, we invoke photohadronic (p$\gamma$) interactions inside the jet to model the spectral energy distribution (SED) and neutrino emission. We analyze the long-term $\gamma$-ray and X-ray data to generate the broadband SED. The temporal light curve indicates that the source was in a high state in optical UV, $\gamma$-ray, and X-ray range during the neutrino detection epoch. In the one-zone lepto-hadronic model, the SSC photons do not provide enough seed photons for p$\gamma$ interactions to explain the neutrino event. However, including an external photon field yields a neutrino event rate of 0.2 in 100 days, for the IceCube detector, for physically motivated values of the magnetic field, an external photon field peaking at optical wavelength, and other jet parameters without exceeding the Eddington luminosity. The radiation from secondary electrons at X-ray energies severely constrains the neutrino flux to a lower value than found in previous studies. In contrast, the pion decay cascade flux at GeV energies is subdominant at the high-energy peak of the SED, suggesting a higher correlation of neutrinos with X-ray emission is plausible.

Accreting binary black holes (BBHs) are multi-messenger sources, emitting copious electromagnetic (EM) and gravitational waves. One of their most promising EM signatures is the lightcurve modulation caused by a strong, unique and extended azimuthal overdensity structure orbiting at the inner edge of the circumbinary disc (CBD), dubbed "lump". In this paper, we investigate the origin of this structure using 2D general-relativistic (GR) hydrodynamical simulations of a CBD in an approximate BBH spacetime. First, we use the symmetric mass-ratio case to study the transition from the natural m = 2 mode to m = 1. The asymmetry with respect to m = 2 grows exponentially, pointing to an instability origin. We indeed find that the CBD edge is prone to a (magneto-)hydrodynamical instability owing to the disc edge density sharpness: the Rossby Wave Instability (RWI). The RWI criterion is naturally fullfilled at the CBD edge and we report the presence of vortices, which are typical structures of the RWI. The RWI is also at work in the asymmetric mass-ratio cases (from 0.1 to 0.5). However, the CBD edge sharpness decreases with a decreasing mass ratio, and so the lump. By proposing a scenario for this lump formation, our work further supports its existence in astrophysical CBDs and potential source for an EM signature of BBHs. Finally, because the RWI is not caused by GR effects, it is also a robust candidate for the lump origin in CBDs around non-compact objects, e.g. binary protostars.

Retrieving the physical parameters from spectroscopic observations of exoplanets is key to understanding their atmospheric properties. Exoplanetary atmospheric retrievals are usually based on approximate Bayesian inference and rely on sampling-based approaches to compute parameter posterior distributions. Accurate or repeated retrievals, however, can result in very long computation times due to the sequential nature of sampling-based algorithms. We aim to amortize exoplanetary atmospheric retrieval using neural posterior estimation (NPE), a simulation-based inference algorithm based on variational inference and normalizing flows. In this way, we aim (i) to strongly reduce inference time, (ii) to scale inference to complex simulation models with many nuisance parameters or intractable likelihood functions, and (iii) to enable the statistical validation of the inference results. We evaluate NPE on a radiative transfer model for exoplanet spectra petitRADTRANS, including the effects of scattering and clouds. We train a neural autoregressive flow to quickly estimate posteriors and compare against retrievals computed with MultiNest. NPE produces accurate posterior approximations while reducing inference time down to a few seconds. We demonstrate the computational faithfulness of our posterior approximations using inference diagnostics including posterior predictive checks and coverage, taking advantage of the quasi-instantaneous inference time of NPE. Our analysis confirms the reliability of the approximate posteriors produced by NPE. The accuracy and reliability of the inference results produced by NPE establishes it as a promising approach for atmospheric retrievals. Amortization of the posterior inference makes repeated inference on several observations computationally inexpensive since it does not require on-the-fly simulations, making the retrieval efficient, scalable, and testable.

We present an homogeneous analysis of all DA stars labeled as magnetic in the Montreal White Dwarf Database (MWDD). Our sample is restricted to almost all known magnetic white dwarf showing clear sign of splitting ($B \gtrsim$ 1-2 MG) that have parallax measurements from the second Gaia data release, photometric data from diverse surveys and spectroscopic data from SDSS or archival data from the Montreal group. We determine the atmospheric parameters (effective temperature, surface gravity, magnetic field strength/geometry) of all objects using state-of-the-art model atmosphere/magnetic synthetic spectra, as well as reclassify many objects that were prematurely labeled as potentially magnetic. Finally, we discuss the atmospheric parameters/field properties distribution as well as the implication on our understanding of magnetic white dwarfs origin and evolution.

We describe new functionality in the GYRE stellar oscillation code for modeling tides in binary systems. Using a multipolar expansion in space and a Fourier-series expansion in time, we decompose the tidal potential into a superposition of partial tidal potentials. The equations governing the small-amplitude response of a spherical star to an individual partial potential are the linear, non-radial, non-adiabatic oscillation equations with an extra inhomogeneous forcing term. We introduce a new executable, gyre_tides, that directly solves these equations within the GYRE numerical framework. Applying this to selected problems, we find general agreement with results in the published literature but also uncover some differences between our direct solution methodology and the modal decomposition approach adopted by many authors. In its present form gyre_tides can model equilibrium and dynamical tides of aligned binaries in which radiative diffusion dominates the tidal dissipation (typically, intermediate and high-mass stars on the main sequence). Milestones for future development include incorporation of other dissipation processes, spin-orbit misalignment, and the Coriolis force arising from rotation.

We present an analysis of all single white dwarf stars known to exhibit spectroscopic signatures of neutral helium line splitting due to the presence of a strong magnetic field. Using state-of-the-art models taking into account the effects of magnetic fields on the synthetic spectra, we determine effective temperatures, surface gravities and masses for the stars in our sample. Our analysis uses data from the second and third Gaia (early) data release, photometric data from diverse surveys such as the Sloan Digital Sky Survey and Pan-STARRS, and archived spectroscopic data. We are able to successfully reproduce the spectra of 8 objects using an offset dipole geometry while several others seem to require either a more complexe geometry or a different chemical composition. We also highlight a group of hot featureless white dwarfs that are most probably highly magnetic objects whose spectra are completely smeared due to the field strength distribution across the surface.

Stellar bars have been found to substantially influence the stellar populations properties in galaxies, affecting their ability of forming stars. While this can be easily seen when studying galaxies in relatively isolated environments, such kind of analysis takes a higher degree of complexity when cluster galaxies are considered, due to the variety of interactions which can potentially occur in these denser environments. We use IFU MUSE data from the GASP survey to study the combined effect of the presence of a stellar bar and of ram pressure, on spatially resolved properties of stellar populations. We have analyzed spatially resolved indicators of both recent SFR and average stellar population ages to check for signatures of anomalous central SF activity, also taking into account for the possible presence of nuclear activity. We found an increase of central SFR in ram pressure affected galaxies when compared with unperturbed ones. The most extreme cases of increase SFR and central rejuvenation occur in barred galaxies that are at advanced stages of ram pressure stripping. For low-mass barred galaxies affected by ram pressure, the combined effect is a systematic enhancement of the star formation activity as opposed to the case of high-mass galaxies that present both enhancement and suppression. Barred galaxies that present a suppression of their star formation activity also present signatures of nuclear activity. Our results indicate that the combined effect of the presence of a bar and a strong perturbation by ram pressure is able to trigger the central SF activity and probably ignite nuclear activity.

The Local Group (LG) of galaxies, modeled as a two body problem, is sensitive to cosmological contributions like those related to the presence of a cosmological constant $\Lambda$ into dynamics. Here we study the LG dynamics in the context of Extended Theories of Gravity like $f(R)$ gravity considered as dark energy and dark matter contributions. In the first approach, we perturb the dark energy effect considering a Yukawa-like interaction that naturally emerges from $f(R)$ gravity in the weak field limit. We assume the mass of LG from simulations and, from this, derive constraints on the Yukawa couplings: $\alpha < 0.581$ and $m_{grav} < 5.095 \cdot 10^{-26} \, eV/c^2$. In the second part, considering a minimal extension of General Relativity, i.e. $f(R) \sim R^{1+\epsilon}$, with $|\epsilon|\ll 1$, we investigate the possibility that it replaces dark matter as a MOND-like theory. We find that there is a value of the parameter $\beta$ (derived starting from $\epsilon$) which gives a minimal value for the LG mass. Moreover, this particular potential allows to calculate the ratio of dark matter and baryonic matter for the LG to be. We show that this ratio could falsify MOND-like theories.

RS Oph is a symbiotic recurrent nova containing a massive white dwarf with heavy mass loss during activity. In August 2021, it underwent its seventh optical eruption since the end of the 19th century. The goal of this work is to analyse the structure of the outflows from the outbursting object. Based on broad-band $U$, $B$, $V$, $R_{\rm C}$, and $I_{\rm C}$ photometry and high-resolution H$\alpha$ spectroscopy obtained at days 11--15 of the outburst, we derived some parameters of the system's components and outflows and their changes during our observation. The effective temperature of a warm shell (pseudophotosphere) produced by the ejected material and occulting the hot component of the system was $T_{eff}=15000\pm1000$ K and the electron temperature of the nebula was $T_{e}=17000\pm3000$ K throughout the observations. The effective radius of the pseudophotosphere was $R_{ eff}=13.3\pm2.0$ R$_{\odot}$ and the emission measure of the nebula $EM=(9.50\pm0.59) $10$^{61}$ cm$^{-3}$ for day 11 and $R_{eff}=10.3\pm1.6$ R$_{\odot}$ and $EM=(5.60\pm0.35)$10$^{61}$ cm$^{-3}$ for day 15. To provide this emission measure, the bolometric luminosity of the outbursting object must exceed its Eddington limit. The mass-loss rate of the outbursting object through its wind is much greater than through its streams. The total rate (from wind + streams) was less than $(4-5)$ 10$^{-5}$ (d/1.6kpc)$^{3/2}$ M$_{\odot}$yr$^{-1}$. The streams are not highly collimated. Their mean outflowing velocities are $\upsilon_{b}=-3680\pm60$ km s$^{-1}$ for the approaching stream and $\upsilon_{r}=3520\pm50$ km s$^{-1}$ for the receding one if the orbit inclination is 50$^\circ$.

Neutral hydrogen ($\rm{HI}$) $21$-cm intensity mapping (IM) offers an efficient technique for mapping the large-scale structures in the universe. We introduce the 'Cross' Tapered Gridded Estimator (Cross TGE), which cross-correlates two cross-polarizations (RR and LL) to estimate the multi-frequency angular power spectrum (MAPS) $C_{\ell}(\Delta\nu)$. We expect this to mitigate several effects like noise bias, calibration errors etc., which affect the 'Total' TGE which combines the two polarizations. Here we apply the Cross TGE on a $24.4 \,\rm{MHz}$ bandwidth uGMRT Band $3$ data centred at $432.8 \,\rm{MHz}$ aiming $\rm{HI}$ IM at $z=2.28$. The measured $C_{\ell}(\Delta\nu)$ is modelled to yield maximum likelihood estimates of the foregrounds and the spherical power spectrum $P(k)$ in several $k$ bins. Considering the mean squared brightness temperature fluctuations, we report a $2\sigma$ upper limit $\Delta_{UL}^{2}(k) \le (58.67)^{2} \, {\rm mK}^{2}$ at $k=0.804 \, {\rm Mpc}^{-1}$ which is a factor of $5.2$ improvement on our previous estimate based on the Total TGE. Assuming that the $\rm{HI}$ traces the underlying matter distribution, we have modelled $C_{\ell}(\Delta\nu)$ to simultaneously estimate the foregrounds and $[\Omega_{\rm{HI}} b_{\rm{HI}}] $ where $\Omega_{\rm{HI}}$ and $b_{\rm{HI}}$ are the $\rm{HI}$ density and linear bias parameters respectively. We obtain a best fit value of $[\Omega_{\rm{HI}}b_{\rm{HI}}]^2 = 7.51\times 10^{-4} \pm 1.47\times 10^{-3}$ which is consistent with noise. Although the $2\sigma$ upper limit $[\Omega_{\rm{HI}}b_{\rm{HI}}]_{UL} \leq 0.061$ is $\sim 50$ times larger than the expected value, this is a considerable improvement over earlier works at this redshift.

We utilize medium-resolution JWST/NIRSpec observations of 164 galaxies at $z=2.0-9.3$ from the Cosmic Evolution Early Release Science (CEERS) survey to investigate the evolution of the excitation and ionization properties of galaxies at high redshifts. Our results represent the first statistical constraints on the evolution of the [OIII]/H$\beta$ vs. [NII]/H$\alpha$, [SII]/H$\alpha$, and [OI]/H$\alpha$ ``BPT'' diagrams at $z>2.7$, and the first analysis of the O32 vs. R23 diagram at $z>4$ with a large sample. We divide the sample into five redshift bins containing 30-40 galaxies each. The subsamples at $z\sim2.3$, $z\sim3.3$, and $z\sim4.5$ are representative of the main-sequence star-forming galaxy population at these redshifts, while the $z\sim5.6$ and $z\sim7.5$ samples are likely biased toward high specific star-formation rate due to selection effects. Using composite spectra, we find that each subsample at $z=2.0-6.5$ falls on the same excitation sequence in the [NII] and [SII] BPT diagrams and the O32-R23 diagram on average, offset from the sequences followed by $z=0$ HII regions in the same diagrams. The direction of these offsets are consistent with high-redshift star-forming galaxies uniformly having harder ionizing spectra than typical local galaxies at fixed nebular metallicity. The similarity of the average line ratios suggests that the ionization conditions of the interstellar medium do not strongly evolve between $z\sim2$ and $z\sim6$. Overall, the rest-optical line ratios suggest the $z=2.7-9.3$ CEERS/NIRSpec galaxies at log($M_*/M_{\odot})\sim7.5-10$ have high degrees of ionization and moderately low oxygen abundances ($\sim0.1-0.3~Z_{\odot}$), but are not extremely metal poor ($<0.1~Z_{\odot}$) even at $z>6.5$.

Backed by advances in digital electronics, signal processing, computation, and storage technologies, aperture arrays, which had strongly influenced the design of telescopes in the early years of radio astronomy, have made a comeback. Amid all these developments, an international effort to design and build the world's largest radio telescope, the Square Kilometre Array (SKA), is ongoing. With its vast collecting area of 1 sq-km, the SKA is envisaged to provide unsurpassed sensitivity and leverage technological advances to implement a complex receiver to provide a large field of view through multiple beams on the sky. Many pathfinders and precursor aperture array telescopes for the SKA, operating in the frequency range of 10-300 MHz, have been constructed and operationalized to obtain valuable feedback on scientific, instrumental, and functional aspects. This review article looks explicitly into the progression of digital-receiver architecture from the Murchison Widefield Array (precursor) to the SKA1-Low. It highlights the technological advances in analog-to-digital converters (ADCs),field-programmable gate arrays (FPGAs), and central processing unit-graphics processing unit (CPU-GPU) hybrid platforms around which complex digital signal processing systems implement efficient channelizers, beamformers, and correlators. The article concludes with a preview of the design of a new generation signal processing platform based on radio frequency system-on-chip (RFSoC).

The low-frequency radio telescope of the Square kilometre Array (SKA) is being built by the international radio astronomical community to i) have orders of magnitude higher sensitivity and ii) be able to map the sky several hundred times faster, than any other existing facilities over the frequency range 50 MHz to 350 MHz. The sensitivity of a radio telescope array is in general dependent upon the number of electromagnetic sensors used to receive the sky signal. The total number of them is further constrained by the effects of mutual coupling between the sensor elements, allowable grating lobes in their radiation patterns etc. The operating frequency band is governed by the desired spatial and spectral response, acceptable sidelobe and backlobe levels, radiation efficiency, polarization purity and calibratability of sensors' response. This paper presents a brief review of several broadband antennas considered as potential candidates by various engineering groups across the globe, for the low frequency radio telescope of SKA covering the frequency range of 50-350 MHz, on the basis of their suitability for conducting primary scientific objectives.

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described phenomena in nature, its core prediction of an Alfvenic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report the first observational evidence for the Alfvenic weak-to-strong transition in MHD turbulence in the terrestrial magnetosheath using the four Cluster spacecraft. The observed transition indicates the universal existence of strong turbulence regardless of the initial level of MHD fluctuations. Moreover, the observations demonstrate that the nonlinear interactions of MHD turbulence play a crucial role in the energy cascade, widening the directions of the energy cascade and broadening the fluctuating frequencies. Our work takes a critical step toward understanding the complete picture of turbulence cascade, connecting the weak and strong MHD turbulence systems. It will have broad implications in star formation, energetic particle transport, turbulent dynamo, and solar corona or solar wind heating.

The standard model of cosmology ($\Lambda$CDM) is facing a serious crisis caused by the inconsistencies in the measurements of some fundamental cosmological parameters (Hubble constant $H_{0}$ and cosmic curvature parameter $\Omega_{k}$ for example). On the other hand a strictly linear evolution of the cosmological scale factor is found to be an excellent fit to a host of observations. Any model that can support such a coasting presents itself as a falsifiable model as far as the cosmological tests are concerned. In this article the observational data of strong gravitational lensing (SGL) systems from SLACS, BELLS, LSD and SL2S surveys has been used to test the viability of linearly coasting cosmology. Assuming the spherically symmetric mass distribution in lensing galaxies, the ratio of angular diameter distance from lens to source and angular diameter distance of the source is evaluated and is used to constrain the power law cosmology. It is found that the linear coasting is consistent with the SGL data within 1-$\sigma$; in agreement with various independent studies.

The density of cosmic rays inside molecular clouds determines the ionization rate in the dense cores where stars form. It is also one of the drivers of astrochemistry leading to the creation of complex molecules. Through Fermi Large Area Telescope observations of nearby giant molecular clouds, we observed deficits (holes) in the gamma-ray residual map when modelling with the expected gamma-ray diffuse emission from uniform cosmic rays interacting with the molecular content. We propose that the deficit is due to the lack of penetration of the low-energy (sub-GeV to GeV) cosmic rays into denser regions or clumps. This differs from the prevailing view of fast cosmic ray transport in giant molecular clouds where the magnetic turbulence is suppressed by neutral-ion damping, as our results require a slow diffusion inside dense molecular clumps. Through modelling we find that while the shielding is negligible on the cloud scale, it becomes important in the denser, parsec-sized regions where the gravitational collapse is already at play, changing the initial condition of star formation and astrochemistry.

Gravitational waves (GWs) from the compact binary coalescences can be used as standard sirens to explore the cosmic expansion history. In the next decades, it is anticipated that we could obtain the multi-band GW standard siren data (from nanohertz to a few hundred hertz), which are expected to play an important role in cosmological parameter estimation. In this work, we give for the first time the joint constraints on cosmological parameters using the future multi-band GW standard siren observations. We simulate the multi-band GW standard sirens based on the SKA-era pulsar timing array (PTA), the Taiji observatory, and the Cosmic Explorer (CE) to perform cosmological analysis. In the $\Lambda$CDM model, we find that the joint PTA+Taiji+CE data could provide a tight constraint on the Hubble constant with a $0.5\%$ precision. Moreover, PTA+Taiji+CE could break the cosmological parameter degeneracies generated by CMB, especially in the dynamical dark energy models. When combining the PTA+Taiji+CE data with the CMB data, the constraint precisions of $\Omega_{\rm m}$ and $H_0$ are $1.0\%$ and $0.3\%$, meeting the standard of precision cosmology. The joint CMB+PTA+Taiji+CE data give $\sigma(w)=0.028$ in the $w$CDM model and $\sigma(w_0)=0.11$ and $\sigma(w_a)=0.32$ in the $w_0w_a$CDM model, which are comparable with or close to the latest constraint results by the CMB+BAO+SN. In conclusion, it is worth expecting to use the future multi-band GW observations to explore the nature of dark energy and measure the Hubble constant.

Studying Doppler shifts provides deeper insights into the flow of mass and energy in the solar atmosphere. We perform a comprehensive measurement of Doppler shifts in the transition region and its center-to-limb variation (CLV) in the strong field regions ($|\textbf{B}| \geq$ 50 G) of 50 active regions (ARs), using the \ion{Si}{4} 1394~{\AA} line recorded by the Interface Region Imaging Spectrometer(IRIS). To locate the ARs and identify strong field regions, we have used the magnetograms obtained by the Helioseismic and Magnetic Imager (HMI). We find that in strong field regions, on average, all the ARs show mean redshifts ranging between 4{--}11~ km/s, which varies with ARs. These flows show a mild CLV, with sizable magnitudes at the limb and substantial scatter at the mid-longitude range. Our observations do not support the idea that redshifts in the lower transition region (T $<\sim$ 0.1 MK) are produced by field-aligned downflows as a result of impulsive heating and warrant alternative interpretation, such as downflow of type-\rm{II} spicules in the presence of a chromospheric wall created by cooler type-\rm{I} spicules.

As humanity has begun to explore space, the significance of space weather has become apparent. It has been established that coronal holes, a type of space weather phenomenon, can impact the operation of aircraft and satellites. The coronal hole is an area on the sun characterized by open magnetic field lines and relatively low temperatures, which result in the emission of the solar wind at higher than average rates. In this study, To prepare for the impact of coronal holes on the Earth, we use computer vision to detect the coronal hole region and calculate its size based on images from the Solar Dynamics Observatory (SDO). We then implement deep learning techniques, specifically the Long Short-Term Memory (LSTM) method, to analyze trends in the coronal hole area data and predict its size for different sun regions over 7 days. By analyzing time series data on the coronal hole area, this study aims to identify patterns and trends in coronal hole behavior and understand how they may impact space weather events. This research represents an important step towards improving our ability to predict and prepare for space weather events that can affect Earth and technological systems.

Filaments are an important part of star-forming interstellar clouds. Theie properties hold clues to their formation mechanisms and role in the star-formation process. We compare the properties of filaments in the Orion Molecular Cloud 3 (OMC-3), as seen in mid-infrared (MIR) absorption and far-infrared (FIR) dust emission. We calculated optical depth maps of the OMC-3 filaments based on the MIR absorption seen in Spitzer data and FIR dust emission observed with Herschel and the ArT\'eMiS instrument. The widths of the selected OMC-3 filament segments are in the range 0.03-0.1 pc, with similar average values seen in both MIR and FIR analyses. Compared to the widths, the individual parameters of the fitted Plummer functions are much more uncertain. The asymptotic power-law index has typically values p~3 but with a large scatter. Modelling shows that the FIR observations can systematically overestimate the filament widths. The effect is potentially tens of per cent at column densities above N(H$_2$) ~ $10^{22}$ cm$^{-2}$ but is reduced in more intense radiation fields, such as the Orion region. Spatial variations in dust properties could cause errors of similar magnitude. In the MIR analysis, dust scattering should generally not be a significant factor, unless there are high-mass stars nearby or the dust MIR scattering efficiency is higher than in the tested dust models. Thermal MIR dust emission can be a more significant source of error, especially close to embedded sources. The analysis of interstellar filaments can be affected by several sources of systematic error, but mainly at high column densities and, in the case of FIR observations, in weak radiation fields. The widths of the OMC-3 filaments were consistent between the MIR and FIR analyses and did not reveal systematic dependence on the angular resolution of the observations.

Radio holographic measurements using the MeerKAT telescope are presented for each of its supported observing bands, namely UHF (544--1087 MHz), L (856--1711 MHz) and S (1750--3499 MHz). Because the UHF-band receiver design is a scaled version of that of the L band, the electromagnetic performance in these two bands are expectedly similar to one another. Despite also being linearly polarized, S-band receivers have an entirely different design and distinct performance characteristics from the lower two bands. As introduced in previous work for the L band, evidence of higher-order waveguide mode activation also appears in S-band measurements but there are differences in its manifestation. Frequency-dependent pointing (beam squint), beam width, beam ellipticity, errorbeam, instrumental polarization and cross-polarization power measurements are illustrated for each of MeerKAT's observational bands in a side-by-side style to facilitate the comparison of features. The derivation of collimation errors and main reflector surface errors from measurements made at these relatively low observation frequencies is also discussed. Results include elevation and ambient temperature effects on collimation, as well as the signatures of collimation degrading over time. The accompanying data release includes a snapshot of full Jones matrix primary beam patterns for all bands and antennas, with corresponding derived metrics.

Aims. Our aim is to derive a fast and accurate method for computing the gravitational potential of astrophysical objects with high contrasts in density, for which nested or adaptive meshes are required. Methods. We present an extension of the convolution method for computing the gravitational potential to the nested Cartesian grids. The method makes use of the convolution theorem to compute the gravitational potential using its integral form. Results. A comparison of our method with the iterative outside-in conjugate gradient and generalized minimal residual methods for solving the Poisson equation using nonspherically symmetric density configurations has shown a comparable performance in terms of the errors relative to the analytic solutions. However, the convolution method is characterized by several advantages and outperforms the considered iterative methods by factors 10--200 in terms of the runtime, especially when graphics processor units are utilized. The convolution method also shows an overall second-order convergence, except for the errors at the grid interfaces where the convergence is linear. Conclusions. High computational speed and ease in implementation can make the convolution method a preferred choice when using a large number of nested grids. The convolution method, however, becomes more computationally costly if the dipole moments of tightly spaced gravitating objects are to be considered at coarser grids.

We present the observations and analysis of a high-magnification microlensing planetary event, KMT-2022-BLG-0440, for which the weak and short-lived planetary signal was covered by both the KMTNet survey and follow-up observations. The binary-lens models with a central caustic provide the best fits, with a planet/host mass ratio, $q = 0.75$--$1.00 \times 10^{-4}$ at $1\sigma$. The binary-lens models with a resonant caustic and a brown-dwarf mass ratio are both excluded by $\Delta\chi^2 > 70$. The binary-source model can fit the anomaly well but is rejected by the ``color argument'' on the second source. From Bayesian analyses, it is estimated that the host star is likely a K or M dwarf located in the Galactic disk, the planet probably has a Neptune-mass, and the projected planet-host separation is $1.9^{+0.6}_{-0.7}$ or $4.6^{+1.4}_{-1.7}$ au, subject to the close/wide degeneracy. This is the third $q < 10^{-4}$ planet from a high-magnification planetary signal ($A \gtrsim 65$). Together with another such planet, KMT-2021-BLG-0171Lb, the ongoing follow-up program for the KMTNet high-magnification events has demonstrated its ability in detecting high-magnification planetary signals for $q < 10^{-4}$ planets, which are challenging for the current microlensing surveys.

We present electron densities $n_{\rm e}$ in the inter-stellar medium (ISM) of star-forming galaxies at $z=4-9$ observed by the JWST/NIRSpec GLASS, ERO, and CEERS programs. We carefully evaluate line-spread functions of the NIRSpec instrument as a function of wavelength with the calibration data of a planetary nebula taken onboard, and obtain secure [OII]$\lambda\lambda$3726,3729 doublet fluxes for 14 galaxies at $z=4.02-8.68$ falling on the star-formation main sequence with the NIRSpec high and medium resolution spectra. We thus derive the electron densities of singly-ionized oxygen nebulae with the standard $n_{\rm e}$ indicator of [OII] doublet, and find that the electron densities of the $z=4-9$ galaxies are $n_{\rm e}\gtrsim 300$ cm$^{-3}$ significantly higher than those of low-$z$ galaxies at a given stellar mass, star-formation rate (SFR), and specific SFR. Interestingly, typical electron densities of singly ionized nebulae increase from $z=0$ to $z=1-3$ and $z=4-9$, which is approximated by the evolutionary relation of $n_{\rm e}\propto(1+z)^{p}$ with $p\sim 1-2$. Although it is not obvious that the ISM property of $n_{\rm e}$ is influenced by global galaxy properties, these results may suggest that nebula densities of high-$z$ galaxies are generally high due to the compact morphologies of high-$z$ galaxies evolving by $r_{\rm e}$ approximately proportional to $(1+z)^{-1}$ ($r_{\rm vir} \propto (1+z)^{-1}$) for a given stellar (halo) mass whose inverse square corresponds to the $p\sim 2$ evolutionary relation. The $p\sim 1-2$ evolutionary relation can be explained by a combination of the compact morphology and the reduction of $n_{\rm e}$ due to the high electron temperature of the high-$z$ metal poor nebulae.

The young massive-star-forming region M17 contains optically visible massive pre-main-sequence stars that are surrounded by circumstellar disks. Such disks are expected to disappear when these stars reach the main sequence. The physical and dynamical structure of these remnant disks are poorly constrained, especially the inner regions where accretion, photo-evaporation, and companion formation and migration may be ongoing. We aim to constrain the physical properties of the inner parts of the circumstellar disks of massive young stellar objects B243 (6 Msun) and B331 (12 Msun), two systems for which the central star has been detected and characterized previously despite strong dust extinction. Two-dimensional radiation thermo-chemical modelling with ProDiMo of double-peaked hydrogen lines of the Paschen and Brackett series observed with X-shooter was used to probe the properties of the inner disks. Additionally, the dust structure was studied by fitting the optical and near-infrared spectral energy distribution. B243 features a hot gaseous inner disk with dust at the sublimation radius at 3 AU. The disk appears truncated at roughly 6.5 AU; a cool outer disk of gas and dust may be present, but it cannot be detected with our data. B331 also has a hot gaseous inner disk. A gap separates the inner disk from a colder dusty outer disk starting at up to 100 AU. In both sources the inner disk extends to almost the stellar surface. Chemistry is essential for the ionization of hydrogen in these disks. The lack of a gap between the central objects and these disks suggests that they accrete through boundary-layer accretion. This would exclude the stars having a strong magnetic field. Their structures suggest that both disks are transitional in nature, that is to say they are in the process of being cleared, either through boundary-layer accretion, photo-evaporation, or through companion activity.

We exploit the increased sensitivity of the recently installed AO SOUL at the LBT to obtain new high-spatial-resolution NIR images of the massive young stellar object IRAS20126+4104 and its outflow. We aim to derive the jet proper motions and kinematics, as well as to study its photometric variability by combining the novel performances of SOUL together with previous NIR images. We used both broad-band ($K_{s}$, $K'$) and narrow-band (Br$\gamma$, H2) observations from a number of NIR cameras (UKIRT/UFTI,SUBARU/CIAO,TNG/NICS,LBT/PISCES,and LBT/LUCI1) to derive maps of the continuum and the H$_2$ emission in the 2.12 $\mu$m line. Three sets of images, obtained with AO systems (CIAO,2003; FLAO,2012; SOUL,2020), allowed us to derive the proper motions of a large number of H$_2$ knots along the jet. Photometry from all images was used to study the jet variability. We derived knot proper motions in the range of 1.7-20.3 mas yr$^{-1}$ (i.e. 13-158 km s$^{-1}$ at 1.64 kpc, avg. outflow tangential velocity $\sim$ 80 km s$^{-1}$). The derived knot dynamical age spans a $\sim$ 200-4000 yr interval. A ring-like H$_2$ feature near the protostar location exhibits peculiar kinematics and may represent the outcome of a wide-angle wind impinging on the outflow cavity. Both H$_2$ geometry and velocities agree with those inferred from proper motions of the H$_2$O masers, located at a smaller distance from the protostar. Although the total H$_2$ line emission from the knots does not exhibit time variations at a $\widetilde{>}$ 0.3 mag level, we have found a clear continuum flux variation (radiation scattered by the dust in the cavity opened by the jet) which is anti-correlated between the blue-shifted and red-shifted lobes and may be periodic (with a period of $\sim$ 12-18 yr). We suggest that the continuum variability might be related to inner-disc oscillations which have also caused the jet precession.

Context. The existence of meteor clusters has long since been a subject of speculation and so far only seven events have been reported, among which two involve less than five meteors, and three were seen during the Leonid storms. Aims. The 1995 outburst of Comet 73P/Schwassmann-Wachmann was predicted to result in a meteor shower in May 2022. We detected the shower, proved this to be the result of this outburst, and detected another meteor cluster during the same observation mission. Methods. The {\tau}-Herculids meteor shower outburst on 31 May 2022 was continuously monitored for 4 hours during an airborne campaign. The video data were analyzed using a recently developed computer-vision processing chain for meteor real-time detection. Results. We report and characterize the detection of a meteor cluster involving 38 fragments, detected at 06:48 UT for a total duration of 11.3 s. The derived cumulative size frequency distribution index is relatively shallow: s = 3.1. Our open-source computer-vision processing chain (named FMDT) detects 100% of the meteors that a human eye is able to detect in the video. Classical automated motion detection assuming a static camera was not suitable for the stabilized camera setup because of residual motion. Conclusions. From all reported meteor clusters, we crudely estimate their occurrence to be less than one per million observed meteors. Low heliocentric distance enhances the probability of such meteoroid self-disruption in the interplanetary space.

Since 2014, the University Campus of the Hellenic Open University (HOU) hosts the Astroneu array which is dedicated to the detection of Extensive Air Showers (EAS) induced by high energy Cosmic Rays (CR). The Astroneu array incorporates 9 large particle scintillation detectors and 6 antennas sensitive in the Radio Frequency (RF) range 1-200 MHz. The detectors are adjusted in three autonomous stations operating in an environment with strong electromagnetic background. As shown by previous studies, EAS radio detection in such environments is possible using innovative noise rejection methods, as well as advanced analysis techniques. In this work, we present the analysis of the collected radio data corresponding to an operational period of approximately four years. We present the performance of the Astroneu radio array in reconstructing the EAS axis direction using different RF detector geometrical layouts and a technique for the estimation of the shower core by comparing simulation and experimental data. Moreover, we measure the relative amplitudes of the two mechanisms that give rise to RF emission (Askaryan effect and Geomagnetic emission) and show that they are in good agreement with previous studies as well as with the simulation predictions.

One of the main scientific goals of the TESS mission is the discovery of transiting small planets around the closest and brightest stars in the sky. Here, using data from the CARMENES, MAROON-X, and HIRES spectrographs, together with TESS, we report the discovery and mass determination of a planetary system around the M1.5 V star GJ 806 (TOI-4481). GJ 806 is a bright (V=10.8 mag, J=7.3 mag) and nearby (d=12 pc) M dwarf that hosts at least two planets. The innermost planet, GJ 806 b, is transiting and has an ultra-short orbital period of 0.93 d, a radius of 1.331+-0.023 Re, a mass of 1.90+-0.17 Me, a mean density of 4.40+-0.45 g/cm3, and an equilibrium temperature of 940+-10 K. We detect a second, non-transiting, super-Earth planet in the system, GJ 806c, with an orbital period of 6.6 d, a minimum mass of 5.80+-0.30 Me, and an equilibrium temperature of 490+-5 K. The radial velocity data also shows evidence for a third periodicity at 13.6 d, although the current dataset does not provide sufficient evidence to unambiguously distinguish between a third super-Earth mass (Msin(i)=8.50+-0.45 Me) planet or stellar activity. Additionally, we report one transit observation of GJ 806 b taken with CARMENES in search for a possible extended atmosphere of H or He, but we can only place upper limits to its existence. This is not surprising as our evolutionary models support the idea that any possible primordial H/He atmosphere that GJ 806 b might have had, would long have been lost. However, GJ 806b's bulk density makes it likely that the planet hosts some type of volatile atmosphere. In fact, with a transmission spectroscopy metrics (TSM) of 44 and an emission spectroscopy metrics (ESM) of 24, GJ 806 b the third-ranked terrestrial planet around an M dwarf suitable for transmission spectroscopy studies, and the most promising terrestrial planet for emission spectroscopy studies.

The edge-on debris disk of the nearby young star Beta Pictoris shows an unusual brightness asymmetry in the form of a clump. The clump has been detected in both the mid-IR and CO and its origin has so far remained uncertain. Here we present new mid-IR observations of Beta Pic to track any motion of the dust clump. Together with previous observations, the data span a period of 12 years. We measured any projected displacement of the dust clump over the 12-yr period to be $0.2^{+1.3}_{-1.4}$ au away from the star based on the median and 1$\sigma$ uncertainty, and constrain this displacement to be <11 au at the 3$\sigma$ level. This implies that the observed motion is incompatible with Keplerian motion at the 2.8$\sigma$ level. It has been posited that a planet migrating outwards may trap planetesimals into a 2:1 resonance, resulting in the observed clump at pericentre of their orbits that trails the planet. The observed motion is also incompatible with such resonant motion at the 2.6$\sigma$ level. While Keplerian motion and resonant motion is still possible, the data suggest that the dust clump is more likely stationary. Such a stationary dust clump could originate from the collision or tidal disruption of a planet-sized body, or from secular perturbations due to a planet that create regions with enhanced densities in the disk.

Microwave Kinetic Inductance Detectors (MKIDs) are cryogenic photon detectors and are attractive because they permit simultaneous time, energy and spatial resolution of faint astronomical sources. We present a cost-effective alternative to dedicated (e.g. analogue) electronics for prototyping readout of single-pixel Optical/NIR MKIDs by repurposing existing and well-known ROACH-1 boards. We also present a pipeline that modernises previously-developed software and data frameworks to allow for extensiblity to new applications and portability to new hardware (e.g. Xilinx ZCU111 or 2x2 RFSoC boards).

Impacts by small solar system bodies (meteoroids, asteroids, comets and transitional objects) are characterized by a combination of energy dynamics and chemical modification on both terrestrial and small solar system bodies. In this context, the discovery of glycine amino acid in meteorites and comets has led to a hypothesis that impacts by astronomical bodies could contribute to delivery and polymerization of amino acids in the early Earth to generate proteins as essential molecules for life. Besides the possibility of abiotic polymerization of glycine, its decomposition by impacts could generate reactive groups to form other essential organic biomolecules. In this study, the high-pressure torsion (HPT) method, as a new platform for simulation of impacts by small solar system bodies, was applied to glycine. In comparison with high-pressure shock experiments, the HPT method simultaneously introduces high pressure and deformation strain. It was found that glycine was not polymerized in the experimental condition assayed, but partially decomposed to ethanol under pressures of 1 and 6 GPa and shear strains of <120 m/m. The detection of ethanol implies the inherent availability of remaining nitrogen-containing groups, which can incorporate to the formation of other organic molecules at the impact site. In addition, this finding highlights a possibility of the origin of ethanol previously detected in comets.

At least one dimensionless physical constant (i.e., a physically observable) must change for the cosmic time to make the varying speed of light (VSL) models phenomenologically feasible. Various physical constants and quantities also should be functions of cosmic time to satisfy all known local laws of physics, including special relativity, thermodynamics, and electromagnetism. Adiabaticity is another necessary condition to keep the homogeneity and isotropy of three-dimensional space [1]. To be a self-consistent theory, one should consider cosmic evolutions of physical constants and quantities when one derives Einstein's field equations and their solutions. All these conditions are well satisfied in the so-called minimally extended varying speed of light (meVSL) model [2]. Unlike other VSL models, we show that the redshift-drift formula of the meVSL model is the same as standard model one. Therefore, we cannot use this as an experimental tool to verify the meVSL. Instead, one can still use the cosmological chronometers (CC) as a model-independent test of the meVSL. The current CC data cannot distinguish meVSL from the standard model (SM) when we adopt the best-fit present values of the Hubble parameter and matter density contrast from the Planck mission. However, the CC data prefer the meVSL when we choose Pantheon's best-fit values of them.

There has recently been growing evidence that broad-line regions (BLRs) in active galactic nuclei have regular substructures, such as spiral arms, which are supported by the fact that the radii of BLRs measured by RM observations are generally consistent with the self-gravitating regions of accretion disks. We have shown in Paper I that the spiral arms excited by the gravitational instabilities in these regions may exist in some disk-like BLRs. As the second paper of the series, we investigate the loosely wound spiral arms excited by gravitational instabilities in disk-like BLRs and present their observational characteristics. Following the treatments of Adams et al. (1989), we solve the governing integro-differential equation by a matrix scheme. The emission-line profiles, velocity-delay maps, and velocity-resolved lags of the BLR spiral arms are calculated. We find that the spiral arms can explain some phenomena in observations: (1) the emission-line profiles in the mean and rms spectra have different asymmetries, (2) some velocity-delay maps, e.g., NGC 5548, have complex sub-features (incomplete ellipse), (3) the timescales of the asymmetry changes in emission-line profiles (rms spectra) are short. These features are attractive for modeling the observed line profiles and the properties of reverberation, and for revealing the details of the BLR geometry and kinematics.

We have recently suggested that the combination of the scalar virial theorem ($M_s \sim R_e \sigma^2$) and the $L=L'_0 \sigma^\beta(t)$ law, with L'_0 and $\beta$ changing from galaxy to galaxy (and with time), can provide a new set of equations valid for investigating the evolution of early-type galaxies (ETGs) (Donofrio & Chiosi, 2022). These equations are able to account for the tilt of the Fundamental Plane (FP) and to explain the observed distributions of ETGs in all its projections. In this paper we analyze the advantages offered by those equations, derive the $\beta$ and $L'_0$ parameters for real and simulated galaxies, and demonstrate that, according to the value of $\beta$, galaxies can move only along some permitted directions in the FP projections. Then, we show that simple galaxy models that grow in mass by infall of gas and form stars with a star formation rate depending on the stellar velocity dispersion nicely reproduce the observed distributions of ETGs in the FP projections and yield $\beta$s that agree with the measured ones. We derive the mutual relationships among the stellar mass, effective radius, velocity dispersion, and luminosity of ETGs as a function of $\beta$ and calculate the coefficients of the FP. Then, using the simple infall models, we show that the star formation history of ETGs is compatible with the $\sigma$-dependent star formation rate, and that both positive and negative values of $\beta$ are possible in a standard theory of galaxy evolution. The parameter $\beta(t)$ offers a new view of the evolution of ETGs. In brief, i) it gives a coherent interpretation of the FP and of the motions of galaxies in its projections; ii) it is the fingerprint of their evolution; iii) it measures the degree of virialization of ETGs; iv) and finally it allows us to infer their evolution in the near past.

There is a class of binary post-asymptotic giant branch (post-AGB) stars that exhibit remarkable near-infrared (NIR) excess. These stars are surrounded by disks with Keplerian or quasi-Keplerian dynamics and outflows composed of gas escaping from the rotating disk. Depending on the dominance of these components, there are two subclasses of binary post-AGB stars: disk-dominated and outflow-dominated. We aim to properly study the hourglass-like structure that surrounds the Keplerian disk around 89 Her. We present total-power on-the-fly maps of $^{12}$CO and $^{13}$CO $J$=2-1 emission lines in 89 Her. Previous studies are known to suffer from flux losses in the most extended components. We merge these total-power maps with previous NOEMA maps. The resulting combined maps are expected to detect the whole nebula extent of the source. Our new combined maps contain the entirety of the detectable flux of the source and at the same time are of high spatial resolution thanks to the interferometric observations. We find that the hourglass-like extended outflow around the rotating disk is larger and more massive than suggested by previous works. The total nebular mass of this very extended nebula is 1.8E-2 solar masses, of which 65% comes from the outflow. The observational data and model results lead us to classify the envelope around 89 Her as an outflow-dominated nebula, together with R Sct and IRAS 19125+0343 (and very probably AI CMi, IRAS 20056+1834, and IRAS 18123+0511). The updated statistics on the masses of the two post-AGB main components reveal that there are two distinct subclasses of nebulae around binary post-AGB stars depending on which component is the dominant one. We speculate that the absence of an intermediate subclass of sources is due to the different initial conditions of the stellar system and not because both subclasses are in different stages of the post-AGB evolution.

We have observed the three low-mass protostars, IRAS 15398$-$3359, L1527 IRS and TMC-1A, with the ALMA 12-m array, the ACA 7-m array, and the IRAM-30m and APEX telescopes in the C$^{18}$O $J=2$-1 emission. Overall, the C$^{18}$O emission shows clear velocity gradients at radii of $\sim$100-1000 au, which likely originate from rotation of envelopes, while velocity gradients are less clear and velocity structures are more perturbed on scales of $\sim$1000-10,000 au. IRAS 15398$-$3359 and L1527 IRS show a break at radii of $\sim$1200 and $\sim$1700 au in the radial profile of the peak velocity, respectively. The peak velocity is proportional to $r^{-1.38}$ or $r^{-1.7}$ within the break radius, which can be interpreted as indicating a rotational motion of the envelope with a degree of contamination of gas motions on larger spatial scales. The peak velocity follows $v_\mathrm{peak} \propto r^{0.68}$ or $v_\mathrm{peak} \propto r^{0.46}$ outside the break radius, which is similar to the $J/M$-$R$ relation of dense cores. TMC-1A exhibits the radial profile of the peak velocity not consistent with the rotational motion of the envelope nor the $J/M$-$R$ relation. The origin of the relation of $v_\mathrm{peak} \propto r^{0.46\operatorname{-}0.68}$ is investigated by examining correlations of the velocity deviation ($\delta v$) and the spatial scale ($\tau$) in the two sources. Obtained spatial correlations, $\delta v \propto \tau^{\sim0.6}$, are consistent with the scaling law predicted by turbulence models, which may suggest the large-scale velocity structures originate from turbulence.

The Search for Extraterrestrial Intelligence (SETI) has traditionally been conducted at radio wavelengths, but optical searches are well-motivated and increasingly feasible due to the growing availability of high-resolution spectroscopy. We present a data analysis pipeline to search Automated Planet Finder (APF) spectroscopic observations from the Levy Spectrometer for intense, persistent, narrow bandwidth optical lasers. We describe the processing of the spectra, the laser search algorithm, and the results of our laser search on 1983 spectra of 388 stars as part of the Breakthrough Listen search for technosignatures. We utilize an empirical spectra-matching algorithm called SpecMatch-Emp to produce residuals between each target spectrum and a set of best-matching catalog spectra, which provides the basis for a more sensitive search than previously possible. We verify that SpecMatch-Emp performs well on APF-Levy spectra by calibrating the stellar properties derived by the algorithm against the SpecMatch-Emp library and against Gaia catalog values. We leverage our unique observing strategy, which produces multiple spectra of each target per night of observing, to increase our detection sensitivity by programmatically rejecting events which do not persist between observations. With our laser search algorithm we achieve a sensitivity equivalent to the ability to detect an 84 kW laser at the median distance of a star in our dataset (78.5 ly). We present the methodology and vetting of our laser search, finding no convincing candidates consistent with potential laser emission in our target sample.

We propose a direct imaging method for the detection of exoplanets based on a combined low-rank plus structured sparse model. For this task, we develop a dictionary of possible effective circular trajectories a planet can take during the observation time, elements of which can be efficiently computed using rotation and convolution operation. We design a simple alternating iterative hard-thresholding algorithm that jointly promotes a low-rank background and a sparse exoplanet foreground, to solve the non-convex optimisation problem. The experimental comparison on the $\beta$-Pictoris exoplanet benchmark dataset shows that our method has the potential to outperform the widely used Annular PCA for specific planet light intensities in terms of the Receiver operating characteristic (ROC) curves.

We present a dynamical analysis of a quasar-host galaxy at $z\simeq 6$ (SDSS J2310+1855) using a high-resolution ALMA observation of the [CII] emission line. The observed rotation curve is fitted with mass models that consider the gravitational contribution of a thick gas disc, a thick star-forming stellar disc, and a central mass concentration, likely due to a combination of a spheroidal component (i.e. a stellar bulge) and a supermassive black hole (SMBH). The SMBH mass of $5\times 10^9\ \rm M_{\odot}$, previously measured using the CIV and MgII emission lines, is not sufficient to explain the high velocities in the central regions. Our dynamical model suggests the presence of a stellar bulge with a mass of $\rm M_{bulge}\sim 10^{10}\ \rm M_{\odot}$ in this object, when the Universe was less than 1 Gyr old. To end up on the local $M_{\rm SMBH}-M_{\rm bulge}$ relation, the bulge mass should increase by a factor of $\sim$40 from $z=6$ to 0, while the SMBH mass should grow at most by a factor of 4, pointing to asynchronous galaxy-BH co-evolution. Imaging with JWST will allow us to validate this scenario.

Most structural and evolutionary properties of galaxies strongly rely on the stellar initial mass function (IMF), namely the distribution of the stellar mass formed in each episode of star formation. As the IMF shapes the stellar population in all stellar systems, it turns out to become one of the most fundamental concepts of modern astronomy. Both constant and variable IMFs across different environments have been claimed despite a large number of theoretical and observational efforts. However, the measurement of the IMF in Galactic stellar populations has been limited by the relatively small number of photometrically observed stars, leading to high uncertainties. Here we report a star-counting result based on ~93,000 spectroscopically observed M-dwarf stars, an order of magnitude more than previous studies, in the 100--300 parsec (pc) Solar neighbourhood. We find unambiguous evidence of a variable IMF that depends on both metallicity and stellar age. Specifically, the stellar population formed at the early time contains fewer low-mass stars compared to the canonical IMF, independent of stellar metallicities. In present days, on the other hand, the proportion of low-mass stars increases with stellar metallicity. The variable abundance of low-mass stars in our Milky Way establishes a powerful benchmark for models of star formation and can heavily impact results in Galactic chemical enrichment modelling, mass estimation of galaxies, and planet formation efficiency.

We report the detection of large-amplitude, rapid radial velocity (RV) variations and line-splitting in high-resolution Keck/NIRSPEC spectra of the M9 dwarf LP 413-53. We attribute these features to binary motion. Analyzing data spanning 9 months, we infer an orbital period of 0.852725$^{+0.000002}_{-0.000003}$~day, an eccentricity of 0.080$^{+0.020}_{-0.013}$, a primary RV semi-amplitude of 24.2$^{+1.8}_{-1.4}$ km~s$^{-1}$, and a secondary RV semi-amplitude of 29.4$^{+2.2}_{-1.7}$ km~s$^{-1}$, implying a system mass ratio $M_\mathrm{secondary}$/$M_\mathrm{primary}$ = 0.822$^{+0.009}_{-0.008}$. These measurements identify LP 413-53 as the shortest-period ultracool binary discovered to date, and one of the smallest separation main sequence binaries known. The position and velocity of the system rules out previously reported membership in the Hyades Moving Group, and indicate that this is likely a pair of evolved (age $\gtrsim$ 1 Gyr), very-low-mass stars. Assuming masses consistent with evolved late-M and L dwarfs, we estimate an orbital separation of 0.0093-0.0095~au or 19-22 stellar radii, and an orbital inclination angle of 27$\pm$2 deg, making it unlikely that this system exhibits eclipse events. The larger radii of these stars at young ages would have put them in contact at the system's current separation, and we speculate that this system has undergone dynamical evolution, either through orbital angular momentum loss or ejection of a third component followed by tidal circularization.

Optical Very-Long-Baseline Interferometers (VLBI), widely used in astronomy, require phase-stable optical links across stations, which impose a limit on baseline distances, and, in turn, limits measurement precision. Here we describe a novel type of two-photon quantum-assisted interferometer, which may allow improvements in precision by orders of magnitude benefiting numerous fields in cosmology and astrophysics. We tested a tabletop version of the interferometer and unambiguously observe correlated behavior in detections of photon pairs from two thermal light sources, in agreement with theoretical predictions. This work opens new possibilities in astronomical measurements.

This paper studies both $\Lambda$CDM and CDM models under the \"uber gravity theory, named \"u$\Lambda$CDM and \"uCDM respectively. We report bounds over their parameter phase-space using several cosmological data, in particular, the recent Pantheon+ sample. Based on the joint analysis, the best fit value of the \"uber characteristic parameter is $z_\oplus = 0.028^{+0.036}_{-0.020}$ and $z_\oplus = 0.960^{+0.031}_{-0.030}$ at 68\% confidence level for \"u$\Lambda$CDM and \"uCDM respectively. Furthermore, we find that the $\mathbb{H}0(z)$ diagnostic suggests the $H_0$ tension is not alleviated. Finally, both models are statistically compared with $\Lambda$CDM through the Akaike and Bayesian information criteria, which suggest that there is a modest evidence against for \"u$\Lambda$CDM, a strongest evidence for \"uCDM, against for the joint analysis, but both \"uber gravity models and $\Lambda$CDM are equally preferred for most of the single samples.

We describe new JWST/NIRSpec observations of galaxies at $z\gtrsim7$ taken as part of the CEERS survey of the EGS field. Previous observations of this area have revealed associations of Ly$\alpha$ emitters at redshifts ($z=7.5$, $7.7$, $8.7$) where the IGM is thought mostly neutral, leading to suggestions that these systems are situated in large ionized bubbles. We identify 21 $z\gtrsim7$ galaxies with robust redshifts in the CEERS dataset, including 10 in the Ly$\alpha$ associations. Their spectra are indicative of very highly ionized and metal poor gas, with line ratios (O32 $=17.84$ and Ne3O2 $=0.89$) and metallicity ($12+\log{[\rm{O/H}]}=7.84$) that are rarely seen at lower redshifts. We find that the most extreme spectral properties are found in the six $z\gtrsim7$ Ly$\alpha$ emitting galaxies in the sample. Each have hard ionizing spectra indicating that their visibility is likely enhanced by efficient ionizing photon production. Ly$\alpha$ velocity offsets are found to be very large ($\gtrsim300$ km s$^{-1}$), likely also contributing to their detectability. We find that Ly$\alpha$ in $z\gtrsim7$ galaxies is $6-12\times$ weaker than in lower redshift samples with matched rest-frame optical spectral properties. If the bubbles around the Ly$\alpha$ emitters are relatively small ($\lesssim0.5-1$ pMpc), we may expect such significant attenuation of Ly$\alpha$ in these ionized regions. We discuss several other effects that may contribute to weaker Ly$\alpha$ emission at $z\gtrsim7$. Deep spectroscopy of fainter galaxies in the vicinity of the Ly$\alpha$ emitters will better characterize the physical scale of the ionized bubbles in this field.

We develop a numerical method for directly computing the dissipative dynamical tidal response of rapidly rotating, oblate stars and gaseous planets with realistic internal structures. Applying these calculations to neutrally and stably stratified polytropes, we identify the most relevant resonances in models with rotation rates up to nearly the mass-shedding limit. We then compute the dynamical tidal response for Jupiter interior models including both stably stratified and convective regions. These calculations show that resonances involving waves with both gravito-inertial and purely inertial character are capable of explaining a discrepancy between observations and hydrostatic calculations of Jupiter's response to tidal forcing by Io. This result contrasts with recent work that excluded Jupiter's rotational flattening, and opens the door to resonances involving a wider range of internal oscillation modes than previously considered.

A unified model independent effective description of cosmological perturbations is derived in terms of two effective quantities, playing the role of effective propagation speeds of curvature perturbations and gravitational waves, encoding the effects of the interaction of perturbations at any order, and inducing a modification of the friction term of the perturbations propagation equation. The approach can be applied to dark energy, modified gravity, dark matter, for fields of arbitrary number and spin, and is particularly suitable for model independent analysis of observational data. The structure of the effective actions and equations is the same for scalar and tensor perturbations. The effective actions can be written as the Klein-Gordon action in terms of an appropriately defined effective metric, dependent on the effective speed. In this geometrical interpretation, the effective metric emerges as the result of the interaction and self-interaction of perturbations, hinting to possible connections with emergent gravity. As an example of an application we find that the effective speeds of curvature perturbations and gravitational waves can be frequency and polarization dependent even for a minimally coupled scalar field in general relativity, when higher order terms effects are computed, going beyond the quadratic action, or in axion inflation. We discuss the relation between the effective speed and quantum correlators.

We study scalar field theories invariant under transverse diffeomorphisms in cosmological contexts. We show that in the geometric optics approximation, the corresponding particles move along geodesics and contribute with the same active mass (energy) to the gravitational field as in Diff invariant theories. However, for low-frequency modes, the contributions to the energy-momentum tensor differ from that of Diff invariant theories. This opens up a wide range of possibilities for cosmological model building. As an example, we show that the simplest TDiff invariant scalar field theory with only kinetic term could drive inflation and generate a red-tilted spectrum of density fluctuations. We also present a detailed analysis of cosmological perturbations and show that the breaking of full Diff invariance generically induces new non-adiabatic pressure perturbations. A simple scalar field dark matter model based on a purely kinetic term that exhibits the same clustering properties as standard cold dark matter is also presented.

The accurate sky localization of gravitational wave (GW) sources is an important scientific goal for space-based GW detectors. Due to the effects of gravity on three spacecrafts, it is hard to maintain the equality of the arm length, so the time-delay interferometry (TDI) method is needed to cancel out the laser frequency noise for space-based GW detectors. By considering the first-generation TDI combination, we employ the Fisher information matrix to study the accuracy of sky localizations for future space-based GW detectors and their combined network. The main difference between future space-based GW detectors includes the time-changing orientation of the detector plane, the arm length, the orbital period of spacecrafts and the noise curve. We study the effects of these factors on the accuracy of source localization at different frequencies. We find that the amplitude modulation caused by the rotation of the detector plane can help LISA and Taiji not only to improve the accuracy of source localization but also to enlarge the sky coverage at frequencies below 1 mHz. As the frequency of monochromatic GWs increases, the Doppler modulation becomes dominate and the equatorial pattern appears in the sky map. The effect of arm length on the angular resolution mainly comes from the noise curve and it is almost the same for both heliocentric and geocentric constellations. The orbital period of the spacecrafts has little effect on the angular resolutions. The improvement on the angular resolutions by the network of combined detectors is small compared with a single detector and the angular resolutions are almost the same with and without the TDI combination.

Suspended optics in gravitational wave (GW) observatories are susceptible to alignment perturbations and, in particular, to slow drifts over time due to variations in temperature and seismic levels. Such misalignments affect the coupling of the incident laser beam into the optical cavities, degrade both circulating power and optomechanical photon squeezing, and thus decrease the astrophysical sensitivity to merging binaries. Traditional alignment techniques involve differential wavefront sensing using multiple quadrant photodiodes, but are often restricted in bandwidth and are limited by the sensing noise. We present the first-ever successful implementation of neural network-based sensing and control at a gravitational wave observatory and demonstrate low-frequency control of the signal recycling mirror at the GEO 600 detector. Alignment information for three critical optics is simultaneously extracted from the interferometric dark port camera images via a CNN-LSTM network architecture and is then used for MIMO control using soft actor-critic-based deep reinforcement learning. Overall sensitivity improvement achieved using our scheme demonstrates deep learning's capabilities as a viable tool for real-time sensing and control for current and next-generation GW interferometers.

We consider Kerr black holes (BHs) surrounded by perfect dark fluid matter (PFDM), with an additional parameter ($k$) because of PFDM, apart from mass ($M$) and rotation parameter ($a$) -- the rotating PFDM BHs. We analyze the photon orbits around PFDM BHs and naked singularities (NSs) and emphasise the effect of PFDM on photon boomerangs. Interestingly, the azimuthal oscillations first increase and then decrease for retrograde orbits, whereas they first decrease and then increase for prograde orbits, with increasing $k$. Unlike in the Kerr NSs, photon boomerangs can form around rotating PFDM NSs. We use the Event Horizon Telescope (EHT) observational results for Schwarzschild shadow deviations of M87* and Sgr A*, $\delta_{M87^*}=-0.01\pm0.17$ and $\delta_{Sgr A^*} = -0.08^{+0.09}_{-0.09}~\text{(VLTI)},-0.04^{+0.09}_{-0.10}~\text{(Keck)}$, to report the upper bounds on the PFDM parameter: $0\leq k\leq 0.0792M$ and $k^{max}\in[0.0507M, 0.0611M]$ respectively. Together with the EHT bounds on the shadows of Sgr A$^*$ and M87$^*$, our analysis concludes that a substantial part of the rotating PFDM BH parameter space agrees with the EHT observations. Thus, one must consider the possibility of the rotating PFDM BHs being strong candidates for the astrophysical BHs.

A radiative decaying Big Bang relic with a mass $m_a\simeq 5-25 \,\rm eV$, which we dub "blue axion", can be probed with direct and indirect observations of the cosmic optical background (COB). The strongest bounds on blue-axion cold dark matter come from the Hubble Space Telescope (HST) measurements of COB anisotropies at $606$ nm. We suggest that new HST measurements at higher frequencies ($336$ nm and $438$ nm) can improve current constraints on the lifetime up to an order of magnitude, and we show that also thermally produced and hot relic blue axions can be competitively probed by COB anisotropies. We exclude the simple interpretation of the excess in the diffuse COB detected by Long Range Reconnaissance Imager (LORRI) as photons produced by a decaying hot relic. Finally, we comment on the reach of upcoming line intensity mapping experiments, that could detect blue axions with a lifetime as large as $10^{29}\,\rm s$ or $10^{27}\,\rm s$ for the cold dark matter and the hot relic case, respectively.

The stability of a class of electrically charged fluid spheres under radial perturbations is studied. Among these spheres there are regular stars, overcharged tension stars, regular black holes, quasiblack holes, and quasinonblack holes, all of which have a Reissner-Nordstr\"om exterior. We formulate the dynamical perturbed equations by following the Chandrasekhar approach and investigate the stability against radial perturbations through numerical methods. It is found that (i) under certain conditions that depend on the adiabatic index of the radial perturbation, there are stable charged stars and stable tension stars; (ii) also depending on the adiabatic index there are stable regular black holes; (iii) quasiblack hole configurations formed by, e.g., charging regular pressure stars or by discharging regular tension stars, can be stable against radial perturbations for reasonable values of the adiabatic index; (iv) quasinonblack holes are unstable against radial perturbations.

Light massive preheat fields acquire a non-vanishing dispersion during parametric resonance from their quantum particle production. This in turn will modify the inflaton potential, which in some cases can induce a transient period of acceleration. We illustrate this phenomenon in the setup of non-supersymmetric non-minimal M-flation (non-$\mathbb{M}$-flation) which has some motivations from the brane compactifications in string theory. Implementing a lattice simulation by the LATTICEEASY code, we compute the potential correction term in our scenario and show that the modified term indeed causes the universe to make a transition from the decelerated expansion to a temporary phase of acceleration. The correction term reduces to some extent the number density of the particles generated during preheating, but the efficiency of preheating remains still enough to have successful particle production after inflation. We also compute the spectrum of the gravitational waves (GWs) generated during preheating in our setup by using the LATTICEEASY code. Although the peak frequency remains almost the same, the inclusion of the correction term reduces the amplitude of the gravitational spectrum by almost one order of magnitude.

By the twelfth century, northern European scholars gradually embraced Arabic innovations in science and technology. England naturally developed into a significant centre of the new learning in western Europe. Hereford, and specifically its cathedral school, played a particularly important role in the transition of English scholarship to the new learning. Hereford cathedral developed into a focal point for high-level scholarship, attracting numerous scholars from across the continent. Roger of Hereford stands out among his peers as an enlightened scholar who made more practical use than most of the full astronomical and astrological knowledge base available in England at the time. A significant body of recent scholarship focuses on twelfth-century ecclesiastical developments, including those relating to Roger of Hereford's Computus. However, much less scholarly emphasis is placed on Roger's astronomical calculations, particularly those which allowed him to establish an important reference meridian at Hereford. Those aspects are the focus of this paper.

We demonstrate how to quantify the frequency-domain amplitude and phase accuracy of waveform models, $\delta A$ and $\delta \phi$, in a form that could be marginalized over in gravitational-wave inference using techniques currently applied for quantifying calibration uncertainty. For concreteness, waveform uncertainties affecting neutron-star inspiral measurements are considered, and post-hoc error estimates from a variety of waveform models are made by comparing time-domain and frequency-domain analytic models with multiple-resolution numerical simulations. These waveform uncertainty estimates can be compared to GW170817 calibration envelopes or to Advanced LIGO and Virgo calibration goals. Signal-specific calibration and waveform uncertainties are compared to statistical fluctuations in gravitational-wave observatories, giving frequency-dependent modeling requirements for detectors such as Advanced LIGO Plus, Cosmic Explorer, or Einstein Telescope. Finally, the distribution of waveform error for GW170817 over the parameters of the low-spin posterior is computed from tidal models and compared to the constraints on $\delta \phi$ or $\delta A$ from GWTC-1 by Edelman et. al. In general, $\delta \phi$ and $\delta A$ can also be interpreted in terms of unmodeled astrophysical energy transfer within or from the source system.

We formulate a Bayesian framework to analyze ringdown gravitational waves from colliding binary black holes and test the no-hair theorem. The idea hinges on mode cleaning -- revealing subdominant oscillation modes by removing dominant ones using newly proposed ${\it rational~filters}$. By incorporating the filter into Bayesian inference, we construct a likelihood function that depends only on the mass and spin of the remnant black hole (no dependence on mode amplitudes and phases) and implement an efficient pipeline to constrain the remnant mass and spin without Markov chain Monte Carlo (MCMC). We test ringdown models by cleaning combinations of different modes and evaluating the consistency between the residual data and pure noise. The model evidence and Bayes factor are used to demonstrate the presence of a particular mode and to infer the mode starting time. In addition, we design a hybrid approach to estimate the remnant black hole properties exclusively from a single mode using MCMC after mode cleaning. We apply the framework to GW150914 and demonstrate more definitive evidence of the first overtone by cleaning the fundamental mode. This new framework provides a powerful tool for black hole spectroscopy in future gravitational-wave events.

The Brazilian gravitational-wave detector Mario Schenberg was conceived in the early 2000s and operated until 2016 when it was dismantled. A straight path to evaluate the viability of the reassembly of the Schenberg antenna is to verify the possibility of detecting gravitational wave (GW) signals within its design sensitivity features. The eventual identification of significant signals would operate as motivation for the Schenberg rebuild. As the antenna was dismantled, we can get some indication from the third observing run (O3) data of the LIGO detectors. It is based on the similarity between Schenberg sensitivity and the sensitivity of the interferometers in the O3 [3150-3260] Hz band. We search for signals with milliseconds to a few seconds without making assumptions about their morphology, polarization, and arrival sky direction. The data were analyzed with the coherent WaveBurst pipeline (cWB) with frequencies from 512 Hz to 4096 Hz and the search targets only signals with bandwidth overlapping the Schenberg frequency band. No statistically significant evidence of GW bursts during O3 was found. The null result was used to feature the search efficiency in identifying different simulated signal morphologies and establish upper limits on the GW burst event rate as a function of its strain amplitude. The present search, and consequently Schenberg, is sensitive to sources emitting isotropically 5 x 10e(-6) M_sun c^2 in GWs from a distance of 10 kpc with 50% detection efficiency and with a false alarm rate of 1/100 years. The feasibility of detecting f-modes of neutron stars excited by glitches was also investigated. The Schenberg antenna would need at least 5.3 years of observation run to get a single detection of the f-mode signal, given E_(glitch) approx 10e(-10) M_sun c^2.

Both ultralight dark matter and exploring the quantum nature of black holes are all topics of great interest in gravitational wave astronomy at present. The superradiant instability allows an exotic compact object (ECO) to be surrounded by an ultralight boson cloud, which leads to the emission of gravitational waves and further triggers rich dynamical effects. In this paper, we study the gravitational effects of superradiant instabilities by calculating the energy fluxes of gravitational waves emitted from ultralight scalar dark matter fields by solving the Teukolsky equation in the background of a massive ECO phenomenologically described by a Kerr geometry with a reflective boundary condition at its physical boundary. We find that both the amplitude and phase of the reflectivity will either suppress or enhance the energy flux of GWs by several orders of magnitude if $M\mu \gtrsim 0.5$ where $M$ and $\mu$ are the mass of ECO and boson, respectively. However, the modifications to energy flux are negligible if $M \mu \lesssim 0.5$. Our results suggest that reflectivity will play a significant role in the near-horizon physics of ECO.

We report how to alleviate both the $H_0$ and $\sigma_8$ tensions simultaneously within $f(T)$ gravity. In particular, we consider the parametrization $f(T)=-T-2\Lambda/M_P^2+\alpha T^\beta$, where two out of the three parameters are independent. This model can efficiently fit observations solving the two tensions. To our knowledge, this is the first time where a modified gravity theory can alleviate both $H_0$ and $\sigma_8$ tensions simultaneously, hence, offering an additional argument in favor of gravitational modification.

We study the LISA sources that arise from isolated binary evolution, and how these depend on age and metallicity, using model stellar populations from BPASS. We calculate the combined GW spectrum of all the binaries within these populations, including all types of compact binaries as well as those with living stars, and use these results to investigate the detectability of star clusters with LISA. We find at late times the dominant sources are WD-WD binaries, but at times before approx. $10^9$ years we find a significant population of NS-WD and BH-WD binaries, which is related to the treatment of mass transfer and common envelope events in BPASS, wherein mass transfer is relatively likely to be stable. Metallicity also has an effect on the GW spectrum and on the relative dominance of different types of binaries. Our results suggest that nearby star clusters might produce GWs detectable by LISA or future missions throughout most of their evolution.

We present a model of a slowly rotating Tolman VII (T-VII) fluid sphere, at second order in the angular velocity. The structure of this configuration is obtained by integrating the Hartle-Thorne equations for slowly rotating relativistic masses. We model a sequence in adiabatic and quasi-stationary contraction, by varying the tenuity parameter $R/R_{\mathrm{S}}$, where $R$ is the radius of the configuration and $R_{\mathrm{S}}$ is its Schwarzschild radius. We determined the moment of inertia $I$, mass quadrupole moment $Q$, and the ellipticity $\varepsilon$, for various configurations. Similar to previous results for Maclaurin and polytropic spheroids, in slow rotation, we found a change in the behaviour of the ellipticity when the tenuity reaches a certain critical value. We compared our results of $I$ and $Q$ for the T-VII model with those predicted by the universal fittings proposed for realistic neutron stars. For the relevant range of compactness, we found that relative errors are within $10\%$, thus suggesting the T-VII solution as a very good approximation for the description of the interior of neutron stars.

We derive exact and closed-form expressions for a large class of two-point and three-point inflation correlators with the tree-level exchange of a single massive particle. The intermediate massive particle is allowed to have arbitrary mass, spin, chemical potential, and arbitrary nonderivative or derivative couplings to external inflaton modes. We also allow the coupling coefficients to have arbitrary complex power dependences on the conformal time. Our results feature closed-form expressions involving only familiar special functions and without any infinite sums. This is achieved by an improved bootstrap method with a suitable change of variables. Our results cover a wide range of cosmological collider models and can be directly used for future phenomenological studies. Our results can also be used as basic building blocks for constructing more complicated inflation correlators.

Dense neutrino gases form in extreme astrophysical sites, and the flavor content of the neutrinos likely has an important impact on the subsequent dynamical evolution of their environment. Through coherent forward scattering among neutrinos, the flavor content of the gas evolves under a time-dependent potential which can be modeled in a quantum many-body formalism as an all-to-all coupled spin-spin interaction. This two-body potential generically introduces entanglement and greatly complicates the study of these systems. In this work we study the evolution of the quantum many-body problem as well as the typically employed mean-field approximation to it for a small number of neutrinos ($N = 16$). We consider randomly chosen one- and two-body couplings in the Hamiltonian, and the resulting evolution of several initial product states. We subsequently compare many-body and mean-field predictions for one-body observables, and we consider one- and two-body entanglement to assess under what conditions the many-body and mean-field predictions are likely to disagree. Except for a special category of prototypical initial conditions, we find that the typically employed mean-field approximation is insufficient to capture the evolution of one-body operators in the systems we consider. We also observe a loss of coherence in one- and two-body trace-reduced subsystems which suggests that the evolution may be well approximated as a classical mixture of separable states.

Hypothetical particles called heavy neutral leptons (HNLs) can be produced in large quantities in the cores of supernovae during the first seconds of the explosion. These particles then decay, producing secondary energetic neutrinos that can be detected by neutrino detectors. In this paper, I identify a region of the HNL parameter space that could be tested using this method, assuming a supernova explosion at distances from 0.2 to 10 kpc. The range of HNLs masses $m_N \sim 160-700$ MeV and lifetimes of $\tau_N \gtrsim 0.02$ seconds can be probed using the Hyper-Kamiokande neutrino detector. This region of the parameter space is complementary to existing bounds from primordial nucleosynthesis and to the expected sensitivity of the future SHiP experiment, thus covering a gap in our current knowledge of HNLs up to masses of $m_N \simeq 400$ MeV.