Subregion complexity in holographic thermalization with dS boundary

Shao-Jun Zhang

arXiv:1905.10605·hep-th·September 4, 2019

Subregion complexity in holographic thermalization with dS boundary

Shao-Jun Zhang

PDF

TL;DR

This paper investigates the evolution of holographic subregion complexity in a Vaidya spacetime with a de Sitter boundary, revealing stages of growth, saturation, and dependence on subregion size and horizon proximity.

Contribution

It provides a detailed analysis of holographic subregion complexity dynamics in dS boundary conditions, extending previous flat boundary studies and deriving explicit relations for saturation time.

Findings

01

HSC exhibits four distinct evolution stages similar to flat boundary cases.

02

The linear growth rate depends mainly on the mass parameter.

03

Saturation time scales with subregion size as tt_{sat} = anh^{-1} ( ilde{R})".

Abstract

We study the time evolution of holographic subregion complexity (HSC) in Vaidya spacetime with dS boundary. The subregion on the boundary is chosen to be a sphere within the cosmological horizon. It is found that the behaviour of HSC is similar to that in cases with flat boundary. The whole evolution can be divided into four stages: First, it grows almost linearly, then the growth slows down; After reaching a maximum it drops down quickly and gets to saturation finally. The linear growth rate in the first stage is found to depend almost only on the the mass parameter. As the subregion size approaches the cosmological horizon, this stage is expected to last forever with the subsequent three stages washed out. The saturation time $\tilde{t}_{s a t}$ depends almost only on the subregion size $\tilde{R}$ as $\tilde{t}_{s a t} = tanh^{- 1} (\tilde{R})$ which is linear in $\tilde{R}$ when…

Equations68

C_{V} (Σ) = \partial B = Σ max (\frac{Vol ( B )}{G _{N} R}),

C_{V} (Σ) = \partial B = Σ max (\frac{Vol ( B )}{G _{N} R}),

C (γ_{A}) = \frac{Vol ( Γ _{A} )}{8 π G _{d + 1} L},

C (γ_{A}) = \frac{Vol ( Γ _{A} )}{8 π G _{d + 1} L},

S = \frac{1}{16 π G _{N}} \int d^{d + 1} x - g (R - 2Λ),

S = \frac{1}{16 π G _{N}} \int d^{d + 1} x - g (R - 2Λ),

R_{μν} - \frac{1}{2} g_{μν} R + Λ g_{μν} = 0.

R_{μν} - \frac{1}{2} g_{μν} R + Λ g_{μν} = 0.

d s^{2} = \frac{L ^{2}}{z ^{2}} (g_{μν} (z, x) d x^{μ} d x^{ν} + d z^{2}),

d s^{2} = \frac{L ^{2}}{z ^{2}} (g_{μν} (z, x) d x^{μ} d x^{ν} + d z^{2}),

d s^{2} = \frac{L ^{2}}{z ^{2}} (- f (r) g (z) d t^{2} + f^{- 1} (r) g (z) d r^{2} + r^{2} g (z) d Ω_{d - 2}^{2} + d z^{2}),

d s^{2} = \frac{L ^{2}}{z ^{2}} (- f (r) g (z) d t^{2} + f^{- 1} (r) g (z) d r^{2} + r^{2} g (z) d Ω_{d - 2}^{2} + d z^{2}),

f (r) = 1 - H^{2} r^{2}, g (z) = (1 - \frac{H ^{2} z ^{2}}{4})^{2} .

d s_{Σ}^{2} = - f (r) d t^{2} + f^{- 1} (r) d r^{2} + r^{2} d Ω_{d - 2}^{2},

d s_{Σ}^{2} = - f (r) d t^{2} + f^{- 1} (r) d r^{2} + r^{2} d Ω_{d - 2}^{2},

d s^{2} = \frac{L ^{2}}{z ^{2}} (- h (z) d t^{2} + \frac{d z ^{2}}{h ( z )} + \frac{H ^{2} L ^{2}}{f ( r ) ^{2}} d r^{2} + \frac{H ^{2} L ^{2}}{f ( r )} r^{2} d Ω_{d - 2}^{2}),

d s^{2} = \frac{L ^{2}}{z ^{2}} (- h (z) d t^{2} + \frac{d z ^{2}}{h ( z )} + \frac{H ^{2} L ^{2}}{f ( r ) ^{2}} d r^{2} + \frac{H ^{2} L ^{2}}{f ( r )} r^{2} d Ω_{d - 2}^{2}),

h (z) = 1 - \frac{z ^{2}}{L ^{2}} - \frac{m z ^{d}}{L ^{2 (d - 1)}} .

m = \frac{L ^{2 (d - 1)}}{z _{+}^{d}} (1 - \frac{z _{+}^{2}}{L ^{2}}) .

m = \frac{L ^{2 (d - 1)}}{z _{+}^{d}} (1 - \frac{z _{+}^{2}}{L ^{2}}) .

T_{H} = \frac{L ^{2} d - ( d - 2 ) z _{+}^{2}}{4 π L ^{2} z _{+}} .

