Papers for Tuesday, Feb 14 2023

This paper corresponds to an invited oral contribution to the session 5A organised by the IAU inter-commission B2-B5 working group (WG) "Laboratory Astrophysics Data Compilation, Validation and Standardization: from the Laboratory to FAIR usage in the Astronomical Community" at the IAU 2022 General Assembly (GA). This WG provides a platform where to discuss the Findability, Accessibility,Interoperability, Reuse (FAIR) usage of laboratory Atomic and Molecular (A&M) data in astronomy and astrophysics. A&M data play a key role in the understanding of the physics and chemistry of processes in several research topics, including planetary science and interdisciplinary research in particular the atmospheres of planets and planetary explorations, etc. Databases, compilation of spectroscopic parameters, and facility tools are used by computer codes to interpret spectroscopic observations and simulate them. In this talk I presented existing A&M databases of interest to the planetary community focusing on access, organisation, infrastructures, limitations and issues, etc.

This is the second data release (DR2) of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 21 new systems with observations completed before March 2022. For each system, we release four one-dimensional (1D) spectra to the ESO Science Archive, one spectrum for each arm of the X-Shooter spectrograph. In this paper, we focus on the line-of-sight velocity distribution, measuring integrated stellar velocity dispersions from the spectra, and assessing their robustness and the associated uncertainties. For each of the 21 new systems, we systematically investigated the effect of the parameters and set-ups of the full spectral fitting on the stellar velocity dispersion ($\sigma$) measurements. In particular, we tested how $\sigma$ changes when several parameters of the fit as well as the resolution and spectral coverage of the input spectra are varied. We found that the effect that causes the largest systematic uncertainties on $\sigma$ is the wavelength range used for the fit, especially for spectra with a lower signal-to-noise ratio (S/N $\leq$ 30). When using blue wavelengths (UVB arm) one generally underestimates the velocity dispersion (by $\sim$15 km/s). The values obtained from the near-IR (NIR) arm present a larger scatter because the quality of the spectra is lower. We finally compared our results with those in literature, finding a very good agreement overall. Joining results obtained in DR1 with those presented here, INSPIRE contains 40 ultra-compact massive galaxies, corresponding to 75% of the whole survey. By plotting these systems in a stellar mass-velocity dispersion diagram, we identify at least four highly reliable relic candidates among the new systems. Their velocity dispersion is larger than that of normal-sized galaxies of similar stellar mass.

RR Lyrae variables are excellent population II distance indicators thanks to their well-defined period-luminosity relations (PLRs) at infrared wavelengths. We present results of near-infrared (NIR) monitoring of Galactic globular clusters to empirically quantify the metallicity dependence of NIR PLRs for RR Lyrae variables. Our sample includes homogeneous, accurate, and precise photometric data for 964 RR Lyrae variables in 11 globular clusters covering a large metallicity range ($\Delta\textrm{[Fe/H]}\sim2$~dex). We derive $JHK_s$ band period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations anchored using 346 Milky Way field RR Lyrae stars with {\it Gaia} parallaxes, and simultaneously solved for independent distances to globular clusters. We find a significant metallicity dependence of $\sim0.2$~mag/dex in $JHK_s$ band PLZ and PWZ relations for RR Lyrae stars independent of the adopted metallicity scale. The metallicity coefficients and the zero-points of the empirical PLZ and PWZ relations are in excellent agreement with the predictions from the horizontal branch evolution and pulsation models. Furthermore, RR Lyrae based distances to our sample of globular clusters are also statistically consistent with other independent measurements in the literature. Our recommended empirical $JHK_s$ band PLZ relations are also provided for RR Lyrae based distance measurements.

Gamma-ray bursts (GRBs) have long been considered a possible source of high-energy neutrinos. While no correlations have yet been detected between high-energy neutrinos and GRBs, the recent observation of GRB 221009A - the brightest GRB observed by Fermi-GBM to date and the first one to be observed above an energy of 10 TeV - provides a unique opportunity to test for hadronic emission. In this paper, we leverage the wide energy range of the IceCube Neutrino Observatory to search for neutrinos from GRB 221009A. We find no significant deviation from background expectation across event samples ranging from MeV to PeV energies, placing stringent upper limits on the neutrino emission from this source.

We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $\Omega_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider the three dimensional case, with misaligned magnetic and rotation axes. We discuss applications of our results to the spindown of normal stars, highly-irradiated exoplanets, and to nascent highly-magnetic and rapidly-rotating neutron stars born in massive star core collapse.

We study the bar pattern speeds and corotation radii of 225 barred galaxies, using IFU data from MaNGA and the Tremaine-Weinberg method. Our sample, which is divided between strongly and weakly barred galaxies identified via Galaxy Zoo, is the largest that this method has been applied to. We find lower pattern speeds for strongly barred galaxies than for weakly barred galaxies. As simulations show that the pattern speed decreases as the bar exchanges angular momentum with its host, these results suggest that strong bars are more evolved than weak bars. Interestingly, the corotation radius is not different between weakly and strongly barred galaxies, despite being proportional to bar length. We also find that the corotation radius is significantly different between quenching and star forming galaxies. Additionally, we find that strongly barred galaxies have significantly lower values for R, the ratio between the corotation radius and the bar radius, than weakly barred galaxies, despite a big overlap in both distributions. This ratio classifies bars into ultrafast bars (R < 1.0; 11% of our sample), fast bars (1.0 < R < 1.4; 27%) and slow bars (R > 1.4; 62%). Simulations show that R is correlated with the bar formation mechanism, so our results suggest that strong bars are more likely to be formed by different mechanisms than weak bars. Finally, we find a lower fraction of ultrafast bars than most other studies, which decreases the recently claimed tension with {\Lambda}CDM. However, the median value of R is still lower than what is predicted by simulations.

We present a comprehensive catalog of observations and stellar population properties for 23 highly secure host galaxies of fast radio bursts (FRBs). Our sample comprises six repeating FRBs and 17 apparent non-repeaters. We present 82 new photometric and eight new spectroscopic observations of these hosts. Using stellar population synthesis modeling and employing non-parametric star formation histories (SFHs), we find that FRB hosts have a median stellar mass of $\approx 10^{9.8}\,M_{\odot}$, mass-weighted age of $\approx 5.2$ Gyr, and ongoing star formation rate $\approx 1.3\,M_{\odot}$ yr$^{-1}$ but span wide ranges in all properties. Classifying the hosts by degree of star formation, we find that 91\% (21/23 hosts) are star-forming, one is transitioning, and one is quiescent. The majority trace the star-forming main sequence of galaxies, but at least two FRBs originate in less active environments, both of which are apparent non-repeaters. Across all modeled properties, we find no statistically significant distinction between the hosts of repeaters and non-repeaters. However, the hosts of repeating FRBs generally extend to lower stellar masses, and the hosts of non-repeaters arise in more optically luminous galaxies. Further, the four galaxies with the most clear and prolonged rises in their SFHs all host repeating FRBs, demonstrating heightened star formation activity in the last $\lesssim 100$ Myr. Our results support the young magnetar model for most FRBs in which their progenitors formed through core-collapse supernovae, but the presence of some FRBs in less active environments suggests that a fraction form through more delayed channels.

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5<z<4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam (m=30.6-30.9, 5sigma) to discover galaxies to z>12, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products.

Interstellar dust captures a significant fraction of elements heavier than helium in the solid state and is an indispensable component both in theory and observations of galaxy evolution. Dust emission is generally the primary coolant of the interstellar medium (ISM) and facilitates the gravitational collapse and fragmentation of gas clouds from which stars form, while altering the emission spectrum of galaxies from ultraviolet (UV) to far-infrared wavelengths through the reprocessing of starlight. However, the astrophysical origin of various types of dust grains remains an open question, especially in the early Universe. Here we report direct evidence for the presence of carbonaceous grains from the detection of the broad UV absorption feature around $2175 \, \mathring{\rm A}$ in deep near-infrared spectra of galaxies up to the first billion years of cosmic time, at a redshift ($z$) of $\sim 7$. This dust attenuation feature has previously only been observed spectroscopically in older, more evolved galaxies at redshifts of $z < 3$. The carbonaceous grains giving rise to this feature are often thought to be produced on timescales of hundreds of millions of years by asymptotic giant branch (AGB) stars. Our results suggest a more rapid production scenario, likely in supernova (SN) ejecta.

We study how well void-finding algorithms identify cosmic void regions and whether we can quantitatively and qualitatively describe their biases by comparing the voids they find with dynamical information from the underlying matter distribution. Using the ORIGAMI algorithm to determine the number of dimensions along which dark matter particles have undergone shell-crossing (crossing number) in $N$-body simulations from the AbacusSummit simulation suite, we identify dark matter particles which have undergone no shell crossing as belonging to voids. We then find voids in the corresponding halo distribution using two different void-finding algorithms: VoidFinder and V$^2$, a ZOBOV-based algorithm. The resulting void catalogs are compared to the distribution of dark matter particles to examine how their crossing numbers depend on void proximity. While both algorithms' voids have a similar distribution of crossing numbers near their centers, we find that beyond 0.25 times the effective void radius, voids found by VoidFinder exhibit a stronger preference for particles with low crossing numbers than those found by V$^2$. We examine two possible methods of mitigating this difference in efficacy between the algorithms. While we are able to partially mitigate the ineffectiveness of V$^2$ by using distance from the void edge as a measure of centrality, we conclude that VoidFinder more reliably identifies dynamically-distinct regions of low crossing number.

Constraints on the formation and evolution of the Milky Way Galaxy require multi-dimensional measurements of kinematics, abundances, and ages for a large population of stars. Ages for luminous giants, which can be seen to large distances, are an essential component of studies of the Milky Way, but they are traditionally very difficult to estimate precisely for a large dataset and often require careful analysis on a star-by-star basis in asteroseismology. Because spectra are easier to obtain for large samples, being able to determine precise ages from spectra allows for large age samples to be constructed, but spectroscopic ages are often imprecise and contaminated by abundance correlations. Here we present an application of a variational encoder-decoder on cross-domain astronomical data to solve these issues. The model is trained on pairs of stars observed by APOGEE and Kepler in order to reduce the dimensionality of the APOGEE spectra in a latent space while removing abundance information. The low dimensional latent representation of these spectra can then be trained to predict age with just $\sim$ 1,000 precise seismic ages. We demonstrate that this model produces more precise spectroscopic ages ($\sim$ 22% overall, $\sim$ 11% for red-clump stars) than previous data-driven spectroscopic ages while being less contaminated by abundance information (in particular, our ages do not depend on [$\alpha$/M]). We create a public age catalog for the APOGEE DR17 data set and use it to map the age distribution and the age-[Fe/H]-[$\alpha$/M] distribution across the radial range of the Galactic disk.

