Weak cosmic censorship in Born–Infeld electrodynamics and bound on charge-to-mass ratio

We construct a class of counterexamples to cosmic censorship in four dimensional Einstein–Born–Infeld theory with asymptotically anti-de Sitter boundary conditions, and investigate the effect of the Born–Infeld parameter b in comparison with the counterpart of Einstein–Maxwell theory. When a charged massive scalar field is included into the action, we find that this class of counterexamples to cosmic censorship would be removed if the charge of scalar fields is above the bound of charge qw\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_{w}$$\end{document}. In particular, the bound of charge required to preserve cosmic censorship increases with the increasing of Born–Infeld parameter. Meanwhile, we also show the lower bounds on charge-to-mass ratio with the different values of Born–Infeld parameter.


Introduction
The study of black hole singularities has been an interesting subject since the original work on the weak cosmic censor- ship conjecture (WCCC) [1,2], which states naked singularities arising in the solutions of Einstein's equations [3] must be hidden within event horizons of black hole, and could not be observed from future null infinity. Although the concrete proof for the validity of WCCC has not been established, there exist lots of works of testing the validity of WCCC through the gedanken experiment [4], which was firstly proposed by Wald to destroy an extremal Kerr-Newman black hole by throwing a massive particle. Moreover, by considering the gravitational collapse of a scalar field, the authors in [5] firstly found that the final fate of a continual collapse would be a naked singularity. See [6] for a review. It's worth noting that time evolution of the matter field leads to formation of a naked singularity in finite time. Obviously, it is an interesting question to ask whether there exists a region of arbitrarily large curvature that is time-independent and observed to distant observers. It is well known that because in higher dimensions the event horizon of a class of black holes is not necessarily the case of topologically spherical surface, black holes would be unstable under gravitational perturbations in the region of the Gregory-Laflamme instability [7] and could produce a naked singularity, which implies the violation of cosmic censorship [8][9][10]. Besides, by analysing the non-linear evolution of black holes [11,12], one found that a naked singularity should form in finite time. In four dimensions, it is more difficult to construct the naked singularity. By investigating the model of Einstein gravity coupling to a Maxwell field with a negative cosmological constant [13,14], one obtained a class of counterexamples to cosmic censorship with asymptotically anti-de Sitter (AdS) boundary for the first time. Furthermore, considering the above model with a vanishing Maxwell field, the authors in Ref. [15] constructed smooth stationary solutions with differential rotation to the boundary metric, which also provides a possible vacuum counterexample to weak cosmic censorship in AdS spacetime.
Recently, by introducing a charged scalar field to the Einstein-Maxwell solutions [16], the authors found that when the charge of the scalar field was sufficient large, the static solution of Einstein-Maxwell action would become unstable and a new stable solution could appear with nontrivial charged scalar field, which is analogous to the instability of a charged black hole to develop scalar condensation in the study of holographic superconductor [17][18][19]. Meanwhile, there exists a minimum value of charge required to remove the counterexamples and preserve cosmic censorship. It is surprising that the minimum value of charge appears to agree precisely with that proposed in the weak gravity conjecture [20], which states that any consistent quantum theory of gravity must have a stable state whose charge-to-mass ratio q/m is equal to or larger than that of an extremal black hole. In the case of Einstein-Maxwell model, the bound of q/m is equal to 1. To further test the weak gravity -cosmic censorship connection, the authors in [21] also studied the static solutions in Einstein theory with a dilaton field and the multi-charged scalar field case.
Considering that Maxwell theory is only the theory of linear electrodynamics, we would like to know whether or not there exists a class of counterexamples to cosmic censorship and the charge-to-mass ratio bound in four dimensional Einstein-nonlinear electrodynamics theory. The Born-Infeld electrodynamics [22][23][24] is a nonlinear generalization of the Maxwell's theory, which can remove the divergence of self-energy of a point-like charge in Maxwell electrodynamics. Moreover, the Lagrangian of Born-Infeld can arise from the low-energy effective theory describing electromagnetism [25,26]. Many works on the black holes solutions in Einstein-Born-Infeld theory have been studied in [27][28][29][30][31][32][33][34].
In the present paper, we show how to construct a class of counterexamples to cosmic censorship in four dimensional Einstein-Born-Infeld theory with asymptotically anti-de Sitter boundary. Furthermore, comparing with the counterpart of Einstein-Maxwell theory, we investigate the counterexample of Einstein-Born-Infeld theory for several values of Born-Infeld parameter b. Besides, introducing a charged massive scalar field into Einstein-Born-Infeld theory, we find that this class of counterexamples to cosmic censorship would be removed, and we also study the value of the lower bound on the scalar field charge required to preserve cosmic censorship in the case of Born-Infeld action.
The paper is organized as follows. In Sect. 2, we introduce the model of Einstein-Born-Infeld coupling to a complex, charged scalar field and the numerical DeTurck method. In Sect. 3, we explore the ansatz of metric and matter field, and analyze the boundary conditions. Numerical results of a class of counterexamples to cosmic censorship and static solutions with charged scalar condensation are shown in Sect. 4. The conclusion and discussion are given in the last section.

