Papers for Friday, Dec 15 2023

We are conducting a comprehensive investigation of V815 Her using photometric and spectroscopic data to understand the origin of the activity and what influences it in the short and long term. Using TESS photometry we performed light curve modeling in order to derive astrophysical and orbital parameters for the eclipsing binary subsystem V815 Her B. Using archival photometric data covering a century we carried out a time frequency analysis. Spectral synthesis was applied to determine the basic astrophysical parameters of the rapidly rotating primary using high-resolution STELLA spectra recorded in 2018. Photometric analysis revealed multiple cycles on timescales between ~6.5 and ~26 years. From TESS photometry we obtained orbital solution for the V815 Her B subsystem. The STELLA spectra covering the 200 day-long observing season enabled to create 19 time-series Doppler images, which revealed a constantly changing spotted surface. From the consecutive image pairs we measured a weak solar-type surface differential rotation of the spotted star. We found evidence that the V815 Her B component previously apostrophized as a third body is actually an eclipsing close binary subsystem of two M dwarfs with a period of 0.5 d, i.e., V815 Her is a 2+2 hierarchical quadruple system. The system is apparently young, only a few times ten million years old, consistent with the spotted primary V815 Her Aa being a zero-age main sequence star. Spot activity on the primary was found to be vivid. Fast starspot decay suggests that convective-turbulent erosion plays a more significant role in such a rapidly rotating star. The weak differential rotation of V815 Her Aa is presumably confined by tidal forces of the close companion V815 Her Ab. The slowly increasing photometric cycle of 6.5 years on average is interpreted as a spot cycle of V815 Her Aa, which is probably modulated by the eccentric wide orbit.

Observations of $z \gtrsim 6$ quasars provide information on the early phases of the most massive black holes (MBHs) and galaxies. Current observations, able to trace both gas and stellar properties, show that most MBHs at high redshift seem overmassive compared to the local population, in line with the elliptical galaxy population, or even above, thus implying a very rapid growth of these objects. To assess the physical conditions for such a rapid growth and explain the existence of a population of already mature MBHs when the Universe was less than 1~Gyr old, we here explore whether episodes of accretion above the Eddington limit can occur across cosmic epochs. By employing state-of-the-art high-resolution cosmological zoom-in simulations of a $z\sim 7$ quasar, where different accretion regimes are consistently included, together with their associated radiative and kinetic feedback, we show that super-Eddington phases can be sustained for relatively long time-scales (tens of Myr) and discuss how they affect the growth of MBHs. We also show, by means of a semi-analytic evolution, that the MBH spin remains relatively low during super-Eddington phases, and this would result in a lower feedback efficiency, hence a potentially faster growth that might then explain the overmassiveness of high-redshift MBHs.

The giant spiral galaxy M101 is host to the nearest recent Type Ia Supernova (SN 2011fe) and thus has been extensively monitored in the near-infrared to study the late-time lightcurve of the supernova. Leveraging this existing baseline of observations, we derive the first Mira-based distance to M101 by discovering and classifying a sample of 211 Miras with periods ranging from 240 to 400 days in the supernova field. Combined with new HST WFC3/IR channel observations, our dataset totals 11 epochs of F110W (HST $YJ$) and 13 epochs of F160W (HST $H$) data spanning $\sim$2900 days. We adopt absolute calibrations of the Mira Period-Luminosity Relation based on geometric distances to the Large Magellanic Cloud and the water megamaser host galaxy NGC 4258, and find $\mu_{\rm M101} = $ 29.10 $\pm$ 0.06 mag. This distance is in 1$\sigma$ agreement with most other recent Cepheid and Tip of the Red Giant Branch distance measurements to M101. Including the previous Mira-SNIa host, NGC 1559 and SN 2005df, we determine the fiducial SN Ia peak luminosity, $M^0_B = -19.27 \pm 0.09$ mag. With the Hubble diagram of SNe Ia, we derive $H_0 = 72.37 \pm 2.97 $ km s$^{-1}$Mpc$^{-1}$, a $4.1\%$ measurement of $H_0$ using Miras. We find excellent agreement with recent Cepheid distance ladder measurements of $H_0$ and confirm previous indications that the local universe value of $H_0$ is higher than the early-universe value at $\sim$ $95\%$ confidence. Currently, the Mira-based $H_0$ measurement is still dominated by the statistical uncertainty in the SN Ia peak magnitude.

In this work, we look for evidence of a non-unity mass ratio binary dwarf galaxy merger in the Sagittarius stream. Simulations of such a merger show that, upon merging with a host, particles from the less-massive galaxy will often mostly be found in the extended stream and less-so in the central remnant. Motivated by these simulations, we use APOGEE DR17 chemical data from approximately 1100 stars in both the Sagittarius remnant and stream to look for evidence of contamination from a second dwarf galaxy. This search is initially justified by the idea that disrupted binary dwarf galaxies provide a possible explanation of the Sagittarius bifurcation, and the location of the massive, chemically peculiar globular cluster NGC 2419 found within the stream of Sagittarius. We separate the Sagittarius data into its remnant and stream and compare the [Mg/Fe] content of the two populations. In particular, we select [Mg/Fe] to search for hints of unique star formation histories among our sample stars. Comparing the stream and remnant populations, we find regions have distinct [Mg/Fe] distributions for fixed [Fe/H], in addition to distinct chemical tracks in [Mg/Fe] -- [Fe/H] abundance space. We show that there are large regions of the tracks for which the probability of the two samples being drawn from the same distribution is very low (p < 0.05). Furthermore, we show that the two tracks can be fit with unique star formation histories using simple, one zone galactic chemical evolution models. While more work must be done to discern whether the hypothesis presented here is true, our work hints at the possibility that Sagittarius may consist of two dwarf galaxy progenitors.

Measurements of chemical abundances in high-$z$ star-forming (SF) galaxies place important constraints on the enrichment history of massive stars and the physical conditions in the early universe. JWST is beginning to enable direct chemical abundance measurements in galaxies at $z$$>$2 via the detection of the faint T$_e$-sensitive auroral line [O III]$\lambda$4363. However, direct abundances of other elements (e.g., S and Ar) in high-$z$ galaxies remain unconstrained due to a lack of T$_e$ data and wavelength coverage. Here, we present multiple direct abundances in D40, a galaxy at $z$$\sim$3 observed with JWST/NIRSpec as part of the CECILIA program. We report the first simultaneous measurement of T$_e$[O III] and T$_e$[S III] in a high-$z$ galaxy, finding good agreement with the temperature trends in local SF systems and photoionization models. We measure a gas-phase metallicity of 12+log(O/H) $=8.16\pm0.05$ dex, and an N/O abundance indicative of primary nucleosynthesis. The S/O abundance in D40 is slightly sub-solar but consistent with local H II regions. In contrast, the log(Ar/O) in D40 is $-$2.71$\pm$0.09 dex, sub-solar by $>$2$\sigma$. The [Ar III]$\lambda$7135 intensities of D40 and other CECILIA galaxies are similar to those of local SF systems with Ar-deficient interstellar media, suggesting that low gas-phase Ar abundance is common in high-$z$ galaxies. Recent nucleosynthesis models find that Ar is produced in Type Ia supernovae: if low Ar/O is consistently observed in high-$z$ SF galaxies, it would provide further observational evidence that early galaxies are enriched primarily by core-collapse supernovae, which produce enhanced O relative to Ar and Fe.

Binary evolution codes are essential tools to help in understanding the evolution of binary systems. They contain a great deal of physics, for example stellar evolution, stellar interactions, mass transfer, tides, orbital evolution. Since many of these processes are difficult to account for in detail, we often rely on prescriptions obtained in earlier studies. We highlight that the impact of the dynamical tides with radiative damping has been implemented inconsistently with respect to its original theoretical formulation in many studies. We derive a new analytical solution for the evolution toward synchronization in the case of circular orbits and propose turnkey equations for the case of eccentric orbits that can be used in population synthesis studies. We compare the strength of the tidal torque obtained with this new formula with respect to that obtained with the formula generally used in literature by studying how the evolution toward synchronization of main sequence stellar models is affected. We conclude that by using an incorrect formula for the tidal torque, as has been done in many binary codes, the strength of the dynamical tides with radiative damping is over- or underestimated depending on whether the star is close to or far from synchronization.

We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC~4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas temperatures and number densities peak in the central regions and decrease along the outflows. These profiles are inconsistent with a spherical, adiabatically-expanding wind model, suggesting the need to include mass loading and/or a non-spherical outflow geometry. We estimate the mass outflow rate of the hot wind to be $2.1\:M_{\odot}~\rm{yr}^{-1}$. Emission from charge exchange is detected in the northern outflow, and we estimate it contributes 12\% to the emitted, broad-band ($0.5-7$~keV) X-ray flux.

A characteristic feature of quasars is the observed non-linear relationship between their monochromatic luminosities at rest-frame 2500 {\AA} and 2 keV. This relationship is evident across all redshifts and luminosities and, due to its non-linearity, can be implemented to estimate quasar distances and construct a Hubble Diagram for quasars. Historically, a significant challenge in the cosmological application of this relation has been its high observed dispersion. Recent studies have demonstrated that this dispersion can be reduced by excluding biased objects from the sample. Nevertheless, the dispersion remains considerable ($\delta \sim 0.20$ dex), especially when compared to the Phillips relation for supernovae Ia. Given the absence of a comprehensive physical model for the relation, it remains unclear how much of the remaining dispersion is intrinsically tied to the relation itself and how much can be attributed to observational factors not addressed by the sample selection and by the choice of X-ray and UV indicators. Potential contributing factors include (i) the scatter produced by using X-ray photometric results instead of spectroscopic ones, (ii) the intrinsic variability of quasars, and (iii) the inclination of the accretion disc relative to our line of sight. In this study, we thoroughly examine these three factors and quantify their individual contributions to the observed dispersion. Based on our findings, we argue that the intrinsic dispersion of the X-ray/UV luminosity relation is likely below 0.06 dex. We also discuss why high-redshift subsamples can show a significantly lower dispersion than the average one.