Radio spectroscopy provides a unique inspection perspective for solar and space weather research, which can reveal the plasma and energetic electron information in the solar corona and inner heliosphere. However, Radio-Frequency Interference (RFI) from human activities affects sensitive radio telescopes, and significantly affects the quality of observation. Thus, RFI detection and mitigation for the observations is necessary to obtain high quality, science-ready data. The flagging of RFI is particularly challenging for the solar and space weather observations at low frequency, because the solar radio bursts can be brighter than the RFI, and may show similar temporal behavior. In this work, we investigate RFI flagging methods for solar and space weather observations, including a strategy for AOFlagger, and a novel method that makes use of a morphology convolution. These algorithms can effectively flag RFI while preserving solar radio bursts.

We recently introduced a new method for simulating collisional gravitational N-body systems with approximately linear time scaling with $N$, based on the Multi-Particle Collision (MPC) scheme, previously applied in Plasma Physics. We simulate globular clusters with a realistic number of stellar particles (at least up to several times $10^6$) on a standard workstation. We simulate clusters hosting an intermediate mass black hole (IMBH), probing a broad range of BH-cluster and BH-average-star mass ratios, unrestricted by the computational constraints affecting direct N-body codes. We use either single mass models or models with a Salpeter mass function, with the IMBH initially sitting at the centre. The force exerted by and on the IMBH is evaluated with a direct scheme. We measure the evolution of the Lagrangian radii and core density and velocity dispersion over time. In addition, we study the evolution of the velocity anisotropy profiles. We find that models with an IMBH undergo core collapse at earlier times, the larger the IMBH mass the shallower, with an approximately constant central density at core collapse. The presence of an IMBH tends to lower the central velocity dispersion. These results hold independently of the mass function. For the models with Salpeter MF we observe that equipartition of kinetic energies is never achieved. Orbital anisotropy at large radii appears driven by energetic escapers on radial orbits. We measure the wander radius. Among the results we obtained, which mostly confirm or extend previously known trends that had been established over the range of parameters accessible to direct N-body simulations, we underline that the leptokurtic nature of the IMBH wander radius distribution might lead to IMBHs presenting as off-center more frequently than expected, with implications on observational IMBH detection.

The goal of this work is to search for Active Galactic Nuclei (AGN) in the Galactic disc at very low latitudes with |b| $<$ 2$^\circ$. For this, we studied the five sources from the VVV near-infrared galaxy catalogue that have also WISE counterparts and present variability in the VIrac VAriable Classification Ensemble (VIVACE) catalogue. In the near-infrared colour-colour diagrams, these objects have in general redder colours compared to the rest of the sources in the field. In the mid-infrared ones, they are located in the AGN region, however there is a source that presents the highest interstellar extinction and different mid-IR colours to be a young stellar object (YSO). We also studied the source variability using two different statistical methods. The fractional variability amplitude $\sigma_{rms}$ ranges from 12.6 to 33.8, being in concordance with previous results found for type-1 AGNs. The slopes of the light curves are in the range (2.6$-$4.7) $\times 10^{-4}$ mag/day, also in agreement with results reported on quasars variability. The combination of all these results suggest that four galaxies are type-1 AGN candidates whereas the fifth source likely a YSO candidate.

We systematically inspect the microlensing data acquired by the KMTNet survey during the previous seasons in order to find anomalous lensing events for which the anomalies in the lensing light curves cannot be explained by the usual binary-lens or binary-source interpretations. From the inspection, we find that interpreting the three lensing events OGLE-2018-BLG-0584, KMT-2018-BLG-2119, and KMT-2021-BLG-1122 requires four-body (lens+source) models, in which either both the lens and source are binaries (2L2S event) or the lens is a triple system (3L1S event). Following the analyses of the 2L2S events presented in \citet{Han2023}, here we present the 3L1S analysis of the KMT-2021-BLG-1122. It is found that the lens of the event KMT-2021-BLG-1122 is composed of three masses, in which the projected separations (normalized to the angular Einstein radius) and mass ratios between the lens companions and the primary are $(s_2, q_2)\sim (1.4, 0.53)$ and $(s_3, q_3) \sim (1.6, 0.24)$. By conducting a Bayesian analysis, we estimate that the masses of the individual lens components are $(M_1, M_2, M_3)\sim (0.47\,M_\odot, 0.24\,M_\odot, 0.11\,M_\odot)$. The companions are separated in projection from the primary by $(a_{\perp,2}, a_{\perp,3})\sim (3.5, 4.0)$~AU. The lens of KMT-2018-BLG-2119 is the first triple stellar system detected via microlensing.

Using the Effective Field Theory (EFT) framework of single field inflation, we investigate the possibility of the formation of Primordial Black Holes (PBHs) in the Slow Roll (SR) to Ultra Slow Roll (USR) transition. We demonstrate that, due to one loop correction to the power spectrum, causality is violated ($c_s>1$) for the mass range of PBHs, $M_{\rm PBH}>10^{2}{\rm gm}$ created during the said transition. We find that non-canonical features with $c_s<1$ worsen the predictions of canonical framework of single field inflation.

In this study, we reported the results of high-resolution (0.14 arcsec) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 225 GHz dust continuum and CO molecular emission lines from the transitional disk around SY Cha. Our high-resolution observations clearly revealed the inner cavity and the central point source for the first time. The radial profile of the ring can be approximated by a bright narrow ring superimposed on a fainter wide ring. Furthermore, we found that there is a weak azimuthal asymmetry in dust continuum emission. For gas emissions, we detected $\rm{}^{12}CO$(2$-$1), $\rm{}^{13}CO$(2$-$1) and $\rm{}C^{18}O$(2$-$1), from which we estimated the total gas mass of the disk to be $2.2\times10^{-4}M_\odot$, assuming a CO/H$_2$ ratio of $10^{-4}$. The observations showed that the gas is present inside the dust cavity. The analysis of the velocity structure of the $\rm{}^{12}CO$(2$-$1) emission line revealed that the velocity is distorted at the location of the dust inner disk, which may be owing to warping of the disk or radial gas flow within the cavity of the dust disk. High-resolution observations of SY Cha showed that this system is composed of a ring and a distorted inner disk, which may be common, as indicated by the survey of transitional disk systems at a resolution of $\sim$0.1~arcsec.

The $\Lambda$CDM model provides a good fit to most astronomical observations but harbors large areas of phenomenology and ignorance. With the improvements in the precision and number of observations, discrepancies between key cosmological parameters of this model have emerged. Among them, the most notable tension is the 4$\sigma$ to 6$\sigma$ deviation between the Hubble constant ($H_{0}$) estimations measured by the local distance ladder and the cosmic microwave background (CMB) measurement. In this review, we revisit the $H_{0}$ tension based on the latest research and sort out evidence from solutions to this tension that might imply new physics beyond the $\Lambda$CDM model. The evidence leans more towards modifying the late-time universe.

In this chapter we examine how our knowledge of present day Venus can inform terrestrial exoplanetary science and how exoplanetary science can inform our study of Venus. In a superficial way the contrasts in knowledge appear stark. We have been looking at Venus for millennia and studying it via telescopic observations for centuries. Spacecraft observations began with Mariner 2 in 1962 when we confirmed that Venus was a hothouse planet, rather than the tropical paradise science fiction pictured. As long as our level of exploration and understanding of Venus remains far below that of Mars, major questions will endure. On the other hand, exoplanetary science has grown leaps and bounds since the discovery of Pegasus 51b in 1995, not too long after the golden years of Venus spacecraft missions came to an end with the Magellan Mission in 1994. Multi-million to billion dollar/euro exoplanet focused spacecraft missions such as JWST, and its successors will be flown in the coming decades. At the same time, excitement about Venus exploration is blooming again with a number of confirmed and proposed missions in the coming decades from India, Russia, Japan, the European Space Agency and the National Aeronautics and Space Administration. In this chapter, we review what is known and what we may discover tomorrow in complementary studies of Venus and its exoplanetary cousins.

We discuss a cosmological scenario with a stochastic background of gravitational waves sourced by the tensor perturbation due to a hybrid inflationary model with cubic potential. The tensor-to-scalar ratio for the present hybrid inflationary model is obtained as $r \approx 0.0006$. Gravitational wave spectrum of this stochastic background, for large-scale CMB modes, $10^{-4}Mpc^{-1}$ to $1Mpc^{-1}$ is studied. The present-day energy spectrum of gravitational waves $\Omega_0^{gw}(f)$ is sensitively related to the tensor power spectrum and r which is, in turn, dependent on the unknown physics of the early cosmos. This uncertainty is characterized by two parameters: $\hat{n_t(f)}$ logarithmic average over the primordial tensor spectral index and $\hat{w(f)}$ logarithmic average over the effective equation of state parameter. Thus, exact constraints in the $\hat{w(f)}$, $\hat{n_t(f)}$ plane can be obtained by comparing theoretical constraints of our model on r and $\Omega_0^{gw}(f)$. We obtain a limit on $\hat{w}(10^{-15}Hz)$<$0.33$ around the modes probed by CMB scales.

