Researcher profile

Carlos A. R. Herdeiro

Carlos A. R. Herdeiro contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2026arXiv

There and back again: Outspiraling motion in non-Kerr compact objects

In Keplerian dynamics, a test body orbiting a point particle in circular motion has a monotonically increasing frequency, with decreasing radius. If a dissipative channel is introduced, such as gravitational wave (GW) emission, under the quadrupole approximation, the corresponding GW strain has an ever-increasing frequency with time. A similar statement holds for equatorial motion of a test particle on the Kerr manifold, except such inspiral is cut off at the ISCO, wherein stable circular orbits cease to exist and a plunge is expected. We analyze circular timelike orbits in generic spinning spacetimes and study the conditions in which exotic motion can occur, due to the presence of non-Kerr features. In particular, we derive conditions under which an inspiral towards a compact object is naturally followed by an outspiral motion, and give concrete examples as well as the corresponding GW phenomenology. This analysis serves both as a theoretical exploration of non-Kerrness as well as an example of a concrete smoking gun of exotic spacetimes.

preprint2022arXiv

Black holes: on the universality of the Kerr hypothesis

To what extent are all astrophysical, dark, compact objects both black holes (BHs) and described by the Kerr geometry? We embark on the exercise of defying the universality of this remarkable idea, often called the "Kerr hypothesis". After establishing its rationale and timeliness, we define a minimal set of reasonability criteria for alternative models of dark compact objects. Then, as proof of principle, we discuss concrete, dynamically robust non-Kerr BHs and horizonless imitators, that 1) pass the basic theoretical, and in particular dynamical, tests, 2) match (some of the) state of the art astrophysical observables and 3) only emerge at some (macroscopic) scales. These examples illustrate how the universality (at all macroscopic scales) of the Kerr hypothesis can be challenged.

preprint2022arXiv

Constraining spontaneous black hole scalarization in scalar-tensor-Gauss-Bonnet theories with current gravitational-wave data

We examine the constraining power of current gravitational-wave data on scalar-tensor-Gauss-Bonnet theories that allow for the spontaneous scalarization of black holes. In the fiducial model that we consider, a slowly rotating black hole must scalarize if its size is comparable to the new length scale $λ$ that the theory introduces, although rapidly rotating black holes of any mass are effectively indistinguishable from their counterparts in general relativity. With this in mind, we use the gravitational-wave event GW190814$\,\unicode{x2014}\,$whose primary black hole has a spin that is bounded to be small, and whose signal shows no evidence of a scalarized primary$\,\unicode{x2014}\,$to rule out a narrow region of the parameter space. In particular, we find that values of ${λ\in [56, 96]~M_\odot}$ are strongly disfavored with a Bayes factor of $0.1$ or less. We also include a second event, GW151226, in our analysis to illustrate what information can be extracted when the spins of both components are poorly measured.

preprint2022arXiv

Deconstructing scaling virial identities in General Relativity: spherical symmetry and beyond

Derrick-type virial identities, obtained via dilatation (scaling) arguments, have a variety of applications in field theories. We deconstruct such virial identities in relativistic gravity showing how they can be recast as self-evident integrals of appropriate combinations of the equations of motion. In spherical symmetry, the appropriate combination and gauge choice guarantee the geometric part can be integrated out to yield a master form of the virial identity as a non-trivial energy-momentum balance condition, valid for both asymptotically flat black holes and self-gravitating solitons, for any matter model. Specifying the matter model we recover previous results obtained via the scaling procedure. We then discuss the more general case of stationary, axi-symmetric, asymptotically flat black hole or solitonic solutions in General Relativity, for which a master form for their virial identity is proposed, in a specific gauge but regardless of the matter content. In the flat spacetime limit, the master virial identity for both the spherical and axial cases reduces to a balance condition for the principal pressures, discussed by Deser.

preprint2022arXiv

Generating gravitational waveform libraries of exotic compact binaries with deep learning

