Magnetic-induced spontaneous scalarization in dynamical Chern–Simons gravity

In the framework of the dynamical Chern–Simons gravity, we study the scalar field perturbations of the Reissner–Nordström–Melvin spacetime, which describes a charged black hole permeated by a uniform magnetic field. In the presence of the magnetic field, the scalar field acquires an effective mass whose square takes negative value in the half domain of the angular direction. This inevitably introduces the tachyonic instability and associated spontaneous scalarization as long as the coupling constant between the scalar field and the Chern–Simons invariant exceeds a threshold value. We study the object pictures of the time evolutions of the scalar field perturbations at the linear level, and find that the presence of the magnetic field will dramatically change the waveforms and associated ringdown modes. Nonlinear evolutions for the unstable perturbations are also performed in the decoupling limit, which demonstrate the scalar cloud as the final fate. Influences of the coupling constant and the black hole charge on the wave dynamics are also studied.


I. INTRODUCTION
In recent years there is a lot of activity trying to understand the plethora of astrophysical observations.First of all, the LIGO and Virgo collaborations reported the detection of Gravitational Waves (GW) from a Binary Black Hole (BH) Merger [1][2][3][4][5].More recently the First M87 Event Horizon Telescope reviled the shadow of a supermassive BH [6].Also, the recent results of the Event Horizon Telescope are very interesting detecting a magnetic field structure near the event horizon of a BH [7].General Relativity (GR) describes accurately the formation and the properties of BHs.A BH is the end state of the collapse of a very massive cloud [8]- [13].The generated structure is characterized by only three parameters, the mass, charge and angular momentum, obeying the powerful no-hair theorems [14][15][16].
In the same time many astrophysical observations showed that many strong radio sources take the form of two emitting regions situated on opposite sites of a galaxy.A theory to account for this is that the magnetic fields and high energy particles responsible for the synchrotron radiation were blown out of the galactic halo in a giant explosion.It was speculated that such explosions could be generated from gravitational collapse.A simple model was presented [17], which was composed by a configuration that contains only electromagnetic field in a form of a collection of parallel magnetic lines hold together by its own gravitational attraction.Again GR can describe such a theory as a rigorous static cylindrically-symmetric solution of the combined sourceless Einstein-Maxwell equations [18]- [20] [For another possibility, see [21,22] and references therein].
The recent observational results on dark matter and on dark energy require a generalization of GR in an attempt to have a viable cosmological theory of Gravity on short and large distances [23][24][25][26][27].We expect that these modified gravity theories (MOGs) will provide important information on the properties and structure of the compact objects which are described by these theories and they are consistent with the observational signatures, which they introduce.In these theories matter sources parameterize by a scalar field is directly coupled to second order algebraic curvature invariants.One of the most well studied MOGs theories are the scalar-tensor theories in which matter is interacting with black holes like the scalar-tensor Horndeski theories [28].Exploring the strong field regime of gravity with the aim to detect gravitational waves and black hole shadows the effects of higher-order curvature terms become significant.The most well known high dimensional theories are the scalar-Einstein-Gauss-Bonnet theory (sEGB) [29]- [51] and the dynamical Chern-Simons gravity theory (dCSG) [52]- [84].
In dSCG, a dynamical scalar field is introduced to couple to the gravitational Chern-Simons invariant non-minimally [70].The motivation for studying dSCG is that it is a gravity theory that contain higher powers of curvature which are consistent Lanczos-Lovelock gravity theories in higher than four dimensions.The Chern-Simons black holes are special solutions of these theories resulting in second order field equations for the metric with well defined AdS asymptotic solutions.
The action for CS modified gravity is defined by the sum of the Einstein-Hilbert action and a new parity violating, four-dimensional correction.Interest in the model spiked when it was found that string theory unavoidably requires such a correction to remain mathematically consistent [72,73].
In general, such a correction arises in the presence of Ramond-Ramond scalars because of the presence of duality symmetries and these parity violating scalar fields can dress a slowly rotating black hole with scalar hair.
One of the first work of a hairy slowly rotating black hole, which resulted because the scalar field was coupled to a Lorentz CS term, was presented in [74,75].This result suggested that non-minimal gravitational couplings may produce interesting new effects in black hole backgrounds.These works were further extended in [76] where it was found that the scalar hair are characterized by the mass, angular momentum and gauge charges of the background rotating black hole.Allowing the scalar field to have a dynamical behavior a solution describing a rotating black hole in the small coupling slow rotation limit was found in [59].Static and rotating black string solutions in dynamical CS modified gravity were also studied in [77,78].
The coupling of a scalar field to high curvature terms allows us to evade the no-hair theorems and obtained more general hairy black hole configurations.In particular, for certain classes of the coupling function it was shown in the case of GB high curvature terms, that we have spontaneous scalarization of black holes [32]- [36].It was found in regions of strong curvature, that below a n critical mass the Schwarzschild black hole develops instabilities and then when the scalar field backreacts to the metric, new branches of scalarized black holes are generated at certain masses as solutions in the theory [32,33,37].Actually, this novel phenomenon has been observed long time ago in neutron stars but there the instability is triggered by the surrounding matter instead of the curvature [85].The mechanism of spontaneous scalarization can also be applied to the Chern-Simons theory.More precisely, on the Kerr BH background, the massless scalar field perturbations will acquire an effective mass due to its coupling to the CS invariant.The square of the effective mass becomes negative in half domain of angular directions resulting in the tachyonic instability, and thus leading to the spontaneous scalarization.However, most of the studies on spontaneous scalarization have considered that the background BHs live in vacuum.In the realistic universe, BHs are usually not isolated but always live in complicated astrophysical environments, such as magnetic fields, accretion disks, dark matter halos, cosmological expansion, etc.These environments may have considerable influences on dynamics of BHs under perturbations, dramatically change the waveforms and associated ringdown modes [86,87].Then one can ask how much of the impact does the realistic astrophysical environments have on the phenomenon of spontaneous scalarization?
In this work, we will focus on the influences of one factor of these environments, the magnetic field.It is believed to be ubiquitous in nature and pervade our universe on various scales with different amplitudes, ranging from order of ∼ 10 −4 Gauss at the center of our galaxy [88] to order of ∼ 10 16 Gauss at surface of some magnetars [89].Also, existences of BHs permeated in magnetic backgrounds are also supported by astrophysical observations of recent years.For example, the presence of the magnetar SGR J1745-2900 orbiting the supermassive BH Sagittarius A * [89][90][91][92], and the presence of strong magnetic field in vicinity of the event horizon of M 87 * [7].Moreover, it is believed that strong magnetic fields around BH play a central role in some of the most energetic events of our universe, such as emission of relativistic jets via, for example, the Blandford-Zanjek process [93].
In general, full comprehension of the interactions between BH and the surrounding magnetic field is a rather involved problem.As a simple model, we consider stationary magnetized BH solutions in which a qualitative picture can be drawn.Actually, a class of exact solutions of the Einstein-Maxwell equations has been known, describing BHs immersed in a uniform magnetic field aligned along the symmetry axis, the Kerr-Newmann-Melvin (KNM) BHs [94][95][96][97].For a recent review of this subject, we refer readers to [98,99].Different from the well-known Kerr-Newmann BH in GR, the KNM BH spactime is not asymptotically flat but resembles the magnetic Melvin universe [17].The thermodynamical properties of these BHs have been studied extensively in recent years [100][101][102][103][104].
In the KNM model, effects of the magnetic field on the spontaneous scalarization is studied preliminary in sEGB theory in [105,106].With this by analyzing the behavior of effective mass square of the scalar field perturbation near the horizon directly, it was found that the presence of the magnetic field pushes up the threshold dimensionless spin of the BH for spin-induced tachyonic instability.However, a complete object picture of the times evolution of the scalar field perturbations to show the influences of the magnetic field on the dynamics is still lack.Interestingly, it is found in [107] that even in the absence of compact objects (like BHs), a real massive scalar field can still condense in the Melvin magnetic universe when coupled non-minimally either to the magnetic field or to the curvature.Another type of mechanism generating BH hairs -the superra-diant instability in these BHs has also been studied in [87,108].By considering massless neutral scalar field perturbations, they found that the presence of magnetic fields enacts both requirements simultaneously for the occurrence of the superradiant instability and formation of stationary scalar clouds, namely an ergoregion to trigger superradiance and a confinement mechanism to trap the scalar field.In [109] a magnetized compact object, like a neutron star, was studied in scalar-tensor theories.
Inspired by these works, we will study the effects of magnetic fields on spontaneous scalarization in dCSG.We take Reissner-Nordstrom-Melvin (RNM) BH as the background.This BH is characterized by three parameters, the BH mass M , the BH charge Q and the magnetic field B. It reduces to the RN BH in the limit B → 0, while reduces to the Schwarzschild-Melvin solution in the limit Q → 0 [95].It should be noted that in the absence of magnetic fields, the CS invariant (and thus the effective mass square that the scalar field acquires) vanishes identically.So in this case, there is no tachyonic instability and spontaneous scalarization.However, when the magnetic field is present, situation changes in two aspects.On the one hand, the scalar field acquires an effective mass whose square may become negative thus would trigger tachyonic instability.On the other hand, the presence of magnetic fields will change the asymptotical structure of the spactime which then may influence the waveforms of the scalar field perturbations.Therefore the presence of the magnetic field introduces high symmetries making the system more involved compared to a system of astrophysical black holes which are nearly asymptotically flat.Then natural questions arise: How do the presence of the magnetic field affect the stability of the BH?And how do the presence of the magnetic field influence the wave dynamics of scalar perturbations?These are main questions we would like to answer in this work.
The paper is organized as follows.In Sec.II, we will give a brief introduction of dCSG theory and the RNM BH and establish the scalar field perturbation equation.In Sec.III, we describe our numerical method to solve the scalar field perturbation equation in detail.In Sec.IV, we report our numerical results.The final section is devoted to summary and discussions.