Aim. The satellite systems of Milky-Way-like galaxies offer a useful means by which to study the galaxy formation process in the cosmological context. It has been suggested that the currently observed anisotropic distribution of the satellites in such galaxy systems is inconsistent with the concordant $\Lambda CDM$ cosmology model on the galactic scale if the observed satellites are random samples of the dark matter (DM) sub-halos that are nearly isotropically distributed around the central galaxy. Here we examine such systems in the context of a simulation of the Local Group. Methods. We use high resolution hydrodynamical simulations of the formation and evolution of Milky-Way-like galaxy systems to investigate the spatial distribution of the luminous satellites and DM sub-halos around the central galaxy. Results. We find that the satellites and DM sub-halos are indeed anisotropically distributed with a direction roughly along the local DM filament. The observed alignment of the satellites in these systems is likely the result of preferred accretion along the local DM filaments in the galaxy formation process. For the dynamical properties of the satellites, we find that the direction of the angular momentum vector of the whole satellite system is different from the normal direction of the fitted "disk of satellites" (DOS) and the normal direction of the velocity dispersion of the system. Hence, the fitted DOS appears to be infalling and is not rotationally supported in this study.

MaGiV-1 is a candidate triple ellipsoidal star system in Vulpecula at coordinates RA(J2000) 19:52:19.13 DEC(J2000) +23:29:59.7 classified as ELL+ELL, number 2344411 in AAVSO VSX database. Through photometry from the TESS Space Telescope, two significant periods describing the orbital times of the components were identified using the Fourier transform. The analysis led to determining P(A-BC) = 4.269d the orbital period of A-BC pair, the primary component with the secondary component described by another pair, and P(BC) = 0.610d the orbital period of B-C pair, the inner ellipsoidal system. However, it cannot be completely ruled out that the shorter period can be explained by pulsations of one of the two components (e.g. by the GDOR type).

We present a computationally efficient expectation-maximization framework for multi-frame image deconvolution and super-resolution. Our method is well adapted for processing large scale imaging data from modern astronomical surveys. Our Tensorflow implementation is flexible, benefits from advanced algorithmic solutions, and allows users to seamlessly leverage Graphical Processing Unit (GPU) acceleration, thus making it viable for use in modern astronomical software pipelines. The testbed for our method is a set of $4$K by $4$K Hyper Suprime-Cam exposures, which are closest in terms of quality to imaging data from the upcoming Rubin Observatory. The preliminary results are extremely promising: our method produces a high-fidelity non-parametric reconstruction of the night sky, from which we recover unprecedented details such as the shape of the spiral arms of galaxies, while also managing to deconvolve stars perfectly into essentially single pixels.

The effects of a scalar field, known as the "assistant field," which nonminimally couples to gravity, on single-field inflationary models are studied. The analysis provides analytical expressions for inflationary observables such as the spectral index ($n_s$), the tensor-to-scalar ratio ($r$), and the local-type nonlinearity parameter ($f_{\rm NL}^{(\rm local)}$). The presence of the assistant field leads to a lowering of $n_s$ and $r$ in most of the parameter space, compared to the original predictions. In some cases, $n_s$ may increase due to the assistant field. This revives compatibility between ruled-out single-field models and the latest observations by Planck-BICEP/Keck. The results are demonstrated using three example models: loop inflation, power-law inflation, and hybrid inflation.

HS 2325+8205 is a long-period eclipsing dwarf nova with an orbital period above the period gap (Porb>3 h) and is reported to be a Z Cam-type dwarf nova. Based on the photometry of the Transiting Exoplanet Survey Satellite (TESS), the light variation and the quasi-periodic oscillation (QPOs) of HS 2325+8205 are studied. Using Continuous Wavelet Transform (CWT), Lomb-Scargle Periodogram (LSP), and sine fitting methods, we find for the first time that there is a QPOs of ~ 2160s in the long outburst top light curves of HS 2325+8205. Moreover, we find that the oscillation intensity of the QPOs of HS 2325+8205 is related to the orbital phase, and the intensity in orbital phases 0.5-0.9 are stronger than in orbital phases 0.1-0.5. Therefore, the relationship between the oscillation intensity of QPOs and the orbital phase may become a research window for the origin of QPOs. In addition, we use the LSP to correct the orbital period of HS 2325+8205 as 0.19433475(6) d.

Only three of the dozen central compact objects (CCOs) in supernova remnants (SNRs) show thermal X-ray pulsations due to non-uniform surface temperature (hot-spots). The absence of X-ray pulsations from several unpulsed CCOs has motivated suggestions that they have uniform-temperature carbon atmospheres (UTCAs), which adequately fit their spectra with appropriate neutron star (NS) surface areas. This is in contrast to the two-temperature blackbody or hydrogen atmospheres that also fit well. Here we investigate the applicability of UTCAs to CCOs. We show the following: (i) The phase-averaged spectra of the three pulsed CCOs can also be fitted with a UTCA of the appropriate NS area, despite pulsed CCOs manifestly having non-uniform surface temperature. A good spectral fit is therefore not strong support for the UTCA model of unpulsed CCOs. (ii) An improved spectrum of one unpulsed CCO, previously analyzed with a UTCA, does not allow an acceptable fit. (iii) For two unpulsed CCOs, the UTCA does not allow a distance compatible with the SNR distance. These results imply that, in general, CCOs must have hot, localized regions on the NS surface. We derive new X-ray pulse modulation upper limits on the unpulsed CCOs, and constrain their hot spot sizes and locations. We develop an alternative model that accounts for both the pulsed and unpulsed CCOs: a range of angles between hot spot and rotation axes consistent with an exponential distribution with scale factor $\lambda \sim 20^{\circ}$. We discuss physical mechanisms that could produce such small angles and small hot-spots.

We present the {\it XMM-Newton} X-ray analysis of 19 X-ray luminous galaxy clusters of low-to-mid redshift ($< 0.4$) selected from the MCXC cluster catalogue in the Hyper Supri%survey as the first work in our series paper. We derive the hydrostatic equilibrium mass and study scaling relations using i) the whole sample, ii) only relaxed clusters and iii) only disturbed clusters. When considering the whole sample, the $Y_{\rm X}$-$M_{\rm tot}$ and $M_{\rm gas}$-$M_{\rm tot}$ relations agree with self-similarity. In terms of morphology, relaxed clusters show a flatter relation in $L_{\rm X,ce}$-$M_{\rm tot}$, $L_{\rm X,bol}$-$M_{\rm tot}$, $L_{\rm X,ce}$-$T$, $L_{\rm bol,ce}$-$T$, $M_{\rm gas}$-$M_{\rm tot}$ and $Y_{\rm X}$-$M_{\rm tot}$. The $L_{\rm bol,ce}$-$M_{\rm tot}$, $L_{\rm X,ce}$-$M_{\rm tot}$ $L_{\rm bol,ce}$-$T$ and $L_{\rm X,ce}$-$T$ relations show a slope $\sim$3$\sigma$ steeper. The residuals in the $M_{\rm gas}$-$M_{\rm tot}$ and $T$-$M_{\rm tot}$ relations and the intrinsic covariance between $M_{\rm gas}$ and $T$ show hints of positive correlation, casting doubt on whether the $Y_{\rm X}$ parameter is a truly low scatter mass proxy. The $M_{\rm gas}$-$M_{\rm tot}$ and $T$-$M_{\rm tot}$ plots color-coded with the offset of the $L_{\rm X,ce}$-$M_{\rm tot}$ relation show these two relations to be brightness dependent but not the $L_{\rm X,ce}$-$T$ relation, suggesting relations involving $M_{\rm tot}$ are biased due to sample selection based on luminosity. Following \citet{2017A&A...606A..24A} and combining our result with literature studies, we find the $M_{\rm tot}$ derived not using mass proxies deviate from $L_{\rm X}$ $\propto$ $M_{\rm gas}^{2}M_{\rm tot}^{-1}$ and $M_{\rm tot}$ based on hydrostatic equilibrium are more massive than what is expected by their relation using caustic masses. This indicates mass bias plays an important role in scaling relations.

With the well-established correlation between the relative stabilities of isomers and their interstellar abundances coupled with the prevalence of isomeric species among the interstellar molecular species, isomerization remains a plausible formation route for isomers in the interstellar medium. The present work reports an extensive investigation of the isomerization energies of 246 molecular species from 65 isomeric groups using the Gaussian-4 theory composite method with atoms ranging from 3 to 12. From the results, the high abundances of the most stable isomers coupled with the energy sources in interstellar medium drive the isomerization process even for barriers as high as 67.4 kcal/mol. Specifically, the cyanides and their corresponding isocyanides pairs appear to be effectively synthesized via this process. The following potential interstellar molecules; CNC, NCCN, c-C5H, methylene ketene, methyl Ketene, CH3SCH3, C5O, 1,1-ethanediol, propanoic acid, propan-2-ol, and propanol are identified and discussed. The study further reaffirms the importance of thermodynamics in interstellar formation processes on a larger scale and accounts for the known isomeric species. In all the isomeric groups, isomerization appears to be an effective route for the formation of the less stable isomers (which are probably less abundant) from the most stable ones that are perhaps more abundant.

This article presents the early results of synchronous multiwavelength observations of one of the brightest gamma-ray bursts (GRBs) GRB 160625B with the detailed continuous fast optical photometry of its optical counterpart obtained by MASTER and with hard X-ray and gamma-ray emission, obtained by the Lomonosov and Konus-Wind spacecraft. The detailed photometry led us to detect the quasi-periodical emission components in the intrinsic optical emission. As a result of our analysis of synchronous multiwavelength observations, we propose a three-stage collapse scenario for this long and bright GRB. We suggest that quasiperiodic fluctuations may be associated with forced precession of a self-gravitating rapidly rotating superdense body (spinar), whose evolution is determined by a powerful magnetic field. The spinar's mass allows it to collapse into a black hole at the end of evolution.

We consider planets composed of water ice and rock, located far from a central star. In an earlier study, computing the growth of planets by continuous accretion, we found that a large fraction of the ice evaporates upon accretion, creating a water vapor atmosphere. Here we consider accretion as a discrete series of planetesimal impacts (of order $10^8$), at the same time-dependent accretion rate, and investigate the fate of the vapor, as a result of its interaction with the accreting planetesimals. We find that a large fraction of the vapor escapes. The remaining fraction may form an outer layer of ice after the termination of accretion and cooling of the surface. The escaped water mass may significantly alter the ice-to-rock ratio of the planet. We investigate the effect of different choices of parameters such as the ice-to-rock ratio, the planetesimal size distribution, and the impact velocities. We find that the planetesimal size distribution has a negligible effect and explain why. By contrast, the ice-to-rock ratio and impact velocities affect the fraction of retained water masses considerably.