Current gravitational wave (GW) detections rely on the existence of libraries of theoretical waveforms. Consequently, finding new physics with GWs requires libraries of non-standard models, which are computationally demanding. We discuss how deep learning frameworks can be used to generate new waveforms "learned" from a simulation dataset obtained, say, from numerical relativity simulations. Concretely, we use the WaveGAN architecture of a generative adversarial network (GAN). As a proof of concept we provide this neural network (NN) with a sample of ($>500$) waveforms from the collisions of exotic compact objects (Proca stars), obtained from numerical relativity simulations. Dividing the sample into a training and a validation set, we show that after a sufficiently large number of training epochs the NN can produce from 12\% to 25\% of the synthetic waveforms with an overlapping match of at least 95\% with the ones from the validation set. We also demonstrate that a NN can be used to predict the overlapping match score, with 90\% of accuracy, of new synthetic samples. These are encouraging results for using GANs for data augmentation and interpolation in the context of GWs, to cover the full parameter space of, say, exotic compact binaries, without the need of intensive numerical relativity simulations.

preprint2022arXiv

Null and timelike circular orbits from equivalent 2D metrics

The motion of particles on spherical $1 + 3$ dimensional spacetimes can, under some assumptions, be described by the curves on a 2-dimensional manifold, the optical and Jacobi manifolds for null and timelike curves, respectively. In this paper we resort to auxiliary 2-dimensional metrics to study circular geodesics of generic static, spherically symmetric, and asymptotically flat $1 + 3$ dimensional spacetimes, whose functions are at least $C^2$ smooth. This is done by studying the Gaussian curvature of the bidimensional equivalent manifold as well as the geodesic curvature of circular paths on these. This study considers both null and timelike circular geodesics. The study of null geodesics through the optical manifold retrieves the known result of the number of light rings (LRs) on the spacetime outside a black hole and on spacetimes with horizonless compact objects. With an equivalent procedure we can formulate a similar theorem on the number of marginally stable timelike circular orbits (TCOs) of a given spacetime satisfying the previously mentioned assumption

preprint2022arXiv

Spin-induced scalarization and magnetic fields

In the presence of certain non-minimal couplings between a scalar field and the Gauss-Bonnet curvature invariant, Kerr black holes can scalarize, as long as they are spinning fast enough. This provides a distinctive violation of the Kerr hypothesis, occurring only for some high spin range. In this paper we assess if strong magnetic fields, that may exist in the vicinity of astrophysical black holes, could facilitate this distinctive effect, by bringing down the spin threshold for scalarization. This inquiry is motivated by the fact that self-gravitating magnetic fields, by themselves, can also promote "spin-induced" scalarization. Nonetheless, we show that in the \textit{vicinity of the horizon} the effect of the magnetic field $B$ on a black hole of mass $M$, up to $BM\lesssim 1$, works \textit{against} spin-induced scalarization, requiring a larger dimensionless spin $j$ from the black hole. A geometric interpretation for this result is suggested, in terms of the effects of rotation $vs.$ magnetic fields on the horizon geometry.

preprint2021arXiv

Black holes, stationary clouds and magnetic fields

As the electron in the hydrogen atom, a bosonic field can bind itself to a black hole occupying a discrete infinite set of states. When (i) the spacetime is prone to superradiance and (ii) a confinement mechanism is present, some of such states are infinitely long-lived. These equilibrium configurations, known as stationary clouds, are states "synchronized" with a rotating black hole's event horizon. For most, if not all, stationary clouds studied in the literature so far, the requirements (i)-(ii) are independent of each other. However, this is not always the case. This paper shows that massless neutral scalar fields can form stationary clouds around a Reissner-Nordström black hole when both are subject to a uniform magnetic field. The latter simultaneously enacts both requirements by creating an ergoregion (thereby opening up the possibility of superradiance) and trapping the scalar field in the black hole's vicinity. This leads to some novel features, in particular, that only black holes with a subset of the possible charge to mass ratios can support stationary clouds.

preprint2021arXiv

Spin-induced scalarized black holes

It was recently shown that a scalar field suitably coupled to the Gauss-Bonnet invariant $\mathcal{G}$ can undergo a spin-induced linear tachyonic instability near a Kerr black hole. This instability appears only once the dimensionless spin $j$ is sufficiently large, that is, $j \gtrsim 0.5$. A tachyonic instability is the hallmark of spontaneous scalarization. Focusing, for illustrative purposes, on a class of theories that do exhibit this instability, we show that stationary, rotating black hole solutions do indeed have scalar hair once the spin-induced instability threshold is exceeded, while black holes that lie below the threshold are described by the Kerr solution. Our results provide strong support for spin-induced black hole scalarization.