II. THE MODEL
In this work, we utilize the units c = G = 4π 0 = 1, where c, G, 0 are the the speed of light in vacuum, the Newton gravitational constant and the vacuum permittivity respectively.
We consider dCSG with an additional electromagnetic field.The general action is [79] where the real scalar field Φ is non-minimally coupled to the CS invariant through the coupling constant α and the coupling function h(Φ).V (Φ) is the scalar self-interaction potential.The explicit definition of the CS invariant is By varying the action, one can derive the equations of motion In this work, we consider the scalar field to be massless without self-interaction so that V (Φ) = 0, and take the form of the coupling function as where β is a constant.With these choices, the theory admits GR electro-vacuum solutions with vanishing scalar hair Φ = 0, among which to our interest in this work is the Reissner-Nordström-Melvin (RNM) solution [95,97] where which describes a RN black hole with charge Q permeated by a uniform magnetic field with strength B aligned along the symmetry axis.Note that this BH solution has cylindrical symmetry [22], and is not asymptotically flat but resembles the magnetic Melvin Universe [17].The event horizon is It reduces to the RN black hole solution in the limit B → 0, while reduces to the Schwarzschild-Melvin solution in the limit Q → 0 [95].On this background, we would like to study the wave dynamics of scalar field perturbations and study possible instabilities.
Then the scalar field perturbation equation ( 2) becomes where the effective mass square m 2 eff = −α * RR is position-dependent and depends on parameters (M, Q, B).The full explicit expression of the CS invariant * RR, valued in the background, is rather involved and uninspiring and will not be displayed here.One can get a preliminary feeling of the influence of the magnetic field on the CS invariant from its small-B expansion It should be noted that throughout this work we will perform the computations with the full form of the CS invariant rather than its small-B expansion.Several comments on the scalar field perturbation equation ( 9) now are in order: • For vanishing magnetic field B = 0 (vanishing BH charge Q = 0), the background reduces to RN (Schwarzschild-Melvin) BH, and the CS invariant and the effective mass square vanish identically.In this case, the above equation describes wave propagation of a massless neutral scalar field in the RN (Schwarzschild-Melvin) background which has already been studied thoroughly and no instability is observed [87,110,111].
• The CS invariant is odd under the transformation So in the scalar field perturbation equation ( 9), the sign of α can always be absorbed into the CS invariant by redefining the θ-coordinate.Taking into account this symmetry, we will only consider α > 0 in the following.Moreover, it also means that, for non-zero magnetic field B/BH charge Q, the effective mass square m 2 eff will always take negative value in half interval of θ domain, implying the inevitable tachyonic instability as long as the coupling constant α is large enough.
In Fig. 1, profiles of the CS invariant are shown for some typical values of parameters.From the figure, one can see that the CS invariant * RR takes a finite non-zero value near the horizon while approaches to zero at infinity.This can be understood because the CS term is a high curvature term, so the strong gravitational field is expected to occur only near the horizon of the BH while it is negligible far away from the BH.Moreover, increasing the magnetic field B/BH charge Q will make the CS invariant more positive (thus the effective mass square m 2 eff more negative for α > 0) in the half interval of θ domain near the horizon, hence triggering more violent instability.
In the following sections, we will study carefully the time evolution of the massless scalar field perturbations and obtain object pictures of the influences of the magnetic field B, the coupling constant α and also the BH charge Q on wave dynamics.