Weak gravitational lensing convergence peaks, the local maxima in weak lensing convergence maps, have been shown to contain valuable cosmological information complementary to commonly used two-point statistics. To exploit the full power of weak lensing for cosmology, we must model baryonic feedback processes because these reshape the matter distribution on non-linear and mildly non-linear scales. We study the impact of baryonic physics on the number density of weak lensing peaks using the FLAMINGO cosmological hydrodynamical simulation suite. We generate ray-traced full-sky convergence maps mimicking the characteristics of a Stage IV weak lensing survey. We compare the number densities of peaks in simulations that have been calibrated to reproduce the observed galaxy mass function and cluster gas fraction or to match a shifted version of these, and that use either thermally driven or jet AGN feedback. We show that the differences induced by realistic baryonic feedback prescriptions (typically $5-30$% for $\kappa = 0.1-0.4$) are smaller than those induced by reasonable variations in cosmological parameters ($20-60$% for $\kappa = 0.1-0.4$) but must be modeled carefully to obtain unbiased results. The reasons behind these differences can be understood by considering the impact of feedback on halo masses, or by considering the impact of different cosmological parameters on the halo mass function. Our analysis demonstrates that the baryonic suppression is insensitive to changes in cosmology up to $\kappa \approx 0.4$ and that the higher $\kappa$ regime is dominated by Poisson noise and cosmic variance.

The supermassive black holes (SMBHs) that power active galactic nuclei found at $z\geq 6$ were formed during the epoch of reionization. Because reionization is a spatially inhomogeneous process, where different regions of the Universe become reionized at different times, the physical properties of SMBH host galaxy environments will vary spatially during reionization. We construct a semi-analytic model to estimate the impact of reionization on SMBH growth due to reduced gas accretion onto dark matter halos. Using a series of merger trees, reionization models, and black hole growth models, we find that early reionization can reduce an SMBH's mass by up to [50, 70, 90] % within dark matter halos of mass [$10^{12}$, $10^{11}$, $10^{10}$] M$_{\odot}$ by $z$ = 6. Our findings also suggest that the redshift range in which black hole growth is impacted by reionization is strongly dependent on whether the Eddington accretion rate can be exceeded. If so, we find that black hole masses are significantly suppressed principally during the early phases of reionization ($z$ $\gtrsim$ 10), while they are more readily suppressed across the full redshift range if super-Eddington growth is not allowed. We find that the global average impact of reionization may be to reduce the masses of black holes residing in $\lesssim$ 10$^{11}$ M$_{\odot}$ halos by a factor of $\gtrsim$ 2. The census of SMBHs that the James Webb Space Telescope is uncovering provides a promising means by which to test these predictions.

A neutron star harbors of order $10^{56}$ electrons in its core, and almost the same number of muons, with muon decay prohibited by Pauli blocking. However, as macroscopic properties of the star such as its mass, rotational velocity, or magnetic field evolve over time, the equilibrium lepton abundances (dictated by the weak interactions) change as well. Scenarios where this can happen include spin-down, accretion, magnetic field decay, and tidal deformation. We discuss the mechanisms by which a star disrupted in one of these ways re-establishes lepton chemical equilibrium. In most cases, the dominant processes are out-of-equilibrium Urca reactions, the rates of which we compute for the first time. If, however, the equilibrium muon abundance decreases, while the equilibrium electron abundance increases (or decreases less than the equilibrium muon abundance), outward diffusion of muons plays a crucial role as well. This is true in particular for stars older than about 10,000 yrs whose core has cooled to $\lesssim 20$ keV. The muons decay in a region where Pauli blocking is lifted, and we argue that these decays lead to a flux of $\mathcal{O}$(10 MeV) neutrinos. Realistically, however, this flux will remain undetectable for the foreseeable future.

Located at the core of the Galactic Centre, the S-star cluster serves as a remarkable illustration of chaos in dynamical systems. The long-term chaotic behaviour of this system can be studied with gravitational $N$-body simulations. By applying a small perturbation to the initial position of star S5, we can compare the evolution of this system to its unperturbed evolution. This results in the two solutions diverging exponentially, defined by the separation in position space $\delta_{r}$, with an average Lyapunov timescale of $\sim$420 yr, corresponding to the largest positive Lyapunov exponent. Even though the general trend of the chaotic evolution is governed in part by the supermassive black hole Sagittarius $\rm A^{*}$ (Sgr $\rm A^{*}$), individual differences between the stars can be noted in the behaviour of their phase-space curves. We present an analysis of the individual behaviour of the stars in this Newtonian chaotic dynamical system. The individuality of their behaviour is evident from offsets in the position space separation curves of the S-stars and the black hole. We propose that the offsets originate from the initial orbital elements of the S-stars, where Sgr $\rm A^{*}$ is considered in one of the focal points of the Keplerian orbits. Methods were considered to find a relationship between these elements and the separation in position space. Symbolic regression turns out to provide the clearest diagnostics for finding an interpretable expression for the problem. Our symbolic regression model indicates that $\left\langle\delta_r\right\rangle \propto e$, implying that the time-averaged individual separation in position space is directly proportional to the initial eccentricity of the S-stars.

The variations in Ly$\alpha$ forest opacity observed at $z>5.3$ between lines of sight to different background quasars are too strong to be caused by fluctuations in the density field alone. The leading hypothesis for the cause of this excess variance is a late, ongoing reionization process at redshifts below six. Another model proposes strong ionizing background fluctuations coupled to a short, spatially varying mean free path of ionizing photons, without explicitly invoking incomplete reionization. With recent observations suggesting a short mean free path at $z\sim6$, and a dramatic improvement in $z>5$ Ly$\alpha$ forest data quality, we revisit this latter possibility. Here we apply the likelihood-free inference technique of approximate Bayesian computation to jointly constrain the hydrogen photoionization rate $\Gamma_{\rm HI}$ and the mean free path of ionizing photons $\lambda_{\rm mfp}$ from the effective optical depth distributions at $z=5.0$-$6.1$ from XQR-30. We find that the observations are well-described by fluctuating mean free path models with average mean free paths that are consistent with the steep trend implied by independent measurements at $z\sim5$-$6$, with a concomitant rapid evolution of the photoionization rate.

Line intensity mapping (LIM) serves as a potent probe in astrophysics, relying on the statistical analysis of integrated spectral line emissions originating from distant star-forming galaxies. While LIM observations hold the promise of achieving a broad spectrum of scientific objectives, a significant hurdle for future experiments lies in distinguishing the targeted spectral line emitted at a specific redshift from undesired line emissions originating at different redshifts. The presence of these interloping lines poses a challenge to the accuracy of cosmological analyses. In this study, we introduce a novel approach to quantify line-line cross-correlations (LIM-LLX), enabling us to investigate the true signal amidst instrumental noise and interloping emissions. For example, at a redshift of approximately $z\sim3.7$, we observed that the measured auto-power spectrum of [CII] exhibited substantial bias, from interloping line emission. However, cross-correlating [CII] with CO(6-5) lines using a FYST-like experiment yielded a promising result, with a Signal-to-noise ratio (SNR) of $\sim 10$. This measurement is notably unbiased. Additionally, we explore the extensive capabilities of cross-correlation by leveraging various CO transitions to probe the tomographic Universe at lower redshifts through LIM-LLX. We further demonstrate that incorporating low-frequency channels, such as 90 GHz and 150 GHz, into FYST's EoR-Spec-like experiment can maximize the potential for cross-correlation studies, effectively reducing the bias introduced by instrumental noise and interlopers.

The cross-correlation of cosmic voids with the lensing convergence ($\kappa$) map of the CMB fluctuations offers a powerful tool to refine our understanding of the dark sector in the consensus cosmological model. Our principal aim is to compare the lensing signature of our galaxy data set with simulations based on the concordance model and characterize the results with an $A_{\kappa}$ consistency parameter. In particular, our measurements contribute to the understanding of the "lensing-is-low" tension of the $\Lambda$CDM model. We selected luminous red galaxies from the WISE-Pan-STARSS data set, allowing an extended 14,200 deg$^2$ sky area, that offers a more precise measurement compared to previous studies. We created 2D and 3D void catalogs to cross-correlate their locations with the Planck lensing map and studied their average imprint signal using a stacking methodology. Applying the same procedure, we also generated a mock catalog from the WebSky simulation for comparison. The 2D void analysis revealed good agreement with the standard cosmological model with $A_{\kappa}\approx1.06 \pm 0.08$, i.e. $S/N=13.3$, showing a higher $S/N$ than previous studies using voids detected in the Dark Energy Survey data set. The 3D void analysis exhibited a lower $S/N$ and demonstrated worse agreement with our mock catalog than the 2D voids. These deviations might be attributed to limitations in the mock catalog, such as imperfections in the LRG selection, as well as a potential asymmetry between the North and South patches of the WISE-Pan-STARSS data set in terms of data quality. Overall, we present a significant detection of a CMB lensing signal associated with cosmic voids, largely consistent with the concordance model. Future analyses using even larger data sets also hold great promise of further sharpening these results, given their complementary nature to large-scale structure analyses.