We report on MASTER optical observations of an afterglow-like optical and X-ray transient AT2021lfa/ZTF21aayokph. We detected the initial steady brightening of the transient at 7{\sigma} confidence level. This allowed us to use smooth optical self-similar emission of GRBs model to constrain the explosion time to better than 14 min as well as to estimate its initial Lorentz factor {\Gamma}0 = 20 +/- 10. Taking into consideration the low {\Gamma}0 and non-detection in gamma-rays, we classify this transient as the first failed GRB afterglow.

We present our study of stability of differentially rotating, axisymmetric neutron stars described by a polytropic equation of state with $\Gamma = 2$. We focus on quasi-toroidal solutions with a degree of differential rotation $\widetilde A=1$. Our results show that for a wide range of parameters hypermassive, quasi-toroidal neutron stars are dynamically stable against quasi-radial perturbations, which may have implications for newly born neutron stars and binary neutron stars mergers.

We presented a study of outburst activity in the BL Lacertae source OJ~287, observed extensively with the X-ray telescope (XRT) onboard Neil Gehrels Swift Observatory. We demonstrated that the results of our analysis of X-ray flaring activity using the Swift/XRT data allow to refine the key characteristics of the OJ~287 secondary (its nature and mass). We discovered that the energy spectra in all spectral states can be fitted using the XSPEC Bulk Motion Comptonization (BMC) model. As a result we found that the X-ray photon index of the BMC model, Gamma correlates with the mass accretion rate, dot. We established that Gamma increases monotonically with dot from the low-hard state, Gamma~ 1.5 to the high-soft state, Gamma~2.8 and finally saturates. The index behavior was similar to that in a number of black hole (BH) candidates in which we showed that its saturation was an observational evidence of the presence of a BH. Based on this correlation, we applied a scaling method and determined that a secondary BH mass in OJ~287 is about ~1.25x10^8 solar masses, using the well-studied X-ray BH binaries XTE~1550--564, H~1743--322, 4U~1630--47, GRS~1915+105 as well as extragalactic BHs ESO~243--49 and M101~ULX--1, as reference sources. Also using the power spectrum analysis we inferred the size of the Compton cloud L_{CC}~ 10^{13} cm where X-ray spectra were formed. Using this value of L_{CC} we confirmed that a BH mass of the secondary in OJ~287 was of order of 10^8 solar masses as we derived using the index, Gamma-correlation (the scaling method) with respect of the mass accretion rate.

We use the Subaru Hyper Suprime-Cam (HSC) data to study structural parameters and systemic proper motion of the Sextans dwarf spheroidal galaxy at the heliocentric distance of 86 kpc, which is one of the most important targets for studies of dark matter nature and galaxy formation physics. Thanks to the superb image quality and wide area coverage of the Sextans field, the HSC data enables a secure selection of member star candidates based on the colour-magnitude cut, yielding about 10,000 member candidates at magnitudes down to $i\sim 24$. We use a likelihood analysis of the two-dimensional distribution of stars to estimate the structural parameters of Sextans taking into account the contamination of foreground halo stars in the Milky Way, and find that the member star distribution is well-fitted by an elliptical King profile with ellipticity $\epsilon \simeq 0.25$ and the core and tidal radii of $R_c=(368.2\pm 8.7)$ pc and $R_t=(2.54\pm 0.045)$ kpc, respectively. Then using the two HSC datasets of 2.66 years time baseline on average, we find the systemic proper motions of Sextans to be $(\mu_\alpha, \mu_\delta)=(-0.508^{+0.071}_{-0.071},0.083^{+0.071}_{-0.072})$ $\mathrm{mas}\ \mathrm{yr}^{-1}$, which is consistent with some of the previous works using the $Gaia$ data of relatively bright member stars in Sextans. Thus, our results give a demonstration that a ground-based, large-aperture telescope data which covers a wide solid angle of the sky and has a long time baseline, such as the upcoming LSST data, can be used to study systemic proper motions of dwarf galaxies.