III. NUMERICAL METHOD
We will apply the numerical method discussed in [112][113][114][115][116] to solve the scalar field perturbation equation (9).With the RNM metric Eq. ( 6), the scalar field perturbation equation ( 9) becomes After introducing the tortoise coordinate x defined as dx ≡ dr/f to map the radial domain r ∈ (r + , +∞) to x ∈ (−∞, +∞), the above equation becomes In the case of scalar field perturbations in the Kerr BH [113], it has already been noted that the azimuthal coordinate ϕ will introduce unphysical pathologies near the horizon, which can be removed by adopting the Kerr azimuthal coordinate [113].In the current case, a similar problem will appear and can be solved by introducing a Kerr-like azimuthal coordinate φ defined as Then the equation becomes Due to the axial symmetry of the RNM spacetime, the scalar field perturbation can be decomposed as with m being the azimuthal number.Then, by introducing an auxiliary variable Π ≡ ∂ t Ψ, the perturbation equation can be cast into a form of two coupled first-order partial differential equations which are suitable for the method presented in [112].Precisely, we apply the fourth-order Runnge-Kutta integrator to perform the time evolution, while the finite difference scheme for the spatial derivatives.
To solve the perturbation equation, we also need to impose physical boundary conditions.At the horizon, ingoing wave condition is implemented following [117].At the outer boundary, the magnetic field will behave like an infinite "wall" at x ∼ 1/B to confine the perturbation [87], so it is sensible to impose a Dirichlet boundary condition Ψ = 0 there.Also, at the poles θ = 0 and π, we should impose physical boundary conditions Ψ| θ=0,π = 0 for m = 0 while ∂ θ Ψ| θ=0,π = 0 for m = 0 [115].

