We study black hole (BH) formation from the nonlinear growth and collapse of primordial perturbations during the matter-dominated era. Modelling cold dark matter (CDM) as pressureless dust, we describe the collapse in a fully nonlinear relativistic framework using the Lemaître-Tolman-Bondi (LTB) and quasi-spherical Szekeres solutions as exact perturbations of a spatially-flat Friedmann-Lemaître-Robertson-Walker (FLRW) $Λ$CDM background. At first order in relativistic scalar perturbation theory, the growing mode of any relevant quantity can be expressed in terms of the conserved gauge-invariant curvature perturbation $\mathcal{R}_c$, which acts as a potential for the 3-curvature of hypersurfaces orthogonal to the matter 4-velocity. We use this result to express the active gravitational mass and curvature functions of the LTB and Szekeres models in terms of the initial values of $\mathcal{R}_c$ and its spatial derivatives. From these initial curvature data we derive: (i) the turn-around, collapse, and apparent-horizon formation times, and (ii) the regularity conditions required for BH formation. We show that sinusoidal and Gaussian profiles do not provide viable BH-forming channels, whereas broad compensated curvature peaks, naturally predicted by peak theory, do. We then estimate the formation times of $10^{3}-10^{6}~\mathrm{M}_\odot$ massive BH seeds produced by the direct collapse of primordial CDM curvature peaks, finding full BH formation at redshifts $z>5$, with core collapse beginning at $10 \lesssim z \lesssim 16$. Finally, we characterize the local dynamics and singularity type of the collapse (point-like, cigar-like, or pancake-like) directly from the initial comoving curvature data, clarifying the role of the initial shear in selecting the collapse end-state. less

We extend the frequency-domain analysis of quasinormal modes in a dynamical, spherically symmetric black hole spacetime undergoing constant-rate mass evolution. In particular, we report a novel feature of the spectrum: the presence of purely imaginary eigenvalues in addition to the usual light-ring modes. We study the frequencies of these modes both analytically and numerically. The analytical calculation uses a novel formalism based on recent advances in connection coefficients of Heun functions. We then compute the frequencies numerically using a spectral method on hyperboloidal slices and find excellent agreement between the two approaches. Finally, we validate the frequency-domain results against an independent set of time-domain simulations. Our analysis shows that the purely imaginary modes govern the late-time signal through exponentially decaying tails. In the Schwarzschild limit, both frequency- and time-domain studies consistently show that the purely imaginary modes give rise to the familiar Schwarzschild power-law tail. less

The possibility that Dark Matter (DM) is partially or totally constituted by stable Planckian remnants of light Primordial Black Holes (PBHs), suggested for instance by Loop Quantum Gravity (LQG), is investigated. Distinct phenomenological regimes are identified, including scenarios that trigger an early matter-dominated epoch. New constraints are derived on the initial PBH and final remnant abundances. We show that a significant initial abundance of PBHs lighter than $10^3$ kg would overproduce Planckian relics, implying that any observational evidence for such PBHs would challenge models with quasi-stable remnants. Conversely, the products of Hawking radiation from PBHs with masses between $10^3$ and $10^{12}$ kg impose that Planckian relics could only be a highly subdominant DM component. We identify a PBH mass around $10^3$ kg for which Hawking evaporation naturally reheats the Universe while the remnants entirely constitute the present-day DM. Such a scenario does not require fine-tuning the initial abundance of PBH of this mass, which could range from $10^{-10}$ to order one. These early-Universe cosmologies yield distinct observational signatures: scalar-induced gravitational waves sourced by primordial or Poisson fluctuations that are amplified by the early PBH-dominated era. Current and future observations of LIGO/Virgo/KAGRA, the Einstein Telescope and LISA, as well as probes of the effective number of relativistic degrees of freedom, can be used to probe and constrain the initial PBH abundance and the present-day abundance of Planckian relics. less