Set up
Let us begin with the action of the Born-Infeld field and a charged complex scalar field in the four-dimensional Einstein gravity spacetime with a negative cosmological constant, which is written as with the field strength of the U(1) gauge field F = F ab F ab , and L is the radius of asymptotic AdS spacetime. The constants m and q represent the mass and the charge of the complex scalar field , respectively. The constant b is the Born-Infeld parameter and the Born-Infeld field will reduce to the Maxwell case when b → 0. Where D a = ∇ a − iq A a is the gauge covariant derivative with respect to A a . Note that the values of m 2 must satisfy the Breitenlohner-Freedman (BF) bound m 2 ≥ −9/4 [35] for the (3+1)-dimensional spacetime.
The motion equations can derived from Eq. (2.1) with the energy-momentum tensor of matter field where F ≡ ∂L B I ∂ F , and it equals to -1 in Maxwell condition. If the complex, massive scalar field ψ vanishes, the solution of Einstein equations (2.2), which can describe the asymptotically spherically black hole with charge, is the well-known Born-Infeld AdS black hole. In terms of spherical coordinates, the metric of Born-Infeld AdS black hole is the following form where M and Q are the ADM mass and the electric charge of BI AdS black hole, respectively, and 2 F 1 is a hypergeometric function [37]. The gauge potential is When there exists a non-trivial configuration of the complex scalar field, it is obvious that we should solve the equations of motion (2.2) numerically instead of seeking the analytical solutions. We will use DeTurk method to solve these equations, which provides a good tool for solving Einstein equations in these papers. By adding a gauge fixing term to Einstein equation (2.2a), we could obtain a set of elliptic equations, which are known as Einstein-DeTurk equation is the Levi-Civita connection associated with a reference metricg, which should be choose to be as same boundary and horizon structure as g.

Ansatz and boundary conditions
Motived by the coordinate systems used in [13,16,21,36], we also consider Poincaré coordinates for pure AdS spacetime, When we use the new coordinate system, the metric (3.2) can become the following form: where x and y are polar-like coordinates. The Poincaré horizon is now located at y = 1, and the axis of rotation is located at x = 0. The conformal boundary is located at x = 1, and y = 0 denotes the intersection of the conformal boundary with the axis of symmetry. In order to construct static, axisymmetric solutions with a timelike Killing vector and an axisymmetric Killing vector, we adopt the following ansatz of metric 3), which is the reference metricg we take in DeTurk method.
Considering the above metric (3.4), an ansatz of matter fields should be described as below where the function U 6 and U 7 are the function of x and y. According to AdS/CFT duality, we take O 2 = (1− y 2 ) 2 U 7 to describe the scalar condensation. To simplify, we only choose m 2 = −2 in our paper. Next, we will discuss the boundary conditions. At conformal boundary, located at x = 1, the metric must reduce to pure AdS spacetime, so we must take and the boundary condition of scalar field is which can be obtained by expanding the equations of motion near x = 1 as a power series in (1 − x). For the gauge field at the boundary, we take where the constant a is the amplitude. This kind of profile for gauge field has been studied in [13,16,21,36]. The asymptotic behaviors of the equations of motion near x = 0 give the conditions at x = 0: At y = 0, i.e. the intersection of the conformal boundary with the axis of symmetry, we set the boundary conditons (3.11) and at Poincaré horizon y=1, Dirichlet boundary conditions are imposed on

Numerical results
In this section, we will numerically solve the coupled system of nonlinear partial differential equations (2.2) with the ansatzs (3.4) and (3.5). We employ finite element methods in the integration region 0 ≤ x ≤ 1 and 0 ≤ y ≤ 1 defined on non-uniform grids, allowing the grids to be more finer grid points near the boundary of y = 0 and y = 1. Our iterative process is the Newton-Raphson method. The relative error for the numerical solutions in this work is estimated to be below 10 −5 . In order to keep good agreement with the aforementioned error, the grid size has to be increased and typically a 120 × 200 to 120 × 300 grid was used. We would study the following two cases with scalar field = 0 and = 0, respectively. For convenience, in the following results, we will set L = 1.