T_{H} = \frac{L ^{2} d - ( d - 2 ) z _{+}^{2}}{4 π L ^{2} z _{+}} .

z_{+}^{e x t} = L \frac{d}{d - 2}, m^{e x t} = - \frac{2 L ^{2 (d - 1)}}{( d - 2 ) ( z _{+}^{e x t} ) ^{d}}

z_{+}^{e x t} = L \frac{d}{d - 2}, m^{e x t} = - \frac{2 L ^{2 (d - 1)}}{( d - 2 ) ( z _{+}^{e x t} ) ^{d}}

d s^{2} = \frac{L ^{2}}{z ^{2}} (- h (v, z) d v^{2} - 2 d v d z + \frac{H ^{2} L ^{2}}{f ( r ) ^{2}} d r^{2} + \frac{H ^{2} L ^{2}}{f ( r )} r^{2} d Ω_{d - 2}^{2}),

d s^{2} = \frac{L ^{2}}{z ^{2}} (- h (v, z) d v^{2} - 2 d v d z + \frac{H ^{2} L ^{2}}{f ( r ) ^{2}} d r^{2} + \frac{H ^{2} L ^{2}}{f ( r )} r^{2} d Ω_{d - 2}^{2}),

h (v, z) = 1 - \frac{z ^{2}}{L ^{2}} - \frac{m ( v ) z ^{d}}{L ^{2 (d - 1)}} .

R_{μν} - \frac{1}{2} g_{μν} R + Λ g_{μν} = 8 π G_{N} T_{μν}^{e x t},

R_{μν} - \frac{1}{2} g_{μν} R + Λ g_{μν} = 8 π G_{N} T_{μν}^{e x t},

8 π G_{N} T_{μν}^{e x t} = \frac{( d - 1 ) z ^{d - 1}}{2 L ^{2 (d - 1)}} \frac{d m ( v )}{d v} δ_{μv} δ_{ν v},

m (v) = \frac{M}{2} [1 + tanh (\frac{v}{v _{0}})],

m (v) = \frac{M}{2} [1 + tanh (\frac{v}{v _{0}})],

z (\tilde{R}) = ϵ, v (\tilde{R}) = \tilde{t},

z (\tilde{R}) = ϵ, v (\tilde{R}) = \tilde{t},

z^{'} (0) = v^{'} (0) = 0, z (0) = z_{*}, v (0) = v_{*},

z^{'} (0) = v^{'} (0) = 0, z (0) = z_{*}, v (0) = v_{*},

d s^{2} = \frac{L ^{2}}{z ^{2}} (\frac{L ^{2}}{( 1 - r ~ ^{2} ) ^{2}} - h (v, z) v^{'2} - 2 v^{'} z^{'}) d \tilde{r}^{2} + \frac{L ^{4}}{z ^{2}} \frac{r ~ ^{2}}{1 - r ~ ^{2}} d Ω_{d - 2}^{2} .

d s^{2} = \frac{L ^{2}}{z ^{2}} (\frac{L ^{2}}{( 1 - r ~ ^{2} ) ^{2}} - h (v, z) v^{'2} - 2 v^{'} z^{'}) d \tilde{r}^{2} + \frac{L ^{4}}{z ^{2}} \frac{r ~ ^{2}}{1 - r ~ ^{2}} d Ω_{d - 2}^{2} .

S = \frac{L ^{2 d - 3}}{4 G _{N}} Ω_{d - 2} \int_{0}^{\tilde{R}} \frac{d r ~}{z ^{d - 1}} Q P^{d - 2},

S = \frac{L ^{2 d - 3}}{4 G _{N}} Ω_{d - 2} \int_{0}^{\tilde{R}} \frac{d r ~}{z ^{d - 1}} Q P^{d - 2},

Q \equiv \frac{L ^{2}}{( 1 - r ~ ^{2} ) ^{2}} - h (v, z) v^{'2} - 2 v^{'} z^{'}, P \equiv \frac{r ~}{1 - r ~ ^{2}} .

h^{2} (\tilde{r}^{2} - 1)^{2} v^{'3} ((d - 1) (\tilde{r}^{2} - 1) \tilde{r} z^{'} + (d - 2) z)

h^{2} (\tilde{r}^{2} - 1)^{2} v^{'3} ((d - 1) (\tilde{r}^{2} - 1) \tilde{r} z^{'} + (d - 2) z)

+ v^{'} ((\tilde{r}^{2} - 1) z^{'} ((\tilde{r}^{2} - 1) (2 z^{'} ((d - 1) (\tilde{r}^{2} - 1) \tilde{r} z^{'} + (d - 2) z) - \tilde{r} (\tilde{r}^{2} - 1) z z^{''}) - (d - 1) h L^{2} \tilde{r})

+ h L^{2} z (- d + 2 \tilde{r}^{2} + 2))

+ 3 h (\tilde{r}^{2} - 1)^{2} v^{'2} z^{'} ((d - 1) (\tilde{r}^{2} - 1) \tilde{r} z^{'} + (d - 2) z) - (d - 1) L^{2} \tilde{r} (\tilde{r}^{2} - 1) z^{'2}

+ L^{2} z (- d + 2 \tilde{r}^{2} + 2) z^{'} + v^{''} (h L^{2} \tilde{r} (\tilde{r}^{2} - 1) z + \tilde{r} (\tilde{r}^{2} - 1)^{3} z z^{'2}) + L^{2} (\tilde{r}^{2} - 1) \tilde{r} z z^{''} = 0,

(d - 1) h^{2} \tilde{r} (\tilde{r}^{2} - 1)^{4} v^{'4} - 2 (d - 1) h L^{2} \tilde{r} (\tilde{r}^{2} - 1)^{2} v^{'2}

(d - 1) h^{2} \tilde{r} (\tilde{r}^{2} - 1)^{4} v^{'4} - 2 (d - 1) h L^{2} \tilde{r} (\tilde{r}^{2} - 1)^{2} v^{'2}

+ z^{'} (3 (d - 1) h \tilde{r} (\tilde{r}^{2} - 1)^{4} (v^{'})^{3} - 3 (d - 1) L^{2} \tilde{r} (\tilde{r}^{2} - 1)^{2} v^{'}

+ z (\tilde{r} (\tilde{r}^{2} - 1)^{4} v^{'} v^{''} - 2 (d - 2) (\tilde{r}^{2} - 1)^{3} v^{'2}))

