We study physical (near-kernel of constraints) states of 4-d Euclidean loop quantum gravity in Smolin's weak coupling limit on the complete graph $K_5$ using variational Monte Carlo with neural network quantum states. We investigate the Hamilton constraint $\hat{H}$ in the ordering proposed by Thiemann, as well as $\hat{H}^\dagger$ and $\hat{H}+\hat{H}^\dagger$. We find that the variational optimisation selects distinct solution families for $\hat{H}$ and $\hat{H}^\dagger$ across several considered cutoffs on the kinematical degrees of freedom. The solution family of $\hat{H}$ is flat on all minimal loops and has non-vanishing volume expectation values. Its edge-charge marginals delocalise with increasing cutoff, which indicates they are approximations to solutions that are non-normalisable in the kinematical inner product. The solution family for $\hat{H}^\dagger$ is normalisable, shows non-trivial charge correlations, lies in the kernel of volume and is not flat. $\hat{H}+\hat{H}^\dagger$ turns out to be much harder to solve and yields quasi-solutions combining features of both previous families. We characterise all solutions using chromaticity 1- and 2-point functions, minimal loop holonomies, geometric area and volume observables and show that the two families can be interpreted as, on the one hand, a family of states close to the Ashtekar-Lewandowski vacuum and the Dittrich-Geiller vacuum with some numerical noise on the other hand. We also present some results that link solutions of the truncated theory to solutions of the continuum theory. less

The oscillation modes of stars play an important role in observations, and on the understanding of stellar stability properties. The role of viscosity in the oscillation modes of compact stars has been so far understood very loosely only, in absence of a well posed framework. We use recent breakthroughs in the formulation of a causal and stable theory of relativistic hydrodynamics, to study oscillation modes of neutron stars. We characterize the axial spectrum of compact stars and uncover new, viscosity-driven families of modes, without a perfect fluid counterpart. Our results show mode avoidance in some of these families, and a spectrum of long-lived modes, whose role in astrophysical, dynamical processes is yet to be understood. less

Black hole shadow images are primarily determined by the properties of photon spheres and can exhibit degeneracies across different spherically symmetric spacetime geometries. We show that time delay observables associated with higher-order images of transient sources provide a robust probe to break such degeneracies in spacetimes admitting multiple photon spheres. Adopting a model-independent, parametrized, static, spherically symmetric framework that captures the generic features of double-peaked effective potentials, we investigate photon geodesics and quantify them in terms of angular deflection, travel time, and the order of the image. We identify distinctive signatures of trajectories probing the region between the unstable photon spheres. In particular, we find that these trajectories are characterized by the nontrivial temporal behavior, including a minimum travel time, a minimum angular deflection, and a characteristic triplet structure of higher-order images with a specific arrival sequence. We further show that the influence of the depth of the potential well, between the two photon spheres, on the observed time delays provides a direct handle on otherwise inaccessible regions of the spacetime. Our results highlight that time-domain lensing observables encode information beyond static shadow images and offer a promising avenue for probing the structure of compact objects and the strong-field regime of gravity. less

Multiple potential wells for massive test particles, allowing distinct families of bound orbits to coexist, are a characteristic feature of certain exotic compact objects beyond general relativity. Taking the dyonic black hole as a representative example, we demonstrate that such multi-well geometries generically support multiple coexisting branches of bound orbits, in contrast to the single-branch behavior observed in the Schwarzschild spacetime. Crucially, the periodic orbits sharing identical rational rotation number, and hence identical topological indices can nevertheless produce \emph{radiatively distinct} gravitational waves in a representative extreme-mass-ratio inspirals: their amplitude modulation and harmonic content differ because each branch spans different regions of spacetime curvature. These ``topologically equivalent yet waveform-distinguishable'' signatures provide a direct observational probe of strong field gravitational dynamics beyond general relativity, potentially accessible to future space-based gravitational wave detectors. less

Transient noise artifacts, commonly referred to as glitches, pose a major challenge to parameter inference for space-based gravitational-wave (GW) observations. We develop a glitch-robust amortized inference framework for massive black hole binaries in the Taiji detector configuration by combining conditional normalizing flows, a time-frequency multimodal fusion encoder, and contrastive learning. To enable large-scale training on contaminated data, we further introduce a neural glitch generator that produces high-fidelity synthetic transients at substantially reduced computational cost. Systematic experiments show that, under glitch contamination, the proposed method yields more accurate and better-calibrated posteriors than a conventional Markov Chain Monte Carlo baseline. In ablation studies, the full time-frequency model with contrastive learning performs best overall and remains robust to variations in glitch duration and merger-relative timing. We further show that standard coverage diagnostics alone are insufficient to fully assess posterior fidelity. We therefore complement them with the continuous ranked probability score, which provides a stricter assessment of global distributional agreement in non-ideal GW data. Taken together, these results establish deep-learning-based amortized inference as a promising framework for fast and robust Bayesian parameter estimation in future space-based GW observations. less

