Charged BTZ black holes cannot be destroyed via the new version of the gedanken experiment

The singularity at the center of charged Bañados–Teitelboim–Zanelli (BTZ) black holes is called a conical singularity. Unlike the canonical singularity in typical black holes, a conical singularity does not destroy the causality of spacetime. Due to the special property of the conical singularity, we examine the weak cosmic censorship conjecture (WCCC) using the new version of the gedanken experiment proposed by Sorce and Wald. A perturbation process wherein the spherically symmetric matter fields pass through the event horizon and fall into the black holes is considered. Assuming that the cosmological constant is obtained by the matter fields, it therefore can be seen as a dynamical variable during the process. From this perspective, according to the stability condition and the null energy condition, the first- and second-order perturbation inequalities are derived. Based on the first-order optimal condition and the second-order perturbation inequality, we show that the nearly extremal charged BTZ black hole cannot be destroyed in the above perturbation process. The result also implies that even if the singularity at the center of the black hole is conical, it still should be surrounded by the event horizon and hidden inside the black hole.


Introduction
A gravitational singularity at the center of black holes is usually mathematically ill-defined because it can cause the curvature to diverge. Since the divergence can affect the validity of the causality law of spacetime, the singularity should always be hidden inside the black hole and is not allowed to appear in spacetime. In other words, an observer at infinity cannot detect any information from the singularity. To ensure well-define spacetime and the validity of the causality law, a e-mail: xinyangwang@bnu.edu.cn b e-mail: jiejiang@mail.bnu.edu.cn (corresponding author) Penrose [1] proposed the weak cosmic censorship conjecture (WCCC). This conjecture postulates that the singularity should be surrounded by the event horizon and can never be exposed to spacetime.
Although the WCCC is suggested for any black hole, the general demonstration for the WCCC is notoriously difficult. Moreover, the validity of the conjecture also depends on the research technique. A notable gedanken experiment attempting to destroy Kerr black holes was proposed by Wald [2]. For Kerr black holes, the angular momentum should be bounded by the mass as J ≤ M 2 when the event horizon of the black hole is guaranteed. To destroy Kerr black holes, a process whereby a test particle with sufficient angular momentum falls into Kerr black holes is considered. If the particle successfully drops into the black hole, the value of the angular momentum can exceed the boundary, and the existence of the event horizon will no longer be guaranteed. In this situation, a naked singularity will emerge in spacetime, and the WCCC for Kerr black holes is no longer valid. Fortunately, it was shown that when the particle carries the angular momentum sufficient to over-spin the Kerr black hole, it cannot be captured by the black hole due to the centrifugal force. This means that the WCCC for Kerr black holes cannot be violated under this process. Since this method was introduced, it has been used to examine the WCCC for other kinds of black holes [3][4][5][6][7]. However, the method has an inherent defect because the interaction between the particle and the background spacetime is not considered. Moreover, Hubeny [8] proposed that using the above method, violation of the WCCC might occur if one suitably adds the particle to a slightly non-extremal black hole. Therefore, to solve the defect, Sorce and Wald [9] proposed a new version of the gedanken experiment to examine the WCCC for Kerr-Newman black holes. In this experiment, the interaction between the black hole spacetime and matter fields outside the black hole was sufficiently considered. In this way, the black hole and matter fields are regarded as a complete dynamical system, while the process by which matter fields pass through the event horizon to perturb the black hole can be treated as a dynamical evolution process of the system. From this perspective, it was shown that the WCCC for near-extremal Kerr-Newman black holes cannot be as easily violated under the second-order approximation of the matter field perturbation. Furthermore, based on the new version of the gedanken experiment, the WCCC has also been demonstrated to be valid for other kinds of black holes [10][11][12][13].
For three-dimensional Bañados-Teitelboim-Zanelli (BTZ) black holes, the singularity at the center of the black hole is called a conical singularity. However, a conical singularity is not similar to a canonical singularity because it does not cause the spacetime curvature to diverge [14]. This means that a conical singularity does not influence the well-defined spacetime and the validity of the causality law. Moreover, an assumption was proposed in Ref. [15] which states that if black holes do not have a canonical singularity, the black hole will not need to obey the requirement of the WCCC. Due to the specific property of the conical singularity, investigating the WCCC for charged BTZ black holes is the best way to examine the assumption. Therefore, according to Ref. [14], the WCCC for rotated BTZ black holes was examined using the new version of the gedanken experiment. It was shown that after perturbation of the matter fields, the WCCC for rotated BTZ black holes cannot be violated. However, in this investigation, the Maxwell field is not contained in the spacetime, and the cosmological constant is not treated as a dynamical variable. Therefore, following our previous work [16], the cosmological constant is regarded as a portion of matter fields, while black holes and matter fields can be seen as a complete dynamical system using the new version of the gedanken experiment. In this case, the cosmological constant should be regarded as a dynamical variable. From this perspective, we will examine the WCCC for charged BTZ black holes under the second-order approximation of the perturbation to comprehensively check the assumption.
The paper is organized as follows. In Sect. 2, we discuss the spacetime geometry of charged BTZ black holes under the perturbation of matter fields. In Sect. 3, based on the Iyer-Wald formalism, we derive the first-order and the second-order perturbation inequalities of the black hole. In Sect. 4, using the first-order optimal option and the secondorder inequality, we examine the WCCC for charged BTZ black holes under the second-order approximation of the perturbation. The paper ends with conclusions in Sect. 5.

