Papers for Tuesday, Feb 07 2023

We present GLASS, the Generator for Large Scale Structure, a new code for the simulation of galaxy surveys for cosmology, which iteratively builds a light cone with matter, galaxies, and weak gravitational lensing signals as a sequence of nested shells. This allows us to create deep and realistic simulations of galaxy surveys at high angular resolution on standard computer hardware and with low resource consumption. GLASS also introduces a new technique to generate transformations of Gaussian random fields (including lognormal) to essentially arbitrary precision, an iterative line-of-sight integration over matter shells to obtain weak lensing fields, and flexible modelling of the galaxies sector. We demonstrate that GLASS readily produces simulated datasets with per cent-level accurate two-point statistics of galaxy clustering and weak lensing, thus enabling simulation-based validation and inference that is limited only by our current knowledge of the input matter and galaxy properties.

The flux of $\gamma$ rays is measured with unprecedented accuracy by the $Fermi$ Large Area Telescope from 100 MeV to almost 1 TeV. In the future, the Cherenkov Telescope Array will have the capability to measure photons up to 100 TeV. To accurately interpret this data, precise predictions of the production processes, specifically the cross section for the production of photons from the interaction of cosmic-ray protons and helium with atoms of the ISM, are necessary. In this study, we determine new analytical functions describing the Lorentz-invariant cross section for $\gamma$-ray production in hadronic collisions. We utilize the limited total cross section data for $\pi^0$ production channels and supplement this information by drawing on our previous analyses of charged pion production to infer missing details. In this context, we highlight the need for new data on $\pi^0$ production. Our predictions include the cross sections for all production channels that contribute down to the 0.5% level of the final cross section, namely $\eta$, $K^+$, $K^-$, $K^0_S$ and $K^0_L$ mesons as well as $\Lambda$, $\Sigma$, and $\Xi$ baryons. We determine the total differential cross section $d\sigma(p+p\rightarrow \gamma+X)/dE_{\gamma}$ from 10 MeV to 100 TeV with an uncertainty of 10% below 10 GeV of $\gamma$-ray energies, increasing to 20% at the TeV energies. We provide numerical tables and a script for the community to access our energy-differential cross sections.

The possibility of an expanding decelerating Universe in the distant future is investigated in the context of a quintessence scalar field cosmology [1]. Such a conceivable evolution is tested against SNe Ia and $H(z)$ data, and also through a model independent method based on Gaussian Processes. The scalar field model is an extension of the exponential Ratra-Peebles quintessence whose potential [2] depends on a pair of parameters ($\alpha, \lambda)$ and predicts a decelerated expansion in the future. The model analysis is updated with the most recent SNe Ia and $H(z)$ data thereby obtaining $H_0 =68.6\pm3.7$ km/s/Mpc, $\Omega_{\Phi0} = 0.735^{+0.083}_{-0.069} $, $\alpha < 6.56$ and $\lambda< 0.879 $ (at $2\sigma$ c.l.). It is found that the model allows for a future deceleration both for $H(z)$ and SNe Ia data. In the (model-independent) Gaussian Processes analysis, however, future deceleration is allowed only in the case of $H(z)$ data.

We reconstruct the extra-galactic gamma-ray source-count distribution, or $dN/dS$, of resolved and unresolved sources by adopting machine learning techniques. Specifically, we train a convolutional neural network on synthetic 2-dimensional sky-maps, which are built by varying parameters of underlying source-counts models and incorporate the Fermi-LAT instrumental response functions. The trained neural network is then applied to the Fermi-LAT data, from which we estimate the source count distribution down to flux levels a factor of 50 below the Fermi-LAT threshold. We perform our analysis using 14 years of data collected in the $(1,10)$ GeV energy range. The results we obtain show a source count distribution which, in the resolved regime, is in excellent agreement with the one derived from catalogued sources, and then extends as $dN/dS \sim S^{-2}$ in the unresolved regime, down to fluxes of $5 \cdot 10^{-12}$ cm$^{-2}$ s$^{-1}$. The neural network architecture and the devised methodology have the flexibility to enable future analyses to study the energy dependence of the source-count distribution.

Measurement of black hole mass for low-$z$ ($z<0.8$) Active Galactic Nuclei (AGNs) is difficult due to the strong contribution from host galaxy stellar light necessitating detailed spectral decomposition to estimate the AGN luminosity. Here, we present an empirical relation to estimate host galaxy stellar luminosity from the optical spectra of AGNs at $z\leq 0.8$. The spectral data were selected from the fourteenth data release of the Sloan Digital Sky Survey (SDSS-DR14) quasar catalog having a signal-to-noise ratio at 5100 \AA (SNR$_{5100}$) $>$10 containing 11415 quasars. The median total luminosity (log ($L_\text{total}$/[erg s$^{-1}$])), stellar luminosity (log ($L_\text{star}$/[erg s$^{-1}$])), and AGN continuum luminosity (log ($L_\text{cont}$/[erg s$^{-1}$])) in our sample are 44.52, 44.06, and 44.30, respectively. We fit the AGN power-law continuum, host galaxy, and iron blend contribution, simultaneously over the entire available spectrum. We found the host galaxy fraction to anti-correlate with continuum luminosity and can be well-represented by a polynomial function, which can be used to correct the stellar light contribution from AGN spectra. We also found anti-correlation between host galaxy fraction and iron strength, Eddington ratio, and redshift. The empirical relation gives comparable results of host-fraction with the image decomposition method.

We present new VLA observations, between 6.8mm and 66mm, of the edge-on Class~I disk IRAS04302+2247. Observations at 6.8mm and 9.2mm lead to the detection of thermal emission from the disk, while shallow observations at the other wavelengths are used to correct for emission from other processes. The disk radial brightness profile transitions from broadly extended in previous ALMA 0.9mm and 2.1mm observations to much more centrally brightened at 6.8mm and 9.2mm, which can be explained by optical depth effects. The radiative transfer modeling of the 0.9mm, 2.1mm, and 9.2mm data suggests that the grains are smaller than 1cm in the outer regions of the disk and allows us to obtain the first lower limit for the scale height of grains emitting at millimeter wavelengths in a protoplanetary disk. We find that the millimeter dust scale height is between 1au and 6au at a radius 100au from the central star, while the gas scale height is estimated to be about 7au, indicating a modest level of settling. The estimated dust height is intermediate between less evolved Class 0 sources, that are found to be vertically thick, and more evolved Class II sources, which show a significant level of settling. This suggests that we are witnessing an intermediate stage of dust settling.

The orbits of asteroids from the MPC catalogue of May 31, 2020 with absolute magnitudes H < 16m, in the 3:1, 5:2 and 2:1 mean motion resonances (MMRs) with Jupiter were selected. The number of the orbits in the 2:1 MMR is dozens times greater than in the other two resonances. There are fragments of parent bodies of neighbour asteroid families, in particular the Themis family, among bodies in the 2:1 MMR. Numerical calculations were performed to evaluate the evolution of the selected orbits over hundreds of thousand years. Perturbations from all eight major planets and the relativistic effects of orbital perihelion displacement were taken into account. For all orbits in the 3:1 and 5:2 MMRs an increase in the orbit eccentricities, which are sufficient for the bodies to approach Mars, was obtained. In the 2:1 MMR, a sufficient increase in the orbit eccentricities was not detected. An increase in orbit eccentricities in this resonance can occur due to the action of non-gravitational effects (NGEs). The action of the Yarkovsky effect can explain the exit of an asteroid with a size of 5 km from the 2:1 MMR over a period about 1 billion years or more. More than 2 billion years ago, there were dozens of bodies over 50 km in size in the 2:1 gap. To form the gap in the 2:1 resonance at the very beginning, the physical conditions in the asteroid belt had to be significantly different from the today ones. In particular, the intensity of the solar radiation in the early Solar system could be much higher as compared to the today one.

Artificial intelligence methods show great promise in increasing the quality and speed of work with large astronomical datasets, but the high complexity of these methods leads to the extraction of dataset-specific, non-robust features. Therefore, such methods do not generalize well across multiple datasets. We present a universal domain adaptation method, \textit{DeepAstroUDA}, as an approach to overcome this challenge. This algorithm performs semi-supervised domain adaptation and can be applied to datasets with different data distributions and class overlaps. Non-overlapping classes can be present in any of the two datasets (the labeled source domain, or the unlabeled target domain), and the method can even be used in the presence of unknown classes. We apply our method to three examples of galaxy morphology classification tasks of different complexities ($3$-class and $10$-class problems), with anomaly detection: 1) datasets created after different numbers of observing years from a single survey (LSST mock data of $1$ and $10$ years of observations); 2) data from different surveys (SDSS and DECaLS); and 3) data from observing fields with different depths within one survey (wide field and Stripe 82 deep field of SDSS). For the first time, we demonstrate the successful use of domain adaptation between very discrepant observational datasets. \textit{DeepAstroUDA} is capable of bridging the gap between two astronomical surveys, increasing classification accuracy in both domains (up to $40\%$ on the unlabeled data), and making model performance consistent across datasets. Furthermore, our method also performs well as an anomaly detection algorithm and successfully clusters unknown class samples even in the unlabeled target dataset.

Solar flares are efficient particle accelerators, with a substantial fraction of the energy released manifesting as non-thermal particles. While the role that non-thermal electrons play in transporting flare energy is well studied, the properties and importance of non-thermal protons is rather less well understood. This is in large part due to the paucity of diagnostics, particularly at the lower-energy (deka-keV) range of non-thermal proton distributions in flares. One means to identify the presence of deka-keV protons is by an effect originally described by \cite{1976ApJ...208..618O}. In the Orrall-Zirker effect, non-thermal protons interact with ambient neutral hydrogen, and via charge exchange produce a population of energetic neutral atoms (ENAs) in the chromosphere. These ENAs subsequently produce an extremely redshifted photon in the red wings of hydrogen spectral lines. We revisit predictions of the strength of this effect using modern interaction cross-sections, and numerical models capable of self-consistently simulating the flaring non-equilibrium ionization stratification, and the non-thermal proton distribution (and, crucially, their feedback on each other). We synthesize both the thermal and non-thermal emission from \lya\ and \lyb, the most promising lines that may exhibit a detectable signal. These new predictions are are weaker and more transient than prior estimates, but the effects should be detectable in fortuitous circumstances. We degrade the \lyb\ emission to the resolution of the Spectral Imaging of the Coronal Environment (SPICE) instrument on board Solar Orbiter, demonstrating that though likely difficult, it should be possible to detect the presence of non-thermal protons in flares observed by SPICE.

Galaxies show different halo scaling relations such as the Radial Acceleration Relation, the Mass Discrepancy Acceleration Relation (MDAR) or the dark matter Surface Density Relation (SDR). At difference with traditional studies using phenomenological $\Lambda$CDM halos, we analyze the above relations assuming that dark matter (DM) halos are formed through a Maximum Entropy Principle (MEP) in which the fermionic (quantum) nature of the DM particles is dully accounted for. For the first time a competitive DM model based on first physical principles, such as (quantum) statistical-mechanics and thermodynamics, is tested against a large data-set of galactic observables. In particular, we compare the fermionic DM model with empirical DM profiles: the NFW model, a generalized NFW model accounting for baryonic feedback, the Einasto model and the Burkert model. For this task, we use a large sample of 120 galaxies taken from the Spitzer Photometry and Accurate Rotation Curves (SPARC) data-set, from which we infer the DM content to compare with the models. We find that the Radial Acceleration Relation and MDAR are well explained by all the models with comparable accuracy, while the fits to the individual rotation curves, in contrast, show that cored DM halos are statistically preferred with respect to the cuspy NFW profile. However, very different physical principles justify the flat inner halo slope in the most favored DM profiles: while generalized NFW or Einasto models rely on complex baryonic feedback processes, the MEP scenario involves a quasi-thermodynamic equilibrium of the DM particles.

Glycolaldehyde (HOCH2CHO) is the most straightforward sugar detected in the Interstellar Medium (ISM) and participates in the formation pathways of molecules fundamental to life, red such as ribose and derivatives. Although detected in several regions of the ISM, its formation route is still debated and its abundance cannot be explained only by reactions in the gas phase. This work explores a new gas-phase formation mechanism for glycolaldehyde and compares the energy barrier reduction when the same route happens on the surface of amorphous ices. The first step of the mechanism involves the formation of a carbon-carbon bond between formaldehyde (H2CO) and the formyl radical (HCO), with an energy barrier of 27 kJ mol-1 (gas-phase). The second step consists of barrierless hydrogen addition. Density functional calculations under periodic boundary conditions were applied to study this reaction path on 10 different amorphous ice surfaces through an Eley-Rideal type mechanism. It was found that the energy barrier is reduced on average by 49 per cent, leading in some cases to a 100 per cent reduction. The calculated adsorption energy of glycolaldehyde suggests that it can be promptly desorbed to the gas phase after its formation. This work, thus contributes to explaining the detected relative abundances of glycolaldehyde and opens a new methodological framework for studying the formation routes for Complex Organic Molecules (COMs) in interstellar icy grains.

We present an estimate of the bulk flow in a volume of radii $150-200h^{-1}$Mpc using the minimum variance (MV) method with data from the CosmicFlows-4 (CF4) catalog. The addition of new data in the CF4 has resulted in an increase in the estimate of the bulk flow in a sphere of radius $150h^{-1}$Mpc relative to the CosmicFlows-3 (CF3). This bulk flow has less than a $0.03\%$ chance of occurring in the Standard Cosmological Model ($\Lambda$CDM) with cosmic microwave background derived parameters. Given that the CF4 is deeper than the CF3, we were able to use the CF4 to accurately estimate the bulk flow on scales of $200h^{-1}$Mpc (equivalent to 266 Mpc for Hubble constant $H_o=75$ km/s/Mpc) for the first time. This bulk flow is in even greater tension with the Standard Model, having less than $0.003\%$ probability of occurring. To estimate the bulk flow accurately, we introduce a novel method to calculate distances and velocities from distance moduli that is unbiased and accurate at all distances. Our results are completely independent of the value of $H_o$.

We measure the 3D kinematic structures of the young stars within the central 0.5 parsec of our Galactic Center using the 10 m telescopes of the W.~M.~Keck Observatory over a time span of 25 years. Using high-precision measurements of positions on the sky, and proper motions and radial velocities from new observations and the literature, we constrain the orbital parameters for each young star. Our results show two statistically significant sub-structures: a clockwise stellar disk with 18 candidate stars, as has been proposed before, but with an improved disk membership; a second, almost edge-on plane of 10 candidate stars oriented East-West on the sky that includes at least one IRS 13 star. We estimate the eccentricity distribution of each sub-structure and find that the clockwise disk has <$e$> = 0.39 and the edge-on plane has <$e$> = 0.68. We also perform simulations of each disk/plane with incompleteness and spatially-variable extinction to search for asymmetry. Our results show that the clockwise stellar disk is consistent with a uniform azimuthal distribution within the disk. The edge-on plane has an asymmetry that cannot be explained by variable extinction or incompleteness in the field. The orientation, asymmetric stellar distribution, and high eccentricity of the edge-on plane members suggest that this structure may be a stream associated with the IRS 13 group. The complex dynamical structure of the young nuclear cluster indicates that the star formation process involved complex gas structures and dynamics and is inconsistent with a single massive gaseous disk.

The progenitors for many types of supernovae (SNe) are still unknown, and an approach to diagnose their physical origins is to investigate the light curve brightness and shape of a large set of SNe. However, it is often difficult to compare and contrast the existing sample studies due to differences in their approaches and assumptions, for example in how to eliminate host galaxy extinction, and this might lead to systematic errors when comparing the results. We therefore introduce the Hybrid Analytic Flux FittEr for Transients (haffet), a Python-based software package that can be applied to download photometric and spectroscopic data for transients from open online sources, derive bolometric light curves, and fit them to semi-analytical models for estimation of their physical parameters. In a companion study, we have investigated a large collection of SNe Ib and Ic observed with the Zwicky Transient Facility (ZTF) with haffet, and here we detail the methodology and the software package to encourage more users. As large-scale surveys such as ZTF and LSST continue to discover increasing numbers of transients, tools such as haffet will be critical for enabling rapid comparison of models against data in statistically consistent, comparable and reproducable ways. Additionally, haffet is created with a Graphical User Interface mode, which we hope will boost the efficiency and make the usage much easier.

The distribution of matter that is measured through galaxy redshift and peculiar velocity surveys can be harnessed to learn about the physics of dark matter, dark energy, and the nature of gravity. To improve our understanding of the matter of the Universe, we can reconstruct the full density and velocity fields from the galaxies that act as tracer particles. We use a convolutional neural network, a V-net, trained on numerical simulations of structure formation to reconstruct the density and velocity fields. We find that, with detailed tuning of the loss function, the V-net could produce better fits to the density field in the high-density and low-density regions, and improved predictions for the amplitudes of the velocities. We also find that the redshift-space distortions of the galaxy catalogue do not significantly contaminate the reconstructed real-space density and velocity field. We estimate the velocity field $\beta$ parameter by comparing the peculiar velocities of mock galaxy catalogues to the reconstructed velocity fields, and find the estimated $\beta$ values agree with the fiducial value at the 68\% confidence level.

We report the results of intensive monitoring of the variability in the H$\alpha$ line for two F-type stars, $\tau$ Boo and $\upsilon$ And, during the last four years 2019-2022, in order to investigate their stellar magnetic activity. The 4-year H$\alpha$ line intensity data taken with the 1.88-m reflector at Okayama Branch Office, Subaru Telescope, shows the existence of a possible $\sim$ 123-day magnetic activity cycle of $\tau$ Boo. The result of the H$\alpha$ variability as another tracer of the magnetic activity on the chromosphere is consistent with previous studies of the Ca II H&K line and suggests that the magnetic activity cycle is persisted in $\tau$ Boo. For $\upsilon$ And, we suggest a quadratic long-term trend in the H$\alpha$ variability. Meanwhile, the short-term monitoring shows no significant period corresponding to specific variations likely induced by their hot-Jupiter in both cases ($\approx$ 3.31 and 4.62 days, respectively). In this H$\alpha$ observation, we could not find any signature of the Star-Planet Magnetic Interaction. It is speculated that the detected magnetic activity variability of the two F-type stars is related to the stellar intrinsic dynamo.

We briefly review progress in developing a pathway to nonlinear astereoseismology, both from theoretical and observational aspects. As predicted by the theory of weak nonlinear interactions between resonant modes, their amplitude and frequency can be modulated according to various kinds of patterns. However, those subtle modulations could hardly be well characterized from ground-based photometric monitoring. The {\sl Kepler} spacecraft offered a new window to find clear-cut evidence of well-determined amplitude and frequency modulations, leading to the first discoveries of such variations in pulsating white dwarf and hot B subdwarf stars. Following that direction, a systematic survey of oscillation mode properties in compact pulsators monitored by {\sl Kepler} suggests that mode variability is likely a common phenomenon, which remain unaccounted for by standard linear non-radial pulsation theory. To reach this conclusion firmly, the survey has now been extended to a larger context including compact stars observed by K2 and TESS. We expect that this extended survey will help to constrain key parameters governing weak nonlinear effects in stellar oscillations.

Context. To date, stellar activity is one of the main limitations in detecting small exoplanets via transit photometry. Since this activity is enhanced in young stars, traditional filtering algorithms may severely under-perform in detecting such exoplanets. Aims. This paper aims to compare the relative performances of four algorithms developed by independent research groups specifically for the filtering of activity in the light curves (LCs) of young active stars, prior to the search for planetary transit signals: Notch and LOCoR(N&L), Young Stars Detrending(YSD), K2 Systematics Correction(K2SC) and VARLET. We include in the comparison also the two best-performing algorithms implemented in Wotan, namely the Tukey's biweight and the Huber Spline. Methods. We performed a series of injection-retrieval tests of planetary transits of different types, from Jupiter down to Earth-sized planets, moving both on circular and eccentric orbits. The tests were carried out over 100 simulated LCs of both quiet and active solar-like stars that will be observed by the ESA space telescope PLATO. Results. We found that N&L is the best choice in many cases, since it misses the lowest number of transits. However, it under-performs if the planetary orbital period closely matches the stellar rotation period, especially in the case of small planets for which the biweight and VARLET algorithms work better. For LCs with a large number of data, the combined results of YSD and Huber Spline yield the highest recovery percentage. Filtering algorithms allow us to get a very precise estimate of the orbital period and the mid-transit time of the detected planets, while the planet-to-star radius is under-estimated most of the time, especially in the case of grazing transits or eccentric orbits. A refined filtering taking into account the presence of the planet is compulsory for a proper planetary characterization.

Some Miras -- long-period variables in late evolutionary stages -- have meandering pulsation periods and light curve asymmetries, the causes of which are still unclear. We aim to understand better the origin of these phenomena by investigating a sample of solar-neighbourhood Miras. We characterised this group of stars and related their variability characteristics to other stellar parameters. We analysed observations from several databases to obtain light curves with maximum time span and temporal coverage for a sample of 548 Miras. We determined their pulsation period evolution over a time span of many decades, searched for changes in the periods, and determined the amplitude of the period change. We also analysed the Fourier spectra with respect to possible secondary frequency maxima. The sample was divided into two groups with respect to the presence of light curve bumps. IR colours and indicators of the third dredge-up were collected to study the sample stars' mass-loss and deep mixing properties. Our analysis revealed one new star, T~Lyn, with a continuously changing period. The group of Miras with meandering period changes is exclusively made up of M-type stars. The Fourier spectra of the meandering period Miras have no prominent additional peaks, suggesting that additional pulsation modes are not the cause of the meandering periods. We confirm that bumps are more common among S and C Miras and show, for the first time, that Miras with bumps have lower mass-loss rates than those with regular, symmetric light curves. Also Miras with meandering period changes have relatively little mass loss. We conclude that Miras with strongly changing periods or asymmetries in their light curves have relatively low dust mass-loss rates. Meandering period changes and light curve asymmetries could be connected to He-shell flashes and third dredge-up episodes.