preprint2021arXiv

Spontaneous scalarization of a conducting sphere in Maxwell-scalar models

We study the spontaneous scalarization of a standard conducting charged sphere embedded in Maxwell-scalar models in flat spacetime, wherein the scalar field $ϕ$ is nonminimally coupled to the Maxwell electrodynamics. This setup serves as a toy model for the spontaneous scalarization of charged (vacuum) black holes in Einstein-Maxwell-scalar (generalized scalar-tensor) models. In the Maxwell-scalar case, unlike the black hole cases, closed-form solutions exist for the scalarized configurations. We compute these configurations for three illustrations of nonminimal couplings: one that \textit{exactly} linearizes the scalar field equation, and the remaining two that produce nonlinear continuations of the first one. We show that the former model leads to a runaway behaviour in regions of the parameter space and neither the Coulomb nor the scalarized solutions are stable in the model; but the latter models can heal this behaviour producing stable scalarized solutions that are dynamically preferred over the Coulomb one. This parallels reports on black hole scalarization in the extended-scalar-Gauss-Bonnet models. Moreover, we analyse the impact of the choice of the boundary conditions on the scalarization phenomenon. Dirichlet and Neumann boundary conditions accommodate both (linearly) stable and unstable parameter space regions, for the scalar-free conducting sphere; but radiative boundary conditions always yield an unstable scalar-free solution and preference for scalarization. Finally, we perform numerical evolution of the full Maxwell-scalar system, following dynamically the scalarization process. They confirm the linear stability analysis and reveal that the scalarization phenomenon can occur in qualitatively distinct ways.

preprint2021arXiv

Ultralight bosons for strong gravity applications from simple Standard Model extensions

We construct families, and concrete examples, of simple extensions of the Standard Model that can yield ultralight {real or} complex vectors or scalars with potential astrophysical relevance. Specifically, the mass range for these putative fundamental bosons ($\sim 10^{-10}-10^{-20}$ eV) would lead dynamically to both new non-black hole compact objects (bosonic stars) and new non-Kerr black holes, with masses of $\sim M_\odot$ to $\sim 10^{10} M_\odot$, corresponding to the mass range of astrophysical black hole candidates (from stellar mass to supermassive). For each model, we study the properties of the mass spectrum and interactions after spontaneous symmetry breaking, discuss its theoretical viability and caveats, as well as some of its potential and most relevant phenomenological implications {linking them to the} physics of compact objects.

preprint2020arXiv

Black holes with synchronised Proca hair: linear clouds and fundamental non-linear solutions

Recent studies have made key progress on the black hole/solitonic solutions of the Einstein-Proca system. Firstly, fully non-linear dynamical evolutions of the Kerr black hole superradiant instability, triggered by a Proca field, have shown the formation of a new equilibrium state, a spinning black hole with synchronised Proca hair. Secondly, non-linear evolutions of spinning Proca stars have established that they are dynamically stable, unlike their scalar cousins. Thirdy, separability of the Proca equation on the Kerr background has been achieved. Motivated by these results, in this paper we reconsider Kerr black holes with synchronised Proca hair. The separability of the Proca equation on the Kerr background allows us to examine the stationary Proca clouds in greater detail, in particular their dependence on the different quantum numbers. These stationary clouds occur at a set of existence lines in the Kerr parameter space, from which the black holes with synchronised Proca hair bifurcate. We construct the domain of existence of these black holes, comparing the fundamental states missed in the original study with the first excited states and with the cousin scalar model, giving illustrative examples of Kerr-like and non-Kerr-like BHs. In the vanishing event horizon limit, these hairy black holes connect to the fundamental states of spinning Proca stars, which include the dynamically stable solutions.

preprint2020arXiv

EHT constraint on the ultralight scalar hair of the M87 supermassive black hole