Perturbed geometry of charged BTZ black holes
For the three-dimensional Einstein-Maxwell-AdS gravitational theory, the Lagrangian three-form is where R is the Ricci scalar, is the cosmological constant with a negative value, F = d A is the strength of the electromagnetic field, A is the gauge potential of the electromagnetic field, and is the volume element of the three-dimensional spacetime. From the Lagrangian, a class of static spherically symmetric solutions describing charged BTZ black holes is given as where the blackening factor f (r ) is The parameters M and Q in f (r ) correspond to the mass and the electric charge of the black hole. The radius of the event horizon r h is the largest root of the equation f (r ) = 0. According to the radius of the event horizon, the area, surface gravity, and electric potential for the event horizon are further given as Subsequently, we should first consider a process to examine the WCCC of charged BTZ black holes. In this process, matter fields pass through the event horizon and fall into the black hole. We further suggest that the cosmological constant can be regarded as an effective parameter determined by the matter source coupling to the Einstein-Maxwell gravity. This implies that the cosmological constant can be treated as a portion of matter fields. Moreover, the spacetime geometry of the black hole and matter fields can be seen as a complete dynamical system. From Eq. (1), the Lagrangian of this dynamical system can be rewritten as where L mt is the Lagrangian of matter fields. From the Lagrangian, the Einstein-Maxwell-AdS gravitational theory can be efficiently derived again when there is a static solution of matter fields such that where T ab is the stress-energy tensor of matter fields. Since the cosmological constant is contained in matter fields, when matter fields fall into the black hole, their value can vary with the process. This suggests that the cosmological constant should be considered a dynamical variable in this situation.
To simplify the calculation and discussion, we only consider the case in which the configuration of matter fields is spherically symmetric. When the spacetime of the black hole and matter fields is treated as the complete dynamical system, the configurations of the metric g ab , the gauge potential of the electromagnetic field A, and matter fields can be uniformly described by a symbol α. The variation of the configurations with the process can be labeled a one-parameter family, i.e., α(λ). For the case of λ = 0, α(0) represents the configuration of the dynamical fields on the background spacetime which is a charged BTZ black hole. When the change in λ with the process is small enough, the process can be treated as a perturbation. Therefore, under the perturbation, the equation of motion in this one-parameter family can be written as where j a is the current of the electric charge, and T EM ab is the stress-energy tensor of the electromagnetic field, which can be expressed as The metric of spacetime during the perturbation can be generally written as where μ(v, r, λ) is an arbitrary function. When f (v, r, 0) = f (r ) and μ(v, r, 0) = 1, the metric will degenerate into the case of the background spacetime. Following the train of thought proposed by Sorce and Wald, we also suggest that the perturbation should satisfy two conditions. The first is that the bifurcation surface B is not affected by the perturbation. The second is that the perturbation should satisfy the stability condition. This condition means that after perturbation of the matter fields, the spacetime geometry can also be described by the class of charged BTZ solutions. In other words, the solutions of the equation of motion after the perturbation are still described by Eq. (2), and the parameters M, Q, and are replaced by M(λ), Q(λ), and (λ), respectively, i.e., where the blackening factor f (r, λ) is In addition, when the stability condition is satisfied, the stress-energy tensor of matter fields at sufficiently late time should have a similar form as Eq. (6), i.e., The currents with the energy, the electric charge, and the cosmological constant are supposed to pass through a finite portion of the future event horizon and drop into the black hole to ensure that the two conditions are satisfied. So we can always choose a hypersurface = H ∪ 1 on spacetime. For the hypersurface, H is a part of the event horizon which starts from the bifurcation surface B and continues along the future event horizon until the section of the event horizon B 1 at a sufficiently late time. After that, the hypersurface becomes spacelike and goes to asymptotic infinity through the isochronous surface 1 . Since the hypersurface is independent of the parameter λ, r h is always the radius of the event horizon of the background spacetime. In addition, due to the stability condition, the spacetime geometry on the hypersurface 1 can be described by Eq. (10) directly.

