Papers for Wednesday, Mar 15 2023

The cm-wavelength radio flares on Cygnus X-3 have been studied for many years. Our recent paper (Spencer et al., 2022) looked again at the minor flares (flux density S of a few 100 mJy) and compared their properties with those of a sample of major flares (S > 1 Jy). We find that the minor flares have rise times and duration of ~ 1 hour, as opposed to ~ days for the major flares. Minor flares show more rapid expansion of the synchrotron radiation emitting material than in the strong flares. They also appear closer to the binary, whereas the large flares form a more developed jet, i.e. the jets formed in minor flares are short and wide, those in major flares are long and thin. We used the results of Fender & Bright (2019) to calculate the magnetic field and expansion velocity as a fraction B of the speed of light under minimum energy conditions when the source is optically thick for samples of minor and major flares. The minimum power in the source was found using the rise time of the flares. The minor flares have lower minimum power but have larger velocities and energy densities than the major flares. Minor flares can occur while a major flare is in progress, suggesting an indirect coupling between them. The spectral evolution of the minor flares can be explained by either an expanding synchrotron source or a shock model. Further investigation requires high resolution VLBI observations at the 1 mas level if we wish to understand the development of the source. The problem is that Cygnus X-3 is strongly scattered by the interstellar medium so high frequencies in the several 10s of GHz are required for the resolution needed. The minor flares are rapid [...] and hence only short snapshot VLBI observations can capture the structure. Large numbers of telescopes are required which is a problem at the highest frequencies. We discuss the VLBI possibilities and trade-offs for this awkward object.

The first meter-scale interstellar meteor (IM1) was detected by US government sensors in 2014, identified as an interstellar object candidate in 2019, and confirmed by the Department of Defense (DoD) in 2022. We use data from a nearby seismometer to localize the fireball to a $1 \mathrm{\; km^2}$ region within the $10^2 \mathrm{\; km^2}$ zone allowed by the precision of the DoD-provided coordinates. The improved localization is of great importance for a forthcoming expedition to retrieve the meteor fragments.

Decaying or annihilating dark matter and other exotic energy injections can modify the spectrum of the universe's photon bath, resulting in e.g. new contributions to spectral distortions of the cosmic microwave background blackbody spectrum and modifications to the temperature and ionization history of the universe. Here, we present an improved version of the $\texttt{DarkHistory}$ code, which is now capable of consistently calculating the spectrum of low-energy photons by properly treating the interactions of these photons with the levels of hydrogen atoms. Other changes to the code include a more detailed treatment of energy deposition by low-energy electrons, and spectral distortions from heating of the intergalactic medium. All of the improvements we have made to $\texttt{DarkHistory}$ are publicly available.

A black hole candidate orbiting a luminous star in the Large Magellanic Cloud young cluster NGC 1850 ($\sim100$Myr) has recently been reported based on radial velocity and light curve modelling. Subsequently, an alternative explanation has been suggested for the system: a bloated post-mass transfer secondary star (M$_{\rm initial} \sim 4-5M_{\odot}$, M$_{\rm current} \sim 1-2M_{\odot}$) with a more massive, yet luminous companion (the primary). Upon reanalysis of the MUSE spectra, we found that the radial velocity variations originally reported were underestimated ($K_{\rm 2,revised} = 176\pm3$km/s vs $K_{\rm 2,original} = 140\pm3$km/s) because of the weighting scheme adopted in the full-spectrum fitting analysis. The increased radial velocity semi-amplitude translates into a system mass function larger than previously deduced ($f_{\rm revised}$=2.83$M_{\odot}$ vs $f_{\rm original}$=1.42$M_{\odot}$). By exploiting the spectral disentangling technique, we place an upper limit of 10\% of a luminous primary source to the observed optical light in NGC1850 BH1, assuming that the primary and secondary are the only components contributing to the system. Furthermore, by analysing archival near-infrared data, we find clues to the presence of an accretion disk in the system. These constraints support a low-mass post-mass transfer star but do not provide a definitive answer whether the unseen component in NGC1850 BH1 is indeed a black hole. These results predict a scenario where, if a primary luminous source of mass M $\ge 4.7M_{\odot}$, is present in the system (given the inclination and secondary mass constraints), it must be hidden in a optically thick disk to be undetected in the MUSE spectra.

We calculate the post-recombination contribution to the Cosmic Microwave Background (CMB) spectral distortion due to general exotic energy injections, including dark matter (DM) decaying or annihilating to Standard Model particles. Upon subtracting the background distortion that would be present even without such energy injections, we find residual distortions that are still potentially large enough to be detectable by future experiments such as PIXIE. The distortions also have a high-energy spectral feature that is a unique signature of the injection of high-energy particles. We present a calculation of the global ionization history in the presence of decaying dark matter with sub-keV masses, and also show that previous calculations of the global ionization history in the presence of energy injection are not significantly modified by these additional spectral distortions. Our improved treatment of low-energy electrons allows us to extend calculations of the CMB anisotropy constraints for decaying DM down to arbitrarily low masses. We also recast these bounds as constraints on the coupling of axion-like particles to photons.

Intermediate-mass black holes (IMBHs) are the missing link between stellar-mass and supermassive black holes, widely believed to reside in at least some dense star clusters, but not yet observed directly. Tidal disruptions of white dwarfs (WDs) are luminous only for black holes less massive than $\sim 10^5\,M_{\odot}$, therefore providing a unique smoking gun that could finally probe the existence of IMBHs beyond any reasonable doubt. Here, we investigate the tidal captures of WDs by IMBHs in dense star clusters, and estimate a typical rate of $\sim 1\,{\rm Myr}^{-1}$ for galactic nuclei and $\sim 0.01\,{\rm Myr}^{-1}$ for globular clusters. Following the capture, the WD inspirals onto the IMBH producing gravitational waves detectable out to $\sim100$ Mpc by LISA for $\sim 10^4\,M_{\odot}$ IMBHs. The subsequent tidal stripping/disruption of the WD can also release bright X-ray and gamma-ray emission with luminosities of at least $\gtrsim10^{40}\,\rm{erg\,s^{-1}}$, detectable by \textit{Chandra}, \textit{Swift}, and upcoming telescopes, such as the \textit{Einstein Probe}.

We have measured the mass of metals in 13 submillimetre galaxies at z~4 in which the gas, based on previous observations, lies in a cold rotating disk. We measured the metal masses using either the submillimetre line or continuum emission from three tracers of the overall metal content - carbon atoms, carbon monoxide molecules and dust grains - using the first calibration of this technique that treats all three tracers simultaneously (Dunne et al. 2022). We obtain very similar mass estimates from the different tracers, which are similar to the entire metal content of a present-day massive early-type galaxy. We used the dynamical masses of these galaxies to set an upper limit on the mass of the interstellar medium in each galaxy, allowing us to set a lower limit on the metal abundance of the ISM, finding values for many of the galaxies well above the solar value. The only caveat to this result is that we cannot be certain for all galaxies that the emission from the tracer comes from within the region covered by the dynamical analysis. Our high values for the metal abundance are supported by a recent JWST study of one galaxy which found a super-solar metal abundance approximately 3 times greater than our lower limit. Using two chemical evolution models, we show that the high metal masses and metal abundances are what is expected shortly after the formation of a galaxy for a top-heavy IMF. We suggest a scenario for galaxy evolution in which massive galaxies reach a high metal abundance during their formation phase, which is then gradually reduced by dry mergers with lower mass galaxies. We use the chemical-evolution models to show that the metals in the outflows from massive early-type galaxies in their formation phase can quantitatively explain the long-standing puzzle that approximately 75% of the metals in clusters of galaxies is in the intracluster gas rather than in the galaxies.

The fraction of ionizing photons that escape their host galaxies to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era. In this paper we use the Meraxes semi-analytic galaxy formation model to infer the mean ionizing photon escape fraction and its dependence on galaxy properties through joint modelling of the observed high redshift galaxy population and existing constraints on the reionization history. Using a Bayesian framework, and under the assumption that escape fraction is primarily related to halo mass, we find that the joint constraints of the UV luminosity function, CMB optical depth, and the Ly$\alpha$ forest require an escape fraction of $(18\pm5)\%$ for galaxies within haloes of $M\lesssim10^{9}$M$_\odot$ and $(5\pm2)\%$ for more massive haloes. When considered with respect to galaxy properties, we find that the transition from values of escape fraction from $\sim18\%$ to $\sim5\%$ occurs at stellar masses of $M_\star\sim10^7$M$_\odot$, nearly independent of redshift. As a function of redshift, we find that reionization is dominated by the smaller $M_\star\lesssim10^7$M$_\odot$ galaxies with high escape fractions at $z\gtrsim6$, and by the larger $M_\star\gtrsim10^7$M$_\odot$ galaxies with lower escape fractions at $z\lesssim6$. Galaxies with star formation rates of $10^{-2.5}$M$_\odot$yr$^{-1}$ to $10^{-1.5}$M$_\odot$yr$^{-1}$ provide the dominant source of ionizing photons throughout reionization. Our results are consistent with recent direct measurements of a $\sim5\%$ escape fraction from massive galaxies at the end of reionization, and support recent findings using both hydrodynamic simulations and direct observations which indicate that the dominant sources of escaping ionizing photons during reionization are low mass galaxies.