Hypothetical ultralight bosonic fields will spontaneously form macroscopic bosonic halos around Kerr black holes, via superradiance, transferring part of the mass and angular momentum of the black hole into the halo. Such process, however, is only efficient if resonant: when the Compton wavelength of the field approximately matches the gravitational scale of the black hole. For a complex-valued field, the process can form a stationary, bosonic field-black hole equilibrium state - a black hole with synchronised hair. For sufficiently massive black holes, such as the one at the centre of the M87 supergiant elliptic galaxy, the hairy black hole can be robust against its own superradiant instabilities, within a Hubble time. Studying the shadows of such scalar hairy black holes, we constrain the amount of hair which is compatible with the Event Horizon Telescope (EHT) observations of the M87 supermassive black hole, assuming the hair is a condensate of ultralight scalar particles of mass $μ\sim 10^{-20}$ eV, as to be dynamically viable. We show the EHT observations set a weak constraint, in the sense that typical hairy black holes that could develop their hair dynamically, are compatible with the observations, when taking into account the EHT error bars and the black hole mass/distance uncertainty.

preprint2020arXiv

Einstein-Maxwell-scalar black holes: the hot, the cold and the bald

The phenomenon of spontaneous scalarisation of charged black holes (BHs) has recently motivated studies of various Einstein-Maxwell-scalar models. Within these models, different classes of BH solutions are possible, depending on the non-minimal coupling function $f(ϕ)$, between the scalar field and the Maxwell invariant. Here we consider the class wherein both the (bald) electrovacuum Reissner-Nordström (RN) BH and new scalarised BHs co-exist, and the former are never unstable against scalar perturbations. In particular we examine the model, within this subclass, with a quartic coupling function: $f(Φ) = 1+αΦ^4$. The domain of existence of the scalarised BHs, for fixed $α$, is composed of two branches. The first branch (cold scalarised BHs) is continuously connected to the extremal RN BH. The second branch (hot scalarised BHs) connects to the first one at the minimum value of the charge to mass ratio and it includes overcharged BHs. We then assess the perturbative stability of the scalarised solutions, focusing on spherical perturbations. On the one hand, cold scalarised BHs are shown to be unstable by explicitly computing growing modes. The instability is quenched at both endpoints of the first branch. On the other hand, hot scalarised BHs are shown to be stable by using the S-deformation method. Thus, in the spherical sector this model possesses two stable BH local ground states (RN and hot scalarised). We point out that the branch structure of BHs in this model parallels the one of BHs in five dimensional vacuum gravity, with [Myer-Perry BHs, fat rings, thin rings] playing the role of [RN, cold scalarised, hot scalarised] BHs.

preprint2020arXiv

Electromagnetic dual Einstein-Maxwell-scalar models

Electromagnetic duality is discussed in the context of Einstein-Maxwell-scalar (EMS) models including axionic-type couplings. This family of models introduces two non-minimal coupling functions $f(ϕ)$ and $g(ϕ)$, depending on a real scalar field $ϕ$. Interpreting the scalar field as a medium, one naturally defines constitutive relations as in relativistic non-linear electrodynamics. Requiring these constitutive relations to be invariant under the $SO(2)$ electromagnetic duality rotations of Maxwell's theory, defines 1-parameter, closed $\textit{duality orbits}$ in the space of EMS models, connecting different electromagnetic fields in "dual" models with different coupling functions, but leaving both the scalar field and the spacetime geometry invariant. This mapping works as a solution generating technique, extending any given solution of a specific model to a (different) solution for any of the dual models along the whole duality orbit. We illustrate this technique by considering the duality orbits seeded by specific EMS models wherein solitonic and black hole solutions are known. For dilatonic models, specific rotations are equivalent to $S$-duality.

preprint2020arXiv

Rotating Axion Boson Stars

We construct and study rotating axion boson stars (RABSs). These are the spinning generalisations of the spherical gravitating solitons recently introduced in arXiv:1909.05515. RABSs are asymptotically flat, stationary, axially symmetric, everywhere regular solutions of the Einstein-Klein-Gordon theory, in the presence of a periodic scalar potential arising in models of axion-like particles. The potential is characterised by two parameters: the mass of the scalar field $m_a$ and the decay constant $f_a$. We present an overview of the solution space, for different values of $f_a$, and analyse some of their phenomenological properties. For large decay constants the solutions become identical to the standard spinning mini boson stars. For small decay constants, on the other hand, the solutions develop distinctive features. In particular, we analyse their compactness, the emergence of ergoregions, light rings and the distribution of stable and unstable equatorial timelike circular orbits, including the innermost stable circular orbit. We also establish the occurrence of violations of the strong energy condition for physical observers, for some RABSs. We observe some analogy between RABSs and spinning gravitating $Q$-balls.