Perturbation inequalities for charged BTZ black holes
Starting with this section, we would like to examine the WCCC for charged BTZ black holes under perturbation of the matter fields. Before examining the WCCC, the firstorder and second-order perturbation inequalities should be derived. According to the Iyer-Wald formalism [17], we focus mainly on the Lagrangian three-form L of threedimensional Einstein-Maxwell gravitational theory and use the off-shell variation to obtain the two perturbation inequalities. For the theory of gravity, since the dynamical fields consist of the Lorentz signature metric g ab and the electromagnetic field A, we also use the unified symbol φ to represent the configuration of the dynamical field operators, i.e., φ = (g ab , A). When we consider the perturbation of the matter fields, the changing behavior of the configuration can also be described by the one-parameter family λ, i.e., φ(λ), while the variation of the quantity η related to the dynamical fields φ is defined by Using the off-shell variation, the first-order variation of the Lagrangian L is formally given as where E φ = 0 is the equation of motion of the on-shell fields related to the Lagrangian L, and is called the symplectic potential two-form which is locally constructed out of φ and its derivatives. In Einstein-Maxwell gravity, from the Lagrangian the equation of motion is formally given as with where T ab and j a are the stress-energy tensor and the electric current of matter fields, respectively. For the background spacetime, the form of the stress-energy tensor is similar to Eq. (6). The total symplectic potential two-form (φ, δφ) in Eq. (14) can be decomposed into two parts which represent the part of gravity and the part of the electromagnetic field, respectively, Using the symplectic potential, the symplectic current two-form can be defined as Since the symplectic potential can be linearly decomposed as the gravitational part and the electromagnetic part, the symplectic current can also be decomposed as the two parts, i.e., The specific expression of the two parts is respectively given as For the gravitational part, we denote with We set ζ a as an infinitesimal generator of the diffeomorphism. Replacing δ with L ζ in Eq. (14), one can define the Noether current two-form J ζ associated with ζ a as In addition, it is shown that the Noether current can also be represented as [2] J where Q ζ is called the Noether charge, and C ζ = ζ · C are the constraints of the theory. If C a = 0, the dynamical fields will satisfy the equations of motion. In the theory of Einstein-Maxwell gravity, the Noether charge can be decomposed into two parts as well that represent the conservation charge of the gravity and the electric charge of the electromagnetic field, i.e., where the specific expression of the Q GR ζ and Q EM ζ can respectively be given as Meanwhile, the constraints C abc in Eq. (26) are defined as To investigate the perturbation that comes from the spherically symmetric matter fields, the diffeomorphism generated by the static Killing vector field ξ a = (∂/∂v) a on the background spacetime is involved. Due to the diffeomorphism, a gauge condition that makes the coordinates (v, r, φ) fixed under the variation can be chosen. This means that under the gauge condition, the Killing vector is invariable under the variation, i.e., δξ a = 0. Therefore, taking the first-order variations of Eqs. (25) and (26), the first-order variational identity can be expressed as Furthermore, taking the variation on the above identity again, the second-order variational identity can also be obtained as where we have used the fact that L ξ φ = 0 for the background dynamical fields.
In the following, according to Eqs. (30) and (31), we will derive the first-order and second-order perturbation inequalities, respectively. Furthermore, based on the two inequalities, we will examine the WCCC for charged BTZ black holes under the second-order approximation of the matter fields perturbation.