We have performed low-resolution ground-based spectroscopy of HATS-46 b in transmission, using the EFOSC2 instrument on the ESO New Technology Telescope (NTT). HATS-46 b is a highly-inflated exoplanet that is a prime target for transmission spectroscopy, having a Jupiter-like radius (0.95 R$_\textrm{Jup}$) but a much lower mass (0.16 M$_\textrm{Jup}$). It orbits a G-type star with a 4.7 d period, giving an equilibrium temperature of 1100 K. We observed one transit of HATS-46 b with the NTT, with the time-series spectra covering a wavelength range of 3900 - 9000 Angstrom at a resolution of $R \sim 380$. We achieved a remarkably precise transmission spectrum of 1.03 $\times$ photon noise, with a median uncertainty of $357$ ppm for $\sim 200$ Angstrom wide bins, despite the relative faintness of the host star with $V_{\mathrm{mag}} = 13.6$. The transmission spectrum does not show strong absorption features and retrievals favour a cloudy model, ruling out a clear atmosphere with $3.0\sigma$ confidence. We also place a conservative upper limit on the sodium abundance under the alternative scenario of a clear atmosphere. This is the eighth planet in the LRG-BEASTS survey, which uses 4m-class telescopes such as the NTT to obtain low-resolution transmission spectra of hot Jupiters with precisions of around one atmospheric scale height.

We compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a 3-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientations (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfv\'enic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above $\sim$$4 \times 10^{3}$ cm$^{-3}$. This transition is recovered in most of our synthetic observations of optically thin molecular tracers, however for $^{12}$CO it does not occur and the PRS remains in parallel alignment across the whole observer-space. We calculate the physical depth of the optical depth $\tau = 1$ surface and find that for $^{12}$CO it is largely located in front of the cloud midplane, suggesting that $^{12}$CO is too optically thick and instead mainly probes low volume density gas. In our super-Alfv\'enic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend toward no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfv\'enic simulations when viewed with the magnetic field in the plane-of-the-sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field.

With tantalizing evidence of the recent e-Rosita mission, re-discovering very soft X-rays and EUV radiation from a cluster of galaxies or its environment, the question of the origin of cluster EUV excess is revisited in this work. It will be shown that the gas temperature, density, and frozen-in magnetic field of the intracluster medium, collectively support the emission and propagation of coherent \u{C}erenkov radiation, which is low frequency and large amplitude radiation capable of accelerating charged particles to relativistic speeds. Owing to the spectrum of \u{C}erenkov radiation, most of the incipient relativistic electrons undergo inverse-Compton scattering with the cosmic microwave background. It turns out the scattered radiation has observable ramifications only in the EUV band, of photon energy $70 -- 100$~eV, having a luminosity $\approx 10^{44}$~ergs~s$^{-1}$. This luminosity is on par with the EUV excess level detected from Abell 1795 and the Coma cluster. It should be stressed, as {\it caveat emptor}, that although the main subject is the putative large amplitude coherent \u{C}erenkov modes which are highly nonlinear, the results presented were derived using a quasi-linear approach to highlight the observable features of the phenomenon, namely the EUV emission.

We use JWST NIRCam observations from 1.5 to 4.44 micron of the massive lensing cluster field A2744 to show that extreme red selected galaxies (f(F444W) > 1 uJy and f(F444W)/f(F150W) > 3.5) pick out all 9 >4.5-sigma ALMA 1.1 or 1.2 mm sources and 17 of the 19 >5-sigma SCUBA-2 850 micron sources in the covered areas. Next, we use the red selected galaxies as priors to probe deeper in the SCUBA-2 850 micron image, identifying a sample of 44 >3-sigma SCUBA-2 850 micron sources with accurate positions, photometric redshifts, and magnifications. We show that there is a strong drop (a factor of around 5) in the ratio of the far-infrared luminosity to the bolometric luminosity at rest-frame 5000 A between z<4 and z>4, which we argue is due to the high-redshift sources having much lower extinction than the lower redshift sources.

With infrared flux contrasts larger than typically seen in hot Jupiter, tidally-locked white dwarf-brown dwarf binaries offer a superior opportunity to investigate atmospheric processes in irradiated atmospheres. NLTT5306 is such a system, with a M$_{BD}$ = 52 $\pm$ 3 M$_{Jup}$ brown dwarf orbiting a T$_{eff}$ = 7756 $\pm$ 35 K white dwarf with an ultra-short period of $\sim$102 min. We present HST/WFC3 spectroscopic phase curves of NLTT5306, consisting of 47 spectra from 1.1 to 1.7 microns with an average S/N=65 per wavelength. We extracted the phase-resolved spectra of the brown dwarf NLTT5306B, finding a small (<100 K) day/night temperature difference ($\sim$5 percent of the average day-side temperature). Our best-fit model phase curves revealed a complex wavelength-dependence on amplitudes and relative phase offsets, suggesting longitudinal-vertical atmospheric structure. The night-side spectrum was well-fit by a cloudy, non-irradiated atmospheric model while the day-side was best-matched by a cloudy, weakly irradiated model. Additionally, we created a simple radiative energy redistribution model of the atmosphere and found evidence for efficient day-to-night heat redistribution and a moderately high Bond albedo. We also discovered an internal heat flux much higher than expected given the published system age, leading to an age reassessment that resulted in NLTT5306B most likely being much younger. We find that NLTT5306B is the only known significantly irradiated brown dwarf where the global temperature structure is not dominated by external irradiation, but rather its own internal heat. Our study provides an essential insight into the drivers of global circulation and day-to-night heat transport as a function of irradiation, rotation rate, and internal heat.

Eccentric mergers are a signature of the dynamical formation channel of binary black holes (BBHs) in dense stellar environments and hierarchical triple systems. Here, we investigate the production of eccentric mergers via binary-single interactions, by means of $2.5\times10^{5}$ direct $\textit{N}$-body simulations. Our simulations include post-Newtonian terms up to the 2.5th order and are designed to reflect the environmental conditions of young (YSCs), globular (GCs), and nuclear star clusters (NSCs). We find that $0.6\%$ ($1\%$) of the mergers in NSCs (GCs) reach the coalescence with eccentricity ${>0.1}$, while in YSCs only $0.08\%$ of the mergers is eccentric. Approximately $\sim63\%$ of these events are produced by chaotic, resonant interactions where temporary binaries are continuously formed and destroyed, while $\sim31\%$ arise from prompt interactions where two black holes (BHs) merge almost radially. Lastly, $\sim 6\%$ of these eccentric mergers occur in temporary hierarchical triples. Binaries undergoing a flyby generally develop smaller tilt angles with respect to exchanges, with the distribution that peaks at $\sim15^{\circ}$ for the former, and $\sim90^{\circ}$ for the latter. This result challenges the idea that dynamics produces perfectly isotropic spin orientations. We further study the impact of the host star cluster on the population of BBH mergers. The environment dramatically affects BH retention: $0.8\%$, $12.5\%$, and $36.4\%$ of all the remnant BHs are retained in YSCs, GCs, and NSCs, respectively. The fraction of massive BHs also depends on the host cluster properties, with pair-instability ($60\leq\,$M$_{\rm BH}$/M$_{\odot}\leq$100) and intermediate-mass (M$_{\rm BH}\geq$100$\,$M$_{\odot}$) BHs accounting for approximately $\sim44\%$ and $1.6\%$ of the mergers in YSCs, $\sim33\%$ and $0.7\%$ in GCs, and $\sim28\%$ and $0.4\%$ in NSCs.

We present the list of variable stars we found in the Kepler superstamp data covering approximately 9 arcminutes from the central region of NGC 6791. We classified the variable stars based on the variability type and we established their cluster membership based on the available Gaia Early Data Release 3 astrometry, by means of the Bayesian Gaussian mixture models. In total we found 278 variable objects, among which 17 binaries, 45 pulsators, 62 rotational and five unclassified variables are cluster members. The remaining 28 binaries, 25 pulsators, 83 rotational, four unclassified and nine unidentified variables are either not members or their membership is not established. In the case of eclipsing binaries we calculated the mid-times of eclipses and derived ephemerides. We searched for eclipse timing variation by means of the observed minus calculated diagrams. Only three objects show significant orbital period variation. Independently of a report published just recently by Colman et al(2022) we found 119 new variables. We used isochrones calculated within the MIST project and derived the age (8.91 Gyr), average distance (4134 pc) and iron content [Fe/H] (0.26-0.28), of NGC 6791. Using the cluster members with membership probabilities greater than 0.9, we calculated the distance to the cluster of 4123(31) pc, which agrees with the result from our isochrone fitting.

Dark matter energy injection in the early universe modifies both the ionization history and the temperature of the intergalactic medium. In this work, we improve the CMB bounds on sub-keV dark matter and extend previous bounds from Lyman-$\alpha$ observations to the same mass range, resulting in new and competitive constraints on axion-like particles (ALPs) decaying into two photons. The limits depend on the underlying reionization history, here accounted self-consistently by our modified version of the publicly available {\tt DarkHistory} and {\tt CLASS} codes. Future measurements such as the ones from the CMB-S4 experiment may play a crucial, leading role in the search for this type of light dark matter candidates

We present results on the properties of extreme gas outflows in massive ($\rm M_* \sim$10$^{11} \ \rm M_{\odot}$), compact, starburst ($\rm SFR \sim$$200 \, \rm M_{\odot} \ yr^{-1}$) galaxies at z = $0.4-0.7$ with very high star formation surface densities ($\rm \Sigma_{SFR} \sim$$2000 \,\rm M_{\odot} \ yr^{-1} \ kpc^{-2}$). Using optical Keck/HIRES spectroscopy of 14 HizEA starburst galaxies we identify outflows with maximum velocities of $820 - 2860$ \kmps. High-resolution spectroscopy allows us to measure precise column densities and covering fractions as a function of outflow velocity and characterize the kinematics and structure of the cool gas outflow phase (T $\sim$10$^4$ K). We find substantial variation in the absorption profiles, which likely reflects the complex morphology of inhomogeneously-distributed, clumpy gas and the intricacy of the turbulent mixing layers between the cold and hot outflow phases. There is not a straightforward correlation between the bursts in the galaxies' star formation histories and their wind absorption line profiles, as might naively be expected for starburst-driven winds. The lack of strong \mgii \ absorption at the systemic velocity is likely an orientation effect, where the observations are down the axis of a blowout. We infer high mass outflow rates of $\rm \sim$50 $-$ 2200 $\rm M_{\odot} \, yr^{-1}$, assuming a fiducial outflow size of 5 kpc, and mass loading factors of $\eta\sim$5 for most of the sample. %with $\eta\sim$20 for two galaxies. While these values have high uncertainties, they suggest that starburst galaxies are capable of ejecting very large amounts of cool gas that will substantially impact their future evolution.

We build a bijective mapping between different physical fields from hydrodynamic CAMELS simulations. We train a CycleGAN on three different setups: translating dark matter to neutral hydrogen (Mcdm-HI), mapping between dark matter and magnetic fields magnitude (Mcdm-B), and finally predicting magnetic fields magnitude from neutral hydrogen (HI-B). We assess the performance of the models using various summary statistics, such as the probability distribution function (PDF) of the pixel values and 2D power spectrum ($P(k)$). Results suggest that in all setups, the model is capable of predicting the target field from the source field and vice versa, and the predicted maps exhibit statistical properties which are consistent with those of the target maps. This is indicated by the fact that the mean and standard deviation of the PDF of maps from the test set is in good agreement with those of the generated maps. The mean and variance of $P(k)$ of the real maps agree well with those of generated ones. The consistency tests on the model suggest that the source field can be recovered reasonably well by a forward mapping (source to target) followed by a backward mapping (target to source). This is demonstrated by the agreement between the statistical properties of the source images and those of the recovered ones.

In order to study the initial conditions of massive star formation, we have previously built a sample of 463 high-mass starless clumps (HMSCs) across the inner Galactic plane covered by multiple continuum surveys. Here, we use $^{13}$ CO(2-1) line data from the SEDIGISM survey, which covers 78$^{\circ}$ in longitude ($-60^{\circ}<l<18^{\circ}$, $\vert b\vert<0.5^{\circ}$) with 30$^{\prime \prime}$ resolution, to investigate the global dynamical state of these parsec-scale HMSCs (207 sources with good quality data, mass $10^{2}\sim 10^{5}\ \rm M_{\odot}$, size $0.1\sim3.6$ pc). We find that most HMSCs are highly turbulent with a median Mach number $\mathcal{M_{S}}\sim$ 8.2, and 44\%$\sim$55\% of them are gravitationally bound (with virial parameter $\alpha_{\rm vir} \lesssim 2$) if no magnetic fields were present. A median magnetic field strength of 0.33$\sim$0.37 mG would be needed to support these bound clumps against collapse, in agreement with previous observations of magnetic fields in massive star formation regions. Luminosity-to-mass ratio, an important tracer for evolutionary stage, is strongly correlated with dust temperature. Magnetic field strength is also corrected with density. The Larson linewidth-size scaling does not hold in HMSCs. This study advances our understanding of global properties of HMSCs, and our high-resolution ALMA observations are on the way to study the resolved properties.

We investigate the possibility of oscillon formation during the preheating phase of asymmetric inflationary potentials. We analytically establish the existence of oscillon-like solutions for the Klein-Gordon equation for a polynomial potential of the form $V(\phi)=\frac{1}{2}\phi^2+A\phi^3+B\phi^4$ using the small amplitude analysis, which naturally arises as a Taylor expansion of the $\alpha$-attractor E-model for $\phi\ll M_\text{pl}$ and $\alpha\sim\mathcal{O}(1)$. We perform a detailed numerical analysis to study the formation of nonlinear structures in the $\alpha$-attractor E-model using the publicly available lattice simulation code \textsf{$\mathcal{C}\text{osmo}\mathcal{L}\text{attice}$} for parameters in the range $10^{-5}\lesssim\alpha\lesssim 5\times 10^{-4}$. We find the backreaction of the field fluctuations onto the evolution of the homogeneous inflaton condensate to be significant for $\alpha\lesssim 2\times 10^{-4}$ for which we observe the formation of highly nonlinear structures with average equation of state $w\simeq 0$. These nonlinear structures maybe interpreted as \textit{oscillons}, providing evidence that they can form during the inflaton oscillations around an asymmetric potential and are found to be present for the entirety of the runtime of our simulations, comprising $\gtrsim 40\%$ of the total energy density.

Solar flares result from the sudden release of energy deposited by sub-photospheric motions into the magnetic field of the corona. The deposited energy accumulates secularly between events. One may interpret the observed event statistics as resulting from a state-dependent Poisson process, in which the instantaneous flare rate is a function of the stress in the system, and a flare becomes certain as the stress approaches a threshold set by the micro-physics of the flare trigger. If the system is driven fast, and if the threshold is static and uniform globally, a cross-correlation is predicted between the size of a flare and the forward waiting time to the next flare. This cross-correlation is broadly absent from the \emph{Geostationary Operational Environmental Satellite} (\emph{GOES}) soft X-ray flare database. One also predicts higher cross-correlations in active regions where the shapes of the waiting time and size distributions match. Again there is no evidence for such an association in the \emph{GOES} data. The data imply at least one of the following: i) the threshold at which a flare is triggered varies in time; ii) the rate at which energy is driven into active regions varies in time; iii) historical flare catalogs are incomplete; or iv) the description of solar flares as resulting from a build-up and release of energy, once a threshold is reached, is incomplete.