Non-singular gravitational theories are expected to be relevant in the early universe. In this paper, we derive a set of effective Friedmann equations describing the dynamics of matter shells in the presence of a gravitational regulator $\ell$. We find that such a regulator induces a primordial black hole mass gap such that below a certain mass $M_\text{gap}(\ell, R_H)$ no black holes can form. The order of magnitude of this mass gap is set by the regulator $\sim c^2\ell/G$, with subleading dependence on the horizon radius at time of formation $R_H$. Finally, we show that over a wide range of equation of state parameters $ω= 0 \dots 1/3$, the mass gap implies a Carr criterion of the form $δ_H > 2G M_\text{gap}/R_H - 1$. If the horizon size is of the same order of the regulator, $R_H \sim \ell$, this new criterion is stronger than the traditional Carr criterion for primordial black hole formation. This connects the primordial black hole abundance directly to the presence of gravitational regulators. less

We investigate thin-disk imaging of Kerr black holes with synchronized scalar hair, focusing on backreacted parity-odd excited states of a complex scalar field minimally coupled to Einstein gravity. The spacetime displays a core-double-torus lensing structure, with a central black hole surrounded by two scalar clouds. We study the dependence of the images on hair strength and viewing angle, identifying a weak-hair regime close to Kerr. With increasing hair, the photon ring and shadow region shrink and become more distorted. In the strong-hair regime, gravitational lensing produces new features, including multiple disconnected shadow components, crescent-shaped structures, and signatures of chaotic lensing. For nearly edge-on viewing angles, repeated equatorial crossings generate nested ring-like patterns. These results highlight possible geometric signatures of black holes with excited scalar hair. less

A new two-parameter asymptotically flat (AF) toric gravitational instanton is identified as a special case of the Euclidean double Kerr-NUT solution, by imposing certain symmetry and regularity conditions on its rod structure. These conditions are solved explicitly, except for one which takes the form of a fifth-order polynomial. This gravitational instanton has Euler number $χ=4$ and Hirzebruch signature $τ=0$, and its global topology is $\mathbb{C}P^2\#\overline{\mathbb{C}P^2}$ with a circle $S^1$ removed appropriately. It is the third of an infinite sequence of AF toric gravitational instantons that was proved to exist by Li and Sun, the first two being the Kerr and Chen--Teo instantons. It is also the first known example of a Ricci-flat gravitational instanton that is not Hermitian. less

With the multitude of gravitational wave observations that have been made in the past ten years, probing the dynamical and nonlinear nature of strong gravity is becoming more and more feasible. One promising way to test the nonlinear nature of Einstein's theory of general relativity (GR) is through the gravitational wave null memory effect: a nonlinear prediction of GR which corresponds to initially comoving observers being permanently displaced due to a burst of gravitational radiation. Previous studies have shown that, while it is unlikely that the memory effect will be observed in a single event by the LIGO-Virgo-KAGRA (LVK) detectors, evidence for memory in the population of LVK events should be attainable after $\sim$2,000 gravitational wave detections. These works, however, largely relied on Bayes factors to perform their memory analyses: an approach that can depend sensitively on the analysis priors and, when naively multiplied across many events, can even favor incorrect conclusions. In this work, using the GWTC-4.0 catalog of binary black hole observations, we instead perform hierarchical Bayesian inference -- which is not subject to the issues associated with Bayes factors -- to measure the evidence for memory in current LVK observations. We find that we can constrain what we call the memory enhancement factor -- the constant appearing in front of the contribution to the strain from the supermomentum flux -- to $0.32^{+6.30}_{-5.12}$ (with $\pm$ values denoting the 68% credible interval), consistent with its GR value of 1. We forecast that $\sim$2,500 detections will be needed to constrain the memory enhancement factor away from zero at the $1σ$ level. less

The configuration space of causal sets is vast. It is a critical goal to map out this space. Here, we take a practical step towards this goal. We investigate nine classes of causal sets, most of them not studied before. These include manifoldlike causal sets with inhomogeneous Ricci curvature, both topologically trivial and nontrivial. We also study classes of non-manifoldlike causal sets, including lattices, layered orders as well as Lorentzian quasicrystals. Finally, we study classes of causal sets that are not manifoldlike, but are expected to become manifoldlike under a suitable coarse-graining process. We use this broad range of distinct classes of causal sets as a testbed for observables. Rather than focusing on continuum-geometry inspired observables, such as curvature invariants, which often exhibit large fluctuations and are computationally very expensive, we focus on graph observables, including some observables that constitute subgraph statistics and some that are global. We find that three observables, namely the link degree distribution, the eigenvalues of the graph Laplacian of the symmetrized Hasse diagram and the recently proposed abundance of causal intervals, can distinguish between the distinct classes of causal sets. This is made possible by the small fluctuations that these observables have in most classes. less

