<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2012.329161</article-id><article-id pub-id-type="publisher-id">JMP-23110</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Approximation Method for the Relaxed Covariant Form of the Gravitational Field Equations for Particles
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>manuel</surname><given-names>Gallo</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Osvaldo</surname><given-names>M. Moreschi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>FaMAF, Universidad Nacional de Córdoba, Instituto de Fsica Enrique Gaviola (IFEG), CONICET, Ciudad Universitaria, Córdoba, Argentina</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>emanuelgg@gmail.com(MG)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>28</day><month>09</month><year>2012</year></pub-date><volume>03</volume><issue>09</issue><fpage>1247</fpage><lpage>1254</lpage><history><date date-type="received"><day>July</day>	<month>7,</month>	<year>2012</year></date><date date-type="rev-recd"><day>August</day>	<month>6,</month>	<year>2012</year>	</date><date date-type="accepted"><day>August</day>	<month>14,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  We present a study of the so called relaxed field equations of general relativity in terms of a decomposition of the metric; which is designed to deal with the notion of particles. Several known results are generalized to a coordinate free covariant discussion. We apply our techniques to the study of a particle up to second order.
 
</p></abstract><kwd-group><kwd>General Relativity; Approximation Methods; Particles</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The notion of particle is fundamental to the Newtonian mechanics framework; in fact, this whole theoretical framework can be constructed in terms of the notion of test particles and massive particles. It is then natural to ask whether this notion can be translated to other frameworks, as is general relativity.</p><p>Within general relativity one understands Newtonian mechanics as the limit of weak field and slow motion. So we know that one can regain the notion of particle in this regime. Also in general relativity, the concept of test particle is a natural one, which allows to discuss several physically interesting situations.</p><p>At first sight it is not at all clear that one can extend the notion of particles (non-test) to the realm of general relativity. To begin with, if one imagines a process in which one shrinks the sizes of an object to obtain a point like object, one knows that at some moment in the process one would end up with the formation of a black hole, which has a characteristic size. However, the post-Newtonian approach to compact objects is frequently constructed in terms of the notion of particles; although post-Newtonian systems are normally required to have weak fields and slowly moving objects.</p><p>It is interesting to note that the most simple black hole, namely the one describing a vacuum spherically symmetric spacetime, can be expressed in terms of the so called Kerr-Schild decomposition. In this way, the Schwarzschild black hole, whose maximal analytic extension is described in terms of the well known causal conformal diagrams, when expressed in the Kerr-Schild decomposition shows a point like description in terms of the flat reference metric of the Kerr-Schild form.</p><p>This indicates that it might be possible to give a particle notion to a compact object in general relativity when expressed with respect to background reference metrics.</p><p>If one intends to study the problem of a systems composed of several compact objects, it appears as an appealing strategy to use approximation techniques for solving the field equations. Several problems are related to this.</p><p>In building approximation schemes for the study of the field equations in general relativity it is often useful to recur to the relaxed form of the field equations; that we recall below. Also, it frequently useful to decompose the physical metric in terms of a background metric. In this work we plan to study both techniques.</p><p>In the process of decomposing the metric a key issue is the notion of gauge, since in general one has more than one way to decompose the physical metric. In order to study this issue we bring the techniques used by Friedrich in his study of the hyperbolic nature of the gravitational field equations. We will present here a generalization of Friedrich’s results that is convenient for our discussion.</p><p>Although we work with coordinate independent expressions, we also relate our work with the widely used harmonic gauge condition; and take the opportunity to restate Anderson’s result in a coordinate independent fashion.</p><p>An approximation scheme is suggested in which the previous studies are taking into account.</p><p>We apply our techniques to the problem of a single particle up to the second order.</p></sec><sec id="s2"><title>2. The Decomposition of the Metric</title><p>Let us express the metric <img src="18-7500906\4c8fce1d-41c4-41d3-9682-9f9981e2f8e6.jpg" /> of the spacetime <img src="18-7500906\be61e77e-f31d-44af-b026-8b322e2899ca.jpg" /> in terms of a reference metric<img src="18-7500906\96122a31-0ebe-49a6-b69f-10ee70be3f78.jpg" />, such that</p><disp-formula id="scirp.23110-formula46142"><label>(1)</label><graphic position="anchor" xlink:href="18-7500906\8b02ad53-7565-4ff6-be63-343cdfa0253e.jpg"  xlink:type="simple"/></disp-formula><p>Let <img src="18-7500906\c03fd32c-6640-4aa9-8cb1-29e98af71cfc.jpg" /> denote the torsion free metric connection of <img src="18-7500906\1e18f983-d65a-410f-8d7b-891461c63c34.jpg" /> and <img src="18-7500906\313ede22-b95f-4d91-a1de-19e64c91f975.jpg" /> the torsion free metric connection of<img src="18-7500906\46e42a08-66f1-4059-bb20-29218508d870.jpg" />; then one can express the covariant derivative of an arbitrary vector <img src="18-7500906\b2b92ecc-9d03-4763-af42-0cb7a1f9b335.jpg" /> by</p><disp-formula id="scirp.23110-formula46143"><label>(2)</label><graphic position="anchor" xlink:href="18-7500906\89ff33cf-8db6-4b7f-b448-1cc4c4ea0119.jpg"  xlink:type="simple"/></disp-formula><p>and one can prove that</p><disp-formula id="scirp.23110-formula46144"><label>(3)</label><graphic position="anchor" xlink:href="18-7500906\50024a8b-161b-4d2b-ab73-64c4f267e94d.jpg"  xlink:type="simple"/></disp-formula><p>Let us observe that</p><disp-formula id="scirp.23110-formula46145"><label>(4)</label><graphic position="anchor" xlink:href="18-7500906\1e2a1ba1-e1e2-4343-a5f9-f3a7988ff39f.jpg"  xlink:type="simple"/></disp-formula><p>The relation between <img src="18-7500906\a55f6b71-85c3-40ea-9b54-094ba32025f1.jpg" /> and the curvature tensor can be calculated from</p><disp-formula id="scirp.23110-formula46146"><label>(5)</label><graphic position="anchor" xlink:href="18-7500906\7d273017-4bdb-4577-81f7-1272cf290a14.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="18-7500906\6a6c17a1-b813-4562-8edb-29db27bf8b5a.jpg" /> is the curvature of the <img src="18-7500906\44e9c203-a89d-41fc-b835-698547af4b4c.jpg" /> connection. Then the Ricci tensor can be calculated from</p><disp-formula id="scirp.23110-formula46147"><label>(6)</label><graphic position="anchor" xlink:href="18-7500906\367b29e8-a412-4aba-b501-784e2abdc3c2.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="18-7500906\6a9b0bfc-2214-4297-a8c1-b3c5bce64564.jpg" /> is the Ricci tensor of the connexion<img src="18-7500906\d264425b-2cec-480e-9bd5-9ce586ddc9c1.jpg" />.</p></sec><sec id="s3"><title>3. Auxiliary Functions or Gauge Vector</title><p>Let us consider four independent auxiliary functions<img src="18-7500906\7eb44135-6aa1-4384-92c6-8f715ad13769.jpg" />, with<img src="18-7500906\582ff9dc-aa2e-4db9-994f-814bba43029a.jpg" />. Then let us observe that</p><disp-formula id="scirp.23110-formula46148"><label>(7)</label><graphic position="anchor" xlink:href="18-7500906\0a0a707e-341d-47aa-9452-f9be1075318a.jpg"  xlink:type="simple"/></disp-formula><p>Then, if <img src="18-7500906\84b3311f-d059-4189-99a8-cb59e3c5bc04.jpg" /> denotes the inverse of<img src="18-7500906\15907562-d1f9-4c9e-8653-d8d12e4ca1f4.jpg" />, which exists by assumption of the independence of the set<img src="18-7500906\0044c2b9-4e6a-4b2d-abac-6032b6669fed.jpg" />, one has</p><disp-formula id="scirp.23110-formula46149"><label>(8)</label><graphic position="anchor" xlink:href="18-7500906\313dcda7-0f58-4c82-ae6a-7c8498827fa7.jpg"  xlink:type="simple"/></disp-formula><p>where we are using</p><disp-formula id="scirp.23110-formula46150"><label>(9)</label><graphic position="anchor" xlink:href="18-7500906\614ef6df-c1b6-491a-abb2-c4d793bfe7e8.jpg"  xlink:type="simple"/></disp-formula><p>Alternatively, let us define the gauge vector <img src="18-7500906\1361b991-9487-4f92-9f8c-a1a6efd1ed8e.jpg" /></p><disp-formula id="scirp.23110-formula46151"><label>(10)</label><graphic position="anchor" xlink:href="18-7500906\1d4c8329-6791-479b-8f74-11330b93f60f.jpg"  xlink:type="simple"/></disp-formula><p>which implies</p><disp-formula id="scirp.23110-formula46152"><label>(11)</label><graphic position="anchor" xlink:href="18-7500906\6d8fdcc6-2ffa-43f3-bd02-43ca4d1bfc1a.jpg"  xlink:type="simple"/></disp-formula><p>so that one also has</p><disp-formula id="scirp.23110-formula46153"><label>(12)</label><graphic position="anchor" xlink:href="18-7500906\ea12c1fc-d9fd-4d92-9cc3-f60207b9ed6b.jpg"  xlink:type="simple"/></disp-formula><p>These equations show the relation that exist between working with a coordinate system, given by the set of functions<img src="18-7500906\b3e276c1-66f3-4990-a7a6-e3fbc54d111f.jpg" />, and the gauge vector<img src="18-7500906\ce9162c3-2df8-4c67-b9bc-550aee76c89d.jpg" />; which does not need any reference to coordinate systems at all. In what follows we will try to use the covariant approach that employs the use of the gauge vector<img src="18-7500906\6593c54b-7337-47c0-8184-786193824f99.jpg" />. We emphasize that Latin indices are abstract; and therefore our expressions are coordinate independent and covariant.</p><p>Then, the Ricci tensor can be expressed by</p><disp-formula id="scirp.23110-formula46154"><label>(13)</label><graphic position="anchor" xlink:href="18-7500906\ccc9fb67-1586-473f-9cee-5f220eae7b0f.jpg"  xlink:type="simple"/></disp-formula><p>Let us note that if the vector field <img src="18-7500906\4b7a0adf-c08b-49bc-9331-1edf661f53a5.jpg" /> is given by (12), then for any function <img src="18-7500906\60af249d-3098-4417-ac5c-73d9c39a2164.jpg" /> one has</p><disp-formula id="scirp.23110-formula46155"><label>(14)</label><graphic position="anchor" xlink:href="18-7500906\02f78aae-a8c8-470c-9e4a-bfe8a2145480.jpg"  xlink:type="simple"/></disp-formula><p>In the standard studies on approximations to the solution of the field equations, one frequently finds the choice of harmonic coordinates for the set of the<img src="18-7500906\da772be7-31c5-421f-881a-5a97142630e9.jpg" />’s; however, in Equation (13) one can see that only the vector field <img src="18-7500906\16ff13fa-5932-4c51-813a-58a4b69df578.jpg" /> appears, without any reference to a choice of auxiliary functions. Therefore one could just refer to the gauge vector<img src="18-7500906\6ccacefe-d0ce-4116-8895-0f93e935e2d2.jpg" />.</p></sec><sec id="s4"><title>4. The Field Equations in Relaxed Covariant Form</title><p>Previous to the discussion of the relaxed covariant field equations, we would like to refer to the work of Friedrich [<xref ref-type="bibr" rid="scirp.23110-ref1">1</xref>] and its extension to this coordinate independent discussion.</p><sec id="s4_1"><title>4.1. Friedrich’s Theorem without the Use of Coordinates</title><p>The field equations are</p><disp-formula id="scirp.23110-formula46156"><label>(15)</label><graphic position="anchor" xlink:href="18-7500906\373cfcb6-c800-43be-838b-0c20994b2ae6.jpg"  xlink:type="simple"/></disp-formula><p>Equation (3.22) in reference [<xref ref-type="bibr" rid="scirp.23110-ref1">1</xref>] can be obtained from (13) by expressing it in a coordinate frame and neglecting the <img src="18-7500906\ee060eeb-4437-4c65-8042-8cbf56309349.jpg" /> terms. In this way, one would obtain the analogous expression where all the appearance of <img src="18-7500906\9db94e5b-2d70-4979-b211-343c1a122aa2.jpg" /> derivatives are replaced by coordinate derivatives<img src="18-7500906\6fc1d455-2f9c-4f8e-9e94-1418e672e250.jpg" />, the tensors <img src="18-7500906\18b91980-db5a-4860-b716-47ea00f95490.jpg" /> are replaced by the Chrirstoffel symbols and one uses<img src="18-7500906\548519ae-211a-4e30-9e10-21f5e97e3f64.jpg" />; namely:</p><disp-formula id="scirp.23110-formula46157"><label>(16)</label><graphic position="anchor" xlink:href="18-7500906\40ee149b-d9bb-42f0-8718-c96049ee6b8e.jpg"  xlink:type="simple"/></disp-formula><p>Friedrich has studied [<xref ref-type="bibr" rid="scirp.23110-ref1">1</xref>] this system introducing the notion of “coordinate gauge source’’</p><disp-formula id="scirp.23110-formula46158"><label>(17)</label><graphic position="anchor" xlink:href="18-7500906\058fabf8-1825-4516-9102-1c074b9f0133.jpg"  xlink:type="simple"/></disp-formula><p>Subsequently, Friedrich studied the case in which <img src="18-7500906\5e772316-f7a3-4e2b-80d3-001fa639f008.jpg" /> is given arbitrarily.</p><p>Then, we can rephrase Friedrich’s theorem in the following form:</p><p>Theorem 4.1 Let <img src="18-7500906\e6b09c73-c484-44d7-9848-e97a2fdcae9a.jpg" /> be four independent functions that are used as a coordinate system. If <img src="18-7500906\e293106c-cc42-45f7-ba32-4d725a220f74.jpg" /> is a solution of (16) together with the matter equations such that on the initial surface one has<img src="18-7500906\8f64e6d2-ed9e-4bb7-84c3-2dbdb3d60558.jpg" />, <img src="18-7500906\7a99ef37-b5ea-45e0-b0a4-cd66dd54ab28.jpg" />, then <img src="18-7500906\3e19017d-17bc-408f-918d-aa3882851614.jpg" /> is in fact a solution of Einstein’s field equations.</p><p>&#160;This theorem can be understood in two ways. In one of them, we think that the four coordinates <img src="18-7500906\860cbaf8-4d89-411c-9d24-93e97820c364.jpg" /> are given and then the theorem checks whether the<img src="18-7500906\a52650b7-5622-4110-b2db-9a53f8731c64.jpg" />’s satisfy the above equations. In the other way, one think that the<img src="18-7500906\5139498d-318f-4763-813a-55987f6d316e.jpg" />’s are given and then the theorem checks whether there exists a coordinate system of<img src="18-7500906\31222698-775f-4521-87ac-67ef3aacdcff.jpg" />’s such that the equations in the theorem are satisfied.</p><p>From the fact that<img src="18-7500906\27135328-9f81-47f4-9e2f-38e5a38ec7e6.jpg" />, one deduces, using the same techniques as in [<xref ref-type="bibr" rid="scirp.23110-ref1">1</xref>], that the generalized Friedrich’s theorem holds, namely, consider the four functions <img src="18-7500906\4b5ce672-e2a6-4674-91e5-26c3d3ff0c06.jpg" /> as given a priori, then:</p><p>Theorem 4.2 If <img src="18-7500906\cc190c64-91e4-4bd4-80c8-33de2be43d17.jpg" /> is a solution of (15), with the decomposition of the metric as in (1) and with the Ricci tensor as given by (13) with<img src="18-7500906\e66f5e3d-5352-4019-b30e-4af2292451b1.jpg" />, together with the matter equations such that on the initial surface one has<img src="18-7500906\b865d45b-1573-4819-840c-1d9363fd2ea6.jpg" />,</p><p><img src="18-7500906\16601cc5-ef1c-4183-ab51-170b32370f1c.jpg" />, where <img src="18-7500906\43dd24a1-9ae2-4381-b19f-6ab6f26efc5e.jpg" /> are four independent scalars, then <img src="18-7500906\5e9d30b4-4edf-40bb-ad9c-a055a336af62.jpg" /> is in fact a solution of Einstein’s field equations.</p><p>This result gives great freedom in the problem of finding solutions of the field equations in terms of a reference metric. Suppose that one solves (15) for a given vector field<img src="18-7500906\a4f060d8-eb94-4e42-b77d-1bde3707b99a.jpg" />. Also assume that one can solve for the functions <img src="18-7500906\2e00f645-c749-4040-82f2-ab0dd45013d0.jpg" /> such that<img src="18-7500906\2a2c6e55-38f0-4573-a0ee-017ba101ddd3.jpg" />. Then, let us build a flat metric <img src="18-7500906\da30ef4a-5ef4-45bb-8828-62b4cb8c0255.jpg" /> so that<img src="18-7500906\ec2ba8f3-5ce9-4372-9272-b6e6f16b668f.jpg" />; which in particular can be satisfied if the <img src="18-7500906\ac69e1cd-a898-48b6-a7ed-e1242d63485a.jpg" />s are thought as Cartesian coordinates of<img src="18-7500906\b993e1ce-93a4-412b-95fd-c2794b3857e9.jpg" />. In this way one would obtain<img src="18-7500906\a9e43bd6-34ba-438e-9326-63d707698550.jpg" />, and so have a solution of the field equations.</p><p>It also might be of interest to researchers in numerical relativity, since it provides the possibility to use any coordinate system; i.e., not necessarily an harmonic one.</p><p>Instead, one could have a proposition that does not refer to the auxiliary functions whatsoever; namely</p><p>Theorem 4.3 If <img src="18-7500906\1caa88da-9c71-4052-a915-be7a0eda9a21.jpg" /> is a solution of (15), with the decomposition of the metric as in (1) and with the Ricci tensor as given by (13), together with the matter equations such that on the initial surface one has<img src="18-7500906\39ae4a6e-e3ed-4be1-b1a3-c77f5e483ac3.jpg" />, <img src="18-7500906\5c68b506-4316-4324-a0cc-63bb7d974a6d.jpg" />, then <img src="18-7500906\9fe4f97d-ca96-4117-8edb-dabe198ad222.jpg" /> is in fact a solution of Einstein’s field equations.</p><p>This theorem can be understood in two ways. In one of them, we think that the metric <img src="18-7500906\4d2ccda4-74a1-4c70-ac53-269a0b9ae7a1.jpg" /> is given and then the theorem checks whether the vector <img src="18-7500906\9cd45ae6-b752-4d65-bfa2-714fe82290a7.jpg" /> satisfies the above equations. In the other way, one think that the <img src="18-7500906\c1187c9c-0469-4d04-8179-cb271769024a.jpg" /> is given and then the theorem checks whether there exists a metric <img src="18-7500906\04a543db-5289-414d-b326-f33d4e06c76b.jpg" /> such that the equations in the theorem are satisfied.</p></sec><sec id="s4_2"><title>4.2. Relaxed Covariant Form of the Field Equations and a Generalization of Friedrich’s Theorem</title><p>Alternatively one can use the form of the field equations in terms of a slight different logic.</p><p>When we use the expression of the Ricci tensor as given by (13) in (15), without assuming that <img src="18-7500906\71898a64-5c63-4c1e-8986-c07ced7e9162.jpg" /> is<img src="18-7500906\c541d580-183c-4825-96bb-29e5ff577706.jpg" />, namely</p><disp-formula id="scirp.23110-formula46159"><label>(18)</label><graphic position="anchor" xlink:href="18-7500906\f634c1d1-9a2d-48ce-a3e0-8f794906e43b.jpg"  xlink:type="simple"/></disp-formula><p>we will refer to these as the relaxed field equations [<xref ref-type="bibr" rid="scirp.23110-ref2">2</xref>].</p><p>Using the standard harmonic gauge technique, one would say: solve the relaxed field equation in the coordinate frame, with<img src="18-7500906\63c9a1a5-8845-4f73-ae9d-e7adbbb1e5be.jpg" />, and then require the equation</p><disp-formula id="scirp.23110-formula46160"><label>(19)</label><graphic position="anchor" xlink:href="18-7500906\6ec56f15-b145-4454-aa63-9fb6af8bd9a5.jpg"  xlink:type="simple"/></disp-formula><p>In the standard approach one makes use of coordinate basis; therefore the previous statement would be the complete story. However in our case, <img src="18-7500906\5055a026-7137-4197-b9af-2407cb91b2bf.jpg" />has a second term where two covariant derivatives of <img src="18-7500906\30cc2360-18c1-44eb-a275-7bea4f8c017b.jpg" /> with respect to the metric <img src="18-7500906\c3fa203a-0891-4353-9816-7a03f04a47bb.jpg" /> appears. At this point it is important to notice that if we have the solutions <img src="18-7500906\96d0df33-6719-41d7-9766-c310ee099b53.jpg" /> from (19) then, on constructing <img src="18-7500906\bbfc5d4a-41dc-4888-9a44-1b628020432a.jpg" /> with this as a Cartesian coordinate system, one would obtain<img src="18-7500906\de41c751-880d-4272-ac05-009a0c33260f.jpg" />.</p><p>In some occasions it is preferable to work with a different set of equations. In this respect, several authors have indicated that actually to request Equation (19) is equivalent [2-4] to demand</p><disp-formula id="scirp.23110-formula46161"><label>(20)</label><graphic position="anchor" xlink:href="18-7500906\036021be-4531-44ac-aa88-52bb6c7c9292.jpg"  xlink:type="simple"/></disp-formula><p>When dealing with Einstein equations in the relaxed form, and treating the vacuum case, Equation (20) is understood as the condition that the divergence of the Einstein tensor must be zero (which of course is identically zero in the non relaxed form).</p><p>Let us study the relation between the divergence of the energy-momentum tensor and the vector<img src="18-7500906\99b6adac-89a8-4f10-998a-4ad283931167.jpg" />. One can write the relaxed field Equations (18) in the usual form in which on the right hand side we have just the standard term<img src="18-7500906\4eef32d6-b0d8-45f9-be65-a34e49159c82.jpg" />; and so on the left hand side, the terms involving <img src="18-7500906\c85b431f-f745-47e1-b1ee-f0133333e199.jpg" /> would be</p><disp-formula id="scirp.23110-formula46162"><label>(21)</label><graphic position="anchor" xlink:href="18-7500906\a1de2920-30cc-4f22-8770-deec52ac2e60.jpg"  xlink:type="simple"/></disp-formula><p>where we have used that the term <img src="18-7500906\23146589-a3e3-412c-baff-53ebb9665d9f.jpg" /> contributes with the term</p><disp-formula id="scirp.23110-formula46163"><label>(22)</label><graphic position="anchor" xlink:href="18-7500906\1c160137-2b7c-43e1-9e98-b1b8a1f6e23f.jpg"  xlink:type="simple"/></disp-formula><p>Then in taking its divergence, on the left hand side, the terms involving <img src="18-7500906\a8cd3326-9fba-46dc-af1b-cd183c399c85.jpg" /> are</p><disp-formula id="scirp.23110-formula46164"><label>(23)</label><graphic position="anchor" xlink:href="18-7500906\fb687adb-f98e-41f0-b9b6-fb7f5c83b823.jpg"  xlink:type="simple"/></disp-formula><p>If we replace <img src="18-7500906\1bca43fd-2314-45c7-a376-65964feec9fd.jpg" /> by<img src="18-7500906\15c36c57-b39f-4396-b251-1a5072a50ca0.jpg" />, the divergence of the left hand side would be identically zero, since the Einstein tensor has divergence zero. Therefore we conclude that the divergence of the stress energy-momentum tensor is</p><disp-formula id="scirp.23110-formula46165"><label>(24)</label><graphic position="anchor" xlink:href="18-7500906\c0883fe6-4717-4833-922b-4f912bfb1044.jpg"  xlink:type="simple"/></disp-formula><p>Therefore, the stress energy-momentum is conserved if and only if</p><disp-formula id="scirp.23110-formula46166"><label>(25)</label><graphic position="anchor" xlink:href="18-7500906\cc99f060-d7dd-4a40-8f54-bb7ae44ed654.jpg"  xlink:type="simple"/></disp-formula><p>Working out the relations, one finds that the previous equation can be expressed as</p><disp-formula id="scirp.23110-formula46167"><label>(26)</label><graphic position="anchor" xlink:href="18-7500906\5a758d6e-0c2a-4504-887f-9faa10605ab2.jpg"  xlink:type="simple"/></disp-formula><p>Which coincides with Friedrich calculation.</p><p>It follows that if one solves Equation (26) such that on an initial hypersurface <img src="18-7500906\0a5f08bb-821b-4309-82bd-36ddbc898c14.jpg" /> and<img src="18-7500906\9a87ba3a-1175-48a9-bfc7-f87c537d301f.jpg" />, then the energy-momentum tensor will be conserved in the evolution of the system. Furthermore, one also deduces that:</p><p>Theorem 4.4 If, given the metric<img src="18-7500906\70d92f22-b28f-49e5-b8b7-153b430b906b.jpg" />, one solves the relaxed field equations for <img src="18-7500906\80061bfc-42be-4854-b8ca-551360a7a0c7.jpg" /> together with the matter equations, which include the conservation of the energymomentum tensor, such that <img src="18-7500906\75735780-4f46-47a4-a5ad-451b5ecae094.jpg" /> and <img src="18-7500906\944d4bc2-418d-4efd-82f4-89189ae12ccb.jpg" /> on an initial hypersurface, then <img src="18-7500906\97415152-f202-4cc0-a4de-98b5f2a3d456.jpg" /> is a solution of Einstein equations.</p><p>This is a rephrasing of Friedrich’s theorem applied to a decomposition of the metric and to its general relaxed covariant form of the field equations.</p><p>It is interesting to remark that Anderson [<xref ref-type="bibr" rid="scirp.23110-ref3">3</xref>], using a retarded integral expression for<img src="18-7500906\eaf6f860-301a-46ee-8579-1f9a670fcf7d.jpg" />, was able to prove the equivalence between the conservation of the energymomentum tensor with the harmonic gauge condition. In relation to this let us remark that if the set of functions <img src="18-7500906\907b4000-08f6-4e8e-9b83-9bb037549c80.jpg" /> is obtained from the solutions of (19); and one uses them as harmonic coordinates of the metric<img src="18-7500906\e5b2eef4-c7b7-4684-8230-54ec9b6e2085.jpg" />, then one deduces that<img src="18-7500906\2d1a0570-7bfa-4b48-9506-59d1d471d61c.jpg" />. And also, if<img src="18-7500906\72bf5f37-f6ec-45fd-ae2d-512f28fe650b.jpg" />, then Cartesian coordinates of <img src="18-7500906\d2a9dfa0-ea73-4469-9eeb-aba2b6efa87e.jpg" /> are harmonic coordinates of<img src="18-7500906\9584c320-49f4-45d5-a59b-202ea4566c49.jpg" />. This means that we can state Anderson’s result in a coordinate independent way, namely:</p><p>Theorem 4.5 Let <img src="18-7500906\cfcb9d6b-3536-42fc-8cf0-245ef0abd559.jpg" /> be the retarded solution, with respect to a flat metric <img src="18-7500906\c16ec1d5-d3cf-4608-a994-ffd8f0fc1af5.jpg" /> of the relaxed field equations together with the matter equations of state, such that<img src="18-7500906\97dfbadf-22f9-4b5c-8040-6c02e1fd8e15.jpg" />, then the conservation of the energy-momentum tensor implies that <img src="18-7500906\b8ab4a7a-1b05-419b-aa3a-62f4125d23b2.jpg" /> is a solution of Einstein equations.</p></sec></sec><sec id="s5"><title>5. The Approximation Method and the Treatment of Particles</title><p>The approximation method that we introduce below, is adapted to the treatment of particles; therefore, it is convenient to begin by treating the problem of one single particle in the context of linearized gravity, in order to clarify some of the techniques.</p><sec id="s5_1"><title>5.1. The Gravitational Field from One Particle in Linearized Gravity</title><sec id="s5_1_1"><title>5.1.1. The Description of a Particle</title><p>Let us consider a massive point particle with mass <img src="18-7500906\32520227-d13e-43d8-951d-98f7421472cd.jpg" /> describing, in a flat space-time<img src="18-7500906\110a3bf5-5ec6-465e-a9d7-fe982b8f389c.jpg" />, a curve C which in a Cartesian coordinate system <img src="18-7500906\31c26d0c-ac2e-4c5b-912e-64b86a6be05d.jpg" /> reads</p><disp-formula id="scirp.23110-formula46168"><label>(27)</label><graphic position="anchor" xlink:href="18-7500906\61c35c8c-21da-48db-a7d0-8fa78c51b72b.jpg"  xlink:type="simple"/></disp-formula><p>with <img src="18-7500906\2f53a551-f6de-48ba-aa50-bbf2cc0dbd99.jpg" /> meaning the proper time of the particle along<img src="18-7500906\c836cc97-7e98-4ae4-9b80-a47012eacef5.jpg" />.</p><p>The unit tangent vector to<img src="18-7500906\7ce04c61-b927-4bf2-b2bb-cebf11b34cf6.jpg" />, with respect to the flat background metric is</p><disp-formula id="scirp.23110-formula46169"><label>(28)</label><graphic position="anchor" xlink:href="18-7500906\04e405cb-8597-42ee-8291-db67a69b0f43.jpg"  xlink:type="simple"/></disp-formula><p>that is,<img src="18-7500906\2b8ac31c-ffcc-431c-97e1-b965d499e56a.jpg" />. Now, for each point <img src="18-7500906\32b2a7ea-d77a-411d-b385-19c9287674d6.jpg" /> of<img src="18-7500906\7cc31403-c015-4587-af01-1e334add8538.jpg" />, we draw a future null cone <img src="18-7500906\0d9af050-b1a2-4da8-9717-218470a0095d.jpg" /> with vertex in<img src="18-7500906\91c0ddf1-0b31-46d1-9f2b-de3603587418.jpg" />. If we call <img src="18-7500906\d8f679a6-b383-4a84-b5ce-d04f02cd7452.jpg" /> the Minkowskian coordinates of an arbitrary point on the cone<img src="18-7500906\846aa866-efa3-4d03-a7fe-47929f4bfffb.jpg" />, then we can define the retarded radial distance on the null cone by</p><disp-formula id="scirp.23110-formula46170"><label>(29)</label><graphic position="anchor" xlink:href="18-7500906\c7873491-8345-4c58-a926-e4f903551243.jpg"  xlink:type="simple"/></disp-formula><p>The energy momentum tensor <img src="18-7500906\1ac0c6dc-f6af-484b-a625-b2e2fefa384d.jpg" /> of a point particle is proportional to<img src="18-7500906\5d02ed0b-a37e-4d8b-8ae8-f06f8bdf77f2.jpg" />; where <img src="18-7500906\c148e297-8167-4a94-9d7a-a7a4624a16cb.jpg" /> is the mass and <img src="18-7500906\c02f2e79-47f6-4946-9d52-49a1b17a3a94.jpg" /> its four velocity. We are distinguishing between the unit tangent vector <img src="18-7500906\f27f6689-7dab-4da1-bda2-56046467be20.jpg" /> and the four velocity vector<img src="18-7500906\fc195fc3-926a-478e-b5f2-ecfdef38016c.jpg" />, because in future works we would like to consider the possibility to normalize the vector <img src="18-7500906\ca8037ab-87f8-4614-bdf9-46335d6ed817.jpg" /> with respect to a different metric. In order to represent a point particle <img src="18-7500906\06413bf5-d383-47bb-a6cb-6111e5570004.jpg" /> must also be proportional to a three dimensional delta function that has support on the world line of the particle.</p><p>We will suppose that the particle does not have multipolar structure. Then, given an arbitrary Minkowskian frame<img src="18-7500906\4d9088dd-6049-4faf-b1fa-f5f8b9bf1ddd.jpg" />, we will express the energy momentum by</p><disp-formula id="scirp.23110-formula46171"><label>(30)</label><graphic position="anchor" xlink:href="18-7500906\d5aa34f8-62be-4921-9780-edc7f39d830d.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s5_1_2"><title>5.1.2. The First Order Solution</title><p>The retarded solution, in terms of Green functions, for the relaxed field Equations (18) for particle<img src="18-7500906\86b94d89-7a53-4d96-9f52-75d9f761fae6.jpg" />, in which we take <img src="18-7500906\a09d7646-a97c-4635-858a-9e05c397f92f.jpg" /> and <img src="18-7500906\9c11a12b-d86e-4128-8ccc-b45ed273d1b1.jpg" /> a flat metric, is</p><disp-formula id="scirp.23110-formula46172"><label>(31)</label><graphic position="anchor" xlink:href="18-7500906\e233baf4-9484-4f07-86ac-54457879c2a5.jpg"  xlink:type="simple"/></disp-formula><p>so that in general</p><disp-formula id="scirp.23110-formula46173"><label>(32)</label><graphic position="anchor" xlink:href="18-7500906\6b96fa83-8ff0-4cf4-bf78-a5c2423db979.jpg"  xlink:type="simple"/></disp-formula><p>In these equations we have considered the definition</p><disp-formula id="scirp.23110-formula46174"><label>(33)</label><graphic position="anchor" xlink:href="18-7500906\ec427925-4eb3-4740-89aa-9fe9ac79fbcd.jpg"  xlink:type="simple"/></disp-formula><p>It is interesting to realize that the exact inverse of this metric is</p><disp-formula id="scirp.23110-formula46175"><label>(34)</label><graphic position="anchor" xlink:href="18-7500906\30b1806f-1dbd-46b2-a585-c8f5b6f5e99c.jpg"  xlink:type="simple"/></disp-formula><p>Note that one can solve for <img src="18-7500906\c0a9b26b-3996-464c-9865-bb39daa2783c.jpg" /> for an arbitrary motion of the particle; however, the complete solution of the problem involves having to set also<img src="18-7500906\bb2828ee-48b9-4169-b111-50445c8f1278.jpg" />; which in terms of a coordinate frame treatment is equivalent to the harmonic condition. Then, recalling, as mentioned previously, that Anderson has proved [<xref ref-type="bibr" rid="scirp.23110-ref3">3</xref>] the equivalence between the harmonic condition and the divergence free condition on the energy-momentum tensor; one deduces from this, that for the case of the energy-momentum tensor of a particle it implies its geodesic motion.</p></sec></sec><sec id="s5_2"><title>5.2. Iterative Approximation Method</title><p>Now we present a general iterative method to solve the relaxed field equations.</p><p>First of all, let us note that given the decomposition (1) and defining the tensor <img src="18-7500906\733d0cc4-7b71-42e3-b5e8-40952e6977ce.jpg" /> from</p><disp-formula id="scirp.23110-formula46176"><label>(35)</label><graphic position="anchor" xlink:href="18-7500906\5099a8d5-98f7-4044-9c1f-9d14ef4c47c4.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.23110-formula46177"><label>(36)</label><graphic position="anchor" xlink:href="18-7500906\565f0aff-85e9-43f9-b9be-bb0b4f6dd5ca.jpg"  xlink:type="simple"/></disp-formula><p>that is <img src="18-7500906\77c191ca-315d-41ef-8db5-378b42916c4d.jpg" /> is the inverse of<img src="18-7500906\c6214c9b-988c-40a3-b3ef-f6235017e52d.jpg" />, one can always express the inverse <img src="18-7500906\d1ca009a-c85c-4df7-83e3-9ec0b930aa16.jpg" /> in the form</p><disp-formula id="scirp.23110-formula46178"><label>(37)</label><graphic position="anchor" xlink:href="18-7500906\8006d75b-7fc4-4ceb-913d-5edaa572e210.jpg"  xlink:type="simple"/></disp-formula><p>Then making the contraction</p><disp-formula id="scirp.23110-formula46179"><label>(38)</label><graphic position="anchor" xlink:href="18-7500906\503c2809-7e6a-4b48-88f6-63366aac83d6.jpg"  xlink:type="simple"/></disp-formula><p>one finds</p><disp-formula id="scirp.23110-formula46180"><label>(39)</label><graphic position="anchor" xlink:href="18-7500906\2a7fa432-ceb8-4705-9179-7510420c8ae1.jpg"  xlink:type="simple"/></disp-formula><p>which can be considered an implicit equation for<img src="18-7500906\ffed443a-7739-4324-bc10-f1e68e46996f.jpg" />; but it also shows explicitly that <img src="18-7500906\d6a41653-6739-4d3c-98cf-9a8fc39b2514.jpg" /> is quadratic in terms of<img src="18-7500906\5796d72d-557c-4078-bfa8-56970fa2d5a8.jpg" />.</p><p>This suggests the natural series <img src="18-7500906\d4b27eb6-d377-49b9-92b3-3b05a1c25807.jpg" /> defined by</p><disp-formula id="scirp.23110-formula46181"><label>(40)</label><graphic position="anchor" xlink:href="18-7500906\0009fc23-64d6-4f59-a19a-23c9f49fa325.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46182"><label>(41)</label><graphic position="anchor" xlink:href="18-7500906\fd2bd9d8-37e8-4080-ad14-fbbe86c0eccd.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46183"><label>(42)</label><graphic position="anchor" xlink:href="18-7500906\567c59d9-836f-4ebd-8c17-456e63d882cf.jpg"  xlink:type="simple"/></disp-formula><p>for natural numbers<img src="18-7500906\59e57d8d-6594-45a6-ab78-75c00146be5e.jpg" />. It is clear that <img src="18-7500906\0cc0a36e-a82f-4883-8a47-1649fee58001.jpg" /> is order<img src="18-7500906\390dbd43-7897-40f4-9236-9fdaf0b3c88d.jpg" />.</p><p>However, we have seen that in the first order solution for a single particle, the inverse of the metric has a term which is conformal to the flat metric<img src="18-7500906\0de64dd5-bfbc-4b88-a5ca-fd76ae1a7aea.jpg" />; which it will be convenient to take into account. For this reason we propose the following method of approximation where this issue is considered.</p><p>The idea is to express (18) and eventually (19) in the form</p><disp-formula id="scirp.23110-formula46184"><label>(43)</label><graphic position="anchor" xlink:href="18-7500906\315826d9-ad5c-48fd-a29c-220cc4405a15.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="18-7500906\cb89758a-8f9a-4245-a261-44f13f546095.jpg" /> is the term proportional to <img src="18-7500906\7b7d7983-e0de-4719-bad8-f04eb8d5db16.jpg" /> that is contained in<img src="18-7500906\2877fec3-8eab-451c-9cee-804725681020.jpg" />; while the general case would be to consider just <img src="18-7500906\62acc9b7-7810-49c0-a3a9-a95bf20141e3.jpg" /> for the left hand side. This equation can also be expressed by</p><disp-formula id="scirp.23110-formula46185"><label>(44)</label><graphic position="anchor" xlink:href="18-7500906\3a69a66b-72dd-48d0-baa2-2e9f459704ab.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.23110-formula46186"><label>(45)</label><graphic position="anchor" xlink:href="18-7500906\74a0a45d-f071-4feb-8a9e-677e90c0e6e6.jpg"  xlink:type="simple"/></disp-formula><p>Now one would like to solve Equation (44) by iterations.</p><p>Let us define the sets <img src="18-7500906\86252be3-f492-4880-8df9-6ecd146f75a3.jpg" /> such that for<img src="18-7500906\c354dd30-436c-4179-821f-06c4acef0df1.jpg" />, one takes<img src="18-7500906\ff4d7e88-7ca8-4e64-8400-e85fffa88352.jpg" />, <img src="18-7500906\361bcf0c-5ad1-4852-94a8-fce9d1efc2ae.jpg" />to be harmonic functions of the metric <img src="18-7500906\5dcea4e2-c8e3-4a80-a5ae-743f58bd17a5.jpg" /> and<img src="18-7500906\5b19f341-ef78-46ab-adf3-c0c7107bf33e.jpg" />; and for<img src="18-7500906\c1643970-d05d-4bf1-b281-592edc1d1467.jpg" />, <img src="18-7500906\b6248f35-ed24-4b16-878a-63b225648b62.jpg" />is the solution of</p><disp-formula id="scirp.23110-formula46187"><label>(46)</label><graphic position="anchor" xlink:href="18-7500906\22217632-8d75-463d-a9c4-d339894023ac.jpg"  xlink:type="simple"/></disp-formula><p>using the retarded Green function. As we have seen above, <img src="18-7500906\66eebe28-d7f4-4d0e-b06d-7698756b6e95.jpg" />clearly arises in the first order calculation.</p><p>The application of this method to the first order, for a single particle, reproduces the calculation explained in Subsection 5.1.2. Next we study this case at second order.</p></sec><sec id="s5_3"><title>5.3. The Second Order Solution</title><p>Let us remark that the first order solution is stationary and spherically symmetric. This structure transports to the second order solution.</p><p>The equation for <img src="18-7500906\b919f7ee-8f61-4c1b-8fd9-669602b815c0.jpg" /> is</p><disp-formula id="scirp.23110-formula46188"><label>(47)</label><graphic position="anchor" xlink:href="18-7500906\43cb2d9e-e5e8-48dc-847c-5d570fd398c5.jpg"  xlink:type="simple"/></disp-formula><p>We will call the right hand side, the tensor<img src="18-7500906\315b9076-ee11-4c26-9657-7db286e5cde5.jpg" />; which has the structure</p><disp-formula id="scirp.23110-formula46189"><label>(48)</label><graphic position="anchor" xlink:href="18-7500906\6ff350fc-7a6b-4475-abc6-e0b1fdbed98d.jpg"  xlink:type="simple"/></disp-formula><p>where we are using the three dimentional notation <img src="18-7500906\ed19a894-0cc9-4e89-9844-42bf6bc07394.jpg" /> and where</p><disp-formula id="scirp.23110-formula46190"><label>(49)</label><graphic position="anchor" xlink:href="18-7500906\dd106b48-7525-40d4-bd05-1ab37fbd974d.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46191"><label>(50)</label><graphic position="anchor" xlink:href="18-7500906\fa09fde7-abcd-44c9-b698-2a6a7f914a19.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46192"><label>(51)</label><graphic position="anchor" xlink:href="18-7500906\d66b4ed8-734f-41bd-abe7-cfb7af6b458a.jpg"  xlink:type="simple"/></disp-formula><p>Therefore one assumes for <img src="18-7500906\c5203c52-decf-44a2-b114-f95e391a2780.jpg" /> the same form, namely</p><disp-formula id="scirp.23110-formula46193"><label>(52)</label><graphic position="anchor" xlink:href="18-7500906\5b97328f-8b49-4394-806b-cc7188f4b855.jpg"  xlink:type="simple"/></disp-formula><p>In this way one has</p><disp-formula id="scirp.23110-formula46194"><label>(53)</label><graphic position="anchor" xlink:href="18-7500906\0388bcc9-e217-4d3f-ba2e-8da8121ec7f1.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46195"><label>(54)</label><graphic position="anchor" xlink:href="18-7500906\a7756332-09a5-4544-9a4d-d2ed72cde844.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46196"><label>(55)</label><graphic position="anchor" xlink:href="18-7500906\bca915db-49dd-4c32-ab3e-e2eb7ae12080.jpg"  xlink:type="simple"/></disp-formula><p>where the index <img src="18-7500906\d67daa12-7f60-495f-949f-a1adfe8a7ff9.jpg" /> denote spatial coordinates.</p><p>One can see then that the equations to solve are</p><disp-formula id="scirp.23110-formula46197"><label>(56)</label><graphic position="anchor" xlink:href="18-7500906\912b4273-55f7-4a51-864d-3f578979558a.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46198"><label>(57)</label><graphic position="anchor" xlink:href="18-7500906\e3134fc1-19f1-4c88-8a2d-942fb225d91b.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46199"><label>(58)</label><graphic position="anchor" xlink:href="18-7500906\7e2439fe-c0b6-4429-8ad4-75b25b009c87.jpg"  xlink:type="simple"/></disp-formula><p>where we are using the symbol <img src="18-7500906\d98b9dc7-f173-4a1e-8bc2-b5ee2f2aef09.jpg" /> to denote<img src="18-7500906\241754aa-d691-400d-ad76-660602a4309c.jpg" />.</p><p>We can solve these equations in two ways, either using Green function techniques, or, recalling the stationary nature of the solution, just integrating the Laplace operator. For this presentation we choose the second option. Let us note that for any function <img src="18-7500906\593445cf-4c96-4ea4-b7b1-8c31a95565dc.jpg" /> one has that</p><disp-formula id="scirp.23110-formula46200"><label>(59)</label><graphic position="anchor" xlink:href="18-7500906\14bbc893-6538-4da8-98a4-4a2e741e5ce2.jpg"  xlink:type="simple"/></disp-formula><p>Therefore one can find <img src="18-7500906\6ef0afce-2ee0-403e-afaf-bae6cc04c0af.jpg" /> by two consecutive integrations, obtaining</p><disp-formula id="scirp.23110-formula46201"><label>(60)</label><graphic position="anchor" xlink:href="18-7500906\cb5779d8-03a7-4465-9464-9930be93f2cd.jpg"  xlink:type="simple"/></disp-formula><p>Similarly one can see that the function <img src="18-7500906\f6817381-d076-4a38-a17c-1a7e13e33af1.jpg" /> satisfies</p><disp-formula id="scirp.23110-formula46202"><label>(61)</label><graphic position="anchor" xlink:href="18-7500906\bb3a54bc-576d-4065-b18f-a26ad1c49e31.jpg"  xlink:type="simple"/></disp-formula><p>which after integration gives</p><disp-formula id="scirp.23110-formula46203"><label>(62)</label><graphic position="anchor" xlink:href="18-7500906\3f1ad785-4602-4f13-ab2a-cfa6f1453bce.jpg"  xlink:type="simple"/></disp-formula><p>Then the function <img src="18-7500906\5f078625-09d3-4aea-9d64-d4a9ed76799f.jpg" /> is given by</p><disp-formula id="scirp.23110-formula46204"><label>(63)</label><graphic position="anchor" xlink:href="18-7500906\7806d6af-1dae-4f04-a19c-da8f0e4ea3e8.jpg"  xlink:type="simple"/></disp-formula><p>Our choice for the integration constants is:<img src="18-7500906\6308c563-d8ad-4152-b57b-465c7f5f0adf.jpg" />,</p><p><img src="18-7500906\3792ed78-5cdd-4baf-a647-30fdf48fb629.jpg" />, <img src="18-7500906\1d2eecd8-862f-4e3d-9d52-93fe16d3f037.jpg" />, <img src="18-7500906\196e53bd-41df-47c4-b81d-dd23b7e7528a.jpg" />, <img src="18-7500906\a0d2bde4-750b-4d38-a0a1-6a8e788098c6.jpg" />and</p><p><img src="18-7500906\7ce53268-dbe7-486d-b220-950eb133fd9c.jpg" />. This choice is made taking into consideration the exact solution and the integration of the solution coming from a Green function approach; that we will not discuss here due to considerations of space. By the exact functions we mean the metric components of the Schwarzschild spacetime in harmonic coordinates; which are</p><disp-formula id="scirp.23110-formula46205"><label>(64)</label><graphic position="anchor" xlink:href="18-7500906\75ec092f-d80d-4aff-94f3-43476a2ecc7f.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46206"><label>(65)</label><graphic position="anchor" xlink:href="18-7500906\b3520435-d845-485b-baa1-fefb2f693977.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23110-formula46207"><label>(66)</label><graphic position="anchor" xlink:href="18-7500906\8ac448f8-4bf0-45f1-88fe-eb2d598529c4.jpg"  xlink:type="simple"/></disp-formula><p>A graphical comparison with the exact functions of the Schwarzschild solution in harmonic coordinates are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>It can be observed that in second order one obtains anexcellent comparison of the solution with the exact values of the metric components; even for very small values of the radial coordinate. Although this comparison has limited value, it is in any case remarkable that it is only necessary to go only to second order to obtain such a good approximation.</p></sec></sec><sec id="s6"><title>6. Final Comments</title><p>We have presented an study of an approach to the gravitational field equation through the relaxed covariant form of them. The whole approach is intended to deal with the notion of compact objects.</p><p>The relaxed field equations was studied using Friedrich approach to the problem and we have also refer to Anderson’s result in the field of harmonic conditions.</p><p>We have generalized Friedrich results to a covariant formulation in terms of a decomposition of the metric.</p><p>Anderson’s result has been restated in a form that does not make reference to coordinate conditions.</p><p>We have presented an approximation method that can be applied to the notion of particles in general relativity; and which is successful in second order for the case of a solitary compact body.</p><p>It is our intention to apply these techniques to the problem of a binary system in general relativity.</p></sec><sec id="s7"><title>7. Acknowledgements</title><p>We acknowledge support from CONICET and SeCyTUNC.</p></sec><sec id="s8"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.23110-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">H. Friedrich, “On the Hyperbolicity of Einstein’s and Other Gauge Field Equations,” Communications in Mathematical Physics, Vol. 100, No. 4, 1985, pp. 525-543.  
Hdoi:10.1007/BF01217728</mixed-citation></ref><ref id="scirp.23110-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">M. Walker and C. M. Will, “The Approximation of Radiative Effects in Relativistic Gravity-Gravitational Radiation Reaction and Energy Loss in Nearly Newtonian Systems,” Astrophysical Journal, Vol. 242, 1980, pp. L129-L133. Hdoi:10.1086/183417</mixed-citation></ref><ref id="scirp.23110-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">J. L. Anderson, “Satisfaction of deDonder and Trautman Conditions by Radiative Solutions of the Einstein Field Equations,” General Relativity and Gravitation, Vol. 4, No. 4, 1973, pp. 289-297. Hdoi:10.1007/BF00759848</mixed-citation></ref><ref id="scirp.23110-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">A. Einstein, L. Infeld and B. Hoffmann, “The Gravitational Equations and the Problem of Motion,” Annals of Mathematics, Vol. 39, No. 1, 1938, pp. 65-100. 
Hdoi:10.2307/1968714</mixed-citation></ref></ref-list></back></article>