We present results of [CII]$\,158\,\rm{\mu m}$ emission line observations, and report the spectroscopic redshift confirmation of a strongly lensed ($\mu\sim20$) star-forming galaxy, MACS0308-zD1 at $z=6.2078\pm0.0002$. The [CII] emission line is detected with a signal-to-noise ratio $>6$ within the rest-frame UV bright clump of the lensed galaxy (zD1.1) and exhibits multiple velocity components; the narrow [CII] has a velocity full-width-half-maximum (FWHM) of $110\pm20\,\rm{km/s}$, while broader [CII] is seen with an FWHM of $230\pm20\,\rm{km/s}$. The broader [CII] component is blueshifted ($-80\pm20\,\rm{km/s}$) with respect to the narrow [CII] component, and has a morphology which extends beyond the UV-bright clump. We find that while the narrow [CII] emission is most likely associated with zD1.1, the broader component is possibly associated with outflowing gas. Based on the non-detection of $\lambda_{\rm 158\,\mu m}$ dust continuum, we find that MACS0308-zD1's star-formation activity occurs in a dust-free environment with the stringent upper limit of infrared luminosity $\lesssim9\times10^{8}\,{\rm L_{\odot}}$. Targeting this strongly lensed faint galaxy for follow-up ALMA and JWST observations will be crucial to characterize the details of typical galaxy growth in the early Universe.

The Square Kilometre Array (SKA) infrastructure will consist of two radio telescopes that will be the most sensitive telescopes on Earth. The SKA community will have to process and manage near exascale data, which will be a technical challenge for the coming years. In this respect, the SKA Global Network of Regional Centres plays a key role in data distribution and management. The SRCNet will provide distributed computing and data storage capacity, as well as other important services for the network. Within the SRCNet, several teams have been set up for the research, design and development of 5 prototypes. One of these prototypes is related to data management and distribution, where a data lake has been deployed using Rucio. In this paper we focus on the tasks performed by several of the teams to deploy new storage endpoints within the SKAO data lake. In particular, we will describe the steps and deployment instructions for the services required to provide the Rucio data lake with a new Rucio Storage Element based on StoRM and WebDAV within the Spanish SRC prototype.

We study the impact of kinks on the cosmic microwave background (CMB) anisotropies generated by cosmic string networks. To do so, we extend the Unconnected Segment Model to describe the stress-energy tensor of a network of cosmic strings with kinks and implement this extension in CMBACT to compute the CMB anisotropies generated by these wiggly string networks. Our results show that the inclusion of kinks leads, in general, to an enhancement of the temperature and polarization angular power spectra, when compared to those generated by cosmic string networks without small-scale structure with the same energy density, on scales corresponding to the distance between kinks. This enhancement, that is more prominent in the case of the temperature anisotropies, is essentially caused by a significant increase of the vector-mode anisotropies, since kinks, due to their shape, generate vortical motions of matter -- a phenomenon that is not taken into account when resorting to an effective description of wiggly cosmic strings.