Actively accreting young stellar objects show H$\alpha$ emission line in their spectra. We present the results of survey for H$\alpha$ emission-line stars in the star-forming region IC~5070 taken with 2-m Himalyan {\it Chandra} Telescope. Based on the H$\alpha$ slitless spectroscopy data, we identified 131 emission-line stars in $\sim$ 0.29 square degrees area of the IC~5070 region. Using Gaia Early Data Release 3, we estimated the mean proper motion and parallax of the emission-line stars. We also estimated the mean distance and reddening toward the region using the emission-line stars, that are $\sim$ 833 pc and $\sim$2 mag, respectively. By examining the locations of these stars in the color-magnitude diagrams constructed using Gaia and PanSTARRS1 data, we found that a majority of the H$\alpha$ emitters are young low-mass ($<$ 1.5 $M\odot$) stars. We also compared our catalog of emission-line stars with the available young stellar catalogs and found that most of them are Class~{II}/ flat spectrum sources with the spectral type ranging from K to M. Based on the proper motion/ parallax values and locations on the color-magnitude diagrams, about 20 emission-line stars are flagged as non-members. The relative proper motion of the emission-line stars with respect to the ionizing source suggest the possibility of the `rocket effect' scenario in the remnant cloud (BRC~31).

The jet composition of gamma-ray bursts (GRBs), as well as how efficiently the jet converts its energy to radiation, are long-standing problems in GRB physics. Here, we reported a comprehensive temporal and spectral analysis of the TeV-emitting bright GRB 190114C. Its high fluence ($\sim$ 4.436$\times$10$^{-4}$ erg cm$^{-2}$) allows us to conduct the time-resolved spectral analysis in great detail and study their variations down to a very short time-scale ($\sim$0.1 s) while preserving a high significance. Its prompt emission consists of three well-separated pulses. The first two main pulses ($P_1$ and $P_2$) exhibit independently strong thermal components, and starting from the third pulse ($P_3$) and extending to the entire afterglow, the spectra are all non-thermal, the synchrotron plus Compton up-scattering model well interpret the observation. By combining the thermal ($P_1$ and $P_2$) and the non-thermal ($P_3$) observations based on two different scenarios (global and pulse properties) and following the method described in \cite{Zhang2021}, we measure the fireball parameters and GRB radiative efficiency with little uncertainties for this GRB. A relevantly high GRB radiative efficiency is obtained based on both the global and pulse properties, suggesting that if GRBs are powered by fireballs, the efficiency can be high sometimes. More interestingly, though the observed parameters are individually different (e.g., the amount of mass loading $M$), the radiative efficiency obtained from $P_1$ ($\eta_\gamma=36.0\pm6.5\%$) and $P_2$ ($\eta_\gamma=41.1\pm1.9\%$) is roughly the same, which implies that the central engine of the same GRB has some common properties.

We investigate dispersive and kinetic effects on the evolution of a two-dimensional kinked Alfv\'en wave packet by comparing results from MHD, Hall-MHD and hybrid simulations of a low-$\beta$ plasma. We find that the Hall term determines the overall evolution of the wave packet over a characteristic time $\tau^*=\tau_a\ell/d_i$ in both fluid and hybrid models. Dispersion of the wave packet leads to the conversion of the wave energy into internal plasma energy. When kinetic protons are considered, the proton internal energy increase has contributions from both plasma compressions and phase space mixing. The latter occurs in the direction parallel to the guiding mean magnetic field, due to protons resonating at the Alfv\'en speed with a compressible mode forced by the wave packet. Implications of our results for switchbacks observations and solar wind energetics are discussed.

We propose a simple analytic fit for the power spectrum of scalar (curvature) perturbations during inflation, in order to describe slow roll of inflaton and formation of primordial black holes in the early universe, in the framework of single-field models. Our fit is given by a sum of the power spectrum in the slow-roll approximation, needed for a viable description of the cosmic microwave background radiation in agreement with Planck/BICEP/Keck measurements, and the log-normal (Gaussian) fit for the power spectrum enhancement (peak) needed for efficient production of primordial black holes. We use the T-type $\alpha$-attractor models in order to describe slow-roll inflation. Demanding the location and height of the peak to yield the masses of primordial black holes in the asteroid-size window allowed for the whole (current) dark matter to be composed of the primordial black holes, we find the restrictions on the remaining parameters and, most notably, on the width of the peak.

Aims. Using the direct imaging technique we searched for low mass companions around the star AF Lep that presents a significant proper motion anomaly (PMa) signal obtained from the comparison of Hipparcos and Gaia eDR3 catalogs. Methods. We observed AF Lep in two epochs with VLT/SPHERE using its subsystems IFS and IRDIS in the near-infrared (NIR) covering wavelengths ranging from the Y to the K spectral bands (between 0.95 and 2.3 {\mu}m). The data were then reduced using the high-contrast imaging techniques angular differential imaging (ADI) and spectral differential imaging (SDI) to be able to retrieve the signal from low mass companions of the star. Results. A faint companion was retrieved at a separation of ~0.335" from the star and with a position angle of ~70.5 deg in the first epoch and with a similar position in the second epoch. This corresponds to a projected separation of ~9 au. The extracted photometry allowed us to estimate for the companion a mass between 2 and 5 MJup. This mass is in good agreement with what is expected for the dynamic mass of the companion deduced using astrometric measures (5.2-5.5 MJup). This is the first companion with a mass well below the deuterium burning limit discovered coupling direct imaging with PMa measures. Orbit fitting done using the orvara tool allowed to further confirm the companion mass and to define its main orbital parameters.

Dwarf galaxies are characterised by a very low luminosity and low mass. Because of significant accretion and ejection activity of massive black holes, some dwarf galaxies also host low-luminosity active galactic nuclei (AGNs). In a few dwarf AGNs, very long baseline interferometry (VLBI) observations have found faint non-thermal radio emission. SDSS J090613.77+561015.2 is a dwarf AGN owning an intermediate-mass black hole (IMBH) with a mass of $M_{BH} = 3.6^{+5.9}_{-2.3} \times 10^5 M_{sun}$ and showing a rarely-seen two-component radio structure in its radio nucleus. To further probe their nature, i.e. the IMBH jet activity, we performed additional deep observations with the European VLBI Network (EVN) at 1.66 GHz and 4.99 GHz. We find the more diffuse emission regions and structure details. These new EVN imaging results allow us to reveal a two-sided jet morphology with a size up to about 150 mas (projected length $\sim$140 pc) and a radio luminosity of about $3\times10^{38}$ erg s$^{-1}$. The peak feature has an optically thin radio spectrum and thus more likely represents a relatively young ejecta instead of a jet base. The EVN study on SDSS J090613.77+561015.2 demonstrates the existence of episodic, relatively large-scale and powerful IMBH jet activity in dwarf AGNs. Moreover, we collected a small sample of VLBI-detected dwarf AGNs and investigated their connections with normal AGNs. We notice that these radio sources in the dwarf AGNs tend to have steep spectra and small linear sizes, and possibly represent ejecta from scaled-down episodic jet activity.

Studies of the time-domain structure of fast radio bursts (FRBs) require an accurate estimate of the FRB dispersion measure in order to recover the intrinsic burst shape. Furthermore, the exact DM is itself of interest when studying the time-evolution of the medium through which multiple bursts from repeating FRBs propagate. A commonly used approach to obtain the dispersion measure is to take the value that maximizes the FRB structure in the time domain. However, various authors use differing methods to obtain this structure parameter, and do not document the smoothing method used. Furthermore, there are no quantitative estimates of the error in this procedure in the FRB literature. In this letter, we present a smoothing filter based on the discrete cosine transform, and show that computing the structure parameter by summing the squares of the derivatives and taking the square root immediately lends itself to calculation of uncertainty of the structure parameter. We illustrate this with FRB181112 and FRB210117 data, which were detected by the Australian Square Kilometre Array Pathfinder, and for which high-time-resolution data is available.

NaCl is a diatomic molecule with a large dipole moment, which allows for its detection even at relatively small abundances. It has been detected towards several evolved stars, among which is the AGB star IK Tau, around which it is distributed in several clumps that lie off-center from the star. We aim to study the three-dimensional distribution of NaCl around the AGB star IK Tau, and to obtain the abundance of NaCl relative to H$_2$ for each of the clumps. First, a new value for the maximum expansion velocity is determined. The observed ALMA channel maps are then deprojected to create a three-dimensional model of the distribution of NaCl. This model is then used as input for the radiative transfer modelling code magritte, which is used to obtain the NaCl abundances of each of the clumps by comparing the observations with the results of the magritte simulations. Additionally, the rotational temperature of the clumps is determined using population diagrams. We derive an updated value for the maximum expansion velocity of IK Tau $\upsilon_\mathrm{exp}$ = 28.4 km/s. A spiral-like shape can be discerned in our three-dimensional distribution model of the NaCl. This spiral lies more or less in the plane of the sky. The distribution is also flatter in the line-of-sight direction than in the plane of the sky. We find clump abundances between $9 \times 10^{-9}$ and $5 \times 10^{-8}$ relative to H$_2$, where the relative abundance is typically lower for clumps closer to the star. For the first time, we used deprojection to understand the three-dimensional environment of an AGB star and calculated the fractional abundance of NaCl in clumps surrounding the star.

We provide novel constraints on the parameters defining the universal pressure profile (UPP) within clusters of galaxies, and explore their dependence on the cluster mass and redshift, from measurements of Sunyaev-Zel'dovich Compton-$y$ profiles. We employ both the $\textit{Planck}$ 2015 MILCA and the ACT-DR4 $y$ maps over the common $\sim 2,100\,\text{deg}^2$ footprint. We combine existing cluster catalogs based on KiDS, SDSS and DESI observations, for a total of 23,820 clusters spanning the mass range $10^{14.0}\,\text{M}_{\odot}<M_{500}<10^{15.1}\,\text{M}_{\odot}$ and the redshift range $0.02<z<0.98$. We split the clusters into three independent bins in mass and redshift; for each combination we detect the stacked SZ cluster signal and extract the mean $y$ angular profile. The latter is predicted theoretically adopting a halo model framework, and MCMCs are employed to estimate the UPP parameters, the hydrostatic mass bias $b_{\rm h}$ and possible cluster miscentering effects. We constrain $[P_0,c_{500},\alpha,\beta]$ to $[5.9,2.0,1.8,4.9]$ with $\textit{Planck}$ and to $[3.8,1.3,1.0,4.4]$ with ACT using the full cluster sample, in agreement with previous findings. We do not find any compelling evidence for a residual mass or redshift dependence, thus expanding the validity of the cluster pressure profile over much larger $M_{500}$ and $z$ ranges; this is the first time the model has been tested on such a large (complete and representative) cluster sample. Finally, we obtain loose constraints on the hydrostatic mass bias in the range 0.2-0.3, again in broad agreement with previous works.

We aim to observationally and statistically constrain the influence of flybys in the formation and evolution of debris disks. We compiled a sample of 254 debris disks with ages between 2 Myr and 8 Gyr that are either part of an association or isolated, drawing the binary and planetary companions of the systems mainly from the literature. Using the Gaia eDR3 astrometric data and radial velocities of our sample, as well as all the sources in a specific region of the sky, we reconstructed the relative linear motions in the last 5 Myr and made predictions for the next 2 Myr. Relating the Hill radius of each debris disk system and the closest distances reached by the two sources, we defined the flyby events in terms of position and time. We find that in the period between the last 5 Myrs and the next 2 Myrs, 90% of the analyzed systems have experienced at least a close flyby, while 7% of them have experienced flybys at distances greater than 0.5R Hill. In particular, 75% of them have experienced at least one past close encounter and 36% multiple past close encounters. From the sub-sample of resolved debris disk (41 out of 94), 80% of the analyzed systems experience at least an encounter within 0.8 pc. From the subsample of 10 debris disks with planets, half of these systems do show misalignments between disk and planet, stirring, or asymmetries. Systems with a misalignment between the planetary orbit and the disk do indeed experience at least one flyby event. In particular, when the planet orbits have a difference with the disk inclination higher than about 20 degree, as in the case of HD 38529, we find that multiple close encounters have taken place in the last 5 Myr, as theoretically predicted. The high incidence of encounters, particularly close encounters, experienced by the systems in the last 5 Myr suggests the fundamental impact of flybys on the evolution of debris disks.

The three-dimensional extinction and structure are studied for the Taurus, Orion, Perseus and California molecular clouds based on the LAMOST spectroscopy. Stellar color excess is calculated with the intrinsic color index derived from the atmospheric parameters in the LAMOST DR8 catalog and the observed color index in the Gaia EDR3 and the 2MASS PSC. In combination with the distance from the Gaia EDR3 parallax, the three-dimensional dust extinction maps are retrieved in the color excesses $E_{\rm{G_{BP},G_{RP}}}$ and $E_{\rm{J,K_{S}}}$ with an uncertainty of $\sim$0.03mag and $\sim$0.07mag respectively. The extinction maps successfully separate the clouds that overlap in the sky area and manifest the structure of the individual cloud. Meanwhile, a bow-like structure is found with a distance range from 175pc to 250pc, half of which is a part of the Per-Tau Shell in similar coordinates and distance while the other half is not. Three low-extinction rings are additionally discovered and briefly discussed.

(abridged) Expanding the sample of directly imaged companions to nearby, young stars that are amenable to detailed astrometric and spectroscopic studies is critical for the continued development and validation of theories of their evolution and atmospheric processes. The recent release of the {\it Gaia} astrometric catalogue allows us to efficiently search for these elusive companions by targeting those stars that exhibit the astrometric reflex motion induced by an orbiting companion. The nearby (27 pc), young (24 Myr) star AF Leporis (AF Lep) was targeted because of its astrometric acceleration, consistent with a wide-orbit planetary companion detectable with high-contrast imaging. We used the SPHERE instrument on the VLT to search for faint substellar companions in the immediate vicinity of AF Lep. We used observations of a nearby star interleaved with those of AF Lep to efficiently subtract the residual point spread function. This provided sensitivity to faint planetary-mass companions within 1 arcsec ($\sim$30 au) of the star. We detected the companion AF Lep b at a separation of 339 mas (9 au), within the inner edge of its unresolved debris disk. The measured $K$-band contrast and the age of the star yield a model-dependent mass of 4--6 $M_{\rm Jup}$, consistent with the mass derived from an orbital fit of $4.3_{-1.2}^{+2.9}$ $M_{\rm Jup}$. The near-infrared SED of the planet is consistent with an object at the L--T transition, but under-luminous with respect to field-gravity objects. AF Lep b joins a growing number of substellar companions imaged around stars in the young $\beta$ Pic moving group. With a mass of between 3--7 $M_{\rm Jup}$, it occupies a gap in this isochronal sequence between the hotter, more massive companions like PZ~Tel~B and $\beta$~Pic~b, and the cooler 51~Eri~b which is sufficiently cool for the formation of methane within its photosphere.

Photonic Integrated Circuits (PIC) are best known for their important role in the telecommunication sector, e.g. high speed communication devices in data centers. However, PIC also hold the promise for innovation in sectors like life science, medicine, sensing, automotive etc. The past two decades have seen efforts of utilizing PIC to enhance the performance of instrumentation for astronomical telescopes, perhaps the most spectacular example being the integrated optics beam combiner for the interferometer GRAVITY at the ESO Very Large Telescope. This instrument has enabled observations of the supermassive black hole in the center of the Milky Way at unprecedented angular resolution, eventually leading to the Nobel Price for Physics in 2020. Several groups worldwide are actively engaged in the emerging field of astrophotonics research, amongst them the innoFSPEC Center in Potsdam, Germany. We present results for a number of applications developed at innoFSPEC, notably PIC for integrated photonic spectrographs on the basis of arrayed waveguide gratings and the PAWS demonstrator (Potsdam Arrayed Waveguide Spectrograph), PIC-based ring resonators in astronomical frequency combs for precision wavelength calibration, discrete beam combiners (DBC) for large astronomical interferometers, as well as aperiodic fiber Bragg gratings for complex astronomical filters and their possible derivatives in PIC.

I review the scientific and technical history of the Search for Extraterrestrial Intelligence (SETI), discuss the impact of the political involvement, and speculate on the nature of a successful detection and its potential social and cultural impact. Emphasis is on the development of SETI in the United States and the complementary progress in the Former Soviet Union.

This article is based on the tutorial we gave at the hands-on workshop of the ICRANet-ISFAHAN Astronomy Meeting. We first introduce the basic theory of machine learning and sort out the whole process of training a neural network. We then demonstrate this process with an example of inferring redshifts from SDSS spectra. To emphasize that machine learning for astronomy is easy to get started, we demonstrate that the most basic CNN network can be used to obtain high accuracy, we also show that with simple modifications, the network can be converted for classification problems and also to processing gravitational wave data.

Recent anisotropy studies of UHECR data at energies $\gtrsim$ 40 EeV, have disclosed a correlation of their angular distribution with the extragalactic local structure, specifically with either local starburst galaxies or AGN. Using Monte Carlo simulations taking into account photo-disintegration processes, we further explore a framework in which these UHECRs were accelerated by Centaurus A in a recent powerful outburst before being scattered by magnetic fields associated with local, Council of Giant, extragalactic structure. We find that the observed intermediate scale anisotropies can be accounted for by the Council of Giant structure imposing a response function on the initial outburst of UHECRs from a single source located at Centaurus A's position. The presence of these local structures create `echoes' of UHECRs after the initial impulse, and focusing effects. The strongest echo wave has a lag of $\sim$ 20 Myr, comparable to the age of synchrotron-emitting electrons in the giant Centaurus A lobes. Through consideration of the composition of both the direct and echo wave components, we find that the distribution of the light (1$<\ln A<$1.5) component across the sky offers exciting prospects for testing the echo model using future facilities such as Auger prime. Our results demonstrate the potential that UHECR nuclei offer, as "composition clocks", for probing propagation scenarios from local sources.