IV. NUMERICAL RESULTS
We consider the scalar field perturbation initially to be a Gaussian wave-packet localized outside the horizon at x = x c which has time symmetry [83,118], where Y m is the θ-dependent part of the spherical harmonic function and σ is the width of the Gaussian wave-packet.
As in the Kerr spacetime, the RNM spacetime we consider here is not spherically symmetric except when B = 0, so mode-mixing phenomenon [119][120][121] will also occur during evolution: a pure initial -multipole will excite other -multipoles with the same m as it evolves with time.
Taking into account this phenomenon and for simplicity, in the following we will only consider axisymmetric perturbations with = m = 0. Also, we set M = 1 so that all quantities are measured in units of M .Without loss of generality, observers are assumed to locate at x = 6M and θ = π 5 .

Object pictures: Linear and nonlinear level
In Fig. • There exists a critical value of the coupling constant α = α c 7.9, above which tachyonic instability occurs.Moreover, with the increase of α, the instability occurs earlier and becomes more violent.This is expected as the effective mass square m 2 eff ∝ α.
• For α ≤ α c , the waveforms show two types of ringdown modes with different frequencies at different stages.This interesting phenomenon has already been observed in [87], where minimally coupled neutral scalar field perturbations of the Schwarzschild-Melvin BH in framework of GR is studied.Following the arguments there, in our case one can expect that at early times, the ringdown modes are being in fact similar to the quasinormal modes (QNMs) of RN BH under massive scalar field perturbations with the additional effective mass µ eff = mB [122][123][124], which are excited near the horizon.After a time of order t ∼ 1 B , the wave is reflected back by the effective "wall" induced by the magnetic field and the so-called "Melvin-like" modes with smaller amplitudes are excited [87].Both types of modes decay with time and are expected to vanish finally.See also [86] for more discussions on similar interesting phenomenon in the context of "dirty" BHs.
• When α = α c , the waveform between the two types of modes approaches a constant indicating criticality.
• For α > α c , unstable tachyonic modes are excited and the Melvin-like modes are greatly suppressed.When α is large enough, the Melvin-like modes are completely over-dominated by the unstable tachyonic modes and thus have not been observed.
To see the final fate of the instability, we have also studied the time evolutions of the scalar field perturbation in the nonlinear level but in the so-called "decoupling limit" [82,125].That is, we take into account the non-linearities in the scalar field coupling while neglecting its back-reaction on the spacetime.This limit has been shown to be able to capture the basic features of the full nonlinear dynamics qualitatively.In this limit, one should consider the scalar field equation with the full coupling Time evolutions of the scalar field perturbation in this limit are shown in Fig. 2 for α = 0.When α = 0, m 2 eff = 0 and there are no nonlinear effects.From the figure, one can see that when α = 5.0 < α c , there is little difference between the linear and nonlinear levels.This is expected as nonlinear effects can be neglected for stable perturbations.However, for α > α c , when taking into account the nonlinear effect, the instability is quenched and the unstable tachyonic modes are suppressed to approach a non-zero constant while the Melvin-like modes thus revive.With the slowly decaying of the Melvin-like modes in time, it is expected that the scalar field will approach a non-zero constant finally to form a scalar cloud around the BH.Moreover, it is also observed that with the increase of α, the scalar field approaches a larger constant finally which can be understood physically.