The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar (${\Sigma}_{\star}$), HI (${\Sigma}_\mathrm{HI}$), and molecular gas (${\Sigma}_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average ${\Sigma}_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower ${\Sigma}_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical ${\Sigma}_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of ${\Sigma}_\mathrm{HI}$ or ${\Sigma}_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of ${\Sigma}_\mathrm{mol}$-to-${\Sigma}_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results.

SX Phoenicis (SXP) variables are short period pulsating stars that exhibit a period-luminosity (PL) relation. We derived the gri-band PL and extinction-free period-Wesenheit (PW) relations, as well as the period-color (PC) and reddening-free period-Q-index (PQ) relations for 47 SXP stars in located in 21 globular clusters using the optical light curves taken from Zwicky Transient Facility (ZTF). These empirically relations were derived for the first time in the gri filters except for the g-band PL relation. We used our gi band PL and PW relations to derive a distance modulus to Crater II dwarf spheroidal which hosts one SXP variable. Assuming that the fundamental and first-overtone pulsation mode for the SXP variable in Crater II, we found distance moduli of $20.03 \pm 0.23$ mag and $20.37 \pm 0.24$ mag, respectively, using the PW relation, where the latter is in excellent agreement with independent RR Lyrae based distance to Crater II dwarf galaxy.

Xtend is a soft X-ray imaging telescope developed for the X-Ray Imaging and Spectroscopy Mission (XRISM). XRISM is scheduled to be launched in the Japanese fiscal year 2022. Xtend consists of the Soft X-ray Imager (SXI), an X-ray CCD camera, and the X-ray Mirror Assembly (XMA), a thin-foil-nested conically approximated Wolter-I optics. The SXI uses the P-channel, back-illuminated type CCD with an imaging area size of 31 mm on a side. The four CCD chips are arranged in a 2$\times$2 grid and can be cooled down to $-120$ $^{\circ}$C with a single-stage Stirling cooler. The XMA nests thin aluminum foils coated with gold in a confocal way with an outer diameter of 45~cm. A pre-collimator is installed in front of the X-ray mirror for the reduction of the stray light. Combining the SXI and XMA with a focal length of 5.6m, a field of view of $38^{\prime}\times38^{\prime}$ over the energy range from 0.4 to 13 keV is realized. We have completed the fabrication of the flight model of both SXI and XMA. The performance verification has been successfully conducted in a series of sub-system level tests. We also carried out on-ground calibration measurements and the data analysis is ongoing.

In this multi-messenger astronomy era, all the observational probes are improving their sensitivities and overall performance. The Focusing on Relativistic universe and Cosmic Evolution (FORCE) mission, the product of a JAXA/NASA collaboration, will reach a 10 times higher sensitivity in the hard X-ray band ($E >$ 10~keV) in comparison with any previous hard X-ray missions, and provide simultaneous soft X-ray coverage. FORCE aims to be launched in the early 2030s, providing a perfect hard X-ray complement to the ESA flagship mission Athena. FORCE will be the most powerful X-ray probe for discovering obscured/hidden black holes and studying high energy particle acceleration in our Universe and will address how relativistic processes in the universe are realized and how these affect cosmic evolution. FORCE, which will operate over 1--79 keV, is equipped with two identical pairs of supermirrors and wideband X-ray imagers. The mirror and imager are connected by a high mechanical stiffness extensible optical bench with alignment monitor systems with a focal length of 12~m. A light-weight silicon mirror with multi-layer coating realizes a high angular resolution of $<15''$ in half-power diameter in the broad bandpass. The imager is a hybrid of a brand-new SOI-CMOS silicon-pixel detector and a CdTe detector responsible for the softer and harder energy bands, respectively. FORCE will play an essential role in the multi-messenger astronomy in the 2030s with its broadband X-ray sensitivity.

Bodies such as planets, moons, and asteroids in our solar system are the brightest objects in the low-frequency radio astronomy at $\lesssim$ 10 GHz. The low-frequency radio emissions from our solar system bodies exhibit various observed characteristics in the spectrum, polarization, periodicity, and flux. The observed characteristics are essential probes for explorations of the bodies' magnetosphere, atmosphere, surface, and even their interior. Generation and propagation theories of the radio emissions associate the characteristics with fundamental physics embedded in the environments: e.g., auroral electron acceleration, betatron acceleration, and atmospheric momentum transfer. Here we review previous studies on the low-frequency radio emissions from our solar system bodies to unveil some outstanding key questions on the dynamics and evolution of the bodies. To address the key questions by the future observations with the Square Kilometre Array (SKA), we made feasibility studies for detection and imaging of the radio emissions. Possible extensions of the solar system observations with SKA to the exoplanets are also proposed in the summary.

Fast transitions between different types of power density spectra (PDS) happening over timescales of several tens of seconds are rare phenomena in black hole X-ray binaries. In this paper, we report a broadband spectral-timing analysis of the fast transitions observed in the 2021 outburst of GX 339-4 using NICER and HXMT observations. We observe transitions between band-limited noise-dominated PDS and type-B quasi-periodic oscillations (QPOs), and their rapid appearance or disappearance. We also make a detailed comparison between the fast transitions in GX 339-4 with those seen in MAXI J1820+070 and MAXI J1348--630. By comparing the spectra of the periods with and without type-B QPOs, we find that the spectral ratios above 10 keV are nearly constant or slightly decreasing, and the values are different between sources. Below 10 keV, the flux change of the Comptonization component is inversely proportional to the flux change of the thermal component, suggesting that the appearance of type-B QPOs is associated with a redistribution of the accretion power between the disc and the Comptonizing emission region. The spectral ratios between the periods with type-B QPO and those with broadband noise are significantly different from that with type-B QPO and without type-B QPO, where the ratios (type-B QPO/broadband noise) show a maximum at around 4 keV and then decrease gradually towards high energies. Finally, we discuss the possible change of the geometry of the inner accretion flow and/or jet during the transitions.

Numerical simulations have shown the occurence of a scenario termed ''super-competitive accretion'', a term that describes a situation where only the central few objects grow supermassive while a larger number of stars compete for the reservoir, with significant accretion flows of $\gtrsim0.1$ M$_\odot$ yr$^{-1}$. This scenario particularly implies that the presence of fragmentation will not necessarily impeed the formation of a central massive object. We here explore this phenomenon using analytical estimates for growth via collisions and accretion, considering accretion due to self-gravity as well as Bondi-Hoyle accretion. Particularly, we explore under what conditions the accretion onto the central massive object breaks down, and derive a criterion that depends on the mass of the most massive object and the mass in fragments. For compact clusters with sizes about $0.1$~pc, we further find that the mass growth by collisions is comparable to the growth via accretion. Our results are validated through the comparison with numerical simulations, and we overall conclude that super-competitive accretion is a valid mechanism for the formation of very massive objects in the early Universe.

The dark age of the universe, when no luminous object had existed, ended with the birth of the first stars, galaxies, and blackholes. This epoch is called cosmic dawn. Cosmic reionization is the major transition of the intergalactic medium (IGM) in the universe driven by ionizing photons emitted from luminous objects. Although the epoch through the dark age to reionization is a milestone in the universe, our knowledge of this epoch has not been sufficient yet. Cosmic 21cm signal, which is emitted from neutral hydrogen, is expected to open a new window for this epoch. In this review paper, we first introduce the basic physics of the 21cm line and how first stars impact on the 21cm line signal. Next, we briefly summarize how we extract astrophysical information from the 21cm line signal by means of statistical and machine learning approaches. We also discuss the synergy between the 21cm line signal and other emission lines. Finally, we summarize the current status of 21cm experiments.

In the most distant reaches of the Universe, the 21-cm hyperfine transition in neutral hydrogen provides one of the only available tracers of large-scale structure. A number of instruments have been working and planned to measure the 21-cm line signals, and in particular, Experiment to Detect the Global EoR Signature (EDGES) recently has reported the first detection of an absorption signal, which corresponds to the 21-cm line global signal at the epoch of reionization (EoR). The future large radio telescope, Square Kilometre Array (SKA) will be able to deliver the high-precision measurement of 21-cm line emission/absorption signals. In this paper, we review the current status for the 21-cm line global and fluctuation signals from EoR to the dark ages, and then summarize the recent studies of how we probe the primordial Universe particularly motivated by the recent EDGES result and future observations by SKA. We focus on two applications for constraining cosmology with the EDGES result: constraints on the primordial magnetic fields and those on the primordial power spectrum. We also discuss the potential of future SKA for probing the inflationary Universe, by discussing expected constraints on the primordial power spectrum, its adiabaticity, and primordial non-Gaussianities from future observations of 21-cm fluctuations.

The solar system's Zodiacal Cloud is visible to the unaided eye, yet the origin of its constituent dust particles is not well understood, with a wide range of proposed divisions between sources in the asteroid belt and Jupiter Family comets. The amount of dust contributed by Oort Cloud comets is uncertain. Knowledge of the Zodiacal Cloud's structure and origins would help with NASA's aim of characterizing potentially Earth-like planets around nearby stars, since the exo-Earths must be studied against the light scattered from extrasolar analogs of our cloud. As the only example where the parent bodies can be tracked, our own cloud is critical for learning how planetary system architecture governs the interplanetary dust's distribution. Our cloud has been relatively little-studied in the near-ultraviolet, a wavelength range that is important for identifying potentially-habitable planets since it contains the broad Hartley absorption band of ozone. We show through radiative transfer modeling that our cloud's shape and size at near-UV wavelengths can be measured from Earth orbit by mapping the zodiacal light's flux and linear polarization across the sky. We quantify how well the cloud's geometric and optical properties can be retrieved from a set of simulated disk observations, using a Markov chain Monte Carlo analysis. The results demonstrate that observations with sufficient precision, covering a set of fields distributed along the ecliptic and up to the poles, can be used to determine the division between asteroidal, Jupiter Family, and Oort Cloud dust components, primarily via their differing orbital inclination distributions. We find that the observations must be repeated over a time span of several months in order to disentangle the zodiacal light from the Galactic background using the Milky Way's rotation across the sky.

We report a $z=2.30$ galaxy protocluster (COSTCO-I) in the COSMOS field, where the Lyman-$\alpha$ forest as seen in the CLAMATO IGM tomography survey does not show significant absorption. This departs from the transmission-density relationship (often dubbed the fluctuating Gunn-Peterson approximation; FGPA) usually expected to hold at this epoch, which would lead one to predict strong Ly$\alpha$ absorption at the overdensity. For comparison, we generate mock Lyman-$\alpha$ forest maps by applying FGPA to constrained simulations of the COSMOS density field, and create mocks that incorporate the effects of finite sightline sampling, pixel noise, and Wiener filtering. Averaged over $r=15\,h^{-1}\,\mathrm{Mpc}$ around the protocluster, the observed Lyman-$\alpha$ forest is consistently more transparent in the real data than in the mocks, indicating a rejection of the null hypothesis that the gas in COSTCO-I follows FGPA ($p=0.0026$, or $2.79 \sigma$ significance). It suggests that the large-scale gas associated with COSTCO-I is being heated above the expectations of FGPA, which might be due to either large-scale AGN jet feedback or early gravitational shock heating. COSTCO-I is the first known large-scale region of the IGM that is observed to be transitioning from the optically-thin photoionized regime at Cosmic Noon, to eventually coalesce into an intra-cluster medium (ICM) by $z=0$. Future observations of similar structures will shed light on the growth of the ICM and allow constraints on AGN feedback mechanisms.

We present the first polarized dust emission measurements of the Horsehead Nebula, obtained using the POL-2 polarimeter on the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera on the James Clerk Maxwell Telescope (JCMT). The Horsehead Nebula contains two sub-millimeter sources, a photodissociation region (PDR; SMM1) and a starless core (SMM2). We see well-ordered magnetic fields in both sources. We estimated plane-of-sky magnetic field strengths of 56$\pm$9 and 129$\pm$21 $\mu$G in SMM1 and SMM2, respectively, and obtained mass-to-flux ratios and Alfv\'en Mach numbers of less than 0.6, suggesting that the magnetic field can resist gravitational collapse and that magnetic pressure exceeds internal turbulent pressure in these sources. In SMM2, the kinetic and gravitational energies are comparable to one another, but less than the magnetic energy. We suggest a schematic view of the overall magnetic field structure in the Horsehead Nebula. Magnetic field lines in SMM1 appear have been compressed and reordered during the formation of the PDR, while the likely more-embedded SMM2 may have inherited its field from that of the pre-shock molecular cloud. The magnetic fields appear to currently play an important role in supporting both sources.

Compact objects (COs) can exist and evolve in an active galactic nuclei (AGN) disk, triggering a series of attractive CO-related multi-messenger events around a supermassive black hole. To better understand the nature of an embedded CO and its surroundings and to investigate CO-related events more accurately, in this paper, we study the specific accretion process of a CO in an AGN disk and explore the role of outflow feedback. We show that the asymptotically isotropic outflow generated from the CO hyper-Eddington accretion would truncate the circum-CO disk and push out its surrounding gas, resulting in recurrent formation and refilling of an outflow cavity to intermittently stop the accretion. Applying this universal cyclic process to black holes (BHs) and neutron stars (NSs), we find that, even if it is above the Eddington rate, the mass rate accreted onto a BH is dramatically reduced compared with the initial gas captured rate and thus consumes few mass of the AGN disk; outflow feedback on a NS is generally similar, but possesses complexities on the existence of a stellar magnetic field and hard surface. We demonstrate that although outflow feedback itself may be unobservable, it remarkably alters the CO evolution via reducing its mass growth rate, and the AGN disk can survive from the otherwise drastic CO accretion overlooking outflow. In addition, we discuss the potential influence of underdense cavity on CO-related events, which embodies the significant role of outflow feedback as well.

This paper proposes a new approach to separate the $\mu$ spectral distortions of the cosmic microwave background from foregrounds with poorly defined spectral shapes. The idea is based on finding the optimal response to the observed signal. This response is weakly sensitive to foregrounds with parameters that are within some certain limits of their possible variations and, at the same time, very sensitive to the amplitude of $\mu$ distortion. The algorithm described in this paper is stable, easy to implement, and simultaneously minimizes the response to foregrounds and photon noise.

We present photometric and spectroscopic observations of the Type Icn supernova (SN) 2021ckj. Spectral modeling of SN 2021ckj reveals that its composition is dominated by oxygen, carbon and iron group elements, and the photospheric velocity at peak is ~10000 km/s. From the light curve (LC) modeling applied to SNe 2021ckj, 2019hgp, and 2021csp, we find that the ejecta and CSM properties of Type Icn SNe are diverse. SNe 2021ckj and 2021csp likely have two ejecta components (an aspherical high-energy component and a spherical standard-energy component) with a roughly spherical CSM, while SN 2019hgp can be explained by a spherical ejecta-CSM interaction alone. The ejecta of SNe 2021ckj and 2021csp have larger energy per ejecta mass than the ejecta of SN 2019hgp. The density distribution of the CSM is similar in these three SNe, and is comparable to those of Type Ibn SNe. This may imply that the mass-loss mechanism is common between Type Icn (and also Type Ibn) SNe. The CSM masses of SN 2021ckj and SN 2021csp are higher than that of SN 2019hgp, although all these values are within the diversity seen in Type Ibn SNe. The early spectrum of SN 2021ckj shows narrow emission lines from C II and C III, without a clear absorption component, in contrast with that observed in SN 2021csp. The similarity of the emission components of these lines implies that the emitting regions of SNe 2021ckj and 2021csp have similar ionization states, and thus suggests that they have similar properties of the ejecta and CSM, which is inferred also from the LC modeling. Taking into account the difference in the strength of the absorption features, this heterogeneity may be attributed to viewing angle effects in otherwise common aspherical ejecta.

Faraday tomography is a new method of the study of cosmic magnetic fields enabled by broadband low-frequency radio observations. By Faraday tomography, it is possible to obtain the Faraday dispersion function which contains information on the line-of-sight distributions of magnetic fields, thermal electron density, and cosmic-ray electron density by measuring the polarization spectrum from a source of synchrotron radiation over a wide band. Furthermore, by combining it with 2-dimensional imaging, Faraday tomography allows us to explore the 3-dimensional structure of polarization sources. The application of Faraday tomography has been active in the last 20 years, when broadband observation has become technically feasible. However, the Faraday dispersion function is mathematically the Fourier transform of the polarization spectrum, and since the observable band is finite, it is impossible to obtain a complete Faraday dispersion function by performing Fourier transform. In addition, the Faraday dispersion function does not directly reflect the distribution of magnetic field, thermal-electron density, and cosmic-ray electron density in the physical space, and its physical interpretation is not straightforward. Despite these two difficult problems, Faraday tomography is attracting much attention because it has great potential as a new method for studying cosmic magnetic fields and magnetized plasmas. In particular, the next-generation radio telescope SKA (Square Kilometre Array) is capable of polarization observation with unprecedented sensitivity and broad bands, and the application of Faraday tomography is expected to make dramatic progress in the field of cosmic magnetic fields. In this review, we explain the basics of Faraday tomography with simple and instructive examples. Then representative algorithms to realize Faraday tomography are introduced and finally some applications are shown.

The Serpens Molecular Cloud is one of the most active sites of ongoing star formation at a distance of about 300 pc, and hence is very well-suited for studies of young low-mass stars and sub-stellar objects. In this paper, for the Serpens star forming region, we find potential members of the Young Stellar Objects population from the Gaia DR3 data and study their kinematics and distribution. We compile a catalog of 656 YSOs from available catalogs ranging from X-ray to the infrared. We use this as a reference set and cross-match it to find 87 Gaia DR3 member stars to produce a control sample with revised parameters. We queried the DR3 catalog with these parameters and found 1196 stars. We then applied three different density-based machine learning algorithms (DBSCAN, OPTICS and HDBSCAN) to this sample and found potential YSOs. The three clustering algorithms identified a common set of 822 YSO members from Gaia DR3 in this region. We also classified these objects using 2MASS and WISE data to study their distribution and the progress of star formation in Serpens.

It has been thought for decades that rotating black holes (BHs) power the energetic gamma-ray bursts (GRBs) and active galactic nuclei (AGNs), but the mechanism that extracts the BH energy has remained elusive. We here show that the solution to this problem arises when the BH is immersed in an external magnetic field and ionized low-density matter. For a magnetic field parallel to the BH spin, the induced electric field accelerates electrons outward and protons inward at spherical polar angles $-60^\circ\lesssim \theta \lesssim 60^\circ$ (hereafter polar region), and vice versa at $60^\circ\lesssim \theta \lesssim 120^\circ$ (hereafter equatorial region). For an antiparallel magnetic field, protons and electrons exchange their roles. The particles that are accelerated outward radiate off energy and angular momentum to infinity. The BH powers the process by reducing its energy and angular momentum by capturing polar protons and equatorial electrons with net negative energy and angular momentum. The electric potential allows for negative energy states outside the BH ergosphere, so the latter does not play any role in this electrodynamical BH energy extraction process.

To detail the development of RS Ophiuchi and the other Galactic Symbiotic-like Recurrent Novae throughout their outburst and quiescence, with a particular emphasis on the propagation of the shock wave during the outburst of the binaries. The spectral analysis has been performed using archival data according to the features of the individual datasets. Swift grism spectra were reduced and extracted using a combination of the pre-existing UVOTPY Python routine and newly written pipelines in Matlab. Other datasets were directly available in reduced form, already corrected for instrumental or background contamination, calibrated in wavelength and flux or intensity. The work on these was done through pipelines suited for reading the data and elaborating them to extract quantities of interest for the analysis. We find striking similarities in different outbursts of the same object and for different novae. For example, RS Oph 2021 was almost identical to the 2006 outburst, despite having occurred at a different orbital phase with the observations made from a different line of sight through the red giant wind. Despite the intrinsically different properties of the binaries, striking similarities are found for different systems of the same class, for instance, that the trend of the electron density over time during outburst appears to follow a general temporal development.

The formation of supermassive stars is believed to be an essential intermediate step for the formation of the massive black hole seeds that become the supermassive black holes powering the quasars observed in the early Universe. Numerical simulations have shown that supermassive stars can form in atomic-cooling halos when protostars reach accretion rates higher than $\sim10^{-2}$ M$_\odot$ yr$^{-1}$ and fragmentation is suppressed on pc scales. It is however still uncertain if a supermassive star still emerges when fragmentation occurs at smaller scales and a cluster of stars is formed instead. In this work we explore the problem of massive object formation due to the interplay of collisions and accretion in star clusters at low metallicity. We model a small embedded cluster of accreting protostars following sub-parsec scale fragmentation during the collapse of a primordial gas cloud and follow its evolution by performing $N$-body plus hydrodynamical simulations. Our results show that supermassive stars with 10$^3$ and 10$^4$ M$_\odot$ are always formed due to the interplay of collisions and accretion, and in some cases these objects are part of a binary system. The resulting supermassive star is surrounded by tens of smaller stars with typical masses in the range $1$-$100$ M$_\odot$.

We present the first detailed photometric analysis of ATO J108.6991+27.8306 (hereinafter as J108). The short-period close binary J108 was observed by the Nanshan 1 m Wide Field Telescope of the Xinjiang Astronomical Observatory. The obtained BVRI-band light curves were used to determine the photometric solution by using the 2003 version of the Wilson-Devinney code. J108 is a typical deep ( f > 50%), low mass ratio (q < 0.25) overcontact binary system with a mass ratio of q = 0.1501 and a fill-out factor of f = 50.1 %, suggesting that it is in the late evolutionary stage of contact binary systems. We found the target to be a W-type W UMa binary and provided evidence for the presence of starspots on both components. From the temperature-luminosity diagram, the main component is the evolved main sequence star with an evolutionary age of about 7.94 Gyr.

A problem of constructing the trajectory of a spacecraft flight to Venus within the framework of a mission including landing of a lander in a given region of the planet's surface is being considered. A new celestial mechanics related method based on the use of gravity assist maneuver near Venus is proposed to transfer the spacecraft to a heliocentric orbit resonant with the orbit of Venus, so that, at the next approach the planet, the given region of the surface becomes attainable for landing. It is shown that the best resonant orbit in terms of the cost of the characteristic velocity is an orbit with a 1:1 ratio of the period to the orbital period of Venus. A procedure for choosing one of possible resonant orbits depending on coordinates of the desired landing point on the surface and the launch date of the mission is described. An example of calculating the flight trajectory that ensures landing in the Vellamo-South region of the Venus surface at launch from the Earth in 2031 is considered.

Modern Giant Segmented Mirror Telescopes (GSMTs) like the Extremely Large Telescope, which is currently under construction, depend heavily on Adaptive Optics (AO) systems to correct for atmospheric distortions. However, a residual blur always remains in the astronomical images corrected by Single Conjugate AO (SCAO) systems due to fitting and bandwidth errors, which can mathematically be described by a convolution of the true image with a point spread function (PSF). Due to the nature of the turbulent atmosphere and its correction, the PSF is spatially varying, which is known as aniosplanatic effect. The PSF serves, e.g., as a quality measure for the science images and therefore needs to be known as accurately as possible. In this paper, we present an algorithm for PSF reconstruction in directions apart from the guide star direction in an SCAO system adapted to the needs of GSMTs focused on estimating the contribution of the anisoplanatic and generalized fitting error to the PSF. In particular, the PSF reconstruction algorithm for Single Conjugate Adaptive Optics from (Wagner, 2018) is combined with an algorithm for time-dependent atmospheric tomography from (Niebsch, 2021) to obtain a direction dependent reconstruction of the post-AO PSF. Results obtained in an end-to-end simulation tool show a qualitatively good reconstruction of the PSF compared to the PSF calculated directly from the simulated incoming wavefront as well as a stable performance with respect to imprecise knowledge of atmospheric parameters.

The determination of distance is fundamental in astrophysics. Gamma-ray sources are poorly characterized in this sense, as the limited angular resolution and poor photon-count statistics in gamma-ray astronomy makes it difficult to associate them to a multiwavelength object with known redshift. Taking the 1794 active galactic nuclei (AGNs) with known redshift from the Fermi-LAT latest AGN catalog, 4LAC-DR3, we employ machine learning techniques to predict the distance of the rest of AGNs based on their spectral and spatial properties. The state-of-the-art CatBoost algorithm reaches an average 0.56 R2 score with 0.46 root-mean-squared error (RMSE), predicting an average redshift value of $z_{avg}=0.63$, with a maximum $z_{max}=1.97$. We use the SHAP explainer package to gain insights into the variables influence on the outcome, and also study the extragalactic bakground light (EBL) implications. In a second part, we use this regression model to predict the redshift of the unassociated sample of the latest LAT point-source catalog, 4FGL-DR3, using the results of a previous paper to determine the possible AGNs within them.

We study the ionising photon production efficiency at the end of the Epoch of Reionisation ($z \sim 5.4 - 6.6$) for a sample of 35 bright Lyman-$\alpha$ emitters, this quantity is crucial to infer the ionising photon budget of the Universe. These objects were selected to have reliable spectroscopic redshifts, assigned based on the profile of their Lyman-$\alpha$ emission line, detected in the MUSE deep fields. We exploit medium-band observations from the JWST extragalactic medium band survey (JEMS) to find the flux excess corresponding to the redshifted \ha\ emission line. We estimate the UV luminosity by fitting the full JEMS photometry, along with several HST photometric points, with \texttt{Prospector}. We find a median ultra-violet continuum slope of $\beta = -2.21^{+0.26}_{-0.17}$ for the sample, indicating young stellar populations with little-to-no dust attenuation. Supported by this, we derive $\xi_{ion,0}$ with no dust attenuation and find a median value of log$\frac{\xi_{ion,0}}{\text{Hz erg}^{-1}} = 26.36^{+0.17}_{-0.14}$. If we perform dust attenuation corrections and assume a Calzetti attenuation law, our values are lowered by $\sim 0.1$ dex. Our results suggest Lyman-$\alpha$ emitters at the Epoch of Reionisation have enhanced $\xi_{ion,0}$ compared to previous estimations from literature, in particular, when compared to the non-Lyman-$\alpha$ emitting population. This initial study provides a promising outlook on the characterisation of ionising photon production in the early Universe. In the future, a more extensive study will be performed on the entire dataset provided by the JWST Advanced Deep Extragalactic Survey (JADES). Thus, for the first time, allowing us to

The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets.

Unusually high N/O abundances were recently reported for a very compact, intensively star-forming object GN-z11 at z=10.6 from JWST/NIRSpec observations. We present an empirical comparison with the C, N, and O abundance ratios in Galactic globular clusters (GCs) over a large metallicity range. We show that hot hydrogen-burning nucleosynthesis within supermassive stars (SMS) formed through runaway collisions can consistently explain the observed abundances ratio in GN-z11 and in GCs. This suggests that a proto-globular cluster hosting a SMS could be at the origin of the strong N-enrichment in GN-z11. Our model predicts the behavior of N/O, C/O, and Ne/O ratios as a function of metallicity, which can be tested if high-z objects similar to GN-z11 are detected with JWST in the future. Further studies and statistics will help differentiate the proto-GC scenario from the Wolf-Rayet scenario that we quantify with a population synthesis model, and shed more light on this peculiar object.

The acoustic glitches' signature present in solar-like stars holds invaluable information. Indeed, it is caused by a sharp variation in the sound speed, therefore carrying localised information. One such glitch is the helium glitch caused by the hydrogen and first and second partial helium ionisation region, allowing us to constrain the surface helium abundance. However, the function adjusted to the glitch signature depends non-linearly on the acoustic depth at which it occurs, He. Retrieving the faint glitch signature and estimating $\tau_{\textrm{He}}$ are difficult but crucial tasks to accurately measure the glitch parameters and, ultimately, accurately infer the helium abundance. In the present paper, we aim at providing a way to estimate $\tau_{\textrm{He}}$ using precise seismic indicators, independent of stellar modelling. Consequently, we aim at improving the WhoSGlAd (Whole Spectrum and Glitches Adjustment) method by automatically providing a model independent measure of the glitch's parameters. We compute the evolution of $T_{\textrm{He}}$, a dimensionless form of the acoustic depth, along a grid of models and adjust an empirical linear relation between $T_{\textrm{He}}$ and the mean large separation and frequency ratio as defined in WhoSGlAd. We further optimise over the value of this estimate to ensure the stability and accuracy of the approach. The proposed approach provides an excellent estimate of the acoustic depth and allows us to swiftly retrieve the glitch signature of observed spectra. We demonstrate that the we can accurately model the helium abundance of four Kepler targets by comparing model (both versions of WhoSGlAd) and literature values.

Primordial black hole (PBH) formation during cosmic phase transitions and annihilation periods, such as the QCD transition or the $e^+e^-$-annihilation, is known to be particularly efficient due to a softening of the equation of state. We present a detailed numerical study of PBH formation during the QCD epoch in order to derive an accurate PBH mass function. We also briefly consider PBH formation during the $e^+e^-$-annihilation epoch. Our investigation confirms that, for nearly scale-invariant spectra, PBH abundances on the QCD scale are enhanced by a factor $\sim 10^3$ compared to a purely radiation dominated Universe. For a power spectrum producing an (almost) scale-invariant PBH mass function outside of the transition, we find a peak mass of $M_{\rm pbh}\approx 1.9 M_{\odot}$ of which a fraction $f\approx 1.5\times 10^{-2}$ of the PBHs have a mass of $M_{\rm pbh} > 10 M_{\odot}$, possibly contributing to the LIGO-Virgo black hole merger detections. We point out that the physics of PBH formation during the $e^+e^-$-annihilation epoch is more complex as it is very close to the epoch of neutrino decoupling. We argue that neutrinos free-streaming out of overdense regions may actually hinder PBH formation.

Polycyclic aromatic hydrocarbons (PAHs) have been detected in numerous circumstellar discs. We propose the continuous processing of PAHs through clustering, adsorption on dust grains, and their reverse-processes as key mechanisms to reduce the emission-capable PAH abundance in protoplanetary discs. This cycle of processing is driven by vertical turbulence in the disc mixing PAHs between the disc midplane and the photosphere. We used a theoretical Monte Carlo model for photodesorption and a coagulation code in the disc midplane to estimate the relevance and timescale of these processes in a Herbig Ae/Be disc environment. By combining these components in a 1D vertical model, we calculated the gas-phase depletion of PAHs that stick as clusters on dust grains. Our results show that the clustering of gas-phase PAHs is very efficient, and that clusters with more than 100 monomers can grow for years before they are able to freeze out in the disc midplane. Once a PAH cluster is frozen on the dust grain surface, the large heat capacity of these clusters prevents them from evaporating off the grains in UV-rich environments such as the photosphere. Therefore, the clustering of PAHs followed by freeze-out can lead to a depletion of gas-phase PAHs in discs. Evaluated over the lifetime of protoplanetary discs, we find a depletion of PAHs by a factor that ranges between 50 and 1000 compared to the standard ISM abundance of PAHs in the inner disc through turbulent processing. Through these processes, we favour PAHs smaller than circumovalene as the major gas-phase emitters of the disc photosphere as larger PAH monomers cannot photodesorb from the grain surface. These gas-phase PAHs co-exist with large PAH clusters sticking on dust grains. We find a close relation between the amount of PAHs frozen out on dust grains and the dust population, as well as the strength of the vertical turbulence.

The first stars, galaxies, star clusters, and direct-collapse black holes are expected to have formed in low-mass ($\sim$$10^{5}-10^{9} ~ M_{\odot}$) haloes at Cosmic Dawn ($z \sim 10 - 30$) under conditions of efficient gas cooling, leading to gas collapse towards the centre of the halo. The halo mass cooling threshold has been analyzed by several authors using both analytical models and numerical simulations, with differing results. Since the halo number density is a sensitive function of the halo mass, an accurate model of the cooling threshold is needed for (semi-)analytical models of star formation at Cosmic Dawn. In this paper the cooling threshold mass is calculated (semi-)analytically, considering the effects of H$_2$-cooling and formation (in the gas phase and on dust grains), cooling by atomic metals, Lyman-$\alpha$ cooling, photodissociation of H$_2$ by Lyman-Werner photons (including self-shielding by H$_2$), photodetachment of H$^-$ by infrared photons, photoevaporation by ionization fronts, and the effect of baryon streaming velocities. We compare the calculations to several high-resolution cosmological simulations, showing excellent agreement. We find that in regions of typical baryon streaming velocities, star formation is possible in haloes of mass $\gtrsim 1-2 \times 10^6 ~ M_{\odot}$ for $z \gtrsim 20$. By $z \sim 8$, the expected Lyman-Werner background suppresses star formation in all minihaloes below the atomic-cooling threshold ($T_{\rm vir} = 10^4 ~ \textrm{K}$). The halo mass cooling threshold increases by another factor of $\sim$$4$ following reionization, although this effect is slightly delayed ($z \sim 4-5$) because of effective self-shielding.

Very metal-poor stars contain crucial information on the Milky Way's infancy. In our previous study \citep{Plotnikova_2022} we derived a mean age of $\sim$ 13.7 Gyr for a sample of these stars in the Sun's vicinity. In this work, we investigate the chemical and kinematics properties of these stars with the goal of obtaining some insights on their origin and their parent population. We did not find any Al-Mg anti-correlation, which lead us to the conclusion that these stars did not form in globular clusters, while the detailed analysis of their orbital parameters reveals that these stars are most probably associated with the pristine Bulge of the Milky Way. We then sketch a scenario for the formation of the Milky Way in which the first structure to form was the Bulge through rapid collapse. The other components have grown later on, with a significant contribution of accreted structures.

In this article we present RUBIS (Rotation code Using Barotropy conservation over Isopotential Surfaces), a fully Python-based centrifugal deformation program available at https://github.com/pierrehoudayer/RUBIS. The code has been designed to calculate the centrifugal deformation of a star or planet resulting from a given cylindrical rotation profile, starting from a spherically symmetric non-rotating model. Furthermore, it can handle models with discontinuities in the density profile. The underlying assumption in RUBIS is that the relationship between density and pressure is preserved during the deformation process. This leads to many procedural simplifications. For instance, RUBIS only needs to solve Poisson' equation, either in spheroidal or spherical coordinates depending on whether the 1D model has discontinuities or not. In this paper, we present the benefits of using RUBIS to deform polytropic models and more complex barotropic structures, thus providing, to a certain extent, insights into baroclinic models. The resulting structures can be used for a wide range of applications, including the seismic study of models. Finally, we illustrate how RUBIS is beneficial specifically in the analysis of Jupiter's gravitational moments, thanks to its ability to handle discontinuous models while retaining a high accuracy compared to current methods.

Supermassive stars (SMS) with masses $M \gtrsim 10^3-10^4 M_{\odot}$, formed by runaway collisions in young, massive, and dense star clusters have been invoked as a possible solution to the puzzles raised by the presence of multiple stellar populations and peculiar abundance patterns observed in globular clusters. However, such objects have not been observed so far. We developed observational strategies to search for SMS hosted within young massive clusters (thought to be the precursors of globular clusters, GCs), which could be applicable in a relatively general fashion, using both photometric and spectroscopic observations. We used theoretical predictions of spectra of SMS and SMS-hosting clusters, together with predictions from standard simple stellar populations to examine their impact on color-color diagrams and on individual optical spectral lines (primarily Hydrogen emission and absorption lines). As a first step, we apply our search strategies to a sample of $\sim 3000$ young star clusters (YSC) from two nearby galaxies with multi-band observations from the HST and optical integral-field spectroscopy obtained with MUSE on the Very Large Telescope. We focus on models for SMS with large radii (corresponding to $ Teff \lesssim 7000$ K), which predict strong Balmer breaks, and construct proper color-color diagrams to select the corresponding SMS-hosting cluster candidates. We show that their spectrophotometric properties are similar to that of normal clusters with ages of a few hundred Myr, which would, however, show signs of composite stellar populations, in particular the presence of nebular lines (H${\alpha}$ and others). Examining the photometry, overall SEDs, and the spectra of $\sim 100$ clusters with strong Balmer breaks, we have found several objects with peculiar SEDs, the presence of emission lines, or other peculiar signatures. [abridged]

We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG-100) cosmological simulations by reconstructing the cosmic web in each snapshot using the DisPerSE framework. We measure the distance from each galaxy with stellar mass log(M*/Msun)>=8 to the nearest node (dnode) and the nearest filament spine (dfil) and study the dependence of both median specific star formation rate (<sSFR>) and median gas fraction (<fgas>) on these distances. We find that <sSFR> of galaxies is only dependent on cosmic web environment at z<2, with the dependence increasing with time. At z<=0.5, 8<=log(M*/Msun)<9 galaxies are quenched at dnode<1 Mpc, and significantly star formation-suppressed at dfil<1 Mpc, trends which are driven mostly by satellite galaxies. At z<=1, in contrast to the monotonic rise in <sSFR> of log(M*/Msun)<10 galaxies with dnode and dfil, log(M*/Msun)>=10 galaxies actually experience an upturn in <sSFR> at dnode<0.2 Mpc (this is caused by both satellites and centrals). Much of this cosmic web-dependence of star formation activity can be explained by the evolution in <fgas>. Our results suggest that in the past ~10 Gyr, low-mass satellites are quenched by rapid gas stripping in dense environments near nodes and gradual gas starvation in intermediate-density environments near filaments, while at earlier times cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as SDSS and DESI, as well as those planned with JWST and Roman are required to test our predictions against observations.

Deep learning has become a popular trend in recent years in the machine learning community and has even occasionally become synonymous with machine learning itself thanks to its efficiency, malleability, and ability to operate free of human intervention. However, a series of hyperparameters passed to a conventional neural network (CoNN) may be rather arbitrary, especially if there is no surefire way to decide how to program hyperparameters for a given dataset. The random hivemind (RH) alleviates this concern by having multiple neural network estimators make decisions based on random permutations of features. The learning rate and the number of epochs may be boosted or attenuated depending on how all features of a given estimator determine the class that the numerical feature data belong to, but all other hyperparameters remain the same across estimators. This allows one to quickly see whether consistent decisions on a given dataset can be made by multiple neural networks with the same hyperparameters, with random subsets of data chosen to force variation in how data are predicted by each, placing the quality of the data and hyperparameters into focus. The effectiveness of RH is demonstrated through experimentation in the predictions of dangerous solar energetic particle events (SEPs) by comparing it to that of using both CoNN and the traditional approach used by ensemble deep learning in this application. Our results demonstrate that RH outperforms the CoNN and a committee-based approach, and demonstrates promising results with respect to the ``all-clear'' prediction of SEPs.

Extracting the CMB blackbody temperature power spectrum -- which is dominated by the primary CMB signal and the kinematic Sunyaev-Zel'dovich (kSZ) effect -- from mm-wave sky maps requires cleaning other sky components. In this work, we develop new methods to use large-scale structure (LSS) tracers to remove cosmic infrared background (CIB) and thermal Sunyaev-Zel'dovich (tSZ) contamination in such measurements. Our methods rely on the fact that LSS tracers are correlated with the CIB and tSZ signals, but their two-point correlations with the CMB and kSZ signals vanish on small scales, thus leaving the CMB blackbody power spectrum unbiased after cleaning. We develop methods analogous to delensing ($\textit{de-CIB}$ or $\textit{de-(CIB+tSZ)}$) to clean CIB and tSZ contaminants using these tracers. We compare these methods to internal linear combination (ILC) methods, including novel approaches that incorporate the tracer maps in the ILC procedure itself, without requiring exact assumptions about the CIB SED. As a concrete example, we use the $\textit{unWISE}$ galaxy samples as tracers. We provide calculations for a combined Simons Observatory and $\textit{Planck}$-like experiment, with our simulated sky model comprising eight frequencies from 93 to 353 GHz. Using $\textit{unWISE}$ tracers, improvements with our methods over current approaches are already non-negligible: we find improvements up to 20% in the kSZ power spectrum signal-to-noise ratio (SNR) when applying the de-CIB method to a tSZ-deprojected ILC map. These gains could be more significant when using additional LSS tracers from current surveys, and will become even larger with future LSS surveys, with improvements in the kSZ power spectrum SNR up to 50%. For the total CMB blackbody power spectrum, these improvements stand at 4% and 7%, respectively. Our code is publicly available at https://github.com/olakusiak/deCIBing.

Recently, the inner main belt asteroid (152830) Dinkinesh was identified as an additional fly-by target for the Lucy mission. The heliocentric orbit and approximate absolute magnitude of Dinkinesh are known, but little additional information was available prior to its selection as a target. In particular, the lack of color spectrophotometry or spectra made it impossible to assign a spectral type to Dinkinesh from which its albedo could be estimated. We set out to remedy this knowledge gap by obtaining visible wavelength spectra with the Keck telescope on 2022 November 23 and with Gemini-South on 2022 December 27. The spectra measured with the Keck I/Low Resolution Imaging Spectrometer (LRIS) and the Gemini South/Gemini Multi-Object Spectrograph South (GMOS-S) are most similar to the average spectrum of S- and Sq-type asteroids. The most diagnostic feature is the $\approx$15$\pm$1$\%$ silicate absorption feature at $\approx$0.9-1.0~micron. Small S- and Sq-type asteroids have moderately high albedos ranging from 0.17-0.35. Using this albedo range for Dinkinesh in combination with measured absolute magnitude, it is possible to derive an effective diameter and surface brightness for this body. The albedo, size and surface brightness are important inputs required for planning a successful encounter by the Lucy spacecraft.

The Sun could produce Axion-Like Particles (ALPs) with masses in the keV-range. A fraction of them would be trapped in the solar gravitational field and accumulate over cosmic times. Their decay into photons would contribute to the observed solar X-ray flux, whose measurements can be used to constrain ALP models. In this paper, we extend previous works through the addition of two mechanisms. First, we include the ALP production in the solar core via the photon coalescence, the dominant production mechanism for trapped ALPs, improving by one order of magnitude the existing limits between $3~\rm{keV}$ and $40~\rm{keV}$ in the parameter space $(g_{a\gamma\gamma}, m)$. Second, we determine for the first time, both analytically and by simulation, the probability for an ALP to be Compton-absorbed while crossing the Sun during its orbit. For $g_{ae}\neq 0$, we demonstrate that Compton absorption counter-balances the non-hadronic production, partially or entirely, resulting in two well-defined regimes in the exclusion limits, with a transition triggered by the coupling to electrons. Out of the transitional region, the solar X-ray constraints on ALPs are exclusively governed by the coupling to photons.