Elastic scattering of dark matter (DM) particles with baryons induce cosmological signals that may be detectable with modern or future telescopes. For DM-baryon scattering cross sections scaling with negative powers of relative velocity, $\sigma_{\chi b}(v) \propto v^{-2}, v^{-4}$, such interactions introduce a momentum-exchange rate that is nonlinear in DM-baryon bulk relative velocities, thus not amenable for inclusion as-is into standard linear cosmological Boltzmann codes. Linear ansatzes have been adopted in past works, but their accuracy is unknown as they do not arise from first-principles derivations. In this work, for the first time, we construct a rigorous framework for computing linear-cosmology observables as a perturbative expansion in $\sigma_{\chi b}$. We argue that this approach is accurate for Cosmic Microwave Background (CMB) angular power spectra when most or all of the DM is scattering with baryons with cross section $\sigma_{\chi b}(v) \propto v^{-2}, v^{-4}$. We derive exact formal expressions for CMB power spectra at linear order in $\sigma_{\chi b}$, and show that they only depend on a specific velocity integral of the momentum-exchange rate. Consequently, we can obtain the exact power spectra at linear order in $\sigma_{\chi b}$ by substituting the original nonlinear momentum-exchange rate with a uniquely specified linear rate. Serendipitously, we find that the exact substitution we derive from first principles precisely coincides with the most widely used linear ansatz, thus placing previous CMB-anisotropy upper bounds on a more solid footing. In addition to finally providing an exact cosmological solution to the DM-baryon scattering problem in a well-defined region of parameter space, the framework we construct opens the way to computing higher-order correlation functions, beyond power spectra, which are promising yet unexplored probes of DM-baryon scattering.

Context. Double eclipsing binaries are gravitationally bound quadruple systems in a 2+2 configuration where both of the binaries are eclipsing. These systems are interesting objects to better understand stellar formation, to investigate the dynamical interaction between the two binary systems or to study certain stages of stellar evolution. Aims. With this work, we aim to determine if double eclipsing binaries can be found using ZTF data and what the difficulties are in doing so. Secondly, we aim to significantly increase the number of known double eclipsing systems and determine how this sample differs from samples of double eclipsing binaries found with other telescopes. Methods. We develop a new method to systematically search for double eclipsing binaries in sparsely sampled light curves. For this we use box-least-squares (BLS) to search for the period of the first binary in the system. We then remove that signal from the light curves, and search the residual light curve again with BLS to find the second period. We applied this method to ZTF light curves of 575 526 eclipsing binaries known in the Gaia eclipsing binary catalogue. Results. We report the discovery of 198 new double eclipsing binary systems. Conclusions. We successfully implemented a method that systematically searches for double eclipsing binary systems in sparsely sampled data. In total 198 new double eclipsing binary systems have been found in 575 526 light curves. The ZTF sample typically contains more short period binaries compared to the TESS sample, but is also able to find systems with longer periods. We expect that at least three to four times more quadruples can be found by applying this method to all ZTF stellar light curves, by increasing the number of data points as a result of longer observations, and by implementing an automatic detection mechanism that replaces visual inspection.

We present a comprehensive analysis of statistical tools for evaluating tensions in cosmological parameter estimates arising from distinct datasets. Focusing on the unresolved Hubble constant ($H_0$) tension, we explore the Pantheon Plus + SH0ES (PPS) compilation, which includes low-redshift Cepheid data from the SH0ES collaboration, and the Cosmic Chronometers (CC) dataset. Employing various tension metrics, we quantitatively assess the inconsistencies in parameter estimates, emphasizing the importance of capturing multidimensional tensions. Our results reveal substantial tension between PPS and Planck 2018 datasets. We highlight the importance of the adoption of these metrics to enhance the precision of future cosmological analyses and facilitate the resolution of existing tensions.

M 39 is a nearby young open cluster hardly studied in the last decades. No giant is known among its members and its chemical composition has never been studied. In order to investigate it we performed high-resolution spectroscopy of 20 expected cluster members with the HARPS and FIES spectrographs. By combining our observations with archival photometry and $Gaia$-DR3 data we searched for evolved members and studied cluster properties such as the radial velocity, extinction and age. For the first time, we provide stellar parameters and chemical abundances for 21 species with atomic numbers up to 56. We have not found any new giant as likely member and notice a negligible reddening along the cluster field, that we place at 300 pc. We obtain a mean radial velocity for M 39 of -5.5$\pm$0.5 km s$^{-1}$ and an isochrone-fitting age of 430$\pm$110 Ma, which corresponds to a MSTO mass of around 2.8 Msol. This value is consistent with the Li content and chromospheric activity shown by its members. Based on main-sequence stars the cluster exhibits a solar composition, [Fe/H]=+0.04$\pm$0.08 dex, compatible with its Galactocentric location. However, it has a slightly subsolar abundance of Na and an enriched content of neutron-capture elements, specially Ba. In any case, the chemical composition of M 39 is fully compatible with that shown by other open clusters that populate the Galactic thin disc

The gravitational microlensing technique is most sensitive to planets in a Jupiter-like orbit and has detected more than 200 planets. However, only a few wide-orbit ($s > 2$) microlensing planets have been discovered, where $s$ is the planet-to-host separation normalized to the angular Einstein ring radius, $\theta_{\rm E}$. Here we present the discovery and analysis of a strong candidate wide-orbit microlensing planet in the event, OGLE-2017-BLG-0448. The whole light curve exhibits long-term residuals to the static binary-lens single-source model, so we investigate the residuals by adding the microlensing parallax, microlensing xallarap, an additional lens, or an additional source. For the first time, we observe a complex degeneracy between all four effects. The wide-orbit models with $s \sim 2.5$ and a planet-to-host mass-ratio of $q \sim 10^{-4}$ are significantly preferred, but we cannot rule out the close models with $s \sim 0.35$ and $q \sim 10^{-3}$. A Bayesian analysis based on a Galactic model indicates that, despite the complicated degeneracy, the surviving wide-orbit models all contain a super-Earth-mass to Neptune-mass planet at a projected planet-host separation of $\sim 6$ au and the surviving close-orbit models all consist of a Jovian-mass planet at $\sim 1$ au. The host star is probably an M or K dwarf. We discuss the implications of this dimension-degeneracy disaster on microlensing light-curve analysis and its potential impact on statistical studies.

We present total intensity and linear polarization images of OJ287 at 1.68GHz, obtained through space-based VLBI observations with RadioAstron on April 16, 2016. The observations were conducted using a ground array consisting of the VLBA and the EVN. Ground-space fringes were detected with a maximum projected baseline length of 5.6 Earth's diameter, resulting in an angular resolution of 530 uas. With this unprecedented resolution at such a low frequency, the progressively bending jet structure of OJ287 has been resolved up to 10 pc of the projected distance from the radio core. In comparison with close-in-time VLBI observations at 15, 43, 86 GHz from MOJAVE and VLBA-BU-BLAZAR monitoring projects, we obtain the spectral index map showing the opaque core and optically thin jet components. The optically thick core has a brightness temperature of 10$^{13}$ K, and is further resolved into two sub-components at higher frequencies labeled C1 and C2. These sub-components exhibit a transition from optically thick to thin, with a SSA turnover frequency estimated to be 33 and 11.5 GHz, and a turnover flux density 4 and 0.7 Jy, respectively. Assuming a Doppler boosting factor of 10, the SSA values provide the estimate of the magnetic field strengths from SSA of 3.4 G for C1 and 1.0 G for C2. The magnetic field strengths assuming equipartition arguments are also estimated as 2.6 G and 1.6 G, respectively. The integrated degree of linear polarization is found to be approximately 2.5 %, with the electric vector position angle being well aligned with the local jet direction at the core region. This alignment suggests a predominant toroidal magnetic field, which is in agreement with the jet formation model that requires a helical magnetic field anchored to either the black hole ergosphere or the accretion disk. Further downstream, the jet seems to be predominantly threaded by a poloidal magnetic field.

The recently discovered population of faint FR0 radiogalaxies has been interpreted as the extension to low power of the classical FRI sources. Their radio emission appears to be concentrated in very compact (pc-scale) cores, any extended emission is very weak or absent and VLBI observations show that jets are already mildly or sub-relativistic at pc scales. Based on these observational properties we propose here that the jets of FR0s are strongly decelerated and disturbed at pc scale by hydrodynamical instabilities. With the above scenario in mind, we study the dynamics of a low-power relativistic jet propagating into a confining external medium, focusing on the effects of entrainment and mixing promoted by the instabilities developing at the jet-environment interface downstream of a recollimation shock. We perform a 3D relativistic hydrodynamical simulation of a recollimated jet by means of the state-of-the-art code PLUTO. The jet is initially conical, relativistic (with initial Lorentz Factor $\Gamma$=5), cold and light with respect to the confining medium, whose pressure is assumed to slowly decline with distance. The magnetic field is assumed to be dynamically unimportant.The 3D simulation shows that, after the first recollimation/reflection shock system, a rapidly growing instability develops, as a result of the interplay between Kelvin-Helmholtz and Richtmyer-Meshkov modes. In turn, the instability promotes strong mixing and entrainment that rapidly lead to the deceleration of the jet and spread its momentum to slowly moving, highly turbulent external gas. We argue that this mechanism could account for the peculiarities of the low-power FR0 jets. For outflows with higher power, Lorentz factor or magnetic field, we expect that the destabilizing effects are less effective, allowing the survival of the jet up to the kpc scale, as observed in FRIs.

Using multi-instrument and multi-wavelength observations, we studied a CME eruption that led to intense geomagnetic storm on 23 April 2023. The eruption occurred on April 21 in solar active region 13283 near the disk-center. The AR was in its decay stage, with fragmented polarities and a pre-existing long filament channel, a few days before the eruption. The study of magnetic field evolution suggests that the flux-rope (filament) has been built up by monotonous helicity accumulation over several days, and further converging and canceling fluxes lead to helicity injection change, resulting in the unstable nature of the magnetic flux-rope (MFR) and its further eruption. Importantly, the CME morphology revealed that the MFR apex underwent a rotation of upto 56\degree~in clockwise-direction owing to its positive helicity. The CME decelerates in the LASCO-FOV and has a plane-of-sky velocity of 1226 km/s at 20\,R$_\odot$. In the Heliospheric Imager FOV, the CME lateral expansion is tracked more than the earthward motion. This implies that the arrival time estimation is difficult to assess. The in-situ arrival of ICME shock was at 07:30 UT on April 23, and a geomagnetic storm commenced at 08:30\,UT. The flux rope fitting to the in-situ magnetic field observations reveals that the MC flux rope orientation is consistent with its near Sun orientation, which has a strong negative Bz-component. The analysis of this study indicates that the near-Sun rotation of the filament during its eruption to the CME is the key to the negative Bz-component and consequently the intense geomagnetic storm.

In recent years, a number of Lyman continuum (LyC) leaker candidates at intermediate redshifts have been found, providing insight into how the Universe was reionised at early cosmic times. Here we identify new LyC leaker candidates at $z\approx 3-4.5$ and compare them to objects from the literature to get an overview of the different observed escape fractions and their relation to the properties of the Lyman $\alpha$ (Ly$\alpha$) emission line. The aim of this work is to test indicators for LyC leakage and to improve our understanding of the kind of galaxies from which LyC radiation can escape. We use data from the Hubble Deep Ultraviolet (HDUV) legacy survey to search for LyC emission based on a sample of $\approx 2000$ Ly$\alpha$ emitters (LAEs) detected previously in two surveys with the Multi-Unit Spectroscopic Explorer (MUSE), MUSE-Deep and MUSE-Wide. Based on their known redshifts and positions, we look for potential LyC leakage in the WFC3/UVIS F336W band of the HDUV. The escape fractions are measured and compared based on spectral energy distribution (SED) fitting performed using the CIGALE software. We add twelve objects to the sample of known LyC leaker candidates, one of which was previously known, and compare their Ly$\alpha$ properties to their escape fractions. We find escape fractions between $\sim 20\%$ and $\sim 90\%$, assuming a high transmission in the intergalactic medium (IGM). We show a method to use the number of LyC leaker candidates we find to infer the underlying average escape fraction of galaxies, which is $\approx 12\%$. Based on their Ly$\alpha$ properties, we conclude that LyC leakers are not very different from other high-z LAEs and suggest that most LAEs could be leaking LyC even if this can not always be detected due to the direction of emission and the transmission properties of the IGM.

Follow-up of gravitational wave alerts has proven to be challenging in the past due to the large uncertainty on the localisation, much larger than the field of view of most instruments. A smart pointing strategy helps to increase the chance of observing the true position of the underlying compact binary merger event and so to detect an electromagnetic counterpart. To tackle this, a software called tilepy has been developed and was successfully used by the H.E.S.S. collaboration to search for very-high energy gamma-ray emission from GWs during the O2 and O3 runs. The optimised tiling strategies implemented in tilepy allowed H.E.S.S. to be the first ground facility to point toward the true position of GW 170817. Here we present the main strategy used by the software to compute an optimal observation schedule. The Astro-COLIBRI platform helps to plan follow-up of a large range of transient phenomena including GW alerts. The integration of tilepy in this tool allow for an easy planning and visualisation of of follow-up of gravitational wave alert helping the astronomer to maximise the chance of detecting a counterpart. The platform also provides an overview of the multi-wavelength context by grouping and visualising information coming from different observatories alongside GW alerts.

We describe how a novel online changepoint detection algorithm, called Poisson-FOCuS, can be used to optimally detect gamma-ray bursts within the computational constraints imposed by miniaturized satellites such as the upcoming HERMES-Pathfinder constellation. Poisson-FOCuS enables testing for gamma-ray burst onset at all intervals in a count time series, across all timescales and offsets, in real-time and at a fraction of the computational cost of conventional strategies. We validate an implementation with automatic background assessment through exponential smoothing, using archival data from Fermi-GBM. Through simulations of lightcurves modeled after real short and long gamma-ray bursts, we demonstrate that the same implementation has higher detection power than algorithms designed to emulate the logic of Fermi-GBM and Compton-BATSE, reaching the performances of a brute-force benchmark with oracle information on the true background rate, when not hindered by automatic background assessment. Finally, using simulated data with different lengths and means, we show that Poisson-FOCuS can analyze data twice as fast as a similarly implemented benchmark emulator for the historic Fermi-GBM on-board trigger algorithms.

We aim to investigate the origin of the HI Br$\gamma$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$\gamma$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of HD58647 with both high spectral and high angular resolution. Thanks to the extensive $uv$ coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Br$\gamma$-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Br$\gamma$-emitting region. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14\degr. The size of the continuum emission at 5-sigma levels is $\sim$ 1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas $\pm$ 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD58647 traced by the Br$\gamma$ ,has been imaged. This allowed us to constrain the size of the Br$\gamma$-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas $\pm$ 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk.

Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainties in foregrounds and high redshift astrophysics. We include the latest publicly available 21-cm power spectrum upper limits from HERA, 21-cm global signal data from SARAS 3, and the synergistic probe of the unresolved X-ray background in our final analysis. This paper thus constitutes the first attempt to use 21-cm power spectrum data to probe cosmic strings. In contrast to previous works, we find no strong constraints can be placed on superconducting cosmic strings from current 21-cm measurements. This is because of uncertainties in the X-ray emission efficiency of the first galaxies, with X-ray emissivities greater than $3 \times 10^{40}$erg s$^{-1}$ M$_{\odot}^{-1}$ yr able to mask the presence of cosmic strings in the 21-cm signal. We conclude by discussing the prospects for future constraints from definitive 21-cm signal measurements and argue that the recently proposed soft photon heating should be cause for optimism due to its potential to break degeneracies that would have otherwise made the signatures of cosmic strings difficult to distinguish from those of astrophysical origin.

AT2018cqh is a unique tidal disruption event (TDE) candidate discovered in a dwarf galaxy. Both the light curve fitting and galaxy scaling relationships suggest a central black hole mass in the range of 5.9<logM_BH/M_sun<6.4. A delayed X-ray brightening was found around 590 days after the optical discovery, but shows unusual long-time rising to peak over at least 558 days, which could be coming from delayed accretion of a newly forming debris disk. We report the discovery of delayed radio flares around 1105 days since its discovery, characterized by an initial steep rise of ~>175 days, a flattening lasting about 544 days, and a phase with another steep rise. The rapid rise in radio flux coupled with the slow decay in the X-ray emission points to a delayed launching of outflow, perhaps due to a transition in the accretion state. However, known accretion models can hardly explain the origins of the secondary radio flare that is rising even more rapidly in comparison with the initial one. If confirmed, AT2018cqh would be a rare TDE in a dwarf galaxy exhibiting optical, X-ray and radio flares. We call for continued multi-frequency radio observations to monitor its spectral and temporal evolution, which may help to reveal new physical processes that are not included in standard TDE models.

We present $\texttt{Miko}$, a catalog-to-cosmology pipeline for general flat-sky field-level inference, which provides access to cosmological information beyond the two-point statistics. In the context of weak lensing, we identify several new field-level analysis systematics (such as aliasing, Fourier mode-coupling, and density-induced shape noise), quantify their impact on cosmological constraints, and correct the biases to a percent level. Next, we find that model mis-specification can lead to both absolute bias and incorrect uncertainty quantification for the inferred cosmological parameters in realistic simulations. The Gaussian map prior infers unbiased cosmological parameters, regardless of the true data distribution, but it yields over-confident uncertainties. The log-normal map prior quantifies the uncertainties accurately, although it requires careful calibration of the shift parameters for unbiased cosmological parameters. We demonstrate systematics control down to the $2\%$ level for both models, making them suitable for ongoing weak lensing surveys.

We report on the timing characteristics of MXB 0656-072 throughout its 2007-2008 type I outbursts utilising RXTE/PCA and Fermi/GBM data. Using pulse timing technique, we explore the spin frequency evolution of the source during this interval. Subsequently, by examining the torque-luminosity relation, we show that the overall frequency evolution is substantially in line with the Ghosh-Lamb model. Furthermore, the residuals of the spin frequencies do not exhibit clear orbital modulations, which possibly indicate that the system is observed on a relatively top view. In the RXTE/PCA observations, the pulsed emission is found to be disappearing below $\sim$$5 \times 10^{36}$ erg s$^{-1}$, whereas the profiles maintain stability above this value within our analysis timeframe. In addition, we incorporate two novel methods along with the conventional Deeter method in order to generate higher-resolution power density spectra (PDS). A red noise pattern in the PDSs is also verified in these new methods, common in disk-fed sources, with a steepness of $\Gamma \sim -2$, reaching saturation at a time-scale of $\sim$150 d. Considering the models for spectral transitions, we discuss the possible scenarios for the dipolar magnetic field strength of MXB 0656-072 and its coherence with deductions from the cyclotron resonance scattering feature (CRSF).

Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which aims to investigate the influence of the starburst environment on the formation of stars and planets, and on the evolution of both low and high mass stars. The primary targets of this project are Westerlund 1 and 2, the closest supermassive star clusters to the Sun. Methods. The project is based primarily on recent observations conducted with the Chandra and JWST observatories. Specifically, the Chandra survey of Westerlund 1 consists of 36 new ACIS-I observations, nearly co-pointed, for a total exposure time of 1 Msec. Additionally, we included 8 archival Chandra/ACIS-S observations. This paper presents the resulting catalog of X-ray sources within and around Westerlund 1. Sources were detected by combining various existing methods, and photon extraction and source validation were carried out using the ACIS-Extract software. Results. The EWOCS X-ray catalog comprises 5963 validated sources out of the 9420 initially provided to ACIS-Extract, reaching a photon flux threshold of approximately 2x10^-8 photons/cm^2/s. The X-ray sources exhibit a highly concentrated spatial distribution, with 1075 sources located within the central 1 arcminute. We have successfully detected X-ray emissions from 126 out of the 166 known massive stars of the cluster, and we have collected over 71000 photons from the magnetar CXO J164710.20-455217

