Papers for Tuesday, Feb 28 2023

Radiation belts are present in all large-scale Solar System planetary magnetospheres: Earth, Jupiter, Saturn, Uranus, and Neptune. These persistent equatorial zones of trapped high energy particles up to tens of MeV can produce bright radio emission and impact the surface chemistry of close-in moons. Recent observations confirm planet-like radio emission such as aurorae from large-scale magnetospheric current systems on very low mass stars and brown dwarfs. These objects, collectively known as ultracool dwarfs, also exhibit quiescent radio emission hypothesized to trace stellar coronal flare activity or extrasolar radiation belt analogs. Here we present high resolution imaging of the ultracool dwarf LSR J1835+3259 demonstrating that this radio emission is spatially resolved and traces a long-lived, double-lobed, and axisymmetric structure similar in morphology to the Jovian radiation belts. Up to 18 ultracool dwarf radii separate the two lobes. This structure is stably present in three observations spanning >1 year. We infer a belt-like distribution of plasma confined by the magnetic dipole of LSR J1835+3259, and we estimate ~15 MeV electron energies that are consistent with those measured in the Jovian radiation belts. Though more precise constraints require higher frequency observations, a unified picture where radio emissions in ultracool dwarfs manifest from planet-like magnetospheric phenomena has emerged.

Galactic nuclei, the densest stellar environments in the Universe, exhibit a complex geometrical structure. The stars orbiting the central supermassive black hole follow a mass segregated distribution both in the radial distance from the center and in the inclination angle of the orbital planes. The latter distribution may represent the equilibrium state of vector resonant relaxation (VRR). In this paper, we build simple models to understand the equilibrium distribution found previously in numerical simulations. Using the method of maximising the total entropy and the quadrupole mean-field approximation, we determine the equilibrium distribution of axisymmetric two-component gravitating systems with two distinct masses, semimajor axes, and eccentricities. We also examine the limiting case when one of the components dominates over the total energy and angular momentum, approximately acting as a heat bath, which may represent the surrounding astrophysical environment such as the tidal perturbation from the galaxy, a massive perturber, a gas torus, or a nearby stellar system. Remarkably, the bodies above a critical mass in the subdominant component condense into a disk in a ubiquitous way. We identify the system parameters where the transition is smooth and where it is discontinuous. The latter cases exhibit a phase transition between an ordered disk-like state and a disordered nearly spherical distribution both in the canonical and in the microcanonical ensembles for these long-range interacting systems.

The nuclear transient AT2019cuk/Tick Tock/SDSS J1430+2303 has been suggested to harbor a supermassive black hole (SMBH) binary near coalescence. We report results from high-cadence NICER X-ray monitoring with multiple visits per day from January-August 2022, as well as continued optical monitoring during the same time period. We find no evidence of periodic/quasi-periodic modulation in the X-ray, UV, or optical bands, however we do observe exotic hard X-ray variability that is unusual for a typical AGN. The most striking feature of the NICER light curve is repetitive hard (2-4 keV) X-ray flares that result in distinctly harder X-ray spectra compared to the non-flaring data. In its non-flaring state, AT2019cuk looks like a relatively standard AGN, but it presents the first case of day-long, hard X-ray flares in a changing-look AGN. We consider a few different models for the driving mechanism of these hard X-ray flares, including: (1) corona/jet variability driven by increased magnetic activity, (2) variable obscuration, and (3) self-lensing from the potential secondary SMBH. We prefer the variable corona model, as the obscuration model requires rather contrived timescales and the self-lensing model is difficult to reconcile with a lack of clear periodicity in the flares. These findings illustrate how important high-cadence X-ray monitoring is to our understanding of the rapid variability of the X-ray corona and necessitate further high-cadence, multi-wavelength monitoring of changing-look AGN like AT2019cuk to probe the corona-jet connection.

We present UV and Ly$\alpha$ radial surface brightness (SB) profiles of Ly$\alpha$ emitters (LAEs) at $z=2.84$ detected with the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The depth of our data, together with the wide field coverage including a protocluster, enable us to study the dependence of Ly$\alpha$ halos (LAHs) on various galaxy properties, including Mpc-scale environments. UV and Ly$\alpha$ images of 3490 LAEs are extracted, and stacking the images yields SB sensitivity of $\sim1\times10^{-20}\mathrm{~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$ in Ly$\alpha$, reaching the expected level of optically thick gas illuminated by the UV background at $z\sim3$. Fitting of the two-component exponential function gives the scale-lengths of $1.56\pm0.01$ and $10.4\pm0.3$ pkpc. Dividing the sample according to their photometric properties, we find that while the dependence of halo scale-length on environment outside of the protocluster core is not clear, LAEs in the central regions of protoclusters appear to have very large LAHs which could be caused by combined effects of source overlapping and diffuse Ly$\alpha$ emission from cool intergalactic gas permeating the forming protocluster core irradiated by active members. For the first time, we identify ``UV halos'' around bright LAEs which are probably due to a few lower-mass satellite galaxies. Through comparison with recent numerical simulations, we conclude that, while scattered Ly$\alpha$ photons from the host galaxies are dominant, star formation in satellites evidently contributes to LAHs, and that fluorescent Ly$\alpha$ emission may be boosted within protocluster cores at cosmic noon and/or near bright QSOs.

In this paper, we compute and analyze synthetic radio images of gamma-ray bursts and kilonova afterglows. For modeling the former, we consider GRB170817A-inspired set of parameters, while for the latter, we employ ejecta profiles from numerical-relativity simulations. We find that the kilonova afterglow sky map has a doughnut-like structure at early times that becomes more ring-like at late times. This is caused by the fact that the synchrotron emission from electrons following Maxwellian distribution function dominates the early, beamed, emission while emissions from electrons following power-law distribution is important at late times. For an on-axis observer, the image flux centroid moves on the image plane initially away from the observer. Flux centroid displacement The image sizes, we find, are the largest for equal mass merger simulations with the soft equation of state. The presence of a kilonova afterglow affects the properties inferred from the source sky map even if the gamma-ray burst afterglow dominates the total flux density. The main effect is the reduction of the mean apparent velocity of the source, and an increase in the source size. Thus, neglecting the presence of the kilonova afterglow may lead to systematic errors in the inference of gamma-ray burst properties from the sky map observations. Notably, at the observing angle inferred for GRB170817A the presence of kilonova afterglow would affect the sky map properties only at very late times $t\gtrsim1500\,$days.

The Pioneer Venus Orbiter (PVO) Gamma-ray burst experiment detected 318 gamma-ray bursts over about 14 years between 1978 and 1992 with near $4\pi$ coverage. This data set complements BATSE by determining the properties of the brightest gamma-ray bursts. PVO places a constrains on the slope of the bright end of the Log N-Log P distribution. The slope is -1.52$\pm 0.15$.

For the first time, the sulphur abundance relative to hydrogen (S/H) in the Narrow Line Regions of a sample of Seyfert 2 nuclei (Sy 2s) has been derived via direct estimation of the electron temperature. Narrow emission line intensities from the SDSS DR17 [in the wavelength range 3000 < $\lambda$ < 9100] and from the literature for a sample of 45 nearby ($z$ < 0.08) Sy 2s were considered. Our direct estimates indicate that Sy 2s have similar temperatures in the gas region where most of the S+ ions are located in comparison with that of star-forming regions (SFs). However, Sy 2s present higher temperature values ($\sim$10000 K) in the region where most of the S++ ions are located relative to that of SFs. We derive the total sulphur abundance in the range of 6.2 < 12 + log(S/H) < 7.5, corresponding to 0.1-1.8 times the solar value. These sulphur abundance values are lower by $\sim$0.4 dex than those derived in SFs with similar metallicity, indicating a distinct chemical enrichment of the ISM for these object classes. The S/O values for our Sy 2 sample present an abrupt ($\sim$0.5 dex) decrease with increasing O/H for the high metallicity regime [12 + log(O/H) > 8.7)], what is not seen for the SFs. However, when our Sy 2 estimates are combined with those from a large sample of star-forming regions, we did not find any dependence between S/O and O/H.

We report a detailed study of features of electron-impact excitation (EIE) of Ca IV for the first time using the relativistic Breit-Pauli R-Matrix method with a large close coupling wavefunction expansion of 54 fine structure levels belonging to n=2,3,4 complexes. Our study predicts presence of a strong 3.2 $\mu$m emission line in IR. The EIE collision strength ($\Omega$) shows extensive resonances with enhanced background resulting in an effective collision strength ($\gamma$) of 2.2 at about 10,000 K that increases to 9.66 around 300,000 K. The present results include collision strength of all 1431 excitations among the 54 levels and effective collision strength for a limited number of transitions of possible interest. We have found extensive resonances in the low energy region, convergence of the resonances, and of the partial waves with the 54 levels wavefunction. At higher energy, the collision strength decreases beyond the resonance region for forbidden transitions, is almost constant or decreases slowly for dipole-allowed transitions with low oscillator strengths, and rises with Bethe-Coulomb behavior of ln(E)to almost a plateau for transitions with high f-values.

The AGN STORM 2 collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including HST, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state ten times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The following analysis of NICER spectra attributes variability in the observed X-ray flux to changes in both the column density of the obscurer by at least one order of magnitude ($N_\mathrm{H}$ ranges from $2.85\substack{+0.48\\ -0.33} \times 10^{22}\text{ cm}^{-2}$ to $25.6\substack{+3.0\\ -3.5} \times 10^{22} \text{ cm}^{-2}$) and the intrinsic continuum brightness (the unobscured flux ranges from $10^{-11.8}$ to $10^{-10.5}$ erg s$^{-1}$ cm$^{-2}$ ). While the X-ray flux generally remains in a faint state, there is one large flare during which Mrk 817 returns to its historical mean flux. The obscuring gas is still present at lower column density during the flare but it also becomes highly ionized, increasing its transparency. Correlation between the column density of the X-ray obscurer and the strength of UV broad absorption lines suggests that the X-ray and UV continua are both affected by the same obscuration, consistent with a clumpy disk wind launched from the inner broad line region.

As neutron stars merge they can approach very high nuclear density. Here, we summarized recent results for the evolution and gravitational wave emission from binary neutron star mergers using a a variety of nuclear equations of state with and without a crossover transition to quark matter. We discuss how the late time gravitational wave emission from binary neutron star mergers may possibly reveal the existence of a crossover transition to quark matter.