The first-order perturbation inequality
We first calculate the integral form of the first-order variational identity to derive the first-order perturbation inequality. Integrating the first-order variational identity on the hypersurface and utilizing the condition L ξ φ = 0, we have Considering the property whereby the hypersurface consists of a portion of the event horizon H and the spacelike hypersurface 1 , and using the Stokes' theorem, the integral form of the first-order variational identity can be decomposed as When the cosmological constant is regarded as a variable, the divergence term will appear as the result of the integral on the spacelike infinity. In order to formally obtain the result of the integral, we will choose a cut-off sphere S c with radius r c to replace the asymptotic infinity boundary of 1 . Moreover, we will take the limitation such that the cut-off sphere S c approaches asymptotic infinity again to obtain the final result. Following a similar consideration as Ref. [9], the gauge condition of the electromagnetic field, ξ a δ A a = 0, on the event horizon H can also be imposed. However, on the hypersurface 1 , the gauge condition cannot be used to calculate quantities related to δ A and δ 2 A because the specific expression of the gauge potential, does not satisfy the gauge condition on the hypersurface 1 . Therefore, we need to find another method to calculate the relevant quantities. Since the strength of electromagnetic F(λ) is gauge-invariant, it is reasonable for us to calculate the quantities which contain the gauge potential on hypersurface 1 using the electromagnetic strength. In the following calculation, we use only the specific expression of F(λ) and neglect the expression of the gauge potential A(λ) on the hypersurface 1 .
The volume element on the hypersurface 1 can be written as Since the expression of does not contain the parameter λ, this implies that the variation of the element volume is vanishing, i.e., δ = 0. According to the equation of motion, the stress-energy tensor T ab and electric current j a on 1 can be expressed as Firstly, we consider the first term in Eq. (33). As mentioned above, the Noether charge and the symplectic potential are both decomposed as the gravity part and the electromagnetic field part. This means that the first term can also be decomposed into the two parts, For the gravity part, from Eqs. (19) and (28), using the specific expression of the metric in Eq. (10), we can obtain For the electromagnetic field part, using Eqs. (19) and (28), the integrand can be written as Substituting the specific expression of the strength of the electromagnetic field into Eq. (39), we find the following relation between the second term and the third term This implies that the last two terms in the integrand can cancel each other exactly, and only the integrand of the electromagnetic part remains as the first term. Hence, the result of the integral can be expressed as Combining Eq. (38) with Eq. (41), the integral result of the first term in Eq. (33) is According to the condition wherein the perturbation cannot influence the bifurcation surface B, the second term in Eq. (33) can be directly neglected.
Secondly, we will calculate the integrals on the hypersurface 1 . For the third term in Eq. (33), substituting Eqs. (16) and (36) into it, we have From the expression of the metric, we can obtain the condition g ab δg ab = 0. This condition implies that the integral of the third term is equal to zero. In order to evaluate the fourth term of Eq. (33), we should calculate the constraints on 1 firstly. According to Eq. (29), the expression of the Killing vector contracting with the constraints on the hypersurface 1 can be specifically written as Taking the variation on the above expression and integrating it on the hypersurface 1 , the fourth term can be obtained as Finally, we turn to calculate the integrals on the event horizon H. For the fifth term in Eq. (33), the integral result is directly equal to zero because the Killing vector contracting with the volume element vanishes on the event horizon H. For the sixth term of Eq. (33), according to the definition of the constraints and j (0) = 0 on the background spacetime, using the expression of the gauge potential A, we have where we have denoted the volume element of the event horizon as˜ = dv ∧ˆ , andˆ = r dϕ is the volume element of a cross-section of the event horizon. According to the electromagnetic part of the equation of motion, we can then obtain Therefore, Eq. (46) can be reduced as Substituting Eqs. (42), (45), and (48) into Eq. (33), the integral form of the first-order variational identity can be written as To derive the first-order perturbation inequality, we should determine the connection between the right-hand side of Eq. (49) and the null energy condition. During the perturbation, since we consider the process wherein the spherically symmetric matter fields fall into charged BTZ black holes, a null vector field can be chosen as where Using the null vector field, the null energy condition during the perturbation process can generally be expressed as It can be demonstrated that the null energy condition should satisfy the following relation Utilizing the fact that β(0) = 0 for the background spacetime, the null energy condition under the first-order approximation can be obtained as Substituting Eq. (54) into Eq. (49), the integral form of the first-order variational identity can be reduced as This inequality is called the first-order perturbation inequality.
Since the main objective of our investigation is to examine the WCCC for charged BTZ black holes under a second-order approximation to check the assumption that if a black hole does not have a canonical singularity, it does not need to satisfy the requirement of the WCCC. When the first-order perturbation inequality is satisfied, this indicates that the WCCC cannot be violated under the first-order approximation of the perturbation, while the higher-order approximation can be largely neglected. However, if a first-order perturbation inequality is chosen as an optimal option, i.e., the WCCC cannot be examined under the first-order approximation, and the second-order approximation needs to be further considered in this situation. In addition, the optimal option also implies that the energy flux through the event horizon vanishes under the first-order approximation of perturbation.