We conducted a deep XMM$\unicode{x2013}$Newton observing campaign on the 23.5-s radio pulsar PSR J0250+5854 in order to better understand the connection between long-period, radio-emitting neutron stars and their high-energy-emitting counterparts. No X-ray emission was detected resulting in an upper limit in the bolometric luminosity of PSR J0250+5854 of $<$10$^{31}$ erg s$^{-1}$ for an assumed blackbody with a temperature of 85 eV, typical of an X-ray Dim Isolated Neutron Star (XDINS). We compared the upper limit in the bolometric luminosity of PSR J0250+5854 with the known population of XDINSs and found that the upper limit is lower than the bolometric luminosity of all but one XDINS. We also compared PSR J0250+5854 with SGR 0418+5729, the magnetar with low dipole magnetic field strength, where the upper limit suggests that if PSR J0250+5854 has a thermal hot spot like SGR 0418+5729, it would have a blackbody temperature of $<$200 eV, compared to 320 eV of the magnetar.

We generate a model description of the solar wind based on an explicit wave-turbulence-driven heating mechanism, and constrain our model with observational data. We included an explicit coronal heating source term in the general 3D magnetohydrodynamic code NIRVANA to simulate the properties of the solar wind. The adapted heating mechanism is based on the interaction and subsequent dissipation of counter-propagating Alfv\'en waves in the solar corona, accounting for a turbulent heating rate Q_p. The solar magnetic field is assumed to be an axisymmetric dipole with a field strength of 1 G. Our model results are validated against observational data taken by the Parker Solar Probe (PSP). Our NIRwave solar wind model reconstructs the bimodal structure of the solar wind with slow and fast wind speeds of 410 km/s and 650 km/s respectively. The global mass-loss rate of our solar wind model is 2.6e-14 solar masses per year. Despite implementing simplified conditions to represent the solar magnetic field, the solar wind parameters characterising our steady-state solution are in reasonable agreement with previously established results and empirical constraints. The number density from our wind solution is in good agreement with the derived empirical constraints, with larger deviations for the radial velocity and temperature. In a comparison to a polytropic wind model generated with NIRVANA, we find that our NIRwave model is in better agreement with the observational constraints that we derive.

The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution.

Inner regions of protoplanetary disks host many complex physical processes such as star-disk interactions, magnetic fields, planet formation, and the migration of new planets. To directly study this region requires milli-arcsecond angular resolution, beyond the diffraction limit of the world's largest optical telescopes and even too small for the mm-wave interferometer ALMA. However, we can use infrared interferometers to image the inner astronomical unit. Here, we present new results from the CHARA and VLTI arrays for the young and luminous Herbig Be star HD 190073. We detect a sub-AU cavity surrounded by a ring-like structure that we interpret as the dust destruction front. We model the shape with 6 radial profiles, 3 symmetric and 3 asymmetric, and present a model-free image reconstruction. All the models are consistent with a near face-on disk with inclination $\lesssim 20^\circ$, and we measure an average ring radius of 1.4 $\pm 0.2$ mas (1.14 AU). Around $48\%$ of the total flux comes from the disk with ~$15\%$ of that emission appearing to emerge from inside the inner rim. The cause of emission is still unclear, perhaps due to different dust grain compositions or gas emission. The skewed models and the imaging point to an off-center star, possibly due to binarity. Our image shows a sub-AU structure, which seems to move between the two epochs inconsistently with Keplerian motion and we discuss possible explanations for this apparent change.

We have investigated the impact of circumplanetary rings consisting of spherical micrometer-sized particles on the net scattered light polarization of extrasolar gas giants. Using the three-dimensional Monte Carlo radiative transfer code POLARIS, we studied the impact of the macroscopic parameters that define the ring, such as its radius and inclination, and the chemical composition of the ring particles on the net scattered polarization. For the spherical ring particles, we applied the Mie scattering theory. We studied the flux and polarization of the scattered stellar radiation as a function of planetary phase angle and wavelength from the optical to the near-infrared. For the chosen grain size distribution, the dust particles in the ring show strong forward scattering at the considered wavelengths. Thus, the reflected flux of the planet dominates the total reflected and polarized flux at small phase angles. However, the scattered and polarized flux of the ring increase at large phase angles and exceeds the total reflected planetary flux. For large rings that contain silicate particles, the total reflected flux is dominated by the radiation scattered by the dust in the ring at all phase angles. As a result, the orientation of polarization is parallel to the scattering plane at small phase angles. In contrast, for a ring that contains water ice particles, the orientation of polarization is parallel to the scattering plane at large phase angles. Depending on the ring inclination and orientation, the total reflected and polarized flux show a specific distribution as well. Large particles show a strong polarization at large phase angles compared to smaller particles. For a Jupiter-like atmosphere that contains methane and aerosols, methane absorption features are missing in the spectrum of a ringed planet.

Single field inflationary models are investigated within Palatini quadratic gravity represented by $R+\alpha R^2$ along with a non-minimal coupling of the form $f(\phi) R$ between the inflaton field $\phi$ and the gravity. The treatment is performed in the Einstein frame, where the minimal coupling to gravity is recovered through conformal transformation. We consider various limits of the model with different inflationary scenarios characterized as canonical slow-roll inflation in the limit $\alpha \dot{\phi}^2\ll (1+f(\phi)) $, constant-roll k-inflation for $\alpha \ll 1$, and slow-roll K-inflation for$ \alpha \gg 1$ . A cosine and exponential potential are examined with the limits mentioned above and different well-motivated non-minimal couplings to gravity. We compare the theoretical results, exemplified by the tensor-to-scalar $r$ ratio and spectral index $n_s$, with the recent observational results of Planck 2018 $\&$ BICEP/Keck . Furthermore, we include the results of a new study forecast precision with which $n_s$ and $r$ can be constrained by currently envisaged observations, including CMB (Simons Observatory, CMB-S4, and LiteBIRD)

The distribution of the spin directions of spiral galaxies in the Sloan Digital Sky Survey has been a topic of debate in the past two decades, with conflicting conclusions reported even in cases where the same data was used. Here we follow one of the previous experiments by applying the SpArcFiRe algorithm to annotate the spin directions in original dataset of Galaxy Zoo 1. The annotation of the galaxy spin directions is done after a first step of selecting the spiral galaxies in three different manners: manual analysis by Galaxy Zoo classifications, by a model-driven computer analysis, and with no selection of spiral galaxies. The results show that when spiral galaxies are selected by Galaxy Zoo volunteers, the distribution of their spin directions as determined by SpArcFiRe is not random, which agrees with previous reports. When selecting the spiral galaxies using a model-driven computer analysis or without selecting the spiral galaxies at all, the distribution is also not random. Simple binomial distribution analysis shows that the probability of the parity violation to occur by chance is lower than 0.01. Fitting the spin directions as observed from Earth to cosine dependence exhibits a dipole axis with statistical strength of 2.33$\sigma$ to 3.97$\sigma$. These experiments show that regardless of the selection mechanism and the analysis method, all experiments show similar conclusions. These results are aligned with previous reports using other methods and telescopes, suggesting that the spin directions of spiral galaxies as observed from Earth exhibit a dipole axis formed by their spin directions. Possible explanations can be related to the large-scale structure of the Universe, or to internal structure of galaxies. The catalogs of annotated galaxies generated as part of this study are available.

In the last decade, growing evidence has emerged supporting a non-universal stellar Initial Mass Function (IMF) in massive galaxies, with a larger number of dwarf stars with respect to the Milky-Way (bottom-heavy IMF). However, a consensus about the mechanisms that cause IMF variations is yet to be reached. Recently, it has been suggested that stars formed early-on in cosmic time, via a star formation burst, could be characterised by a bottom-heavy IMF. A promising way to confirm this is to use relics, ultra-compact massive galaxies, almost entirely composed by these "pristine" stars. The INSPIRE Project aims at assembling a large sample of confirmed relics, that can serve as laboratory to investigate on the conditions of star formation in the first 1-3 Gyr of the Universe. In this third INSPIRE paper, we build a high signal-to-noise spectrum from five relics and one from five galaxies with similar sizes, masses, and kinematical properties, but characterised by a more extended star formation history (non-relics). Our detailed stellar population analysis suggests a systematically bottom-heavier IMF slope for relics than for non-relics, adding new observational evidence for the non-universality of the IMF at various redshifts and further supporting the above proposed physical scenario.