Effects of B and Q
In Figs. 3 and 4, we plot time evolutions of the scalar field perturbation for various values of B and Q to see their effects on the dynamics of wave propagation.Meanwhile, we also show their From the left panel of Fig. 3, one can see that there exists a critical value of B = B c 0.155, above which tachyonic instability is triggered.For B ≤ B c , increasing B will make the effective "wall" closer to the horizon and thus the Melvin-like modes excited earlier.For B > B c , once again we can see that the Melvin-like modes is suppressed by the unstable tachyonic modes.Increasing B will make the instability to occur earlier and more violent, which can be understood physically from Fig. 1 as larger B will make the effective mass square more negative near the horizon.The "Melvin-like" modes with smaller amplitude are excited due to the reflections of the wave from the effective "wall" induced by the magnetic field [87].Both types of modes decay with time (although very slowly) and are expected to vanish finally.However, when α > α c , unstable tachyonic modes are excited and the Melvin-like modes are suppressed greatly.When taking into account the non-linearities of the scalar field coupling in decoupling limit, the instability is quenched and the unstable tachyonic mode approaches to a non-zero constant to form a scalar cloud.It is expected that when back-reaction of the scalar field to the spacetime geometry is considered, the scalar cloud will become the scalar hair of the BH.
Let us give some comments on the astrophysical relevance of the parametric value of the magnetic field we considered.In this work, typical value of the magnetic field considered is BM ∼ 0.1.
Retaining physical units, we have the following relation where M is the solar mass.So BM ∼ 0.1 corresponds to a magnetic field B ∼ 2.36 × (for example the Sagittarius A * ), B ∼ 10 12 Gauss.In general, for BHs with mass M > 10 2 M , the magnetic field considered will be smaller than the ever-measured strongest magnetic field B ∼ 10 16 Gauss [89].Moreover, as have been observed, even for smaller magnetic field, the same tachyonic instability and spontaneous scalarization can still be triggered as long as the coupling constant α is large enough.
Of course, a more precise description of the formation of the hair requires nonlinear evolutions of the full system including the spacetime, which will be rather involved technically.It will also be interesting to construct the final spontaneously scalarized black hole solutions, which is still an open question in dCSG theory, whether in our considered model in this work or in the original model without magnetic field.Moreover, it will also be interesting to extend current studies to the rotating case by considering the background BH to be the Kerr-Melvin BH, or other models of spontaneous scalarizations, like the Einstein-Maxwell-scalar model [126], to see the effects of the environmental magnetic field.We would like to leave these questions for further investigations.

FIG. 1 :
FIG. 1: (color online) The CS invariant * RR as a function of r and θ.In the left panel, Q = 0.5 is fixed while B is varied; In the right panel, B = 0.1 is fixed while Q is varied.

2 ,
time evolutions of the scalar field perturbations are shown for B = 0.1, Q = 0.5 and various values of the coupling constant α in linear level.For other values of parameters, we have similar phenomena.Several properties are observed:

FIG. 2 :
FIG. 2: (color online) Time evolutions of the scalar field perturbation for B = 0.1, Q = 0.5 and various values of α both in the linear and nonlinear levels (with the nonlinear parameter β = 10.0).From left to right and top to bottom, α increases from 0 to 20.0.Initial data of the perturbation is Ψ(t = 0, x) ∼ e − (x−xc ) 2 2σ 2 with

FIG. 3 :FIG. 4 :
FIG. 3: (color online) Left: time evolutions of the scalar field perturbation in the linear level for α = 5.0, Q = 0.5 and various values of B. Initial data of the perturbation is Ψ(t = 0, x) ∼ e − ) 2 2σ 2 with x c = 6.0 and

M
Gauss whose explicit value is inversely proportional to the mass of the BH; For stellarmass BHs with M ∼ 10M , B ∼ 10 17 Gauss, while for supermassive BHs with M ∼ 10 6 M