The second-order perturbation inequality
To consider the second-order approximation of the perturbation, we should sequentially derive the second-order perturbation inequality. Integrating the second-order variational identity on the hypersurface , Eq. (31) can be written as According to Stokes' theorem, and utilizing the property of the hypersurface , the integral expression can be decomposed as where In the integral form of the second-order variational identity, we can see that, except for the last two terms, any term in Eq. (58) only takes the variation again on the corresponding term in the first-order variational identity. Therefore, we can use the integral result directly in the first-order variational identity to evaluate the integral of the second-order perturbation inequality.
As with the first-order variational identity, the first term in Eq. (58) can also be decomposed into the gravity part and the electromagnetic field part. From Eqs. (38) and (41), the integrand of the two parts can be obtained as and Integrating the above two equations on the hypersurface S c , respectively, and summing the integral results of the two parts, the result of the first term in Eq. (58) can be given as Based on the definition of the constraints (29) and j (0) = 0 on the background spacetime, the second-order variation of the Killing vector contracting with the constraints can be expressed as Integrating Eq. (64) on the event horizon H, using the result of Eq. (47) and the gauge condition ξ a δ A a = 0, the integral result of the third term in Eq. (58) is From the result of Eq. (44), the fourth term can be directly calculated as Based on the condition g ab (λ)δg ab (λ) = 0 and the fact that ξ a contracting with the volume element is zero on the event horizon H, we have For the seventh term, since the symplectic current can be decomposed as the gravitational part and the electromagnetic part, its integral can be decomposed as the two parts as well, i.e., For the integral of the gravitational part, according to the specific expression of the metric in Eq. (9), the integral can be directly calculated as where we have used the optimal option of the first-order approximation.
To evaluate the integral of the electromagnetic part, the specific expression of the symplectic current of the electromagnetic part should first be considered. Based on the definition, the symplectic current of the electromagnetic field can be expressed as The expressions of the electromagnetic field strength in Eq. (10) indicates that if the index c in the volume element is chosen as the component of r , i.e., (dr ) c , after contracting (dr ) c with F cd , a proportional relation between the result of the contraction and the Killing vector can be obtained as (dr ) c F cd ∝ −ξ d . According to the proportional relation and the gauge condition ξ a δ A a = 0 on the event horizon, we can find that the last two terms in Eq. (70) both vanish, and Eq. (70) can be reduced as For the first term in Eq. (71), it will not appear in the final expression of the integral because it only contains a boundary term after the integration, and the boundary term will not contribute to E H . Using the gauge condition ξ a δ A a = 0 again on the event horizon, the expression of Eq. (68) can finally be simplified as On the other hand, based on the stress-energy tensor of the electromagnetic field, the result of Eq. (72) can be rewritten as Therefore, combining Eqs. (62), (65), (66), (67), and (73), the second-order variational identity can be obtained as Next we will evaluate the last two terms E 1 (φ, δφ) and Y (φ, δφ) in Eq. (74) to obtain the complete expression of the second-order variational identity. From the expressions of E 1 (φ, δφ) and Y (φ, δφ), we can see that the two quantities depend only on the first-order variation of the dynamical fields, and the integral is implemented only on the hypersurface 1 . Following the method in Ref. [9], an auxiliary spacetime should be introduced to calculate the two terms. Since the auxiliary spacetime involves only the dynamical fields and their first-order variation, the first-order variation of the dynamical fields δφ can be labeled δφ BTZ , where φ BTZ represents the dynamical fields in the auxiliary spacetime. Therefore, we can replace δφ with δφ BTZ in the expression of the two quantities, i.e., According to the stability condition, the spacetime geometry still belongs to the class of charged BTZ solutions after the perturbation, and the changing behavior of the configurations of the dynamical fields can be described by the one-parameter family λ. Therefore, after the perturbation process, the solutions of the equation of motion in the auxiliary spacetime can be directly written as where the blackening factor f BTZ (r, λ) is Since only the first-order variation of the dynamical fields is contained in the auxiliary spacetime, and the higher-order variation does not exist, the parameters M BTZ (λ), Q BTZ (λ), and BTZ (λ) under the first-order approximation can be expanded as where δ M, δ Q, and δ are chosen to agree with the values of the first-order approximation. This implies that the relation δ 2 M = δ 2 Q = δ 2 = E H (φ, δφ) = 0 will be given naturally in the auxiliary spacetime. Based on the above discussions, the two terms E 1 (φ, δφ) and Y (φ, δφ) can be directly evaluated in the auxiliary spacetime.
Using Stokes' theorem, the integral of the second-order variational identity can be written as For the third and fourth terms in Eq. (79), because δ 2 M = δ 2 Q = δ 2 = 0, we can easily demonstrate that Therefore, Eq. (79) can be reduced as Using the result of Eq. (62) and the gauge condition ξ a δ A BTZ a = 0 on the event horizon, the first and second terms of Eq. (81) can be obtained as and Substituting Eqs. (82) and (83) into Eq. (81), we have Finally, substituting Eq. (84) into Eq. (74), the complete expression of the second-order variational identity can be obtained as According to Eq. (53) and T ab (0) = 0 for the background spacetime, using the optimal option of the first-order approximation, the null energy condition under the second-order approximation can be given as Substituting Eq. (86) into Eq. (85), the second-order variational identity becomes an inequality, i.e., which is called the second-order perturbation inequality. Thus far, the first-order and second-order perturbation inequalities have both been obtained. In the following, based on the two perturbation inequalities, we will examine the WCCC for charged BTZ black holes under the second-order approximation of the perturbation to check whether the conical singularity should still be surrounded by the event horizon and hidden inside the black hole.