The Burst Observer and Optical Transient Exploring System (BOOTES) was first designed as an asset of autonomous telescopes that started to be deployed in 1998, taking 24 years to be fully developed around the Earth. Nowadays BOOTES has became a global network of robotic telescopes, being the first one present in all continents, as of 2022. Here we present the details of the network and review its achievements over the last two decades regarding follow-up observations of high-energy transient events. Moreover, considering the recent operations of neutrino and gravitational wave detectors, some hot-topic expectations related to robotic astronomy are discussed within the framework of multi-wavelength astrophysics.

We propose a novel set of Poisson Cluster Process models to detect Ultra-Diffuse Galaxies (UDGs), a recently-discovered class of galaxies that are challenging to detect and are of substantial interests in modern astrophysics. We construct an improved spatial birth-death-move MCMC algorithm to make inferences about the locations of these otherwise un-observable galaxies. Our novel models significantly out-perform existing approaches based on the Log-Gaussian Cox Process; the novel marked point process we propose can also improve the detection performance for UDGs in noisy environments. We find evidence of a potential new "dark galaxy" that was not detected using previous methods.

We present an revised table of 390 Galactic radio supernova remnants (SNRs) and their basic parameters. Statistical analyses are performed on SNR diameters, ages, spectral indices, Galactic heights and spherical symmetries. Furthermore, the accuracy of distances estimated using the $\Sigma$-D relation is examined. The arithmetic mean of the Galactic SNR diameters is $30.5$ pc with standard error $1.7$ pc and standard deviation $25.4$ pc. The geometric mean and geometric standard deviation factor of Galactic SNR diameters is $21.9$ pc and $2.4$, respectively. We estimate ages of 97 SNRs and find that the supernova (SN) birth rate to be lower than, but within $2\sigma$ of currently accepted values for SN birth rate. The mean spectral index of shell-type SNRs is $-0.51 \pm 0.01$ and no correlations are found between spectral indices and the SNR parameters of molecular cloud (MC) association, SN type, diameter, Galactic height and surface brightness. The Galactic height distribution of SNRs is best described by an exponential distribution with a scale height of $48 \pm 4$ pc. The spherical symmetry measured by the ovality of radio SNRs is not correlated to any other SNR parameters considered here or to explosion type.

Spiral density waves can arise in galactic disks as linear instabilities of the underlying stellar distribution function. Such instabilities grow exponentially in amplitude at some fixed growth rate $\beta$ before saturating nonlinearly. However, the mechanisms behind nonlinear saturation, and the resulting saturated spiral amplitude, have not received much attention. Here we argue that the most important nonlinear saturation mechanism is likely trapping of stars near the spiral's corotation resonance. Under these circumstances, we show analytically that an $m$-armed spiral will saturate when the libration frequency of resonantly trapped orbits reaches $\omega_\mathrm{lib} \sim 3\, m^{1/2} \beta$. For a galaxy with a flat rotation curve this implies a maximum relative spiral surface density $\vert \delta\Sigma/\Sigma_0\vert \sim \mathrm{a\,\,few} \times (\beta/\Omega_\mathrm{p})^2 \cot \alpha$, where $\Omega_\mathrm{p}$ is the spiral pattern speed and $\alpha$ is its pitch angle. This result is in good agreement with recent N-body simulations of isolated stellar disks, and is consistent with observed trends from spectroscopic surveys.

In the present work, a consistent Lagrangian model that encapsulates fully kinetic ions and gyrokinetic electrons for solar wind electromagnetic turbulence is formulated. Using a consistent method, where both electrons and protons are treated with the same mathematical formalism, we derive and implement a model in which high frequency waves and kinetic electrons effects are described in a computationally cost-efficient way. To that aim, higher order Lie-transform perturbation methods applied to Hamiltonian formulation of guiding center motion are used in order to describe the dynamics of particles and fields. Furthermore, the use of a Hamiltonian formulation allow us to introduce an abelian and gauge invariant electromagnetic field theory for the closure of the system.

We present here a concept of measuring local ties between collocated GNSS and VLBI stations using the microwave technique that effectively transforms a GNSS receiver to an element of a VLBI network. This is achieved by modifying the signal chain that allows to transfer voltage of the GNSS antenna to a digitizer via a coaxial cable. We discuss the application of this technique to local tie measurement. We have performed observations with a GNSS antenna and FD-VLBA radiotelescope and detected a strong interferometric signal from both radiogalaxies and GNSS satellites.

The arrival time prediction of Coronal mass ejections (CMEs) is an area of active research. Many methods with varying levels of complexity have been developed to predict CME arrival. However, the mean absolute error (MAE) of predictions remains above 12 hours, even with the increasing complexity of methods. In this work we develop a new method for CME arrival time prediction that uses magnetohydrodynamic simulations involving data-constrained flux-rope-based CMEs, which are introduced in a data-driven solar wind background. We found that, for 6 CMEs studied in this work, the MAE in arrival time was ~8 hours. We further improved our arrival time predictions by using ensemble modeling and comparing the ensemble solutions with STEREO-A&B heliospheric imager data. This was done by using our simulations to create synthetic J-maps. A machine learning (ML) method called the lasso regression was used for this comparison. Using this approach, we could reduce the MAE to ~4 hours. Another ML method based on the neural networks (NNs) made it possible to reduce the MAE to ~5 hours for the cases when HI data from both STEREO-A&B were available. NNs are capable of providing similar MAE when only the STEREO-A data is used. Our methods also resulted in very encouraging values of standard deviation (precision) of arrival time. The methods discussed in this paper demonstrate significant improvements in the CME arrival time predictions. Our work highlights the importance of using ML techniques in combination with data-constrained magnetohydrodynamic modeling to improve space weather predictions.

Scintillator detector response modelling has become an essential tool in various research fields such as particle and nuclear physics, astronomy or geophysics. Yet, due to the system complexity and the requirement for accurate electron response measurements, model inference and calibration remains a challenge. Here, we propose Compton edge probing to perform non-proportional scintillation model (NPSM) inference for inorganic scintillators. We use laboratory-based gamma-ray radiation measurements with a NaI(Tl) scintillator to perform Bayesian inference on a NPSM. Further, we apply machine learning to emulate the detector response obtained by Monte Carlo simulations. We show that the proposed methodology successfully constrains the NPSM and hereby quantifies the intrinsic resolution. Moreover, using the trained emulators, we can predict the spectral Compton edge dynamics as a function of the parameterized scintillation mechanisms. The presented framework offers a novel way to infer NPSMs for any inorganic scintillator without the need for additional electron response measurements.

We report the existence of novel static spherical black-hole solutions in a vector-tensor gravitational theory called the bumblebee gravity model which extends the Einstein-Maxwell theory by allowing the vector to nonminimally couple to the Ricci curvature tensor. A test of the solutions in the strong-field regime is performed for the first time using the recent observations of the supermassive black-hole shadows in the galaxy M87 and the Milky Way from the Event Horizon Telescope Collaboration. The parameter space is found largely unexcluded and more experiments are needed to fully bound the theory.

Astronomical observations suggest that the Universe may be anisotropic on the largest scales. In order to model this situation, we develop a new approach to cosmology that allows for large-scale anisotropy to emerge from the growth of non-linear structure. This is achieved by decomposing all relevant fields with respect to a preferred space-like direction, and then averaging the resulting scalar quantities over spatial domains. Our approach allows us to derive a set of large-scale effective field equations that govern the dynamics of any emergent large-scale anisotropy, and which (up to back-reaction terms) take the form of the field equations of the locally rotationally symmetric Bianchi cosmologies. We apply our approach to the dust-filled Farnsworth solutions, which are an interesting set of exact cosmological models that allow for both anisotropic expansion and large-scale bulk flow.

The Single electron Sensitive Read Out (SiSeRO) is a novel on-chip charge detector output stage for charge-coupled device (CCD) image sensors. Developed at MIT Lincoln Laboratory, this technology uses a p-MOSFET transistor with a depleted internal gate beneath the transistor channel. The transistor source-drain current is modulated by the transfer of charge into the internal gate. At Stanford, we have developed a readout module based on the drain current of the on-chip transistor to characterize the device. In our earlier work, we characterized a number of first prototype SiSeROs with the MOSFET transistor channels at the surface layer. An equivalent noise charge (ENC) of around 15 electrons root mean square (RMS) was obtained. In this work, we examine the first buried-channel SiSeRO. We have achieved substantially improved noise performance of around 4.5 electrons root mean square (RMS) and a full width half maximum (FWHM) energy resolution of 132 eV at 5.9 keV, for a readout speed of 625 kpixel/s. We also discuss how digital filtering techniques can be used to further improve the SiSeRO noise performance. Additional measurements and device simulations will be essential to further mature the SiSeRO technology. This new device class presents an exciting new technology for the next-generation astronomical X-ray telescopes requiring fast, low-noise, radiation-hard megapixel imagers with moderate spectroscopic resolution.

We consider the back-reaction of cosmological fluctuations on the local expansion rate averaged over a space-like hypersurface of constant value of a clock field. We show that in the infrared limit, the fluctuations lead to a decrease in the average expansion rate, measured at a fixed value of the clock field, compared to what would be obtained in a homogeneous universe. We work in the context of Einstein gravity coupled to perfect fluid matter.

\item[Background] The nuclear responses for antineutrinos associated with double beta decays (DBDs) and astro-antineutrino interactions are studied by measuring ordinary muon capture (OMC) rates. \item[Purpose]The experimental studies of absolute OMC rates and their mass number dependence for $^{100}$Mo and the natural Mo are currently of interest in astro-antineutrinos and DBDs. \item[Method]The OMC rates were obtained experimentally by measuring the time spectrum of the trapped muon's decay into electrons to obtain the half-lives of the trapped muons. \item[Results]The OMC rate for the enriched isotope of $^{100}$Mo is $\Lambda$($^{100}$Mo)=(7.07$\pm$0.32)$\times10^{6}$ s$^{-1}$, while that for the natural Mo is $\Lambda$($^{\rm nat}$Mo)=(9.66$\pm$0.44)$\times10^{6}$ s$^{-1}$, i.e., $\Lambda$($^{100}$Mo) is about 27$\%$ of $\Lambda$($^{\rm nat}$Mo), reflecting the blocking effect of the excess neutrons for the proton-to-neutron transformation in OMC. The present experimental observation is consistent with the predictions using Goulard-Primakoff's (GPs) and Primakoff's (Ps) empirical equations. \item[Conclusions] The absolute OMC rates for $^{100}$Mo and $^{\rm nat}$Mo were measured. The large neutron excess in $^{100}$Mo gives a much lower OMC rate than $^{\rm nat}$Mo. On both $^{100}$Mo and $^{\rm nat}$Mo, consistent OMC rates with the GP and P values are observed.

