Gravitational waves at the first post-Newtonian order with the Weyssenhoff fluid in Einstein–Cartan theory

The generation of gravitational waves from a post-Newtonian source endowed with a quantum spin, modeled by the Weyssenhoff fluid, is investigated in the context of Einstein–Cartan theory at the first post-Newtonian level by resorting to the Blanchet–Damour formalism. After having worked out the basic principles of the hydrodynamics in Einstein–Cartan framework, we study the Weyssenhoff fluid within the post-Newtonian approximation scheme. The complexity of the underlying dynamical equations suggests to employ a discrete description via the point-particle limit, a procedure which permits the analysis of inspiralling spinning compact binaries. We then provide a first application of our results by considering binary neutron star systems.


Introduction
Nowadays, gravitational-wave (GW) astronomy represents a fundamental mean to investigate gravity at extreme regimes and offers valuable insights into the physics of compact objects [1,2], as the best GW candidates are represented by black holes (BHs) and neutron stars (NSs) [3]. GWs manifest as perturbations of the spacetime and their theoretical description lays solid roots in general relativity (GR). The GW theory is intimately intertwined with the two-body problem. A variety of analytical and numerical techniques has been developed to foretell approximately the dynamics and the corresponding waveforms of compact-object binary systems during their inspiral, plunge, merger, and ringdown stages [4,5]. The motion and the radiation of post-Newtonian (i.e., slowly moving, weakly stressed, and weakly self-gravitating) isolated sources during their early a e-mail: emmanuele.battista@univie.ac.at; emmanuelebattista@gmail.com (corresponding author) b e-mail: vittorio.defalco-ssm@unina.it inspiralling stage can be tackled via the Blanchet-Damour scheme. This framework is built on the pioneering works by Bonnor and collaborators [6][7][8][9] and Thorne [10], and exploits two approximation strategies: the multipolar-post-Minkowskian (MPM) and the post-Newtonian (PN) methods [11,12]. Other two fundamental GW generation formalisms are the Will-Wiseman-Pati approach, which extends the pattern first developed by Epstein and Wagoner [13], and the gravitational self-force (GSF) model. The former reckons with the direct integration of the relaxed Einstein equations (DIRE) and differs from the aforementioned MPM-PN program in the definition of the source multipole moments [14][15][16][17]. GSF is based on BH perturbation theory and explores the dynamics and the radiative phenomena of extrememass-ratio inspirals [18][19][20][21]. The effective-one-body (EOB) framework provides a highly-accurate description of the motion and the gravitational amplitude of coalescing binaries in their late evolution phases [22][23][24][25][26][27]. A valid support, especially for the demanding task of solving the Einstein field equations in the most extreme regimes, is provided by numerical relativity (NR). Indeed, NR simulations are firmly harnessed to predict waveforms through the merger and the ringdown, and validate the approaches used to study binary systems [28][29][30][31][32]. Furthermore, NR is largely exploited together with several phenomenological patterns (Phenom) to perform the fitting and the parameter estimation of the GW data [33][34][35][36][37].
The state of the art on the conservative dynamics consists in the complete set of the equations of motion of pointparticle nonspinning binaries at the fourth post-Newtonian (4PN) order, where the underlying calculations have been undertaken within three different patterns: the Arnowitt-Deser-Misner (ADM) Hamiltonian formulation of GR [38,39], the Fokker-action approach in harmonic coordinates [40,41], and the effective-field-theory model [42,43]; recently, the 5PN and 6PN levels have been worked out modulo a small number of unknowns [44][45][46][47]. In addition, this field is currently investigated via new methodologies making use of tools stemming from effective field theory and modern scattering amplitude programs [48][49][50]. Regarding the GW emission aspects, the 3PN-accurate waveform has been determined [51] and a great deal of efforts is being made in the literature to extend our present knowledge regarding GW templates of inspiralling compact binaries up to 4PN level [52][53][54]. In the case of spinning binary systems, both the radiative aspects and equations of motion have been worked out with high PN accuracy [55][56][57][58][59][60][61][62][63][64][65].
The huge amount of complementary observational data, triggered by the fast advance of the technological progress, is increasing our capacity to acquire more and more accurate information on gravitational sources and gravity itself. These scientific achievements have motivated us to explore the interplay between quantum and GR effects in GW phenomena via the Einstein-Cartan (EC) theory. Indeed, in a previous paper [66] we have solved, in the context of EC model, the GW generation problem at 1PN level by exploiting the Blanchet-Damour formalism, which permits relating the outgoing radiative gravitational field to the structure and the motion of a spinning PN isolated source. The final solution is encoded by the 1PN-accurate relations for the radiative moments, which are given in the form of well-defined (compact-support) integral expressions over the stress-energy distribution of the matter field. This result can be obtained after a detailed analysis regarding the gauge condition, the coordinates covering the interior and the exterior zone of the source, and the structure of the torsion tensor, which is supposed to admit a vanishing trace. Furthermore, the invariance of the Riemann tensor in the context of the linearized EC framework must be invoked in order to match the internal and the external fields and write the final expressions of the radiative moments involving the new spin contributions. Hereafter, "spin" is intended as the quantum (microscopic) angular momentum of elementary particles [78].
In this paper, we apply the abovementioned Blanchet-Damour scheme to a GW source described by the Weyssenhoff model of an ideal fluid with spin [67,68]. This choice is motivated by the fact that it configures as a natural extension of the GR perfect fluid and permits to simplify the demanding calculations framed in EC theory. The hydrodynamical picture of the Weyssenhoff treatment can be obtained by considering first a perfect fluid, and then by assigning to each "fluid element", which contains a set of microscopic spin configurations, a value of the spin density tensor via an average procedure [69,70]. The Weyssenhoff fluid has been considered both in cosmological models and astrophysical problems. In the first case, it has been proved that spin interactions in the early stages of the universe could bring about significant results, such as the avoidance of the big-bang singularity, as well as the reproduction of cosmic inflation and dark energy mechanisms [69][70][71][72]. In astrophysical settings, the spin may avert the spacetime singularity caused by the gravitational collapse of a star [73][74][75][76][77].
In our analysis, the dynamics of the Weyssenhoff fluid is investigated through the PN method, which permits obtaining a hierarchy of PN mathematical problems to be solved perturbatively order by order starting from the 0PN level. However, this approach still leads to a set of partial and integro-differential dynamical equations, which in general can be tackled through NR. Therefore, in order to obtain analytical results, we resort to the point-particle procedure, which is a valid pattern widely exploited in the literature [17]. In this way, the bodies can be treated as point-like objects and their dynamics is described in terms of ordinary differential equations. This prepares the ground for the analysis of inspiralling spinning compact binaries, which can represent the high-energy astrophysical testbed of our model.
The article is organized as follows: in Sect. 2, we briefly recall the fundamental results of Ref. [66]; in Sect. 3, after having set out the main concepts of the hydrodynamics in EC theory, we deal with the Weyssenhoff model; then, we pass from the continuous description to the discrete picture of the fluid by exploiting the point-particle limit, which is worked out in Sect. 4; in Sect. 5, we apply our theoretical apparatus to NS binaries, and provide a first estimate of EC corrections by studying the ensuing gravitational flux and waveform; eventually, in Sect. 6, we give a summary of the paper and draw the conclusions.

Einstein-Cartan theory
The EC theory is an extension of GR where both the spin and the mass of matter play a dynamical role. This model is defined on a spacetime M , endowed with a symmetric metric tensor g αβ and the most general metric-compatible affine connection whereΓ λ μν denotes the Christoffel symbols and dubbed the Cartan torsion tensor. Hereafter, a hat symbol refers to quantities framed in GR.
Given the matter Lagrangian density, it is possible to introduce the metric energy-momentum tensor T αβ , the spin angular momentum tensor τ βα γ , and the spin energy potential μ βα γ [66,78]. In particular, the tensors μ βα γ and τ βα γ are related in the following way: Another fundamental object is the total energy-momentum tensor of matter T αβ , defined as [78,79] where ∇ * is the modified covariant derivative operator, whose action on a generic tensor field of type (1, 1) is The EC field equations are (χ ≡ 16π G/c 4 ) [78] where the combined energy-momentum tensor Θ αβ satisfieŝ The matter source's dynamical equations can be obtained via the generalized conservation laws of energy-momentum and angular momentum, which read as, respectively, 1 [78,79] the Riemann tensor being given by

The mathematical problem
The GW generation in EC theory is described by the following well-posed mathematical problem [66]: where h αβ ≡ √ −gg αβ − η αβ , ≡ η αβ ∂ α ∂ β , and is the effective stress-energy pseudotensor encompassing both the matter fields, described by Θ αβ , and the effective gravitational source term Λ αβ = Λ αβ (h), which includes all the nonlinearities of EC field equations. The harmonic gauge (9b) can be imposed provided we further require [66] S αμ μ = 0.
The gravitational source is supposed to be confined to the region Ω, defined by withd the typical size of the source. This condition ensures, on the one hand, that T αβ , τ βα γ , and Θ αβ are smooth functions in R 4 having a spatially compact support in Ω, and, on the other, that the torsion tensor vanishes outside Ω.
We deal with PN sources in EC theory, for which the spatial domain R 3 can be decomposed as R 3 = D e ∪ D i , where the set D i is the near zone and covers entirely the source, while D e encompasses the external weak-field region of the source and is called exterior zone. For a PN source these two domains intersect in the overlapping region D o . Finally, the spatial region where a detector apparatus is located is known as wave zone.

Resolution method
The (approximate) solution of problem (9) can be worked out as follows [66,81]. First of all, two preliminary steps need to be carried out in the overlapping region D o : (1) calculating the multipole (re)expansion of the 1PN series of the inner metric, which can be expressed in terms of the source multipole moments I L and J L 2 ; (2) evaluating the 1PN (re)expansion of the MPM external metric, which depends on the canonical multipole moments M L and S L . Eventually, I L , J L and M L , S L are related to the radiative moments U L , V L (representing the physical observables in the wave zone) by exploiting the matching procedure. At 1PN order, the relation between I L , J L and U L , V L is given by (with the superscript (l) denoting the l-th time derivative with respect to the variable u, and I L (u) = d 3 y y L σ ( y, u) where and y L = q y L stands for the symmetric-trace-free (STF) projection of y L .

Asymptotic gravitational waveform and radiated power
Given a set X μ = (cT, X) of radiative coordinates, the external metric can be put in the so-called radiative form, where its coefficients admit an asymptotic expansion in powers of R −1 at future null infinity (i.e., R ≡ |X| → ∞ with U ≡ T −R/c and N ≡ X/R fixed) [12,[80][81][82][83]. The asymptotic waveform H TT i j (X ) describing the outgoing radiation is defined as the transverse-traceless (TT) projection of the leading R −1 term of such an expansion. At 1PN order, it reads as where P i jkl (N ) being the TT projection operator onto the plane orthogonal to N . Starting from the results contained in Refs. [17,83,84], we have proved that the standard GR formula of the total radiated power F (also called, in the astrophysics literature, total gravitational luminosity or flux of the source) is valid also in EC theory. Therefore, at 1PN order, the total energy radiated per unit time, expressed as a function of the retarded time U, reads as (1) U i j (U) (1) U i jk (U)

The semiclassical spin fluid and its post-Newtonian approximation
Having set out the main aspects of the Blanchet-Damour approach in EC theory, we consider, as a first approach to the description of spin effects inside matter, the class of semiclassical spin fluid models. Before getting to the heart of the discussion, in Sect. 3.1 we first present the main principles of the hydrodynamics in EC theory by making use of the hypothesis (11). In Sect. 3.2, we consider the Weyssenhoff fluid, which represents one of the most common frameworks studied in the literature. Finally, in Sect. 3.3, we investigate the Weyssenhoff fluid dynamics within the PN approximation scheme.
3.1 Hydrodynamics in Einstein-Cartan theory under the hypothesis S αβ β = 0 In this section, we introduce the basic concepts underlying the hydrodynamics in EC theory. The analysis of a fluid in EC framework turns out to be more complex than in the GR case, but the assumption (11) entails a great simplification. Indeed, the modified covariant derivative (4) becomes and the EC covariant divergence of a generic vector A μ assumes the same form as in GR, namely The spacetime position of each fluid element is labelled in a system of coordinates by x μ = (ct, x). The description of the kinematic properties of the fluid can be performed in terms of two fundamental quantities: the timelike fourvelocity u μ = dx μ dλ (where λ denotes the proper time of an observer comoving with the fluid) and the four-acceleration It is useful to introduce the projection operator on the spatial hypersurface orthogonal to the timelike four-velocity u μ and the EC substantial derivative of the vector ξ μ along u μ [67] ξ μ ≡ u ν ∇ ν ξ μ .
The latter is the natural extension of the concept used in GR and in classical continuous mechanics, which allows to define the EC time derivative for densities [67] Dξ μ ≡ ∇ ν u ν ξ μ =ξ μ + ξ μ∇ ν u ν .
The above derivative operator is applied to objects which are densities (e.g., the spin density tensor, which will be introduced in Sect. 3.2) and has a precise physical meaning. Indeed, its definition relies on the fact that the densities must be not only transported along the fluid worldlines, but it must be also taken into account how the volumes transform during the dynamical evolution 3 (see Ref. [67], for more details). 3 By defining the expansion tensor ϑ αβ ≡ P μ α P ν β ∇ (ν u μ) , it is easy to prove that the expansion scalar ϑ ≡ ϑ α α assumes the same form as in GR [85], i.e., ϑ = 2∇ α u α . This quantity describes how the fluid volume changes during the motion.
The metric stress-energy tensor modeling the relativistic perfect fluid which includes also spin contributions, can be generally written as where Φ αβ is a symmetric tensor containing torsion terms and whose explicit form depends on the chosen spin model, P denotes the isotropic fluid pressure, and the total energy density, which, in turn, depends on the restmass density ρ and the internal energy density ε. We note that e can be functionally split in the sum of ρc 2 and ε due to the perfect fluid hypothesis. Upon introducing the rest-mass density current we can write the conservation equation for the rest mass as (cf. Eq. (20)) Furthermore, Eq. (19) permits simplifying the expressions of the canonical stress-energy tensor (3) and of the combined stress-energy tensor (5b), which read as, respectively, Therefore, from Eq. (7a), we have The dynamics of the perfect fluid can be examined via Eqs. (27) and (29), supplemented by Eq. (7b). In this way, we obtain a system of highly-non-linear differential equations, dubbed cardinal equations of the hydrodynamics in EC theory, which can be written as Equation (30b) can be evaluated in the directions orthogonal (by means of the projection operator (21)) and parallel to the fluid four-velocity. The former components give the Euler equation in EC theory, while the latter the energy-balance law. Furthermore, Eq. (30c) leads to the rotational equations of the fluid motion.
The system (30) comprises, in general, 11 independent equations and at most 26 unknowns, which are represented by: the 3 components of the fluid four-velocity (the fourth one is constrained by the normalization condition u μ u μ = −c 2 ), the rest-mass density ρ, the internal energy density ε, the pressure P, and (at most) the 20 independent components of the spin angular momentum tensor τ μν λ . Indeed, the independent components of the torsion tensor S μν λ have been lowered to 20 by virtue of the hypothesis (11) [66] and hence also τ μν λ has (at most) 20 independent components. This is in agreement with the physical content of the EC theory, where the torsion tensor is the geometrical counterpart of the spin of matter, encoded by the tensor τ μν λ . In the worst case, S μν λ can be determined once the 20 components of τ μν λ are known, but, as we will see in the next section, this number can be drastically reduced by assigning a suitable functional form for τ μν λ . It is clear that the system (30) is in general not closed, since there is no balance between the number of independent equations and unknowns. However, this not a drawback of the model, which, on the contrary, can reveal a rich dynamical structure. In order to supply the missing equations and correctly characterize the structure of the fluid under investigation, constitutive equations must be added to the system (30), e.g., P = P (ρ, ε), P = P(ε), ρ = ρ(ε).
A fundamental and general aspect of our approach relies on the possibility of describing the spin effects inside the fluid by exploiting all kinds of (relativistic) spin models. Indeed, the heart of this pattern consists in providing the functional form of Φ αβ and τ αβγ , which is equivalent to assign T αβ and μ αβγ . In general, these two unknowns play a crucial physical and geometrical role in the EC gravity framework. Indeed, we can distinguish the following two classes: (1) spin-geometry, represented by either τ αβγ or μ αβγ (skewsymmetric quantities), which rule the torsion field S μν λ ; (2) mass-geometry, described by either Φ αβ or T αβ (symmetric quantities), which shape the metric tensor g μν .

The Weyssenhoff fluid
This section is devoted to the description of the Weyssenhoff fluid. In Sect. 3.2.1, we retrace the historical ideas behind Weyssenhoff approach, since they are useful for its full comprehension. In Sect. 3.2.2, we present the model and the related dynamical equations. We will see that it can be readily analyzed within our general framework put forth in Sect.
3.1. A short digression on the first thermodynamic law is contained in Sect. 3.2.3.

Historical introduction to the fluid with spin
The study of matter in its microphysical aspects can be tackled following different strategies framed either in quantum models or from the standpoint of classical vision. Due to the arguments developed in this paper, it is more suitable to follow the latter approach, where the structure of physical systems is derived starting from relativistic theories. We consider what in the early literature was dubbed free spin-particle, namely a material particle endowed with spin and on which no force, apart from the gravitational pull, acts [67,86,87]. We pursue the route paved by the "Krakow school" in the years 1937-1947, characterized by the prolific scientific activity of several Polish physicists, like M. Mathisson, J. Lubański, J. Weyssenhoff, and A. Raabe. In 1927, Einstein and Gromer derived the dynamical equations of a free particle from the equations of the gravitational field as the singularities in this field [88]. In 1937, inspired by these ideas, Mathisson and Lubański deduced the equations of motion of a free spin-particle through a variational principle in the case of a linearized gravitational field [86,87,89,90]. This new dynamics exhibits two peculiarities: (1) in the classical limit, it does not reduce to the Newton laws of motion, since it keeps an additional term depending on the internal angular momentum or spin of the particle; (2) the presence of spin generates dynamical equations of the third order. However, Mathisson did not immediately realize that his results were equivalent to those previously published in 1926 by the Russian physicist Frenkel [91]. Indeed, the latter author, after a discussion with Pauli, who showed him a letter by Thomas on the famous "precession factor 1/2" [92], became interested in deriving the equations of motion for a spinning electron, which coincide exactly with those of a free spinparticle only if the terms depending on the electromagnetic field vanish. In 1947, after the Second World War, Weyssenhoff and Raabe proposed a third different method to obtain the equations of a free spin-particle. These authors showed that their results agree with those derived by Frenkel and Mathisson [67]. However, their approach is more rigorous, since it is constructed by exploiting fundamental concepts from continuous mechanics and GR. In addition, the final equations are presented in a much simpler form thanks to the introduction of the linear energy-momentum four-vector [67].
The last approach is the most used in the literature to derive the equations of the incoherent spinning fluid. For this reason, it is also called Weyssenhoff(-Raabe) fluid. Having defined s αβ as the spin density per unit rest-volume, we split it into two three-dimensional vectors The fundamental hypothesis behind the Weyssenhoff fluid is that the vector q vanishes in the rest system of the fluid, which translates in having the covariant relation [67,91] since the fluid velocity u α = (u 0 , u) has vanishing spatial components when evaluated in the fluid rest frame. This hypothesis, also known in the literature as Frenkel condition since it was first employed by Frenkel, arises from the following practical needs: (1) simplifying and closing the set of differential equations underlying the dynamics by matching the numbers of unknowns and equations; (2) letting s be the vector which encodes the real physical degrees of freedom, since q can be gauged away by an appropriate choice of a coordinate system. Equation (32) can be written in three-dimensional form as from which it follows The vectors s and q have the same transformation properties of the magnetic and the electric fields, respectively. Therefore, the Frenkel condition physically tells us that s αβ is "purely magnetic", meaning that its electric component vanishes in the rest-frame coordinate system.
Following the general approach devised in Sect. 3.1, a specific fluid can be characterized by assigning the spin angular momentum tensor τ αβ γ and the symmetric tensor Φ αβ . The former, is given by [69] s αβ = s [αβ] being the spin density tensor, which is constrained, due to the hypothesis (11), to satisfy i.e., the Frenkel condition (cf. Eq. (32)). This assumption permits several considerable simplifications. First of all, the torsion and contortion tensors can be written as respectively; furthermore, the EC time derivative of s μν does not involve contortion terms and reads as (cf. Eq. (23)) finally, we can write the useful identity where the acceleration a μ is the same as in GR, i.e., The tensor Φ αβ of the Weyssenhoff fluid is given by [69] where s 2 ≡ s αβ s αβ is the spin density scalar, satisfying the condition Ds 2 = 0. Therefore, from Eq. (41), jointly with Eqs. (5c), (24), and (35), we have and hence the explicit expressions of the canonical and the combined energy-momentum tensors are, respectively, (cf. Eq. (28)) The general form of the four-momentum density is where α is a spacelike vector orthogonal to u α . An inspection of Eq. (44a) suggests α = −2a σ s σ α /c 2 and hence upon exploiting the relation e = −p α u α [69]. By substituting Eq. (46) in Eq. (44a), we recover the usual form of the canonical energy-momentum tensor for the Weyssenhoff fluid [78], i.e., The internal structure of the fluid is described in terms of the rest-mass density ρ and the internal energy density ε. The total energy density e can be written as (cf. Eq. (25)) where s denotes the specific entropy (i.e., the entropy per unit mass). The functional dependence of ε partially resembles that of classical physics, apart from the presence of s μν , which embodies the new contribution due to quantum mechanical effects and is specific of the EC theory. We suppose that the fluid is adiabatic and that the rest mass of the system is conserved. The first hypothesis can be written as [85,97] whereas the second one is represented by Eq. (30a). The translational equations pertaining to the fluid dynamics can be obtained from Eq. (30b), along with Eqs. (2), (35), (37), and (41). As pointed out before, the projection along the fluid four-velocity u μ gives the energy-balance law (or energy-conservation equation), which upon exploiting Eq. (30a) reads aṡ whereas the projection onto the hypersurface orthogonal to u μ gives the Euler equation, which after some algebra can be written as We note that, for s μν = 0, Eq. (51) reduces to the Euler equation of GR [17,85]. Bearing in mind Eqs (30c), (35), (37), and (44a), we find that the rotational fluid motion is ruled by where we have exploited Eq. (46). Since the tensor s μν has in general six independent components, the Frenkel condition (36) is crucial to correctly account for the three true dynamical degrees of freedom associated with the spin of a particle. For this reason, Eq. (52) gives rise to three independent equations.
As discussed in Sect. 3.1, the system of the cardinal hydrodynamic equations can be closed once the needed constitutive equations are provided. Therefore, the full set of differential equations governing the Weyssenhoff fluid dynamics is (see Eqs. (27) and (49)-(52)) The above system contains ten independent equations and exactly ten unknowns, i.e., the three components of the fluid velocity u μ , the rest-mass density ρ, the internal energy density ε, the specific entropy s, the isotropic pressure P, and the three independent components of s μν . It is also worth noting that it explicitly includes an equation describing the behaviour of s (cf. Eq. (53e)). This is due to the fact that, unlike GR [17,85], in EC theory Eqs. (53a) and (53d) do not imply that the fluid evolves adiabatically (see Sect. 3.2.3).
In conclusion, the closed self-consistent system of gravitymatter equations is represented by EC field equations (5) supplemented by the set (53).

The first thermodynamic law
In the study of the hydrodynamics, it is important to investigate the implications of the energy-conservation equation (53d) at thermodynamic level [69,85]. However, we note that this analysis does not provide additional constraining equations to the system (53), but it allows to better understand how the fluid behaves in EC theory.
The first thermodynamic law for a spinning fluid can be written as [69,93,96,98] where is the specific internal energy, and θ being the temperature and ω μν the thermodynamic variable conjugated to s μν coinciding with the microscopic angular velocity of the fluid. Starting from the continuity Eq. (53a) and the energybalance law (53d), we obtaiṅ which, once compared with Eq. (54), gives Bearing in mind Eq. (53e), it follows that At this point, some remarks are in order. First of all, as pointed out before, it is known that in GR a perfect fluid evolves adiabatically [17,85]. This result can be proved by following the same procedure as the one we have adopted in this section, i.e., by exploiting the fluid continuity equation jointly with the energy-balance law. On the contrary, in EC theory the same calculations do not imply the adiabatic-flow hypothesis as long as the spin contribution in Eq. (58) does not vanish. This shows that in EC hydrodynamics the adiabatic condition must be explicitly imposed via Eq. (53e), which in fact has been employed to obtain the final relation (59).

Post-Newtonian analysis of the Weyssenhoff fluid dynamics
In this section, we perform the PN investigation of the Weyssenhoff fluid dynamics, which is useful both for the calculations carried out in Sect. 4

The computational program
The EC hydrodynamics, although equipped with the simplifying constraint S αβ β = 0, gives rise to a set of coupled nonlinear differential equations (cf. Eq. (30)). Also in the more specific case, where we consider the Weyssenhoff fluid endowed with the Frenkel condition (36), the nonlinear nature and the differential structure still endure (cf. Eq. (53)), and we are, in general, not able to solve the mathematical problem exactly. An analytical solution may be determined by imposing symmetry requirements, like time independence and spatial isotropy (as it occurs, for example, in GR for the Schwarzschild metric). However, this approach reveals to be very limited in its range of applications, where time variability and spatial anisotropies occur frequently.
Since the EC gravity framework is geometrically more tangled than GR theory, our goal is therefore not to search for further exact solutions, but rather to take advantage of a comprehensive and solid approximation method, which goes beyond the symmetries and the particular functional form of the starting problem. The mathematical formulation of the GW generation problem naturally suggests us to resort to the perturbation theory (see Refs. [99][100][101], for further details and examples). This approach permits to determine a solution under the form of a power series of a small parameter (represented, in our case, by = 1/c), where expressions of higher powers of become smaller as the order increases. The approximate perturbed solution is obtained by truncating such a series at a certain established order, suggested either by the observational sensitivity of the phenomenon under investigation or, from a theoretical point of view, by the computational complexity (i.e., gradually increasing the PN order calculations).
The strategy to build up the solution via the PN formalism relies on expanding all the involved quantities, and then extracting the coefficients of the expansions which are useful to construct the differential equations at each PN level. In this way, we obtain a series of PN problems, where the lowest order permits obtaining the classical equations with eventual spin corrections, and the higher terms provide the successive perturbative corrections. The method is therefore iteratively based on first determining (either analytically or numerically) the solution of the 0PN problem, and then proceeding order by order to solve (either analytically or numerically) the successive PN problems. Indeed, to compute the N th PN order, we need to use the parameters determined in the previous iterations. Finally, gathering all these quantities together we can provide an approximate solution to the original problem up to the desired PN level.
It is important to remark that we apply this computational program in the derivation of the 0PN and 1PN orders of Eq. (53) without solving them. As we will see, the obtained results will be useful in Sect. 4.

Preliminary: expansions of basic quantities
Given the harmonic coordinates x μ = (ct, x) (cf. Eq. (9b)), the 1PN series of the metric tensor, expressed in terms of the retarded potentials V and V i , reads as [66] where (cf. Eq. (15)) The expansion of Eq. (61) for small retardation effects gives the potentials U , U i , and the superpotential X being given by respectively. It should be noted that the term of order c −3 in Eq. (62a) is a function of time only. Stating from Eq. (60), one easily obtains the PN structure of the Christoffel symbols (see Ref. [99], for details.) The spin angular momentum tensor τ μν λ admits, in general, the following PN series [66]: where (n) τ μν λ ∼M d 2v n c n−2 (v andM being some typical internal velocity and mass of the source, respectively). Therefore, the tensor (5c) has the PN structure Bearing in mind the definition of the spin angular momentum of a test particle in EC theory (see Eq. (40) of Ref. [66]), for the spin density tensor we find where (n) s μν indicates a factor going likeMv n d 2 c n−1 . If we write the fluid four-velocity equivalently as (with γ ≡ u 0 c and v ≡ dx dt ) its PN expansion can be easily constructed by means of Eq. (60) and the normalization condition. Therefore, the Frenkel condition (36) leads to where (n) u μ denotes a term of orderv n c n−1 . We note that Eq. (68a) results from a linear combination of Eq. (68b) with coefficients v j /c (recalling that (n) u j = (n−1) u 0 c v j with n = 1, 3). This implies that the latter gives rise to three independent relations which can be used to gauge away the components s 0i of the spin density tensor at different PN orders in the equations of motion (cf. Sect. 3.2.1 and the discussion below Eq. (52)). Indeed, from Eq. (68b) we obtain to leading order and to next-to-leading order, respectively, or, equivalently, From the above equations it is clear that no components s 0i are present at leading order, whereas they can be written in terms of (1) s ik at next-to-leading order. Furthermore, we see that Eq. (68a) simply states that the three-vectors s 0i and u i are perpendicular, but gives no information about the form of s 0i . On the contrary, Eq. (68b) allows to write explicitly, at each PN level, an expression for s 0i which turns out to be both orthogonal to u i and vanishing in the fluid rest frame. For our forthcoming calculations, it is useful to express Eq. (53a) in an equivalent form. Indeed, bearing in mind Eq. (67), and upon introducing the rescaled mass density ρ ≡ γ √ −gρ [17], Eq. (53a) yields the (exact) equation the relation between ρ and the proper rest-mass density ρ being The PN expansion of the first thermodynamic law (54) reads as (cf. Eq. (55)) where d dt and we have exploited Eq. (72), the relations ω 0i = O c −1 and ω i j = O c 0 (see Refs. [69,109] for details), and we have supposed that both θ and s are O(c 0 ) quantities.

Equations of motion: 0PN expansion
The results contained in Sect. 3.3.2 allow us to study the hydrodynamic equations (53) within the PN formalism. In our investigation, the rest-mass conservation Eq. (53a) will be replaced by its equivalent expression (71) and will not be expanded. Furthermore, by means of Eq. (67), Eq. (53e) assumes the exact form and the constitutive Eq. (53f) will be simply written in terms of ρ (cf. Eq. (72)). Therefore, by computing the 0PN expansion of the remaining equations of the system (53) and by exploiting the Frenkel condition (70a), we end up, after a lengthy calculation, with the following system pertaining to the dynamics of a Weyssenhoff fluid with 0PN accuracy: Equation (76b) originates from the 0PN expansion of Eq. (53b) and represents the Euler equation of the Newtonian theory. Equation (76c) is the leading-order piece of Eq. (53c). (1) s i j . Finally, Eq. (76d) is the 0PN-accurate energy-conservation Eq. (53d).

The point-particle limit of the Weyssenhoff fluid
The final outcome of the PN treatment contained in Sect. 3 is represented by Eqs. (76) and (80), which completely determine the behaviour of the Weyssenhoff fluid at 0PN and 1PN level, respectively. In general, these two sets comprise both partial and integro-differential equations, whose resolution, either analytically or via numerical means, turns out to be highly challenging. This makes, as a consequence, the evaluation of the radiative moments parametrizing the asymptotic waveform and the radiated power (cf. Sect. 2.4) rather demanding.
A widely used approach in the literature which offers a way around this issue consists in employing the so-called point-particle procedure, where the fluid configuration is supposed to be described by a collection of separated pointlike components, usually refereed to as "bodies". Within this pattern, the equations underlying the fluid dynamics become less complex as they are turned into ordinary differential equations. For these reasons, in Sect. 4.1, we introduce the point-particle procedure and work out the 1PN-accurate expressions of the radiative moments. As we will see, the results of Sect. 3 will be crucial, since they allow us to go from the fine-grained description to the coarse-grained picture of the Weyssenhoff fluid. Finally, in Sect. 4.2, we deal with the special case of binary systems, which are known to be the main candidates of GW events in high-energy astrophysics.

The point-particle procedure
The evaluation of the 1PN-accurate total power of emission (18) requires the knowledge of the 1PN radiative mass quadrupole moment U i j as well as the 0PN mass octupole U i jk and current quadrupole V i j ; on the other hand, the 1PN asymptotic waveform (16) can be worked out once the 0PN mass 2 4 -pole U i jkl and current octupole V i jk are also known. Having solved the 1PN GW generation problem [66], we know that these are related to the source multipole moments via Eqs. (13) and (14). Therefore, in order to compute the point-particle limit of the required radiative moments, we first need to determine their fine-grained expression, i.e., the form assumed in the case of a continuous distribution of matter represented by the Weyssenhoff fluid, see Sect. 4.1.1. Then, after some premises outlined in Sects. 4.1.2 and 4.1.3, the point-particle limit is performed in Sect. 4.1.4.

The fine-grained form of the radiative moments
In this section, we calculate the fine-grained expression of the radiative moments occurring in the total power of emission and the wave amplitude for the specific case of the Weyssenhoff fluid. This entails, first of all, the computation of the PN expansion of metric stress-energy tensor (42) and the tensor (43). With the help of the results of Sect. 3.3.2, the former yields whereas the latter where we have exploited Eqs. (62a) and (70) where σ and σ i have been defined in Eq. (15), while σ i j ≡ Θ i j [66], and, inspired by Ref. [102], with (see Eq. (63a)) It should be noticed that the spin contributions in Eq. (83) appear both explicitly, via the terms involving s i j , and implicitly, through, e.g., the specific internal energy Π (see Eq. (55)). For our forthcoming analysis, it is useful to define a new set of STF mass-type radiative moments I rad L and currenttype radiative moments J rad L according to Bearing in mind the above equations jointly with Eqs. (13), (14), and (83), the 1PN fine-grained form of I rad i j , and the leading-order expressions of I rad i jk , I rad i jkl , J rad i j , and J rad i jk are given by, respectively, where we have employed the Gauss theorem to discard the integrals containing a total divergence (recall that, in our model, both the spin angular momentum tensor (35) and the metric energy-momentum tensor (42) have compact support) and, in Eq. (91), the identity

The ADM mass and the center of mass
In Ref. [66], we have demonstrated that the mass monopole moment I , which can be read off Eq. (14a) with l = 0, gives rise to a generalized notion of ADM mass (or total massenergy) of the fluid system and can be written as and, in addition, satisfies In the case of the Weyssenhoff fluid, we find, with the help of Eq. (83), where we have exploited Gauss theorem to discard total derivatives occurring in σ and σ ii (cf. Eqs. (83a) and (83c)). For this reason, Eq. (96) resembles formally the GR expression [17]. The ADM mass (96) assumes the equivalent form where and denote the total material mass, the translational kinetic energy, the total gravitational potential energy, and the total internal energy of the fluid system, respectively. It should be noted that the contribution of the spin kinetic energy is included in Eq. (99d). Furthermore, Eqs. (71), (76b), and which can be obtained through Eqs. (57) and (72), assure the constancy of the ADM mass (97), in agreement with Eq. (95). The case l = 1 of Eq. (14a) yields the definition of the dipole moment, which for the Weyssenhoff fluid becomes where we have exploited the Gauss theorem.
The position of the center of mass of the fluid system, which is defined as R i ≡ I i /M ADM [17,80], reads as and satisfies the relation P i being the conserved ADM three-momentum of the system (see Ref. [66] for further details). It should be noted that an explicit spin correction term occurs in Eq. (102). This is similar to the contribution appearing in the analysis of spinning bodies within GR [17,103]. However, we recall the basic difference between EC model and Einstein theory, as in the former the spin refers to the intrinsic quantum angular momentum of elementary particles, while in the latter the spin is related to a macroscopic rotation [78].

The coarse-grained description
We are now ready to apply the results of the previous sections to the framework where the fine-grained description of the Weyssenhoff fluid is replaced by the coarse-grained picture involving separated fluid elements, which are characterized by a small number of variables [17,102]. Therefore, we consider a setup where the fluid distribution breaks up into a collection of N weakly self-gravitating, slowly moving, widely separated, and isolated spinning components surrounded by vacuum regions of space. Each component, referred to as body, is assigned a label A = 1, 2, . . . , N , and, in our hypotheses, the ratio α ≡ λ/δ 1 (λ and δ being the typical size of A and the typical separation between the bodies, respectively). The coordinate fluid density can be expressed as ρ = A ρ A , where the sum extends over each body, and ρ A vanishes everywhere except within the volume occupied by the body A. Furthermore, the (conserved) material mass of A is m mat A = A d 3 x ρ , the domain of integration being a time-independent portion of the three-dimensional space which extends slightly beyond the volume occupied by the body A and does not include nor intersect another body within the system (recall that ρ = ρ A within such a domain).
Similarly, the spin density tensor s i j is zero in the vacuum exterior of the bodies, whereas s i j = s A i j inside the volume of the body A. The spin vector s A of A is defined by and, owing to Eq. (76c), is conserved modulo O c −2 corrections, i.e., d dt Driven by Eqs. (102) and (96), the center of mass worldline r i A (t) of the body A is defined as where is the total mass-energy of the body, which, along the same lines as for M ADM , can be shown to be conserved. In Eqs. (106) and (107), Π A denotes the specific internal energy of A; furthermore, u A is the internal (self-gravity) potential, which is readily obtained from the Newtonian potential (86); finally, upon introducing the vector y i A ≡ x i − r i A (t), measuring the position of a fluid element relative to the center of mass r i , and represents the velocity of this fluid element relative to the body velocity v i The above definitions permit renormalizing the internal selfgravity of A into the mass m A .

The point-particle limit of the radiative moment
At this stage, we have all the ingredients needed to perform the point-particle limit of the radiative moments (88)- (92). For this reason we employ the same techniques as in Ref. [102].
In EC theory, the following crucial identities hold: where P A is the pressure within body A and all functions are supposed to depend on t and the position y A + r A (t). Equation (108a) is a consequence of Eqs. (106) and (107); Eq. (108b) stems from Eq. (108a) and, in turn, Eq. (108c) can be obtained by evaluating the time derivative of Eq. (108b); finally, Eq. (108d) derives from the virial theorem, which reads as (cf. Eq. (83)) 1 2 By supposing that all bodies are spherically symmetric and in static equilibrium, we find that the point-particle counterpart of the radiative mass quadrupole moment (88) reads as moreover, for the the mass octupole and current quadrupole, the point-particle procedure yields, respectively, the mass 2 4 -pole and the current octupole give, respectively, Equations (110)-(112) are required for the computation of the instantaneous luminosity (18), whereas the whole set of radiative moments (110)-(114) appears in the asymptotic waveform (16) (cf. Eq. (87)). It is clear that EC theory brings in explicit corrections due to the spin s A of the body A in I rad i j , J rad i j , and J rad i jk . We recall that, in deriving Eqs. (110)-(114), we have neglected terms O α 2 .

Binary systems
In this section, we set out the features of the special case of binary systems, i.e., a collection of N = 2 weakly selfgravitating, slowly moving, and widely separated spinning fluid bodies having masses m 1 , m 2 (with m 1 ≥ m 2 ), position vectors r 1 , r 2 , velocities v 1 , v 2 , and spin vectors s 1 , M , and ν ≡ μ M denote the total mass, the reduced mass, and the symmetric mass ratio of the system, respectively. As in the Newtonian framework [17,104], the study of the dynamics is simplified once the origin of the coordinate frame is attached to the barycenter R i of the system and the position of each body is determined in terms of their separation vector. Accordingly, we introduce the instantaneous relative position vector R and the instantaneous relative velocity vector V 4 of the two objects Moreover, starting from Eq. (102), we find that, within the coarse-grained description of the fluid and resorting to the same techniques as in the last sections, the position of the barycenter of a binary system is defined by (we discard, like before, O α 2 corrections) By employing a post-Galilean transformation (i.e., a particular subclass of general PN transformations [17]), it is always possible to define the center of mass frame of the system by setting R i = 0 and P i = 0. Therefore, by means of Eq.
(116), we find that in a mass-centered coordinate system the motion of the bodies is related to their relative motion by the following relations: If we replace Eq. (117) in Eqs. (110)-(114), we find that the general form of the radiative moments for binaries of spinning objects reads as The above equations completely determine the 1PN generation of GWs from binary systems in EC theory.

First application to binary neutron star systems
The theoretical pattern developed in the previous sections will be now applied to the study of binary NSs. A full analysis requires the knowledge of the 1PN dynamics in EC the-ory, which, at the moment, is not at our disposal. Despite that, we can provide a first estimation by following a hybrid approach. We exploit the EC definition of center of mass and the general expression of the EC radiative moments obtained in Sect. 4.2, along with the conservation equation of the spin vector (cf. Eqs. (76c) and (105)); furthermore, we consider the quasi-elliptic 1PN-accurate GR motion of a binary system determined by Damour and Deruelle [104] (see Sect. 5.1).
We note that, in the full EC framework, only the time derivatives of the radiative mass quadrupole moment I rad i j will contain new O(c −2 ) spin-dependent terms, whereas the other moments will remain unaffected. Indeed, as pointed out before, the EC 0PN-accurate translational motion coincides with the Newtonian Euler equation owing to the Frenkel condition (cf. Eq. (76b)).
In this hybrid setup, we obtain the explicit expressions of the flux and the gravitational waveform (see Sect. 5.2). We conclude the section with a numerical estimate concerning the EC corrections by examining binary NS systems (see Sect. 5.3).
Any effect of GW back-reaction on the source dynamics will be neglected. This hypothesis, which is valid to a good approximation in some astrophysical GW sources (see e.g. Refs. [105,106]), permits to derive a first model for the description of GW phenomena in EC theory. Hereafter, we no longer mention O c −n terms.

The Damour-Deruelle solution
We consider, within the hybrid approach set forth above, a binary system consisting of two PN widely separated spinning bodies having masses m 1 , m 2 and spin vectors s 1 , s 2 . As pointed out before, we assume that their dynamics is governed by the Damour-Deruelle solution, which we now briefly outline.
The problem of solving the motion of the binary system can be reduced to the simpler equivalent task of determining the relative motion in the PN center of mass frame (which can be defined by means of the results of Sect. 4.2). For this reason, given a harmonic coordinate system, we define, in the same way as before, the instantaneous relative position vector R and the instantaneous relative velocity vector V (cf. Eq. (115)).
The total energy E and the total angular momentum J of the system (which are conserved in the Damour-Deruelle dynamics) read as, respectively, The motion takes place in the plane orthogonal to J, which is henceforth supposed to be directed along the z-axis. We can thus introduce polar coordinates (R, ϕ) and write the equation of the PN relative orbit as where ϕ in is the initial angle and the following orbital parameters have been introduced:

The flux and the gravitational waveform
It follows from the definition (87) that the instantaneous luminosity (18) attains the equivalent form (restoring, for a while, the O(c −n ) terms) I rad i j + whereas the asymptotic amplitude (16) becomes I rad kl (u) J rad l)a (u) I rad klab (u) where u ≡ t − |x|/c, n ≡ x/|x|, and the radiative moments have been derived in Eqs. (118)-(122). We also note that in Eqs. (127) and (128) we have exploited the fact that, at this order, there is no difference between the harmonic and the radiative coordinates [80,81]. At this stage, we can compute F and H TT 11 by exploiting the Damour-Deruelle solution.
The instantaneous luminosity (127) can be written as where the expression of the GR flux can be found in Ref. [102], while the EC contribution can be obtained after a lengthy calculation by exploiting the results of Sects. 4.2 and 5.1 . In the hypotheses that the motion occurs in the x y-plane (i.e., R z = V z = 0) and the spins of the two bodies are aligned with the total angular momentum J (i.e., s x1 = s x2 = s y1 = s y2 = 0), we find for the EC luminosity The general form of the gravitational waveform can be obtained after a long calculation by considering the time derivatives of the radiative moments appearing in Eq. (128). In this section, we give the component H TT 11 of the asymptotic amplitude in the hypotheses that the GW propagates along the direction n = (0, 0, 1) and, like before, the spins of the two bodies are aligned and orthogonal to the x y-plane of motion. The resulting expression can be written as the sum of the GR and the EC contributions, i.e., (in order to ease the notation, henceforth we write H 11 ≡ H TT 11 ) where d so being the (constant) distance to the astrophysical source.

Numerical estimates
In this section, we deal with binary NS systems and provide some numerical estimates of the contributions introduced by EC model in the flux and the waveform. In Sect. 5.3.1, we set the parameters which are necessary to perform the numerical computations. These are then discussed in Sect. 5.3.2.
Hereafter, a dot signifies a differentiation with respect to the t variable.

Parameter setting
The set of initial conditions characterizing our numerical investigation is represented by (R in , ϕ in ,Ṙ in ,φ in ). The initial radius R in is expressed in terms of R g , with R g ≡ G M/c 2 ; for the initial angle and radial velocity, we assume ϕ in = 0 andṘ in = 0, respectively;φ in is a fraction of the Keplerian velocity, i.e., where for β = 1 the Newtonian eccentricity e 0 vanishes (leading to circular orbits), whereas in the limiting case β → 0 we have e 0 → 1. Therefore, given these premises, the initial conditions are specified once we assign M, R in , and β. A crucial point of our analysis regards the spins of the NSs. These are modeled as follows where, following Ref. [66], n = 10 44 m −3 is estimated as the inverse of the nucleon volume. Therefore, if the masses m 1 and m 2 are known, then the spin components s z1 and s z2 can be immediately calculated.
In order to gain useful information about the binary system's dynamics, we determine the minimum, average, and maximum values of the relative radius (i.e., R min , R av , R max ). Furthermore, to perform some consistency checks, we define a set of parameters, which must be less than 1 due to the hypotheses underlying our model; first of all, the slow-motion condition demands that we compute the maximum values v max 1 /c, v max 2 /c attained by the ratios v 1 /c, v 2 /c, respectively (the velocities v 1 and v 2 of the two bodies can be obtained starting from Eq. (117)); to verify whether the two bodies remain widely separated, we calculate finally, we monitor the strength of the gravitational field through the factor The values of the aforementioned variables, along with other quantities characterizing the binary NS system to be investigated in Sect. 5.3.2, are listed in Table 1.

Discussion of the results
We consider a gravitational system consisting of two NSs, whose parameters can be found in Table 1.
The order of magnitude of the spin components s z1 and s z2 (in units ofh), physically representing the number of neutrons inside the NSs, is consistent with the values reported in the literature (which are of the order of 10 57 neutrons) [107]; moreover, the magnitude of the parameters v max 1 , v max 2 , α 1 , α 2 , γ confirms that the slow-motion, wide-separation, and weakfield hypotheses are fulfilled.
In order to estimate the EC contributions to the GR flux and waveform, we define (cf. Eqs. (129)-(133)) The above quantities, along with the function R(t) representing the relative distance of the NSs, are shown in Fig. 1. From the plot of E F , we see that the spin effects become more significant at the closest point of approach between the objects, where the gravitational field becomes more intense. This agrees with the spirit of EC theory, whose importance is expected to increase in the strong-gravity regime. In our example, the average contributions predicted by EC theory are smaller than GR ones by a factor of 10 −23 . This difference is consistent with the fact that the bodies are widely separated  Table 1). during their dynamical evolution. The EC corrections can be also figured out starting from the trend of E H . As shown in Fig. 1, it goes from its minimum to its maximum values when the NSs get closer. Moreover, this function vanishes as soon as H 11 (t) = 0 (cf. Eqs. (131) and (138b)). These points indicate when GR and EC effects become comparable and in our example we have |H GR 11 | = |H EC 11 | ∼ 10 −52 . The same information, expressed in terms of ϕ, can be inferred from the zeroes of H 11 (ϕ), which occur at ϕ ≈ 1 4 π, 3 4 π, 7 4 π, 11 4 π (see Fig. 2; the functional form of H 11 (ϕ) can be promptly deduced from Eqs. (131)-(133)).

Conclusions
This work configures as a natural continuation of the research program started out in Ref. [66], where we have solved the GW generation problem in EC model at 1PN level by resort- Fig. 2 Gravitational waveform H 11 as a function of the azimuthal angle ϕ. The ϕ-intercepts occurr at ϕ = 0.31π, 0.81π, 1.34π, 1.79π ing to the Blanchet-Damour formalism. This general treatment finds an explicit application here, where the matter source is described by the Weyssenhoff fluid.
The structure of the paper is sketched in Fig. 3. In Sect. 2, we have summarized the key steps of the previous article. In this framework, the spinning PN source is supposed to be a generic hydrodynamical fluid system. For this reason, in Sect. 3.1, we have introduced the fundamental pillars of the EC hydrodynamics verifying the simplifying hypothesis that the torsion tensor has a vanishing trace, i.e., S αβ β = 0. Subsequently, we have modeled the spin effects inside matter by employing the Weyssenhoff model of a semiclassical ideal spinning fluid supplemented by the Frenkel condition (see Sect. 3.2). The study of the Weyssenhoff fluid within the PN approximation scheme, representing a fundamental tool of the Blanchet-Damour formalism, is contained in Sect. 3.3. Both at 0PN and at 1PN level, the dynamics is ruled by a system of partial and integro-differential equations, whose resolution is extremely demanding. A less involved pattern can be obtained if we employ the point-particle procedure, which allows to characterize the fluid dynamics in terms of ordinary differential equations by going from a continuous picture to a discrete description of the system (see Sect. 4). We have then derived, within EC theory and for the particular case of binary systems, the 1PN formula (116) of the center of mass position and the general expressions (118)-(122) of the radiative multipole moments. Starting from these results and the conservation law (105) of the spin vector, we have resorted to the Damour-Deruelle solution in GR (which has been briefly discussed in Sect. 5.1) to set up a hybrid approach for dealing with binaries of spinning PN NSs, where we have provided some numerical estimates of the EC contributions to the flux and the waveform (see Sects. 5.2 and 5.3).
This paper contains some new theoretical results, which can be summarized as follows: Fig. 3 Flowchart of the paper. "Previous paper" refers to Ref. [66] (1) development of a general pattern (subject to the hypothesis S αβ β = 0) for the hydrodynamics in EC theory, where the spin effects are modeled through the tensors Φ αβ and τ μν λ (cf. Eq. (30) The calculation of the parameter n appearing in the spin formula (135) has led to established results for the number of neutrons inside a NS. Moreover, our investigation has revealed that EC contributions become more important when the gravitational field strength grows (see Fig. 1). Therefore, we could in principle extend our hybrid scheme to binary BH systems and evaluate the corrections to the radiated power and the asymptotic amplitude foretold by EC model. If we suppose that, similarly to Eq. (135), the BH spin can be written as (m being the BH mass) then we obtain that the EC effects are, for known astrophysical masses (i.e., 6M M 10 10 M , with M the total mass of the system), between 23 and 13 orders of magnitude lower than GR ones. These differences are justified by the fact that our approach is restricted to compact binaries in their early inspiralling stage. Despite that, our analysis permits to infer that EC corrections can fulfill a relevant role in the later evolution phases. This topic might be useful for testing quantum phenomena in the strong-gravity regime, and it is worth examining it in a separate paper. Furthermore, GW phenomena could lead to interesting implications in the context of generalized EC theories. In particular, Ref. [108] puts forth a model where the gravitational Lagrangian has an additional quadratic-torsion term whose coupling constant differs from the Newtonian gravitational constant G. This new contribution changes the interaction strength between torsion and matter fields and hence novel results can be found also in the GW framework.
In the context of the GR pole-dipole approximation, the lowest-order "classical spin" contributions to the mass-type and current-type radiative moments emerge at 1.5PN and 0.5PN level, respectively [80,110]. Therefore, in EC theory we recover formally the same expressions as in GR, see Eqs.
(118)-(122). This is a consequence of the Frenkel condition (70) and the ensuing PN expansions of the tensors T μν and S μν (see Eqs. (81) and (82)). If we define the Kerr angular momentum as J = Gm 2 c a with a ∈ (0, 1), we can provide an estimate of the effects introduced by EC theory on a single body by evaluating the ratio s/J by means of Eqs. In conclusion, this paper prepares the ground for a systematic study of spinning PN binaries in EC theory. This entails a comprehensive investigation of the 1PN dynamics of binary systems in EC theory, which deserves consideration in a separate paper.