We study the Tolman-Oppenheimer-Volkoff equation in the presence of a cosmological constant for general thermodynamically consistent equations of state, without imposing regularity at the center. Formulating the problem as an initial value system integrated from an outer boundary inwards, we obtain a general classification of solutions and show that singular configurations dominate the solution space. We demonstrate that all singular solutions share a universal geometric structure and give rise to spacetimes that are bounded-acceleration complete, indicating that the associated singularities are comparatively mild. Our results extend the classification previously obtained for Λ=0 and reveal qualitatively new features for $Λ\neq 0$. For $Λ< 0$, we identify solutions with approximate horizon structures that mimic black holes in equilibrium with their Hawking radiation. For $Λ> 0$, we find four distinct classes of solutions with cosmological horizons, distinguished by the behavior of their temperature gradients. less

We investigate the stability of ultra-compact stellar configurations in the context of an interacting vacuum component. By extending the Tolman-Oppenheimer-Volkoff equations to include a covariant energy exchange between the fluid and vacuum sectors, we examine how the classical Buchdahl stability limit is modified. We analyze two phenomenological interaction models: a coupling to the matter energy density gradient and a direct coupling to the spacetime curvature. Numerical integration reveals that while standard General Relativity predicts a central pressure divergence as the compactness approaches the Buchdahl threshold, the interaction term $Q_ν$ relaxes the pressure gradient and maintains a finite, well-behaved central pressure for proper domains of the coupling parameter. These results demonstrate that an interacting vacuum provides a physical mechanism to bypass classical geometric bounds, potentially supporting ultra-compact objects in regimes previously considered singular. less

The emergence of precision gravity simulators in quantum and fluid systems is opening new avenues for probing curved-spacetime physics and black-hole phenomenology under controlled laboratory conditions. In parallel, advances in understanding how fundamental physics can be probed in the spectral signatures of black holes and exotic compact objects motivate the development of modern spectroscopic techniques within analogue-gravity experiments. In this work, we model the spectral properties of analogue black holes sourced by broadband stochastic noise, a crucial aspect in realistic experiments that poses substantial challenges for established data-analysis techniques. Using simulation-based inference, we demonstrate that the physical parameters encoded in noisy spectra can be reliably extracted, showing that these techniques provide a powerful tool for studying both spacetime properties and boundary effects in gravity simulators. less

We introduce pyEFPEHM, a post-Newtonian (PN) inspiral waveform model for eccentric and spin-precessing compact binaries that includes higher-order modes and matter effects. Accurate and efficient waveform models capturing these effects are essential for probing compact-binary formation channels and exploiting current and future gravitational-wave (GW) observations. pyEFPEHM extends pyEFPE, significantly improving its physical content and accuracy. In particular, we show that above 2.5PN order the quasi-circular contributions to the orbital phasing dominate at each PN order, and incorporate all available higher-order quasi-circular PN corrections to the phasing, including adiabatic tidal effects. We generalize the multiple-scale analysis solution of the spin-precession equations, extending it to higher PN orders and including all available quasi-circular corrections. Finally, we add eccentric corrections up to 1PN order in the waveform amplitudes, including the GW multipoles $(l,|m|)=(2,2),(2,1),(2,0),(3,3),(3,2),(3,1),(3,0),(4,4),(4,2),(4,0)$. We validate pyEFPEHM against analytical waveform models and numerical relativity simulations, showing that it provides a robust and computationally efficient description of the inspiral, with good agreement across a broad region of parameter space and up to close to merger. The accuracy degrades in the late inspiral for systems with very unequal masses ($m_2/m_1 \lesssim 0.1$), significant spins aligned with the orbital angular momentum ($|χ_\mathrm{eff}| \gtrsim 0.5$), and high eccentricities ($e \gtrsim 0.6$), where the PN expansion is expected to break down. pyEFPEHM represents a significant step toward physically complete and efficient waveform modeling of eccentric and precessing binaries, providing a foundation for future extensions including higher-order corrections, calibration to numerical relativity, and merger ringdown modeling. less

We investigate the cosmological implications of non-polynomial quasi-topological gravity (NPQTG), a novel class of modified gravitational theories in which the background dynamics is encoded in a single function of the Hubble parameter. This framework provides a minimal and theoretically consistent extension of general relativity, incorporating higher-curvature effects while preserving second-order field equations and avoiding higher-derivative instabilities. We first establish the general conditions for cosmological viability and construct explicit realizations, including polynomial, quartic, power-law and non-polynomial models, demonstrating how different functional forms lead to distinct expansion histories. Focusing on the quartic and power-law cases, we show that the resulting cosmological evolution reproduces the standard thermal history of the Universe and gives rise to an effective dark-energy sector of geometric origin, with dynamical equation-of-state behavior that can lie in the quintessence or phantom regime. We then confront the models with observational data from Type Ia Supernovae, Cosmic Chronometers, and Baryon Acoustic Oscillations, using a Bayesian MCMC analysis. We find that both models provide an excellent fit to the data, remaining fully compatible with current constraints and statistically competitive with $Λ$CDM. Our results demonstrate that NPQTG offers a simple and efficient framework for describing late-time cosmic acceleration with dynamical dark energy, while maintaining theoretical consistency and observational viability. less