preprint2020arXiv

Spherical electro-vacuum black holes with resonant, scalar $Q$-hair

The asymptotically flat, spherical, electro-vacuum black holes (BHs) are shown to support static, spherical configurations of a gauged, self-interacting, scalar field, minimally coupled to the geometry. Considering a $Q$-ball type potential for the scalar field, we dub these configurations $Q$-clouds, in the test field approximation. The clouds exist under a resonance condition, at the threshold of (charged) superradiance. This is similar to the stationary clouds supported by Kerr BHs, which exist for a synchronisation condition, at the threshold of (rotational) superradiance. In contrast with the rotating case, however, $Q$-clouds require the scalar field to be massive and self-interacting; no similar clouds exist for massive but free scalar fields. First, considering a decoupling limit, we construct $Q$-clouds around Schwarzschild and Reissner-Nordström BHs, showing there is always a mass gap. Then, we make the $Q$-clouds backreact, and construct fully non-linear solutions of the Einstein-Maxwell-gauged scalar system describing spherical, charged BHs with resonant, scalar $Q$-hair. Amongst other properties, we observe there is non-uniqueness of charged BHs in this model and the $Q$-hairy BHs can be entropically preferred over Reissner-Nordström, for the same charge to mass ratio; some $Q$-hairy BH solutions can be overcharged. We also discuss how some well known no-hair theorems in the literature, applying to electro-vacuum plus minimally coupled scalar fields, are circumvented by this new type of BHs.

preprint2020arXiv

Spinning black holes in shift-symmetric Horndeski theory

We construct spinning black holes (BHs) in shift-symmetric Horndeski theory. This is an Einstein-scalar-Gauss-Bonnet model wherein the (real) scalar field couples linearly to the Gauss-Bonnet curvature squared combination. The BH solutions constructed are stationary, axially symmetric and asymptotically flat. They possess a non-trivial scalar field outside their regular event horizon; thus they have scalar hair. The scalar "charge" is not, however, an independent macroscopic degree of freedom. It is proportional to the Hawking temperature, as in the static limit, wherein the BHs reduce to the spherical solutions found by Sotirou and Zhou. The spinning BHs herein are found by solving non-perturbatively the field equations, numerically. We present an overview of the parameter space of the solutions together with a study of their basic geometric and phenomenological properties. These solutions are compared with the spinning BHs in the Einstein-dilaton-Gauss-Bonnet model and the Kerr BH of vacuum General Relativity. As for the former, and in contrast with the latter, there is a minimal BH size and small violations of the Kerr bound. Phenomenological differences with respect to either the former or the latter, however, are small for illustrative observables, being of the order of a few percent, at most.

preprint2020arXiv

Stationary black holes and light rings

The ringdown and shadow of the astrophysically significant Kerr Black Hole (BH) are both intimately connected to a special set of bound null orbits known as Light Rings (LRs). Does it hold that a generic equilibrium BH must possess such orbits? In this letter we prove the following theorem. A stationary, axi-symmetric, asymptotically flat black hole spacetime in 1+3 dimensions, with a non-extremal, topologically spherical, Killing horizon admits, at least, one standard LR outside the horizon for each rotation sense. The proof relies on a topological argument and assumes $C^2$-smoothness and circularity, but makes no use of the field equations. The argument is also adapted to recover a previous theorem establishing that a horizonless ultra-compact object must admit an even number of non-degenerate LRs, one of which is stable.

preprint2020arXiv

Stationary scalar and vector clouds around Kerr-Newman black holes

Massive bosons in the vicinity of Kerr-Newman black holes can form pure bound states when their phase angular velocity fulills the synchronisation condition, i.e. at the threshold of superradiance. The presence of these stationary clouds at the linear level is intimately linked to the existence of Kerr black holes with synchronised hair at the non-linear level. These configurations are very similar to the atomic orbitals of the electron in a hydrogen atom. They can be labeled by four quantum numbers: $n$, the number of nodes in the radial direction; $\ell$, the orbital angular momentum; $j$, the total angular momentum; and $m_j$, the azimuthal total angular momentum. These synchronised configurations are solely allowed for particular values of the black hole's mass, angular momentum and electric charge. Such quantization results in an existence surface in the three-dimensional parameter space of Kerr-Newman black holes. The phenomenology of stationary scalar clouds has been widely addressed over the last years. However, there is a gap in the literature concerning their vector cousins. Following the separability of the Proca equation in Kerr(-Newman) spacetime, this work explores and compares scalar and vector stationary clouds around Kerr and Kerr-Newman black holes, extending previous research.

preprint2020arXiv

Tidal Love numbers of Proca stars

Proca stars (arXiv:1508.05395) are everywhere regular, asymptotically flat self-gravitating solitons composed of complex, massive vector bosons, forming a macroscopic, star-sized, Bose-Einstein condensate. They have been suggested as a classical, dynamical model of exotic compact objects, black hole mimickers and dark matter candidates. In spherical symmetry, they are qualitatively similar to their scalar cousins, the standard (scalar) boson stars in all respects studied so far. Here we study the tidal deformability and the quadrupolar tidal Love numbers, both electric-type and magnetic-type, of spherically symmetric Proca stars. The equations for the perturbations are derived and the numerical values of the Love numbers for some concrete background solutions are computed. We observe that both Love numbers are qualitatively similar to the ones of the scalar boson stars; in particular, the electric-type (magnetic-type) quadrupolar Love numbers are positive (negative) for Proca stars. Quantitatively, for the same compactness, the electric and magnetic Love numbers in the Proca case are closer (in magnitude) than in the scalar case, with the electric Love numbers having a slightly larger magnitude than the magnetic ones.

preprint2019arXiv

A class of solitons in Maxwell-scalar and Einstein-Maxwell-scalar models

Recently, no-go theorems for the existence of solitonic solutions in Einstein-Maxwell-scalar (EMS) models have been established in arXiv:1902.07721. Here we discuss how these theorems can be circumvented by a specific class of non-minimal coupling functions between a real, canonical scalar field and the electromagnetic field. When the non-minimal coupling function diverges in a specific way near the location of a point charge, it regularises all physical quantities yielding an everywhere regular, localised lump of energy. Such solutions are possible even in flat spacetime Maxwell-scalar models, wherein the model is fully integrable in the spherical sector, and exact solutions can be obtained, yielding an explicit mechanism to de-singularise the Coulomb field. Considering their gravitational backreaction, the corresponding (numerical) EMS solitons provide a simple example of self-gravitating, localised energy lumps.

preprint2019arXiv

Einstein and Eddington and the eclipse in Principe: Celebration and science 100 years after

On May 29, 1919, at Roça Sundy, Principe island, Eddington confirms Einstein's general relativity theory for the first time by photographing stars behind the obscured Sun during a total eclipse. History was made. At Sobral, Eddington's astronomer colleagues photograph the same eclipse and also conclude that light from distant stars suffers a deflection when passing by the gravitational field of the Sun, in accordance with general relativity. With the confirmation of general relativity, a theory of gravitation at a fundamental level, physics became, once and for all, relativist and its future was outlined. The first world war had finished a few months before and the deep wounds between nations were yet to heal. Science wanted to be above it all showing that people could be united by a common goal. This year, the 100 years of this achievement was commemorated with the scientific conference "From Einstein and Eddington to LIGO: 100 years of gravitational light deflection" in Principe to celebrate this landmark event. We report here on this conference of celebration.

preprint2019arXiv

On the inexistence of self-gravitating solitons in generalised axion electrodynamics

Building upon the methods used recently, we establish the inexistence of self-gravitating solitonic solutions for both static and strictly stationary asymptotically flat spacetimes in generalised axion electrodynamics. This is an Einstein-Maxwell-axion model, where the axion field $θ$ is non-minimally coupled to the electromagnetic field. Considering the standard QCD axion coupling, we first present an argument for the absence of static axionic solitons, $i.e.$ localised energy axionic-electromagnetic configurations, yielding an everywhere regular, horizonless, asymptotically flat, static spacetime. Then, for generic couplings $f(θ)$ and $g(θ)$ (subject to mild assumptions) between the axion field and the electromagnetic field invariants, we show there are still no solitonic solutions, even when dropping the staticity assumption and merely requiring a strictly stationary spacetime, regardless of the spatial isometries.