Gedanken experiment to examine nearly extremal charged BTZ black holes
In the following, based on the first-order and second-order perturbation inequalities, we would like to examine the WCCC for charged BTZ black holes under the second-order approximation of the perturbation. According to the stability condition, the spacetime geometry on the hypersurface 1 still belongs to the class of charged BTZ solutions. Therefore, testing the WCCC is equivalent to testing whether the event horizon exists after the perturbation. If the event horizon exists, the geometry on 1 can be described by charged BTZ black holes. This means that the WCCC cannot be violated under the perturbation of the matter fields. Figure 1 shows that the event horizon of the black hole exists when the minimal value of the blackening factor f (r ) is negative. If the sign of the minimal value of f (r ) becomes positive, the event horizon will disappear, while the naked singularity will be exposed to spacetime. This implies that examining the WCCC for charged BTZ black holes is equivalent to checking the sign of the minimal value of f (r, λ) in Eq. (11). Therefore, we define a function to represent the minimal value of the blackening factor, where r m (λ) is defined as the position of the minimal value of the function f (r, λ), and its value can be determined by Meanwhile, Eq. (89) also gives the following relation Considering the second-order approximation of the matter fields perturbation, the expression of h(λ) can be expanded with respect to the parameter λ at λ = 0. The specific expression of h(λ) under the second-order approximation is where we have used Eq. (90) to simplify the expression. Taking the first-order variation of Eq. (89) and using Eq. (90) as well, the first-order variation of r m can be expressed as Substituting Eq. (92) into Eq. (91), the second-order expansion of h(λ) can be simplified as Following a similar idea as in Ref. [9], a small parameter which should agree with the first-order approximation of the matter fields perturbation is introduced. For nearly extremal charged BTZ black holes, the minimal value of the blackening factor r m and the radius of the event horizon r h satisfy the relation r m = (1 − )r h . Using this relation, Eq. (89) can be written as under the first-order approximation of . The minimal value of the blackening factor f (r m , λ) under the second-order approximation of is expressed as Substituting the specific expression of f (r ) into the result of Eq. (95), we have Utilizing the first-order optimal option and the second-order perturbation inequality, combining the above results, the expression of h(λ) under the second-order approximation is reduced as In addition, under the zero-order approximation of , we have the following relation Combining Eq. (98) with Eq. (97), the function h(λ) can finally be written as Equation (99) shows that the value of h(λ) is negative under the second-order approximation. This illustrates that the event horizon still exists and that the WCCC for a charged BTZ black hole cannot be violated during the perturbation of matter fields. This result also implies that even if the singularity at the center of the black hole is conical, it should also be surrounded by the event horizon of the black hole and cannot be exposed to spacetime.