With charged scalar field = 0
In this subsection, we will analyze the solutions with charged scalar field = 0 and present the evidence for counterexample to cosmic censorship. When one obtains a solution of Einstein equations, it is important to investigate whether or not the spacetime of solution is regular. In general, one of the most useful ways is to check for the finiteness of the Kretschmann scalar, which is also called Riemann tensor squared and written as where R αβγ δ is the Riemann curvature tensor. Therefore, Kretschmann scalar is a sum of squares of tensor components and a quadratic invariant. Numerical results are presented in Fig. 1. In the left panel we present the distribution of Kretschmann scalar K as a function of x and y coordinate with b = 0.25, n = 8 and a = 6. It is obvious that the spacetime is not flat, and the value of the Kretschmann scalar in the purple area is about to 24/L 4 , which means it recovers the result of Ad S 4 spacetime. In addition, the maximum of the Kretschmann scalar appears   [16]. When the amplitude reaches to a maximum value, there will appear a singularity, which means producing arbitrarily large curvature in spacetime. This maximum depends on both n and b. We show our results in the following table.
From Table 1, we could see that the maximum amplitude a max increases with the increasing of n. Meanwhile, it decreases with the increasing of the Born-Infeld parameter b. When Maxwell field is extended to Born-Infeld model, the static solutions could still become singular at a finite value of a.

With charged scalar field = 0
In the last subsection, we obtain a family of Einstein-Born-Infeld solutions with charged scalar field = 0. We would like to know whether or not there exist the charged static solutions with no-vanished scalar field, which is analogous to the case of holographic superconductor studied in [17,19,38,39]. When fixed the parametes b and a, there is a minimal value of q, below which the scalar field vanishes and the solution is the simple planar. Above this value, the charged static solution becomes unstable to develop scalar hair. To obtain the minimal value of the charge q, we could use the method in [16,40] to treat the equation of motion of perturbatively. Therefore, before numerically solving the full dynamic equations of motion including scalar field in Eqs.
(2.2), we solve the time-independent scalar field equations to search for zero mode at the fixed background constructed in Sect. 4.1, which one could recognize as an eigenvalue problem with the eigenfunction and eigenvalue q 2 . We will use q min to denote the smallest eigenvalue of zero-mode, which is the minimal charge and is also a function of the amplitude a.
The results are shown in Fig. 2. In left panel, with the fixed n = 8, we plot minimal charge q min as a function of the amplitude a for the corresponding values of b = 1, 0.5, 0.25 and 0, represented by cyan, blue, red and black lines, respectively. The vertical dashed gridlines indicate the maximum amplitude a max = 8.06 (black), 6.56 (red), 5.94 (blue) and 5.20 (cyan), respectively, which corresponds to the values of n = 8 shown in Table 1. For each curve, the minimal charge q min decreases with the increasing of amplitude a. We see that for the fixed amplitude a, the critical charge q min increases with the increasing of Born-Infeld parameter b. Furthermore, when the amplitude a decreases, the minimal charge q min tends to the same value for several values of b.
In the right of Fig. 2, with the fixed b = 0.25, we plot minimal charge q min as a function of the amplitude a for the corresponding values of n = 10, 8, 6 and 4, represented by cyan, blue, red and black lines, respectively. For each curve, the minimal charge q min decreases with the increasing of amplitude a, which is similar to the behavior in the left panel. For the fixed amplitude a, the minimal charge q min increases with the increasing of n.
Next, we could solve the coupled equations of motion (2.2) to find the solutions for the above critical charge q min . We show the profiles of the expectation value for the operator dual to as a function of the boundary radial coordinate r in Fig. 3. In the left panel, with the fixed n = 8 and q = 2.7, the curves from top to bottom correspond to b = 0, 0.5, 1, represented by black, red and blue lines, respectively. For each curve, the scalar condensation O 2 decreases with the increasing amplitude of r . When the amplitude r increases, the scalar condensation tends to zero value for several values of b. Moreover, at the origin r = 0, there exists the largest value of the scalar condensation. Furthermore, for the fixed amplitude r , the scalar condensation O 2 decreases with the increasing value of Born-Infeld parameter b.
In the right of Fig. 3, with the fixed b = 0.25 and q = 2.22, we plot the scalar condensation O 2 as a function of r for several values of a. The curves from top to bottom correspond to a = 10, 9, 8, 7 and 5.8, represented by black, red, blue, cyan and purple lines, respectively. For each curve, the scalar condensation decreases with the increasing of r , which is similar to the behavior in the left panel. For the fixed amplitude b, the scalar condensation O 2 increases with the increasing value of the parameter a.
In order to further understand how increasing the amplitude a affects the maximum of the condensation with the different values of b, we show in the left panel of Fig. 4 the condensate O 2 at r = 0 as a function of the amplitude a with the fixed parameters n = 8, q = 2.7. From top to bottom, the curves correspond to the values of b = 0, b = 0.5, and b = 1, respectively. In addition, the vertical dashed gridlines indicate the critical amplitude a c = 4.74 (black), a c = 4.82 (red) and a c = 4.98 (blue), respectively. As one can see that there exists a critical amplitude a c above which the condensate appears, then rises as the system is imposed with much larger amplitudes a. This behaviour of After obtaining the numerical solution with charged scalar condensate, we want to know what happens to the solutions with = 0 as the decreasing of charge q when a > a max . We present in the right panel of Fig. 4 the typical maximal value of Kretschmann scalar K as a function of the charge q with b = 0.5, n = 8 and a = 6 > a max , and the vertical black dashed gridlines indicate the q min = 2.276. From this plot, we can see that the full nonlinear solutions with the scalar condensate are in the range of q > q min . It is obvious that there exists the growth of the maximum of Kretschmann scalar with the decreasing of q, which indicates the formation of a curvature singularity outside Poincaré horizon at some minimal values q = q min , and is similar to the case of Maxwell case.
To study the properties of the bound on charge-to-mass ratio, we exhibit the phase diagram of the minimal charge q min versus the amplitude a with several values of Born-Infeld parameter in Fig. 5, where we add the curves of minimal charge with a > a max to previous plot for a < a max in the left panel of Fig. 2. The Born-Infeld parameters b = 0.25, 0.5, 1 are represented by the red, blue and orange lines, respectively. The black line represents the curve in the Maxwell model which has been discussed in [16]. We can see that there exists an asymptotical value q w of q min as a increases. The asymptotical value q w increases with the increasing of Born-Infeld parameter b. More details are shown in the inset, and the horizonal black, red, blue and orange dashed gridlines indicate the values of q w = 2, 2.23, 2.43, 2.708, respectively. If we consider scalar field with q < q w , there still exists a singularity outside the Poincaré horizon as a increases to a finite value. But for keeping q > q w , there would not exist a singularity when one increases the amplitude a to an arbitrarily large value.
It is noteworthy that for the curve of Maxwell model in the solutions with = 0, the minimal charge q min decreases firstly with the decreasing of the amplitude a, and then it reaches a minimal point. Further increasing a to a max , the value of q min continues to decrease, and a second branch with lower q min is obtained. When the curve with lower q min moves toward a max , the numerical error begins to increase and a finer mesh is required to calculate. It is difficult to handle numerical calculation near the value of a max , so we could not get the curve connecting the solutions of = 0 and = 0. However, for larger value of b, it is relatively easy to obtain the connected curve near the value of a max .

Conclusions
In this paper, we have presented a class of counterexamples to cosmic censorship in four-dimensional Einstein-Born-Infeld theory. The maximum of Kretschmann scalar goes to infinite when the amplitude a increases to a max , which means there appears a singularity outside Poincaré horizon. If we consider the charged scalar field in this model, we found there exists a minimal value of charge q min , below which the static solution with zero scalar condensation is stable. However, when the charge of the scalar field was above q min , the static solution of Born-Infeld action will become unsta- ble and a new stable solution with nontrivial charged scalar field could appear, which is analogous to the instability of a charged static solution to develop scalar condensation in the case of Maxwell theory [16]. Meanwhile, there exists the asymptotical value of charge q w required to remove the counterexamples and preserve cosmic censorship. Comparing with the Einstein-Maxwell model, the bound of chargeto-mass ratio q/m is larger than 1 and increases with the Born-Infeld parameter b.
In the case of Maxwell theory [16], it is surprising to find that the bound of the charge to mass ratio which is necessary to preserve weak cosmic censorship is precisely equal to the weak gravity bound. For the nonlinear Born-Infeld theory, we also obtain the the bound of the charge to mass ratio necessary to preserve weak cosmic censorship, but, because now the weak gravity conjecture in Born-Infeld electrodynamics is still not clear, one could not determine whether bound of charge-to-mass ratio is equal to weak gravity conjecture. However, it is sure that a class of static solution with charged scalar condensate could prevent the violation of the weak cosmic censorship, and the bound of charge-to-mass ratio could also be obtained in the case of Born-Infeld electrodynamic.
It will be interesting to investigate the several further researches. First, since we have studied the static solution with a free scalar condensate in Born-Infeld electrodynamic, we would like to investigate how self-interactions of the scalar field to prevent the violation of the weak cosmic censorship and the bound of charge-to-mass ratio. The second extension of our study is to construct generalized multiscalar hair configurations, such as two coexisting states of the charge scalar field presented in Born-Infeld electrodynamic. Finally, we are planning to study the model of the Einstein-Born-Infeld-vector model and construct the static solution with charged vector hair necessary to preserve weak cosmic censorship in future work.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The datasets generated during the current study are available from the corresponding author on reasonable request.] 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 .