We show that extremal Kerr black holes are sensitive probes of new physics. Stringy or quantum corrections to general relativity are expected to generate higher-curvature terms in the gravitational action. We show that in the presence of these terms, asymptotically flat extremal rotating black holes have curvature singularities on their horizon. Furthermore, near-extremal black holes can have arbitrarily large tidal forces for infalling observers. In addition, we consider five-dimensional extremal charged black holes and show that higher-curvature terms can have a large effect on the horizon geometry.

We investigate the production of a spectator scalar dark matter field that is directly coupled to the inflaton during inflation and reheating. We consider two specific inflationary potentials, namely the Starobinsky and T-model of inflation, which satisfy the constraints on the scalar tilt, $n_s$, and tensor-to-scalar ratio, $r$, measured by the Planck satellite. Excitation of light scalar dark matter during inflation may result in large isocurvature perturbations, which can be avoided by inducing a sizable effective dark matter mass during the inflationary phase. For purely gravitational production, the Planck isocurvature constraints require the dark matter mass to be larger than the Hubble scale at horizon exit, with $m_{\chi} \gtrsim 0.5H_*$. For small bare dark matter masses $m_{\chi} \ll H_*$, these constraints translate into a lower bound on the dark matter coupling to the inflaton. We argue that these constraints can applied to a wide class of single-field slow-roll inflation models. We also derive isocurvature, dark matter abundance, and Lyman-$\alpha$ constraints on the direct coupling and bare dark matter mass.