Conclusions
Based on the new version of the gedanken experiment proposed by Sorce and Wald, we examine the WCCC for charged BTZ black holes under a second-order approximation of the perturbation of the matter fields to check the assumption that if a black hole does not have a canonical singularity, it does not need to satisfy the requirement of the WCCC. In our investigation, we consider a process wherein the spherically symmetric matter fields pass through the event horizon to perturb charged BTZ black holes. During the perturbation, the cosmological constant is regarded as a portion of the matter fields, while the spacetime geometry of the black hole and matter fields can be seen as a complete dynamical system. In this situation, the cosmological constant should be considered a dynamical variable. In addition, a stability condition is proposed to examine the WCCC, which states that the spacetime geometry still belongs to the class of charged BTZ solutions after the perturbation. According to the stability condition and the null energy condition, the first-order and second-order perturbation inequalities are derived. Then, using the optimal option of the first-order approximation and the second-order perturbation inequality, based on the stability condition and the null energy condition, the WCCC for charged BTZ black holes under the second-order approximation of the perturbation is examined. The result shows that the event horizon of charged BTZ black holes still exists after the perturbation of the matter fields, and the WCCC for charged BTZ black holes cannot be violated. It also illustrates that even if the conical singularity does not cause the spacetime curvature to diverge, it still should be surrounded by the event horizon and hidden inside the black hole. or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated or analysed during the current study.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .