Papers for Friday, Apr 11 2025

The International Astronomical Union definitions for Planet and Dwarf Planet both require that a body has sufficient mass to overcome rigid body forces and self gravitate into a nearly round shape. However, quantitative standards for determining when a body is sufficiently round have been lacking. Previously published triaxial ellipsoid solutions for asteroids, satellites, and Dwarf Planets in the radius range 135 to 800 km are examined to identify a minimum mass above which the entire population, regardless of composition, is round. From this data, the minimum mass to meet the roundness criterion is 5.0 x 10E20 kg. The triaxial shape data suggests three radius ranges: (1) bodies with a radius less than 160 km are non-spheroidal, (2) bodies with a radius in the range 160 to 450 km are transitional in shape or nearly round, (3) bodies with a radius greater than 450 km are spheroidal. Bodies orbiting the Sun with a mass greater than 5.0 x 10E20 kg are Planets or Dwarf Planets. Arguments are presented for including the 16 spheroidal moons of the Solar System as a third dynamical class that can be identified as Satellite Planets. Definitions are proposed that expand upon the taxonomy started in 2006 with the International Astronomomical Union Planet and Dwarf Planet classes.

Quasi-Periodic Eruptions (QPEs) are recurring bursts of soft X-ray emission from supermassive black holes (SMBHs), which a growing class of models explains via extreme mass-ratio inspirals (EMRIs). QPEs exhibit blackbody-like emission with significant temperature evolution, but the minimal information content of their almost pure-thermal spectra has limited physical constraints. Here we study the recently discovered QPEs in ZTF19acnskyy (``Ansky''), which show absorption-like features evolving dramatically within eruptions and correlating strongly with continuum temperature and luminosity, further probing the conditions underlying the emission surface. The absorption features are well-described by dense ionized plasma of column density $N_{\rm H}\gtrsim 10^{21}$ cm$^{-2}$, blueshift $0.06\lesssim v/c \lesssim 0.4$, and either collisional or photoionization equilibrium. With high-resolution spectra, we also detect ionized blueshifted emission lines suggesting a nitrogen over-abundance of $21.7^{+18.5}_{-11.0}\times$ solar. We interpret our results with orbiter-disk collisions in an EMRI system, in which each impact drives a shock that locally heats the disk and expels X-ray emitting debris undergoing radiation pressure-driven homologous expansion. We explore an analytical toy model that links the rapid change in absorption lines to the evolution of the ionization parameter and the photosphere radius, and suggest that $\sim 10^{-3}M_\odot$ ejected per eruption with expansion velocities up to $v_{\rm max}\sim 0.15c$, can reproduce the absorption features. With these assumptions, we show a P Cygni profile in a spherical expansion geometry qualitatively matches the observed line profiles. Our work takes a first step towards extending existing physical models for QPEs to address their implications for spectral line formation.

Quasi-periodic eruptions (QPEs) are rapid, recurring X-ray bursts from supermassive black holes, believed to result from interactions between accretion disks and surrounding matter. The galaxy SDSS1335+0728, previously stable for two decades, exhibited an increase in optical brightness in December 2019, followed by persistent Active Galactic Nucleus (AGN)-like variability for 5 years, suggesting the activation of a $\sim$10$^6\;M_\odot$ black hole. From February 2024, X-ray emission has been detected, revealing extreme $\sim$4.5-day QPEs with the highest fluxes and amplitudes, longest time scales, largest integrated energies, and a $\sim$25-day super-period. Low-significance UV variations are reported for the first time in a QPE host, likely related to the long timescales and large radii from which the emission originates. This discovery broadens the possible formation channels for QPEs, suggesting they are not linked solely to tidal disruption events but more generally to newly formed accretion flows, which we are witnessing in real time in a turn-on AGN candidate.

The X-ray polarization of Sco X-1 was measured for the first time with very high significance by the Imaging X-ray Polarimeter Explorer (IXPE). A polarization degree of 1.0% $\pm$ 0.2% at a PA of 8° $\pm$ 6° at 90% confidence level (CL) in the 2-8 keV energy band is obtained, while the source was in its soft state with short flaring periods. The source state was determined by a strictly simultaneous X-ray observation campaign jointly with IXPE, which involved NICER, NuSTAR, and Insight-HXMT, allowing for a broad-band spectrum characterization and study of quasiperiodic oscillations. The spectropolarimetric analysis yielded a polarization of <3.2% for the accretion disk and a polarization of 1.3% $\pm$ 0.4% for the hard Comptonized component. A constraint on the polarization of the reflection component, modeled using relxillNS, is obtained. All the results about the polarization degree match the theoretical expectations, while the polarization angle of 8° +/- 6° at 90% CL shows a rotation of 46° $\pm$ 9° with respect to the measured position angle of the radio jet and previous marginal results by PolarLight. This may suggest a variation in the polarization angle related to the source state, which is linked to the variation of corona geometry as reported by IXPE observations of Z sources, or possibly to relativistic precession.

The IRAS01425+2902 wide binary system was recently reported to have both a young planet and a puzzling geometric arrangement, where the planet and binary both orbit edge-on, but misaligned by 60 deg to the circumprimary disc. This is the youngest transiting planet yet to be detected but its misalignment to the disc is difficult to explain. In this paper we explore the dissolution of an unstable triple system as a potential mechanism to produce this system. We simulate the effects of an ejection interaction in models using a highly inclined, retrograde flyby centred on the primary star of IRAS01425. The escaping star of ~0.35 solar masses inclines both the disc and binary orbits such that they have a relative misalignment of greater than 60 deg, as inferred from observations. The planet orbit also becomes inclined relative to the disc, and our interpretation predicts that the binary should have a highly eccentric orbit (e > 0.5 from our simulations). We additionally demonstrate that despite the high relative misalignment of the disc it is unlikely to be vulnerable to von Zeipel-Kozai-Lidov oscillations.

We measure galaxy sizes from $2 < z < 10$ using COSMOS-Web, the largest-area JWST imaging survey to date, covering $\sim$0.54 deg$^2$. We analyze the rest-frame optical (~5000A) size evolution and its scaling relation with stellar mass ($R_e\propto M_*^\alpha$) for star-forming and quiescent galaxies. For star-forming galaxies, the slope $\alpha$ remains approximately 0.20 at $2 < z < 8$, showing no significant evolution over this redshift range. At higher redshifts, the slopes are $-0.13 \pm 0.15$ and $0.37 \pm 0.36$ for $8 < z < 9$ and $9 < z < 10$, respectively. At fixed galaxy mass, the size evolution for star-forming galaxies follows $R_e \propto (1+z)^{-\beta}$, with $\beta = 1.21 \pm 0.05$. For quiescent galaxies, the slope is steeper $\alpha\sim 0.5$-$0.8$ at $2 < z < 5$, and $\beta=0.81\pm0.26$. We find that the size-mass relation is consistent between UV and optical at $z < 8$ for star-forming galaxies. However, we observe a decrease in the slope from UV to optical at $z > 8$, with a tentative negative slope in the optical at $8 < z < 9$, suggesting a complex interplay between intrinsic galaxy properties and observational effects such as dust attenuation. We discuss the ratio between galaxies' half-light radius, and underlying halos' virial radius, $R_{vir}$, and find the median value of $R_e/R_{vir}=2.7\%$. The star formation rate surface density evolves as $\log\Sigma_\text{SFR} = (0.20\pm0.08)\,z+(-0.65\pm0.51)$, and the $\Sigma_\text{SFR}$-$M_*$ relation remains flat at $2<z<10$. Lastly, we identify a threshold in stellar mass surface density $\log\Sigma_e\sim9.5$-$10\, M_{\odot}/kpc^2$ marking the transition to compact, quenched galaxies from extended, star-forming progenitors. In summary, our findings show that the extensive COSMOS-Web dataset at $z > 3$ provides new insights into galaxy size and related properties in the rest-frame optical.

Direct imaging campaigns executed with James Webb Space Telescope (JWST) will enable the study of the faintest observable exoplanets yet. To assist observers in the JWST proposal process, we present an in-depth exploration of the effects of moderate contrast binary (visual and physical) systems on JWST NIRCam coronagraphic Reference Differential Imaging (RDI) methods. All work in this paper is based on simulation data generated using the Python package PanCAKE, in addition to a variety of custom scripts which we have made publicly available on GitHub. Presenting both contrast curves and more involved 'heatmaps' of sensitivity loss, we present quantifiable measurements for how a binary companion will impact contrast both totally and locally, as a function of magnitude, separation, and position angle. Observers can use results in this work to estimate the impact of a known binary, and in some cases will find that PSF subtraction can still be reliably performed. We have found several scenarios where JWST NIRCam coronagraphic RDI PSF subtraction can be viably performed using a binary reference and make several suggestions. The brightest binary companions analyzed, with a relative brightness of 1e-3, resulted in the worst local sensitivity loss of 3.02 units of magnitude. The faintest binary companions looked at, 1e-6 relative brightness, have almost no effect on local sensitivity. Changing position angle impacts sensitivity loss by 0.5 to 0.3 depending on companion flux. Binary companion separation considerations should be on a case-by-case science goal basis. This work also discusses the trade space for these suggestions in detail.

The Moon is the closest celestial gamma-ray emitting object. Its gamma-ray emission arises from interactions between Galactic cosmic rays (CRs) and the lunar surface. While the lunar GeV gamma-ray spectrum is dominated by a continuum from hadronic decay processes, the MeV emission exhibits both continuum and distinctive spectral lines from nuclear de-excitation and radioactive decay processes. Using Geant4 Monte Carlo particle simulations, we model the lunar gamma-ray spectrum. Our results demonstrate its consistency with Fermi-LAT observations, and predict that next-generation MeV gamma-ray instruments will detect both the lunar MeV continuum and several key spectral line features, notably the $1.779~\mathrm{MeV}$ line from $\mathrm{^{28}Si}$ de-excitation enhanced by the lunar surface composition, the $e^+e^-$ annihilation line, and radioactive decay lines from $\mathrm{^{22}Na}$ ($\tau\approx3.75\,\mathrm{yr}$) and long-lived $\mathrm{^{26}Al}$ ($\tau\approx1\,\mathrm{Myr}$). These gamma-ray lines are sensitive to CRs with energies $\lesssim1\,\mathrm{GeV\,nuc^{-1}}$, offering unique temporal probes of CR activity over different timescales. Observations of the lunar MeV gamma-ray spectrum will therefore open a new window to study the current irradiation of the solar-terrestrial environment by low-energy CRs and its long-term temporal evolution.

We present the HST ACS G800L grism spectroscopy observation of the faint active galactic nuclei (AGN) candidates in the COSMOS field at redshift of 6 selected by the point-source morphology and the photometry drop-off at 8000Å. Among the sample of 7 objects, only one is detected by multiple bands, and has similar shape of spectral energy distribution as the so-called ``little red dots'' JWST selected AGN candidates, but our object is 3 magnitude brighter than the JWST sample. We draw the upper limit of the AGN luminosity function $\Phi=1.1\times 10^{-7}$Mpc$^3$ mag$^{-1}$ for $M_{UV}$=$-21$ at redshift of 6. The rest of the sample shows inconsistent flux density when comparing magnitudes of HST ACS F814W to the Subaru $i$-band and $z$-band magnitudes combined. The HST ACS G800L grism observation shows that this inconsistency cannot be created from an emission line. Therefore, we speculate that these objects are transients with the light curve decay timescale at most 6 years in observed frame.

The LIGO-Virgo-KAGRA catalog has been analyzed with an abundance of different population models due to theoretical uncertainty in the formation of gravitational-wave sources. To expedite model exploration, we introduce an efficient and accurate variational Bayesian approach that learns the population posterior with a normalizing flow and serves as a drop-in replacement for existing samplers. With hardware acceleration, inference takes just seconds for the current set of black-hole mergers and readily scales to larger catalogs. The trained posteriors provide an arbitrary number of independent samples with exact probability densities, unlike established stochastic sampling algorithms that otherwise match with Jensen-Shannon divergences below 0.1 nats in our 14-dimensional parameter space, while requiring up to three orders of magnitude fewer likelihood evaluations and as few as $\mathcal{O}(10^3)$. Provided the posterior support is covered, discrepancies can be addressed with smoothed importance sampling, which quantifies a goodness-of-fit metric for the variational approximation while also estimating the evidence for Bayesian model selection. Neural variational inference thus enables interactive development, analysis, and comparison of population models, making it a useful tool for astrophysical interpretation of current and future gravitational-wave observations.

Be X-ray binaries (BeXRBs) may show strong X-ray and optical variability, and can exhibit some of the brightest outbursts that break through the Eddington limit. Major X-ray outbursts are often accompanied by strong optical flares that evolve parallel to the X-ray outburst. Our goal is to provide a simple quantitative explanation for the optical flares with an application to a sample of the brightest outbursts of BeXRBs in the Magellanic clouds and the Galactic Ultra-luminous X-ray (ULX) pulsar Swift J0243.6+6124. We constructed a numerical model to study X-ray irradiation in a BeXRB system. We then conducted a parametric investigation of the model parameters and made predictions for the intensity of the optical flares based on geometric and energetic constraints. From our modeling we found that the optical emission during major outbursts is consistent with being the result of X-ray irradiation of the Be disk. For individual systems, if this method is combined with independent constraints of the geometry of the Be disk, the binary orbital plane, and the plane of the observer, it can provide estimates of the Be disk size during major outbursts. Moreover, we computed a semi-analytical relation between optical flare luminosity and X-ray luminosity that is consistent with both model predictions and observed properties of flares.

We study the timing and spectral properties of the Be/X-ray binary pulsar SXP 138 using four NuSTAR observations spanning 2016 to 2017. Analysis of the light curves using the Lomb-Scargle periodogram shows an increase in the spin period of SXP 138 from 140.69 to 140.85 seconds, indicating that the source is in the propeller regime. We calculate the associated rate of spin period change and characterize its non-linearity with a quadratic fit. Pulse profiles obtained by folding the light curves at the spin period show two primary high peaks and two secondary peaks for all the observations. Such features in the pulse profile can result from the combined effect of pencil and fan beam emissions from the two antipodal hotspots accompanied by the relativistic bending of photons. The energy spectra fit with both the blackbody and the power-law spectral model. The best-fit models of all observations show an overall increasing trend in the temperature of the source from 1.8 keV to 2.5 keV, during the first three observations, and a possible decreasing trend in the photon index from 1.8 to 1.7. The power-law component is dominant in the later observations, which is associated with an increase in Compton scattering and the accretion rate. SXP 138 has not been studied in detail, and this work bridges the gap by providing the first comprehensive timing and spectral analysis of this source. These findings improve our understanding of its accretion processes and emission mechanisms, placing SXP 138 in the broader context of Be/X-ray binaries.

Rotational evolution of stellar radiative zones is an old puzzle. We argue that angular momentum (AM) transport by turbulent processes induced by differential rotation is insufficient, and propose that a key role is played by ``magnetic webs." We define magnetic webs as stable magnetic configurations that enforce corotation of their coupled mass shells. Stable magnetic configurations naturally form through relaxation of helical magnetic fields deposited in parts of radiative zones. We discuss the conditions for a magnetic configuration to be sufficiently sturdy to prevent the build up of differential rotation, and conclude that these conditions are easily met in stellar interiors. Low mass stars on the red giant branch (RGB) likely have their compact cores coupled to the lower part of their extended radiative mantle by a magnetic web that was deposited by the receding zone of core convection on the main sequence. This results in moderate core rotation that is broadly consistent with asteroseismic observations, as we illustrate with a stellar evolution model with mass $1.6M_\odot$. Evolving massive stars host more complicated patterns of convective zones that may leave behind many webs, transporting AM towards the surface. Efficient web formation likely results in most massive stars dying with magnetized and slowly rotating cores.

Dark matter could be composed of macroscopic objects with large masses and geometric cross-sections spanning many decades. We investigate the potential interaction of such `stuff-sized' dark matter by considering its interactions with asteroids, planetary rings, and terrestrial bodies. This hail of dark matter could catastrophically destroy these Solar System objects, evaporate them from their orbits, or cause substantial cratering. We estimate these effects and use them to place competitive bounds on a wide, previously-unconstrained swathe of the dark matter parameter space.

Searching for planets analogous to Earth in terms of mass and equilibrium temperature is currently the first step in the quest for habitable conditions outside our Solar System and, ultimately, the search for life in the universe. Future missions such as PLATO or LIFE will begin to detect and characterise these small, cold planets, dedicating significant observation time to them. The aim of this work is to predict which stars are most likely to host an Earth-like planet (ELP) to avoid blind searches, minimises detection times, and thus maximises the number of detections. Using a previous study on correlations between the presence of an ELP and the properties of its system, we trained a Random Forest to recognise and classify systems as 'hosting an ELP' or 'not hosting an ELP'. The Random Forest was trained and tested on populations of synthetic planetary systems derived from the Bern model, and then applied to real observed systems. The tests conducted on the machine learning (ML) model yield precision scores of up to 0.99, indicating that 99% of the systems identified by the model as having ELPs possess at least one. Among the few real observed systems that have been tested, 44 have been selected as having a high probability of hosting an ELP, and a quick study of the stability of these systems confirms that the presence of an Earth-like planet within them would leave them stable. The excellent results obtained from the tests conducted on the ML model demonstrate its ability to recognise the typical architectures of systems with or without ELPs within populations derived from the Bern model. If we assume that the Bern model adequately describes the architecture of real systems, then such a tool can prove indispensable in the search for Earth-like planets. A similar approach could be applied to other planetary system formation models to validate those predictions.

Main-sequence stars follow a well-defined rotation-activity relation. There are two primary regimes: saturated, where the fractional X-ray luminosity $\log(L_{\rm X}/L_*)$ is approximately constant, and unsaturated, where the fractional X-ray luminosity decreases with increasing Rossby number (or decreasing rotation rate). Pre-main sequence (PMS) stars have a larger scatter in $\log(L_{\rm X}/L_*)$ than main-sequence stars, are observed to have saturated levels of X-ray emission, and do not follow the rotation-activity relation. We investigate how PMS stars evolve in the rotation-activity plane and the timescale over which the X-ray rotation-activity relation emerges. Using observational data of $\sim$600 stars from four PMS clusters, stellar internal structure models, a rotational evolution model, and observed X-ray luminosity trends with age, we simulate the evolution of the PMS stars in the rotation-activity plane up to ages of 100 Myr. Our model reproduces the rotation-activity relation found for main-sequence stars, with higher-mass stars beginning to form the unsaturated regime from around 10 Myr. After $\sim$25 Myr, the gradient of the unsaturated regime matches that found for main-sequence stars. For stars of mass greater than 0.6 M$_{\odot}$, the maximum age by which a star has left the saturated regime correlates with when the star leaves the PMS. We find that an intra-cluster age spread is a key factor in contributing to the observed scatter in $\log(L_{\rm X}/L_*)$, particularly for ages < 10 Myr.

The subluminous red nova (SLRN) ZTF SLRN-2020 is the most compelling direct detection of a planet being consumed by its host star, a scenario known as a planetary engulfment event. We present JWST spectroscopy of ZTF SLRN-2020 taken +830 d after its optical emission peak using the NIRSpec fixed-slit $3-5$ $\mu$m high-resolution grating and the MIRI $5-12$ $\mu$m low-resolution spectrometer. NIRSpec reveals the $^{12}$CO fundamental band ($\nu=1-0$) in emission at $\sim4.7$ $\mu$m, Brackett-$\alpha$ emission, and the potential detection of PH$_3$ in emission at $\sim4.3$ $\mu$m. The JWST spectra are consistent with the claim that ZTF SLRN-2020 arose from a planetary engulfment event. We utilize DUSTY to model the late-time $\sim1-12$ $\mu$m spectral energy distribution (SED) of ZTF SLRN-2020, where the best-fit parameters indicate the presence of warm, $720^{+80}_{-50}$ K, circumstellar dust with a total dust mass of Log$\left(\frac{M_\mathrm{d}}{\mathrm{M}_\odot}\right)=-10.61^{+0.08}_{-0.16}$ M$_\odot$. We also fit a DUSTY model to archival photometry taken +320 d after peak that suggested the presence of a cooler, T$_\mathrm{d}=280^{+450}_{-20}$ K, and more massive, Log$\left(\frac{M_\mathrm{d}}{\mathrm{M}_\odot}\right)=-5.89^{+0.29}_{-3.21}$, circumstellar dust component. Assuming the cool component originates from the ZTF SLRN-2020 ejecta, we interpret the warm component as fallback from the ejecta. From the late-time SED model we measure a luminosity of L$_* = 0.29^{+0.03}_{-0.06}$ L$_\odot$ for the remnant host star, which is consistent with a $\sim0.7$ M$_\odot$ K-type star that should not yet have evolved off the main sequence. If ZTF SLRN-2020 was not triggered by stellar evolution, we suggest that the planetary engulfment was due to orbital decay from tidal interactions between the planet and the host star.

In this paper, written in memory of Alexei Starobinsky, we discuss the observational viability of the Ph-$\Lambda_{\rm s}$CDM model - a dynamical dark energy scenario based on a phantom scalar field undergoing an anti-de Sitter (AdS) to de Sitter (dS) transition - and revisit the Sahni-Shtanov braneworld model in light of updated BAO Ly-$\alpha$ data at $z \sim 2.3$. Both models are able to remain consistent with Planck CMB data while offering potential resolutions to the $H_0$ tension. In both cases, the expansion rate $H(z)$ is suppressed relative to Planck-$\Lambda$CDM at high redshift and enhanced at low redshift, while remaining consistent with the comoving distance to recombination as estimated by Planck-$\Lambda$CDM. Comparing model predictions with BAO-inferred values of $H(z)$, we find that SDSS Ly-$\alpha$ data at $z \approx 2.33$ mildly favor such dynamical models, whereas the recent DESI Ly-$\alpha$ measurements agree more closely with $\Lambda$CDM. Although current high-redshift BAO data do not decisively favor one model over another, our findings illustrate how frameworks originally developed to address earlier anomalies - such as the braneworld scenario - may gain renewed relevance in confronting today's cosmological tensions.

We present the discovery of 15 well-resolved giant radio galaxies (GRGs) with angular sizes >5 arcmin and physical sizes >1 Mpc in wide-field Phased Array Feed 944 MHz observations on the Australian Square Kilometre Array Pathfinder (ASKAP). We identify their host galaxies, examine their radio properties as well as their environment, and classify their morphologies as FR I (4), FR II (8), intermediate FR I/II (2), and hybrid (1). The combined ~40 deg^2 ASKAP image of the Sculptor field, which is centred near the starburst galaxy NGC 253, has a resolution of 13" and an rms sensitivity of >10 microJy/beam. The largest GRGs in our sample are ASKAP J0057-2428 (zphot = 0.238), ASKAP J0059-2352 (zphot = 0.735) and ASKAP J0107-2347 (zphot = 0.312), for which we estimate linear projected sizes of 2.7, 3.5 and 3.8 Mpc, respectively. In total we catalog 232 extended radio galaxies of which 77 (33%) are larger than 0.7 Mpc and 35 (15%) are larger than 1 Mpc. The radio galaxy densities are 5.8 deg^-2 (total) and 0.9 (1.9) deg^-2 for those larger than 1 (0.7) Mpc, similar to previous results. Furthermore, we present the ASKAP discovery of a head-tail radio galaxy, a double-lobe radio galaxy with a spiral host, and radio emission from several galaxy clusters. As the ASKAP observations were originally conducted to search for a radio counterpart to the gravitational wave detection GW190814 (z ~ 0.05), we highlight possible host galaxies in our sample.

We present spatially resolved Alma Band-9 observations of the [CII] 158 $\mu$m fine structure line from an optically selected quasar, SDSS J100038.01+020822.4 (J1000), at z=1.8275. By utilizing [OI] 63 $\mu$m line observations from Herschel/PACS and constructing a detailed dust SED using Herschel and Spitzer archival imaging data, we show that the [CII] line emission is well explained by a photodissociation region (PDR) model, in which the emission arises from the surfaces of molecular clouds exposed to far-UV radiation fields $\sim 5\cdot10^3$ times the local interstellar radiation field (G$_0$). We find a factor of 30 variation in spatially resolved [CII]/Far-IR continuum across the source which is explained by the reduced fraction of cooling via [CII] line emission at such high far-UV field strengths. By matching derived PDR parameters to the observed far-IR line and continuum intensities we derive cloud size-scales and find that typical cloud radii in J1000 are $\sim$3.5 pc perhaps indicating an ISM that is highly fractured due to intense star formation activity. We model the galaxy dynamically and find that the [CII] emission is contained within a compact, dynamically cold disk with v/$\sigma$=6.2, consistent with cosmological simulations. We also report the discovery of a companion galaxy to j1000 confirmed by the detection of [CII] and use recently obtained JWST/NirCAM imaging of the system to argue for J1000 being an interacting system. With total stellar mass $\sim 1.5 \times 10^{10}$ M$_\odot$ and main-component dynamical mass $\gtrsim 10^{11}$ M$_\odot$, the J1000 system is a progenitor to the most massive galaxies seen in the local Universe.

We present the results obtained from timing and spectral studies of 15 thermonuclear X-ray bursts from 4U 1820-30 observed with the Neutron Star Interior Composition Explorer (NICER) during its five years of observations between 2017-2022. All bursts showed clear signs of photospheric radius expansion, where the neutron star (NS) photosphere expanded more than 50 km above the surface. One of the bursts produced a super-expansion with a blackbody emission radius of 902 km for the first time with NICER. We searched for burst oscillations in all 15 bursts and found evidence of a coherent oscillation at 716 Hz in a burst, with a 2.9$\sigma$ detection level based on Monte Carlo simulations. If confirmed with future observations, 4U 1820-30 would become the fastest-spinning NS known in X-ray binary systems. The fractional rms amplitude of the candidate burst oscillation was found to be 5.8% in the energy range of 3-10 keV. Following the variable persistent model from burst time-resolved spectroscopy, an anti-correlation is seen between the maximum scaling factor value and the (pre-burst) persistent flux. We detected a low value of ionization at the peak of each burst based on reflection modeling of burst spectra. A partially interacting inner accretion disk or a weakly ionized outer disk may cause the observed ionization dip during the photospheric radius expansion phase.

We study 15 thermonuclear X-ray bursts from 4U 1820--30 observed with the Neutron Star Interior Composition Explorer (NICER). We find evidence of a narrow emission line at 1.0 keV and three absorption lines at 1.7, 3.0, and 3.75 keV, primarily around the photospheric radius expansion phase of most bursts. The 1.0 keV emission line remains constant, while the absorption features, attributed to wind-ejected species, are stable but show slight energy shifts, likely due to combined effects of Doppler and gravitational redshifts. We also examine with NICER the ``aftermath'' of a long X-ray burst (a candidate superburst observed by MAXI) on 2021 August 23 and 24. The aftermath emission recovers within half a day from a flux depression. During this recovery phase, we detect two emission lines at 0.7 and 1 keV, along with three absorption lines whose energies decreased to 1.57, 2.64, and 3.64 keV. Given the nature of the helium white-dwarf companion, these absorption lines during the aftermath may originate from an accretion flow, but only if the accretion environment is significantly contaminated by nuclear ashes from the superburst. This provides evidence of temporary metal enhancement in the accreted material due to strong wind loss. Moreover, we suggest that the absorption features observed during the short X-ray bursts and in the superburst aftermath share a common origin in heavy nuclear ashes enriched with elements like Si, Ar, Ca, or Ti, either from the burst wind or from an accretion flow contaminated by the burst wind.

Before solar eruptions, a short-term slow-rise phase is often observed, during which the pre-eruption structure ascends at speeds much greater than the photospheric motions but much less than those of the eruption phase. Numerical magnetohydrodynamic (MHD) simulations of the coronal evolution driven by photospheric motions up to eruptions have been used to explain the slow-rise phase, but their bottom driving speeds are much larger than realistic photospheric values. Therefore, it remains an open question how the excessively fast bottom driving impacts the slow-rise phase. Here we modelled the slow-rise phase before eruption initiated from a continuously sheared magnetic arcade. In particular, we performed a series of experiments with the bottom driving speed unprecedentedly approaching the photospheric value of around $1$ km s$^{-1}$. The simulations confirmed that the slow-rise phase is an ideal MHD process, i.e., a manifestation of the growing expansion of the sheared arcade in the process of approaching a fully open field state. The overlying field line above the core flux has a slow-rise speed modulated by the driving speed's magnitude but is always over an order of magnitude larger than the driving speed. The core field also expands with speed much higher than the driving speed but much lower than that of the overlying field. By incrementally reducing the bottom-driving speed to realistic photospheric values, we anticipate better matches between the simulated slow-rise speeds and some observed ones.

Using galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) along with haloes from the dark matter only constrained ELUCID (Exploring the Local Universe with the reConstructed Initial Density field) simulation, we examine the properties of galaxies and haloes with respect to their distance to cosmic filaments, determined by the medial-axis thinning technique of the COsmic Web Skeleton (COWS) method. Our findings suggest that galaxies or subhaloes grow in mass as they approach these filaments. Galaxies exhibit a redder colour and diminished specific star formation rates as they approach these filaments. Additionally, older subhaloes tend to be more common near the central regions of these filaments. Elliptical galaxies are more frequently found than spiral galaxies in the central regions of the filaments. Lower-mass galaxies typically display reduced sizes in proximity to filaments, whereas higher-mass galaxies tend to exhibit increased sizes when close to filaments. Moreover, the concentration and spin of the haloes grow as they approach the filaments. These findings support the notion that the large-scale structure of the universe, characterized by cosmic web structures, plays a vital role in shaping galaxy and halo properties.

Recently, polarization angle (PA) orthogonal jumps over millisecond timescale were discovered from three bursts of a repeating fast radio burst source FRB 20201124A by the FAST telescope. We investigate the physical implications of this phenomenon. In general, PA jumps can arise from the superposition of two electromagnetic waves, either coherently or incoherently, as the dominance of the two orthogonal modes switches. In the coherent case, PA jumps occur when linear polarization reaches a minimum and circular polarization peaks, with the total polarization degree conserved. However, incoherent superposition can lead to depolarization. The observations seem to be more consistent with incoherent superposition. The amplitudes of the two orthogonal modes are required to be comparable when jumps occur, placing constraints on the intrinsic radiation mechanisms. We provide general constraints on FRB emission and propagation mechanisms based on the data. Physically, it is difficult to produce PA jumps by generating two orthogonal modes within millisecond timescales, and a geometric effect due to sweeping line-of-sight is a more plausible reason. This requires the emission region to be within the magnetosphere of a spinning central engine, likely a magnetar. The two orthogonal modes may be produced by intrinsic radiation mechanisms or Alfvén-O-mode transition. Plasma birefringence is not easy to achieve when the plasma is moving relativistically. Curvature radiation predicts $|E_{\rm X}/E_{\rm O}|\gtrsim1$, and is difficult to produce jumps; whereas inverse Compton scattering can achieve the transition amplitude ratio $|E_{\rm X}/E_{\rm O}|=1$ to allow jumps to occur under special geometric configurations.

The cosmic growth rate, which is related to peculiar velocity and is a primary scientific objective of galaxy spectroscopic surveys, can be inferred from the Redshift Space Distortion effect and the kinetic Sunyaev-Zel'dovich effect. However, the reconstruction noise power spectrum of the radial velocity field in kSZ is significantly dependent on the measurement of the small-scale galaxy-electron power spectrum $P_{ge}$. In this study, we thoroughly discuss the enhancement of cosmic growth rate measurements facilitated by Fast Radio Bursts, which probe the electron density of the universe along their propagation paths to provide crucial additional information on $P_{ge}$. Subsequently, we utilize future spectroscopic surveys from the Chinese Space Station Telescope and the CMB-S4 experiment, combined with FRB dispersion measures, to achieve precise measurements of the cosmic growth rate at redshifts $z_g = 0.15,0.45,0.75$. Employing Fisher matrix forecasting analysis, we anticipate that constraints on $f\sigma_8$ will reach a precision of 0.001 with a sample size of $10^6$ FRBs. Furthermore, we perform a global analysis using Markov Chain Monte Carlo methods to constrain key parameters of three distinct dark energy models and a modified gravity model based on cosmic growth rate measurements. The results demonstrate that these refined $f\sigma_8$ measurements considerably enhance the constraints on relevant cosmological parameters compared to those obtained from Planck. As the number of observed FRBs increases, alongside more precise galaxy surveys and next-generation CMB observations, new opportunities will arise for constraining cosmological models using the kSZ effect and for developing novel cosmological applications of FRBs.

Polycyclic aromatic hydrocarbon (PAH) ions are crucial intermediates in interstellar chemistry and may play a key role in the infrared emission features observed in space. Here, we investigate the infrared spectra of the indenyl (C$_9$H$_7^-$) and fluorenyl (C$_{13}$H$_9^-$) anions and the indenyl cation (C$_9$H$_7^+$) using infrared pre-dissociation (IRPD) spectroscopy. The experiments were performed in a cryogenic 22 pole ion trap at the FELion beamline of the tunable free-electron laser FELIX. Spectral analysis of the two anionic PAHs, in combination with density functional theory (DFT) computations, revealed key vibrational modes near 1300 cm$^{-1}$, making these ions potential carriers of the 7.7 {\mu}m PAH emission band seen in many astronomical objects. The feature-rich spectrum of cationic indenyl could not be entirely explained by modeling through time-independent anharmonic DFT calculations. Although a better match has been achieved through molecular dynamics simulations, we cannot completely rule out the presence of multiple cationic isomers of the H-loss fragments of indene in the experiments.

The number of resident space objects is rising at an alarming rate. Mega-constellations and breakup events are proliferating in most orbital regimes, and safe navigation is becoming increasingly problematic. It is important to be able to track RSOs accurately and at an affordable computational cost. Orbital dynamics are highly nonlinear, and current operational methods assume Gaussian representations of the objects' states and employ linearizations which cease to hold true in observation-free propagation. Monte Carlo-based filters can provide a means to approximate the a posteriori probability distribution of the states more accurately by providing support in the portion of the state space which overlaps the most with the processed observations. Moreover, dynamical models are not able to capture the full extent of realistic forces experienced in the near-Earth space environment, and hence fully deterministic propagation methods may fail to achieve the desired accuracy. By modeling orbital dynamics as a stochastic system and solving it using stochastic numerical integrators, we are able to simultaneously estimate the scale of the process noise incurred by the assumed uncertainty in the system, and robustly track the state of the spacecraft. In order to find an adequate balance between accuracy and computational cost, we propose three algorithms which are capable of tracking a space object and estimating the magnitude of the system's uncertainty. The proposed filters are successfully applied to a LEO scenario, demonstrating the ability to accurately track a spacecraft state and estimate the scale of the uncertainty online, in various simulation setups.

A sizeable fraction of gamma-ray burst (GRB) light curves (LCs) features a sequence of peaks, which holds information on the unknown way energy is dissipated into gamma-rays over time. Traditional searches for periodic signals in GRB LCs turned out to be inconclusive, partly because they are challenging as a consequence of the short-lived, coloured-noise, and non-stationary nature of the LCs themselves. Yet, recent claims have revived the issue. We searched for periodic components in GRB LCs through a new approach to GRBs, which escapes most of the issues faced by traditional techniques. We identified peaks through a well tested algorithm and selected GRBs with at least 10 peaks out of 5 GRB catalogues (Swift/BAT, CGRO/BATSE, Fermi/GBM, Insight-HXMT, BeppoSAX/GRBM). Each GRB was simply treated as a discrete point process, whose realisation coincides with the sequence of peak times. We searched for possible periodic recurrences based on the multinomial distribution, after accounting for the clustering of peaks due to the non-stationarity of the GRB signals. The best candidate has a p-value of 3e-4 that there is no periodic recurrence. However, accounting for the multiple trials of 555 searched GRBs, its statistical significance is demoted to 17%. The overall distribution of the p-values obtained for all GRBs is compatible with a uniform distribution in [0,1]. We found no robust evidence for multi-peaked GRBs with periodic recurrences. We can exclude that a sizeable fraction (>~ 0.75) of peaks of each GRB with at least 10 peaks are periodic. While our result does not necessarily clash with claimed periodicities based on Fourier techniques, it constrains the putative recurrent behaviour, which would not manifest itself through the sequence of peaks, but, evidently, in a more elusive way.

We provide an updated test for modifications of gravity from a sample of wide-binary stars from GAIA DR3, and their sky-projected relative velocities. Here we extend on our earlier 2023 study, using several updated selection cuts aimed at reducing contamination from triple systems with an undetected third star. We also use improved mass estimates from FLAMES, and we add refinements to previous modelling of the triple and other populations and the model-fitting. We fit histograms of observed vs Newtonian velocity differences to a flexible mixture of binary + triple populations with realistic eccentricity distributions, plus unbound flyby and random-chance populations. We find as before that Newtonian models provide a significantly better fit than MOND, though improved understanding of the triple population is necessary to make this fully decisive.

Classical Be stars are well known to eject mass, but the details governing the initial distribution and evolution of this matter into a disk are poorly constrained by observations. By combining high-cadence spectroscopy with contemporaneous space photometry from TESS, we have sampled about 30 mass ejection events in 13 Be stars. Our goal is to constrain the geometrical and kinematic properties of the ejecta, facilitating the investigation into the initial conditions and evolution, and understanding its interactions with preexisting material. The photometric variability is analyzed together with measurements of the rapidly changing emission features to identify the onset of outburst events and obtain information about the geometry of the ejecta and its evolution. All Be stars observed with sufficiently high cadence exhibit rapid oscillations of line asymmetry with a single frequency in the days following the start of the event. The emission asymmetry cycles break down after roughly 5 - 10 cycles, with the emission line profile converging toward approximate symmetry. In photometry, several frequencies typically emerge at relatively high amplitude at some point during the mass ejection process. In all observed cases, freshly ejected material was initially within a narrow azimuthal range, indicating it was launched from a localized region on the star. The material orbits the star with a frequency consistent with the near-surface Keplerian orbital frequency. This material circularizes into a disk configuration after several orbital timescales. This is true whether or not there was a preexisting disk. We find no evidence for precursor phases prior to the ejection of mass in our sample. The several photometric frequencies that emerge during outburst are at least partially stellar in origin. (Abstract abridged)

The halo mass function (HMF) is fundamental for interpreting the number counts of galaxy clusters, serving as a pivotal theoretical tool in cosmology. With the advent of high-precision surveys such as LSST, eROSITA, DESI, and Euclid, accurate HMF modeling becomes indispensable to avoid systematic biases in cosmological parameter estimation from cluster cosmology. Moreover, these surveys aim to shed light on the dark sector and uncover dark energy's puzzling nature, necessitating models that faithfully capture its features to ensure robust parameter inference. We aim to construct a model for the HMF in dynamical dark energy cosmologies that preserves the accuracy achieved for the standard $\Lambda (\nu)$CDM model of cosmology, while meeting the precision requirements necessary for future cosmological surveys. Our approach models the HMF parameters as functions of the deceleration parameter at the turnaround, a quantity shown to encapsulate essential information regarding the impact of dynamical dark energy on structure formation. We calibrate the model using results from a comprehensive suite of $N$-body simulations spanning various cosmological scenarios, ensuring sub-percent systematic accuracy. We present an HMF model tailored for dynamical dark energy cosmologies. The model is calibrated following a Bayesian approach, and its uncertainty is characterized by a single parameter controlling its systematic error, which remains at the sub-percent level. This ensures that theoretical uncertainties from our model are subdominant relative to other error sources in future cluster number counts analyses.

Aims. We analyze the clustering of galaxy clusters in a large contiguous sample, the Constrain Dark Energy with X-ray (CODEX) sample. We construct a likelihood for cosmological parameters by comparing the measured clustering signal and a theoretical prediction, and use this to obtain parameter constraints. Methods. We measured the three multipole moments (monopole, quadrupole, and hexadecapole, $\ell = 0, 2, 4$) of the power spectrum of a subset of the CODEX clusters. To fully model cluster clustering, we also determined the expected clustering bias of the sample using estimates for the cluster masses and a mass-to-bias model calibrated using N-body simulations. We estimated the covariance matrix of the measured power spectrum multipoles using a set of simulated dark-matter halo catalogs. Combining all these ingredients, we performed a Markov chain Monte Carlo sampling of cosmological parameters $\Omega_m$ and $\sigma_8$ to obtain their posterior. Results. We found the CODEX clustering signal to be consistent with an earlier X-ray selected cluster sample, the REFLEX II sample. We also found that the measured power spectrum multipoles are compatible with the predicted, bias-scaled linear matter power spectrum when the cosmological parameters determined by the Planck satellite are assumed. Furthermore, we found the marginalized parameter constraints of $\Omega_m = 0.24^{+0.06}_{-0.04}$ and $\sigma_8 = 1.13^{+0.43}_{-0.24}$. The full 2D posterior is consistent, for example, with the Planck cosmology within the 68% confidence region.

Many neutron star low-mass X-ray binaries (NS LMXBs) with short orbital periods (~hours) cycle between outburst and quiescent phases, and thus provide an excellent way to study the accretion process. The cause of such outbursts is believed to be thermal-viscous instability in the accretion disc. However, some of these transient sources show unusually long outbursts. For example, EXO 0748-676 remained in outburst for at least 23 years before entering a quiescence, only to re-emerge 16 years later. We aim to investigate if such long outbursts could be due to the usual disc instability, or if any other mechanism is required. In order to address this question, we systematically compare various properties of long outburst and short outburst NS LMXBs. For this, we analyze the long-term X-ray light curves of many short orbital period (hours) NS LMXBs, examining the outburst duration and the inferred accretion rate, and estimate the accretion disc mass. Our study shows that long outburst sources are well-separated from the short outburst ones in parameter spaces involving accretion rate, disc mass, outburst duration, etc. in four ways. This implies that the thermal-viscous instability in the disc cannot explain the long outbursts, but could explain the short ones. Moreover, we discuss that both donor star related and disc related models have difficulties in explaining long outbursts. Our finding will be crucial to understanding the accretion process of transiently accreting neutron stars and black holes.

Time-delay cosmography (TDC) using multiply-lensed quasars (QSOs) by galaxies has recently emerged as an independent and competitive tool to measure the value of the Hubble constant. Lens galaxy clusters hosting multiply-imaged QSOs, when coupled with an accurate and precise knowledge of their total mass distribution, are equally powerful cosmological probes. However, less than ten such systems have been identified to date. Our study aims to expand the limited sample of cluster-lensed QSO systems by identifying new candidates within rich galaxy clusters. Starting from a sample of ~$10^5$ galaxy cluster candidates (Wen & Han, 2022), built from Dark Energy Survey and Wide-field Infrared Survey Explorer imaging data, and a highly-pure catalogue of over one million QSOs, based on Gaia DR3 data, we cross-correlate them to identify candidate lensed QSOs near the core of massive galaxy clusters. Our search yielded 3 lensed double candidates over an area of ~$5000$ sq. degree. In this work, we focus on the best candidate consisting of a double QSO with Gaia-based redshift of 1.35, projected behind a moderately rich cluster (WHJ0400-27) at $z_{phot}=0.65$. Based on a first spectroscopic follow-up study, we confirm the two QSOs at $z=1.345$, with indistinguishable spectra, and a brightest cluster galaxy at $z=0.626$. These observations seem to support the strong lensing nature of this system, although some tension emerges when the cluster mass from a preliminary lens model is compared with that from other mass proxies. We also discuss the possibility that such system is a rare physical association of two distinct QSOs with a projected physical distance of ~$150$ kpc. If further spectroscopic observations confirm its lensing nature, such a rare lens system would exhibit one of the largest image separations observed to date ($\Delta\vartheta=17.8''$), opening interesting TDC applications.

The recent observational evidence of deviations from the $\Lambda$-Cold Dark Matter ($\Lambda$CDM) model points towards the presence of evolving dark energy. The simplest possibility consists of a cosmological scalar field $\varphi$, dubbed quintessence, driving the accelerated expansion. We assess the evidence for the existence of such a scalar field. We find that, if the accelerated expansion is driven by quintessence, the data favour a potential energy $V(\varphi)$ that is concave, i.e., $m^2=d^2V/d\varphi^2<0$. Furthermore, and more significantly, the data strongly favour a scalar field that is non-minimally coupled to gravity (Bayes factor $\log(B) = 7.34 \pm 0.60$), leading to time variations in the gravitational constant on cosmological scales, and the existence of fifth forces on smaller scales. The fact that we do not observe such fifth forces implies that either new physics must come into play on non-cosmological scales or that quintessence is an unlikely explanation for the observed cosmic acceleration.

Photometric redshifts of galaxies obtained by multi-wavelength data are widely used in photometric surveys because of its high efficiency. Although various methods have been developed, template fitting is still adopted as one of the most popular approaches. Its accuracy strongly depends on the quality of the Spectral Energy Distribution (SED) templates, which can be calibrated using broadband photometric data from galaxies with known spectroscopic redshifts. Such calibration is expected to improve photometric redshift accuracy, as the calibrated templates will align with observed photometric data more closely. The upcoming China Space Station Survey Telescope (CSST) is one of the Stage IV surveys, which aiming for high precision cosmological studies. To improve the accuracy of photometric redshift estimation for CSST, we calibrated the CWW+KIN templates using a perturbation algorithm with broadband photometric data from the CSST mock catalog. This calibration used a training set consisting of approximately 4,500 galaxies, which is 10% of the total galaxy sample. The outlier fraction and scatter of the photometric redshifts derived from the calibrated templates are 2.55% and 0.036, respectively. Compared to the CWW+KIN templates, these values are reduced by 34% and 23%, respectively. This demonstrates that SED templates calibrated with a small training set can effectively optimize photometric redshift accuracy for future large-scale surveys like CSST, especially with limited spectral training data.

High order corrections to the perihelion precession are obtained in non-Newtonian central potentials, via complex analysis techniques. The result is an exact series expansion whose terms, for a perturbation of the form $\delta V=\frac{\gamma}{r^{s}}$, are calculated in closed form. To validate the method, the series is applied to the specific case of s=3, and the results are compared with those presented in literature, which are relate to the Schwarzschild metric. As a further test, a numerical simulation was carried out for the case where s=4. The algebraic calculations and numerical simulations were carried out via software with symbolic capabilities.

High-frequency gravitational-wave (GW) radiation has been detected by LIGO-Virgo-KAGRA in the merger of compact stars. However, two GW events, GW190814 and GW200210, the mass of one companion object falls into the mass region of $(2.2-3)\rm~M_\odot$, and how to identify such object (e.g., as a low-mass black hole (BH) or a massive neutron star (NS)) remains an open question. In this paper, we propose a method to identify the mystery compact object (MCO) with the mass region of $(2.2-3)\rm~M_\odot$ in a binary system via the possible electromagnetic (EM) radiations before and after the mergers. A multi-band EM emission can be produced with $L\propto(-t)^{7/4}$ (or $L\propto(-t)^{-5/4}$) during the inspiral phase due to the BH battery (or interaction magnetospheres) mechanism, and a bright (or dark) kilonova emission is powered by radioactive decay with ejecta mass ratio $q>1.7$ (or $q<1.7$) during the post-merge state when MCO is as a low-mass BH (or massive NS) to merger with NS. Moreover, by considering the merger system between MCO and a BH when MCO is a massive NS, we find that it requires the BH with high spin (e.g., $a\sim0.8-0.99$) to make sure the tidal disruption event (TDE) occurred, and a multi-band precursor emission and bright kilonova emission can also be produced during the inspiral phase and post-merge state, respectively. In any case, no matter which mechanism we adopt, such precursor emissions are too weak to be detected by most current telescopes unless the distance is close enough.

The extreme ultraviolet (EUV) and X-ray radiation emitted during solar flares has been shown to significantly increase the electron density of the Earth's ionosphere. During flares, quasi-periodic pulsations (QPPs) in X-ray flux originating in the corona have previously been linked to subsequent pulsations in the Earth's ionospheric D-region. Similar pulsations have been detected in chromospheric EUV emission, although their impact on the Earth's ionosphere has not previously been investigated. Here, for the first time, synchronous pulsations were detected in solar EUV emission and ionospheric Total Electron Content (TEC) measurements. Using wavelet and periodogram analysis, we detect QPPs with approximately 85 second periods in chromospheric EUV emission lines (He II 304 Å, C III 977 Å and H I 972 Å) from the Solar Dynamics Observatory Extreme Ultraviolet Variability Experiment (SDO/EVE) during the impulsive phase of an X5.4 flare on March 7, 2012. These lines contribute to ionization in the ionospheric E- and F-regions, resulting in subsequent variations of electron density with the same periodicity, which was detected in TEC measurements. This work demonstrates that the Earth's ionosphere is responsive to fine-scale fluctuations in EUV emission during flares, with a time delay of approximately 30 seconds found. These findings may have applications in atmospheric modelling and solar-terrestrial studies, including the calculation of ionospheric recombination rates.

Low-thrust, many-revolution transfers between near-rectilinear halo orbits and low lunar orbits are challenging due to the many-revolutions and is further complicated by three-body perturbation. To address these challenges, we extend hybrid differential dynamic programming by enhancing with a continuation of dynamical system. The optimization begins with the Sundman-transformed two-body problem and gradually transitions to the Sundman-transformed circular restricted three-body problem expressed in the moon-centered inertial frame. Numerical examples demonstrate the robust convergence of our method, where optimal transfers from low lunar orbit to near-rectilinear halo orbit are obtained with a poor initial guess of low lunar orbit.

Active Galactic Nuclei (AGN) host accreting supermassive black holes (SMBHs). The accretion can lead to the formation of a hot, X-ray emitting corona close to the SMBH capable of accelerating relativistic electrons. Observations in the millimetre (mm) band can probe its synchrotron emission. We provide a framework to derive physical information of SMBH coronae by modelling their spectral energy distribution (SED) from radio to far infrared frequencies. We also explore the possibilities of deriving additional information from mm observations, such as the SMBH mass, and studying high-redshift lensed sources. We introduce a corona emission model based on a one-zone spherical region with a hybrid thermal and non-thermal plasma. We investigate in detail how the corona SED depends on different parameters such as size, opacity, and magnetic field strength. Other galactic emission components from dust, ionised gas and diffuse relativistic electrons are also included in the SED fitting scheme. We apply our code consistently to a sample of radio-quiet AGN with strong indications of a coronal component in the mm. The detected mm emission from SMBH coronae is consistent with having a non-thermal relativistic particle population with an energy density that is ~0.5-10% of that in the thermal plasma. This requires magnetic energy densities close to equipartition with the thermal gas, and corona sizes of 60-250 gravitational radii. The model can also reproduce the observed correlation between mm emission and SMBH mass when accounting for uncertainties in the corona size. The mm band offers a unique window into the physics of SMBH coronae, enabling the study of highly dust-obscured sources and high-redshift lensed quasars. Gaining a deeper understanding of the relativistic particle population in SMBH coronae can provide key insights into their potential multiwavelength and neutrino emission.

Fast moving celestial objects are characterized by velocities across the celestial sphere that significantly differ from the motions of background stars. In observational images, these objects exhibit distinct shapes, contrasting with the typical appearances of stars. Depending on the observational method employed, these celestial entities may be designated as near-Earth objects or asteroids. Historically, fast moving celestial objects have been observed using ground-based telescopes, where the relative stability of stars and Earth facilitated effective image differencing techniques alongside traditional fast moving celestial object detection and classification algorithms. However, the growing prevalence of space-based telescopes, along with their diverse observational modes, produces images with different properties, rendering conventional methods less effective. This paper presents a novel algorithm for detecting fast moving celestial objects within star fields. Our approach enhances state-of-the-art fast moving celestial object detection neural networks by transforming them into physical-inspired neural networks. These neural networks leverage the point spread function of the telescope and the specific observational mode as prior information; they can directly identify moving fast moving celestial objects within star fields without requiring additional training, thereby addressing the limitations of traditional techniques. Additionally, all neural networks are integrated using the mixture of experts technique, forming a comprehensive fast moving celestial object detection algorithm. We have evaluated our algorithm using simulated observational data that mimics various observations carried out by space based telescope scenarios and real observation images. Results demonstrate that our method effectively detects fast moving celestial objects across different observational modes.

This paper presents a study of the dependence of the simulated intensities of recombination lines from hydrogen and helium atoms on the number of $n\ell$-resolved principal quantum numbers included in the calculations. We simulate hydrogen and helium emitting astrophysical plasmas using the code Cloudy and show that, if not enough $n\ell$-resolved levels are included, recombination line intensities can be predicted with significant errors than can be more than 30\% for H~I IR lines and 10\% for He~I optical lines ($\sim$20\% for He~I IR recombination lines) at densities $\sim1\text{cm}^{-3}$, comparable to interstellar medium. This can have consequences in several spectroscopic studies where high accuracy is required, such as primordial helium abundance determination. Our results indicate that the minimum number of resolved levels included in the simulated hydrogen and helium ions of our spectral emission models should be adjusted to the specific lines to be predicted, as well as to the temperature and density conditions of the simulated plasma.

Recent DESI DR2 observations indicate that dark energy has crossed from phantom to quintessence regime, a behavior known as the quintom-B realization. In this work we constrain dynamical dark energy and modified gravity using the swampland Trans-Planckian Censorship Conjecture (TCC), which forbids eternal acceleration since in this case any trans-Planckian quantum fluctuation would eventually stretch beyond the Hubble radius, breaking the applicability of any effective field theory and cosmological techniques. By combining DESI DR2 data with the TCC criterion, we impose tight constraints on the dark energy equation of state and its parameter space in scenarios such as the Chevallier-Polarski-Linder, Barboza-Alcaniz, Jassal-Bagla-Padmanabhan, EXP and LOG parameterizations, significantly constraining the quintom-A behavior. Also we examine models within the framework of $f(T)$ and $f(Q)$ modified gravity theories, demonstrating that TCC is very powerful to constrain or exclude them, a result that indicates the necessity to consider infrared modifications on General Relativity apart from the usual ultraviolet ones. Our findings imply that viable dynamical dark energy scenarios must asymptotically transit to deceleration, shedding light on new physics consistent with both cosmological observations and quantum gravity principles.

Being one of the first exoplanets observed by the James Webb Space Telescope (JWST), WASP-39 b has become an iconic target and many transit spectra recorded with different instruments (NIRISS, NIRCAM, NIRSpec G395H, NIRSpec PRISM and MIRI) are currently available, allowing in-depth studies of its atmosphere. We present here a novel approach to interpret WASP-39 b's transit spectroscopic data, consisting of a multi-step process where ab initio equilibrium chemistry models and blind retrievals are used iteratively to find physically robust, optimal solutions. Following this approach, we have identified a new scenario to explain WASP-39 b's atmospheric composition, in which silicon-based chemistry plays a major role. In this scenario, SiO may explain the spectral absorption at 4.1 $\mu$m, currently interpreted as being due to SO$_2$. SiO and the other gas species identified by the retrieval models, i.e. H$_2$O, CO$_2$, Na and K, are consistent with an atmosphere in chemical equilibrium with a temperature-pressure profile constrained by H$_2$O and CO$_2$ absorption bands. In addition, silicate clouds and hazes can produce the spectral features observed by MIRI in the spectral window 5-12 $\mu$m. While we advocate the need for more data, possibly at higher spectral resolution, to confirm our results for WASP-39 b's atmospheric composition, we highlight a refined atmospheric retrieval strategy with pre-selection and post-reconstruction to guide the next generation of transit spectroscopy.

Affleck-Dine (AD) baryogenesis can produce the baryon asymmetry of the Universe through the $CP$-violating dynamics of AD field. The field generally fragments into Q-balls, whose rapid decay induces enhanced gravitational waves. In this Letter, we investigate the anisotropies in this gravitational wave background as a new essential observable for AD baryogenesis. The evolution of AD field causes non-Gaussian baryonic isocurvature perturbations, and the non-Gaussianity modulates the spatial distribution of Q-balls on large scales, resulting in large-scale anisotropies in the Q-ball-induced gravitational wave background. We present that the anisotropies can be significantly large with a reduced angular power spectrum $\sim 10^{-2}$, and can be detected by future experiments like LISA. Moreover, these anisotropies universally reveal the $CP$-violating dynamics of AD field, opening a novel road to explore the longstanding baryon asymmetry puzzle.

Binary black hole systems (BBHs) have become a vivid reality in astrophysics as stellar-mass black hole mergers are now detected through their related gravitational wave emission during the merger stage. If many studies were recently dedicated to the last stages of BBH where black holes are surrounded by a circumbinary disk (CBD), the structure of these systems prior to the formation of the CBD remains mostly unexplored. The aim of the present article is to investigate the potential modifications induced by the presence of a secondary black hole onto the structure of the accretion disk surrounding the primary black hole. We performed 2D classical hydrodynamical simulations of an accretion disk surrounding the primary black hole while taking into account all gravitational effects induced by both the primary black hole and the secondary black hole orbiting on circular orbits around the center of mass of the this http URL report three main effects of the presence of a secondary black hole orbiting a circular orbit beyond the outer edge of the accretion disk: 1/ the outer radius of the accretion disk is significantly reduced and its ratio to the black hole separation is directly linked to only the mass ratio of the black holes; 2/ two spiral arms are visible in the gas density structure of the disk and 3/ the outer edge of the accretion disk exhibits an elliptical shape that mainly depends on the mass ratio of the black holes. Our results show that an accretion disk orbiting a primary black hole in a pre-CBD BBH exhibits specific features induced by the gravitational force generated by the presence of a secondary black hole beyond its outer edge. Such features, directly linked to the binary separation and mass ratio, has therefore the potential to help in the search and identification of BBH in the pre-CBD stage.

The inflationary diffusion of (pseudo-)scalar fields with discrete symmetries can seed the formation of a gas of closed domain walls after inflation, when the distance between degenerate minima in field space is not too far from the inflationary Hubble scale. Primordial black holes (PBHs) can then be formed once sufficiently heavy domain walls re-enter the Hubble sphere. In this scenario, inflation determines a distinctive PBH mass distribution that is rather flat and can thus lead to a sizable total abundance of PBHs, while avoiding some of the downsides of PBH formation from critical collapse. We show that generic QCD axion models, with decay constant close to the inflationary Hubble scale, can yield up to $1\%$ of the dark matter (DM) today in the form of PBHs, while being compatible with isocurvature constraints from Cosmic Microwave Background observations. This occurs for values of axion decay constants around $f_a\simeq 10^{8}~\text{GeV}$, that is the region targeted by axion helioscopes and partially constrained by astrophysical observations. The resulting PBHs have \textit{stupendously} large masses, above $10^{11}M_\odot$, and their existence can be probed by Large Scale Structure observations. Larger PBH abundances can be generated by axion-like particles. Alternatively, in scenarios where isocurvature constraints can be relaxed, we find that the totality of the DM can be produced by the QCD axion misalignment mechanism, accompanied by a ${\cal O}(10^{-3})$ DM fraction in PBHs of masses $(10^5-10^6)~M_\odot$. These can act as seeds for the formation of massive black holes at large redshifts, as suggested by recent JWST observations.

(abridged) The black hole Sagittarius A* (Sgr A*) is a prime target for next-generation Earth-space very-long-baseline interferometry missions such as the Black Hole Explorer (BHEX), which aims to probe baselines of the order of 20 G$\lambda$. At these baselines, Sgr A* observations will be affected by the diffractive scattering effects from the interstellar medium (ISM). Therefore, we study how different parameter choices for turbulence in the ISM affect BHEX's observational capabilities to probe strong lensing features of Sgr A*. By using a simple geometric model of concentric Gaussian rings for Sgr A*'s photon ring signal and observing at 320 GHz, we find that the BHEX-ALMA baseline has the required sensitivity to observe Sgr A* for a broad range of values of the power-law index of density fluctuations in the ISM and the inner scale of turbulence. For other baselines with moderate sensitivities, a strong need for observations at shorter scales of $\approx$ 13.5 G$\lambda$ is identified. For this purpose, an orbit migration scheme is proposed. It is modeled using both chemical propulsion (CP)-based Hohmann transfers and electric propulsion (EP)-based orbit raising with the result that a CP-based transfer can be performed in a matter of hours, but with a significantly higher fuel requirement as compared to EP, which however requires a transfer time of around 6 weeks. The consequences of these orbits for probing Sgr A*'s spacetime is studied by quantifying the spatial resolution, temporal resolution and the angular sampling of the photon ring signal in the Fourier coverage of each of these orbits. We show that higher orbits isolate spacetime features while sacrificing both, signal lost to scattering and temporal resolution, but gain greater access to the morphology of the photon ring.

Magnetic fields exist in and around galaxies, but the properties of these fields have not been fully explored due to the challenges inherent in observing and modeling them. In this Note, we explore the differences in magnetic field strength of central and satellite galaxies from the magnetohydrodynamic TNG100 simulation. We find that on average, magnetic fields in satellite galaxies are roughly an order of magnitude stronger than those of central galaxies with comparable masses. The difference is greater for satellites that have already approached within $1 R_{200}$ of their host galaxies. These results indicate that magnetic fields in satellite galaxies are amplified by environmental processes as they fall into a host halo.

Short-period Galactic double white dwarf (DWD) systems will be observable both in visible light through photometric monitoring and in mHz-range gravitational waves (GWs) with forthcoming space-based laser interferometry such as LISA. When only photometric variability is used to measure DWD intrinsic properties, there is a degeneracy between the chirp mass and binary tidal interaction, as orbital frequency time derivative is set by both GW radiation and tides. Without expensive radial velocity data from spectroscopic monitoring, this degeneracy may be lifted in principle by directly measuring the second time derivative of the orbital frequency through photometric monitoring over an ultra-long time baseline. Alternatively, the degeneracy can be removed by exploiting information in both photometric variability and the coherent GW waveform. Investigating both approaches, we find that direct measurement of the second time derivative is likely infeasible for most DWDs, while the multi-messenger method will disentangle measurements of the chirp mass and the binary moments of inertia, for a large sample of tidally locked systems. The latter information will enable empirical tests of WD structure models with finite temperature effects.

We dig into the semi-classical description of gravity by studying one-loop corrections to primordial power spectra generated during cosmic inflation from gravitational nonlinear interactions. In the realm of the Effective Field Theory (EFT) of inflationary fluctuations, we renormalize the quadratic Lagrangian dictating the linear dynamics of gauge-invariant perturbations. Since gravity is a non-renormalizable theory, this procedure is performed perturbatively in terms of negative powers of the EFT strong coupling scales. Since the interactions we consider are purely gravitational, they are ubiquitous and independent of the details of the EFT. Our results are thus relevant for a large class of approximately scale-invariant inflationary scenarios, be them driven by a single scalar field with canonical kinetic terms, or with a non-canonical structure \textit{à la} $P(X,\phi)$, or for an effective single-field description at the level of fluctuations only and emerging from a covariant multifield theory. Using dimensional regularization, we show that time-dependent Ultra-Violet (UV) divergences appearing at the loop level can be canceled at all times by an appropriate splitting of the bare Lagrangian into renormalized operators and counterterms. Moreover, we explicitly compute all finite contributions to the loops and we prove that, taking into account backreaction, the final one-loop renormalized power spectra of both the primordial curvature perturbation and of gravitational waves are exactly conserved on super-horizon scales. Conclusions of our work imply that the scalar and tensor propagation speeds are immune to radiative corrections from gravitational nonlinearities. We discuss a first application to multifield inflation.

Decay of helical (hyper)magnetic fields that may have been present in the Universe during the Electroweak epoch can contribute to generation of the baryon asymmetry of the Universe. We revise constraints on the strength and correlation length of such fields from the requirement that their decay does not lead to over-production of the baryon asymmetry. We show that the helical fields with strength down to 1e-5 of the maximal possible strength during the Electroweak epoch should have had their correlation at least ~1e-6 of the Hubble radius during this epoch. For weaker fields this lower bound on the correlation length relaxes proportionally to the square of magnetic field strength. A field with parameters saturating the bound may actually be responsible for the baryon asymmetry observed today. We show that relic of such a field, surviving in the present day Universe in the form of intergalactic magnetic field detectable with Cherenkov Telescope Array Observatory, may have the strength up to 10-100 pG and can have parameters needed to affect the cosmological recombination and relax the Hubble tension. We also show that there is no constraint on the parameters of helical or non-helical magnetic fields stemming from the requirement that the baryon isocurvature perturbations produced by such fields during the Electroweak epoch are within the observational limits.

We investigate scale-dependent modifications to the primordial scalar power spectrum as potential solutions to the Hubble tension. We use the Fisher-bias formalism, recently adapted to examine perturbed recombination solutions to the Hubble tension, and extend its range of validity with an iterative method. We first analyze the Planck cosmic microwave background (CMB) anisotropy data, demonstrating the existence of modifications to the primordial power spectrum capable of fully resolving the tension between Planck and SH0ES. As a proof of concept, we interpret these solutions in terms of small, time-dependent variations in the first slow roll parameter or in the sound speed of curvature perturbations during a stage of primordial inflation. However, these solutions are associated with a low total matter density $\Omega_m$, which makes them inconsistent with baryon acoustic oscillations (BAO) and uncalibrated supernovae (SNIa) data. When incorporating additional BOSS and PantheonPlus data, the solutions that reduce the Hubble tension tend to overfit Planck CMB data to compensate for the worsened fit to BAO and SNIa data, making them less compelling. These findings suggest that modifying the primordial power spectrum alone is unlikely to provide a robust resolution to the tension and highlight how the viability of such data-driven solutions depends on the specific datasets considered, emphasizing the role of future high-precision observations in further constraining possible resolutions to the tension.

Recently, the Alpha Magnetic Spectrometer (AMS-02) Collaboration presented tentative evidence for the detection of cosmic antihelion-3 (${}^3\overline{\rm He}$) events, alongside a comparable number of antideuterons ($\overline{\rm D}$). If confirmed, these observations could revolutionize our understanding of cosmic-ray production and propagation and/or serve as compelling indirect evidence for dark matter. Given that the detection of cosmic $\overline{\rm D}$ is already at the limit of AMS-02 sensitivity, explaining the observation of ${}^3\overline{\rm He}$ even within the standard coalescence framework poses a significant challenge. It has recently been shown that a previously overlooked mechanism within the Standard Model of particle physics-namely, the production of antihelion via the displaced-vertex decay of $\bar{\Lambda}_b^0$ baryons-could substantially enhance the ${}^3\overline{\rm He}$ flux arising from dark matter-induced processes. In light of these challenges, we present a tuning of Pythia that is consistent with LEP data on the fragmentation function of $b$ quarks into $b$-hadrons-a critical factor for determining the $\bar{\Lambda}_b^0$ multiplicity-and with ALICE and ALEPH data for the $\overline{\rm D}$ and ${}^3\overline{\rm He}$ spectra, which we employ to determine our coalescence model. Our refined Pythia tuning, in conjunction with our coalescence model, results in a predicted branching ratio for the production of ${}^3\overline{\rm He}$ from $\bar{\Lambda}_b^0$ decays that is consistent with the recent upper limit measured by LHCb. Furthermore, our prediction indicates that the contribution of $\overline{\rm D}$ and ${}^3\overline{\rm He}$ from beauty-hadron decays is negligible relative to the direct production from hadronization.

We investigate the effect of potentially large distortions of the relic neutrino spectra on cosmological this http URL that end, we consider a phenomenological model of "gray" spectral distributions, described by a single parameter which generalizes the traditional $y$-distortions to possibly large negative values. Implementing these distortions in the primordial nucleosynthesis code PRIMAT, we can constrain the distortion parameter along with the presence of extra radiation, exploiting the complementarity of big bang nucleosynthesis and cosmic microwave background measurements to disentangle gravitational and non-thermal effects. These constraints rule out a distortion where more than $\sim 1/2$ of the neutrinos energy density is replaced by dark radiation. Nonetheless, we find that large distortions, accompanied by extra radiation, are allowed-and even slightly preferred in some cases-by current cosmological observations. As this scenario would require substantial modifications to the physics of neutrino decoupling in the early Universe, these observational constraints call for a renewed attention on the possibility of large deviations from the standard cosmological model in the neutrino sector.

Jupiter's ultraviolet aurora frequently shows a number of arcs between the dusk-side polar region and the main emission, which are denoted as ``bridges''. This work presents a largely automated detection and statistical analysis of bridges over 248 Hubble-Space-Telescope observations, alongside a multi-instrument study of crossings of magnetic field lines connected to bridges by the Juno spacecraft during its first 30 perijoves. Bridges are observed to arise on timescales of $\sim$2 hours, can persist over a full Jupiter rotation, and are conjugate between hemispheres. The appearance of bridges is strongly associated with compression of the magnetosphere by the solar wind. Low-altitude bridge crossings are associated with upward-dominated, broadband electron distributions, consistent with Zone-II aurorae, as well as with plasma-wave bursts observed by Juno-Waves, in agreement with existing theoretical models for the generation of polar-region aurorae. Electron populations associated with crossings of field lines threading the main emission by Juno also become more downward-dominated when separate bridges are present in the nearby aurora. In all, this indicates that bridges are likely Zone-II aurorae that have become spatially separated from the Zone-I aurorae under the influence of the solar wind.

Axion dark matter can satisfy the conditions needed to account for all of the dark matter and solve the strong CP problem. The Axion Dark Matter eXperiment (ADMX) is a direct dark matter search using a haloscope to convert axions to photons in an external magnetic field. Key to this conversion is the use of a microwave resonator that enhances the sensitivity at the frequency of interest. The ADMX experiment boosts its sensitivity using a dilution refrigerator and near quantum-limited amplifier to reduce the noise level in the experimental apparatus. In the most recent run, ADMX searched for axions between 1.10-1.31 GHz to extended Kim-Shifman-Vainshtein-Zakharov (KSVZ) sensitivity. This Letter reports on the results of that run, as well as unique aspects of this experimental setup.

Primordial black holes (PBHs) with masses between $10^{14}$ and $10^{20}$ kg are candidates to contribute a substantial fraction of the total dark matter abundance. When in orbit around the center of a star, which can possibly be a completely interior orbit, such objects would emit gravitational waves, as predicted by general relativity. In this work, we examine the gravitational wave signals emitted by such objects when they orbit typical stars, such as the Sun. We show that the magnitude of the waves that could eventually be detected on Earth from a possible PBH orbiting the Sun or a neighboring Sun-like star within our galaxy can be significantly stronger than those originating from a PBH orbiting a denser but more distant neutron star (NS). Such signals may be detectable by the LISA gravitational-wave detector. In addition, we estimate the contribution that a large collection of such PBH-star systems would make to the stochastic gravitational-wave background (SGWB) within a range of frequencies to which pulsar timing arrays are sensitive.

We investigate a framework for extracting parameters of stochastic gravitational wave background (SGWB) with peak-like templates in the millihertz frequency band, and analyzing transient contamination effects on parameter reconstruction. We present the spectrum and spectrogram under different conditions and provide the results of parameter reconstruction. Using templates from the early universe, we demonstrate that the peak-like templates outperform the broken power law (BPL) templates in power-law exponents recovery and peak frequency localization. The reconstruction results obtained using data from Fast Fourier Transform (FFT) are better than those obtained using data from Short-Time Fourier Transform (STFT) which is based on the spectrogram. For the single-peak template, the estimation accuracy of the exponent and peak frequency surpasses that of the BPL template by an order of magnitude, but demonstrates less precision in amplitude estimation compared to BPL. Regarding the double-peak template, parameter estimation results derived from the FFT methodology consistently outperform those obtained using STFT. Nevertheless, transient signals exhibit a detrimental impact on parameter estimation precision, causing errors to increase by an order of magnitude, particularly in multi-peak scenarios. This framework provides an example for using templates to analyze data from space-based gravitational wave detectors.

We investigate tidal heating associated with the binary inspiral of strange quark stars and its impact on the resulting gravitational wave signal. Tidal heating during the merger of neutron stars composed of nuclear matter may be considered negligible, but it has been demonstrated recently that the presence of hyperons at high densities could significantly enhance the dissipation during inspiral. In this work, we evaluate the bulk viscosity arising from non-leptonic weak processes involving quarks and show that it can be several orders of magnitude higher than the viscosity of nuclear matter at temperatures relevant to the inspiral phase of the merger of strange stars. We model strange quark matter in the normal phase using a non-ideal bag model including electrons and ensure compatibility with astrophysical constraints. By analysing equal-mass binary systems with component masses ranging from 1.4 to 1.8 $\, M_{\odot}$, we find that temperatures close to 0.1 MeV are reached by the end of the inspiral phase. We also estimate the effect on the gravitational waveform and conclude that the additional phase shift could range from $0.1$ to $0.5$ radians for strange quark masses of 200 MeV, making it potentially detectable by next-generation gravitational wave detectors. Given that tidal heating from hyperons is dominant only for very massive neutron stars having masses 1.8 to 2.0 $\, M_{\odot}$, a successful detection of this phase shift during the inspiral of binary systems with relatively low masses of 1.4 to 1.6 $\, M_{\odot}$ could be a smoking gun signature for the existence of strange quark stars.

We explore four-dimensional scalar-tensor theories obtained from well-defined dimensional regularizations of Lovelock invariants. When an infinite tower of corrections is considered, these theories allow for cosmological models in which the Big Bang singularity is replaced by an inflationary phase in the early-universe, and they also admit a specific class of regular black hole solutions.

We clarify how the weak-interaction-driven bulk viscosity $\zeta$ and the bulk relaxation time $\tau_\Pi$ of neutrino-transparent $npe$ matter depend on the nuclear symmetry energy. We show that, at saturation density, the equation-of-state dependence of these transport quantities is fully determined by the experimentally constrained nuclear symmetry energy $S$ and its slope $L$. Variations of $L$ within current experimental uncertainties can change the (frequency-independent) bulk viscosity by orders of magnitude. This suggests that dissipative effects encoded in the gravitational-wave signatures of binary neutron star inspirals may help constrain nuclear symmetry energy properties.

We study spin oscillations of neutrinos in relativistic moving matter inside an accretion disk. These neutrinos are gravitationally scattered off a spinning Kerr black hole surrounded by a thick accretion disk. The disk can co-rotate and counter-rotate with respect to BH spin. We perform numerical simulations of the propagation of a large number of incoming test neutrinos. We briefly discuss our results.

This work conducts an in-depth exploration of exact electrically charged solutions, including traversable wormholes, black holes, and black bounces, within the framework of the scalar-tensor representation of hybrid metric-Palatini gravity (HMPG) with a non-zero scalar potential. By integrating principles from both the metric and Palatini formulations, HMPG provides a flexible approach to addressing persistent challenges in General Relativity (GR), such as the late-time cosmic acceleration and the nature of dark matter. Under the assumption of spherical symmetry, we employ an inverse problem technique to derive exact solutions in both the Jordan and Einstein conformal frames. This method naturally leads to configurations involving either canonical or phantom scalar fields. A thorough examination of horizon structures, throat conditions, asymptotic behaviour, and curvature regularity (via the Kretschmann scalar) reveals the intricate causal structures permitted by this theoretical model. The analysis uncovers a diverse range of geometric configurations, with the phantom sector exhibiting a notably richer spectrum of solutions than the canonical case. These solutions encompass traversable wormholes, black universe models, where the interior of a black hole evolves into an expanding cosmological phase rather than a singularity, as well as black bounce structures and multi-horizon black holes. The results demonstrate that introducing a non-zero scalar potential within HMPG significantly expands the array of possible gravitational solutions, yielding complex causal and curvature properties that go beyond standard GR. Consequently, HMPG stands out as a powerful theoretical framework for modelling extreme astrophysical environments, where deviations from classical gravity are expected to play a crucial role.

We present a formula for the universal anomalous scaling of the multipole moments of a generic gravitating source in classical general relativity. We derive this formula in two independent ways using effective field theory methods. First, we use the absorption of low frequency gravitational waves by a black hole to identify the total multipole scaling dimension as the renormalized angular momentum of black hole perturbation theory. More generally, we show that the anomalous dimension is determined by phase shifts of gravitational waves elastically scattering off generic source multipole moments, which reproduces the renormalized angular momentum in the particular case of black holes. The effective field theory approach thus clarifies the role of the renormalized angular momentum in the multipole expansion. The universality of the point-particle effective description of compact gravitating systems further allows us to extract the universal part of the anomalous dimension, which is the same for any object, including black holes, neutron stars, and binary systems. As an application, we propose a novel resummation of the universal short-distance logarithms (``tails'') in the gravitational waveform of binary systems, which may improve the modeling of signals from current and future gravitational wave experiments.

We study the evolution of domain wall networks and their phenomenological implications in a model of a real scalar $\chi$, where a $Z_2$-symmetry is slightly broken by a potential bias $V_{bias}$. It is demonstrated that the latter triggers domain wall annihilation considerably earlier than previously thought. Namely, we observe that the scaling relation $t_{ann} \propto 1/V^{2/3}_{bias}$ for the annihilation time $t_{ann}$ fits to the simulation data better than a commonly assumed $t_{ann} \propto 1/V_{bias}$. As a result, the energy density of gravitational waves produced by the network of biased domain walls, for a given tiny $V_{bias}$, is suppressed compared to naive expectations. The spectral shape of gravitational waves is similar to that resulting from unbiased domain walls, but with more power in the close-to-maximum ultraviolet part. In the far ultraviolet region, the spectrum of gravitational waves becomes nearly flat; such a plateau has been recognized earlier in the case of unbiased walls. In our investigation we mainly focus on the symmetry breaking potential $V_{breaking} \propto \chi^3$, and argue that no significant modifications of the domain walls evolution take place if one includes higher powers of $\chi$.

We investigate the deflection of photons in the strong deflection limit within static, axisymmetric spacetimes possessing reflection symmetry. As the impact parameter approaches its critical value, the deflection angle exhibits a logarithmic divergence. This divergence is characterized by a logarithmic rate and a constant offset, which we express in terms of coordinate-invariant curvature evaluated at the unstable photon circular orbit. The curvature contribution is encoded in the electric part of the Weyl tensor, reflecting tidal effects, and the matter contribution is encoded in the Einstein tensor, capturing the influence of local energy and pressure. We also express these coefficients using Newman--Penrose scalars. By exploiting the relationship between the strong deflection limit and quasinormal modes, we derive a new expression for the quasinormal mode frequency in the eikonal limit in terms of the curvature scalars. Our results provide a unified and coordinate-invariant framework that connects observable lensing features and quasinormal modes to the local geometry and matter distribution near compact objects.