Subspace diagonalization techniques based on quantum sampling, such as quantum selected configuration interaction (QSCI) and sample-based quantum diagonalization (SQD), have recently emerged as promising quantum-centric approaches for approximating ground-state energies of many-body systems. However, their performance is fundamentally limited by an intrinsic trade-off between sampling efficiency and the sparsity of the ground-state wavefunction, which becomes particularly severe in strongly correlated systems. Here, we introduce a filter-assisted SQD protocol that engineers wavefunction sparsity via a quantum filter, i.e., a unitary transformation of the Hamiltonian designed to concentrate the ground-state weight onto a small number of computational basis states. Using the Gini coefficient as a robust sparsity measure, we establish a quantitative relationship between wavefunction sparsity and the resource requirements of SQD, providing theoretical bounds on the required subspace dimension and sampling cost. To realize the quantum filter, we employ a tensor-network-based circuit-encoding algorithm that maps target states to quantum circuits with controllable fidelity. We benchmark our approach on the quantum Ising model with transverse and longitudinal fields using both numerical simulations and quantum hardware experiments. Our results demonstrate that, compared with standard SQD, the proposed protocol significantly enhances wavefunction sparsity, reduces ground-state energy estimation errors by orders of magnitude, and substantially lowers sampling overhead. These findings establish filter-assisted subspace diagonalization as a powerful and scalable framework for quantum many-body calculations in the strongly correlated regime. less

We construct a detailed microscopic noise model for multi-qubit (MQ) gate operations in the context of trapped ion architecture with all-to-all connectivity. We find that phonon heating and motional dephasing are well captured by effective single- and two-qubit error channels that can, in principle, act between arbitrary pairs of qubits. Nevertheless, the median magnitude of two-qubit errors between uncoupled qubits is substantially smaller than that of errors between gate-coupled qubits. Errors associated with photon scattering are shown to solely propagate to qubits participating in gate operations. Lastly, we combine all noise sources, assigned with experimentally relevant parameters, and explore the scalability of a quantum error correction (QEC) scheme based on the rotated surface code, as a function of error rates and code size. Our analysis bridges device-level physics and QEC performance for MQ gates in trapped-ion architectures. less

Distinct Clifford orbits of magic states can exhibit different stabilizer ranks at small tensor powers. We establish this for qutrits, where the single-qutrit Clifford group has four inequivalent orbits of magic states: Strange, Norrell, Hadamard-eigenstate, and the qutrit T-state, but a nontrivial upper bound on the asymptotic exponent had been pinned down for only the qutrit T-state. For the other three orbits we give explicit stabilizer decompositions, yielding upper bounds on the per-copy asymptotic stabilizer-rank exponent: $γ_S \le \log_3(2)/2 \approx 0.316$ for the Strange state, and $γ_{H_3}, γ_N \le \log_3(4)/3 \approx 0.421$ for the Hadamard-eigenstate and Norrell orbits, all strictly below the prior $γ_{T_3} \le 1/2$ baseline. We also prove the first nontrivial $Ω(m / \log m)$ asymptotic lower bounds for the Hadamard-eigenstate and Norrell orbits, and exhibit two-qutrit Clifford circuits that convert two copies of these states into an injectable phase state with constant success probability, enabling constant-overhead injection of one non-Clifford diagonal gate per orbit. In the case of qubits, we give a closed-form decomposition of the qubit T-type orbit at four copies matching the existing $γ_T \le \log_2(3)/4 \approx 0.396$ exponent via a direct algebraic identity rather than an entangled cat-state construction. An open-source library stabrank accompanies the paper, with Lean 4 proof formalizations of all the decompositions. less