The Cherenkov Telescope Array (CTA) is the next-generation stereoscopic system of Imaging Atmospheric Cherenkov Telescopes (IACTs). In IACTs, the atmosphere is used as a calorimeter to measure the energy of extensive air showers induced by cosmic gamma rays, which brings along a series of constraints on the precision to which energy can be reconstructed. The presence of clouds during observations can severely affect Cherenkov light yield, contributing to the systematic uncertainty in energy scale calibration. To minimize these systematic uncertainties, a calibration of telescopes is of great importance. For this purpose, the influence of cloud transmission and altitude on the CTA-North performance degradation using detailed Monte Carlo simulations has been investigated for the case when no actions are taken to correct the effects of clouds. Variations of instrument response functions in the presence of clouds are presented. In the presence of clouds with low transmission ($\leq 80\%$) the energy resolution is aggravated by 30\% at energies below 1 TeV, and by 10\% at higher energies. For higher transmissions, the energy resolution is degraded by less than 10\% in the whole energy range. The angular resolution varies up to 10\% depending both on the transmission and altitude of the cloud. The sensitivity of the array is most severely reduced at lower energies, even by 60\% at 40 GeV, depending on the clouds' properties. A simple semi-analytical model of sensitivity degradation has been introduced to summarize the influence of clouds on sensitivity and provide useful scaling relations.

Young supernova remnants (SNRs) shocks are believed to be the main sites of galactic cosmic rays production, showing X-ray synchrotron dominated spectra in the vicinity of their shock. While a faint thermal signature left by the shocked interstellar medium (ISM) should also be found in the spectra, proofs for such an emission in Tycho's SNR have been lacking. We perform an extended statistical analysis of the X-ray spectra of five regions behind the blast wave of Tycho's SNR using \textit{Chandra} archival data. We use Bayesian inference to perform extended parameter space exploration and sample the posterior distributions of a variety of models of interest. According to Bayes factors, spectra of all five regions of analysis are best described by composite three-component models taking into account non-thermal emission, ejecta emission and shocked ISM emission. The shocked ISM stands out the most in the Northern limb of the SNR. We find for the shocked ISM a mean electron temperature $kT_{\rm{e}}=1.00^{+1.17}_{-0.42}$ keV for all regions and a mean ionization timescale $n_{\rm e}t=1.96^{+1.18}_{-0.76}\times10^9$ cm$^{-3}$s resulting in a mean ambient density $n_{\rm e}=0.19^{+0.18}_{-0.08}$ cm$^{-3}$ around the remnant. We performed an extended analysis of the Northern limb, and show that the measured synchrotron cutoff energy is not well constrained in the presence of a shocked ISM component. Such results cannot currently be further investigated by analyzing emission lines in the 0.5-1 keV range, because of the low Chandra spectral resolution in this band. We show with simulated spectra that X-IFU future performances will be crucial to address this point.

The operation-planning of satellites, aimed at introducing a certain level of supervised automation during the execution of the operations, poses a great challenge to both designers and operators. From one side, the routine operations for an Earth Observation mission are predictable and typically repeatable; both these aspects are suitable for computerisation. On the other hand, not every non-nominal scenario can be anticipated and correctly formulated in terms of operations. Dealing with contingency presents risks which need to be addressed as early as possible, hopefully already during the operations preparation. It is also possible, however, to intervene at a later operational stage of the mission, optimising the tools already in use to support the operations execution. In this paper, having in mind the idea to improve existing processes in place at EUMETSAT, we present an algorithm able to reschedule the spacecraft's activities in case of anomaly. The main goal is to support the decision-making process while overcoming contingencies both avoiding overloading the spacecraft and planning engineers and ensuring the continuity of the mission, in particular giving the highest priority to the onboard computer memory size and data quality. We tested the method with the data of Sentinel-6, which carries the altimeter POSEIDON-4 operated by EUMETSAT, and the results are hereby presented.

We present updated measurements of the [O III] 88 $\mu$m, [C II] 158 $\mu$m, and dust continuum emission from a star-forming galaxy at $z=7.212$, SXDF-NB1006-2, by utilizing ALMA archival datasets analyzed in previous studies and datasets have not been analyzed before. The follow-up ALMA observation with higher angular resolution and sensitivity reveals a clumpy structure of the [O III] emission on a scale of $0.32-0.85\,\rm{kpc}$. We also combined all the ALMA [O III] ([C II]) datasets and update the [O III] ([C II]) detection to $5.9\,\sigma$ ($3.6\,\sigma-4.5\,\sigma$). The non-detection of [C II] with data from the REBELS large program implies the incompleteness of spectral-scan surveys using [C II] to detect galaxies with high star formation rates (SFRs) but marginal [C II] emission at high-$z$. The dust continuum at 90 $\mu$m and 160 $\mu$m remains non-detection, indicating little dust content of $<3.9\times10^{6}\,\rm{M_\odot}\,(3\sigma)$, and we obtained a more stringent constraint on the total infrared (TIR) luminosity. We update the [O III]/[C II] luminosity ratios to $10.2\pm4.7~(6.1\pm3.5$) and $20\pm12~(9.6\pm6.1$) for $4.5\,\sigma$ and $3.6\,\sigma$ [C II] detections, respectively, where ratios in the parentheses are corrected for the surface brightness dimming (SBD) effect on the extended [C II] emission. We also find a strong [C II] deficit ($0.6-1.3$ dex) between SXDF-NB1006-2 and the mean $L_{\rm{[CII]}}-\rm{SFR}$ relation of galaxies at $0<z<9$.

We have investigated the pitch angle ($\psi$) of the spiral arms of galaxies in the Hubble Space Telescope COSMOS field. The sample consists of 102 face-on galaxies with a two-armed pattern at a mean redshift $\langle z \rangle \approx 0.5$. The typical values of $\psi$ in the spiral arms of distant galaxies are shown to be close to those for nearby spiral galaxies. Within one galaxy the scatter of $\psi$ for different arms is, on average, half the mean pitch angle. In the $z$ range from 1 to 0 we have found a tendency for $\psi$ to decrease. Our analysis of the $\psi$ distributions in galaxies at different redshifts is consistent with the assumption that in most of the galaxies at $z \leq 0.5$ the spiral arms are tidal in origin or they arose from transient recurrent instabilities in their disks.

Determining which between projected local density and distance from the cluster center plays a major role in regulating morphological fractions in clusters is a longstanding debate. Reaching a definitive answer will shed light on the main physical mechanisms at play in the most extreme environments. Here we make use of the data from the OmegaWINGS survey, currently the largest survey of clusters in the local Universe extending beyond 2 virial radii from the cluster cores, to extend previous analysis outside the virial radius. Local density and clustercentric distance seems to play different roles for galaxies of different morphology: the fraction of elliptical galaxies mainly depends on local density, suggesting that their formation was linked to the primordial densities, which now correspond to the cluster cores. Only the fraction of low mass ellipticals shows an anticorrelation with clustercentric distance, suggesting a different origin for these objects. Excluding elliptical galaxies, the relative fraction of S0s and spirals instead depends on local density only far from the cluster cores, while within the virial radius their proportion is regulated by distance, suggesting that cluster specific processes halt the star formation and transform Sp galaxies into S0s. This interpretation is supported by literature results on the kinematical analysis of early and late type galaxies, according to which fast and slow rotators have distinct dependencies on halo mass and local density.

We present the first results for the dust-scattering rings of GRB 221009A, coined as the GRB of the century, as observed by the Neil Gehrels Swift satellite. We perform analysis of both time resolved observations and stacked data. The former approach enable us to study the expansion of the most prominent rings, associate their origin with the prompt X-ray emission of the GRB and determine the location of the dust layers. The stacked radial profiles increase the signal-to-noise ratio of the data and allows detection of fainter and overlapping peaks in the angular profile. We find a total of 15 dust concentrations (with hints of even more) that span about 25 kpc in depth and could be responsible for the highly structured X-ray angular profiles. By comparing the relative scattered fluxes of the five most prominent rings we show that the layer with the largest amount of dust is located at about 1.1 kpc away from us. We finally compare the location of the dust layers with results from experiments that study the 3D structure of our Galaxy via extinction or CO radio observations, and highlight the complementarity of dust X-ray tomography to these approaches.

We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40~hours) NIRCam observations reach a $F277W$ magnitude of 32 at the $5\sigma$ level, more than 2 magnitudes deeper than JWST public datasets already analyzed to find high redshift galaxies. We select a sample of 45 $z>8$ galaxy candidates based on their dropout nature in the $F115W$ and/or $F150W$ filters, a high probability for their photometric redshifts, estimated with three different codes, being at $z>8$, good fits based on $\chi^2$ calculations, and predominant solutions compared to $z<8$ alternatives. We find mild evolution in the luminosity function from $z\sim13$ to $z\sim8$, i.e., only a small increase in the average number density of $\sim$0.2~dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values $\alpha=-2.3\pm0.2$ and $M^*=-20.8\pm0.2$~mag. Comparing our results with the predictions of a wide range of state-of-the-art galaxy evolution models, we find two main results: (1) a slower increase with time in the cosmic star formation rate density compared to a steeper rise predicted by models; (2) nearly a factor of 10 higher star formation activity concentrated in scales around 2~kpc in galaxies with stellar masses $\sim10^8$~M$_\odot$ during the first 350~Myr of the Universe ($z\sim12$), with models matching better the observations $\sim$150~Myr later, by $z\sim9$.

The black hole at the center of our Milky Way Galaxy -- the Galactic Black Hole, or GBH -- is a rather modest representative of its class. With a mass of $4\times10^6$ solar masses, it is well over a thousand times less massive than the most extreme supermassive black holes known to be powering the most luminous quasars. Furthermore, the Galactic Black Hole has a remarkably dim accretion flow, and its luminous energy output is overwhelmed by the dense cluster of bright stars and red giants that surround it, except at radio wavelengths. However, the proximity of the GBH compensates for its restrained activity; being over 100 times closer than the next nearest supermassive black hole in a galactic nucleus, it offers us an unparalleled opportunity to observe its behavior in detail. Consequently, far more observational attention has been paid to the GBH and its entourage of stars and gas than to any other single object outside the solar system. This review covers the history of our recognition of the GBH, its presently known physical characteristics, the manifestations of its current and past activity, and the prospects for refining our knowledge with future research.