We present a new method for modelling the kinematics of galaxies from interferometric observations by performing the optimization of the kinematic model parameters directly in visibility-space instead of the conventional approach of fitting velocity fields produced with the CLEAN algorithm in real-space. We demonstrate our method on ALMA observations of $^{12}$CO (2$-$1), (3$-$2) or (4$-$3) emission lines from an initial sample of 30 massive 850$\mu$m-selected dusty star-forming galaxies with far-infrared luminosities $\gtrsim$$\,10^{12}\,$L$_{\odot}$ in the redshift range $z \sim\,$1.2$-$4.7. Using the results from our modelling analysis for the 12 sources with the highest signal-to-noise emission lines and disk-like kinematics, we conclude the following: (i) Our sample prefers a CO-to-$H_2$ conversion factor, of $\alpha_{\rm CO} = 0.92 \pm 0.36$; (ii) These far-infrared luminous galaxies follow a similar Tully$-$Fisher relation between the circularized velocity, $V_{\rm circ}$, and baryonic mass, $M_{\rm b}$, as more typical star-forming samples at high redshift, but extend this relation to much higher masses $-$ showing that these are some of the most massive disk-like galaxies in the Universe; (iii) Finally, we demonstrate support for an evolutionary link between massive high-redshift dusty star-forming galaxies and the formation of local early-type galaxies using the both the distributions of the baryonic and kinematic masses of these two populations on the $M_{\rm b}\,-\,\sigma$ plane and their relative space densities.

We analyzed TESS photometric data of the flickering-active cataclysmic star AQ Men in 2018--2019. We processed 7 sectors with 14 light curves (LCs) inside them, with a time resolution of 2 min. Aiming to study the "macro-flickering", with quasi periods (QPs) between 10 and 100 hours, we processed LCs after 55 time reduced, with a time resolution of 1.83 hours. The method, developed earlier by us, includes comparing the LCs by their statistical and fractal parameters, as well as revealing QPs by minima of structure functions and relevant maxima of autocorrelation functions. We distinguish the known high state of AQ Men, in the sectors ## 01, 05, 08, 12, 13, as well as the low state, in sectors #19, #20. In the low state the LCs show noticeable eclipses with period 3.4 h, lower average fluxes, higher scatters, and additional QPs. By its statistical and fractal parameters, the macro-flickering of AQ Men in the high state is similar to the ordinary flickering of 3 symbiotic binaries, studied by us earlier (see below). We found 92 QPs in the range of 20--70 h. We reveal 3 QP modes, at 20.9 h, 32.5 h and 54.1 h (1.149, 0.738 and 0.434 c/d; Fig. 7) within a standard error about 10%. The last mode is the most populated one and seems to be a manifestation of the superorbital period. Other 4 QP modes of AQ Men are added from the literature. The regularity of these 7 QP modes follows a power function with a base 1.57 and standard deviation 6.4% (Fig. 8). This power model prognosticates 5 other QP modes: 3 internal and 2 external (Table 2). The bases of the power regularity models for the flickering of the symbiotic binaries RS Oph, T CrB, and MWC 560 (however in the time scale of minutes) are 1.55, 2.0, and 1.34, respectively (Table 1). For unknown reasons in these 4 cases we find (i) regularities with (ii) different bases.

The particle nature of dark matter (DM) has been a long-lasting mystery. Many models suggest that DM could decay or self annihilate into standard model particles, and thus could be a source of gamma rays in the sky. The High Altitude Water Cherenkov (HAWC) observatory has yielded some of the strongest limits in searches of DM decay or annihilation. Building on the flux limits provided by the HAWC collaboration in 2018, we consider the effects of additional components from Galactic secondary Inverse-Compton scatterings and extragalactic DM distributions. We find that these effects can significantly improve the DM constraints, up to an order of magnitude in some cases. This highlights the importance of considering secondary effects in detail in LHAASO-WCDA and SWGO in the future.

We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps, the primary goal of this paper is to provide the methodology to carry out the end-to-end study of their effect for a space-borne CMB polarization experiment, up to the cosmological results in the form of a bias $\delta r$ on the tensor-to-scalar ratio $r$. LiteBIRD is dedicated to target the measurement of CMB primordial $B$ modes by reaching a sensitivity of $\sigma \left( r \right) \leq 10^{-3}$ assuming $r=0$. As a demonstration of our framework, we derive the relationship between the knowledge of the beam far side-lobes and the tentatively allocated error budget under given assumptions on design, simulation and component separation method. We assume no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that $\delta r$ is mostly due to the integrated fractional power difference between the estimated beams and the true beams in the far side-lobes region, with little dependence on the actual shape of the beams, for low enough $\delta r$. Under our set of assumptions, in particular considering the specific foreground cleaning method we used, we find that the integrated fractional power in the far side-lobes should be known at a level as tight as $\sim 10^{-4}$, to achieve the required limit on the bias $\delta r < 1.9 \times 10^{-5}$. The framework and tools developed for this study can be easily adapted to provide requirements under different design, data analysis frameworks and for other future space-borne experiments beyond LiteBIRD.

Understanding the co-evolution of complex life with Earth's geology is an enduring challenge. The rock record evidences remarkable correlations between changes in biology and the wider Earth system, yet cause and effect remain unclear. Here, we link the evolutionary history of eukaryotes with the rise and fall of carbonate rock fraction within continental crust - a key variable in controlling the efficiency of carbon drawdown during weathering, solid Earth degassing rates, and ultimately nutrient supply to life. We use geospatial database analyses to demonstrate a strongly non-linear growth and then collapse in Earth's continental crust carbonate reservoir. Biomineralisers reshaped Earth's surface in their image; armouring continental margins with carbonate platforms, such that the continental carbonate reservoir increased in size by 5-fold in under 100 Myr after the Cambrian Radiation of animal life. This Paleozoic carbonate revolution represents among the most dramatic crustal evolutionary events in Earth's history. The Permo-Triassic extinction event coupled to the rise of open ocean calcifiers initiated a steady decline in continental crustal carbonate content; one that still continues today, which unabated would produce Precambrian-style crustal carbonate distributions in around 500-1000 Myr. Our results demonstrate strongly non-linear crustal evolution after the rise of the complex Phanerozoic biosphere. This outcome suggests that complex life may generate unique biogeochemical trajectories on otherwise geologically similar worlds, posing a new challenge in the hunt for life beyond Earth.

Cosmological gravitational particle production (CGPP) is the creation of particles in an expanding universe due solely to their gravitational interaction. These particles can play an important role in the cosmic history through their connection to various cosmological relics including dark matter, gravitational wave radiation, dark radiation, and the baryon asymmetry. This review explains the phenomenon of CGPP as a consequence of quantum fields in a time-dependent background, catalogs known results for the spectra and cosmological abundance of gravitationally produced particles of various spins, and explores the phenomenological consequences and observational signatures of CGPP.

This chapter investigates the complexities surrounding the iron biogeochemical cycle from the Archean to present, with a focus on assessing the balance between iron sources and sinks during long periods of Earth's history with relatively invariable redox conditions, when steady state can be safely assumed. Currently, the residence time of iron in the ocean may be as short as approximately 5 years. The input flux of iron is highly sensitive to redox cycling in sediments, while its removal primarily occurs through dispersed processes of oxidation and precipitation. In the Archean, we find a significant imbalance between continental and hydrothermal inputs, which collectively contribute between 61,500 to 263,000 Gg/yr of dissolved iron to the oceans, and the most obvious sinks such as iron formations (IFs), which sequester up to ~43,000 Gg/yr of iron. A possible solution to this imbalance involves the dispersed abiotic precipitation and removal of iron as silicates, sulfides, and carbonates in marine basins. Additionally, we calculate the residence time of dissolved iron in the Archean oceans to be between 6 kyr and 3 Myr, which is significantly longer than the ocean mixing timescale. Our estimates indicate that under the anoxic Archean atmosphere, the iron cycle was more protracted than today, and the isotopic compositions and concentrations of dissolved iron were likely more uniform. Distinct water bodies were likely confined to limited areas or specific, dynamic systems with intense iron turnover, such as regions where deep-sea upwelling currents brought hydrothermal iron to photic zones rich in biotic or abiotic oxidants.

Stellar radio emission can measure a star's magnetic field strength and structure, plasma density and dynamics, and the stellar wind pressure impinging on exoplanet atmospheres. However, properly interpreting the radio data often requires temporal baselines that cover the rotation of the stars, orbits of their planets and any longer-term stellar activity cycles. Here we present our monitoring campaign on the young, active M dwarf AU Microscopii with the Australia Telescope Compact Array between 1.1 and 3.1 GHz. With over 250 hours of observations, these data represent the longest radio monitoring campaign on a single main-sequence star to date. We find that AU Mic produces a wide variety of radio emission, for which we introduce a phenomenological classification scheme predicated on the polarisation fraction and time-frequency structure of the emission. The radio emission detected on AU Mic can be broadly categorised into five distinct types of bursts, and broadband quiescent emission. The radio bursts are highly circularly polarised and periodic with the rotation period of the star, implying that the emission is beamed. It is therefore most likely produced by the electron cyclotron maser instability. We present a model to show that the observed emission can be explained with auroral rings on the magnetic poles. The total intensity of the broadband emission is stochastic, but we show that its circular polarisation fraction is also periodic with the rotation of the star. We present a qualitative model to describe the periodicity in the polarisation fraction of the broadband emission using gyromagnetic emission, and infer a magnetic obliquity of at least 20 degrees from the observed variation in polarisation fraction. Finally, we show that the radio emission might be evolving on long timescales, hinting at a potential stellar magnetic activity cycle.