We construct the first four-dimensional multi-black hole solution of general relativity with a positive cosmological constant. The solution consists of two static black holes whose gravitational attraction is balanced by the cosmic expansion. These static binaries provide the first four-dimensional example of non-uniqueness in general relativity without matter.

The antiproton flux measurements from AMS-02 offer valuable information about the nature of dark matter, but their interpretation is complicated by large uncertainties in the modeling of cosmic ray propagation. In this work we present a novel framework to efficiently marginalise over propagation uncertainties in order to obtain robust AMS-02 likelihoods for arbitrary dark matter models. The three central ingredients of this framework are: the neural emulator DarkRayNet, which provides highly flexible predictions of the antiproton flux; the likelihood calculator pbarlike, which performs the marginalisation, taking into account the effects of solar modulation and correlations in AMS-02 data; and the global fitting framework GAMBIT, which allows for the combination of the resulting likelihood with a wide range of dark matter observables. We illustrate our approach by providing updated constraints on the annihilation cross section of WIMP dark matter into bottom quarks and by performing a state-of-the-art global fit of the scalar singlet dark matter model, including also recent results from direct detection and the LHC.

If the temperature of the hot thermal plasma in the Early Universe was within a few orders of magnitude of the Planck scale $M_{\rm Pl}$, then the hoop conjecture predicts the formation of microscopic black holes from particle collisions in the plasma. Although these evaporated instantly, they would have left behind a relic abundance of all stable degrees of freedom which couple to gravity. Here we show that, upon minimal assumptions of a high reheat temperature and semiclassical black hole dynamics, this process could have produced the relic abundance of dark matter observed today for a particle mass anywhere in the range of $100~\mathrm{keV} \lesssim m_{dm} < M_{\rm Pl}$. The production mechanism does not rely on any additional assumptions about non-gravitational dark matter-Standard Model interaction.

Optical satellite-to-ground communication (OSGC) has the potential to improve access to fast and affordable Internet in remote regions. Atmospheric turbulence, however, distorts the optical beam, eroding the data rate potential when coupling into single-mode fibers. Traditional adaptive optics (AO) systems use a wavefront sensor to improve fiber coupling. This leads to higher system size, cost and complexity, consumes a fraction of the incident beam and introduces latency, making OSGC for internet service impractical. We propose the use of reinforcement learning (RL) to reduce the latency, size and cost of the system by up to $30-40\%$ by learning a control policy through interactions with a low-cost quadrant photodiode rather than a wavefront phase profiling camera. We develop and share an AO RL environment that provides a standardized platform to develop and evaluate RL based on the Strehl ratio, which is correlated to fiber-coupling performance. Our empirical analysis finds that Proximal Policy Optimization (PPO) outperforms Soft-Actor-Critic and Deep Deterministic Policy Gradient. PPO converges to within $86\%$ of the maximum reward obtained by an idealized Shack-Hartmann sensor after training of 250 episodes, indicating the potential of RL to enable efficient wavefront sensorless OSGC.

Defining the electric and magnetic field vectors in curved spacetime requires a proper choice of the observer's frame four-vector. Related literature shows that this fundamental issue in physics still needs to be properly resolved. In recent literature on using electromagnetic means to detect gravitational waves, an {\it ad hoc} definition based on regarding $F_{ab}$ with two covariant indices as the special relativistic one is popular. We show that by assigning physical fields to tensor components in that way, one cannot identify the frame four-vector allowing such a choice, thus failing to properly define the external charge and current densities in that frame. We propose the normal-frame as the proper one. In this frame, the weak gravity corrections appear as the effective polarizations and magnetizations in both the homogeneous and inhomogeneous parts of Maxwell equations.

The end state of Hawking evaporation of a black hole is uncertain. Some candidate quantum gravity theories, such as loop quantum gravity and asymptotic safe gravity, hint towards Planck sized remnants. If so, the Universe might be filled with remnants of tiny primordial black holes, which formed with mass $M<10^9\,{\rm g}$. A unique scenario is the case of $M\sim 5\times10^5\,{\rm g}$, where tiny primordial black holes reheat the universe by Hawking evaporation and their remnants dominate the dark matter. Here, we point out that this scenario leads to a cosmological gravitational wave signal at frequencies $\sim 100{\rm Hz}$. Finding such a particular gravitational wave signature with, e.g. the Einstein Telescope, would suggest black hole remnants as dark matter.

Axions are considered as a candidate of dark matter. The axions form coherent clouds which delay and amplify gravitational waves (GWs) at a resonant frequency. That is, the coherent axion clouds produce secondary GWs following a primary wave. All GWs from compact binary mergers detected so far propagate in the Milky Way halo composed of dark matter so that such GWs are followed by the secondary GWs. Since the properties of the secondary GWs depend on the axion mass and coupling with the parity violating sector of gravity, we can search for a characteristic signal produced by axions. In our previous study, we have developed a search method optimized for the axion signals, and obtained about ten times stronger constraint on the coupling than the previous best one for the axion mass range, $[1.7 \times 10^{-13},\,8.5 \times 10^{-12}]\,\mathrm{eV}$. However, in our previous search, we assume that dark matter is homogeneously distributed in the Milky Way halo, which is unrealistic. In this paper, we extend the dark matter profile to a more realistic one called the core NFW profile, and find that the constraint in our previous study is robust to the dark matter profile.

The merger remnant of a binary neutron star coalescence is initially strongly differentially rotating. Some properties of these remnants can be accurately modeled through building equilibrium neutron star models. In the present paper, we study how a modification of general relativity, namely scalar-tensor theory with a massive scalar field, will alter the picture. In contrast to previous studies, we implement a realistic phenomenological differential rotational law which allows for neutron star models to attain maximal angular velocity away from the center. We find that solutions with much higher masses and angular momenta exist in scalar-tensor theory compared to general relativity. They keep their quasi-spherical energy-density distribution for significantly higher values of the angular momentum before transitioning to quasi-toroidal models, in contrast to pure general relativity. Constructing such neutron star solutions is the first step to our final goal that is studying how scalarization alters the stability and gravitational wave emission of post-merger remnants.

Recapitulation of the resonance condition for the fundamental and higher electron cyclotron harmonics in the Electron Cyclotron Maser Instability (ECMI) enables radiation below and confirms the possibility of radiation in a narrow band above harmonics $n>1$. Near $n=1$ resonance on the confined lower X-mode branch, amplification is supported by the decrease of phase and group speeds. Confined slow large-amplitude quasi-electrostatic X-modes nonlinearly modulate the plasma to form cavitons until self-trapped inside them at further increasing wavenumber. They undergo wave-wave interaction, enabling escape to free space in the second harmonic band below $n=2$. At sufficiently large parallel wavenumber (oblique propagation), the fundamental resonance $n=1$ is hyperbolic, a possibility so far missed but vital for an effective ECMI in the upward current region. Here, the resonance hyperbola favourably fits the loss cone boundary, the presumably important ECMI upward-current source-electron distribution, to stimulate ECMI growth at available auroral electron energies.