Development of novel radiation-absorbent materials and devices for millimetre and submillimetre astronomy instruments is a research area of high interest, and with substantial engineering challenges. Alongside low-profile structure and ultra-wideband performance in a wide range of angles of incidence, advanced absorbers in CMB instruments are aimed at reducing optical systematics, notably instrument polarisation, far beyond previous specifications. This paper presents an innovative design of flat thin flexible absorber operating in a wide frequency range of 80-400 GHz. The structure comprises a combination of sub-wavelength metal-mesh capacitive and inductive grids and dielectric layers, making use of the magnetic mirror concept for large bandwidth. The overall stack thickness is a quarter of the longest operating wavelength and is close to the theoretical limit stipulated by Rozanov criterion. The test device is designed to operate at 22.5deg. incidence. The iterative numerical-experimental design procedure of the new absorber is discussed in detail, as well as the practical challenges of its manufacture. A well-established mesh-filter fabrication process has been successfully employed for prototype fabrication, which ensures cryogenic operation of the hot-pressed quasi-optical devices. The final prototype, extensively tested in quasi-optical testbeds using a Fourier-transform spectrometer and a vector network analyser, demonstrated performance closely matching the finite-element analysis simulations, viz., greater than 99% absorbance for both polarisations, with only 0.2% difference, across the frequency band of 80-400 GHz. The angular stability for up to +/-10deg. has been confirmed by simulations. To the best of the authors knowledge, this is the first successful implementation of a low-profile ultrawideband metamaterial absorber for this frequency range and operating conditions.

The Przybylski star spectrum has been studied in order to search for deuterium lines. Since this star is extremely enriched in the s-process elements, which are the product of interaction between free neutrons and seed nuclei, we might as well expect to detect deuterium in this star. However, no visible spectroscopic manifestation of deuterium has been detected. Perhaps, the reason of this result is the convective "destruction" of this isotope.

In this paper, we firstly calibrate the Amati relation (the $E_{\rm p}-E_{\rm iso}$ correlation) of gamma ray bursts (GRBs) at low redshifts ($z<0.8$) via Gaussian process by using the type Ia supernovae samples from Pantheon+ under the philosophy that objects at the same redshift should have the same luminosity distance in any cosmology. As a result, this calibration derives the distance moduli of GRBs at high redshifts ($z>0.8$). For an application of these derived distance modulus of GRBs to cosmology, via Gaussian process again, a series of cosmography parameters, which describe kinematics of our Universe, up to the fifth oder, i.e. the Hubble parameter $H(z)$, the deceleration parameter $q(z)$, the jerk parameter $j(z)$, the snap parameter $s(z)$ and the lerk parameter $l(z)$, are reconstructed from the cosmic observations. The result shows that the current quality of GRBs data points are not good enough to give viable prediction of the kinematics of our Universe at high redshifts.

We present the results of an ALMA survey to identify 183 GHz H$_2$O maser emission from AGN already known to host 22 GHz megamaser systems. Out of 20 sources observed, we detect significant 183 GHz maser emission from 13; this survey thus increases the number of AGN known to host (sub)millimeter megamasers by a factor of 5. We find that the 183 GHz emission is systematically fainter than the 22 GHz emission from the same targets, with typical flux densities being roughly an order of magnitude lower at 183 GHz than at 22 GHz. However, the isotropic luminosities of the detected 183 GHz sources are comparable to their 22 GHz values. For two of our sources -- ESO 269-G012 and the Circinus galaxy -- we detect rich 183 GHz spectral structure containing multiple line complexes. The 183 GHz spectrum of ESO 269-G012 exhibits the triple-peaked structure characteristic of an edge-on AGN disk system. The Circinus galaxy contains the strongest 183 GHz emission detected in our sample, peaking at a flux density of nearly 5 Jy. The high signal-to-noise ratios achieved by these strong lines enable a coarse mapping of the 183 GHz maser system, in which the masers appear to be distributed similarly to those seen in VLBI maps of the 22 GHz system in the same galaxy and may be tracing the circumnuclear accretion disk at larger orbital radii than are occupied by the 22 GHz masers. This newly identified population of AGN disk megamasers presents a motivation for developing VLBI capabilities at 183 GHz.

We study the evolution of various measures of quantumness of the curvature perturbation by integrating out the inaccessible entropic fluctuations in the multi-field models of inflation. In particular, we discuss the following measures of quantumness, namely purity, entanglement entropy and quantum discord. The models being considered in this work are ones that produce large scale curvature power spectra similar to those produced by single-field models of inflation. More specifically, we consider different multi-field models which generate nearly scale invariant and oscillatory curvature power spectrum and compare their quantum signatures in the perturbations with the corresponding single-field models. We find that, even though different models of inflation may produce the same observable power spectrum on large scales, they have distinct quantum signatures arising from the perturbation modes. This may allow for a way to distinguish between different models of inflation based on their quantum signatures. Intriguingly, this result generalizes to bouncing scenarios as well.

BO Cet is renowned cataclysmic variable having an orbital period above the period gap (0.139835 d) and showing both features of a Z Cam/IW And star and an SU UMa star. Using the Asteroid Terrestrial-impact Last Alert System (ATLAS) forced photometry and the All-Sky Automated Survey for Supernovae (ASAS-SN) Sky Patrol data, I found that BO Cet underwent a superoutburst in 2022 October-November after a series of short, normal outbursts with increasing amplitudes. This sequence of outbursts (supercycle) is what is seen in many SU UMa stars and this observation strengthened the suggestion that the accumulating mass and angular momentum in the disk during repeated normal outbursts caused a superoutburst even in the unusual system BO Cet. The outburst just preceding the superoutburst bore characteristics of an IW And-type standstill. This phenomenon reinforces the suggestion that the terminal outburst in IW And stars occurs when the disk radius reaches a certain limit. I consider that this outburst was a failed superoutburst, during which the disk reached the radius of the 3:1 resonance but the outburst faded before superhumps developed. In BO Cet with a mass ratio on the borderline of the stability of the 3:1 resonance, there may have been a competition between the effects of tidal truncation and the 3:1 resonance as the disk radius grew and the latter won in the current case. This finding in BO Cet might suggest that IW And-type and SU UMa-type phenomena are more strongly physically related than have been thought.

We derive and publish data-driven estimates of stellar metallicities [M/H] for ~120 million stars with low-resolution XP spectra published in Gaia DR3. The [M/H] values, along with Teff and logg, are derived using the XGBoost algorithm, trained on stellar parameters from APOGEE, augmented by a set of very metal-poor stars. XGBoost draws on a number of data features: the full set of XP spectral coefficients, narrowband fluxes derived from XP spectra, and broadband magnitudes. In particular, we include AllWISE magnitudes, as they reduce the degeneracy of Teff and dust reddening. We also include the parallax as a data feature, which helps constrain logg and [M/H]. The resulting mean stellar parameter precision is 0.1dex in [M/H], 50K in Teff, and 0.08dex in logg. This all-sky [M/H] sample is two orders of magnitude larger than published samples of comparable fidelity across -3<[M/H]<+0.5. Additionally, we provide a catalog of over 13 million bright (G<16) red giants whose [M/H] are vetted to be precise and pure. We present all-sky maps of the Milky Way in different [M/H] regimes that illustrate the purity of the dataset, and demonstrate the power of this unprecedented sample to reveal the Milky Way's structure from its heart to its disk.

Remote sensing hyperspectral and more generally spectral instruments are common tools to decipher surface features in Earth and Planetary science. While linear mixture is the most common approximation for compounds detection (mineral, water, ice, etc...), the transfer of light in surface and atmospheric medium are highly non-linear. The exact simulation of non-linearities can be estimated at very high numerical cost. Here I propose a very simple non-linear form (that includes the regular linear area mixture) of radiative transfer to approximate surface spectral feature. I demonstrate that this analytical form is able to approximate the grain size and intimate mixture dependence of surface features. In addition, the same analytical form can approximate the effect of Martian mineral aerosols. Unfortunately, Earth aerosols are more complex (water droplet, water ice, soot,...) and are not expected to follow the same trend.

Based on a new HI survey using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), combined with the Pan-STARRS1 images, we identified an isolated HI cloud without any optical counterpart, named FAST J0139+4328. The newly discovered HI cloud appears to be a typical disk galaxy since it has a double-peak shape in the global HI profile and an S-like rotation structure in the velocity-position diagram. Moreover, this disk galaxy has an extremely low absolute magnitude (M_B>-10.0 mag) and stellar mass (<6.9*10^5 Msun). Furthermore, we obtained that the HI mass of this galaxy is 8.3*10^7 Msun, and the dynamical mass to total baryonic mass ratio is 47+-27, implying that dark matter dominates over baryons in FAST J0139+4328. These findings provide observational evidence that FAST J0139+4328 is an isolated dark dwarf galaxy with a redshift of z=0.0083. This is the first time that an isolated dark galaxy has been detected in the nearby universe.

High-mass X-ray binaries (HMXBs) are a particular class of high-energy sources that require multi-wavelength observational efforts to be properly characterised. New identifications and the refinement of previous measurements are regularly published in the literature by independent teams of researchers and might, when they are collected in a catalogue, offer a tool for facilitating further studies of HMXBs. We update previous instances of HMXB catalogues in the Galaxy and provide the community easy access to the most complete set of observables on Galactic HMXBs. In addition to the fixed version that is available in Vizier, we also aim to host and maintain a dynamic version that can be updated upon request from users. Any modification will be logged in this version. Using previous HMXB catalogues supplemented by listings of hard X-ray sources detected in the past 20 years, we produced a base set of HMXBs and candidates by means of identifier and sky coordinate cross matches. We queried in Simbad for unreferenced HMXBs. We searched for as many hard X-ray, soft X-ray, optical, and infrared counterparts to the HMXBs as we could in well-known catalogues and compiled their coordinates. Each HMXB was subjected to a meticulous search in the literature to find relevant measurements and the original reference. We provide a catalogue of 152 HMXBs in the Galaxy with their best known coordinates, the spectral type of the companion star, systemic radial velocities, component masses, orbital period, eccentricity, and spin period when available. We also provide the coordinates and identifiers for each counterpart we found from hard X-rays to the near-infrared, including 111 counterparts from the recent Gaia DR3 catalogue.

Context. Small binary asteroid systems and pairs are thought to form through fission induced by spin up via the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. This process is expected to depend on their structural strength, hence composition. Aims. We aim to determine how taxonomic classes, used as a proxy for composition, distribute amongst binary asteroids and asteroid pairs compared to the general population. Methods. We compare the distribution of taxonomic classes of binary systems and pairs with that of a reference sample of asteroids. We build this sample by selecting asteroids to reproduce the orbital and size distribution of the binaries and pairs to minimize potential biases between samples. Results. A strong deficit of primitive compositions (C, B, P, D types) among binary asteroids and asteroid pairs is identified, as well as a strong excess of asteroids with mafic-silicate rich surface compositions (S, Q, V, A types). Conclusions. Amongst low mass, rapidly rotating asteroids, those with mafic-silicate rich compositions are more likely to form multiple asteroid systems than their primitive counterparts.

Water condensed as ice beyond the water snowline, the location in the Sun's natal gaseous disk where temperatures were below 170 K. As the disk evolved and cooled, the snowline moved inwards. A low temperature in the terrestrial planet-forming region is unlikely to be the origin of water on the planets, and the distinct isotopic compositions of planetary objects formed in the inner and outer disks suggest limited early mixing of inner and outer Solar System materials. Water in our terrestrial planets has rather been derived from H-bearing materials indigenous to the inner disk and delivered by water-rich planetesimals formed beyond the snowline and scattered inwards during the growth, migration, and dynamical evolution of the giant planets.

Aims. New or refined methods for determining instantaneous scalar and vector pattern speeds from a restricted domain are developed, for applications in N-body simulations or in galaxies, for which the stellar coordinates become increasingly better known. Methods. The general feature used throughout follows from the fact that the time-derivative of a function of the coordinates is linearly proportional to its rotation rate and its particle velocities. Knowing these allow in general to retrieve the instantaneous pattern speed vector by linear optimization. Similarly, if an invariant function depends both on the position and velocities, then its instantaneous rotation vectors in space can be retrieved. Knowing the accelerations allows to find rotation in velocity space. Results. The first 3 methods are based on the assumed rotational invariance of functions at each point in space or velocity space: 1) The 6D invariant function method measuring the pattern speed vectors in space and velocity space, 2) The differential/regional 3D Tremaine-Weinberg method evaluated over high signal to noise ratio regions, 3) The 3D Jacobi integral method yielding the potential pattern speed. Extensions to derive the rotation center position, speed and acceleration are introduced in 1) and 3). The last 2 methods are based on the assumed invariance of average functions of the particle coordinates. 4) The moment of inertia methods by using the derivative of the singular value decomposition, 5) The 2D Fourier method (3D for m = 2 mode) giving the mode rotation speeds. Pattern speed accelerations are also derived in 4) and 5). Conclusions. Depending on the available data, the different methods provide a choice of approaches.

We present and compare several methods to mitigate time-correlated (1/f) noise within the HI intensity mapping component of the MeerKAT Large Area Synoptic Survey (MeerKLASS). By simulating scan strategies, the HI signal, white and correlated noise, we assess the ability of various data processing pipelines to recover the power spectrum of HI brightness temperature fluctuations. We use MeerKAT pilot data to assess the level of 1/f noise expected for the MeerKLASS survey and use these measurements to create realistic levels of time-correlated noise for our simulations. We find the time-correlated noise component within the pilot data to be between 16 and 23 times higher than the white noise level at the scale of k = 0.04 Mpc^-1. Having determined that the MeerKAT 1/f noise is partially correlated across all the frequency channels, we employ Singular Value Decomposition (SVD) as a technique to remove 1/f noise but find that over-cleaning results in the removal of HI power at large (angular and radial) scales; a power loss of 20 per cent is seen for a 2-mode SVD clean at the scale of k = 0.04 Mpc^-1. We compare the impact of map-making using weighting by the full noise covariance (i.e. including a 1/f component), as opposed to just a simple unweighted binning, finding that including the time-correlated noise information reduces the excess power added by 1/f noise by up to 30 per cent.

This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We interpret the photometric data within a framework for modelling light curve morphology, and find that the models correctly predict the burst shapes, including their larger amplitudes and later peaks towards longer wavelengths. We are thus able to infer the start and end times of mass loading into the circumstellar disks of these stars. The disk sizes are typically 3-6 times the areas of the central star. The disk temperatures are ~40%, and the disk luminosities are ~10% of those of the central Be star, respectively. The available spectroscopy is consistent with inside-out evolution of the disk. Higher excitation lines have larger velocity widths in their double-horned shaped emission profiles. Our observations and analysis support the decretion disk model for outbursting Be stars.

DE Boo is a unique system, with an edge-on view through the debris disk around the star. The disk, which is analogous to the Kuiper belt in the Solar System, was reported to extend from 74 to 84 AU from the central star. The high photometric precision of the Characterising Exoplanet Satellite (CHEOPS) provided an exceptional opportunity to observe small variations in the light curve due to transiting material in the disk. This is a unique chance to investigate processes in the debris disk. Photometric observations of DE Boo of a total of four days were carried out with CHEOPS. Photometric variations due to spots on the stellar surface were subtracted from the light curves by applying a two-spot model and a fourth-order polynomial. The photometric observations were accompanied by spectroscopic measurements with the 1m RCC telescope at Piszk\'estet\H{o} and with the SOPHIE spectrograph in order to refine the astrophysical parameters of DE Boo. We present a detailed analysis of the photometric observation of DE Boo. We report the presence of nonperiodic transient features in the residual light curves with a transit duration of 0.3-0.8 days. We calculated the maximum distance of the material responsible for these variations to be 2.47 AU from the central star, much closer than most of the mass of the debris disk. Furthermore, we report the first observation of flaring events in this system. We interpreted the transient features as the result of scattering in an inner debris disk around DE Boo. The processes responsible for these variations were investigated in the context of interactions between planetesimals in the system.

Interstellar scattering (ISS) of radio pulsar emission can be used as a probe of the ionised interstellar medium (IISM) and causes corruptions in pulsar timing experiments. Two types of ISS phenomena (intensity scintillation and pulse broadening) are caused by electron density fluctuations on small scales (< 0.01 AU). Theory predicts that these are related, and both have been widely employed to study the properties of the IISM. Larger scales ($\sim$1-100\,AU) cause measurable changes in dispersion and these can be correlated with ISS observations to estimate the fluctuation spectrum over a very wide scale range. IISM measurements can often be modeled by a homogeneous power-law spatial spectrum of electron density with the Kolmogorov ($-11/3$) spectral exponent. Here we aim to test the validity of using the Kolmogorov exponent with PSR~J0826+2637. We do so using observations of intensity scintillation, pulse broadening and dispersion variations across a wide fractional bandwidth (20 -- 180\,MHz). We present that the frequency dependence of the intensity scintillation in the high frequency band matches the expectations of a Kolmogorov spectral exponent but the pulse broadening in the low frequency band does not change as rapidly as predicted with this assumption. We show that this behavior is due to an inhomogeneity in the scattering region, specifically that the scattering is dominated by a region of transverse size $\sim$40\,AU. The power spectrum of the electron density, however, maintains the Kolmogorov spectral exponent from spatial scales of 5$\times10^{-6}$\,AU to $\sim$100\,AU.

Recent studies have shown that 6.7 GHz methanol maser flares can be a powerful tool for verifying the mechanisms of maser production and even the specific signatures of accretion rate changes in the early stages of high-mass star formation. We characterize the spatial structure and evolution of methanol and water masers during a flare of methanol maser emission at 6.7 GHz in the HMYSO G24.33$+$0.14. VLBA was used to image the 6.7 and 12.2 GHz methanol and 22.2 GHz water vapor masers at three epochs guided by monitoring the methanol line with the Torun 32m telescope. The 6.7 GHz maser maps were also obtained with the EVN and LBA during the flare. WISE data were used to find correlations between the 6.7 GHz maser and IR fluxes. The 6.7 GHz methanol maser cloudlets are distributed over $\sim$3500 au, and the morphology of most of them is stable although their brightness varies following the course of the total flux density on a timescale of two months. The 12.2 GHz methanol maser cloudlets cover an area an order of magnitude smaller than that of 6.7 GHz emission, and both transitions emerge from the same masing gas. The 22.2 GHz maser cloudlets lie in the central region and show a systematic increase in brightness and moderate changes in size and orientation, together with the velocity drift of the strongest cloudlet during two months of the VLBI observing period. Time lag estimates imply the propagation of changes in the physical conditions of the masing region with a subluminal speed ($\sim$0.3c). A tight correlation of IR (4.6$\mu$m) and 6.7 GHz flux densities is found, supporting the radiative pumping model. Comparison with the 230 GHz ALMA data indicates that the methanol masers are distributed in the inner part of the rotating disk, whereas the 22.2 GHz emission traces the compact inner component of the bipolar outflow or a jet structure.

We analyse the evolution of the largest ionized region using the topological and morphological evolution of the redshifted 21-cm signal coming from the neutral hydrogen distribution during the different stages of reionization. For this analysis, we use the "Largest Cluster Statistics" - LCS. We mainly study the impact of the array synthesized beam on the LCS analysis of the 21-cm signal considering the upcoming low-frequency Square Kilometer Array (SKA1-Low) observations using a realistic simulation for such observation based on the 21cmE2E-pipeline using OSKAR. We find that bias in LCS estimation is introduced in synthetic observations due to the array beam. This in turn shifts the apparent percolation transition point towards the later stages of reionization. The biased estimates of LCS, occurring due to the effect of the lower resolution (lack of longer baselines) and the telescope synthesized beam will lead to a biased interpretation of the reionization history. This is important to note while interpreting any future 21-cm signal images from upcoming or future telescopes like the SKA, HERA, etc. We conclude that one may need denser $uv$-coverage at longer baselines for a better deconvolution of the array synthesized beam from the 21-cm images and a relatively unbiased estimate of LCS from such images.

The Event Horizon Telescope (EHT) has recently released high resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on accuracy and precision of numerical models of plasma and radiation around compact objects. The goal of this work is to speedup radiative transfer simulations used to create mock images of black holes for comparison with the EHT observations. A ray-tracing code for general relativistic and fully polarized radiative transfer through plasma in strong gravity is ported onto GPU. We describe our GPU implementation and carry out speedup tests using models of optically thin advection dominated accretion flow onto black hole realized semi-analytically and in 3-D GRMHD simulations, low and very high image pixel resolutions, and two different sets of CPU+GPUs. We show that a GPU unit with high double precision compute capability can significantly reduce the image production computational time, with speedup factor up to approximately $1200$. The significant speedup facilitates, e.g., dynamic model fitting to the EHT data, including polariemtric data. The method extension may enable studies of emission from plasma with nonthermal particle distribution functions for which accurate approximate synchrotron emissivities are not available. The significant speedup reduces the carbon footprint of the EHT image libraries generation by at least an order of magnitude.

We introduce a method to obtain the envelopes of eccentric orbits in axially symmetric potentials, $\Phi(R,z)$, endowed with $z$-symmetry of reflection. By making the transformation $z\rightarrow a+\sqrt{a^{2}+ z^{2}}$, with $a>0$, we compute the resulting mass density, referred here as the \emph{effective density} $\rho_{\rm ef}(R,z;a)$, in order to calculate the envelopes $Z(R)$ of orbits in the meridional plane $(R,z)$. We find that they obey the approximated formula $Z(R)\propto [\Sigma_{\rm ef}(R;a\approx 0)]^{-1/3}$, where $\Sigma_{\rm ef}(R;a)$ is the integrated surface density associated with $\rho_{\rm ef}(R,z;a)$. As examples we consider the dynamics in two potentials: the monopole plus quadrupole and the Kalnajs disc.

We applied the newly developed rose diagram overlay method to detect the layered structure of 88 nearby open clusters ($\leq$500~pc) on the three projections after the distance correction of their member stars, based on the catalog in literature. The results show that with the rose diagram overlay method, a total of 74 clusters in our sample have a layered structure, while the remaining clusters are without a clear layered structure. We for the first time defined the layered structure parameters for the sample clusters. Meanwhile, we found that the layered circle core area ($s$) has a strong positive correlation with the number of cluster members, while the kernel instability index ($\eta$) has a strong negative correlation with the number of cluster members. Our study provides a novel perspective for the detection of the layered structure of open clusters.

Galactic short-period close white dwarf binaries (CWDBs) are important objects for space-borne gravitational-wave (GW) detectors in the millihertz frequency bands. Due to the intrinsically low luminosity, only about 25 identified CWDBs are detectable by the Laser Interferometer Space Antenna (LISA), which are also known as verification binaries (VBs). The Gaia Early Data Release 3 (EDR3) provids a catalog containing a large number of CWDB candidates, which also includes parallax and photometry measurements. We crossmatch the Gaia EDR3 and Zwicky Transient Facility public data release 8, and apply period-finding algorithms to obtain a sample of periodic variables. The phase-folded light curves are inspected, and finally we obtain a binary sample containing 429 CWDB candidates. We further classify the samples into eclipsing binaries (including 58 HW Vir-type binaries, 65 EA-type binaries, 56 EB-type binaries, and 41 EW-type binaries) and ellipsoidal variations (209 ELL-type binaries). We discovered four ultrashort period binary candidates with unique light-curve shapes. We estimate the GW amplitude of all of our binary candidates, and calculate the corresponding signal-to-noise ratio (S/N) for TianQin and LISA. We find two (six) potential GW candidates with S/Ns greater than 5 in the nominal mission time of TianQin (LISA), which increases the total number of candidate VBs for TianQin (LISA) to 18 (31).

The nearest stars provide a fundamental constraint for our understanding of stellar physics and the Galaxy. The nearby sample serves as an anchor where all objects can be studied and understood with precise data. This work is an update of the 10 pc sample published by Reyl\'e et al. (2021) that used the unprecedented high precision parallax measurements from the early third data release of the astrometric space mission Gaia. We review this census, all updates being related to close binaries, brown dwarfs and exoplanets. We provide a new catalogue of 541 stars, brown dwarfs, and exoplanets in 336 systems within 10 pc from the Sun. This list is as volume-complete as possible from current knowledge and it provides a list of benchmark stars. We also explore the new products made available in the most recent third Gaia data release.

Supervised machine learning was recently introduced in high-contrast imaging (HCI) through the SODINN algorithm, a convolutional neural network designed for exoplanet detection in angular differential imaging (ADI) data sets. The benchmarking of HCI algorithms within the Exoplanet Imaging Data Challenge (EIDC) showed that (i) SODINN can produce a high number of false positives in the final detection maps, and (ii) algorithms processing images in a more local manner perform better. This work aims to improve the SODINN detection performance by introducing new local processing approaches and adapting its learning process accordingly. We propose NA-SODINN, a new deep learning architecture that better captures image noise correlations by training an independent SODINN model per noise regime over the processed frame. The identification of these noise regimes is based on a novel technique, named PCA-pmaps, which allows to estimate the distance from the star in the image from which background noise starts to dominate over residual speckle noise. NA-SODINN is also fed with local discriminators, such as S/N curves, which complement spatio-temporal feature maps when training the model.Our new approach is tested against its predecessor, as well as two SODINN-based hybrid models and a more standard annular-PCA approach, through local ROC analysis of ADI sequences from VLT/SPHERE and Keck/NIRC-2 instruments. Results show that NA-SODINN enhances SODINN in both the sensitivity and specificity, especially in the speckle-dominated noise regime. NA-SODINN is also benchmarked against the complete set of submitted detection algorithms in EIDC, in which we show that its final detection score matches or outperforms the most powerful detection algorithms, reaching a performance similar to that of the Regime Switching Model algorithm.

Contemporary research suggests that the reionisation of the intergalactic medium (IGM) in the early Universe was predominantly realised by star-forming (proto-)galaxies (SFGs). Due to observational constraints, our knowledge on the origins of sufficient amounts of ionising Lyman continuum (LyC) photons and the mechanisms facilitating their transport into the IGM remains sparse. Recent efforts have thus focussed on the study of local analogues to these high-redshift objects. We used archival spectroscopic SDSS DR12 data to select a sample of low-z He II 4686 emitters and restricted it to a set of SFGs with an emission line diagnostic sensitive to the presence of an AGN, which serves as our only selection criterion. Our final sample consists of eighteen low-mass, low-metallicity dwarf galaxies which appear to be predominantly ionised by stellar sources. We find large O32 ratios and [S II] deficiencies, which provide strong indications for these galaxies to be LyC Emitters (LCEs). At least 40% of these objects are candidates for featuring cosmologically significant LyC escape fractions >10%. Their SFHs exhibit strong similarities and almost all galaxies appear to contain an old (>1 Gyr) stellar component, while also harbouring a young, two-stage (~10 Myr and <1 Myr) starburst, which we speculate might be related to LyC escape. The properties of the compact emission line galaxies presented here align well with those observed in many local LCEs. In fact, our sample may prove as an extension to the rather small catalogue of local LCEs, as the extreme interstellar medium (ISM) conditions we find are assumed to facilitate LyC leakage. Notably, all of our eighteen candidates are significantly closer (z<0.1) than most established LCEs. If the inferred LyC photon loss is genuine, this demonstrates that selecting SFGs from He II 4686 is a powerful selection criterion in the search for LCEs.

The nuclear stellar disc (NSD) is a flat dense stellar structure at the heart of the Milky Way. Recent work shows that analogous structures are common in the nuclei of external spiral galaxies, where there is evidence of an age gradient that indicates that they form inside-out. However, the characterisation of the age of the NSD stellar population along the line of sight is still missing due to its extreme source crowding and the high interstellar extinction towards the Galactic centre. We aim to characterise the age of the stellar population at different average Galactocentric NSD radii to investigate for the first time the presence of an age gradient along the line of sight. We selected two groups of stars at different NSD radii via their different extinction and proper motion distribution. We analysed their stellar population by fitting their de-reddened $K_s$ luminosity functions with a linear combination of theoretical models. We find significant differences in the stellar population at different NSD radii, indicating the presence of an age gradient along the line of sight. Our sample from the closest edge of the NSD contains a significant fraction ($\sim40$ % of its total stellar mass) of intermediate-age stars (2-7 Gyr), that is not present in the sample from stars deeper inside the NSD, in which $\sim90 %$ of the stellar mass is older than 7 Gyr. Our results suggest that the NSD age distribution is similar as the one found in external galaxies and imply that bar-driven processes observed in external galaxies are similarly at play in the Milky Way.

We investigate the effect of observing cadence on the precision of radius ratio values obtained from transit light curves by performing uniform Markov Chain Monte Carlo fits of 46 exoplanets observed by the Transiting Exoplanet Survey Satellite (TESS) in multiple cadences. We find median improvements of almost 50% when comparing fits to 20s and 120s cadence light curves to 1800s cadence light curves, and of 37% when comparing 600s cadence to 1800s cadence. Such improvements in radius precision are important, for example, to precisely constrain the properties of the radius valley or to characterize exoplanet atmospheres. We also implement a numerical Information Analysis to predict the precision of parameter estimates for different observing cadences. We tested this analysis on our sample and found it reliably predicts the effect of shortening observing cadence with errors in the predicted % precision of <0,5% for most cases. We apply this method to 157 TESS object of interest that have only been observed with 1800s cadence to predict the precision improvement that could be obtained by reobservations with shorter cadences and provide the full table of expected improvements. We report the 10 planet candidates that would benefit the most from reobservations at short cadence. Our implementation of the Information Analysis for the prediction of the precision of exoplanet parameters, Prediction of Exoplanet Precisions using Information in Transit Analysis (PEPITA) is made publicly available.

We present a systematic analysis to determine and improve the pulsation periods of 1637 known long-period Mira variables in M33 using $gri$-band light curves spanning $\sim18$~years from several surveys, including M33 variability survey, Panoramic Survey Telescope and Rapid Response System, Palomar Transient Factory (PTF), intermediate PTF, and Zwicky Transient Facility. Based on these collections of light curves, we found that optical band light curves that are as complete as possible are crucial to determine the periods of distant Miras. We demonstrated that the machine learning techniques can be used to classify Miras into O-rich and C-rich based on the $(J-K_s)$ period--color plane. Finally, We derived the distance modulus to M33 using O-rich Miras at maximum light together with our improved periods as $24.67 \pm 0.06$~mag, which is in good agreement with the recommended value given in the literature.

The $\tau$-Herculids (IAU shower number #61 TAH) is a minor meteor shower associated with comet 73P/Schwassmann-Wachmann 3, a Jupiter-Family comet that disintegrated into several fragments in 1995. As a consequence of the nucleus break-up, possible increased meteor rates were predicted for 2022. On May 30-31, observation networks around the world reported two distinct peaks of TAH activity, around solar longitudes 69.02{\deg} and 69.42{\deg}. This work examines the encounter conditions of the Earth with meteoroids ejected from 73P during the splitting event and on previous perihelion passages. Numerical simulations suggest that the main peak observed in 2022 was caused by meteoroids ejected from the splitting nucleus with four times the typical cometary gas expansion speed. High-resolution measurements performed with the Canadian Automated Meteor Observatory indicate that these meteoroids are fragile, with estimated bulk densities of 250 kg/m$^3$. In contrast with the main peak, the first TAH activity peak in 2022 is best modelled with trails ejected prior to 1960. We find that ordinary cometary activity could have produced other TAH apparitions observed in the past, including in 1930 and 2017. The extension of our model to future years predicts significant returns of the shower in 2033 and 2049.

Pre-merger gravitational-wave (GW) sky-localisation of binary neutron star (BNS) and neutron star black hole (NSBH) coalescence events, would enable telescopes to capture precursors and electromagnetic (EM) emissions around the time of the merger. We propose a novel astrophysical scenario that could provide early-warning times of hours to days before coalescence with sub-arcsecond localisation, provided that these events are gravitationally lensed. The key idea is that if the BNS/NSBH is lensed, then so must the host galaxy identified via the EM counterpart. From the angular separation of the lensed host galaxy images, as well as its redshift and the (foreground) lens redshift, we demonstrate that we can predict the time delays assuming a standard lens model. Encouraged by the non-trivial upper limits on the detection rates of lensed BNS/NSBH mergers that we estimate for upcoming observing runs of the LIGO-Virgo-Kagra and third generation networks, we assess the feasibility and benefits of our method. To that end, we study the effect of limited angular resolution of the telescopes on our ability to predict the time delays. We find that with an angular resolution of $0.05''$, we can predict time delays of $> 1$ day with $1\sigma$ error-bar of $\mathcal{O}$(hours) at best. We also construct realistic time delay distributions of detectable lensed BNSs/NSBHs to forecast the early-warning times we might expect in the observing scenarios we consider.

We measure the three-dimensional cross-power spectrum of galaxy density and intrinsic alignment (IA) fields for the first time from the spectroscopic and imaging data of SDSS-III BOSS galaxies, for each of the four samples in the redshift range $0.2 < z < 0.75$. In the measurement we use the power spectrum estimator, developed in our previous work, to take into account the line-of-sight dependent projection of galaxy shapes onto the sky coordinate and the $E/B$-mode decomposition of the spin-2 shape field. Our method achieves a significant detection of the $E$-mode power spectrum with the total signal-to-noise ratio comparable with that of the quadrupole moment of the galaxy density power spectrum, while the measured $B$-mode power spectra are consistent with a null signal to within the statistical errors for all the galaxy samples. We also show that, compared to the previous results based on the two-dimensional projected correlation function, our method improves the precision of the linear shape bias parameter estimation by up to a factor of two thanks to the three-dimensional information. By performing a joint analysis of the galaxy density and IA power spectra in the linear regime, we constrain the isotropic and anisotropic local primordial non-Gaussianities (PNGs) parameters, $f_\mathrm{NL}^{s=0}$ and $f_\mathrm{NL}^{s=2}$, simultaneously, where the two types of PNGs induce characteristic scale-dependent biases at very large scales in the density and IA power spectra, respectively. We do not find any significant detection for both PNGs: the constraints $f^{s=0}_\mathrm{NL}=57^{+30}_{-29}$ and $f^{s=2}_\mathrm{NL} = -67_{-269}^{+285}$ ($68\%$ C.L.), respectively. Our method paves the way for using the IA power spectrum as a cosmological probe for current and future galaxy surveys.

IXPE (Imaging X-ray Polarimetry Explorer) is a Small Explorer mission by NASA and ASI, launched on December 9$^{th}$ 2021, dedicated to investigating X-ray polarimetry allowing angular-, time- and energy-resolved observations in the 2--8 keV energy band. IXPE is in Science Observation phase since January 2022; it comprises of three identical telescopes with grazing-incidence mirrors, each one having in the focal plane a Gas Pixel Detector (GPD). In this paper, we present a possible guideline to obtain an optimal background selection in polarimetric analysis, and a rejection strategy to remove instrumental background. This work is based on the analysis of IXPE observations, aiming to improve as much as possible the polarimetric sensitivity. In particular, the developed strategies have been applied ``as a case study'' to the IXPE observation of the 4U 0142+61 magnetar.

Imaging Science Subsystem (ISS) mounted on the Cassini spacecraft has taken a lot of images, which provides an important source of high-precision astrometry of some planets and satellites. However, some of these images are degraded by trailed stars. Previously, these degraded images cannot be used for astrometry. In this paper, a new method is proposed to detect and compute the centers of these trailed stars automatically. The method is then performed on the astrometry of ISS images with trailed stars. Finally, we provided 658 astrometric positions between 2004 and 2017 of several satellites that include Enceladus, Dione, Tethys, Mimas and Rhea. Compared with the JPL ephemeris SAT427, the mean residuals of these measurements are 0.11 km and 0.26 km in right ascension and declination, respectively. Their standard deviations are 1.08 km and 1.37 km, respectively. The results show that the proposed method performs astrometric measurements of Cassini ISS images with trailed stars effectively.

We revisit constraints on decaying very heavy dark matter (VHDM) using the latest ultrahigh-energy cosmic-ray (UHECR; $E\gtrsim 10^{18}$ eV) data and ultrahigh-energy (UHE) $\gamma$-ray flux upper limits, measured by the Pierre Auger Observatory. We present updated limits on the VHDM lifetime ($\tau_\chi$) for masses up to $\sim10^{15}$~GeV, considering decay into quarks, leptons, and massive bosons. In particular, we consider not only the UHECR spectrum but their composition data that favors heavier nuclei. Such a combined analysis improves the limits at $\lesssim 10^{12}$ GeV because VHDM decay does not produce UHECR nuclei. We also show that the constraints from the UHE $\gamma$-ray upper limits are $\sim10$ more stringent for all of the Standard Model final states we consider.

Active M-type stars are known to often produce superflares on the surface. Radiation from stellar (super-)flares is important for the exoplanet habitability, but the mechanisms are not well understood. In this paper, we report simultaneous optical spectroscopic and photometric observations of a stellar superflare on an active M dwarf YZ CMi with the 3.8-m Seimei telescope and the $Transiting\, Exoplanet\, Survey\, Satellite$. The flare bolometric energy was $1.3^{+1.6}_{-0.6} \times 10^{34} \,\rm{erg}$ and H$\alpha$ energy was $3.0^{+0.1}_{-0.1} \times 10^{32} \,\rm{erg}$. The H$\alpha$ emission line profile showed red asymmetry throughout the flare with a duration of $4.6-5.1 \,\rm{hrs}$. The velocity of the red asymmetry was $\sim 200-500 \,\rm{km\,s^{-1}}$ and line width of H$\alpha$ was broadened up to $34 \pm 14$ $\r{A}$. The redshifted velocity and line width of H$\alpha$ line decayed more rapidly than the equivalent width, and their time evolutions are correlated with that of the white-light emission. This indicates a possibility that the white light, H$\alpha$ red asymmetry, and H$\alpha$ line broadening originate from nearly the same site, i.e., the dense chromospheric condensation region heated by non-thermal electrons. On the other hand, the flux ratio of the redshifted excess components to the central components is enhanced one hour after the flare onset. This may be due to the change of the main source of the red asymmetry to the post-flare loops in the later phase of the flare.

Recovering sharper images from blurred observations, referred to as deconvolution, is an ill-posed problem where classical approaches often produce unsatisfactory results. In ground-based astronomy, combining multiple exposures to achieve images with higher signal-to-noise ratios is complicated by the variation of point-spread functions across exposures due to atmospheric effects. We develop an unsupervised multi-frame method for denoising, deblurring, and coadding images inspired by deep generative priors. We use a carefully chosen convolutional neural network architecture that combines information from multiple observations, regularizes the joint likelihood over these observations, and allows us to impose desired constraints, such as non-negativity of pixel values in the sharp, restored image. With an eye towards the Rubin Observatory, we analyze 4K by 4K Hyper Suprime-Cam exposures and obtain preliminary results which yield promising restored images and extracted source lists.

We investigate the possibility of reducing the number of degrees of freedom (d.o.f.) starting from generic metric theories of gravity by introducing multiple auxiliary constraints (ACs), under the restriction of retaining spatial covariance as a gauge symmetry. Arbitrary numbers of scalar-, vector- and tensor-type ACs are considered a priori, yet we find that no vector- and tensor-type constraints should be introduced, and that scalar-type ACs should be no more than four for the purpose of constructing minimally modified gravity (MMG) theories which propagate only two tensorial d.o.f., like general relativity (GR). Through a detailed Hamiltonian analysis, we exhaust all the possible classifications of ACs and find out the corresponding minimalizing and symmetrizing conditions for obtaining the MMG theories. In particular, no condition is required in the case of four ACs, hence in this case the theory can couple with matter consistently and naturally. To illustrate our formalism, we build a concrete model for this specific case by using the Cayley-Hamilton theorem and derive the dispersion relation of the gravitational waves, which is subject to constraints from the observations.

Bismuth germanium oxide ($\rm Bi_{4} Ge_{3} O_{12}$, BGO) scintillation crystals are widely used as detectors in the fields of particle physics and astrophysics due to their high density, and thus higher efficiency for gamma-ray detection. Owing to their good chemical stability, they can be used in any environment. For rare-event searches, such as dark matter and coherent elastic neutrino-nucleus scattering, BGO crystals are essential to comprehend the response of nuclear recoil. In this study, we have analyzed the events of neutron elastic scattering with oxygen in BGO crystals. Then, we have measured the quenching factor for oxygen recoil energy in the BGO crystal as a function of recoil energy by using a monoenergetic neutron source.

Between the mid-eighteenth and mid-nineteenth centuries, long-haul oceanic voyages of exploration and discovery routinely carried astronomical tent observatories to support land-based longitude determinations using heavy and cumbersome astronomical regulators and transit telescopes. Following James Cook's deployment of a pilot tent observatory on his first voyage to the Pacific in 1768-1771, the tent design was altered by William Bayly for more convenient use on Cook's second and third voyages to the Pacific. Bayly's design became the standard structure of tent observatories assigned to shipboard astronomers during the Age of Sail. By the middle of the nineteenth century, a subtle shift in focus had occurred, with tent observatories now being deployed to observe specific celestial events (such as the 1882 Venus transit or a variety of eclipses), while longitude determinations increasingly relied on the novel, compact and improved box chronometers of the day. A further shift in the application of tent observatories occurred towards the end of the nineteenth century, when astronomical applications largely gave way to a renewed focus on meteorological measurements.

Auroral kilometric radiation (AKR) is the paradigm of intense radio emission from planetary magnetospheres. Being close to the electron gyro frequency and/or its lower harmonics, its observation indicates the non-thermal state of the source plasma. Emission is produced when the plasma enters a state of energetic excitation which results in deformation of the electron distribution function. Under certain conditions this leads to "quasi-coherent" emission. It is believed that the weakly-relativistic electron-cyclotron-maser instability is responsible for this kind of radiation. Since energetically radio radiation normally is not of {primary} importance in the large-scale magnetospheric phenomena, AKR as such has, for the purposes of large-scale magnetospheric physics, become considered a marginal problem. Here this notion is questioned. AKR while applying to the auroral region mainly during magnetospherically disturbed times {carries just a fraction of the total substorm energy. It is, however, of diagnostic power in the physics of the upper auroral ionosphere and Space Weather research}. As a fundamental physical problem of generation of radiation in non-thermal plasmas {it remains not resolved yet. Many questions have been left open even when dealing only with the electron-cyclotron-maser. These can advantageously be studied in the magnetosphere proper both by observation and theory, the only continuously accessible place in space. The most important are listet here with hint on how they should be attacked. Its value is to be sought in the role it should play in application to the other magnetized planets, extra-solar planets, and to strongly magnetized astronomical objects} as an important {tool to diagnose the matter state responsible for radiation in the radio frequency range beyond thermal, shock or synchrotron radiation.

We report constraints on sub-GeV dark matter particles interacting with electrons from the first underground operation of DAMIC-M detectors. The search is performed with an integrated exposure of 85.23 g days, and exploits the sub-electron charge resolution and low level of dark current of DAMIC-M Charge-Coupled Devices (CCDs). Dark-matter-induced ionization signals above the detector dark current are searched for in CCD pixels with charge up to 7 e-. With this data set we place limits on dark matter particles of mass between 0.53 and 1000 MeV/c2, excluding unexplored regions of parameter space in the mass ranges [1.6,1000] MeV/c2 and [1.5,15.1] MeV/c2 for ultra-light and heavy mediator interactions, respectively.

In this work, we study the three-dimensional AdS gravitational vacuum stars (gravastars) in the context of gravity's rainbow theory. Then we extend it by adding the Maxwell electromagnetic field. We compute the physical features of gravastars, such as proper length, energy, entropy, and junction conditions. Our results show that the physical parameters for charged and uncharged states depend significantly on rainbow functions. Besides from charged state, they also depend on the electric field. Finally, we explore the stability of thin shell of three-dimensional (un)charged AdS gravastars in gravity's rainbow. We show that the structure of thin shell of these gravastars may be stable and is independent of the type of matter.

The next generation of spacecraft is anticipated to enable various new applications involving onboard processing, machine learning and decentralised operational scenarios. Even though many of these have been previously proposed and evaluated, the operational constraints of real mission scenarios are often either not considered or only rudimentary. Here, we present an open-source Python module called PASEOS that is capable of modelling operational scenarios involving one or multiple spacecraft. It considers several physical phenomena including thermal, power, bandwidth and communications constraints as well as the impact of radiation on spacecraft. PASEOS can be run both as a high-performance-oriented numerical simulation and/or in a real-time mode directly on edge hardware. We demonstrate these capabilities in three scenarios, one in real-time simulation on a Unibap iX-10 100 satellite processor, another in a simulation modelling an entire constellation performing tasks over several hours and one training a machine learning model in a decentralised setting. While we demonstrate tasks in Earth orbit, PASEOS is conceptually designed to allow deep space scenarios too. Our results show that PASEOS can model the described scenarios efficiently and thus provide insight into operational considerations. We show this in terms of runtime and overhead as well as by investigating the modelled temperature, battery status and communication windows of a constellation. By running PASEOS on an actual satellite processor, we showcase how PASEOS can be directly included in hardware demonstrators for future missions. Overall, we provide the first solution to holistically model the physical constraints spacecraft encounter in Earth orbit and beyond. The PASEOS module is available open-source online together with an extensive documentation to enable researchers to quickly incorporate it in their studies.

Motivated by a version of the ``Dipolar Dark Matter'' model, that aims at a relativistic completion of the phenomenology of MOND, we investigate the gravitational polarization mechanism in the canonical bimetric theory with an effective matter coupling. We explicitly show the fundamental obstacle why such theories cannot achieve a consistent gravitational polarization, and thus fail to recover the MONDian phenomenology at low energies.

Nuclear reactions involving alpha particles play an important role in various astrophysical processes such as the gamma-process of heavy element nucleosynthesis. The poorly known low-energy alpha-nucleus optical (AOMP) potential is a key parameter to estimate the rates of these reactions. The AOMP can be tested by measuring the cross section of alpha-scattering as well as alpha-induced reactions. Low energy elastic alpha-scattering on 144Sm has recently been measured with high precision. The aim of the present work was to complement that work by measuring the (a,n) cross sections on 144Sm at low energies. The experimental data shall be used to constrain the AOMP. From this potential the 144Sm(a,g)148Gd reaction rate can be derived with reduced uncertainties. The 144Sm(a,n)147Gd reaction was studied by bombarding Sm targets with alpha-beams provided by the cyclotron accelerator of Atomki. The cross section was determined using the activation method. The gamma-radiation following the decay of the 147Gd reaction product was measured with a HPGe detector. The experimental data are analyzed within the statistical model. The cross section was measured in the alpha-energy range between 13 and 20 MeV in 1 MeV steps, i.e., from close above the (a,n) threshold. The results were compared with statistical model calculations using various approaches and parametrizations for the AOMP, and excellent agreement was obtained for two recent potentials. However, these potentials cannot reproduce literature data for the 144Sm(a,g)148Gd reaction with the same accuracy. Constraints for the AOMP were derived from an analysis of the new 144Sm(a,n)147Gd data and literature data for 144Sm(a,g)148Gd. These constraints enable a determination of the reaction rate of the 144Sm(a,g)148Gd reaction with significantly reduced uncertainties of less than a factor of two.

This is a homage to the memory of Prof. Steven Weinberg who passed away on 23 July 2021.

Regular black holes are crucially important as approaches to solving the singularity problem, and in this paper, the accretion disk of Bardeen and Hayward models have been studied. For this purpose, we calculated the physical properties of black holes, including radiant energy, luminosity derivative, temperature, and conversion efficiency of accretion mass into radiation. The obtained results show that the non-zero-free parameters of regular black holes cause the radius of the innermost stable circular orbit of the disk to shift to smaller values. As a result of this displacement, we saw an increase in the profiles of radiant energy, luminosity derivative, and temperature. We also find that Bardeen and Hayward's black holes are more efficient in converting mass to radiation than Schwarzschild. Finally, we compared the free parameter of these two black holes with the spin of the rotating black hole and found that the Bardeen and Hayward black holes can mimic the slowly rotating Kerr black hole.

Warm inflation, its different particle physics model implementations and the implications of dissipative particle production for its cosmology are reviewed. First, we briefly present the background dynamics of warm inflation and contrast it with the cold inflation picture. An exposition of the space of parameters for different well-motivated potentials, which are ruled out, or severely constrained in the cold inflation scenario, but not necessarily in warm inflation, is provided. Next, the quantum field theory aspects in realizing explicit microscopic models for warm inflation are given. This includes the derivation of dissipation coefficients relevant in warm inflation for different particle field theory models. The dynamics of cosmological perturbations in warm inflation are then described. The general expression for the curvature scalar power spectrum is shown. We then discuss in details the relevant regimes of warm inflation, the weak and strong dissipative regimes. We also discuss the results predicted in these regimes of warm inflation and how they are confronted with the observational data. We explain how the dissipative dynamics in warm inflation can address several long-standing issues related to (post-) inflationary cosmology. This includes recent discussions concerning the so-called swampland criteria and how warm inflation can belong to the landscape of string theory.

We consider the positivity bounds for WIMP scalar dark matter with effective Higgs-portal couplings up to dimension-8 operators. Taking the superposed states for Standard Model Higgs and scalar dark matter, we show that the part of the parameter space for the effective couplings, otherwise unconstrained by phenomenological bounds, is ruled out by the positivity bounds on the dimension-8 derivative operators. We find that dark matter relic density, direct and indirect detection and LHC constraints are complementary to the positivity bounds in constraining the effective Higgs-portal couplings. In the effective theory obtained from massive graviton or radion, there appears a correlation between dimension-8 operators and other effective Higgs-portal couplings for which the strong constraint from direct detection can be evaded. Nailing down the parameter space mainly by relic density, direct detection and positivity bounds, we find that there are observable cosmic ray signals coming from the dark matter annihilations into a pair of Higgs bosons, $WW$ or $ZZ$.

The elliptical instability is an instability of elliptical streamlines, which can be excited by large-scale tidal flows in rotating fluid bodies, and excites inertial waves if the dimensionless tidal amplitude ($\epsilon$) is sufficiently large. It operates in convection zones but its interactions with turbulent convection have not been studied in this context. We perform an extensive suite of Cartesian hydrodynamical simulations in wide boxes to explore the interactions of the elliptical instability and Rayleigh-B\'enard convection. We find that geostrophic vortices generated by the elliptical instability dominate the flow, with energies far exceeding those of the inertial waves. Furthermore, we find that the elliptical instability can operate with convection, but it is suppressed for sufficiently strong convection, primarily by convectively-driven large-scale vortices. We examine the flow in Fourier space, allowing us to determine the energetically dominant frequencies and wavenumbers. We find that power primarily concentrates in geostrophic vortices, in wavenumbers that are convectively unstable, and along the inertial wave dispersion relation, even in non-elliptically deformed convective flows. Examining linear growth rates on a convective background, we find that convective large-scale vortices suppress the elliptical instability in the same way as the geostrophic vortices created by the elliptical instability itself. Finally, convective motions act as an effective viscosity on large-scale tidal flows, providing a sustained energy transfer (scaling as $\epsilon^2$). Furthermore, we find that the energy transfer resulting from bursts of elliptical instability, when it operates, is consistent with the $\epsilon^3$ scaling found in prior work.

We sample a wide range of relativistic mean-field theories (RMFTs) constrained by chiral effective field theory and properties of symmetric nuclear matter around saturation density and test them against known stellar structure constraints. This includes a relatively new mass and radius measurement of the central compact object in HESS J1731-347, which is reported to have an unusually low mass of $M=0.77^{+0.20}_{-0.17}\,\mathrm{M}_{\odot}$ and a compact radius of $R(1.4\mathrm{M}_{\odot})=10.4^{+0.86}_{-0.78}$ km. We show that none of the sampled nuclear RMFTs meet all stellar structure constraints at the $68\,\%$ credibility level, but that hybrid equations of state with a quark matter inner core and nuclear outer core fulfill all constraints at the $68\,\%$ credibility level. This indicates a tension between astrophysical constraints and low-energy nuclear theory if future measurements confirm the central values reported above.

Parameter inference, i.e. inferring the posterior distribution of the parameters of a statistical model given some data, is a central problem to many scientific disciplines. Posterior inference with generative models is an alternative to methods such as Markov Chain Monte Carlo, both for likelihood-based and simulation-based inference. However, assessing the accuracy of posteriors encoded in generative models is not straightforward. In this paper, we introduce `distance to random point' (DRP) coverage testing as a method to estimate coverage probabilities of generative posterior estimators. Our method differs from previously-existing coverage-based methods, which require posterior evaluations. We prove that our approach is necessary and sufficient to show that a posterior estimator is optimal. We demonstrate the method on a variety of synthetic examples, and show that DRP can be used to test the results of posterior inference analyses in high-dimensional spaces. We also show that our method can detect non-optimal inferences in cases where existing methods fail.