In some quantum gravity (QG) theories, Lorentz symmetry may be broken above the Planck scale. The Lorentz invariance violation (LIV) may induce observable effects at low energies and be detected at high energy astrophysical measurements. The Large High Altitude Air Shower Observatory(LHAASO) has detected the onset, rise, and decay phases of the afterglow of GRB 221009A, covering a wide energy range of photons approximately from $0.2$ to $18$ TeV. This observation provides an excellent opportunity to study the Lorentz invariance violation effect. In this study, we simultaneously utilize the data from the KM2A and WCDA detectors of LHAASO, and apply two event by event methods, namely the pair view method and maximum likelihood method, to investigate LIV. We obtain stringent constraints on the QG energy scale. For instance, through the maximum likelihood method, we determine the 95$\%$ confidence level lower limits to be $E_{QG,1} > 14.7 (6.5)\times 10^{19}$GeV for the subluminal (superluminal) scenario of $n = 1$, and $E_{QG,2} > 12.0 (7.2)\times 10^{11}$GeV for the subluminal (superluminal) scenario of $n = 2$. We find that the rapid rise and slow decay behaviors of the afterglow can impose strong constraints on the subluminal scenario, while the constraints are weaker for the superluminal scenario.

We present results of explicit asymptotic approximations applied to neutrino--electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and timestepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each timestep can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

Modern astronomical surveys have multiple competing scientific goals. Optimizing the observation schedule for these goals presents significant computational and theoretical challenges, and state-of-the-art methods rely on expensive human inspection of simulated telescope schedules. Automated methods, such as reinforcement learning, have recently been explored to accelerate scheduling. However, there do not yet exist benchmark data sets or user-friendly software frameworks for testing and comparing these methods. We present DeepSurveySim -- a high-fidelity and flexible simulation tool for use in telescope scheduling. DeepSurveySim provides methods for tracking and approximating sky conditions for a set of observations from a user-supplied telescope configuration. We envision this tool being used to produce benchmark data sets and for evaluating the efficacy of ground-based telescope scheduling algorithms, particularly for machine learning algorithms that would suffer in efficacy if limited to real data for training.We introduce three example survey configurations and related code implementations as benchmark problems that can be simulated with DeepSurveySim.

We present high-spectral resolution M-band spectra from iSHELL on NASA's Infrared Telescope Facility (IRTF) along the line of sight to the debris disk host star HD 32297. We also present a Gemini Planet Imager (GPI) H-band polarimetric image of the HD 131488 debris disk. We search for fundamental CO absorption lines in the iSHELL spectra of HD 32297 but do not detect any. We place an upper limit on the CO column density of $\sim$6$\times10^{15}$ cm$^{-2}$. By combining the column density upper limit, the CO mass measured with ALMA, and the geometrical properties of the disk, we estimate the scale height of the CO to be $\lesssim$ 2 au across the radial extent of the disk ($\sim$80-120 au). We use the same method to estimate the CO scale height of three other edge-on, CO-rich debris disks that all have CO observed in absorption with HST as well as in emission with ALMA: $\beta$ Pictoris, HD 110058, and HD 131488. We compare our estimated CO scale heights of these four systems to the millimeter dust scale heights and find that, under the assumption of hydrostatic equilibrium, there is a potential correlation between the CO and millimeter dust scale heights. There are multiple factors that affect the gas vertical structure such as turbulence, photodissociation with weak vertical mixing, as well as where the gas originates. One possible explanation for the potential correlation could be that the gas and dust are of a similar secondary origin in these four systems.

One of the most exciting targets of current and future gravitational-wave observations is the angular power spectrum of the astrophysical GW background. This cumulative signal encodes information about the large-scale structure of the Universe, as well as the formation and evolution of compact binaries throughout cosmic time. However, the finite rate of compact binary mergers gives rise to temporal shot noise, which introduces a significant bias in measurements of the angular power spectrum if not explicitly accounted for. Previous work showed that this bias can be removed by cross-correlating GW sky maps constructed from different observing times. However, this work considered an idealised measurement scenario, ignoring detector specifics and in particular noise contributions. Here we extend this temporal cross-correlation method to account for these difficulties, allowing us to implement the first unbiased anisotropic search pipeline for LIGO-Virgo-KAGRA data. In doing so, we show that the existing pipeline is biased even in the absence of shot noise, due to previously neglected sub-leading contributions to the noise covariance. We apply our pipeline to mock LIGO data, and find that our improved analysis will be crucial for stochastic searches from the current observing run (O4) onwards.

The Gaia mission is delivering a large number of astrometric orbits for binary stars. By combining these with spectroscopic orbits for systems with two observable spectra (SB2), it is possible to derive the masses of both components. However, to get masses with a good accuracy requires accurate spectroscopic orbits, which is the major aim of the present paper. A subsidiary aim is to discover SB2 systems hiding among known SB1, and even though this search may often prove unsuccessful, the acquired radial velocities may be used anyway to improve the existing spectroscopic orbits. New radial velocities for 58 binary systems from the Ninth Catalogue of Spectroscopic Binary Orbits (SB9), obtained using the high-resolution HERMES spectrograph installed on the 1.2 m Mercator telescope, were used to possibly identify hitherto undetected SB2 systems. For SB1 systems with inaccurate orbits, these new radial-velocity measurements were used to improve the orbital accuracy. This paper provides 51 orbits (41 SB1 and 10 SB2) that have been improved with respect to the solution listed in the SB9 catalogue, out of the 58 SB9 orbits studied, which belong to 56 stellar systems. Among them, there are five triple and four quadruple systems. Despite the high resolution of HERMES, HIP 115142 A is the only system which we detected as a new SB2 system. The B component of the visual binary HIP 92726 has now been found to be a spectroscopic system as well, which makes HIP 92726 a newly discovered quadruple system (SB1+SB1). The high HERMES resolution allowed us moreover to better isolate the signature of the secondary component of HIP 12390, HIP 73182 and HIP 111170. More accurate masses have thus been derived for them. Among the 30 SB also present in Gaia Data Release 3 (DR3) and with periods shorter than the Gaia DR3 time span (1000 d), only 5 were flagged as binaries by DR3.

In the hunt for Earth-like exoplanets it is crucial to have reliable host star parameters, as they have a direct impact on the accuracy and precision of the inferred parameters for any discovered exoplanet. For stars with masses between 0.35 and 0.5 ${\rm M_{\odot}}$ an unexplained radius inflation is observed relative to typical stellar models. However, for fully convective objects with a mass below 0.35 ${\rm M_{\odot}}$ it is not known whether this radius inflation is present as there are fewer objects with accurate measurements in this regime. Low-mass eclipsing binaries present a unique opportunity to determine empirical masses and radii for these low-mass stars. Here we report on such a star, EBLM J2114-39\,B. We have used HARPS and FEROS radial-velocities and TESS photometry to perform a joint fit of the data, and produce one of the most precise estimates of a very low mass star's parameters. Using a precise and accurate radius for the primary star using Gaia DR3 data, we determine J2114-39 to be a $M_1 = 0.987 \pm 0.059$ ${\rm M_{\odot}}$ primary star hosting a fully convective secondary with mass $M_2~=~0.0986~\pm 0.0038~\,\mathrm{M_{\odot}}$, which lies in a poorly populated region of parameter space. With a radius $R_2 =~0.1275~\pm0.0020~\,\mathrm{R_{\odot}}$, similar to TRAPPIST-1, we see no significant evidence of radius inflation in this system when compared to stellar evolution models. We speculate that stellar models in the regime where radius inflation is observed might be affected by how convective overshooting is treated.

We present findings of 3D filamentary structures in the Smith Cloud, a high-velocity cloud (HVC) located at $l=38^{\circ}$, $b=-13^{\circ}$. We use data from the Galactic Arecibo L-Band Feed Array \ion{H}{i} (GALFA-\ion{H}{i}) along with our new filament detection algorithm, \texttt{fil3d}, to characterize these structures. In this paper, we also discuss how different input parameters affect the output of \texttt{fil3d}. We study filaments in the local ISM and compare them to those found in the Smith Cloud. Based on thermal linewidth estimations we find supporting evidence that the Smith Cloud filaments are part of its warm neutral medium. We also find a relationship between thermal linewidth and the $v_{LSR}$ of the filaments. We study the plane-of-sky magnetic field as traced by Planck 353 GHz polarized dust emission along the line of sight and find the HI filaments in this region are not aligned with the magnetic field. This is likely related to their location close to dynamic processes in the Galactic Plane and/or the low column density of the filaments relative to emission in the Plane. The results show the HI filaments are found in a wide range of Galactic environments and form through multiple processes.

We present the first Swift/XRT long-term monitoring of 2S 0114+650, a wind-fed supergiant X-ray binary for which both orbital and superorbital periods are known (P_orb~11.6d and P_sup~0.8d). Our campaign, summing up to ~ 79ks, is the most intense and complete sampling of the X-ray light curve of this source with a sensitive pointed X-ray instrument, and covers 17 orbital, and 6 superorbital cycles. The combination of flexibility, sensitivity, and soft X-ray coverage of XRT allowed us to confirm previously reported spectral changes along the orbital cycle of the source and unveil the variability in its spectral parameters as a function of the superorbital phase. For completeness, we also report on a similar analysis carried out by exploiting XRT archival data on three additional wind-fed supergiant X-ray binaries IGR J16418-4532, IGR J16479-4514, and IGR J16493-4348. For these sources, the archival data provided coverage along several superorbital cycles but our analysis could not reveal any significant spectral variability.

