Recent work has shown that introducing highercurvature terms to the EinsteinHilbert action causes the approach to a spacelike singularity to unfold as a sequence of Kasner eons. Each eon is dominated by emergent physics at an energy scale controlled by highercurvature terms of a given order transitioning to higherorder eons as the singularity is approached. The purpose of this paper is twofold. First we demonstrate that the inclusion of matter dramatically modifies the physics of eons compared to the vacuum case. We illustrate this by considering a family of quasitopological gravities of arbitrary order minimally coupled to a scalar field. Second we investigate Kasner eons in the interior of black holes with field theory duals and analyze their imprints on holographic observables. We show that the behavior of the thermal afunction twopoint functions of heavy operators and holographic complexity can capture distinct signatures of the eons making them promising tools for diagnosing stringy effects near black hole singularities.
['High Energy Physics - Theory', 'General Relativity and Quantum Cosmology']
Kasner eons with matter holographic excursions to the black hole singularity
{"Elena C'aceres, a ' Angel J. Murcia, b Ayan K. Patra, c Juan F. Pedraza c": "- a Theory Group, Department of Physics, University of Texas, Austin, TX 78712, USA\n- b INFN, Sezione di Padova, Via Francesco Marzolo 8, I-35131 Padova, Italy\n- c Instituto de F'ısica Te'orica UAM/CSIC, Calle Nicol'as Cabrera 13-15, Madrid, E-28049, Spain  \nE-mail: elenac@utexas.edu , angel.murcia@pd.infn.it , a.patra@csic.es ,  \n- j.pedraza@csic.es  \nAbstract: Recent work has shown that introducing higher-curvature terms to the EinsteinHilbert action causes the approach to a space-like singularity to unfold as a sequence of Kasner eons . Each eon is dominated by emergent physics at an energy scale controlled by higher-curvature terms of a given order, transitioning to higher-order eons as the singularity is approached. The purpose of this paper is twofold. First, we demonstrate that the inclusion of matter dramatically modifies the physics of eons compared to the vacuum case. We illustrate this by considering a family of quasi-topological gravities of arbitrary order minimally coupled to a scalar field. Second, we investigate Kasner eons in the interior of black holes with field theory duals and analyze their imprints on holographic observables. We show that the behavior of the thermal a -function, two-point functions of heavy operators, and holographic complexity can capture distinct signatures of the eons, making them promising tools for diagnosing stringy effects near black hole singularities.", '1 Introduction': "General relativity has long been known to predict a complex but manageable behavior near space-like singularities, particularly in the Belinski-Khalatnikov-Lifshitz (BKL) limit [1-3], where spatial derivatives become negligible and ultra-local 'Carrollian' dynamics emerge [4]. The Kasner metric, a key solution in this limit, describes a homogeneous yet anisotropically expanding spacetime with fixed scaling exponents, constrained by specific relations to satisfy Einstein's equations. Under certain simplifying assumptions, this metric is conjectured to characterize the 'final regime' of evolution before general relativity loses its predictive power.  \nMore generally, the BKL limit predicts a sequence of Kasner epochs -periods where the spacetime metric closely follows a Kasner geometry- connected by brief transitions. The dynamics of the BKL limit are particularly intriguing when the scaling exponents evolve through a space of solutions defined by the Kasner relations. This evolution can be influenced by factors like spatial curvature or couplings to specific matter fields. For example, in the context of supergravity, it has been shown that the inclusion of p -forms triggers 'mixmaster' dynamics governed by affine Lie algebras [5]. In this scenario, the evolution of Kasner exponents can be understood as a particle moving on a hyperbolic billiard table, giving rise to the concept of 'cosmological billiards.' The boundaries of this table, which define the symmetry algebra, are determined by the particular setup under consideration.  \nThe long-term evolution of the BKL limit reveals two key stages: Kasner eras and Kasner eons. A Kasner era comprises multiple Kasner epochs, where two spatial directions alternately contract and expand while the others remain fixed [6]. However, numerical evolution through these eras is not so straightforward, due to the emergence of 'spikes,' spatial structures that disrupt expected asymptotic locality [7-9]. Recently, higher-curvature 'stringy' corrections to the Einstein-Hilbert action have been shown to introduce new transitions, causing the approach to a space-like singularity to unfold through a sequence of Kasner eons -extended periods dominated by physics at energy scales defined by the higher-curvature terms [10]. As the singularity nears, the evolution progresses through eons of increasingly higher order.  \nThe analysis in [10] focused on specific higher-order gravities in vacuum. While the emergence of eons is likely a robust feature arising from emergent stringy effects, it is essential to recognize the critical role that matter fields play in any BKL analysis. As space-like singularities are approached, matter fields may undergo infinite growth, resulting in substantial backreaction on the geometry. This instability typically modifies the Kasner relations, which can vary significantly depending on the matter content of the theory [1-3]. Therefore, incorporating matter is crucial for understanding the dynamics near space-like singularities, and should be essential for fully uncovering the physics of Kasner eons. To address this issue, in the first part of this paper we will investigate a family of quasi-topological gravities of arbitrary order coupled to a scalar field, illustrating how higher-order gravitational effects interact with matter fields and influence behavior near singularities. 1 These findings will offer valuable insights into Kasner eons from a purely gravitational perspective.  \nIn the second part of this paper, we shift our focus to Kasner eons within the context of the AdS/CFT correspondence. It is well-known that the singularity inside an AdS-Schwarzschild black hole exhibits Kasner-like behavior. Recently, Frenkel and collaborators [11] studied a class of AdS black holes with scalar hair, uncovering a deformation of the singularity into a more general Kasner universe, consistent with BKL expectations. Following [11], several studies have studied the dynamics of Kasner singularities from a holographic perspective, including scenarios with cosmological bounces and Kasner inversions akin to those observed in mixmaster universes [12-39]. Within the AdS/CFT correspondence, these phenomena are interpreted in terms of Renormalization Group (RG) flows of the dual CFT [22], where various holographic observables capture essential information about the near-singularity geometries. Based on these antecedents, we thus anticipate that we can engineer specific observables to identify the presence of eons from the field theory perspective.  \nThis paper is organized as follows. In section 2 we revisit the concept of Kasner eons, illustrating their emergence in a specific family of quasi-topological gravities that include terms of arbitrary order in curvature and an arbitrary cosmological constant. We investigate Kasner eons in the vacuum of these theories and in the presence of a scalar field, where, as anticipated, we uncover important modifications of the relations that the Kasner exponents  \nmust satisfy. Notably, when matter is included, these theories ensure that any static and plane-symmetric background is fully described by second-order equations of motion, which drastically simplifies our analysis. In section 3 we focus on Kasner eons within black holes with a negative cosmological constant and AdS boundary conditions, which are conducive to a holographic interpretation. We analyze a handful of observables designed to capture distinct signatures of the eons in terms of the dual field theory, including the thermal a -function, two-point functions of heavy operators, and holographic complexity. Finally, we conclude in section 4 with a summary of our findings and potential directions for future research.", '2 Kasner eons in quasi-topological gravities': "In their seminal work in the '70s, Belinski, Khalatnikov and Lifshitz [1-3] showed that the classical approach to a space-like singularity -whether cosmological or within a black holeunfolds as a sequence of Kasner epochs, during which the metric takes the form:  \nd s 2 = -d τ 2 + d ∑ i =1 τ 2 p i d x 2 i , (2.1)  \nwhere we consider a spacetime with D = d + 1 dimensions, and the Kasner exponents p i satisfy the following relations:  \nd ∑ i =1 p i = 1 , d ∑ i =1 p 2 i = 1 . (2.2)  \nThese conditions are needed to fulfill the Einstein equations in vacuum, but are generically altered with the inclusion of matter or higher-curvature terms. Therefore, we will refer to any spacetime of the form (2.1) as Kasner, regardless of whether the Kasner exponents satisfy (2.2). Kasner eras, on the other hand, consist of a sequence of Kasner epochs, during which two spatial directions successively alternate between contracting and expanding.  \nMeanwhile, Kasner eons emerge over longer time scales following the introduction of higher-curvature terms in the gravitational action [10]. In what follows, we will define Kasner eons and demonstrate how a potentially infinite sequence naturally arises as the gravitational action is supplemented with a tower of higher-curvature corrections. First, in subsection 2.1 we will focus on purely gravitational theories without matter, extending the pioneering work of [10] to arbitrary derivative orders and spacetime dimensions D ≥ 5. Next, in subsection 2.2, we will introduce a minimally coupled scalar field and explore how its presence modifies the approach to the singularity by examining the Kasner exponents across different eons.  \nThe computations and derivations in this section are conducted from a purely gravitational perspective. Specifically, Kasner eons are identified by analyzing the bulk equations of motion within the higher-curvature theories under consideration. Later, in section 3, we will consider Kasner eons in the presence of a negative cosmological constant and demonstrate how they can also be probed using CFT observables via the AdS/CFT correspondence.", '2.1 Kasner eons in vacuum': "Let us start by briefly explaining and motivating the emergence of Kasner eons. For now, we will focus on a vacuum scenario and consider a specific theory of gravity in D = d + 1 dimensions that may include an infinite tower of higher-curvature corrections:  \nL = λ 0 ℓ 2 + λ 1 R + ∞ ∑ n =2 λ n ℓ 2( n -1) Z ( n ) . (2.3)  \nHere, ℓ denotes a length scale at which higher-curvature corrections become important, while λ n are dimensionless constants, with λ 0 representing the cosmological constant expressed in units of ℓ . 2 Without loss of generality, we will set λ 1 = 1 unless specified, achieved by rescaling the theory so that λ n → λ n /λ 1 . The terms Z ( n ) refer to specific higher-curvature terms chosen such that the gravitational equations of motion for any static and plane-symmetric background (equivalently, spherically symmetric or hyperbolic symmetric) in vacuum remain first-order in derivatives. These theories constitute a special subset of quasi-topological gravities [40-46], 3 representing the most natural generalization of Lovelock gravities for curvature orders 2 n > D . Interestingly enough, one can prove (following the arguments from [43], although this reference proves something slightly different) that any effective theory of gravity with no derivatives of the curvature may be mapped, order by order in the UV scale ℓ , to a theory of the form given in (2.3). As discussed in [10], examples of these theories -with first-order equations for any static and plane-symmetric background in vacuum- exist at all curvature orders and for dimensions D ≥ 5. For the reader's convenience, we provide in Appendix A specific expressions for Z ( n ) up to fifth order in the curvature for any spacetime dimension D ≥ 5 that satisfies our criteria, and discuss how to construct Z ( n ) for n > 5.  \nThe equations of motion for the theory (2.3) can always be expressed in the following compact form -see, e.g., [47]:  \nP acde R b cde -1 2 g ab L +2 ∇ c ∇ d P acbd = 0 , P abcd ≡ ∂ L ∂R abcd . (2.4)  \nWe will adopt the following static and plane-symmetric ansatz for the spacetime metric:  \nd s 2 = 1 r 2 [ -e -χ ( r ) f ( r )d t 2 + d r 2 f ( r ) +d ⃗x 2 d -1 ] . (2.5)  \nwhere d ⃗x 2 d -1 denotes the ( d -1)-dimensional Euclidean metric. In these coordinates, the singularity is located at r → + ∞ . A key feature of the special higher-curvature theories under consideration (2.3) is that ∇ d P acbd ∣ ∣ χ, f = 0, where ∣ ∣ χ, f denotes evaluation on (2.5).  \nAs a result, the complete set of gravitational equations of motion can be exactly integrated, yielding the following two conditions 4  \n∞ ∑ n =0 λ n ( -1) n f ( r ) n r -d ( d +1 -2 n )( ℓ 2 d ( d +1)) n -1 = C , χ ( r ) = χ 0 . (2.6)  \nHere χ 0 is a constant that can be arbitrarily set to zero through a time reparametrization and C is a constant related to the energy of the spacetime. By solving (2.6), an algebraic equation for f ( r ), the exact solution for the metric can be determined, enabling a detailed analysis of the approach to the singularity -provided it exists. However, an analytical solution for f ( r ) is generally not possible, complicating a direct study of the singularity.  \nThe behavior of the solutions can nevertheless be examined perturbatively as the singularity is approached. In the initial stage, far from the singularity, Einstein gravity dominates. However, as the singularity is approached, higher-curvature corrections become increasingly significant, gradually taking control of the dynamics. First, quadratic-curvature terms surpass the Einstein regime, determining the new Kasner exponents and the corresponding Kasner metric. As one gets even closer, cubic terms dominate, making the quadratic terms negligible and fully controlling the solution. Each phase where a specific curvature order dominates is called a Kasner eon [10], resulting in a sequence of Kasner eons transitioning from one another as the singularity is approached. From (2.6) it follows that, sufficiently far from these transitions, the n -th eon -obtained by setting all λ m = n to zero- is defined by:  \n̸  \nf ( r ) ∼ -r d/n . (2.7)  \nEach Kasner eon can be characterized through the corresponding Ricci scalar derived from (2.5), which diverges as r d/n in the limit r → + ∞ . This indicates that every eon can be associated with a physical curvature singularity, though the severity of the singularity diminishes as n → + ∞ . This observation highlights the necessity of incorporating the complete tower of higher-curvature terms to potentially achieve full singularity resolution [49, 50].  \nLet us now derive the Kasner metric corresponding to the various Kasner eons resulting from (2.6). To do this, we will assume that:  \nf ( r ) = -f n r γ n , (2.8)  \nwhere γ n ≡ d/n and f n > 0. By introducing a new coordinate defined by f n r γ n +2 d τ 2 = d r 2 , (2.5) simplifies to a pure Kasner metric:  \nd s 2 = -d τ 2 + τ 2 p t d t 2 + τ 2 p x d x 2 d -1 , p x = 2 γ n , p t = p x -1 . (2.9)  \nThus, we conclude that the approach to the singularity is governed by a sequence of Kasner eons, characterized by the Kasner exponents satisfying:  \np x = p t +1 = 2 n d . (2.10)  \nThese exponents alter the vacuum Kasner relations of GR (2.2) as follows:  \nd ∑ i =1 p i = 2 n -1 , d ∑ i =1 p 2 i = 4 n ( n -1) d +1 . (2.11)  \nFor pure Lovelock metrics 5 in even dimensions, d +1 = 2 n , the Kasner exponents (2.10) take the values p x = 1, p t = 0. Conversely, in odd dimensions, d = 2 n , they are given by p x = d -1 d , p t = -1 d . In both cases, the Kasner metric exactly satisfies the vacuum equations of motion -cf. [51], where the Kasner solutions for the highest-order non-trivial Lovelock gravity in a given dimension are studied. In fact, we note that (2.9) represents an exact solution of the equations of motion derived from (2.3) for arbitrary n and d , provided λ m = n are set to zero.  \n̸  \nMore generally, we find that the approach to the singularity in our theory of gravity (2.3), which includes potentially an infinite tower of higher-curvature terms, consists of a sequence of Kasner eons (2.9) with Kasner exponents given by (2.10). Specifically, the Kasner exponents p t and p x are shifted by 2 /d as they transition from one Kasner eon to the next one. 6  \nTo identify the distinct Kasner eons encountered during the approach to the singularity, it is useful to introduce the following effective Kasner exponents [10]:  \np eff x = 2 f ( r ) rf ' ( r ) , p eff t = 2 f ( r ) rf ' ( r ) -1 . (2.12)  \nIndeed, whenever f ( r ) corresponds to a proper Kasner eon in a region near the singularity -such that its form is given by (2.8)- the effective exponents p eff x and p eff t are naturally described by (2.10). Consequently, (2.12) serves as a useful indicator for the presence of Kasner eons, identifiable as intervals of the radial variable r where p eff x and p eff t plateau at constant values, which can be easily spotted through a graphical analysis.  \nTo achieve this, the explicit form of f ( r ) must be determined by solving the algebraic equation (2.6). However, beyond cubic order, a closed-form solution for f ( r ) is generally unattainable, necessitating numerical methods. Notably, we find it is numerically more stable to solve the first-order differential equation obtained by differentiating (2.6), which yields:  \n-λ 0 ℓ 2 + ∞ ∑ n =1 λ n ( -1) n ( d +1 -2 n )( dℓ 2 f ( r )( d +1)) n -1 ( nrf ' ( r ) -df ( r ) ) = 0 . (2.13)  \nThe cosmological constant term λ 0 decouples in the near-singularity limit and does not affect the analysis of Kasner eons. However, anticipating the holographic studies of the next section, we set λ 0 to be negative in our numerics, ensuring an asymptotically AdS spacetime. This choice also guarantees the existence of planar black hole solutions, with f ( r ) representing the blackening factor. Without loss of generality, then, we may assume that f ( r ) has one zero,  \n<!-- image -->  \nFigure 1 : Kasner exponents for the higher-curvature gravities (2.3) in d = 4 (left) and d = 5 (right) with specific dimensionless couplings λ n (chosen values: λ 0 = d ( d -1), λ 1 = 1, λ 2 = 10 -8 , λ 3 = -10 -22 , λ 4 = -10 -42 , λ 5 = -10 -67 and λ n> 5 = 0). We have also set r h = ℓ for simplicity. The emergence of Kasner eons is evident, exhibiting the predicted values, although the cubic eon in d = 5 is absent. This absence arises because, in d ≥ 5, cubic quasi-topological gravities behave like cubic Lovelock gravity, which is topological in d = 5.  \n<!-- image -->  \nrepresenting the black hole horizon, which we set to be at some r = r h . Evaluating (2.13) at this point, we find that:  \nλ 0 ( d -1) ℓ 2 + f ' h r h = 0 , (2.14)  \nwhere f ' h = f ' ( r h ). This is equivalent to fixing the black hole temperature T , since 4 πT = -f ' h . For each value of r h and a small δ > 0, we can now proceed with the numerical integration. Starting just inside the horizon at r = r h + δ , we integrate inward towards the singularity; or starting just outside the horizon at r = r h -δ , we integrate outward towards the AdS boundary. Since f ' h is fully fixed by the equations of motion, we can now approximate -with sufficiently good accuracy- the initial value of f ( r ) as f ( r h ± δ ) ≃ ± f ' h δ .  \nIn Figure 1, we present the effective Kasner exponents (2.12) for specific coupling choices λ n in (2.3) across d = 4 and d = 5 spatial dimensions. Each graph clearly displays a sequence of Kasner eons, represented by the flat regions, followed by rapid transitions to the next eon. Additionally, we find that each Kasner eon exhibits the predicted values (2.9) for p t and p x .", '2.2 Kasner eons with matter': "We now aim to investigate how matter fields affect Kasner eons. As an illustrative example to understand the interplay between higher-curvature terms and matter, we will consider a minimally coupled scalar field on top of our gravity theory (2.3):  \nL = λ 0 ℓ 2 + λ 1 R + ∞ ∑ n =2 λ n ℓ 2( n -1) Z ( n ) -1 2 ( ∂ϕ ) 2 -V ( ϕ ) . (2.15)  \nAt this point, it is natural to wonder about including higher-derivative terms in the scalar field ϕ , as well as non-minimal couplings between gravity and the scalar field. However, a technical challenge arises: the resulting equations of motion, even within the highly symmetric configurations described by (2.5), would generally introduce higher derivatives in the fields, complicating the numerical integration. Therefore, in this study, we will limit our analysis to the inclusion of the standard two-derivative kinetic term, which will already introduce significant new physical features, as we will demonstrate below.  \nAfter incorporating the minimally coupled scalar field, the equations of motion for the resulting theory (2.15) are as follows:  \nP acde R b cde -1 2 g ab L +2 ∇ c ∇ d P acbd = 1 2 ∂ a ϕ∂ b ϕ, P abcd ≡ ∂ L ∂R abcd ,  \n(2.16)  \n□ ϕ = V ' ( ϕ ) . (2.17)  \nwhere L now refers to (2.15). As before, we assume a plane-symmetric ansatz for the system, so the metric and the scalar field are written as:  \nd s 2 = 1 r 2 [ -e -χ ( r ) f ( r )d t 2 + 1 f ( r ) d r 2 +d x 2 d -1 ] , ϕ = ϕ ( r ) . (2.18)  \nUnder this assumption, the condition ∇ d P acbd | f,χ = 0 continues to hold. This simplification allows the complete set of equations of motion to be reduced to a system of three ordinary differential equations, each containing no more than two derivatives for f ( r ), χ ( r ) and ϕ ( r ). These equations take the following form:  \nr 2 fϕ '' = 1 2 r ϕ ' ( f ( rχ ' +2( d -1) ) -2 rf ' ) + V ' ( ϕ ( r )) . (2.19)  \nχ ' = r ( ϕ ' ) 2 ∑ ∞ n =1 ˜ λ n ℓ 2( n -1) f n -1 , (2.20)  \nf ' = -2 λ 0 +2 ℓ 2 V ( ϕ ( r )) + 2 d ∑ ∞ n =1 ℓ 2 n ˜ λ n n f n + ℓ 2 fr 2 ( ϕ ' ) 2 2 r ∑ ∞ n =1 ℓ 2 n ˜ λ n f n -1 , (2.21)  \nwhere ˜ λ n are dimensionless constants given by:  \n˜ λ n = ( -1) n +1 n ( d +1 -2 n )( d ( d +1)) n -1 λ n . (2.22)  \nSimilar to the vacuum case, our focus is on determining the Kasner exponents associated with each Kasner eon that emerges as we approach the black hole singularity. To achieve this, we consider solutions in regions where they can be locally approximated by pure Kasner solutions. This approximation holds when f ( r ), χ ( r ) and ϕ ( r ) exhibit the following profiles:  \nf ( r ) = -f n r γ n , χ ( r ) = χ n log r + χ 0 , ϕ ( r ) = { ϕ 1 log r (GR) , ϕ n r ν n (higher-curvature) , (2.23)  \nwhere we explicitly distinguish between the logarithmic profile of the scalar field in Einstein gravity [5] and its power-law behavior in the presence of higher-curvature corrections. In analogy with the vacuum case, we define a new coordinate f n r γ n +2 d τ 2 = d r 2 so that, with the profiles for f ( r ) and χ ( r ) given in (2.23), (2.18) transforms into a pure Kasner metric:  \nd s 2 = -d τ 2 + τ 2 p t d t 2 + τ 2 p x d x 2 d -1 , p x = 2 γ n , p t = p x + χ n γ n -1 . (2.24)  \nNow, observe that any given curvature invariant R ( n ) of order n , when evaluated on a Kasner metric (2.1), diverges in the limit τ → 0 as:  \nR ( n ) ∼ 1 τ 2 n . (2.25)  \nIf we restrict our attention to quadratic potentials of the form V ( ϕ ) = 1 2 m 2 ϕ 2 , then the kinetic term will always diverge faster than V ( ϕ ). Therefore, for the scalar field to influence the Kasner eon, its kinetic terms must match the order of R ( n ) . This implies that the scalar field must diverge as 1 /τ n -1 ∝ r ( n -1) γ n / 2 , so that: 7  \nν n = ( n -1) γ n 2 . (2.26)  \nThe Kasner exponent associated with the scalar, p ϕ , can be naturally defined such that ϕ ( τ ) ∝ τ p ϕ as τ → 0. This definition implies  \np ϕ = -2 γ n ν n = -n +1 . (2.27)  \nNext, let us investigate the values of the remaining Kasner exponents p t and p x in each eon. To do this, consider a regime where the dominant contributions to the equations of motion are the terms of order O ( ℓ 2( n -1) ), for n > 1. This entails imposing (2.23) on the equations of motion and disregarding any subdominant terms. Furthermore, if the potential is quadratic in ϕ , it can be ignored, as its contribution is subleading compared to the kinetic term. Consequently, the scalar equation of motion (2.19) implies that:  \n( n +1) γ n -χ n = 2 d. (2.28)  \nOn the other hand, if we solve for ϕ ' in (2.21), substitute it in (2.20), impose (2.23) and disregard terms arising from curvature invariants of different order than n , we find:  \n-nχ n +2 nγ n = 2 d. (2.29)  \nBy combining the previous two equations, we arrive at:  \nγ n = χ n = 2 d n . (2.30)  \nUsing (2.24), this implies that the exponents p t and p x for the n -th Kasner eon are:  \np t = p x = n d . (2.31)  \nThese Kasner exponents, together with (2.27), modify the Kasner relations as follows: 8  \nd ∑ i =1 p i = n, d ∑ i =1 p 2 i + 1 2 p 2 ϕ = n 2 d + ( n -1) 2 2 . (2.32)  \nLet us now comment on these results. First, note that naively setting n = 1 in (2.32) does not yield the correct relations for Einstein gravity with a minimally coupled scalar field, where ∑ d i =1 p i = 1 and ∑ d i =1 p 2 i + 1 2 p 2 ϕ = 1. This difference is expected, as the scalar field in Einstein gravity exhibits a logarithmic divergence ϕ ( τ ) ∼ p ϕ log τ , in contrast to the power-law behavior in higher-order cases, which results in a distinct definition for p ϕ .  \nSecond, while the Kasner exponents p x and p t were distinct in the absence of the scalar field (2.10), introducing the scalar field dramatically changes the situation. Now, the Kasner exponents for any Kasner eon -except for the unique case of Einstein gravity with matterbecome equal. Interestingly enough, this means that the presence of the scalar field enhances the isotropy of the Kasner solution, effectively rendering the temporal direction (i.e., that associated with the coordinate t ) indistinguishable from any of the spatial directions. One may interpret that matter produces an additional mixing between the Kasner directions that render them identical. The price to pay by that scrambling is that now Kasner exponents are halved compared to the p x exponent in purely gravitational higher-curvature theories. Additionally, the shift δp in Kasner exponents between successive eons is also halved in the presence of the scalar field:  \nδp = 1 d . (2.33)  \nThere is also an important difference between GR and higher-curvature theories when incorporating a minimally coupled scalar field. In the presence of higher-curvature corrections, the Kasner exponents associated with both the spacetime metric and the scalar field are fully determined by (2.27) and (2.31). In contrast, the Kasner exponents are less constrained in Einstein gravity with a minimally coupled scalar field. Specifically, if we examine the equations of motion (2.19)-(2.21) and focus on the Einstein eon, imposing (2.23) led us to:  \nχ 1 = ϕ 2 1 ( d -1) , γ 1 = d + ϕ 2 1 2( d -1) . (2.34)  \nThis implies that the set of Kasner exponents in the GR eon is given by:  \np x = 4( d -1) 2 d ( d -1) + ϕ 2 1 , p t = ϕ 2 1 -2( d -1)( d -2) 2 d ( d -1) + ϕ 2 1 , p ϕ = -p x ϕ 1 . (2.35)  \nNote that the relationships p t + ( d -1) p x = 1 and p 2 t + ( d -1) p 2 x + 1 2 p 2 ϕ = 1 hold true, as expected. Consequently, we find that the Kasner exponents in the GR eon are not fully determined; there is one free parameter that can be identified as either ϕ 1 or p ϕ . In contrast, the set of Kasner exponents for higher-curvature eons is fully determined, which can possibly be attributed to the lack of higher-derivative terms in the scalar field. In this case, we needed to select the appropriate power for the scalar field to ensure it matches the order of the higher-curvature terms -see (2.26)- thereby eliminating the free parameter present in GR.  \nThis concludes the analytical characterization of the various Kasner eons associated with static, plane-symmetric solutions to (2.15). We will now turn our attention to exploring these features through explicit numerical examples.  \nTo begin with, let us provide some details about the numerical treatment of (2.19)-(2.21). As in the vacuum case, we pick a negative λ 0 to guarantee AdS asymptotics and planar black hole solutions, preparing for the holographic studies in section 3. However, we emphasize that the cosmological constant does not play any role in the approach towards the singularity. Now, by assuming regularity at the horizon, r = r h , we may locally expand ϕ , f and χ as:  \nϕ ( r ) = ϕ h + ϕ ' h ( r -r h ) + O [( r -r h ) 2 ] , (2.36)  \nf ( r ) = f ' h ( r -r h ) + O [( r -r h ) 2 ] , (2.37)  \nχ ( r ) = χ h + χ ' h ( r -r h ) + O [( r -r h ) 2 ] , (2.38)  \nwhere the subindex h indicates evaluation at the horizon. Further, by evaluating the field equations at the horizon, we find the following relations between the variables:  \nr 2 h ϕ ' h f ' h = V ' ( ϕ h ) , ( d -1) χ ' h = r h ( ϕ ' h ) 2 , 2 ℓ 2 r h ( d -1) f ' h = -2 λ 0 +2 ℓ 2 V ( ϕ h ) , (2.39)  \nOnce the potential V ( ϕ ) is specified, one can solve for χ ' h , ϕ h and ϕ ' h in terms of r h and f ' h . The latter is related to the black hole temperature T by 4 πT = -f ' h e -χ h / 2 . Notably, the constraints given by (2.39) do not depend on the higher-curvature terms.  \nWe then proceed with the numerical integration by initializing the algorithm at a point r = r h ± δ sufficiently close to the horizon, with δ ≪ 1. Additionally, we set χ h = 0, which can always be accomplished via a time reparametrization. Consequently, by fixing f ' h and r h and specifying a value for δ , we are able to solve for the region between the horizon and the singularity, and between the horizon and the AdS boundary, respectively.  \nTo efficiently visualize the different Kasner eons as we approach the singularity, we now introduce the effective Kasner exponents:  \np eff x = 2 f ( r ) rf ' ( r ) , p eff t = 2 f ( r ) rf ' ( r ) + f ( r ) χ ' ( r ) f ' ( r ) -1 , p eff ϕ ( r ) = -2 χ ' ( r ) ϕ ' ( r ) ϕ ( r ) ( f ' ( r ) f ( r ) ) 2 . (2.40)  \nIndeed, when f ( r ) , χ ( r ) and ϕ ( r ) match a proper Kasner eon near the singularity (2.23), the effective exponents p eff x , p eff t and p eff ϕ are given by (2.31) and (2.27). Furthermore, we explicitly include the factor χ ' ( r ) in p eff ϕ ( r ) since it serves as a clear indicator of the scalar field's presence -whenever χ is constant, leading to χ ' = 0, the scalar field profile becomes  \n-  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n-  \n1.0  \n0  \n5  \n10  \n15  \n20  \n25  \n30  \nFigure 2 : Effective Kasner exponents for higher-curvature gravities (2.3) in d = 4 (left) and d = 5 (right) with the same specific dimensionless couplings λ n used in Figure 1 ( λ 0 = d ( d -1), λ 1 = 1, λ 2 = 10 -8 , λ 3 = -10 -22 , λ 4 = -10 -42 , λ 5 = -10 -67 and λ n> 5 = 0) and m 2 = -10 -4 /ℓ 2 . As reviewed in Appendix B, a negative m 2 is consistent with the Breitenlohner-Freedman stability bound [52, 53]. We have also set r h = ℓ for simplicity. The expected Kasner eons emerge with the predicted values, and the effective Kasner exponent p eff ϕ is rescaled by a factor -1 /d for clarity in the plots.  \n<!-- image -->  \ntrivial, cf. (2.20). These expressions extend the effective Kasner exponents given in [10] to cases where the scalar field is present and χ ( r ) is non-constant.  \nFigure 2 shows the effective Kasner exponents for a specific set of couplings and a massive free scalar field with V ( ϕ ) = 1 2 m 2 ϕ 2 . Notably, the initial Kasner eon is driven by the Einstein gravity eon with a massive scalar field. While this feature cannot be directly observed from Figure 2 -as the GR eon produces a logarithmic profile for the scalar field, which p eff ϕ does not capture- we have confirmed that the numerical solution within the interval defining the GR eon satisfies the conditions (2.35), characteristic of the Einstein gravity eon with a scalar field. Subsequently, the next eon is dominated by a purely gravitational Gauss-Bonnet term, with no significant influence from the scalar field. This behavior suggests that the scalar field requires some time to transition from a logarithmic approach to the singularity to a power-law behavior. Depending on the value of the scalar field at the black hole horizon, this transient period may persist for an extended period, during which the system may only experience transitions between purely gravitational Kasner eons.  \nEventually, the scalar field reaches its characteristic power-law divergence, completing the transition to Kasner eons influenced by the minimally coupled scalar field. In the plots shown in Figure 2, this transition occurs as the purely gravitational quadratic Kasner eon gives way to the next eon. From that point onwards, the scalar field becomes significant enough to impact the Kasner dynamics. The scalar field and higher-curvature terms then become closely intertwined, evolving in accordance with (2.27) and (2.31).", '3 Imprints of Kasner eons on holographic observables': "In this section, we introduce several holographic probes to identify the presence of Kasner eons. Our focus will be on holographic CFTs dual to higher-curvature gravitational theories that include a minimally coupled massive scalar:  \nI = 1 16 πG d +1 ∫ d d +1 x √ -g ( R + d ( d -1) + ∞ ∑ n =2 λ n Z ( n ) -1 2 ( ∂ϕ ) 2 -1 2 m 2 ϕ 2 ) . (3.1)  \nThese theories were examined from a purely gravitational perspective in section 2.2. In fact, the action (3.1) follows from (2.15) with λ 0 = d ( d -1), ℓ = 1 and V ( ϕ ) = 1 2 m 2 ϕ 2 .  \nWe begin by considering the following plane-symmetric ansatz,  \nd s 2 = 1 r 2 ( -f ( r ) e -χ ( r ) d t 2 + d r 2 f ( r ) +d ⃗x 2 d -1 ) , ϕ = ϕ ( r ) . (3.2)  \nwhere d ⃗x 2 d -1 denotes the ( d -1)-dimensional Euclidean metric. In this coordinate system, the conformal boundary is at r = 0, while the singularity lies at r →∞ . The bulk equations of motion, along with their numerical treatment, were thoroughly discussed in section 2.2. It was noted that solving the equations numerically from the horizon towards either the asymptotic boundary or the singularity required setting initial conditions for the scalar field, its derivative, and the metric function χ ( r ), cf. (2.39). In the AdS/CFT correspondence, the scalar field is dual to a scalar operator O with a conformal dimension ∆, satisfying: 9  \nm 2 = ∆(∆ -d ) . (3.3)  \nIn terms of ∆, the initial conditions (2.39) take the form:  \nϕ h = ∓ √ 2 √ d -1 √ -d -f ' h r h √ ∆( d -∆) , ϕ ' h = ± √ 2 √ d -1 √ -d -f ' h r h √ ∆( d -∆) f ' h r 2 h , (3.4)  \nχ ' h = -2( d + f ' h r h ) [∆( d -∆)] f ' 2 h r 3 h . (3.5)  \nIt is important to emphasize that hairy black holes with non-trivial scalar profiles correspond to RG flows in the dual field theory, where the putative CFT is deformed by the operator O . In this context, it is essential to ensure that ∆ < d , so that the dual operator corresponds to a relevant deformation (thereby, affecting the IR of the theory). According to (3.3), this implies we must select a negative m 2 , however, in AdS, there exists an allowable range that complies with the Breitenlohner-Freedman stability bound [52, 53]. 10  \nWith these preliminaries in place, we now turn to examining several holographic quantities that can reveal the presence of Kasner eons from the field theory point of view.", '3.1 Thermal a -function': "In the context of quantum field theories (QFTs), the degrees of freedom decrease due to coarse-graining along the renormalization group flow. The existence of a monotonic function that quantifies the degrees of freedom is well established in two dimensions [54]. This issue has also been explored in higher dimensions [55-59], where the monotonic function is known as the a -function. Via the standard UV/IR relations in holography, it is well known that the radial direction maps to an energy scale in the field theory [60-63]. This indicates that the holographic dual of the a -function should be a bulk quantity that decreases from the boundary (UV) to the horizon (IR). Various proposals for a holographic a -function have been proposed and thoroughly investigated in the literature [64-68]. More recently, a finite temperature generalization of this function, dubbed the thermal a -function, was introduced in [22] to probe and characterize black hole interiors within holographic RG flows. This quantity has been further explored in various contexts in [28, 33, 37, 39, 69].  \nOne effective approach to ensure the monotonicity of the a -function is to require that the matter fields satisfy the Null Energy Condition (NEC),  \nk µ k ν T µν ≥ 0 , (3.6)  \nfor any null vector k α and the corresponding bulk stress-energy tensor T µν . Alternatively, in higher-curvature gravities, we can impose the Null Curvature Condition (NCC), which stipulates that the Ricci curvature of a Lorentzian manifold R µν must not be negative along any arbitrary null vector, 11  \nk µ k ν R µν ≥ 0 . (3.7)  \nThus, (3.6) and (3.7) offer two distinct approaches to defining a monotonic function along the holographic flow. Below, we will propose two monotonic functions corresponding to these conditions and explicitly use them to prove their monotonicity.  \nWorking with the ansatz (3.2) 12 and imposing the NEC (3.6) yields a monotonic function, which we denote as a E T ( r ),  \na E T ( r ) = π d/ 2 Γ ( d 2 ) ( 1 ℓ P ) d -1 exp ( -( d -1) 2 χ -1 2 ∞ ∑ n =2 ˜ λ n ∫ f ( n -1) χ ' dr ) , (3.8)  \nwhere ˜ λ n are the rescaled couplings introduced in (2.22). Alternatively, by imposing the NCC (3.7) we obtain a function that we denote as a C T ( r ),  \na C T ( r ) = π d/ 2 Γ ( d 2 ) ( 1 ℓ P ) d -1 exp ( -( d -1) 2 χ ) . (3.9)  \nMonotonicity of a E T ( r ) and a C T ( t ) can be demonstrated by analyzing their derivative along the entire flow. 13 In particular, we find that  \nda E T dr = -a E T 2 rf ( r ) 2 [ f ( r ) ( T r r -T t t ) ] , (3.10)  \nand  \nrespectively, where  \nda C T dr = -a C T 2 rf ( r ) 2 [ f ( r ) ( R r r -R t t ) ] , (3.11)  \nT r r -T t t = rf ( r ) ( χ ' ( d -1) + ∞ ∑ n =2 ˜ λ n f ( n -1) χ ' ) , (3.12)  \nR r r -R t t = ( d -1) rf ( r ) χ ' . (3.13)  \nOn the other hand, using the null vector ⃗ k = e -χ ( r ) / 2 ∂ t + f ( r ) ∂ r , we can rewrite the NEC and NCC as follows:  \nf ( r ) ( T r r -T t t ) ≥ 0 , f ( r ) ( R r r -R t t ) ≥ 0 , (3.14)  \nimplying da E T /dr ≤ 0 and da C T /dr ≤ 0 everywhere along the flow. Thus, both thermal a -functions monotonically decrease from the AdS boundary to the singularity. We illustrate this behavior in Figures 3a and 3c. 14  \nWe observe several noteworthy properties of these functions. First, when higher-curvature corrections are absent, a E T and a C T coincide. This is expected, as the NEC and NCC are equivalent in Einstein's gravity, which is no longer true when higher-curvature corrections are introduced. Second, both functions effectively capture all matter eons and their transitions through particular combinations of their first and second derivatives. For example, note that  \n∂ z log [ a ' E T ( z ) a E T ( z ) ] = 2 p ' eff ϕ ( z ) p eff ϕ ( z ) + 2( p eff ϕ ( z ) -p eff x ( z ) p ' eff t ( z ) -(1 + p eff t ( z ) -2 p eff x ( z )) p ' eff x ( z ) (1 + p eff t ( z ) -p eff x ( z )) p eff x ( z ) , (3.15) -a '' C T ( z ) a C T ( z ) = ( d -1) ( p eff t ( z ) -p eff x ( z ) + 1 ) 2 -p eff x ( z ) p ' eff t ( z ) + ( p eff t ( z ) + 1 ) p ' eff x ( z ) p eff x ( z ) ( p eff t ( z ) -p eff x ( z ) + 1 ) , (3.16)  \nwhere p eff t ( z ), p eff x ( z ) and p eff ϕ ( z ) are given in (2.40) and we have defined z = log r for convenience. During each Kasner eon, p eff t ( z ), p eff x ( z ) and p eff ϕ ( z ) stabilize at a plateau, indicating that (3.15) and (3.16) must reach constant values as well:  \n∂ z log [ a ' E T ( z ) a E T ( z ) ] → 2 p ϕ p x ( p t -p x +1) , (3.17)  \nFigure 3 : Thermal a -functions in units of π d/ 2 /ℓ d -1 P Γ ( d 2 ) (left) and particular combinations of their first and second derivatives constructed to identify the eons (right) in d = 5 with z = log( r ), λ 2 = 10 -8 , λ 3 = -10 -22 , λ 4 = -10 -42 , λ 5 = -10 -67 , λ n> 5 = 0, and m 2 = -d 2 4 + 1 2 .  \n<!-- image -->  \n-a '' C T ( z ) a C T ( z ) → ( d -1) ( p t -p x +1) p x . (3.18)  \nThis behavior is depicted in Figures 3b and 3d, demonstrating that these monotonic functions serve as a robust diagnostic tool for identifying stringy effects near black hole singularities. Finally, we note that, among the two functions, only a E T defined in (3.8), approaches the zero temperature a -function proposed in [70] for quasi-topological gravities with λ n ≥ 4 = 0. To see this, we switch to the conventional domain wall coordinates,  \nds 2 = e 2 A ( ρ ) ( -h ( ρ ) 2 dt 2 + d⃗x 2 d -1 ) + dρ 2 , (3.19)  \nwhere  \ne 2 A ( ρ ) = 1 r 2 , h ( ρ ) 2 = f ( r ) e -χ ( r ) , dr dρ = -r √ f ( r ) . (3.20)  \nwith ρ ≥ 0. In these coordinates the horizon is situated at ρ = 0, while the conformal boundary is located at ρ → ∞ . This coordinate patch exclusively covers the region outside the black hole. However, as discussed in [22], the black hole's interior can be accessed by  \nanalytically continuing the radial coordinate, i.e., ρ → -iρ . 15 In this patch, our thermal a -function a E T ( ρ ) can be written as follows:  \na E T ( ρ ) = π d/ 2 Γ ( d 2 ) ( 1 ℓ P ) d -1 ( h ( ρ ) A ' ( ρ ) ) d -1 exp [ -∞ ∑ n =2 ˜ λ n ∫ A ' 2( n -1) ( A '' ( ρ ) h ( ρ ) -h ' ( ρ ) A ' ( ρ ) h ( ρ ) A ' ( ρ ) ) dρ ] . (3.21)  \nTo analyze the zero temperature case, we set h = 1. In this case, the integrals in (3.21) can be explicitly carried out, leading to  \na E T = π d/ 2 Γ ( d 2 ) ( 1 ℓ P ) d -1 ( 1 A ' ( ρ ) ) d -1 exp ( -∞ ∑ n =2 ˜ λ n 2 n -2 A ' 2( n -1) ) , (3.22)  \n≈ π d/ 2 Γ ( d 2 ) ( 1 ℓ P ) d -1 ( 1 A ' ( ρ ) ) d -1 ( 1 -∞ ∑ n =2 ˜ λ n 2 n -2 A ' 2( n -1) ) , (3.23)  \nwhere in the last line we only kept the terms that are linear in ˜ λ n ≪ 1. Therefore, our a -function a E T provides a (non-perturbative) generalization of the a -function introduced in [70] for a quasi-topological gravities with an arbitrary number of higher-derivative couplings. 16", '3.2 Two-sided correlators': "In the AdS/CFT correspondence, there is a direct relationship between the length of boundary-anchored geodesics and two-point correlation functions of heavy operators in the boundary CFT [71, 72]. This entry in the holographic dictionary is crucial, as it illustrates how geometric data from the bulk encodes information about the boundary field theory.  \nMore specifically, for two points x 1 and x 2 on the boundary of the AdS, the two-point correlation function G ( x 1 , x 2 ) of a primary operator O ( x ) with large conformal dimension ∆ ≫ 1, is related to the bulk geodesic length L ( x 1 , x 2 ) according to:  \nG ( x 1 , x 2 ) ∼ e -mL ( x 1 ,x 2 ) , (3.24)  \nwhere m ≫ 1 is the mass of the bulk scalar dual to O ( x ). This expression is derived using a 'first quantized' formalism for the correlator, followed by the application of a saddle-point approximation. The geodesic γ connecting the points x 1 and x 2 in the bulk minimizes the distance in the curved AdS space, ensuring it is the shortest path between these points. The length of this geodesic L ( x 1 , x 2 ), which extends deep into the bulk, thus provides a crucial map between geometric data of the AdS bulk and correlation functions in the dual CFT.  \nIn the case of a two-sided black hole, it is well-established that space-like geodesics can traverse the Einstein-Rosen bridge, provided the boundary operators are inserted at the two  \ndisjoint AdS boundaries. In certain limits, these space-like geodesics become null, enabling them to approach arbitrarily close to the singularity. This property was exploited in [73, 74] to investigate how two-point CFT correlators encode information about the singularity, 17 and was further revisited in [11] in the context of holographic RG flows. Building on these pioneering works, we will now focus on calculating traversing geodesics to investigate potential signatures of the Kasner eons encoded in the boundary two-point correlators.  \nTo simplify the calculations, we assume that the boundary operators are symmetrically placed, ensuring that the space-like geodesics anchored at these endpoints exhibit reflection symmetry. The radial geodesics originate from one of the AdS boundaries, travel towards the singularity until they reach a turning point, and then return to the second AdS boundary. Further, they are confined to the ( t, r ) plane, where the induced metric is given by,  \nds 2 = 1 r 2 ( -f ( r ) e -χ ( r ) dt 2 + dr 2 f ( r ) ) . (3.25)  \nWe can parametrize the geodesics using a single function, r = r ( t ), allowing us to express the length functional as:  \nL = ∫ dt 1 r √ -f ( r ) e -χ ( r ) + ˙ r 2 f ( r ) ≡ ∫ dt L g , (3.26)  \nwhere L g represents the 'Lagrangian' of the geodesic, defined implicitly by (3.26).  \nSince L g does not explicitly depend on time, we can define a conserved quantity E , representing the energy associated with the space-like geodesic,  \nE = ˙ r ∂ L g ∂ ˙ r -L g = f ( r ) e -χ ( r ) r √ -f ( r ) e -χ ( r ) + ˙ r 2 f ( r ) . (3.27)  \nThe geodesic length is then obtained by substituting E into L g , yielding the minimal geodesic length anchored at some boundary time slice t = t b ( E ),  \nL = 2 | E | ∫ r t r c e -χ/ 2 dr r 2 √ 1 + f ( r ) e -χ ( r ) r 2 E 2 +2log r c . (3.28)  \nNote that a counterterm has been introduced to subtract the UV divergences arising from the AdS boundary. Here, r c represents the UV cutoff, while r t is the turning point of the geodesic, determined by  \nE 2 = g tt ( r t ) = -f ( r t ) e -χ ( r t ) r 2 t . (3.29)  \nUsing (3.27), we can express boundary time t b as a function of the energy E ,  \n<!-- image -->  \nFigure 4 : Energy of the two-sided geodesics that probe the black hole interior as a function of turning point r t . We consider d = 4 (left) and d = 5 (right) with λ 2 = 10 -8 , λ 3 = -10 -22 , λ 4 = -10 -42 , λ 5 = -10 -67 , λ n> 5 = 0 and m 2 = -d 2 4 + 1 2 . In both cases, the last eon has p t > 0 and the energy has therefore a maximum value. This implies that space-like geodesics get stuck at some finite r and cannot reach arbitrarily close to the singularity.  \n<!-- image -->  \nt b = -P ∫ r t 0 sgn( E ) e χ/ 2 dr f ( r ) √ 1 + f ( r ) e -χ ( r ) r 2 E 2 , (3.30)  \nwhere P denotes the principal value. The relation in (3.29) specifies the depth of these geodesics within the black hole. Real space-like geodesics can be anchored up to a maximum boundary time t c given by  \nt c = -P ∫ r m 0 sgn( E ) e χ/ 2 dr f ( r ) √ 1 + f ( r ) e -χ ( r ) r 2 E 2 , (3.31)  \nwhere r m is the maximum value of the turning point.  \nIn Figures 4 and 5, we plot the energy of the geodesics as a function of the turning point for different scenarios of interest. For a two-sided geodesic to reach the singularity, its energy must diverge, i.e., E →∞ . However, when the singularity resides within a Kasner eon characterized by p t > 0, there is always a finite maximum energy, hindering direct probing of the singularity. This phenomenon was first observed in [12]. 18 Notably, p t > 0 is quite generic for Kasner eons with matter (for n ≥ 2), as in these cases p t = n/d , cf. (2.31). This behavior is depicted in Figure 4. Conversely, for Kasner eons without matter, we have p t = 2 n/d -1 as shown in (2.10). Therefore, for n ≤ d/ 2, we have p t ≤ 0, and space-like geodesics are expected to reach the singularity as E →∞ . 19 This behavior is illustrated in Figure 5.  \nFigure 5 : Left: energy of the two-sided correlator for Gauss-Bonnet theory without matter in d = 5 with λ 2 = 10 -8 and λ n> 2 = 0. Right: the same energy evaluated for quasi-topological gravity theory up to quintic order without matter in d = 12 with λ 2 = 10 -8 , λ 3 = 10 -22 , λ 4 = 10 -42 , λ 5 = 10 -67 and λ n> 5 = 0. In both cases, the last eon has p t < 0, implying that space-like geodesics can reach arbitrarily close to the singularity.  \n<!-- image -->  \nWhether or not the final eon has p t > 0, space-like geodesics can still probe all preceding eons with p t ≤ 0. To effectively diagnose these eons, we find it convenient to perform a time parametrization, t b → τ b ( t b ) = log( r t ( t b )), and define the following combination:  \nL eff ( τ b ) ≡ L '' ( τ b ) L ' ( τ b ) -t '' b ( τ b ) t ' b ( τ b ) = p eff t ( τ b ) p eff x ( τ b ) , (3.32)  \nwhich perfectly captures the p t ≤ 0 Kasner eons, as illustrated in Figure 6. This provides a new means to identify stringy effects within the black hole interior. Crucially, L eff can be expressed entirely in terms of boundary quantities and their derivatives, with L effectively encoding the two-point correlator as described in (3.24). However, a subtle but noteworhy caveat arises: (3.32) presupposes a preferred notion of time, which may be difficult to identify from a purely boundary analysis (in part because the two-point correlators do not cease to exist at t b = t c , while the space-like geodesics do). An intriguing direction for future research would be to explore whether complex geodesics can reveal information about the correlator for t b > t c and, in particular, if they can be used to diagnose the remaining eons with p t > 0. We leave this investigation for future work.", '3.3 Generalized complexity': 'A further probe of the black hole interior that is worth investigating is quantum complexity [76, 77]. In [78, 79], the authors uncovered a broad class of gravitational observables exhibiting features typical of complexity, including late-time linear growth and the switchback effect. This framework, known as complexity=anything or generalized complexity, introduces novel  \n<!-- image -->  \nFigure 6 : L eff and t b as function of τ b in d = 4 (left) and d = 5 (right) with λ 2 = 10 -8 , λ 3 = -10 -22 , λ 4 = -10 -42 , λ 5 = -10 -67 , λ n> 5 = 0 and m 2 = -d 2 4 + 1 2 . The blue dashed line corresponds to the critical time, t c . When this time is reached, the space-like geodesic reaches its maximum depth, r m , situated in the last eon where p t ≤ 0.  \n<!-- image -->  \ngeometric constructs that probe the bulk and have holographic connections to various notions of quantum complexity in the boundary theory. 20  \nIn this section, we aim to investigate whether holographic complexity can be utilized to identify the presence of Kasner eons. Notably, the original complexity-volume (CV) prescription [89] cannot be used to probe the singularity, as extremal volume surfaces are typically trapped at a finite r . However, CV is merely one of many holographic proposals. The relationship between holographic complexities and the behavior near spacetime singularities has been explored in various contexts [15, 22, 28, 37, 39, 90-93]. Moreover, it is known that certain subclasses of generalized complexities can directly probe the singularity [37, 93]. Thus, it is crucial to examine how the Kasner exponents relate to key features of complexity, such as linear growth. In this section, we intend to employ a class of generalized complexities that probe the singularity to clarify their relationship with the Kasner exponents emerging in our particular family of higher derivative theories.  \nLet us start with a brief overview of [37, 93] and discuss how their approach can be adapted to diagnose the eons. The variants of complexity proposed in these studies define the complexity C in terms of a dimensionless scalar function, F , evaluated on a co-dimension one hypersurface, Σ, characterized by constant mean curvature (CMC),  \nC = 1 G N L ∫ Σ d d y √ h F ( y ) , (3.33)  \nwhere L is a length scale. The limit of particular interest occurs when the trace of the extrinsic curvature of the CMC slice diverges, K → ∞ . There are two such surfaces, corresponding  \n<!-- image -->  \nFigure 7 : Left: CMC slices Σ ± in the limit of large extrinsic curvature, at a finite boundary time. The shaded region is the WdW patch. Right: the same situation but at late boundary times. In this case Σ + approaches the τ = 0 slice while Σ -approaches the horizon, r = r h .  \n<!-- image -->  \nto the future and past slices of the Wheeler-deWitt (WdW) patch. See Figure 7 for an illustration. At late boundary times, one of these surfaces approaches the singularity, while the other approaches the horizon. We will focus our attention on the former.  \nWe will assume the following Kasner metric for the eons  \nd s 2 = -c 2 τ d τ 2 + c 2 t τ 2 p t d t 2 + τ 2 p x d ⃗x 2 d -1 . (3.34)  \nThis is equivalent to (2.1). However, we do not rescale t and r such that c r = c t = 1 since the growth of complexity can depend on the choice of these coordinates, and these are naturally inherited from the RG flows. A quick calculation shows that the late-time rate of growth of complexity is given by:  \nd C dt ∣ ∣ ∣ ∣ t →∞ = c t V d -1 G N L F ( τ ) τ p t +( d -1) p x , (3.35)  \nwhere V d -1 is the volume to the ( d -1)-dimensional spaced spanned by ⃗x . Reference [93] considered F = L | K | . Since | K | = p t +( d -1) p x c τ τ , using the Kasner relations of Einstein gravity one finds that the rate of growth is finite but independent of the Kasner exponents. A natural generalization for higher-curvature theories is F = L | K | m . In this case, we find that as τ → 0,  \nd C dt ∣ ∣ ∣ ∣ t →∞ =        0 if 0 ≤ m< 2 n -1 , c t V d -1 c 2 n -1 τ G N (2 n -1) 2 n -1 if m = 2 n -1 , ∞ if m> 2 n -1 , (3.36)  \nfor vacuum eons, while  \nd C dt ∣ ∣ ∣ ∣ t →∞ =        0 if 0 ≤ m<n, c t V d -1 c n τ G N n n if m = n, ∞ if m>n, (3.37)  \nfor matter eons with a scalar field. To arrive at these expressions we have used the linear Kasner relations (2.11) and (2.32), respectively.  \nReference [37] considered an alternative formulation of complexity, where the functional F is constructed from another invariant derived from the extrinsic curvature tensor: F = L √ | K µν K µν | . By recognizing that | K µν K µν | = p 2 t +( d -1) p 2 x c 2 τ τ 2 , and applying the Kasner relations of Einstein gravity, one finds that the growth rate is finite and directly dependent on the Kasner exponents. This version of complexity, therefore, effectively diagnoses the nature of the singularity in Einstein gravity. A natural extension to higher-curvature theories is given by F = L | K µν K µν | m/ 2 , where in the limit τ → 0, we find,  \nd C dt ∣ ∣ ∣ ∣ t →∞ =          0 if 0 ≤ m< 2 n -1 , c t V d -1 c 2 n -1 τ G N ( 4 n ( n -1) d +1 ) 2 n -1 2 if m = 2 n -1 , ∞ if m> 2 n -1 , (3.38)  \nfor vacuum eons, while  \nd C dt ∣ ∣ ∣ ∣ t →∞ =          0 if 0 ≤ m<n, c t V d -1 c n τ G N ( n 2 d ) n 2 if m = n, ∞ if m>n, (3.39)  \nfor matter eons with a scalar field. To arrive at these expressions we have used both the linear and quadratic Kasner relations given in (2.11) and (2.32), respectively.  \nAcouple of comments are in order. First, in the above expressions, n is meant to represent the final eon in which the singularity resides. For example, if the higher-curvature theory is truncated at the Gauss-Bonnet coupling, i.e. λ n> 2 = 0, then n = 2. 21 We will typically consider a situation where the last eon exhibits linear growth in complexity, as this is the expected behavior in the dual CFT. 22 A crucial difference with Einstein gravity, however, is that for higher-curvature theories the Kasner exponents are completely determined in terms of n . Given that the two cases discussed above result in growth rates that explicitly depend on n , we can conclude that both function choices yield equivalent insights into the Kasner scaling within the final eon, making them equally effective probes of the singularity. In principle, we should be able to extract information about the remaining eons and their transitions by examining the full time evolution of complexity and analyzing its t -scaling across different time scales. Since this is computationally more involved, we defer this analysis to future work.', '4 Discussion': "In general relativity, the BKL analysis offers a crucial framework for understanding spacetime dynamics near space-like singularities. However, a general framework for higher-curvature gravity theories is currently lacking. Our work takes a step in this direction by examining a specific class of higher-derivative theories of arbitrary order minimally coupled to a scalar field. In particular, we showed that the presence of matter significantly alters the near-singularity behavior, clearly distinguishing it from the vacuum case.  \nIn the first part of our paper, we considered a family of quasi-topological gravities in vacuum. By incorporating an infinite number of higher-curvature couplings, we explicitly demonstrated the emergence of an infinite sequence of Kasner eons as a space-like singularity is approached. Each eon is governed by emergent physics at energy scales dictated by higher-order curvature terms, with transitions occurring to successively higher-order eons as the singularity nears. This phenomenon arises independently of the cosmological constant, underscoring it as a fundamental characteristic of gravity. We further showed that these higher-order eons result in modified Kasner relations that differ from those established in Einstein gravity. Nevertheless, each higher-order eon can still be represented by a Kasner metric, with the exponents adhering to new constraints based on both the order of the eon and the dimension of spacetime.  \nIntroducing a minimally coupled free massive scalar field into our model significantly altered the near-singularity analysis. We first observed that the dynamics near the singularity remained unaffected by both the mass of the scalar field and the cosmological constant. As in the vacuum case, an infinite sequence of Kasner eons emerged. However, the presence of the scalar field led to a dramatic change in the Kasner exponents. While the vacuum scenario exhibited two distinct Kasner exponents, p t and p x , the inclusion of the scalar field caused these exponents to become equal for all higher-order eons. This indicates that the scalar field enhances the isotropy of the Kasner solution, rendering the temporal direction indistinguishable from the spatial directions. Importantly, we found that the full set of Kasner exponents, including p ϕ , are fully determined by the order of the eon and the spacetime dimensionality, in stark contrast with GR, where the exponents are not fully fixed. We attributed this difference to the absence of higher-derivative corrections in the matter sector.  \nIn the final part of our paper, we focused on cases with a negative cosmological constant and explored the emergence of eons in the interior of black holes with holographic duals. We examined three key holographic observables -the thermal a -function, two-point functions of heavy operators, and holographic complexity- as potential probes for understanding the behavior of spacetime near singularities. We first introduce non-perturbative generalizations of the thermal a -function applicable to hairy black holes in quasi-topological gravities of any curvature order. This generalization was achieved by employing two distinct conditions: the Null Energy Condition (NEC) and the Null Curvature Condition (NCC). While these conditions are equivalent in Einstein gravity, they generally differ in higher-curvature theories. Based on these conditions, we constructed two monotonic functions characterizing holographic  \nRG flows in the presence of a relevant deformation that effectively capture all Kasner eons within the black hole. This demonstrates that these monotonic functions serve as robust diagnostic tools for identifying stringy effects near black hole singularities.  \nNext, we analyzed two-point correlation functions of heavy operators, holographically dual to space-like geodesics anchored at the AdS boundary. Our study revealed that the presence of a scalar field, which induces a positive p t , generally prevents the geodesics from probing regions near the singularity. In contrast, in the absence of a scalar field, the geodesics can reach the singularity if the theory satisfies n ≤ d/ 2, resulting in p t ≤ 0. Regardless of the singularity's characteristics, we demonstrated that space-like geodesics can capture all initial eons with p t < 0, providing, at least partially, a valuable complementary probe of stringy effects in the black hole interior. In the future, it would be interesting to investigate whether complex geodesics can reveal additional aspects of the singularity, particularly their potential to diagnose eons that remain inaccessible through standard geodesics.  \nLastly, we investigated the growth rate of complexity within specific variants of the complexity=anything proposal. Evaluating complexity on a slice of constant mean curvature close to the singularity, we demonstrated that the last eon is effectively captured and can be diagnosed through late-time linear growth of complexity. While we did not delve into this in detail, we anticipate that other eons can be similarly identified by examining different t -scalings in the full-time evolution of complexity, which we plan to explore in future work.  \nOur work represents a crucial first step towards developing a comprehensive BKL-like analysis in higher-curvature gravity, paving the way for future theoretical advancements. Several interesting directions merit exploration. First, it would be interesting to analyze the approach to the singularity in four-dimensional black holes supplemented with highercurvature corrections. Indeed, our study with quasi-topological gravities applies for spacetime dimensions D ≥ 5, leaving the four-dimensional case uncovered. To explore it, one may resort to the so-called generalized quasi-topological gravities [95-97], which exist in four dimensions but whose equations for static and plane-symmetric equations are nevertheless much more involved. Investigation in this area will be conducted in the future. Second, we have focused only on higher-derivative theories coupled to a minimally coupled massive scalar field. We found that even in this relatively simple scenario, the presence of matter alters the GR predictions for the Kasner exponents, without inducing chaotic behavior. A significant future direction would be to explore more complex matter content, potentially giving rise to 'mixmaster' behavior. This would allow for a deeper understanding of the interplay between Kasner epochs, eras, and eons. Third, it would be interesting to investigate the effects of considering higher-curvature corrections in the matter sector, as well as non-minimal couplings. As explained above, we expect such generalizations to give more intricate dynamics and more freedom in the definition of Kasner exponents and Kasner relations associated to the different eons. Finally, while our focus was on space-like singularities within black holes, examining cosmological applications may be of interest too. For instance, we can consider holographic cosmologies from dynamical end-of-the-world (ETW) branes falling into a black hole [98-100]. These branes, whose worldvolume is a big-bang/big-crunch cosmology, end in  \nthe singularity, and thus are sensitive to the presence of eons. Thus, it is natural to ask what are the implications of our work for holographically cosmologies constructed in this way. We hope to return to some of these questions in the near future.", 'Acknowledgments': "We are grateful to Daniel Are'an, Pablo Bueno, Pablo A. Cano, Javier Carballo, Robie Hennigar, Hyun-Sik Jeong, Hanzhi Jiang, Gerben Oling, Le-Chen Qu and Simon Ross for useful discussions and correspondence. EC and ' AM thank the Instituto de F'ısica Te'orica UAMCSIC (IFT UAM-CSIC) for hospitality during the early stages of this project. The work of EC is supported by the National Science Foundation under grants No. PHY-2112725 and No. PHY-2210562. ' AM acknowledges support from the Istituto Nazionale di Fisica Nucleare through the Bando 23590 and the INFN Special Initiative String Theory and Fundamental Interactions . AKP and JFP are supported by the 'Atracci'on de Talento' program (Comunidad de Madrid) grant 2020-T1/TIC-20495, by the Spanish Research Agency via grants CEX2020001007-S and PID2021-123017NB-I00, funded by MCIN/AEI/10.13039/501100011033, and by ERDF A way of making Europe.", 'A Quasi-topological gravities with first-order vacuum equations for any static and plane-symmetric background': 'We consider a static and plane-symmetric ansatz of the form (2.5). For spacetime dimensions d +1 = D ≥ 5, it turns out that there always exists a higher-curvature gravity Z ( n ) which is purely constructed out of curvature invariants of n -th curvature order and whose equations of motion on top of (2.5) contain no more than two derivatives - as a matter of fact, equations are of first-order in derivatives for f ( r ) and χ ( r ). Such theories correspond to special subclasses of quasi-topological gravities and below we present explicit instances of such theories up to order n = 5 and arbitrary dimension. We denote by W abcd the Weyl curvature tensor and by Z ab = R ab -1 D g ab R the traceless part of the Ricci tensor.  \nZ (1) = R,  \nZ (2) = R 2 -4 D ( D -1) ( D -2) 2 Z ab Z ab + D ( D -1) ( D -2)( D -3) W abcd W abcd ,  \nZ (3) = R 3 + 2 D 2 ( D -1) 2 (2 D -3) W ab cd W cd ef W ef ab ( D -2)( D -3)( D (( D -9) D +26) -22) + 24 D 2 ( D -1) 2 W bd ac Z a b Z c d ( D -3)( D -2) 3 -24 D 2 ( D -1) 2 W acde W bcde Z a b ( D -3)( D -4)( D -2) 2 + 16 D 2 ( D -1) 2 Z a b Z b c Z c a ( D -2) 4 + 3( D -1) DRW abcd W abcd ( D -2)( D -3) -12 D ( D -1) Z b a Z a b R ( D -2) 2 ,  \nZ (4) = R 4 + 3 D 3 ( D -1) 2 (3 D -4) ( W abcd W abcd ) ( D -2) 4 ( D -3) 2 -384 D 3 ( D -1) 3 Z ac Z de W bdce Z a b ( D -3)( D -4)( D -2) 4  \n(A.1) (A.2) (A.3) 2  \n- (A.4)  \n+ 8 D 2 ( D -1) 2 (2 D -3) RW ab cd W cd ef W ef ab ( D -2)( D -3)( D (( D -9) D +26) -22) + 48 D 3 ( D -1) 3 ( Z b a Z a b ) 2 ( D -2) 5 ( D -3) -96 D 3 ( D -1) 3 Z b a Z c b Z d c Z a d ( D -2) 5 ( D -3) + 64 D 2 ( D -1) 2 RZ b a Z c b Z a c ( D -2) 4 + 96 D 2 ( D -1) 2 RW bd ac Z a b Z c d ( D -2) 3 ( D -3) -96 D 2 ( D -1) 2 RW acde W bcde Z a b ( D -3)( D -4)( D -2) 2 -48 D 2 ( D -1) 4 W abcd W abcd Z ef Z ef ( D -3)( D -2) 4 -192 D 2 ( D -1) 3 W acbd W cefg W d efg Z ab ( D 2 -5 D +6) 2 ( D -4) + 6 D ( D -1) R 2 W abcd W abcd ( D -2)( D -3) -24 D ( D -1) R 2 Z ab Z ab ( D -2) 2 + 96 D 2 (2 D -1)( D -1) 3 W abcd W aecf Z bc Z ef ( D -2) 3 ( D -3) 2 ,  \nZ (5) = R 5 + 4(4 D -5) D 4 ( D -1) 3 W cd ab W ef cd W ab ef W ghij W ghij ( D -3) 2 ( D -2) 3 ( D (( D -9) D +26) -22)  \n- -640 D 4 ( D -1) 4 Z de Z de Z b a Z c b Z a c ( D -3)( D -4)( D -2) 6 + 15 D 3 ( D -1) 2 (3 D -4) R ( W abcd W abcd ) 2 ( D -3) 2 ( D -2) 4 + 48 D 4 ( D -1) 3 ( 16( D -1)( D -3) Z a b Z b c Z c d Z d f Z f a -5( D -2) 2 W ghij W ghij W acde W bcde Z a b ) ( D -4)( D -3) 2 ( D -2) 6 -1920 D 3 ( D -1) 3 RW bdce Z a b Z ac Z de ( D -3)( D -4)( D -2) 4 + 240 D 3 ( D -1) 3 R ( Z ab Z ab ) 2 ( D -3)( D -2) 5 -480 D 3 ( D -1) 3 RZ b a Z c b Z d c Z a d ( D -3)( D -2) 5 -320( D -1) 4 D 3 ( D ( D +2) -2) W ghij W ghij Z a b Z b c Z c a ( D -4)( D -3)( D -2) 6 ( D +1) + 40 D 4 ( D -1) 3 (9 D 4 -88 D 3 +307 D 2 -440 D +206) W cd ab W ef cd W ab ef Z gh Z gh ( D -2) 3 ( D -3)( D -4)( D 2 -6 D +11)( D (( D -9) D +26) -22) + 480 D 4 ( D -1) 4 Z ef Z ef W bd ac Z a b Z c d ( D -3) 2 ( D -2) 5 + 20 D 2 ( D -1) 2 (2 D -3) R 2 W cd ab W ef cd W ab ef ( D -2)( D -3)( D (( D -9) D +26) -22) + 240 D 2 ( D -1) 2 R 2 W bd ac Z a b Z c d ( D -3)( D -2) 3 + 160 D 2 ( D -1) 2 R 2 Z b a Z c b Z a c ( D -2) 4 -960 D 2 ( D -1) 3 RW acbd W cefg W d efg Z ab ( D -4) ( D 2 -5 D +6) 2 -240 D 2 ( D -1) 2 R 2 W acde W bcde Z a b ( D -3)( D -4)( D -2) 2 -240 D 2 ( D -1) 4 RW abcd W abcd Z ef Z ef ( D -3)( D -2) 4 + 10 D ( D -1) DR 3 W abcd W abcd ( D -2)( D -3) -40 D ( D -1) R 3 Z ab Z ab ( D -2) 2 + 480 D 2 (2 D -1)( D -1) 3 RW abcd W aecf Z bc Z ef ( D -2) 3 ( D -3) 2 -1920 D 3 (3 D -1)( D -1) 4 Z ab W acbd Z ef W c g e f Z dg ( D -2) 4 ( D -3) 2 ( D -4)( D +1) -960 D 4 ( D -1) 4 Z b a Z c b Z cd Z ef W eafd ( D -2) 5 ( D -3) 2 -(A.5)\n- 240 D 4 ( D -1) 3 (3 D -2) Z a b Z b c W daef W efgh W gh dc ( D ( D -6) + 11)( D -2) 3 ( D -3)( D -4) .  \nHaving at disposal quasi-topological gravities with first-order vacuum equations on any static and plane-symmetric background up to quintic order in curvature, it is possible to obtain such theories to arbitrary order in the curvature and any spacetime dimension D via the following recursive formula, adapted from [10, 44]:  \nZ ( n +5) = 3( n +3) R Z ( n +4) n +1 -3( n +4) Z (2) Z ( n +3) n + ( n +3)( n +4) Z (3) Z ( n +2) n ( n +1) , (A.6)  \nwhere n ≥ 1. Surprisingly enough, written in this way the recursive relation (A.6) is dimension-independent.', 'B Holographic dictionary for bulk scalar fields': "In this appendix, we discuss some basic entries of the holographic dictionary for the theory at hand, given by (3.1). By varying the action with respect to the metric g µν and scalar field ϕ , we obtain following equations of motions,  \nϕ '' + ( f ' f -d -1 r -χ ' 2 ) ϕ ' + ∆( d -∆) r 2 f ϕ = 0 , (B.1) χ ' -2 f ' f -1 ∑ ∞ n =1 ˜ λ n f n -1 ( ∆( d -∆) ϕ 2 ( d -1) rf + 2 d ( d -1) rf -2 d r ∞ ∑ n =1 ˜ λ n n f n -1 ) = 0 , (B.2) χ ' -r ( ϕ ' ) 2 ∑ ∞ n =1 ˜ λ n f n -1 = 0 , (B.3)  \nwhere m is expressed in terms of the scaling dimension ∆ of the operator O dual to ϕ [101, 102]:  \nm 2 = ∆(∆ -d ) . (B.4)  \nFrom this relation (B.4), we see that the dual operator is relevant only if ∆ < d , which occurs when m 2 < 0. Using the Breitenlohner-Freedman stability bound [52, 53], we can further determine a lower bound for m 2 , which yields  \n-d 2 4 ≤ m 2 < 0 . (B.5)  \nAdditionally, for holographic applications, we impose the unitarity bound for ∆ in order to have a unitary theory,  \n∆ ≥ d -2 2 . (B.6)  \nIn the standard quantization, this bound further restricts m 2 , so that 23  \n-d 2 4 ≤ m 2 ≤ 1 -d 2 / 4 . (B.7)  \nWe will thus consider massive scalar fields within this range.  \nTo solve the equations in (B.1)-(B.3), we need appropriate boundary conditions for bulk fields. As the fields approach the boundary at r → 0 with ∆ = d/ 2, they behave as [14]:  \n̸  \nϕ ( r ) ∼ ϕ 0 r d -∆ , (B.8)  \nχ ( r ) ∼ d -∆ 2( d -1) ϕ 2 0 r 2( d -∆) , (B.9)  \nf ( r ) ∼ 1 , (B.10)  \nwhere ϕ 0 is the source of the operator O dual to ϕ . Note that for ∆ = d/ 2, the above expressions fail, and need to be modified. For this special case, ∆ = d/ 2, it has been found that they instead behave as [14],  \nϕ ( r ) ∼ ϕ 0 r d/ 2 log r , (B.11)  \nχ ( r ) ∼ ϕ 2 0 4 d ( d -1) r d (2 + 2 d log r + d 2 (log r ) 2 ) , (B.12)  \nf ( r ) ∼ 1 , (B.13)  \nOn the other hand, near the singularity as r →∞ , they generally diverge as  \nf ( r ) = -f n r γ n , χ ( r ) = χ n log r + χ 0 , ϕ ( r ) = { ϕ 1 log r ( n = 1) , ϕ n r ( n -1) γn 2 ( n ≥ 2) , (B.14)  \nwhere χ n = γ n = 2 d n for n ≥ 2, f n > 0 and χ 0 are some arbitrary constants and  \nχ 1 = ϕ 2 1 ( d -1) , γ 1 = d + ϕ 2 1 2( d -1) . (B.15)", 'References': "- [1] V. A. Belinsky, I. M. Khalatnikov and E. M. Lifshitz, Oscillatory approach to a singular point in the relativistic cosmology , Adv. Phys. 19 (1970) 525-573.\n- [2] V. A. Belinsky and I. M. Khalatnikov, Effect of Scalar and Vector Fields on the Nature of the Cosmological Singularity , Sov. Phys. JETP 36 (1973) 591.\n- [3] V. A. Belinsky and I. M. Khalatnikov, On the Influence of Matter and Physical Fields Upon the Nature of Cosmological Singularities , Sov. Sci. Rev. A 3 (1981) 555-590.\n- [4] G. Oling and J. F. Pedraza, Mixmasters in Wonderland: Chaotic dynamics from Carroll limits of gravity , 2409.05836 .\n- [5] T. Damour, M. Henneaux and H. Nicolai, Cosmological billiards , Class. Quant. Grav. 20 (2003) R145-R200 [ hep-th/0212256 ].\n- [6] V. Belinski and M. Henneaux, The Cosmological Singularity . Cambridge Monogr.Math.Phys. Cambridge Univ. Pr., Cambridge, 2017.  \n- [7] B. K. Berger, Numerical approaches to space-time singularities , Living Rev. Rel. 5 (2002) 1 [ gr-qc/0201056 ].\n- [8] J. M. Heinzle, C. Uggla and W. C. Lim, Spike Oscillations , Phys. Rev. D 86 (2012) 104049 [ 1206.0932 ].\n- [9] E. Czuchry, D. Garfinkle, J. R. Klauder and W. Piechocki, Do spikes persist in a quantum treatment of spacetime singularities? , Phys. Rev. D 95 (2017), no. 2 024014 [ 1605.04648 ].\n- [10] P. Bueno, P. A. Cano and R. A. Hennigar, Kasner Epochs, Eras and Eons , 2402.14912 .\n- [11] A. Frenkel, S. A. Hartnoll, J. Kruthoff and Z. D. Shi, Holographic flows from CFT to the Kasner universe , JHEP 08 (2020) 003 [ 2004.01192 ].\n- [12] S. A. Hartnoll, G. T. Horowitz, J. Kruthoff and J. E. Santos, Gravitational duals to the grand canonical ensemble abhor Cauchy horizons , JHEP 10 (2020) 102 [ 2006.10056 ].\n- [13] S. A. Hartnoll, G. T. Horowitz, J. Kruthoff and J. E. Santos, Diving into a holographic superconductor , SciPost Phys. 10 (2021), no. 1 009 [ 2008.12786 ].\n- [14] E. Caceres, A. Kundu, A. K. Patra and S. Shashi, Page curves and bath deformations , SciPost Phys. Core 5 (2022) 033 [ 2107.00022 ].\n- [15] A. Bhattacharya, A. Bhattacharyya, P. Nandy and A. K. Patra, Bath deformations, islands, and holographic complexity , Phys. Rev. D 105 (2022), no. 6 066019 [ 2112.06967 ].\n- [16] L. Sword and D. Vegh, Kasner geometries inside holographic superconductors , JHEP 04 (2022) 135 [ 2112.14177 ].\n- [17] L. Sword and D. Vegh, What lies beyond the horizon of a holographic p-wave superconductor , JHEP 12 (2022) 045 [ 2210.01046 ].\n- [18] Y.-Q. Wang, Y. Song, Q. Xiang, S.-W. Wei, T. Zhu and Y.-X. Liu, Holographic flows with scalar self-interaction toward the Kasner universe , 2009.06277 .\n- [19] S. A. H. Mansoori, L. Li, M. Rafiee and M. Baggioli, What's inside a hairy black hole in massive gravity? , JHEP 10 (2021) 098 [ 2108.01471 ].\n- [20] Y. Liu, H.-D. Lyu and A. Raju, Black hole singularities across phase transitions , JHEP 10 (2021) 140 [ 2108.04554 ].\n- [21] S. Das and A. Kundu, RG flows and thermofield-double states in holography , JHEP 04 (2022) 167 [ 2112.11675 ].\n- [22] E. Caceres, A. Kundu, A. K. Patra and S. Shashi, Trans-IR flows to black hole singularities , Phys. Rev. D 106 (2022), no. 4 046005 [ 2201.06579 ].\n- [23] Y.-S. An, L. Li, F.-G. Yang and R.-Q. Yang, Interior structure and complexity growth rate of holographic superconductor from M-theory , JHEP 08 (2022) 133 [ 2205.02442 ].\n- [24] R. Auzzi, S. Bolognesi, E. Rabinovici, F. I. Schaposnik Massolo and G. Tallarita, On the time dependence of holographic complexity for charged AdS black holes with scalar hair , JHEP 08 (2022) 235 [ 2205.03365 ].\n- [25] M. Mirjalali, S. A. Hosseini Mansoori, L. Shahkarami and M. Rafiee, Probing inside a charged hairy black hole in massive gravity , JHEP 09 (2022) 222 [ 2206.02128 ].  \n- [26] S. A. Hartnoll, Wheeler-DeWitt states of the AdS-Schwarzschild interior , JHEP 01 (2023) 066 [ 2208.04348 ].\n- [27] S. A. Hartnoll and N. Neogi, AdS black holes with a bouncing interior , SciPost Phys. 14 (2023), no. 4 074 [ 2209.12999 ].\n- [28] E. Caceres and S. Shashi, Anisotropic flows into black holes , JHEP 01 (2023) 007 [ 2209.06818 ].\n- [29] Y. Liu and H.-D. Lyu, Interior of helical black holes , JHEP 09 (2022) 071 [ 2205.14803 ].\n- [30] L.-L. Gao, Y. Liu and H.-D. Lyu, Black hole interiors in holographic topological semimetals , JHEP 03 (2023) 034 [ 2301.01468 ].\n- [31] E. Caceres, S. Shashi and H.-Y. Sun, Imprints of phase transitions on Kasner singularities , Phys. Rev. D 109 (2024), no. 12 126018 [ 2305.11177 ].\n- [32] M. J. Blacker and S. Ning, Wheeler DeWitt states of a charged AdS 4 black hole , JHEP 12 (2023) 002 [ 2308.00040 ].\n- [33] E. C'aceres, A. K. Patra and J. F. Pedraza, Shock waves, black hole interiors and holographic RG flows , JHEP 07 (2024) 052 [ 2311.12940 ].\n- [34] M. De Clerck, S. A. Hartnoll and J. E. Santos, Mixmaster chaos in an AdS black hole interior , JHEP 07 (2024) 202 [ 2312.11622 ].\n- [35] R.-G. Cai, M.-N. Duan, L. Li and F.-G. Yang, Towards classifying the interior dynamics of charged black holes with scalar hair , JHEP 02 (2024) 169 [ 2312.11131 ].\n- [36] L.-L. Gao, Y. Liu and H.-D. Lyu, Internal structure of hairy rotating black holes in three dimensions , JHEP 01 (2024) 063 [ 2310.15781 ].\n- [37] D. Are'an, H.-S. Jeong, J. F. Pedraza and L.-C. Qu, Kasner interiors from analytic hairy black holes , 2407.18430 .\n- [38] R.-G. Cai, M.-N. Duan, L. Li and F.-G. Yang, Clarifying Kasner dynamics inside anisotropic black hole with vector hair , 2408.06122 .\n- [39] J. Carballo, A. K. Patra and J. F. Pedraza, Diving inside holographic metals , 2408.07748 .\n- [40] J. Oliva and S. Ray, A new cubic theory of gravity in five dimensions: Black hole, Birkhoff's theorem and C-function , Class. Quant. Grav. 27 (2010) 225002 [ 1003.4773 ].\n- [41] R. C. Myers and B. Robinson, Black Holes in Quasi-topological Gravity , JHEP 08 (2010) 067 [ 1003.5357 ].\n- [42] M. H. Dehghani, A. Bazrafshan, R. B. Mann, M. R. Mehdizadeh, M. Ghanaatian and M. H. Vahidinia, Black Holes in Quartic Quasitopological Gravity , Phys. Rev. D85 (2012) 104009 [ 1109.4708 ].\n- [43] P. Bueno, P. A. Cano, J. Moreno and ' A. Murcia, All higher-curvature gravities as Generalized quasi-topological gravities , JHEP 11 (2019) 062 [ 1906.00987 ].\n- [44] P. Bueno, P. A. Cano and R. A. Hennigar, (Generalized) quasi-topological gravities at all orders , Class. Quant. Grav. 37 (2020), no. 1 015002 [ 1909.07983 ].\n- [45] J. Moreno and ' A. J. Murcia, Classification of generalized quasitopological gravities , Phys. Rev. D 108 (2023), no. 4 044016 [ 2304.08510 ].  \n- [46] J. Moreno and ' A. J. Murcia, Cosmological higher-curvature gravities , Class. Quant. Grav. 41 (2024), no. 13 135017 [ 2311.12104 ].\n- [47] T. Padmanabhan and D. Kothawala, Lanczos-Lovelock models of gravity , Phys. Rept. 531 (2013) 115-171 [ 1302.2151 ].\n- [48] P. Bueno, P. A. Cano, R. A. Hennigar, M. Lu and J. Moreno, Generalized quasi-topological gravities: the whole shebang , Class. Quant. Grav. 40 (2023), no. 1 015004 [ 2203.05589 ].\n- [49] P. A. Cano and ' A. Murcia, Resolution of Reissner-Nordstrom singularities by higher-derivative corrections , Class. Quant. Grav. 38 (2021), no. 7 075014 [ 2006.15149 ].\n- [50] P. Bueno, P. A. Cano and R. A. Hennigar, Regular Black Holes From Pure Gravity , 2403.04827 .\n- [51] X. O. Camanho, N. Dadhich and A. Molina, Pure Lovelock Kasner metrics , Class. Quant. Grav. 32 (2015), no. 17 175016 [ 1505.06675 ].\n- [52] P. Breitenlohner and D. Z. Freedman, Positive Energy in anti-De Sitter Backgrounds and Gauged Extended Supergravity , Phys. Lett. B115 (1982) 197.\n- [53] P. Breitenlohner and D. Z. Freedman, Stability in Gauged Extended Supergravity , Annals Phys. 144 (1982) 249.\n- [54] A. B. Zamolodchikov, Irreversibility of the Flux of the Renormalization Group in a 2D Field Theory , JETP Lett. 43 (1986) 730-732. [Pisma Zh. Eksp. Teor. Fiz.43,565(1986)].\n- [55] Z. Komargodski and A. Schwimmer, On Renormalization Group Flows in Four Dimensions , JHEP 12 (2011) 099 [ 1107.3987 ].\n- [56] Z. Komargodski, The Constraints of Conformal Symmetry on RG Flows , JHEP 07 (2012) 069 [ 1112.4538 ].\n- [57] H. Casini and M. Huerta, A Finite entanglement entropy and the c-theorem , Phys. Lett. B 600 (2004) 142-150 [ hep-th/0405111 ].\n- [58] H. Casini and M. Huerta, On the RG running of the entanglement entropy of a circle , Phys. Rev. D85 (2012) 125016 [ 1202.5650 ].\n- [59] H. Casini, E. Test'e and G. Torroba, Markov Property of the Conformal Field Theory Vacuum and the a Theorem , Phys. Rev. Lett. 118 (2017), no. 26 261602 [ 1704.01870 ].\n- [60] L. Susskind and E. Witten, The Holographic bound in anti-de Sitter space , hep-th/9805114 .\n- [61] A. W. Peet and J. Polchinski, UV / IR relations in AdS dynamics , Phys. Rev. D 59 (1999) 065011 [ hep-th/9809022 ].\n- [62] Y. Hatta, E. Iancu, A. H. Mueller and D. N. Triantafyllopoulos, Aspects of the UV/IR correspondence : energy broadening and string fluctuations , JHEP 02 (2011) 065 [ 1011.3763 ].\n- [63] C. A. Ag'on, A. Guijosa and J. F. Pedraza, Radiation and a dynamical UV/IR connection in AdS/CFT , JHEP 06 (2014) 043 [ 1402.5961 ].\n- [64] D. Z. Freedman, S. S. Gubser, K. Pilch and N. P. Warner, Renormalization group flows from holography supersymmetry and a c theorem , Adv. Theor. Math. Phys. 3 (1999) 363-417 [ hep-th/9904017 ].  \n- [65] R. C. Myers and A. Sinha, Seeing a c-theorem with holography , Phys. Rev. D82 (2010) 046006 [ 1006.1263 ].\n- [66] R. C. Myers and A. Sinha, Holographic c-theorems in arbitrary dimensions , JHEP 01 (2011) 125 [ 1011.5819 ].\n- [67] A. Gonz'alez Lezcano, J. Hong, J. T. Liu, L. A. Pando Zayas and C. F. Uhlemann, c-functions in flows across dimensions , JHEP 10 (2022) 083 [ 2207.09360 ].\n- [68] E. Deddo, J. T. Liu, L. A. Pando Zayas and R. J. Saskowski, c-functions in higher-derivative flows across dimensions , JHEP 08 (2023) 147 [ 2305.18530 ].\n- [69] E. C'aceres, R. Castillo V'asquez, K. Landsteiner and I. Salazar Landea, Holographic a-functions and Boomerang RG flows , JHEP 02 (2024) 019 [ 2310.15983 ].\n- [70] R. C. Myers, M. F. Paulos and A. Sinha, Holographic studies of quasi-topological gravity , JHEP 08 (2010) 035 [ 1004.2055 ].\n- [71] V. Balasubramanian and S. F. Ross, Holographic particle detection , Phys. Rev. D 61 (2000) 044007 [ hep-th/9906226 ].\n- [72] J. Louko, D. Marolf and S. F. Ross, On geodesic propagators and black hole holography , Phys. Rev. D 62 (2000) 044041 [ hep-th/0002111 ].\n- [73] L. Fidkowski, V. Hubeny, M. Kleban and S. Shenker, The Black hole singularity in AdS/CFT , JHEP 02 (2004) 014 [ hep-th/0306170 ].\n- [74] G. Festuccia and H. Liu, Excursions beyond the horizon: Black hole singularities in Yang-Mills theories. I. , JHEP 04 (2006) 044 [ hep-th/0506202 ].\n- [75] G. T. Horowitz, H. Leung, L. Queimada and Y. Zhao, Boundary signature of singularity in the presence of a shock wave , SciPost Phys. 16 (2024), no. 2 060 [ 2310.03076 ].\n- [76] L. Susskind, Computational Complexity and Black Hole Horizons , Fortsch. Phys. 64 (2016) 44-48 [ 1403.5695 ]. [Fortsch. Phys.64,24(2016)].\n- [77] L. Susskind, Entanglement is not enough , Fortsch. Phys. 64 (2016) 49-71 [ 1411.0690 ].\n- [78] A. Belin, R. C. Myers, S.-M. Ruan, G. S'arosi and A. J. Speranza, Does Complexity Equal Anything? , Phys. Rev. Lett. 128 (2022), no. 8 081602 [ 2111.02429 ].\n- [79] A. Belin, R. C. Myers, S.-M. Ruan, G. S'arosi and A. J. Speranza, Complexity equals anything II , JHEP 01 (2023) 154 [ 2210.09647 ].\n- [80] M. A. Nielsen, M. R. Dowling, M. Gu and A. C. Doherty, Quantum Computation as Geometry , Science 311 (2006), no. 5764 1133-1135 [ quant-ph/0603161 ].\n- [81] M. R. Dowling and M. A. Nielsen, The geometry of quantum computation , Quant. Inf. Comput. 8 (2008), no. 10 0861-0899 [ quant-ph/0701004 ].\n- [82] R. Jefferson and R. C. Myers, Circuit complexity in quantum field theory , JHEP 10 (2017) 107 [ 1707.08570 ].\n- [83] S. Chapman, M. P. Heller, H. Marrochio and F. Pastawski, Toward a Definition of Complexity for Quantum Field Theory States , Phys. Rev. Lett. 120 (2018), no. 12 121602 [ 1707.08582 ].\n- [84] E. Caceres, S. Chapman, J. D. Couch, J. P. Hernandez, R. C. Myers and S.-M. Ruan, Complexity of Mixed States in QFT and Holography , JHEP 03 (2020) 012 [ 1909.10557 ].  \n- [85] E. Caceres, R. Carrasco and V. Patil, Lorentzian threads and generalized complexity , JHEP 04 (2024) 010 [ 2312.10606 ].\n- [86] J. F. Pedraza, A. Russo, A. Svesko and Z. Weller-Davies, Lorentzian Threads as Gatelines and Holographic Complexity , Phys. Rev. Lett. 127 (2021), no. 27 271602 [ 2105.12735 ].\n- [87] J. F. Pedraza, A. Russo, A. Svesko and Z. Weller-Davies, Sewing spacetime with Lorentzian threads: complexity and the emergence of time in quantum gravity , JHEP 02 (2022) 093 [ 2106.12585 ].\n- [88] J. F. Pedraza, A. Russo, A. Svesko and Z. Weller-Davies, Computing spacetime , Int. J. Mod. Phys. D 31 (2022), no. 14 2242010 [ 2205.05705 ].\n- [89] D. Stanford and L. Susskind, Complexity and Shock Wave Geometries , Phys. Rev. D90 (2014), no. 12 126007 [ 1406.2678 ].\n- [90] S. Bolognesi, E. Rabinovici and S. R. Roy, On Some Universal Features of the Holographic Quantum Complexity of Bulk Singularities , JHEP 06 (2018) 016 [ 1802.02045 ].\n- [91] P. Caputa, D. Das and S. R. Das, Path integral complexity and Kasner singularities , JHEP 01 (2022) 150 [ 2111.04405 ].\n- [92] G. Katoch, J. Ren and S. R. Roy, Quantum complexity and bulk timelike singularities , JHEP 12 (2023) 085 [ 2303.02752 ].\n- [93] E. Jørstad, R. C. Myers and S.-M. Ruan, Complexity=anything: singularity probes , JHEP 07 (2023) 223 [ 2304.05453 ].\n- [94] L. V. Iliesiu, M. Mezei and G. S'arosi, The volume of the black hole interior at late times , JHEP 07 (2022) 073 [ 2107.06286 ].\n- [95] P. Bueno and P. A. Cano, Einsteinian cubic gravity , Phys. Rev. D94 (2016), no. 10 104005 [ 1607.06463 ].\n- [96] R. A. Hennigar, D. Kubizˇn'ak and R. B. Mann, Generalized quasitopological gravity , Phys. Rev. D95 (2017), no. 10 104042 [ 1703.01631 ].\n- [97] J. Ahmed, R. A. Hennigar, R. B. Mann and M. Mir, Quintessential Quartic Quasi-topological Quartet , JHEP 05 (2017) 134 [ 1703.11007 ].\n- [98] S. Cooper, M. Rozali, B. Swingle, M. Van Raamsdonk, C. Waddell and D. Wakeham, Black hole microstate cosmology , JHEP 07 (2019) 065 [ 1810.10601 ].\n- [99] S. Antonini and B. Swingle, Cosmology at the end of the world , Nature Phys. 16 (2020), no. 8 881-886 [ 1907.06667 ].\n- [100] S. Fallows and S. F. Ross, Constraints on cosmologies inside black holes , JHEP 05 (2022) 094 [ 2203.02523 ].\n- [101] E. Witten, Anti-de Sitter space and holography , Adv. Theor. Math. Phys. 2 (1998) 253-291 [ hep-th/9802150 ].\n- [102] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Gauge theory correlators from noncritical string theory , Phys. Lett. B428 (1998) 105-114 [ hep-th/9802109 ]."}