+ z ((2 - d) h (\tilde{r}^{2} - 1)^{3} v^{'3} + L^{2} (\tilde{r}^{2} - 1) (d - 2 \tilde{r}^{2} - 2) v^{'} - L^{2} \tilde{r} (\tilde{r}^{2} - 1)^{2} v^{''})

+ (d - 1) L^{4} \tilde{r} + 2 (d - 1) \tilde{r} (\tilde{r}^{2} - 1)^{4} v^{'2} z^{'2} - \tilde{r} (\tilde{r}^{2} - 1)^{4} z v^{'2} z^{''} = 0.

\hat{S} \equiv \frac{4 G _{N} ( S _{V ai d y a} - S _{A d S} )}{V _{A}},

\hat{S} \equiv \frac{4 G _{N} ( S _{V ai d y a} - S _{A d S} )}{V _{A}},

d s^{2} =

d s^{2} =

C = \frac{L ^{2 d - 2} Ω _{d - 2}}{8 π G _{N} L} \int_{0}^{z_{*}} d z \int_{0}^{\tilde{r} (z)} d \tilde{r} N P^{d - 2} z^{- d},

C = \frac{L ^{2 d - 2} Ω _{d - 2}}{8 π G _{N} L} \int_{0}^{z_{*}} d z \int_{0}^{\tilde{r} (z)} d \tilde{r} N P^{d - 2} z^{- d},

N \equiv - [\frac{L ^{2}}{( 1 - r ~ ^{2} ) ^{2}} - h v_{\tilde{r}}^{2}] (h v_{z} + 2) v_{z} - (h v_{z} + 1)^{2} v_{\tilde{r}}^{2}, P \equiv \frac{r ~}{1 - r ~ ^{2}},

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

††institutetext: Institute for Theoretical Physics & Cosmology, Zhejiang University of Technology, Hangzhou 310023, China

Subregion complexity in holographic thermalization with dS boundary

Shao-Jun Zhang

[email protected]

Abstract

We study the time evolution of holographic subregion complexity (HSC) in Vaidya spacetime with dS boundary. The subregion on the boundary is chosen to be a sphere within the cosmological horizon. It is found that the behaviour of HSC is similar to that in cases with flat boundary. The whole evolution can be divided into four stages: First, it grows almost linearly, then the growth slows down; After reaching a maximum it drops down quickly and gets to saturation finally. The linear growth rate in the first stage is found to depend almost only on the the mass parameter. As the subregion size approaches the cosmological horizon, this stage is expected to last forever with the subsequent three stages washed out. The saturation time $\tilde{t}_{sat}$ depends almost only on the subregion size $\tilde{R}$ as $\tilde{t}_{sat}=\tanh^{-1}(\tilde{R})$ which is linear in $\tilde{R}$ when $\tilde{R}$ is small but logarithmically divergent as $\tilde{R}$ approaches the cosmological horizon.

1 Introduction

In past decades, with the idea of AdS/CFT or the more generic holographic principle Maldacena:1997re ; Gubser:1998bc ; Witten:1998qj , physicists are trying to build a bridge connecting gravity and other areas of modern theoretical physics, such as condensed matter physics (CMT) Hartnoll:2009sz ; Herzog:2009xv ; McGreevy:2009xe ; Horowitz:2010gk ; Cai:2015cya , QCD Mateos:2007ay ; Gubser:2009md ; CasalderreySolana:2011us , cosmology Banks:2004eb , quantum information theory (QIT) Swingle:2009bg ; Swingle:2012wq ; Qi:2013caa and etc. It is hoped that this bridge may help us get insights into both the strongly coupled problems in the quantum field theory (QFT) side as well as the origin of spacetime in the gravity side. After decades’ efforts, several precise correspondences between the two sides are proposed. Recently, Susskind and his collaborators conjecture that complexity of the boundary QFT may be related to the interior geometry of black hole in the gravity side Susskind:2014rva . In QFT (or QIT), complexity of a target state is an important concept defined as the minimum number of unitary operators (or gates) needed to prepare the state starting from some reference state (for example the vacuum). So far, this conjecture has been refined into two concrete proposals, namely the CV (complexity=volume) and CA (complexity=action) conjectures. In the CV conjecture, complexity of a state living on a time slice $\Sigma$ of the boundary equals to the extremal volume of a codimension-one hypersurface ${\cal B}$ in the bulk ending on $\Sigma$ at the boundary Stanford:2014jda , that is

[TABLE]

where $G_{N}$ is the gravitational constant and $R$ is some typical length scale of the bulk geometry, for example the AdS radius or the horizon radius. While the CA conjecture states that complexity of a state equals to the on-shell gravitational action evaluated on the so-called Wheeler-DeWitt (WDW) patch of the bulk spacetime Brown:2015bva ; Brown:2015lvg . Each conjecture has its own merits and demerits respectively Hashimoto:2018bmb . Inspired by these ideas, an amount of work are raised to study the holographic complexity for various gravity models to check these proposals Momeni:2016ekm ; Cai:2016xho ; Brown:2016wib ; Couch:2016exn ; Yang:2016awy ; Chapman:2016hwi ; Carmi:2016wjl ; Pan:2016ecg ; Brown:2017jil ; Kim:2017lrw ; Cai:2017sjv ; Alishahiha:2017hwg ; Bakhshaei:2017qud ; Tao:2017fsy ; Guo:2017rul ; Zangeneh:2017tub ; Alishahiha:2017cuk ; Abad:2017cgl ; Reynolds:2017lwq ; Hashimoto:2017fga ; Nagasaki:2017kqe ; Miao:2017quj ; Ge:2017rak ; Ghodrati:2017roz ; Qaemmaqami:2017lzs ; Carmi:2017jqz ; Kim:2017qrq ; Cottrell:2017ayj ; Sebastiani:2017rxr ; Moosa:2017yvt ; HosseiniMansoori:2017tsm ; Zhang:2017nth ; Reynolds:2017jfs ; Chapman:2018dem ; Chapman:2018lsv ; Khan:2018rzm ; Caputa:2018kdj ; Feng:2018sqm ; Liu:2019smx ; Jiang:2019qea ; Jiang:2019pgc ; Jiang:2018sqj ; Jiang:2018gft .

The above two conjectures are for the whole boundary system which both are then extended to be defined on subsystem respectively in Refs. Alishahiha:2015rta and Carmi:2016wjl later, and they are now called holographic subregion complexity (HSC). The subregion CV proposal is a natural extension of the well-known Hubney-Rangamani-Takayanagi (HRT) holographic entanglement entropy (HEE) conjecture Ryu:2006bv ; Hubeny:2007xt . Namely, complexity of a subregion ${\cal A}$ of the boundary system equals to the volume of the extremal codimension-one hypersurface $\Gamma_{\cal A}$ enclosed by ${\cal A}$ and the corresponding Hubney-Rangamani-Takayanagi (HRT) surface $\gamma_{\cal A}$ Ryu:2006bv ; Hubeny:2007xt , that is

[TABLE]

where $L$ is the AdS radius. Later studies suggest that it should be dual to the fidelity susceptibility in QIT Alishahiha:2015rta ; MIyaji:2015mia . While in the subregion CA proposal, complexity of subregion ${\cal A}$ is given by the on-shell gravitational action evaluated on the intersection region between WDW patch and the so-called entanglement wedge Czech:2012bh ; Headrick:2014cta . Also, lots of work and effort have been devoted to understand the holographic subregion complexity Caputa2017 ; Caputa2017b ; Czech1706 ; subBenAmi2016 ; subRoy2017 ; subBanerjee2017 ; subBakhshaei2017 ; subSarkar2017 ; subZangeneh2017 ; subMomeni2017 ; subRoy2017b ; subCarmi2017 ; Chen:2018mcc ; Ageev:2018nye ; Ghodrati:2018hss ; Zhang:2018qnt ; Alishahiha:2018lfv ; Alishahiha:2018tep .

On the other hand, in the so-called ”holographic thermalization” topic, the AdS/CFT duality has been applied successfully to study the physics in non-equilibrium processes, especially the thermalization process of hot QCD matter which is strongly coupled and produced in heavy ion collisions at the Relativistic Heavy Ion Collider Gelis:2011xw ; Iancu:2012xa ; Muller:2011ra ; CasalderreySolana:2011us . According to the AdS/CFT dictionary, the thermalization process in the boundary QFT system is dual to a black hole formation process in the bulk which can be modelled simply by a Vaidya-like metric. There are already lots of work on this topic and many interesting results are obtained. For more details on this topic, please refer to the review Balasubramanian:2011ur and references therein.

Complexity in the holographic thermalization process is also studied to investigate its time evolution behaviours under thermal quench. In Refs. Chapman:2018dem ; Chapman:2018lsv , by applying the CV and CA conjectures, it is found that the late time growth of holographic complexity in the Vaidya spacetime is the same as that found for an eternal black hole. In Ref. Chen:2018mcc , the time evolution of subregion complexity is studied in the process with the subregion CV conjecture. And the results show that the subregion complexity is not always a monotonically increasing function of time. Actually, it increases at early time, but after reaching a maximum it decreases quickly and gets to saturation finally. For other related work, please see Refs. Ageev:2018nye ; Ageev:2019fxn ; Ling:2018xpc ; Fan:2018xwf ; Jiang:2018tlu

However, it should be noted that the boundary QFTs considered in the above mentioned work are usually living on the flat Minkowski spacetimes. It would be interesting to generalize the discussions to more realistic situations where QFTs lives on curved spacetimes, which may hep us to understand the extremely hot and condensed physics such as in the very early universe. Several holographic models of the quantum field theory in curved spacetimes (QFTCS) have already been proposed in de Sitter (dS) spacetime and other cosmological backgrounds (please refer to the review Marolf:2013ioa for details). Here we would like to mention the work done in Ref. Marolf:2010tg , where an interesting holographic model was built to relate the QFTs living on the dS boundary to the bulk Einstein gravity. Employing this model, the thermalization process of QFTs in dS spacetime is studied holographically in Ref. Fischler:2013fba . By applying the holographic entanglement entropy as a probe, the whole thermalization is found to be similar to the flat boundary case Liu:2013iza ; Liu:2013qca and can be divided into a sequence of processes. Moreover, the saturation time is found to depend almost only on the entanglement sphere radius. When the radius is small, the saturation time is almost a linear increasing function of the radius, as expected to coincide with the result of the flat boundary case at this time Balasubramanian:2011ur . However, when the radius becomes larger and larger to approach the cosmological horizon, the saturation time blows up logarithmically. Later, the study is extended to include the effect of higher-derivative terms, such as the Gauss-Bonnet correction Zhang:2014cga . And it is found that increasing the Gauss-Bonnet coupling will shorten the saturation time. Please also refer to Refs. Fischler:2014ama ; Fischler:2014tka ; Nguyen:2017ggc for other related work on AdS/CFT with dS boundary.

As the deep connection between holographic entanglement entropy (HEE) and holographic subregion CV (HSCV), it would be interesting to study the time evolution of subregion complexity in the thermalization process of the QFTs living on dS spacetime within the above model. It is natural to ask the following questions: How the existence of the cosmological horizon affects the behaviour of HSCV? Whether the time evolution behaviours of HSCV can be used to describe the the whole thermalization process? Is there any difference between behaviours of HSCV and HEE? The main goal of this work is trying to address these questions.

The work is organized as follows. In the next section, we will give a brief review of the holographic model of QFTs in dS spacetime proposed in Ref. Marolf:2010tg , including the Vaidya-like solution. Then in Sec. III, we study in detail the time evolution of HSCV in the thermalization process. Due to the complication of the equations needed to solve, we rely mainly on numerical calculations. The final section is devoted to discussions and summary.

2 Gravity solutions with dS boundary

In this section, following Refs. Fischler:2013fba ; Zhang:2014cga , we will briefly review the bulk solutions in Einstein gravity with a foliation such that the boundary metric corresponds to a de Sitter spacetime. Three relevant solutions will be presented, including a vacuum AdS, a static AdS black hole and its Vaidya-like cousin.

2.1 Action

We consider $(d+1)$ -dimensional Einstein-Hilbert action as follows

[TABLE]

where $G_{N}$ is the Newton constant and $\Lambda$ negative cosmological constant. The action gives the following equations of motion

[TABLE]

For asymptotically AdS spacetime, the metric can be written in the Fefferman-Graham form Fefferman:1985

[TABLE]

where $L$ is the AdS raduis related to the cosmological constant as $\Lambda=-\frac{d(d-1)}{2L^{2}}$ . The dual quantum field theory lives at the conformal boundary $z=0$ with a metric $ds^{2}_{\Sigma}=g_{\mu\nu}(0,x)dx^{\mu}dx^{\nu}$ . In this paper, we are interested in cases where the boundary metric $ds^{2}_{\Sigma}$ corresponds to a dS spacetime in certain coordinates.

2.2 AdS vacuum solution

The equations of motion (4) admit an AdS vacuum solution as

[TABLE]

The conformal boundary locates at $z=0$ with conformally reduced metric

[TABLE]

which is just the dS spacetime in the static patch with a cosmological horizon at $r=1/H$ , where $H$ denotes the Hubble constant.

The AdS vacuum solution is dual to the vacuum state of the dual QFT with the latter can be taken as the well-known Bunch-Davis or Euclidean vacuum. For a geodesic observer sitting at $r=0$ , the Bunch-Davis vacuum appears to have temperature $T_{dS}=H/2\pi$ natural for the existence of the cosmological horizon.

2.3 AdS black hole solution

The equations of motion (4) also admit an AdS black hole solution with the dS boundary (7)

[TABLE]

The event horizon $z_{+}$ is given by the largest positive root of $h(z)$ . The mass parameter $m$ can be written in terms of the horizon as

[TABLE]

The Hawking temperature of the black hole is

[TABLE]

It should be noted that the zero temperature limit of the black hole solution (2.3) is the not the solution with $m=0$ which is isometric to the AdS vacuum solution (2.2). Actually, the zero temperature limit of the solution has the smallest horizon radius and most ”negative” mass as

[TABLE]

This means that when the mass is negative in the range $0>m>m^{\textrm{e}xt}$ , the black hole still has a regular horizon and reasonable thermodynamics. This is a typical behavior of topological black holes.

Holographically, the black hole solution is dual to the QFT on the static patch of dS spacetime at the temperature given by Eq. (10). Note that this temperature does not have to be the same as the dS temperature $T_{dS}$ . For more discussions on this point, please refer to Ref. Fischler:2013fba .

2.4 Vaidya-like solution

Our aim is to study the holographic thermalization process of the dual QFT under quench. This process can be simply described holographically by a Vaidya-like geometry in the bulk.

Going to the Eddington-Finskelstein coordinates and introducing a time-dependent mass parameter, from the black hole solution one can obtain its Vaidya-like cousin as

[TABLE]

External source should be introduced to maintain the equations of motion

[TABLE]

which implies that the infalling shell is made of null dust. We take the form of the mass function as

[TABLE]

where $M>0$ is the total mass of the dust shell and $v_{0}$ of its thickness. Then the solution describes the collapsing of the null dust shell from the boundary into the bulk to form a black hole. At the QFT side, it corresponds to a sudden global injection of energy into the system and then let it evolve from the Bunch-Davis vacuum to a thermal state with $T>T_{dS}$ .

3 Holographic entanglement entropy and subregion complexity

In this section, by applying the holographic subregion CV (HSCV) (2), we will study the time evolution of holographic subregion complexity in the thermalization process which is described by the Vaidya-like geometry holographically.

On the boundary at time $\tilde{t}$ , taking into account the symmetry of the Vaidya-like metric (2.4), it is convenient to choose the subregion ${\cal A}$ to be a $(d-1)$ -dimensional sphere centred at $\tilde{r}=0$ ( $\tilde{r}\equiv Hr$ ) with raduis $\tilde{R}$ . According to the conjecture (2), the holographic subregion complexity of ${\cal A}$ is given by the extreme volume of the codimension-one hypersurface $\Gamma_{\cal A}$ enclosed by ${\cal A}$ and its corresponding HRT surface $\gamma_{\cal A}$ . So, to study the holographic subregion complexity, we should first find the HRT surface $\gamma_{\cal A}$ whose area gives the holographic entanglement entropy.

3.1 Holographic entanglement entropy

Considering the symmetry, the HRT surface $\gamma_{\cal A}$ in the bulk can be parameterized by functions $z(\tilde{r})$ and $v(\tilde{r})$ , with the boundary conditions

[TABLE]

where $\epsilon$ is an UV cutoff constant. At the tip of the HRT surface, taking into account the symmetry, we have

[TABLE]

where $(z_{\ast},v_{\ast})$ are two parameters labelling the location of the tip and the prime denotes derivative with respect to $\tilde{r}$ . The induced metric on $\gamma_{\cal A}$ is

[TABLE]

The holographic entanglement entropy functional is given by the area of the HRT surface

[TABLE]

To find the extreme value of this functional, we need to solve the two equations of motion, which can be obtained by varying the functional and are rather complicated

[TABLE]

To avoid symbol confusion, we denote the solution of the above two equations as $(v_{0}(\tilde{r}),z_{0}(\tilde{r}))$ which parameterize the HRT surface. The relation between $v$ and $z$ on the HRT surface, denoted as $v_{0}(z)$ , can be obtained by eliminating the parameter $\tilde{r}$ from the two functions.

Generally, the HEE (18) is ultra-divergent. To remove the divergence and for convenience, we define a normalised HEE as

[TABLE]

where ${\cal S}_{AdS}$ is the HEE for the same subregion ${\cal A}$ in pure AdS geometry. And $V_{\cal A}\equiv L^{d-1}\Omega_{d-2}\int_{0}^{\tilde{R}}\frac{\tilde{r}^{d-2}}{(1-\tilde{r}^{2})^{d/2}}d\tilde{r}=L^{d-1}\Omega_{d-2}\frac{\tilde{R}^{d-1}}{d-1}\ _{2}F_{1}\left(\frac{d-1}{2},\frac{d}{2},\frac{d+1}{2},\tilde{R}^{2}\right)$ is the volume of the subregion ${\cal A}$ 111Actually, it should be noted that $V_{\cal A}$ is divergent as $\tilde{R}$ approaches the cosmological horizon to cover the whole boundary space.. So, $\hat{\cal S}$ can be seen as a normalised entanglement entropy density.

3.2 Holographic subregion complexity

Due to the spherical symmetry, the co-dimension one extreme hypersurface $\Gamma_{\cal A}$ , enclosed by ${\cal A}$ and the HRT surface $\gamma_{\cal A}$ , can be parameterized by function $v=v(z,\tilde{r})$ . The induced metric on $\Gamma_{\cal A}$ is

[TABLE]

According to the HSCV proposal (2), the holographic subregion complexity functional of $\Gamma_{\cal A}$ is

[TABLE]

where $v_{z}\equiv\frac{\partial v}{\partial z}$ and $v_{\tilde{r}}\equiv\frac{\partial v}{\partial\tilde{r}}$ . To extremizing the HSCV functional, we need to solve the equation of motion which can be obtained by varying the functional with respect to $v(z,\tilde{r})$

[TABLE]

At first glance, it seems difficult to solve the above equation. However, it is interesting to note that $v(z,\tilde{r})=v_{0}(z)$ is just the solution (Here we would like to emphasize again that $v_{0}(z)$ is just the function giving the relation between $v$ and $z$ on the HRT surface), which can be checked directly by plugging $v_{0}(z)$ into the equation. This simply means that $\Gamma_{\cal A}$ is just formed by dragging the HRT surface $\gamma_{\cal A}$ along the $\tilde{r}$ direction. Similar feature has already been observed in flat boundary case with strip subregion in Ref. Chen:2018mcc .

As HEE, HSC is also ultra-divergent, so we can also define a normalised HSC density as

[TABLE]

where ${\cal C}_{AdS}$ is the HSCV for ${\cal A}$ in pure AdS geometry.

3.3 Numerical results

Having set up the general frame work of HEE and HSCV, now we are ready to study the time evolution of HSC in holographic thermalization. Due to the complication of the equations needed to solve, we rely on numerical method. And for convenience, we set the AdS radius $L=1$ .

3.3.1 General behaviours

In Fig. 1, we plot the time evolution of normalised HSC density $\hat{\cal C}$ for various $\tilde{R}$ with fixed spacetime dimension. From the figure, one can see that the time evolution of $\hat{\cal C}$ is not a monotonically increasing function of the time. Rather, it can be divided into four stages: After quench, firstly it grows quickly and almost linearly, then the growth slows down; After reaching a maximal value $\hat{\cal C}_{max}$ it starts to drop down fast,and shortly after the drop down stops and it saturates to a constant value $\hat{\cal C}_{sat}$ finally. Moreover, it is interesting to note that the final saturation constant may be negative, which means that the final value of the complexity may be smaller than its initial value. These behaviours are very different from that in CV or CA conjectures, where the complexity is always a monotonically increasing function of time Chapman:2018dem ; Chapman:2018lsv . Similar behaviours have been observed in flat boundary cases with strip subregion Chen:2018mcc ; Ling:2018xpc , indicating universality of the behaviours.

From Fig. 1-3, one can see that the maximal value $\hat{\cal C}_{max}$ depends on the subregion size $\tilde{R}$ , the spacetime dimension $(d+1)$ and the mass parameter $M$ . Increasing $\tilde{R}$ or $M$ will yield a bigger $\hat{\cal C}_{max}$ , while increasing the dimension will, on the contrary, lower the maximal value.

Moreover, one can also see that the final saturation constant $\hat{\cal C}_{sat}$ also depends on $(\tilde{R},d,M)$ but in a more complicated way.

3.3.2 Linear growth stage

Let us focus on discussing the first stage when $\hat{\cal C}$ grows almost linearly in time, i.e.,

[TABLE]

where $A$ is the proportional constant which may depend on $(\tilde{R},d,M)$ . From Fig. 1-3, one can see that $A$ is nearly independent of $\tilde{R}$ and $d$ ; While it strongly depends on $M$ . By fitting the numerical data, it is found that $A\approx 0.4M$ .

From Fig. 1, one can also see that larger the subregion size $\tilde{R}$ is, longer time the linear growth stage lasts. It is expected that as $\tilde{R}$ approaches the cosmological horizon $\tilde{r}=1$ to cover the entire boundary space, the linear growth stage will last forever which agrees well with the CV conjecture. We can see this point more clearly in Fig. 4 where we take $d=2$ case as an example. We will give more evidences on this point later.

3.3.3 Saturation time

In Refs. Fischler:2013fba ; Zangeneh:2017tub , one defines the saturation time as the time HEE approaches a constant. Similarly, for the complexity, we can also define a saturation time $\tilde{t}_{sat}$ as the time $\hat{\cal C}$ reaches its saturation constant $\hat{\cal C}_{sat}$ .

In Fig. 5, we plot the time evolution of the two observables, $\hat{\cal S}$ and $\hat{\cal C}$ to make a comparison. From the figure, we can see the well-known fact that $\hat{\cal S}$ is always a monotonically increasing function of time. Moreover, from Fig. 5 and Fig. 2, one can see that $\hat{\cal S}$ and $\hat{\cal C}$ reache their saturation values at almost the same time. And their saturation time $\tilde{t}_{sat}$ is nearly independent of $d$ and $M$ .

In Fig. 6, the saturation time $\tilde{t}_{sat}$ for $\hat{\cal C}$ as a function of the subregion size $\tilde{R}$ is plotted. The numerical results can be well fitted by the function $\tilde{t}_{sat}=\tanh^{-1}(\tilde{R})$ , as for the HEE Fischler:2013fba . It is interesting to note that the $\tilde{t}_{sat}$ is just the time light takes travelling from the origin $\tilde{r}=0$ to the boundary of the subregion $\tilde{r}=\tilde{R}$ .222We thank J.F. Pedraza for pointing out this. From the figure and the fitting, one can easily see that $\tilde{t}_{sat}$ is linear in $\tilde{R}$ when $\tilde{R}$ is small; However, as $\tilde{R}$ approaches the cosmological horizon $\tilde{r}=1$ , the saturation time diverges logarithmically and thus the linear growth stage will also last forever, as we already mentioned above.

4 Summary and Discussions

In this work, we consider the holographic model of thermalization process for QFTs in dS spacetime. By applying the holographic subregion CV conjecture, we study the time evolution of subregion complexity under quench. The subregion ${\cal A}$ is chosen to be a sphere on the boundary time slice. The dual extremal codimension-one hypersurface $\Gamma_{{\cal A}}$ in the bulk, whose volume gives the complexity of ${\cal A}$ , is found to be simply swept out by the HRT surface along the $\tilde{r}$ -direction. The whole time evolution of subregion complexity can be divided into four stages: It first increases almost linearly; Then its growth slows down and after reaching a maximum it starts to drop down quickly, and shortly after the drop down stops and it gets to saturation finally. This picture is similar to that in flat boundary cases but with a strip subregion Chapman:2018lsv ; Ling:2018xpc . This implies that the time evolution behaviours of subregion complexity are very general, and is independent of the subregion shape and the cosmological horizon.

The linear growth rate in the first stage is found to almost only depend on the mass parameter. As the subregion size approaches the cosmological horizon, this stage is expected to last forever, and as the HEE the saturation time is logarithmically divergent. The saturation time is found to depend almost only on the subregion size $\tilde{R}$ , and their relation can be well fitted by the function $\tilde{t}_{sat}=\tanh^{-1}(\tilde{R})$ . It is interesting to note that the $\tilde{t}_{sat}$ is just the time light takes travelling from the origin $\tilde{r}=0$ to the boundary of the subregion $\tilde{r}=\tilde{R}$ . The underlying physical meaning of this fact needs further investigation.

In this work, we only consider the HSCV conjecture. It is interesting to check whether general behaviours of subregion complexity still holds for other conjectures, for example the holographic subregion CA. In Ref. Zhang:2014cga , using HEE as a probe we show that including the Gauss-Bonnet correction will shorten the saturation time. It is also interesting to see how the higher-derivative terms affect the time evolution of subregion complexity. We leave these questions for further investigations.

Acknowledgement

This work was supported by National Natural Science Foundation of China (Nos. 11605155 and 11675144). We thank J.F. Pedraza for helpful comments on this manuscript.

Bibliography105

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1(1) J. M. Maldacena, The Large N limit of superconformal field theories and supergravity , Int. J. Theor. Phys. 38 , 1113 (1999) [Adv. Theor. Math. Phys. 2 , 231 (1998)] [hep-th/9711200].
2(2) S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Gauge theory correlators from noncritical string theory , Phys. Lett. B 428 , 105 (1998) [hep-th/9802109].
3(3) E. Witten, Anti-de Sitter space and holography , Adv. Theor. Math. Phys. 2 , 253 (1998) [hep-th/9802150].
4(4) S. A. Hartnoll, Lectures on holographic methods for condensed matter physics , Class. Quant. Grav. 26 , 224002 (2009) [ar Xiv:0903.3246 [hep-th]].
5(5) C. P. Herzog, Lectures on Holographic Superfluidity and Superconductivity , J. Phys. A 42 , 343001 (2009) [ar Xiv:0904.1975 [hep-th]].
6(6) J. Mc Greevy, Holographic duality with a view toward many-body physics , Adv. High Energy Phys. 2010 , 723105 (2010) [ar Xiv:0909.0518 [hep-th]].
7(7) G. T. Horowitz, Introduction to Holographic Superconductors , Lect. Notes Phys. 828 , 313 (2011) [ar Xiv:1002.1722 [hep-th]].
8(8) R. G. Cai, L. Li, L. F. Li and R. Q. Yang, Introduction to Holographic Superconductor Models , Sci. China Phys. Mech. Astron. 58 , no. 6, 060401 (2015) [ar Xiv:1502.00437 [hep-th]].