Recent studies have shown that the large-scale gas dynamics of protoplanetary disks (PPDs) are controlled by non-ideal magneto-hydrodynamics (MHD), but how this influences dust dynamics is not fully understood. To this end, we investigate the stability of dusty, magnetized disks subject to the Hall effect, which applies to planet-forming regions of PPDs. We find a novel Background Drift Hall Instability (BDHI) that may facilitate planetesimal formation in Hall-effected disk regions. Through a combination of linear analysis and nonlinear simulations, we demonstrate the viability and characteristics of BDHI. We find it can potentially dominate over the classical streaming instability (SI) and standard MHD instabilities at low dust-to-gas ratios and weak magnetic fields. We also identify magnetized versions of the classic SI, but these are usually subdominant. We highlight the complex interplay between magnetic fields and dust-gas dynamics in PPDs, underscoring the need to consider non-ideal MHD like the Hall effect in the broader narrative of planet formation.

The elusive massive black hole (MBH) seeds stand to be revealed by the Laser Space Antenna Interferometer through mergers. As an aftermath of galaxy mergers, MBH coalescence is a vastly multi-scale process connected to galaxy formation. We introduce the "Massive black hole Assembly in Galaxies Informed by Cosmological Simulations" (MAGICS) suite, with galaxy/MBH properties and orbits recovered from large-volume cosmological simulation ASTRID. The simulations include subgrid star formation, supernovae feedback, and MBH accretion/feedback. In this first suite, we extract fifteen representative galaxy mergers with seed MBHs to examine their dynamics at an improved mass and spatial resolution (by $\sim2000$ and $\sim20$) and follow MBH orbits down to $\sim10\,\text{pc}$. We find that the seed MBH energy loss and orbital decay are largely governed by global torques induced by the galaxy merger process on scales resolvable by cosmological simulations. Specifically, pairs sink quickly if their orbits shrink rapidly below $1\,\text{kpc}$ during the first $\sim200\,\text{Myr}$ of pairing due to effective energy loss in major galaxy mergers, whereas MBHs gaining energy in minor galaxy mergers with head-on collisions are likely to stall. High initial eccentricities ($e_\text{init}>0.5$) and high stellar densities at kpc scales ($\rho_\text{star}>0.05\,M_\odot/\text{pc}^3$) also lead to most efficient decays. $\sim50\%$ high-redshift seed MBH pairs experience consecutive galaxy mergers and are more likely to stall at $\sim1\,\text{kpc}$. For a subset of systems, we carry out N-Body re-simulations until binary formation and find that some stalled systems merge at high-z when embedded in sufficient nuclear star clusters.

We study in detail the inner 600 pc of the Seyfert 2 galaxy ESO138-G001 by means of the SOAR Integral Field Spectrograph (SIFS) attached to the SOAR telescope. This source is known for displaying a very rich coronal line spectrum and a blob of high-excitation emission ~3" SE of the active galactic nucleus (AGN). The nature of this emission has not been fully understood yet. The excellent spatial and spectral resolution of SIFS allows us to confirm that the bulk of the coronal line forest emission region is very compact, of ~0.8" in diameter, centred on the AGN and most likely powered by radiation from the AGN. In addition, evidence of a nuclear outflow, restricted to the inner 1" centred at the nucleus is found based on the detection of broad components in the most important emission lines. The gas in the inner few tens of parsecs filters out the AGN continuum so that the NLR is basically illuminated by a modified SED. This scenario is confirmed by means of photoionisation models that reproduce the most important lines detected in the SIFS field-of-view. From the modelling, we also found that the black hole mass M_BH of the AGN is about 10^5.5 solar mass, in agreement with previous X-ray observations. The spectrum of the SE blob is dominated by emission lines of low- to mid-ionisation, with no hints of coronal lines. Our results show that it represents gas in the ionisation cone that is photoionised by the filtered central AGN continuum.

Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the James Webb Space Telescope was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present MIRI imaging mosaics at 5.6 and 21 microns that cover a significant portion of one of M33's spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep ALMA CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new DOLPHOT MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further shows YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass.

We present a catalog of 6266 solar flares detected by the X-Ray Solar Monitor onboard the Chandrayaan-2 lunar orbiter between 1.55 and 12.4 keV (1 and 8 \AA) from 2019 September 12 to 2022 November 4, including 1469 type A flares. The catalog represents the first large sample, including both type A, hot thermal flares, and type B, impulsive flares, with a sub-A class sensitive instrument. We also detect 213 sub-A and 1330 A class flares. Individual flares are fit with an exponentially-modified Gaussian function and multi-flare groups are decomposed into individual flares. We validate our findings with flare catalogs made using visual inspection as well as automatic pipelines on Geostationary Operational Environmental Satellite and Solar Dynamics Observatory data. We find a clear bimodality in the ratio of the width to decay time between type A and B flares. We infer a power-law index of $\alpha_F = 1.92 \pm 0.09$ for the background-subtracted peak flux distribution of XSM flares, which is consistent with the value $\sim 2$ reported in the literature. We also infer $\alpha_F = 1.90 \pm 0.09$ for type B, and $\alpha_F = 1.94 \pm 0.08$ for type A flares, which has previously not been reported in the literature. These comparable values hint at a similarity in their generative processes.

Accretion around black holes is very often characterized by distinctive X-ray reflection features (mostly, iron inner-shell transitions), which arise due to the primary radiation being reprocessed by a dense and relatively colder medium, such as an accretion disk. Most reflection modeling assume that emission stops at the inner-most stable circular orbit (ISCO), and that for smaller radii - in the plunging region - the density drops and the accretion flow is far too ionized for efficient line production. We investigate the spectral features of the reflection in the plunging regions of optically-thick and geometrically-thin accretion disks around black holes. We show that for cases in which the density profile is considered constant (as expected in highly magnetized flows), or in cases in which the disk density is high enough such that the ionization still allows line formation within the ISCO, there is a significant modification of the observed reflected spectrum. Consistent with previous studies, we found that the impact of the radiation reprocessed in the plunging region is stronger the lower the black hole spin, when the plunging region subtends a larger area. Likewise, as for the case of standard reflection modeling, the relativistic broadening of the iron line is more pronounced at low inclination, whereas the blueshift and relativistic beaming effect is dominant at high inclination. We also tested the effects of various prescriptions of the stress at the ISCO radius on the reflection spectrum, and found that several of these cases appear to show line profiles distinct enough to be distinguishable with reasonably good quality observational data.

The JWST is allowing new measurements of gas-phase metallicities in galaxies between cosmic noon and cosmic dawn through observations of multiple rest-frame optical/ultraviolet [OIII], [OII], and hydrogen Balmer series lines. The most robust approach to such measurements uses luminosity ratios between the excited auroral transition, [OIII] 4364A, and the lower [OIII] 5008A/4960A lines to determine the gas temperature. The ratio of the luminosities in the latter transitions to those in hydrogen Balmer series lines then yield relatively clean metallicity estimates. In the absence of detection of the, often faint, [OIII] auroral line, however, the ratios of various [OIII], [OII], [NII], and Balmer lines are instead used to determine metallicities. Here we present a refined approach for extracting metallicities from these ``strong line diagnostics''. Our method exploits empirical correlations between the temperature of OIII/OII regions and gas-phase metallicity, which lie close to theoretical expectations from thermal equilibrium calculations. We then show, from first principles, how to extract metallicities from traditional strong line diagnostics, R2, R3, R23, O3O2, and N2O2. We show that these ratios depend also on ionization correction factors, but that these can be determined self-consistently along with the metallicities. We quantify the success of our method using metallicities derived from galaxies with auroral line determinations and show that it generally works better than previous empirical approaches. The scatter in the observed line ratios and redshift evolution are largely explained by O3O2 variations. We provide publicly available routines for extracting metallicities from strong line diagnostics using our methodology.

This paper presents a method for finding ram-pressure stripped (RPS) galaxy candidates by performing a morphological analysis of galaxy images obtained from the Legacy survey. We consider a sample of about 600 galaxies located in different environments such as groups and clusters, tidally interacting pairs and the field. The sample includes 160 RPS previously classified in the literature into classes from J1 to J5, based on the increasing level of disturbances. Our morphological analysis was done using the {\sc astromorphlib} software followed by the inspection of diagnostic diagrams involving combinations of different parameters like the asymmetry ($A$), concentration ($C$), S\'ersic index ($n$), and bulge strength parameters $F(G,\,M_{20})$. We found that some of those diagrams display a distinct region in which galaxies classified as J3, J4 and J5 decouples from isolated galaxies. We call this region as the morphological transition zone and we also found that tidally interacting galaxies in pairs are predominant within this zone. Nevertheless, after visually inspecting the objects in the morphological transition zone to discard obvious contaminants, we ended up with 33 bonafide new RPS candidates in the studied nearby groups and clusters (Hydra, Fornax, and CLoGS sample), of which one-third show clear evidence of unwinding arms. Future works may potentially further increase significantly the samples of known RPS using such method.

The AMS-02 experiment has observed new properties of primary cosmic rays (CRs) categorized into two groups: He-C-O-Fe and Ne-Mg-Si-S, which are independent of CR propagation. In this study, we investigate the unexpected properties of these nuclei using a spatial propagation model. All nuclei spectra are accurately reproduced and separated into primary and secondary contributions. Our findings include: 1. Primary CR spectra are identical. 2. Our calculations align with AMS-02 results for primary-dominated nuclei within a 10\% difference, but show significant discrepancies for the secondary-dominated nuclei. 3. The primary element abundance is presented for the first time. We anticipate that the DAMPE and future HERD experiments will provide observations of nuclei spectra above TeV energy.

The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu Goldstone boson of global $\mathrm{U}(1)_{\rm PQ}$ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals.