Inspired by recent interest in geometric deep learning, this work generalises the recently developed Slepian scale-discretised wavelets on the sphere to Riemannian manifolds. Through the sifting convolution, one may define translations and, thus, convolutions on manifolds - which are otherwise not well-defined in general. Slepian wavelets are constructed on a region of a manifold and are therefore suited to problems where data only exists in a particular region. The Slepian functions, on which Slepian wavelets are built, are the basis functions of the Slepian spatial-spectral concentration problem on the manifold. A tiling of the Slepian harmonic line with smoothly decreasing generating functions defines the scale-discretised wavelets; allowing one to probe spatially localised, scale-dependent features of a signal. By discretising manifolds as graphs, the Slepian functions and wavelets of a triangular mesh are presented. Through a wavelet transform, the wavelet coefficients of a field defined on the mesh are found and used in a straightforward thresholding denoising scheme.

How the growth of large-scale magnetic fields depends on microphysical transport has long been a focus of magnetic dynamo theory, and helical dynamo simulations have shown that the time to reach saturation in closed systems depends on the magnetic Reynolds number Rm. Because this would be too long for many high-Rm astrophysical systems, here we tackle the long-standing question of how much Rm-independent growth occurs earlier. From modest-Rm numerical simulations, we identify and explain a regime when the large-scale field grows independently of Rm, but to a magnitude that decreases with Rm. For plausible magnetic spectra however, the same analysis predicts the growth in this regime to be Rm-independent and provides a substantial lower bound for the field strength as Rm$\to\infty$. The results provide renewed optimism for the relevance of closed dynamos and pinpoint how modest Rm and hyper-diffusive simulations can cause misapprehension of Rm$\to\infty$ behavior.

We take into account the dynamics of three types of models of rotating galaxies in polar coordinates in a rotating frame. Due to non-axisymmetric potential perturbations, the angular momentum varies with time, and the kinetic energy depends on the momenta and spatial coordinate. The existing explicit force-gradient symplectic integrators are not applicable to such Hamiltonian problems, but the recently extended force-gradient symplectic methods proposed in a previous work are. Numerical comparisons show that the extended force-gradient fourth-order symplectic method with symmetry is superior to the standard fourth-order symplectic method but inferior to the optimized extended force-gradient fourth-order symplectic method in accuracy. The optimized extended algorithm with symmetry is used to explore the dynamical features of regular and chaotic orbits in these rotating galaxy models. The gravity effects and the degree of chaos increase with an increase of the number of the radial terms in the series expansions of the potential. There are similar dynamical structures of regular and chaotical orbits in the three types of models for the same number of the radial terms in the series expansions, energy and initial conditions.

Equatorial ionospheric irregularities have been studied in the past and have produced interesting insights into ionospheric physics and processes. Here, we present the initial results of a long-term study of the ionosphere near the Equatorial Ionization Anomaly (EIA) using Navigation with the Indian Constellation (NavIC). We have characterized the ionospheric irregularities in terms of the power spectral density at different dynamical frequencies. The formalism is similar to as suggested by earlier works using the phase screen modeling of the ionosphere. The observations of the C/N 0 (dB-Hz) variation have been taken by utilizing the L5 (1176.45 MHz) signal of NavIC over Indore located near the northern crest of EIA. We show some initial results as a proof of concept study from a single day (December 4, 2017) of scintillation observations. This is a first-of-its-kind study in this region with NavIC. From the power spectral density analysis, we have demonstrated that NavIC is capable of detecting such irregularities over long periods over this region and has implications for forecasting such events in the future.

Capturing the image of the shadow cast by the event horizon of an illuminated black hole is, at the most basic level, an experiment of extreme light deflection in a strongly curved spacetime. As such, the properties of an imaged shadow can be used to probe the general relativistic Kerr nature of astrophysical black holes. As an example of this prospect it is commonly asserted that a shadow can test the validity of the theory's famous `no hair theorem' for the black hole's mass and spin multipole moments. In this paper we assess this statement by calculating the shadow's equatorial radius in spacetimes with an arbitrary multipolar structure and within a slow rotation approximation. We find that when moments higher than the quadrupole are taken into account, the shadow acquires a high degree of degeneracy as a function of the deviation from the Kerr multipole moments. The results of our analysis suggest that dark objects with strongly non-Kerr multipolar structure could nevertheless produce a Kerr-like shadow with its characteristic quasi-circular shape.

In this paper, preliminary results from the artificial neural network (ANN) based model developed at IIT Indore has been presented. One year's hourly total electron content (TEC) database has been created from the International Reference Ionosphere (IRI) 2016 model. For the first time, a reverse problem has been addressed, wherein the training has been performed for predicting the three indices: 13-month running sunspot number, ionospheric index, and daily solar radio flux also called targets to the network when hourly TEC values are the inputs. The root mean square errors (RMSEs) of these targets have been compared and minimized after several training of the dataset using different sets of combinations. Unknown data fed to the network yielded 0.99%, 3.12%, and 0.90% errors for Rz12, IG12, and F10.7 radio flux, respectively, thus signifying ~97% prediction accuracy of the model.

Intense geomagnetic storms can have a strong impact on the signals (termed ionospheric scintillations) emitted by any global navigation satellite system (GNSS). The paper reports the first studies of scintillations at the Indore region on the NavIC signals due to the impact of the intense geomagnetic storm event reported on September 8, 2017, at 0151 and 1304 UT. The variation of the planetary indices as well as the DST index which dropped to a value of -124 nT on September 8, 2017, indicates the occurrence of an intense geomagnetic storm on September 8, 2017. The observations presented are carried out at Indore, which is located at the equatorial anomaly crest. The S4 index measurements of colocated GNSS receiver showed values of 0.5 or above on the disturbed day between 15 and 18 UT. The analysis presented clearly signifies the degradation of the carrier to noise measurements of the NavIC L5 signal during the same time, which in turn affected the positional accuracy of NavIC, an important consideration for performance.

The scale invariant theory is preserving the fundamental physical properties of General Relativity, while enlarging the group of invariances subtending gravitation theory (Dirac1973; Canuto et al.1977). The Scale Invariant Vacuum (SIV) theory assumes, as gauging condition, that:"The macroscopic empty space is scale invariant, homogeneous and isotropic". Some basic properties in Weyl's Integrable Geometry and cotensor calculus are examined in relation with scalar-tensor theories. Possible scale invariant effects are strongly reduced by matter density, both at the cosmological and local levels. The weak feld limit of SIV tends to MOND, when the scale factor is taken as constant, an approximation valid (<1%) over the last 400 Myr. A better understanding of the a0-parameter is obtained: it corresponds to the equilibrium point of the Newtonian and SIV dynamical acceleration. Parameter a0 is not a universal constant, it depends on the density and age of the Universe. As MOND is doing, SIV theory avoids the call to dark matter, moreover the cosmological models predict accelerated expansion.

Investigations on the resonance in the collisional flavor instability (CFI) of neutrinos, which were reported recently, are reported. We show that the resonance occurs not only for the isotropy-preserving modes as pointed out in the previous work but also for the isotropy-breaking modes and that it enhances the growth rate of CFI by orders of magnitude. Employing the linear analysis and nonlinear numerical simulations in the two-flavor scheme and under the relaxation approximation for the collision term, we discuss the criterion for the resonance, its effect on the nonlinear evolution as well as the influences of homogeneity-breaking (k \ne 0) perturbations as well as of anisotropy in the background on the resonance. We will also touch the cohabitation of the resonance with the fast flavor conversion (FFC).

MICROSCOPE's final results [1] report no violation of the Weak Equivalence Principle (Universality of Free Fall) for Pt and Ti test masses quantified by an Eotvos parameter of about $10^{-15}$, an improvement by about two orders of magnitude over the best ground tests. The measurement is limited by random noise with $1/\sqrt{\nu}$ frequency dependence attributed to thermal noise from internal damping occurring in the grounding wires. From information provided in [2] and the physics of internal damping we calculate the differential acceleration noise spectral density at the signal frequency, and show it varies widely between experiment sessions. Such large variations are inexplicable if translated into physical quantities such as the quality factor. While calibrations interspersed with measurement sessions may cause some such changes, they cannot explain jumps between consecutive sessions without recalibration. A potential explanation is conjectured related to a fluctuating zero depending on measurement initialization errors. The experiment was plagued with "glitches" -- anomalous acceleration spikes related to radiation from the Earth -- injecting significant power at the signal frequency and its harmonics. We argue that the procedure used to eliminate the glitches and fill the gaps with artificially reconstructed data leaves residuals at the critical frequencies that may affect the result of the experiment. We propose that an analysis of the time domain data as a function of the orientation of the sensor would be tolerant of data gaps without affecting the sensitivity. Future experiments aiming to exploit the full potential of space must resolve these issues, rely solely on measured data, and, more generally, readdress the experiment design.

The properties of collisionless shocks are frequently assessed in the magnetohydrodynamics (MHD) model. Yet, in a collisionless plasma, an ambient magnetic field can sustain a stable anisotropy in the upstream or the downstream, resulting in a departure from the MHD predicted behavior. We present a model allowing to derive the downstream anisotropy, hence the shock density jump, in terms of the upstream quantities. For simplicity, the case of a parallel shock in pair plasma is considered. Contrary to previous works where the upstream was assumed isotropic, here the upstream anisotropy $A=T_\perp/T_\parallel$ is a free parameter. The strong sonic shock regime is formally identical to the isotropic upstream case. Yet, for intermediate sonic Mach numbers, a variety of behaviors appear as a result of the anisotropy of the upstream.