The CMB is a powerful probe of early-universe physics but is only observed after passing through large-scale structure, which changes the observed spectra in important model-dependent ways. This is of particular concern given recent claims of significant discrepancies with low redshift data sets when a standard $\Lambda$CDM model is assumed. By using empirical measurements of the CMB lensing reconstruction, combined with weak priors on the smoothness of the lensing spectrum, foregrounds, and shape of any additional integrated Sachs-Wolfe effect, we show how the early-universe parameters can be constrained from CMB observations almost independently of the late-time evolution. This provides a way to test new models for early-universe physics, and measure early-universe parameters, independently of late-time cosmology. Using the empirical measurement of lensing keeps the size of the effect of late-time modelling uncertainty under control, leading to only modest increases in error bars of most early-universe parameters compared to assuming a full evolution model. We provide robust constraints on early-$\Lambda$CDM model parameters using the latest Planck PR4 data and show that with future data marginalizing over a single lensing amplitude parameter is sufficient to remove sensitivity to late-time cosmological model only if the spectral shape matches predictions.

Our vision of galaxies has changed significantly since the era of large galaxy surveys like the Sloan, which gave us extensive statistics with millions of galaxies. The Hubble sequence classification described in Chapter 1 still remains very widely used but has been enriched with broad categories based on color that indicate the recent formation of stars: the red sequence of passive galaxies, consisting solely of old stars, and the blue cloud of galaxies with active star formation. Chapter 3 focused on galaxies with a dominant spheroid, which are generally found on the red sequence. One of the key questions about the evolution of galaxies that remain to be answered is to understand how a galaxy can suddenly pass from one category to another. This is a graduate-student level lecture, not a review article.

Since the 1990s, we have known that there is a super-massive black hole in every galaxy, and that its mass is proportional to the mass of the bulge. To better understand how these black holes were formed, in symbiosis with their galaxies, we will look at their demography, the scaling relations between properties of black holes and host galaxies, and their evolution in a Hubble time. Observations at high angular resolution now allow us to enter the black hole sphere of influence, to see the molecular tori evoked in the AGN unification paradigm, and to understand the feeding processes of black holes. These are often accompanied by feedback processes, which moderate the formation of galaxies. This is a graduate-student level lecture, not a review article.

Radio-loud quasars at high redshift (z > 4) are rare objects in the Universe and rarely observed with Very Long Baseline Interferometry (VLBI). But some of them have flux density sufficiently high for monitoring of their apparent position. The instability of the astrometric positions could be linked to the astrophysical process in the jetted active galactic nuclei in the early Universe. Regular observations of the high-redshift quasars are used for estimating their apparent proper motion over several years. We have undertaken regular VLBI observations of several high-redshift quasars at 2.3 GHz (S band) and 8.4 GHz (X band) with a network of five radio telescopes: 40-m Yebes (Spain), 25-m Sheshan (China), and three 32-m telescopes of the Quasar VLBI Network (Russia) -- Svetloe, Zelenchukskaya, and Badary. Additional facilities joined this network occasionally. The sources have also been observed in three sessions with the European VLBI Network (EVN) in 2018--2019 and one Long Baseline Array (LBA) experiment in 2018. In addition, several experiments conducted with the Very Long Baseline Array (VLBA) in 2017--2018were used to improve the time sampling and the statistics. Based on these 37 astrometric VLBI experiments between 2017 and 2021, we estimated the apparent proper motions of four quasars: 0901+697, 1428+422, 1508+572, and 2101+600.

We aim to measure very precise and accurate model-independent masses and distances of detached binary stars. Precise masses at the $< 1$% level are necessary to test and calibrate stellar interior and evolution models, while precise and independent orbital parallaxes are essential to check for the next Gaia data releases. We combined RV measurements with interferometric observations to determine orbital and physical parameters of ten double-lined spectroscopic systems. We report new relative astrometry from VLTI/GRAVITY and, for some systems, new VLT/UVES spectra to determine the radial velocities of each component. We measured the distance of ten binary systems and the mass of their components with a precision as high as 0.03% (average level 0.2%). They are combined with other stellar parameters (effective temperatures, radii, flux ratios, etc.) to fit stellar isochrones and determine their evolution stage and age. We also compared our orbital parallaxes with Gaia and showed that half of the stars are beyond $1\sigma$ with our orbital parallaxes; although, their RUWE is below the frequently used cutoff of 1.4 for reliable Gaia astrometry. By fitting the telluric features in the GRAVITY spectra, we also estimated the accuracy of the wavelength calibration to be $\sim 0.02$% in high and medium spectral resolution modes. We demonstrate that combining spectroscopic and interferometric observations of binary stars provides extremely precise and accurate dynamical masses and orbital parallaxes. As they are detached binaries, they can be used as benchmark stars to calibrate stellar evolution models and test the Gaia parallaxes.

The cosmic-ray ionization rate (CRIR, $\zeta_2$) is one of the key parameters controlling the formation and destruction of various molecules in molecular clouds. However, the current most commonly used CRIR tracers, such as H$_3^+$, OH$^+$, and H$_2$O$^+$, are hard to detect and require the presence of background massive stars for absorption measurements. In this work, we propose an alternative method to infer the CRIR in diffuse clouds using the abundance ratios of OH/CO and HCO$^+$/CO. We have analyzed the response of chemical abundances of CO, OH, and HCO$^+$ on various environmental parameters of the interstellar medium in diffuse clouds and found that their abundances are proportional to $\zeta_2$. Our analytic expressions give an excellent calculation of the abundance of OH for $\zeta_2$ $\leq$10$^{-15}$ s$^{-1}$, which are potentially useful for modelling chemistry in hydrodynamical simulations. The abundances of OH and HCO$^+$ were found to monotonically decrease with increasing density, while the CO abundance shows the opposite trend. With high-sensitivity absorption transitions of both CO (1--0) and (2--1) lines from ALMA, we have derived the H$_2$ number densities ($n_{\rm H_2}$) toward 4 line-of-sights (LOSs); assuming a kinetic temperature of $T_{\rm k}=50\,{\rm K}$, we find a range of (0.14$\pm$0.03--1.2$\pm$0.1)$\times$10$^2$ cm$^{-3}$}. By comparing the observed and modelled HCO$^+$/CO ratios, we find that $\zeta_2$ in our diffuse gas sample is in the { range of $1.0_{-1.0}^{+14.8}$ $\times$10$^{-16}- 2.5_{-2.4}^{+1.4}$ $\times$10$^{-15}$ s$^{-1}$. This is $\sim$2 times higher than the average value measured at higher extinction, supporting an attenuation of CRs as suggested by theoretical models.

Kilonova spectra imprint valuable information about the elements synthesized in neutron star mergers. In AT2017gfo, the kilonova associated with GW170817, the spectroscopic feature centered around 8000 angstroms has been interpreted as the P-Cygni profile arising from singly ionized Strontium. Recently, Perego et al. (2022) suggested that Helium 10833 line can be an alternative explanation of the feature. Here, we study the line features under non-local thermodynamic equilibrium. We find that the ionization of ejecta by the stopping of radioactive decay product can significantly enhance the ionization states around the line forming region. We compute the kilonova spectrum under the assumption of spherical symmetry and uniform elemental fraction in the line-forming region. We find that 0.2\% (in mass) of Helium in the ejecta can reproduce the P-Cygni feature in the observed spectrum at $1.43$ -- $4.40$ days. Strontium with a mass fraction of $1\%$ is also able to make the absorption feature at $\sim 1.5\,$days, but it gets weaker with time due to ionization by radioactive decay products. The strength of the He line signature depends sensitively on UV strength for the first two epochs. Further modeling of UV line blanketing by $r$-process elements and the optical properties of light $r$-process elements would be crucial to distinguish between Helium and Strontium features. The mass fraction of He is a good indicator for ejecta entropy that allows us to probe the mass ejection mechanism.

Proximate molecular quasar absorbers (PH2) are an intriguing population of absorption systems recently uncovered through strong H2 absorption at small velocity separation from the background quasars. We performed a multi-wavelength spectroscopic follow-up of thirteen such systems with VLT/X-Shooter. Here, we present the observations and study the overall chemical enrichment measured from the HI, H2 and metal lines. We combine this with an investigation of the neutral gas kinematics with respect to the quasar host. We find gas-phase metallicities in the range 2% to 40% of the Solar value, i.e. in the upper-half range of HI-selected proximate damped Lyman-alpha systems, but similar to what is seen in intervening H2-bearing systems. This is likely driven by similar selection effects that play against the detection of most metal and molecular rich systems in absorption. Differences are however seen in the abundance of dust (from [Zn/Fe]) and its depletion pattern, when compared to intervening systems, possibly indicating different dust production or destruction close to the AGN. We also note the almost-ubiquitous presence of a high-ionisation phase traced by NV in proximate systems. In spite of the hard UV field from the quasars, we found no strong overall deficit of neutral argon, at least when compared to intervening DLAs. This likely results from argon being mostly neutral in the H2 phase, which actually accounts for a large fraction of the total amount of metals. We measure the quasar systemic redshifts through emission lines from both ionised gas and CO(3-2) emission, the latter being detected in all 6 cases for which we obtained 3-mm data from complementary NOEMA observations. For the first time, we observe a trend between line-of-sight velocity with respect to systemic redshift and metallicity of the absorbing gas. [truncated]

We investigated the diagnostic potential of the Stokes $V$ profile of the H$\alpha$ line to probe the chromospheric line-of-sight (LOS) magnetic field ($B_{\mathrm{LOS}}$) by comparing the $B_{\mathrm{LOS}}$ inferred from the weak field approximation (WFA) with that of inferred from the multi-line inversions of the Ca II 8542 {\AA}, Si I 8536 {\AA} and Fe I 8538 {\AA} lines using the STiC inversion code. Simultaneous spectropolarimetric observations of a pore in the Ca II 8542 {\AA} and H$\alpha$ spectral lines obtained from the SPINOR at the Dunn Solar Telescope on the 4th of December, 2008 are used in this study. The WFA was applied on the Stokes $I$ and $V$ profiles of H$\alpha$ line over three wavelength ranges viz.: around line core ($\Delta\lambda=\pm0.35$ {\AA}), line wings ($\Delta\lambda=[-1.5, -0.6]$ and $[+0.6, +1.5]$ {\AA}) and full spectral range of the line ($\Delta\lambda=\pm1.5$ {\AA}) to derive the $B_{\mathrm{LOS}}$. We found the maximum $B_{\mathrm{LOS}}$ strengths of $\sim+800$ and $\sim+600$ G at $\log\tau_{\mathrm{500}}$ = $-$1 and $-$4.5, respectively in the pore. The morphological map of the $B_{\mathrm{LOS}}$ inferred from the H$\alpha$ line core is similar to the $B_{\mathrm{LOS}}$ map at $\log\tau_{\mathrm{500}}$ = $-$4.5 inferred from multi-line inversions. The $B_{\mathrm{LOS}}$ map inferred from the H$\alpha$ line wings and full spectral range have a similar morphological structure to the $B_{\mathrm{LOS}}$ map inferred at $\log\tau_{\mathrm{500}}$ = $-$1. The $B_{\mathrm{LOS}}$ estimated from H$\alpha$ using WFA is weaker by a factor of $\approx 0.53$ than that of inferred from the multi-line inversions.

Knowing the conserved quantities that a galaxy's stellar orbits conform to is important in helping us understand the stellar distribution and structures within the galaxy. Isolating integrals of motion and resonances are particularly important, non-isolating integrals less so. We compare the behavior and results of two methods for counting the number of conserved quantities, one based on the correlation integral approach and the other a more recent method using machine learning. Both methods use stellar orbit trajectories in phase space as their only input, and we create such trajectories from theoretical spherical, axisymmetric and triaxial model galaxies. The orbits have known isolating integrals and resonances. We find that neither method is fully effective in recovering the numbers of these quantities, nor in determining the number of non-isolating integrals. From a computer performance perspective, we find the correlation integral approach to be the faster. Determining the algebraic formulae of (multiple) conserved quantities from the trajectories has not been possible due to the lack of an appropriate symbolic regression capability. Notwithstanding the shortcomings we have noted, it may be that the methods are usable as part of a trajectory analysis tool kit.

Solar coronal jets are frequently occurring collimated ejections of solar plasma, originating from magnetically mixed polarity locations on the Sun of size scale comparable to that of a supergranule. Many, if not most, coronal jets are produced by eruptions of small-scale filaments, or minifilaments, whose magnetic field reconnects both with itself and also with surrounding coronal field. There is evidence that minifilament eruptions are a scaled-down version of typical filament eruptions that produce solar flares and coronal mass ejections (CMEs). Moreover, the magnetic processes building up to and triggering minifilament eruptions, which is often flux cancelation, might similarly build up and trigger the larger filaments to erupt. Thus, detailed study of coronal jets will inform us of the physics leading to, triggering, and driving the larger eruptions. Additionally, such studies potentially can inform us of smaller-scale coronal-jet-like features, such as jetlets and perhaps some spicules, that might work the same way as coronal jets. We propose a high-resolution (~0.1 pixels), high-cadence (~5 seconds) EUV-solar-imaging mission for the upcoming decades, that would be dedicated to observations of features of the coronal-jet size scale, and smaller-scale solar features produced by similar physics. Such a mission could provide invaluable insight into the operation of larger features such as CMEs that produce significant Space Weather disturbances, and also smaller-scale features that could be important for coronal heating, solar wind acceleration, and heliospheric features such as the magnetic switchbacks that are frequently observed in the solar wind.

Europa is an icy ocean world, differentiated into a floating ice shell and solid interior, separated by a global ocean. The classical spin-orbit coupling problem considers a satellite as a single rigid body, but in the case of Europa, the existence of the subsurface ocean enables independent motion of the ice shell and solid interior. This paper explores the spin-orbit coupling problem for Europa from a dynamical perspective, yielding illuminating analytical and numerical results. We determine that the spin behavior of Europa is influenced by processes not captured by the classical single rigid body spin-orbit coupling analysis. The tidal locking process for Europa is governed by the strength of gravity-gradient coupling between the ice shell and solid interior, with qualitatively different behavior depending on the scale of this effect. In this coupled rigid model, the shell can potentially undergo large angular displacements from the solid interior, and the coupling plays an outsize role in the dynamical evolution of the moon, even without incorporating the dissipative effects of shell non-rigidity. This work extends classical spin-orbit coupling analysis to the icy ocean worlds, lending a new perspective for dynamically exploring the rotational history of Europa and other similarly differentiated bodies in the outer solar system.

The 'optical depth-linear growth rate' ($\tau_{\rm T}-f$) degeneracy is a long-standing problem in the kinetic Sunyaev Zel'dovich (kSZ) cosmology. It can be spontaneously broken in redshift space, where the velocity field leaves its own distinct imprint on the galaxies' redshift space positions and provides valuable information of the linear growth rate. We validate this idea with the Fisher matrix and Monte Carlo Markov Chain techniques in this work, finding that the level of this degeneracy breaking is further enhanced on non-linear scales due to the non-linear evolution of the density and velocity fields. This result emphasizes the importance of the redshift space analysis of the kSZ effect and its potential as a powerful cosmological probe, especially on non-linear scales. As a by-product, we develop a non-linear model of the redshift space density-weighted pairwise kSZ power spectrum. The fitted $f$ and $\tau_{\rm T}$ values from this model are shown to be accurate within $1-2\sigma$ ranges of the fiducial ones when confronted to the mock galaxies mimicking a DESI+CMB-S4 survey combination, even on small scales of $k\sim 0.5h/{\rm Mpc}$.

Accretion processes in protoplanetary disks produce a diversity of small bodies that played a crucial role in multiple reshuffling events throughout the solar system and in both early and late accretion of planets. Application of thermo-chronometers to meteorites provides precise dating of the formation age of mineralogical components. Nucleosynthetic anomalies indicating a dichotomy between NC and C meteorites and precise parent body (PB) chronology can be combined with thermal evolution models to constrain the timescale of accretion and dynamical processes in the early solar system. Achondrite PBs are considered to have accreted early and mostly in the NC region, while late accretion in the C region produced mostly undifferentiated PBs, such as the CR chondrite PB that formed as late as 4 Ma after CAIs. However, presence of more evolved CR-related meteorites suggests also an earlier accretion timing. We present modeling evidence for a temporally distributed accretion of parent bodies of CR-related meteorite groups that originate from a C reservoir and range from aqueously altered chondrites to partially differentiated primitive achondrites. The PB formation times derived range from <1 Ma to ~4 Ma after solar system formation, with ~3.7 Ma, ~1.5-2.75 Ma, <~0.6 Ma, and <~0.7 Ma for CR, Flensburg, NWA 6704, and NWA 011. This implies that accretion processes in the C reservoir started as early as in the NC reservoir and produced differentiated PBs with carbonaceous compositions in addition to undifferentiated C chondritic PBs. The accretion times correlate inversely with the degree of the meteorites' alteration, metamorphism, or differentiation. Accretion times for CI/CM, Ryugu, and Tafassites PBs of ~3.75 Ma, ~1-3 Ma, and 1.1 Ma, respectively, fit well into this correlation in agreement with the thermal and alteration conditions suggested by the meteorites.

We present a general framework on how the polarization of radiation due to scattering, dichroic extinction, and birefringence of aligned spheroidal dust grains can be implemented and tested in 3D Monte Carlo radiative transfer (MCRT) codes. We derive a methodology for solving the radiative transfer equation governing the changes of the Stokes parameters in dust-enshrouded objects. We utilize the M\"uller matrix, and the extinction, scattering, linear, and circular polarization cross sections of spheroidal grains as well as electrons. An established MCRT code is used and its capabilities are extended to include the Stokes formalism. We compute changes in the polarization state of the light by scattering, dichroic extinction, and birefringence on spheroidal grains. The dependency of the optical depth and the albedo on the polarization is treated. The implementation of scattering by spheroidal grains both for random walk steps as well as for directed scattering (peel-off) are described. The observable polarization of radiation of the objects is determined through an angle binning method for photon packages leaving the model space as well as through an inverse ray-tracing routine for the generation of images. We present paradigmatic examples for which we derive analytical solutions of the optical light polarization by spheroidal dust particles. These tests are suited for benchmark verification of MCpol and other such codes, and allow to quantify the numerical precision reached. We demonstrate that MCpol is in excellent agreement to within 0.1% of the Stokes parameters when compared to the analytical solutions.

We report the discovery of SRGE J075818-612027, a deep stream-eclipsing magnetic cataclysmic variable found serendipitously in SRG/eROSITA CalPV observations of the open cluster NGC~2516 as an unrelated X-ray source. An X-ray timing and spectral analysis of the eROSITA data is presented and supplemented by an analysis of TESS photometry and SALT spectroscopy. X-ray photometry reveals two pronounced dips repeating with a period of $106.144(1)$ min. The 14-month TESS data reveal the same unique period. A low-resolution identification spectrum obtained with SALT displays hydrogen Balmer emission lines on a fairly blue continuum. The spectrum and the stability of the photometric signal led to the classification of the new object as a polar-type cataclysmic variable. In this picture, the dips in the X-ray light curve are explained by absorption in the intervening accretion stream and by a self-eclipse of the main accretion region. The object displays large magnitude differences on long (months) timescales both at optical and X-ray wavelengths, being interpreted as high and low states and thus supporting the identification as a polar. The bright phase X-ray spectrum can be reflected with single temperature thermal emission with 9.7 keV and bolometric X-ray luminosity $L_{\rm X} \simeq 8\times 10^{32}$erg s$^{-1}$ at a distance of about 2.7 kpc. It lacks the pronounced soft X-ray emission component prominently found in ROSAT-discovered polars.

An isolated pulsar is a rotating neutron star possessing a high magnetic dipole moment that generally makes a finite angle with its rotation axis. As a consequence, the emission of magnetic dipole radiation (MDR) continuously takes away its rotational energy. This process leads to a time decreasing angular velocity of the star that is usually quantified in terms of its braking index. While this simple mechanism is indeed the main reason for the spin evolution of isolated pulsars, it may not be the only cause of this effect. Most of young isolated pulsars present braking index values that are consistently lower than that given by the MDR model. Working in the weak field (Newtonian) limit, we take in the present work a step forward in describing the evolution of such a system by allowing the star's shape to wobble around an ellipsoidal configuration as a backreaction effect produced by the MDR emission. It is assumed that an internal damping of the oscillations occurs, thus introducing another form of energy loss in the system, and this phenomenon may be related to the deviation of the braking index from the pure MDR model predictions. Numerical calculations suggest that the average braking index for typical isolated pulsars can be thus simply explained.

Scintillations of pulsar radio signals caused by the interstellar medium can in principle be used for interstellar interferometry. Changes of the dynamic spectra as a function of pulsar longitude were in the past interpreted as having spatially resolved the pulsar magnetosphere. Guided by this prospect we used VLBI observations of PSR B1237+25 with the Arecibo and Green Bank radio telescopes at 324 MHz and analyzed such scintillations at separate longitudes of the pulse profile. We found that the fringe phase characteristics of the visibility function changed quasi-sinusoidally as a function of longitude. Also, the dynamic spectra from each of the telescopes shifted in frequency as a function of longitude. Similar effects were found for PSR B1133+16. However, we show that these effects are not signatures of having resolved the pulsar magnetosphere. Instead the changes can be related to the effect of low-level digitizing of the pulsar signal. After correcting for these effects the frequency shifts largely disappeared. Residual effects may be partly due to feed polarization impurities. In view of our analysis we think that observations with the intend of spatially resolving the pulsar magnetosphere need to be critically evaluated in terms of these constraints on interstellar interferometry.

Weak gravitational lensing of the cosmic microwave background (CMB) carries imprints of the physics operating at redshifts much lower than that of recombination and serves as an important probe of cosmological structure formation, dark matter physics, and the mass of neutrinos. Reconstruction of the CMB lensing deflection field through use of quadratic estimators has proven successful with existing data but is known to be sub-optimal on small angular scales ($\ell > 3000$) for experiments with low noise levels. Future experiments will provide better observations in this regime, but these techniques will remain statistically limited by their approximations. We show that correlations between fluctuations of the large-scale temperature gradient power of the CMB sourced by $\ell < 2000$, and fluctuations to the local small-scale temperature power reveal a lensing signal which is prominent in even the real-space pixel statistics across a CMB temperature map. We present the development of the Small Correlated Against Large Estimator (SCALE), a novel estimator for the CMB lensing spectrum which offers promising complementary analysis alongside other reconstruction techniques in this regime. The SCALE method computes correlations between both the large/small-scale temperature gradient power in harmonic space, and it is able to quantitatively recover unbiased statistics of the CMB lensing field without the need for map-level reconstruction. SCALE can outperform quadratic estimator signal-to-noise by a factor of up to 1.5 in current and upcoming experiments for CMB lensing power spectra $C_{6000<L<8000}^{\phi\phi}$.

The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this RV information remains useful for studying the stellar multiplicity of nearby M~dwarfs. To search for companions around such ``deprioritized" M~dwarfs, we observed 14 IRD-SSP targets using Keck/NIRC2 observations with pyramid wavefront sensing at NIR wavelengths, leading to high sensitivity to substellar-mass companions within a few arcseconds. We detected two new companions (LSPM~J1002+1459~B and LSPM~J2204+1505~B) and two new candidates that are likely companions (LSPM~J0825+6902~B and LSPM~J1645+0444~B) as well as one known companion. Including two known companions resolved by the IRD fiber injection module camera, we detected seven (four new) companions at projected separations between $\sim2-20$~au in total. A comparison of the colors with the spectral library suggests that LSPM~J2204+1505~B and LSPM~J0825+6902~B are located at the boundary between late-M and early-L spectral types. Our deep high-contrast imaging for targets where no bright companions were resolved did not reveal any additional companion candidates. The NIRC2 detection limits could constrain potential substellar-mass companions ($\sim10-75\ M_{\rm Jup}$) at 10~au or further. The failure with Keck/NIRC2 around the IRD-SSP stars having significant RV trends makes these objects promising targets for further RV monitoring or deeper imaging with JWST to search for smaller-mass companions below the NIRC2 detection limits.

We present the first result in exploring the gaseous halo and galaxy correlation using the Dark Energy Spectroscopic Instrument (DESI) survey validation data in the Cosmic Evolution Survey (COSMOS) and Hyper Suprime-Cam (HSC) field. We obtain the multiphase gaseous halo properties in the circumgalactic medium (CGM) by using 115 quasar spectra (S/N > 3). We detect MgII absorption at redshift 0.6 < z < 2.5, CIV absorption at 1.6 < z < 3.6, and HI absorption associated with the MgII and CIV. The CGM is mixed by a higher density phase of detectable MgII and CIV and a lower density of CIV-only phase. By cross-matching the COSMOS2020 catalog, we identify the MgII and CIV host galaxies at 0.9 < z < 3.1 in ten quasar fields. We find that within the impact parameter of 250 kpc, a tight correlation is seen between strong MgII equivalent width and the host galaxy star formation rate. The covering fraction fc of strong MgII selected galaxies, which is the ratio of absorbing galaxy in a certain galaxy population, shows significant evolution in the main-sequence galaxies and marginal evolution in all the galaxy populations within 250 kpc at 0.9 < z < 2.2. The fc increase in the main-sequence galaxies likely suggests the co-evolution of strong MgII absorbing gas and the main-sequence galaxies at the cosmic noon. Furthermore, several MgII and CIV absorbing gas is detected out of the galaxy virial radius, tentatively indicating the feedback produced by the star formation and/or the environmental effects.

The timing follow-up of newly discovered millisecond pulsars (MSPs) is hindered by the larger positional uncertainty (a few tens of arc-minutes) associated with the discovery. In this paper, we present the localization of two MSPs, discovered by the GMRT High Resolution Southern Sky (GHRSS) survey, up to arc-second accuracy using a 33~MHz offline coherently dedispersed gated correlator. This gated correlator is an upgraded version of the earlier 16~MHz design. This new development with a factor of two enhancement in the observing bandwidth offers better sensitivity in the image domain, leading to more precise localization. Aided by the precise position, we followed up these two MSPs with sensitive phased-array (PA) beams of upgraded GMRT from 300 to 1460 MHz. More sensitive observations in the PA mode for these two MSPs yield precise ($\sim$ sub-$\mu s$) time-of-arrivals, with DM uncertainties in the range of $10^{-4}-10^{-5}$ $pc\,cm^{-3}$. We also report the profile evolution of the two MSPs over 300$-$1460 MHz. Finally, we discuss the suitability of these MSPs for the pulsar timing array experiments aimed to detect low-frequency gravitational waves.

The cosmic dust particles found in space are mainly porous aggregates of smaller grains. Theoretically, these aggregates are replicated using fractal geometry, assuming a cluster of spheres. Although, the light scattering response of cosmic dust aggregates has been thoroughly studied using clusters of spherical grains in the past few decades, yet, the effect of irregularities on the surface of each grain in an entire aggregate has mostly been neglected. We, for the first time, introduce a visually realistic cosmic dust model which incorporates a mixture of rough fractal aggregates (RFA) and agglomerated debris (Solids) to replicate the unusual polarization-phase curve observed in case of the interstellar comet 2I/Borisov at multiple wavelengths. The authenticity of the RFA structures has been verified by replicating light scattering results of circumstellar dust analogues from the Granada Amsterdam Light Scattering Database. We demonstrate that the light scattering response from the RFA structures has a very close resemblance with the experimental values. Finally, we model the observed polarization-phase curve of the interstellar comet 2I/Borisov using a mixture of RFA and solid particles. The best-fit data indicates presence of higher percentage of porous RFA structures 80% owing to the fact that the comet carries higher percentage of small and highly porous pristine cosmic dust particles. Further, the model indicates that the unusually steeper polarimetric slope and the high dust-to-gas ratio in relatively newer comets is mainly due to higher porous-to-compact ratio.

We report on Konus-Wind (KW) and Mikhail Pavlinsky ART-XC telescope observations and analysis of a nearby GRB 221009A, the brightest $\gamma$-ray burst (GRB) detected by KW for $>$28 years of observations. The pulsed prompt phase of the burst emission lasts for $\sim 600$~s and is followed by a steady power-law decay lasting for more than 20~ks. From the analysis of the KW and ART-XC light curves and the KW spectral data we derive time-averaged spectral peak energy of the burst $E_p\approx 2.6$~MeV, $E_p$ at the brightest emission peak $E_p\approx 3.0$~MeV, the total 20~keV--10~MeV energy fluence of $\approx0.21$~erg~cm$^{-2}$, and the peak energy flux in the same band of $\approx 0.03$~erg~cm$^{-2}$~s$^{-1}$. The enormous observed fluence and peak flux imply, at redshift $z=0.151$, huge values of isotropic energy release $E_{\mathrm{iso}}\approx1.2\times10^{55}$~erg (or $\gtrsim 6.5$~solar rest mass) and isotropic peak luminosity $L_{\mathrm{iso}}\approx3.4\times10^{54}$~erg~s$^{-1}$ (64-ms scale), making GRB~221009A the most energetic and one the most luminous bursts observed since the beginning of the GRB cosmological era in 1997. The isotropic energetics of the burst fit nicely both `Amati' and `Yonetoku' hardness-intensity correlations for $>$300~KW long GRBs, implying that GRB~221009A is most likely a very hard, super-energetic version of a "normal" long GRB.

In the multi-messenger era, coordinating observations between astronomical facilities is mandatory to study transient phenomena (e.g. Gamma-ray bursts) and is achieved by sharing information with the scientific community through networks such as the Gamma-ray Coordinates Network. The facilities usually develop real-time scientific analysis pipelines to detect transient events, alert the astrophysical community, and speed up the reaction time of science alerts received from other observatories. We present in this work the RTApipe framework, designed to facilitate the development of real-time scientific analysis pipelines for present and future gamma-ray observatories. This framework provides pipeline architecture and automatisms, allowing the researchers to focus on the scientific aspects and integrate existing science tools developed with different technologies. The pipelines automatically execute all the configured analyses during the data acquisition. This framework can be interfaced with science alerts networks to perform follow-up analysis of transient events shared by other facilities. The analyses are performed in parallel and can be prioritised. The workload is highly scalable on a cluster of machines. The framework provides the required services using containerisation technology for easy deployment. We present the RTA pipelines developed for the AGILE space mission and the prototype of the SAG system for the ground-based future Cherenkov Telescope Array observatory confirming that the RTApipe framework can be used to successfully develop pipelines for the gamma-ray observatories, both space and ground-based.

Magnetic reconnection has long been known to be the most important mechanism as quick conversion of magnetic field energy into plasma kinetic energy. In addition, energy dissipation by reconnection has gained attention not only as a plasma heating mechanism, but also as a plasma mechanism for accelerating nonthermal particles. However, the energy partitioning of thermal and nonthermal plasmas during magnetic reconnection is not understood. Here, we studied energy partition as a function of plasma sheet temperature and guide magnetic field. In relativistic reconnection with anti-parallel magnetic field or weak guide magnetic field, it was found that the nonthermal energy density can occupy more than $90 \%$ of the total kinetic plasma energy density, but strengthening the guide magnetic field suppresses the efficiency of the nonthermal particle acceleration. In nonrelativistic reconnection for anti-parallel magnetic field, most dissipated magnetic field energy is converted into thermal plasma heating. For a weak guide magnetic field with a moderate value, however, the nonthermal particle acceleration efficiency was enhanced, but strengthening the guide-field beyond the moderate value suppresses the efficiency.

Gradual solar energetic particle (SEP) events, usually attributed to shock waves driven by coronal mass ejections (CMEs), show a wide variety of temporal behaviors. For example, TO, the >10 MeV proton onset time with respect to the launch of the CME, has a distribution of at least an order of magnitude, even when the source region is not far from the so-called well-connected longitudes. It is important to understand what controls TO, especially in the context of space weather prediction. Here we study two SEP events from the western hemisphere that are different in TO on the basis of >10 MeV proton data from the Geostationary Operations Environmental Satellite, despite similar in the CME speed and longitude of the source regions. We try to find the reasons for different TO, or proton release times, in how the CME-driven shock develops and the Alfv\'en Mach number of the shock wave reaches some threshold, by combining the CME height-time profiles with radio dynamic spectra. We also discuss how CME-CME interactions and active region properties may affect proton release times.

We complete the publication of all microlensing planets (and ``possible planets'') identified by the uniform approach of the KMT AnomalyFinder system in the 21 KMT subprime fields during the 2019 observing season, namely KMT-2019-BLG-0298, KMT-2019-BLG-1216, KMT-2019-BLG-2783, OGLE-2019-BLG-0249, and OGLE-2019-BLG-0679 (planets), as well as OGLE-2019-BLG-0344, and KMT-2019-BLG-0304 (possible planets). The five planets have mean log mass-ratio measurements of $(-2.6,-3.6,-2.5,-2.2,-2.3)$, median mass estimates of $(1.81,0.094,1.16,7.12,3.34)\, M_{\rm Jup}$, and median distance estimates of $(6.7,2.7,5.9,6.4,5.6)\, {\rm kpc}$, respectively. The main scientific interest of these planets is that they complete the AnomalyFinder sample for 2019, which has a total of 25 planets that are likely to enter the statistical sample. We find statistical consistency with the previously published 33 planets from the 2018 AnomalyFinder analysis according to an ensemble of five tests. Of the 58 planets from 2018-2019, 23 were newly discovered by AnomalyFinder. Within statistical precision, half of all the planets have caustic crossings while half do not (as predicted by Zhu et al. 2014), an equal number of detected planets result from major-image and minor-image light-curve perturbations, and an equal number come from KMT prime fields versus subprime fields.

Primordial B-mode detection is one of the main goals of the current and future CMB experiments. However, the weak B-mode signal is overshadowed by several Galactic polarized emissions, such as the thermal dust emission and the synchrotron radiation. Subtracting the foreground components from CMB observations is one of the key challenges in the search for primordial B-mode signal. Here, we construct a deep convolutional neural network (CNN) model, called CMBFSCNN (Cosmic Microwave Background Foreground Subtraction with CNN), which can cleanly remove various foreground components from the simulated CMB observational maps with the sensitivity of the CMB-S4 experiment. The noisy CMB Q (or U) maps are recovered with a mean absolute difference of $0.018 \pm 0.023\ \mu$K (or $0.021 \pm 0.028\ \mu$K). To remove residual instrumental noise in the foreground-cleaned map, inspired by the Needlet Internal Linear Combination method, we divide the whole data into two ``half-split maps'' which share the same sky signal but with uncorrelated noise, and perform the cross-correlation technique to reduce the instrumental noise effect at the power spectrum level. We find that the CMB EE and BB power spectra can be precisely recovered with significantly reduced noise effects. Finally, we apply this pipeline on the current Planck observations. As expected, various foregrounds have been cleanly removed on the Planck observational maps and the recovered EE and BB power spectra are in good agreement with the Planck official results.

The power-law emission and reflection component provide valuable insights into the accretion process around a black hole. In this work, thanks to the broadband spectra coverage of \emph{the Nuclear Spectroscopic Telescope Array}, we study the spectral properties for a sample of low-mass black hole X-ray binaries (BHXRBs). We find that there is a positive correlation between the photon index $\Gamma$ and the reflection fraction $R$ (the ratio of the coronal intensity that illuminates the disk to the coronal intensity that reaches the observer), consistent with previous studies, but except for MAXI J1820+070. It is quite interesting that this source also deviates from the well-known ``V"-shaped correlation between the photon index $\Gamma$ and the X-ray luminosity log$L_{\rm X}$, when it is in the bright hard state. More specifically, the $\Lambda$-shaped correlation between $\Gamma$ and log$L_{\rm X}$ is observed, as the luminosity decreases by a factor of 3 in a narrow range from $\sim 10^{38}$ to $10^{37.5}$ $\rm erg~s^{-1}$. Furthermore, we discover a strong positive correlation between $R$ and the X-ray luminosity for BHXRBs in the hard state, which puts a constraint on the disk-corona coupling and the evolution.

The ASTRI Mini-Array is an international collaboration led by the Italian National Institute for Astrophysics. This project aims to construct and operate an array of nine Imaging Atmospheric Cherenkov Telescopes to study gamma-ray sources at very high energy (TeV) and perform stellar intensity interferometry observations. We describe the software architecture and the technologies used to implement the Online Observation Quality System (OOQS) for the ASTRI Mini-Array project. The OOQS aims to execute data quality checks on the data acquired in real-time by the Cherenkov cameras and intensity interferometry instruments, and provides feedback to both the Central Control System and the Operator about abnormal conditions detected. The OOQS can notify other sub-systems, triggering their reaction to promptly correct anomalies. The results from the data quality analyses (e.g. camera plots, histograms, tables, and more) are stored in the Quality Archive for further investigation and they are summarised in reports available to the Operator. Once the OOQS results are stored, the operator can visualize them using the Human Machine Interface. The OOQS is designed to manage the high data rate generated by the instruments (up to 4.5 GB/s) and received from the Array Data Acquisition System through the Kafka service. The data are serialized and deserialized during the transmission using the Avro framework. The Slurm workload scheduler executes the analyses exploiting key features such as parallel analyses and scalability.

The standardisation of gamma-ray astronomical data emerged in recent years as a necessity for the future generation of gamma-ray observatories. Nevertheless, adopting a common format for gamma-ray instruments can already benefit the current generation of gamma-ray instruments. As the end of their operations approaches, it provides a natural solution for the production of their data legacy. Additionally it eases data combination for multi-instrument analyses, thus enhancing the potential for scientific discovery with the wealth of data so far gathered. In this contribution, we present for the first time the effort to adapt the data of the MAGIC telescopes to the standardised format proposed by the Data Formats for Gamma-ray Astronomy initiative. We validate the data conversion by analysing the standardised data with the open-source software Gammapy and comparing the results obtained against those produced with the MAGIC proprietary software, MARS. For both samples chosen (Crab Nebula and Mrk421 observation), for all the scientific products extracted (spectra and light curves), we observe good agreement between the results of the two software.

On June 30, 1908 at about 0h 14.5m UTC what is known today as the Tunguska Event (TE) occurred, most likely caused by the fall of a small rocky asteroid of about 50-60 meters in diameter over the basin of the Tunguska River (Central Siberia). Unfortunately the first expedition was made by Kulik 19 years after the event and macroscopic meteorites have never been found in epicenter site. After considering the Chelyabinsk event as a guide, we estimated the strewn field of possible macroscopic fragments of the asteroid responsible of the TE: we have reason to believe that there might be fragments with enough strength to survive the airburst and reach the ground. The strewn field, which is located about 15 to 20 km North-West from the epicenter, forms an ellipse with axes from $25 \times 20$ to $40 \times 30$ km at $3\sigma$ level and should be considered for the search of macroscopic bodies, even if the mud and vegetation could have made any trace disappear. Cheko Lake, which by some authors is considered an impact crater, falls about 3 km ($\sim 4\sigma$) outside these areas and, based on our results, it is unlikely that it could be a real impact crater: only if the cosmic body's trajectory had an azimuth of about $160^\circ - 170^\circ$ would be in the strewn field area, but it is not consistent with the most likely trajectory.

This paper describes the polarisation study of a Lynds cloud, LDN 1340, $\alpha$ = 2h32m & $\delta$ = $73^{\circ} 00^\prime$ corresponding to galactic coordinates of $\ell=$ 130$^{\circ}$.07 $b=$ 11$^{\circ}$.6, with emphasis on the RNO 8 area. The cloud has been observed using the 1.2 m telescope at Mt.Abu Infrared Observatory, in the infrared wavelength band using the Near-Infrared Camera, Spectrograph & Polarimeter (NICSPol) instrument. The polarimetric observations were used to map the magnetic field geometry around the region. We combined our measurements with archival data from the 2MASS and WISE surveys. The Gaia EDR3 & DR3 data for the same region were used for distance, proper motion, and other astrophysical information. The analysis of the data reveals areas with ordered polarisation vectors in the region of RNO 8. The position angle measurements reveal polarisation due to dichroic extinction which is consistent with the Galactic magnetic field. The magnetic field strength was calculated for the RNO 8 region using the Chandrashekhar-Fermi method and the value estimated is $\sim$ 42$\mu$G.

\textbf{Context:} White-light continuum observations of solar flares often have coronal counterparts, including the classical ``white-light prominence'' (WLP) phenomenon. \textbf{Aims:} Coronal emissions by flares, seen in white-light continuum, have only rarely been reported previously. We seek to use modern data to understand the morphology of WLP events. \textbf{Methods:} We have identified a set of 14 examples of WLP detected by the HMI (Heliospheric and Magnetic Imager) experiment on board SDO (the Solar Dynamics Observatory satellite), using a new on-line catalogue covering 2011-2017. These invariably accompanied white-light flares (WLF) emission from the lower atmosphere by flares near the limb, as identified by hard X-ray images from RHESSI (the Reuven Ramaty High Energy Spectroscopic Imager). HMI provides full Stokes information, and we have used the linear polarisations (Q~and~U) to distinguish Thomson scattering from cool material, following the analysis pioneered by \cite{2014ApJ...786L..19S}. \textbf{Results:} The event morphologies fit roughly into three categories: Ejection, Loop, and Spike, but many events show multiple phenomena. \textbf{Conclusions:} The coronal white-light continuum, observed by HMI analogously to the observations made by a coronagraph, detects many examples of coronal emission and dynamics. Using the Stokes linear polarisation, we estimate the masses of hot coronal plasma in 11 of the 14 events and find them to be similar to typical CME masses, but not exceeding 10$^{15}$\,g. We note that the HMI observations do not occult the bright solar disk and were not designed for coronal observations, resulting in relatively low signal-to-noise ratios. We therefore believe that future such observations with better optimisation will be even more fruitful.

In this work we have analyzed the variation of speed of sound and trace anomaly inside a neutron star using neural network. We construct family of agnostic equation of state by maintaining thermodynamic stability, speed of sound bound with constraints from recent observation of neutron stars. The data of mass and radius serves an input for the neural network and we obtain speed of sound as an output. The speed of sound shows non-monotonic behaviour and the trained data predicts softer equation of state. The neural network predicts stiff equation of state at the centre of intermediate mass stars and massive stars have relatively softer equation of state at their centre. The trace anomaly shows a non-monotonic behaviour for the untrained data however, when we train our data; with neural network the trace anomaly shows monotonic behaviour and the condition of $\Delta \ge 0$ is maintained.

Context. HESS J1809$-$193 is an unassociated very-high-energy $\gamma$-ray source located on the Galactic plane. While it has been connected to the nebula of the energetic pulsar PSR J1809$-$1917, supernova remnants and molecular clouds present in the vicinity also constitute possible associations. Recently, the detection of $\gamma$-ray emission up to energies of $\sim$100 TeV with the HAWC observatory has led to renewed interest in HESS J1809$-$193. Aims. We aim to understand the origin of the $\gamma$-ray emission of HESS J1809$-$193. Methods. We analysed 93.2 h of data taken on HESS J1809$-$193 above 0.27 TeV with the High Energy Stereoscopic System (H.E.S.S.), using a multi-component, three-dimensional likelihood analysis. In addition, we provide a new analysis of 12.5 yr of Fermi-LAT data above 1 GeV within the region of HESS J1809$-$193. The obtained results are interpreted in a time-dependent modelling framework. Results. For the first time, we were able to resolve the emission detected with H.E.S.S. into two components: an extended component that exhibits a spectral cut-off at $\sim$13 TeV, and a compact component that is located close to PSR J1809$-$1917 and shows no clear spectral cut-off. The Fermi-LAT analysis also revealed extended $\gamma$-ray emission, on scales similar to that of the extended H.E.S.S. component. Conclusions. Our modelling indicates that based on its spectrum and spatial extent, the extended H.E.S.S. component is likely caused by inverse Compton emission from old electrons that form a halo around the pulsar wind nebula. The compact component could be connected to either the pulsar wind nebula or the supernova remnant and molecular clouds. Due to its comparatively steep spectrum, modelling the Fermi-LAT emission together with the H.E.S.S. components is not straightforward. (abridged)

The Hubble tension is investigated taking into account the cosmological look-back time. Specifically, considering a single equation, widely used in standard cosmology, it is possible to recover both values of the Hubble constant $H_0$ reported by the SHOES and Planck collaborations: the former is obtained through cosmological ladder methods (e.g. Cepheids, Supernovae Type IA) and the latter through measurements of the Cosmic Microwave Background. Also, other values obtained in the literature are achieved with the same approach. We conclude that the Hubble tension can be removed if the look-back time is correctly referred to the redshift where the measurement is performed.

The goal of this study is to reconstruct the evolution and the dust formation processes during the final AGB phases of a sample of carbon-rich, post-AGB Galactic stars, with particular attention to the determination of the past mass-loss history. We study the IR excess of sources classified as single stars by means of dust formation modelling where dust grains form and grow in a static wind and expand from the surface of the star. The method is applied to various evolutionary stages of the final AGB phase of stars with different masses and metallicities. The detailed analysis of the SED of the sources investigated, which included the derivation of the luminosities and the dust properties, is used to infer information on mass loss, efficiency of dust formation, and wind dynamics. We confirm previous results that most of the investigated sources descend from low-mass(M<1.5Msun) progenitors that reached the C-star stage. Metal-poor carbon stars are characterised by higher IR excesses with respect to their more metal-rich counterparts of similar luminosity due to a higher surface carbon-to-oxygen excess. This work confirms previous conclusions that more luminous stars descending from higher-mass progenitors are generally more opaque due to shorter evolutionary timescales that place the dust shell closer to the central object. We also find that the mass-loss rate at the tip of the AGB phase of metal-rich low-mass carbon stars is approximately 1-1.5x10^-5Msun/yr, whereas in the metal-poor domain M~4-5x10^-5Msun/yr is required. These results indicate the need for an upwards revision of the theoretical mass-loss rates of low-mass carbon stars in the available literature, which in turn require a revised determination of carbon dust yields by AGB stars.

We present a method for the optical identification of sources detected in wide-field X-ray sky surveys. We have constructed and trained a neural network model to characterise the photometric attributes of the populations of optical counterparts of X-ray sources and optical field objects. The photometric information processing result is used for the probabilistic cross-match of X-ray sources with optical DESI Legacy Imaging Surveys sources. The efficiency of the method is illustrated using the SRG/eROSITA Survey of Lockman Hole. To estimate the accuracy of the method, we have produced a validation sample based on the Chandra and XMM-Newton catalogues of X-ray sources. The cross-match precision in our method reaches 94% for the entire X-ray catalogue and 97% for sources with a flux $F_{\rm x, 0.5-2}>10^{-14}$ erg/s/cm$^2$. We discuss the further development of the optical identification model and the steps needed for its application to the SRG/eROSITA all-sky survey data.

While the inflationary paradigm is consistent with observations, it still suffers from the singularity problem. On the other hand, the classical bouncing scenario, while devoid of this problem, does not conform to observations. Here, we explore scenarios wherein the bouncing phase smoothly transits to an inflationary one, with the pivot scale leaving the Hubble horizon during the latter era, thereby maintaining consistency with observations. Staying within the ambit of Einstein-Hilbert gravity augmented by the inflaton, we achieve the bounce by introducing a second scalar field that helps engineer the requisite violation of the null energy condition. Potential ghost instabilities can be mitigated by invoking a non-trivial coupling between the two scalar fields.

Dwarf galaxies are excellent cosmological probes, because their shallow potential wells make them very sensitive to the key processes that drive galaxy evolution, including baryonic feedback, tidal interactions, and ram pressure stripping. However, some of the key parameters of dwarf galaxies, which help trace the effects of these processes, are still debated, including the relationship between their sizes and masses. We re-examine the Fornax Cluster dwarf population from the point of view of isomass-radius--stellar mass relations (IRSMRs) using the Fornax Deep Survey Dwarf galaxy Catalogue, with the centrally located (among dwarfs) $3.63 \mathcal{M}_{\odot}$~pc$^{-2}$ isodensity radius defining our fiducial relation. This relation is a powerful diagnostic tool for identifying dwarfs with unusual structure, as dwarf galaxies' remarkable monotonicity in light profile shapes, as a function of stellar mass, reduces the relation's scatter tremendously. By examining how different dwarf properties (colour, tenth-nearest-neighbour distance, etc.) correlate with distance from our fiducial relation, we find a significant population of structural outliers with comparatively lower central mass surface density and larger half-light-radii, residing in locally denser regions in the cluster, albeit with similar red colours. We propose that these faint, extended outliers likely formed through tidal disturbances, which make the dwarfs more diffuse, but with little mass loss. Comparing these outliers with ultra-diffuse galaxies (UDGs), we find that the term UDG lacks discriminatory power; UDGs in the Fornax Cluster lie both on and off of IRSMRs defined at small radii, while IRSMR outliers with masses below $\sim 10^{7.5} \mathcal{M}_{\odot}$ are excluded from the UDG classification due to their small effective radii.

We describe the target selection and characteristics of the DESI Peculiar Velocity Survey, the largest survey of peculiar velocities (PVs) using both the fundamental plane (FP) of galaxies and the Tully-Fisher (TF) relationship planned to date. We detail how we identify suitable early-type galaxies for the FP and suitable late-type galaxies for the TF relation using the photometric data provided by the DESI Legacy Imaging Survey DR9. Subsequently, we provide targets for 373 533 early-type galaxies and 118 637 late-type galaxies within the DESI 5-year footprint. We validate these photometric selections using existing morphological classifications. Furthermore, we demonstrate using survey validation data that DESI is able to measure the spectroscopic properties to sufficient precision to obtain PVs for our targets. Based on realistic DESI fiber assignment simulations and spectroscopic success rates, we predict the final DESI Peculiar Velocity Survey will obtain $\sim$133 000 FP-based and $\sim$53 000 TF-based PV measurements over an area of 14 000 $\mathrm{deg^{2}}$. Each of these components will be a factor of 4--5 larger than other recent samples. We forecast the ability of using these data to measure the clustering of galaxy positions and peculiar velocities from the combined DESI PV and Bright Galaxy Surveys (BGS), which allows for cancellation of cosmic variance at low redshifts. With these forecasts, we anticipate a $4\%$ statistical measurement on the growth rate of structure at $z<0.15$. This is over two times better than achievable with redshifts from the BGS alone. The combined DESI PV and Bright Galaxy surveys will enable the most precise tests to date of the time and scale dependence of large-scale structure growth at $z<0.15$.

Thermal X-ray spectra from supernova remnants (SNRs) are dominated by a number of line emission from various elements. Resolving the individual lines is critically important for a variety of scientific topics such as diagnosing high-temperature and low-density non-equilibrium plasmas, identifying spectral features like charge exchange and resonance line scattering, revealing kinematics and elemental abundances of SN ejecta and the circumstellar medium, and studying the interstellar medium or planets' atmospheres from extinction features seen in X-ray spectra of very bright SNRs. This chapter reviews high-resolution X-ray spectroscopy of SNRs obtained so far. Most results were derived with dispersive spectrometers aboard Einstein, Chandra, and XMM-Newton satellites. Because these dispersive spectrometers were slitless, one has to select small objects with angular sizes less than a few arcminutes to successfully perform high-resolution spectroscopy. Despite this limitation, the three satellites delivered fruitful scientific results in the last few decades. Arrays of low-temperature microcalorimeters offer excellent opportunities for high-resolution X-ray spectroscopy of SNRs, as they are non-dispersive spectrometers that work for largely extended sources as well as point-like sources. The microcalorimeter aboard the Hitomi satellite already delivered pioneering results during its short lifetime. The upcoming X-Ray Imaging and Spectroscopy Mission, which is a recovery mission of Hitomi, will truly open the new discovery window to high-resolution X-ray spectroscopy of SNRs.

A short X-ray burst was observed from the radio-loud magnetar 1E 1547.0--5408 in April 2022. Unusually however, the source stopped showing radio pulsations $\gtrsim 3\,$weeks \emph{prior} to the burst. After recovery, radio timing revealed that the object had also undergone a modest glitch. A model for the overall event is constructed where an initially mild perturbation adjusts the magnetic geometry near the polar caps, leading to shallow fractures. Crustal ejecta or particles leaking from a pair-plasma fireball pollute the magnetospheric gaps, shutting off the pulsar mechanism, but the energy release is not yet large enough to noticeably enhance the X-ray flux. This perturbation gradually ramps, eventuating in a large-scale energy redistribution which fuels the burst. The star's mass quadrupole moment changes in tandem, issuing a glitch. Some quantitative estimates for the magnetic reconfiguration under this interpretation are provided, based on a quasi-static model where the fluid evolves through a sequence of hydromagnetic equilibria.

Future space observatories achieve detection of gravitational waves by interferometric measurements of a carrier phase, allowing to determine relative distance changes, in combination with an absolute distance measurement based on the transmission of pseudo-random noise chip sequences. In addition, usage of direct-sequence spread spectrum modulation enables data transmission. Hereafter, we report on the findings of a performance evaluation of planned receiver architectures, performing phase and distance readout sequentially. An analytical model is presented identifying the power spectral density of the chip modulation at frequencies within the measurement bandwidth as the main driver for phase noise. This model, verified by numerical simulations, excludes binary phase-shift keying modulations for missions requiring pico-meter noise levels at the phase readout, while binary offset carrier modulation, where most of the power has been shifted outside the measurement bandwidth, exhibits superior phase measurement performance. Ranging analyses of the delay-locked loop reveal strong distortion of the pulse shape due to the preceding phase tracking introducing ranging bias variations. Numerical simulations show that these variations, however, originating from data transitions, are compensated by the delay tracking loop and are thus not considered critical, irrespective of the modulation type.

In order to investigate the impact of radio jets on the interstellar medium (ISM) of galaxies hosting active galactic nuclei (AGN), we present subarcsecond resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations of the Teacup galaxy. This is a nearby ($D_{\rm L}$=388 Mpc) radio-quiet type-2 quasar (QSO2) with a compact radio jet ($P_{\rm jet}\approx$10$^{43}$ erg s$^{-1}$) that subtends a small angle from the molecular gas disc. Enhanced emission line widths perpendicular to the jet orientation have been reported for several nearby AGN for the ionised gas. For the molecular gas in the Teacup, not only do we find this enhancement in the velocity dispersion but also a higher brightness temperature ratio (T32/T21) perpendicular to the radio jet compared to the ratios found in the galaxy disc. Our results and the comparison with simulations suggest that the radio jet is compressing and accelerating the molecular gas, and driving a lateral outflow that shows enhanced velocity dispersion and higher gas excitation. These results provide further evidence that the coupling between the jet and the ISM is relevant to AGN feedback even in the case of radio-quiet galaxies.

Combining the `time-delay distance' ($D_{\Delta t}$) measurements from galaxy lenses and other distance indicators provides model-independent determinations of the Hubble constant ($H_0$) and spatial curvature ($\Omega_{K,0}$), only based on the validity of the Friedmann-Lema\^itre-Robertson-Walker (FLRW) metric and geometrical optics. To take the full merit of combining $D_{\Delta t}$ measurements in constraining $H_0$, we use gamma-ray burst (GRB) distances to extend the redshift coverage of lensing systems much higher than that of Type Ia Supernovae (SNe Ia) and even higher than quasars, whilst the general cosmography with a curvature component is implemented for the GRB distance parametrizations. Combining Lensing+GRB yields $H_0=71.5^{+4.4}_{-3.0}$~km s$^{-1}$Mpc$^{-1}$ and $\Omega_{K,0} = -0.07^{+0.13}_{-0.06}$ (1$\sigma$). A flat-universe prior gives slightly an improved $H_0 = 70.9^{+4.2}_{-2.9}$~km s$^{-1}$Mpc$^{-1}$. When combining Lensing+GRB+SN Ia, the error bar $\Delta H_0$ falls by 25\%, whereas $\Omega_{K,0}$ is not improved due to the degeneracy between SN Ia absolute magnitude, $M_B$, and $H_0$ along with the mismatch between the SN Ia and GRB Hubble diagrams at $z\gtrsim 1.4$. Future increment of GRB observations can help to moderately eliminate the $M_B-H_0$ degeneracy in SN Ia distances and ameliorate the restrictions on cosmographic parameters along with $\Omega_{K,0}$ when combining Lensing+SN Ia+GRB. We conclude that there is no evidence of significant deviation from a (an) flat (accelerating) universe and $H_0$ is currently determined at 3\% precision. The measurements show great potential to arbitrate the $H_0$ tension between the local distance ladder and cosmic microwave background measurements and provide a relevant consistency test of the FLRW metric.

We present radio coverage of the 2019 outburst of the accreting millisecond X-ray pulsar SAX J1808.4-3658, obtained with MeerKAT. We compare these data to contemporaneous X-ray and optical measurements in order to investigate the coupling between accretion and jet formation in this system, while the optical lightcurve provides greater detail of the outburst. The reflaring activity following the main outburst peak was associated with a radio re-brightening, indicating a strengthening of the jet in this phase of the outburst. We place quasi-simultaneous radio and X-ray measurements on the global radio:X-ray plane for X-ray binaries, and show they reside in the same region of luminosity space as previous outburst measurements, but significantly refine the correlation for this source. We also present upper limits on the radio emission from the accreting millisecond X-ray pulsar MAXI J0911-655 and the transitional Z/Atoll-type transient XTE J1701-462. In the latter source we also confirm that nearby large-scale structures reported in previous radio observations of the source are persistent over a period of ~15 years, and so are almost certainly background radio galaxies and not associated with the X-ray transient.

The Neil Gehrels Swift Observatory (Swift) has been in operation for 18 years. The Ultra-Violet/Optical Telescope (UVOT) onboard Swift was designed to capture the earliest optical/UV emission from gamma-ray bursts (GRBs), spanning the first few minutes to days after the prompt gamma-ray emission. In this article, I provide an overview of the long GRBs (whose prompt gamma-ray duration is >2 s) observed by the Swift/UVOT and review the major discoveries that have been achieved by Swift/UVOT over the last 18 years. I discuss where improvements have been made to our knowledge and understanding of the optical/UV emission, particularly the early optical/UV afterglow.

The Galactic interstellar medium abounds in low-velocity shocks with velocities less than, say, about 70 km/s. Some are descendants of higher velocity shocks, while others start off at low velocity (e.g., stellar bow shocks, intermediate velocity clouds, spiral density waves). Low-velocity shocks cool primarily via Ly-alpha, two-photon continuum, optical recombination lines (e.g., H-alpha), free-bound emission, free-free emission and forbidden lines of metals. The dark far-ultraviolet (FUV) sky, aided by the fact that the two-photon continuum peaks at 1400 angstroms, makes the FUV band an ideal tracer of low-velocity shocks. Recent GALEX FUV images reaffirm this expectation, discovering faint and large interstellar structure in old supernova remnants and thin arcs stretching across the sky. Interstellar bow shocks are expected from fast stars from the Galactic disk passing through the numerous gas clouds in the local interstellar medium within 15 pc of the Sun. Using the best atomic data available to date, we present convenient fitting formulae for yields of Ly$\alpha$, two-photon continuum and H$\alpha$ for pure hydrogen plasma in the temperature range of 10^4 K to 10^5 K. The formulae presented here can be readily incorporated into time-dependent cooling models as well as collisional ionization equilibrium models.

Given their strong stellar winds, Wolf-Rayet (WR) stars exhibit emission line spectra that are predominantly formed in expanding atmospheric layers. The description of the wind velocity field $v(r)$ is therefore a crucial ingredient in the spectral analysis of WR stars, possibly influencing the determination of stellar parameters. In view of this, we perform a systematic study by simulating a sequence of WR-star spectra for different temperatures and mass-loss rates using $\beta$-type laws with $0.5\leq\beta\leq 20$. We quantify the impact of varying $v(r)$ by analysing diagnostic lines and spectral classifications of emergent model spectra computed with the Potsdam Wolf-Rayet (PoWR) code. We additionally cross-check these models with hydrodynamically consistent -- hydro -- model atmospheres. Our analysis confirms that the choice of the $\beta$-exponent has a strong impact on WR-star spectra, affecting line widths, line strengths and line profiles. In some parameter regimes, the entire range of WR subtypes could be covered. Comparison with observed WR stars and hydro models revealed that values of $\beta \gtrsim 8$ are unlikely to be realized in nature, but a range of $\beta$-values needs to be considered in spectral analysis. UV spectroscopy is crucial here to avoid an underestimation of the terminal velocity $v_\infty$. Neither single- nor double-$\beta$ descriptions yield an acceptable approximation of the inner wind when compared to hydro models. Instead, we find temperature shifts to lower $T_{2/3}$ when employing a hydro model. Additionally, there are further hints that round-lined profiles seen in several early WN stars are an effect from non-$\beta$ velocity laws.

High resolution spectroscopy has allowed for unprecedented levels of atmospheric characterization, especially for the hottest gas giant exoplanets known as ultrahot Jupiters (UHJs). High-resolution spectra are sensitive to 3D effects, making complex 3D atmospheric models important for interpreting data. Moreover, these planets are expected to host magnetic fields that will shape their resulting atmospheric circulation patterns, but little modeling work has been done to investigate these effects. In this paper, we generate high-resolution transmission spectra from General Circulation Models for the canonical UHJ WASP-76b with three different magnetic treatments in order to understand the influence of magnetic forces on the circulation. In general, spectra from all models have increasingly blueshifted net Doppler shifts as transit progresses, but we find that the differing temperature and wind fields in the upper atmospheres of these models result in measurable differences. We find that magnetic effects may be contributing to the unusual trends previously seen in transmission for this planet. Our $B=3$ Gauss active drag model in particular shows unique trends not found in the models with simpler or no magnetic effects. The net Doppler shifts are additionally influenced by the dominant opacity sources in each wavelength range considered, as each species probes different regions of the atmosphere and are sensitive to spatial differences in the circulation. This work highlights the ongoing need for models of planets in this temperature regime to consider both 3D and magnetic effects when interpreting high resolution transmission spectra.

IC 10 X-1 is an eclipsing high mass X-ray binary (HMXB) containing a stellar-mass black hole (BH) and a Wolf-Rayet (WR) donor star with an orbital period of P = 34.9 hr. This binary belongs to a group of systems that can be the progenitors of gravitational wave sources, hence understanding the dynamics of systems such as IC 10 X-1 is of paramount importance. The prominent He II 4686 emission line (previously used in mass estimates of the BH) is out of phase with the X-ray eclipse, suggesting that this line originates somewhere in the ionized wind of the WR star or in the accretion disk. We obtained 52 spectra from the GEMINI/GMOS archive, observed between 2001 and 2019. We analyzed the spectra both individually, and after binning them by orbital phase to improve the signal-to-noise ratio. The RV curve from the stacked data is similar to historical results, indicating the overall parameters of the binary have remained constant. However, the He II line profile shows a correlation with the X-ray hardness-ratio values, also, we report a pronounced skewness of the line-profile, and the skewness varies with the orbital phase. These results support a paradigm wherein the He II line tracks structures in the stellar wind that are produced by interactions with the BH's ionizing radiation and the accretion flow. We compare the observable signatures of two alternative hypotheses proposed in the literature: wind irradiation plus shadowing, and accretion disk hotspot; and we explore how the line-profile variations fit into each of these models.

Hotspots of radio galaxies are regions of shock-driven particle acceleration. Multiple hotspots have long been identified as potential indicators of jet movement or precession. Two frequent explanations describe a secondary hotspot as either the location of a prior jet termination point, or a deflected backflow-driven shock: the so-called Dentist's Drill and Splatter Spot models. We created high-resolution simulations of precessing jets with a range of parameters. In addition to the existing mechanisms, our results show three additional mechanisms for multiple hotspot formation: (1) the splitting of a large terminal hotspots into passive and active components; (2) jet stream splitting resulting in two active hotspots; (3) dynamic multiple hotspot complexes that form as a result of jet termination in a turbulent cocoon, linked here to rapid precession. We show that these distinct types of multiple hotspots are difficult to differentiate in synthetic radio maps, particularly hotspot complexes which can easily be mistaken for the jet itself. We discuss the implication for hypothesised binary supermassive black hole systems where jet precession is a key component of the morphology, and show a selection of potential precession candidates found using the LOFAR Two-Metre Sky Survey Data Release 2 (LoTSS DR2).

GRB 221009A has been referred to as the Brightest Of All Time (the BOAT). We investigate the veracity of this statement by comparing it with a half century of prompt gamma-ray burst observations. This burst is the brightest ever detected by the measures of peak flux and fluence. Unexpectedly, GRB 221009A has the highest isotropic-equivalent total energy ever identified, while the peak luminosity is at the $\sim99$th percentile of the known distribution. We explore how such a burst can be powered and discuss potential implications for ultra-long and high-redshift gamma-ray bursts. By geometric extrapolation of the total fluence and peak flux distributions GRB 221009A appears to be a once in 10,000 year event. Thus, while it almost certainly not the BOAT over all of cosmic history, it may be the brightest gamma-ray burst since human civilization began.

Very-high-energy (VHE) $\gamma$-rays ($\gtrsim 0.1\rm~TeV$) and neutrinos are crucial for identifying accelerators of ultrahigh-energy cosmic rays (UHECRs), but this is challenging especially for UHECR nuclei. In this work, we develop a numerical code to solve the transport equation for UHECRs and their secondaries, where both nuclear and electromagnetic cascades are taken into account self-consistently, considering steady UHECR accelerators such as radio galaxies. In particular, we focus on Centaurus A, which has been proposed as one of the most promising UHECR sources in the local universe. Motivated by observations of extended VHE $\gamma$-ray emission from its kiloparsec-scale jet by the H.E.S.S. telescope, we study interactions between UHECRs accelerated in the large-scale jet and various target photon fields including blazar-like beamed core emission, and present a quantitative study on VHE $\gamma$-ray signatures of UHECR nuclei, including the photodisintegration and Bethe-Heitler pair-production processes. We show that VHE $\gamma$-rays from UHECR nuclei could be detected by the ground-based $\gamma$-ray telescopes given that the dominant composition of UHECRs consists of intermediate-mass (such as oxygen) nuclei

We study the dynamics of a collisionless kinetic gas whose particles follow future-directed timelike and spatially bound geodesics in the exterior of a sub-extremal Kerr black hole spacetime. Based on the use of generalized action-angle variables, we analyze the large time asymptotic behavior of macroscopic observables associated with the gas. We show that, as long as the fundamental frequencies of the system satisfy a suitable non-degeneracy condition, these macroscopic observables converge in time to the corresponding observables determined from an averaged distribution function. In particular, this implies that the final state is characterized by a distribution function which is invariant with respect to the full symmetry group of the system, that is, it is stationary, axisymmetric and Poisson-commutes with the integral of motion associated with the Carter constant. As a corollary of our result, we demonstrate the validity of the strong Jeans theorem in our setting, stating that the distribution function belonging to a stationary state must be a function which is independent of the generalized angle variables. An analogous theorem in which the assumption of stationarity is replaced with the requirement of invariance with respect to the Carter flow is also proven. Finally, we prove that the aforementioned non-degeneracy condition holds. This is achieved by providing suitable asymptotic expansions for the energy and Carter constant in terms of action variables for orbits having sufficiently large radii, and by exploiting the analytic dependency of the fundamental frequencies on the integrals of motion.

We argue that, in the same way that in a black hole space-time VECROs will form in order to cancel the gravitational effects of a collapsing mass shell and prevent the formation of a singularity, in a contracting universe a gas of VECROs will form to hold up the contraction, prevent a Big Crunch singularity, and lead to a nonsingular cosmological bounce.

As the European maritime powers expanded their reach beyond north Atlantic coastal waters to distant lands as far away as the East Indies, access to a practical means of maritime navigation in the southern hemisphere became imperative. The first few voyages undertaken by the Dutch East India Company and its predecessor explicitly aimed at compiling star charts and constellations that were only visible south of the Equator, as practical navigation aids. The oldest known star atlas of southern constellations was published in 1603 by Frederick de Houtman. Controversies have plagued de Houtman's astronomical credentials from their inception, however, with contemporaries variously attributing the early southern star charts to Pieter Dirkszoon Keyser, de Houtman, or even to their tutor Petrus Plancius. The balance of available evidence suggests that Keyser initially led the astronomical observing campaign, ably assisted by de Houtman. Upon Keyser's untimely death, de Houtman embraced a leading role in compiling astronomical observations for maritime navigation purposes, whereas Plancius most probably led the delineation of the 12 new southern constellations that soon became part and parcel of the nautical consciousness.

The hypothesis that the mass of BHs increases with time according to the same law as the volume of the part of the Universe containing it and therefore the population of BHs is similar to dark energy in its action was recently proposed. We demonstrate the reasons why it cannot be accepted, even if all the assumptions on which this hypothesis is based are considered true.

This is the first in a series of papers where we study the dynamics of a bubble wall beyond usual approximations, such as the assumptions of spherical bubbles and infinitely thin walls. In this paper, we consider a vacuum phase transition. Thus, we describe a bubble as a configuration of a scalar field whose equation of motion depends only on the effective potential. The thin-wall approximation allows obtaining both an effective equation of motion for the wall position and a simplified equation for the field profile inside the wall. Several different assumptions are involved in this approximation. We discuss the conditions for the validity of each of them. In particular, the minima of the effective potential must have approximately the same energy, and we discuss the correct implementation of this approximation. We consider different improvements to the basic thin-wall approximation, such as an iterative method for finding the wall profile and a perturbative calculation in powers of the wall width. We calculate the leading-order corrections. Besides, we derive an equation of motion for the wall without any assumptions about its shape. We present a suitable method to describe arbitrarily deformed walls from the spherical shape. We consider concrete examples and compare our approximations with numerical solutions. In subsequent papers, we shall consider higher-order finite-width corrections, and we shall take into account the presence of the fluid.

We developed a fast and modular deep learning algorithm to search for lookalike signals of interest in radio spectrogram data. First, we trained an autoencoder on filtered data returned by an energy detection algorithm. We then adapted a positional embedding layer from classical Transformer architecture to a frequency-based embedding. Next we used the encoder component of the autoencoder to extract features from small (~ 715,Hz with a resolution of 2.79Hz per frequency bin) windows in the radio spectrogram. We used our algorithm to conduct a search for a given query (encoded signal of interest) on a set of signals (encoded features of searched items) to produce the top candidates with similar features. We successfully demonstrate that the algorithm retrieves signals with similar appearance, given only the original radio spectrogram data.

We study the momentum-space entanglement between the sub- and super-Hubble modes of a spectator scalar field, with a cubic $\lambda \phi^3$ interaction, in de Sitter space. Momentum-space entanglement has some universal properties for any interacting quantum field theory, and we examine them for this specific curved background using the Hubble scale as a natural delimiter to define UV/IR separation. We show that there are several new subtleties when generalising flat space results due to having a time-dependent interaction term and a non-trivial vacuum state. Our main finding is that the momentum-space entanglement entropy in de Sitter space grows very rapidly, supporting previous similar results for cosmological perturbations [1], which leads to interesting new questions.

The radio signal transmitted by the Mars Express (MEX) spacecraft was observed regularly between the years 2013-2020 at X-band (8.42 GHz) using the European Very Long Baseline Interferometry (EVN) network and University of Tasmania's telescopes. We present a method to describe the solar wind parameters by quantifying the effects of plasma on our radio signal. In doing so, we identify all the uncompensated effects on the radio signal and see which coronal processes drive them. From a technical standpoint, quantifying the effect of the plasma on the radio signal helps phase referencing for precision spacecraft tracking. The phase fluctuation of the signal was determined for Mars' orbit for solar elongation angles from 0 - 180 deg. The calculated phase residuals allow determination of the phase power spectrum. The total electron content (TEC) of the solar plasma along the line of sight is calculated by removing effects from mechanical and ionospheric noises. The spectral index was determined as $-2.43 \pm 0.11$ which is in agreement with Kolomogorov's turbulence. The theoretical models are consistent with observations at lower solar elongations however at higher solar elongation ($>$160 deg) we see the observed values to be higher. This can be caused when the uplink and downlink signals are positively correlated as a result of passing through identical plasma sheets.

Compact boson star binaries are hypothetical sources for ground-based and space gravitational-wave detectors. Their signal would be a messenger for novel fundamental fields and could shed light on the dark matter. In this work, we further develop our analysis in Phys. Rev. D 102, 083002 (2020), aimed at constraining the properties of these objects with future observations. We use a coherent waveform template for the inspiral stage of boson star binaries with large quartic self interactions, including tidal deformability and the nonlinear dependence of the quadrupole moments on the spin in terms of the fundamental couplings of the scalar field theory. Performing a Bayesian analysis, we investigate the ability of a third-generation gravitational-wave detector such as the Einstein Telescope to distinguish these exotic sources from black holes and infer constraints on the fundamental couplings of the model.