We present a novel approach, $\textit{Metric pErTuRbations wIth speCtral methodS}$ (METRICS), to calculate the gravitational metric perturbations and the quasinormal-mode frequencies of rotating black holes of any spin without decoupling the linearized field equations. We demonstrate the method by applying it to perturbations of Kerr black holes of any spin, simultaneously solving all ten linearized Einstein equations in the Regge-Wheeler gauge through purely algebraic methods and computing the fundamental (co-rotating) quadrupole mode frequency at various spins. We moreover show that the METRICS approach is accurate and precise, yielding (i) quasinormal mode frequencies that agree with Leaver's, continuous-fraction solution with a relative fractional error smaller than $10^{-5}$ for all dimensionless spins below up to 0.95, and (ii) metric perturbations that lead to Teukolsky functions that also agree with Leaver's solution with mismatches below $1\%$ for all spins below 0.9. By not requiring the decoupling or the angular separation of the linearized field equations, the METRICS approach has the potential to be straightforwardly adapted for the computation of the quasinormal-mode frequencies of rotating black holes of any spin beyond general relativity or in the presence of matter.

Binary neutron stars mergers widely accepted as potential progenitors of short gamma-ray bursts. After the remnant of the merger has collapsed to a black hole, a jet is powered and may breakout from the the matter expelled during the collision and the subsequent wind emission. The interaction of the jet with the ejecta may affect its dynamics and the resulting electromagnetic counterparts. We here examine how an inhomogeneous and anisotropic distribution of ejecta affects such dynamics, dictating the properties of the jet-ejecta cocoon and of the afterglow radiated by the jet upon deceleration. More specifically, we carry out general-relativistic hydrodynamical simulations of relativistic jets launched within a variety of geometrically inhomogeneous and anisotropic distributions of ejected matter. We find that different anisotropies impact the variance of the afterglow light-curves as a function of the jet luminosity and ejected mass. A considerable amount of the jet energy is deposited in the cocoon through the jet-ejecta interaction with a small but important dependence on the properties of the ejecta. Furthermore, all configurations show a two-component behaviour for the polar structure of the jet, with a narrow core at large energies and Lorentz factors and a shallow segment at high latitudes from the jet axis. Hence, afterglows measured on off-axis lines of sight could be used to deduce the properties of the ejected matter, but also that the latter need to be properly accounted for when modelling the afterglow signal and the jet-launching mechanisms.

The spectra of gravitational waves from black hole evaporation generically peak at frequencies of order the Hawking temperature, making this signal ultra-high frequency for primordial black holes evaporating in the early universe. This motivates us to consider small black holes in theories with large extra dimensions, for which the peak frequency can be lowered substantially, since the true bulk Planck scale $M_*$ can be much smaller than the effective $M_{\rm Pl}$. We study the emission of brane-localized gravitons during the Hawking evaporation of ultra-light primordial black holes in the context of theories with large extra dimensions, with the ultimate goal of computing the contribution to the stochastic gravitational wave background. To accurately model black hole evolution, we compute greybody factors for all particle species emitted on the brane and in the bulk, presuming the majority of emission proceeds during the Schwarzschild phase. We then compute the power spectrum and present day spectral density parameter for brane-localized gravitons contributing to a gravitational wave signal. We find that for an optimal choice of parameters, the peak frequency plateaus in the sub-MHz regime, within range of planned high-frequency gravitational wave detectors, making this scenario a target for detection once their sensitivity exceeds $\Delta N_{\rm eff}$ bounds.

The automatic identification of planetary feature names in astronomy publications presents numerous challenges. These features include craters, defined as roughly circular depressions resulting from impact or volcanic activity; dorsas, which are elongate raised structures or wrinkle ridges; and lacus, small irregular patches of dark, smooth material on the Moon, referred to as "lake" (Planetary Names Working Group, n.d.). Many feature names overlap with places or people's names that they are named after, for example, Syria, Tempe, Einstein, and Sagan, to name a few (U.S. Geological Survey, n.d.). Some feature names have been used in many contexts, for instance, Apollo, which can refer to mission, program, sample, astronaut, seismic, seismometers, core, era, data, collection, instrument, and station, in addition to the crater on the Moon. Some feature names can appear in the text as adjectives, like the lunar craters Black, Green, and White. Some feature names in other contexts serve as directions, like craters West and South on the Moon. Additionally, some features share identical names across different celestial bodies, requiring disambiguation, such as the Adams crater, which exists on both the Moon and Mars. We present a multi-step pipeline combining rule-based filtering, statistical relevance analysis, part-of-speech (POS) tagging, named entity recognition (NER) model, hybrid keyword harvesting, knowledge graph (KG) matching, and inference with a locally installed large language model (LLM) to reliably identify planetary names despite these challenges. When evaluated on a dataset of astronomy papers from the Astrophysics Data System (ADS), this methodology achieves an F1-score over 0.97 in disambiguating planetary feature names.

To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulations, we confirm that the relative fundamental quasinormal mode amplitudes of precessing systems are related to those of non-precessing systems by a simple rotation, and that additional structure in the spectrum is connected to the system's kick velocity and other asymmetries in the orbital dynamics. We find that the ringdowns of precessing systems need not be dominated by the ${(\ell,m)=(2,\pm 2)}$ quasinormal modes. These results build upon previous works on waveform modeling, and are consistent with a recent ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521.

The heliosphere is permeated with highly structured solar wind originating from the sun. One of the primary science objectives of Parker Solar Probe (PSP) is to determine the structures and dynamics of the plasma and magnetic fields at the sources of the solar wind. However, establishing the connection between in situ measurements and structures and dynamics in the solar atmosphere is challenging: most of the magnetic footpoint mapping techniques have significant uncertainties in the source localization of a plasma parcel observed in situ, and the PSP plasma measurements suffer from a limited field of view. Therefore it lacks a universal tool to self-contextualize the in situ measurements. Here we develop a novel time series visualization method named Jensen-Shannon Scalogram. Utilizing this method, by analyzing the magnetic magnitude data from both PSP and Ulysses, we successfully identify in situ remnants of solar atmospheric and magnetic structures spanning more than seven orders of magnitude, from years to seconds, including polar and mid-latitude coronal holes, as well as structures compatible with super-granulation, ``jetlets'' and very small scale flaring activity. Furthermore, computer simulations of Alfv\'enic turbulence support key features of the observed magnetic magnitude distribution. Building upon these discoveries, the Jensen-Shannon Scalogram therefore not only enables us to reveal the fractal fine structures in the solar wind time series from both PSP and decades-old data archive, but will also serve as a general-purpose data visualization method applicable to all time series.

The present paper is a modest attempt to initiate the research program outlined in this abstract. We propose that general relativity and relativistic MOND (RelMOND) are analogues of the broken electroweak symmetry. That is, $SU(2)_R \times U(1)_{YDEM} \rightarrow U(1)_{DEM}$ (DEM stands for dark electromagnetism), and GR is assumed to arise from the broken $SU(2)_R$ symmetry, and is analogous to the weak force. RelMOND is identified with dark electromagnetism $U(1)_{DEM}$, which is the remaining unbroken symmetry after spontaneous symmetry breaking of the darkelectro-grav sector $SU(2)_R \times U(1)_{YDEM}$. This sector, as well as the electroweak sector, arise from the breaking of an $E_8 \times E_8$ symmetry, in a recently proposed model of unification of the standard model with pre-gravitation, this latter being an $SU(2)_R$ gauge theory. The source charge for the dark electromagnetic force is square-root of mass, motivated by the experimental fact that the square-roots of the masses of the electron, up quark, and down quark, are in the ratio 1:2:3, which is a flip of their electric charge ratios 3:2:1 The introduction of the dark electromagnetic force helps understand the weird mass ratios of the second and third generation of charged fermions. We also note that in the deep MOND regime, acceleration is proportional to square-root of mass, which motivates us to propose the relativistic $U(1)_{DEM}$ gauge symmetry as the origin of MOND. We explain why the dark electromagnetic force falls inversely with distance, as in MOND, and not as the inverse square of distance. We conclude that dark electromagnetism is a good mimicker of cold dark matter, and the two are essentially indistinguishable in those cosmological situations where CDM is successful in explaining observations, such as CMB anisotropies, and gravitational lensing.

We suggest the procedure of building the maps of noctilucent clouds (NLC) zonal and meridional velocity, mean altitude and particle size based on three-color photometry by identical all-sky RGB-cameras separated by 115 km in a close-meridional direction. The procedure is applied to the bright NLC event on July 3, 2023. The altitude is measured by precise triangulation technique, and effective particle radius is estimated by comparison of each definite NLC fragment intensity and colors at different scattering angles as it is registered from different observation sites. The results are compared with existing photometric methods for average altitude and particle size measurements. A significant difference in evening and morning twilight NLC parameters is found, which is discussed in comparison with existing analysis of diurnal NLC variations.

A small component of dark matter (DM) that is strongly interacting with the standard model sector is consistent with various experimental observations. Despite the small abundance, strongly-interacting DM can lead to pronounced signals in DM direct detection experiments. We study Belle II sensitivity on strongly-interacting DM that has a MeV-GeV mass and couples with electrons. By taking into account the substantial interactions between DM and electrons within detectors, we compute the ceiling of the mono-photon signature at Belle II, beyond which the mono-photon channel loses its sensitivity, and visible ECL clusters due to DM scatterings assume significance. We study two ECL signatures for strongly-interacting DM: the mono-cluster and the di-cluster channels. To carry out detailed calculations and to compare with other constraints, we consider DM models with light mediators, as they naturally lead to sizable interaction cross sections. We compute exclusion regions for the di-cluster, mono-cluster, and mono-photon channels. We find that Belle II (with currently accumulated data of 362 fb$^{-1}$) can rule out a significant portion of the parameter space above the ceilings of the constraints from various DM direct detection and neutrino experiments, for the vector mediator case with mass $\gtrsim$10 MeV. Belle II also offers superior constraints on new light particles compared to PBH for the scalar mediator with mass $\gtrsim$10 MeV.