Gaussian gravitational constant
The Gaussian gravitational constant (symbol Template:Mvar) is a parameter used in the orbital mechanics of the Solar System. It relates the orbital period to the orbit's semi-major axis and the mass of the orbiting body in Solar masses.
The value of Template:Mvar historically expresses the mean angular velocity of the system of Earth+Moon and the Sun considered as a two body problem, with a value of about 0.986 degrees per day, or about 0.0172 radians per day. As a consequence of the law of gravitation and Kepler's third law, Template:Mvar is directly proportional to the square root of the standard gravitational parameter of the Sun, and its value in radians per day follows by setting Earth's semi-major axis (the astronomical unit, au) to unity, Template:Mvar:(rad/d) = (Template:MvarTemplate:Solar mass)0.5·au−1.5.
A value of Template:Mvar = Template:Val rad/day was determined by Carl Friedrich Gauss in his 1809 work Theoria Motus Corporum Coelestium in Sectionibus Conicis Solem Ambientum ("Theory of the Motion of the Heavenly Bodies Moving about the Sun in Conic Sections").<ref name="Gauss">Template:Cite book</ref> Gauss's value was introduced as a fixed, defined value by the IAU (adopted in 1938, formally defined in 1964), which detached it from its immediate representation of the (observable) mean angular velocity of the Sun–Earth system. Instead, the astronomical unit now became a measurable quantity slightly different from unity. This was useful in 20th-century celestial mechanics to prevent the constant adaptation of orbital parameters to updated measured values, but it came at the expense of intuitiveness, as the astronomical unit, ostensibly a unit of length, was now dependent on the measurement of the strength of the gravitational force.
The IAU abandoned the defined value of Template:Mvar in 2012 in favour of a defined value of the astronomical unit of Template:Val exactly, while the strength of the gravitational force is now to be expressed in the separate standard gravitational parameter Template:MvarTemplate:Solar mass, measured in SI units of m3⋅s−2.<ref name="Smart53">Template:Cite book</ref>
DiscussionEdit
Gauss's constant is derived from the application of Kepler's third law to the system of Earth+Moon and the Sun considered as a two-body problem, relating the period of revolution (Template:Mvar) to the major semi-axis of the orbit (Template:Mvar) and the total mass of the orbiting bodies (Template:Mvar). Its numerical value was obtained by setting the major semi-axis and the mass of the Sun to unity and measuring the period in mean solar days:
- Template:Mvar = 2Template:Pi Template:Sqrt / (Template:Mvar Template:Sqrt ) ≈ 0.0172021 [rad], where:
- Template:Mvar ≈ 365.256 [days], Template:Mvar = (Template:Solar mass+Template:Earth mass+Template:Lunar mass) ≈ 1.00000304 [[[:Template:Solar mass]]], and Template:Mvar = 1 by definition.
The value represents the mean angular motion of the Earth-Sun system, in radians per day, equivalent to a value just below one degree (the division of the circle into 360 degrees in Babylonian astronomy was likely intended as approximating the number of days in a solar year<ref>David H. Kelley, Eugene F. Milone, Exploring Ancient Skies: A Survey of Ancient and Cultural Astronomy (2011), p. 219</ref>). The correction due to the division by the square root of Template:Mvar reflects the fact that the Earth–Moon system is not orbiting the Sun itself, but the center of mass of the system.
Isaac Newton himself determined a value of this constant which agreed with Gauss's value to six significant digits.<ref>"The numerical value of the Gaussian constant was determined by Newton himself 120 years prior to Gauss. It agrees with the modern value to six significant figures. Hence the name 'Gaussian constant' should be regarded as a tribute to Gauss' services to celestial mechanics as a whole, instead of indicating priority in determining the numerical value of the gravitational constant used in celestial mechanics, as is sometimes considered in referring to his work." Sagitov (1970:713). [This claim is questionable since Sagitov does not give a citation to where Newton computed this value.]</ref> Gauss (1809) gave the value with nine significant digits, as 3548.18761 arc seconds.
Since all involved parameters, the orbital period, the Earth-to-Sun mass ratio, the semi-major axis and the length of the mean solar day, are subject to increasingly refined measurement, the precise value of the constant would have to be revised over time. But since the constant is involved in determining the orbital parameters of all other bodies in the Solar System, it was found to be more convenient to set it to a fixed value, by definition, implying that the value of Template:Mvar would deviate from unity. The fixed value of Template:Mvar = 0.01720209895 [rad] was taken to be the one set by Gauss (converted from degrees to radian), so that Template:Mvar = 4Template:Pi2:(Template:Mvar2 Template:Mvar2 Template:Mvar) ≈ 1.<ref name=Sagitov>Sagitov, M. U., "Current Status of Determinations of the Gravitational Constant and the Mass of the Earth", Soviet Astronomy, Vol. 13 (1970), 712–718, translated from Astronomicheskii Zhurnal Vol. 46, No. 4 (July–August 1969), 907–915.</ref>
Gauss's 1809 value of the constant was thus used as an authoritative reference value for the orbital mechanics of the Solar System for two centuries. From its introduction until 1938 it was considered a measured quantity, and from 1938 until 2012 it was used as a defined quantity, with measurement uncertainty delegated to the value of the astronomical unit. The defined value of Template:Mvar was abandoned by the IAU in 2012, and the use of Template:Mvar was deprecated, to be replaced by a fixed value of the astronomical unit, and the (measured) quantity of the standard gravitational parameter Template:MvarTemplate:Solar mass.
Role as a defining constant of Solar System dynamicsEdit
Gauss himself stated the constant in arc seconds, with nine significant digits, as Template:Mvar = Template:Gaps. In the late 19th century, this value was adopted, and converted to radian, by Simon Newcomb, as Template:Mvar = Template:Gaps.<ref name="Clemence65"> Template:Cite journal</ref> and the constant appears in this form in his Tables of the Sun, published in 1898.<ref> "The adopted value of the Gaussian constant is that of Gauss himself, namely: Template:Mvar = Template:Gaps = Template:Gaps". Template:Cite book</ref>
Newcomb's work was widely accepted as the best then available<ref>Template:Cite journal</ref> and his values of the constants were incorporated into a great quantity of astronomical research. Because of this, it became difficult to separate the constants from the research; new values of the constants would, at least partially, invalidate a large body of work. Hence, after the formation of the International Astronomical Union in 1919 certain constants came to be gradually accepted as "fundamental": defining constants from which all others were derived. In 1938, the VIth General Assembly of the IAU declared,
However, no further effort toward establishing a set of constants was forthcoming until 1950.<ref>Template:Cite journal</ref> An IAU symposium on the system of constants was held in Paris in 1963, partially in response to recent developments in space exploration.<ref name=" Clemence65"/> The attendees finally decided at that time to establish a consistent set of constants. Resolution 1 stated that
Resolution 4 recommended
Included in the list of fundamental constants was
These resolutions were taken up by a working group of the IAU, who in their report recommended two defining constants, one of which was
For the first time, the Gaussian constant's role in the scale of the Solar System was officially recognized. The working group's recommendations were accepted at the XIIth General Assembly of the IAU at Hamburg, Germany in 1964.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>
Definition of the astronomical unitEdit
Gauss intended his constant to be defined using a mean distance<ref group=note>Historically,Template:Citation needed the term mean distance was used interchangeably with the elliptical parameter the semi-major axis. It does not refer to an actual average distance.</ref> of Earth from the Sun of 1 astronomical unit precisely.<ref name=" Clemence65"/> With the acceptance of the 1964 resolutions, the IAU, in effect, did the opposite: defined the constant as fundamental, and the astronomical unit as derived, the other variables in the definition being already fixed: mass (of the Sun), and time (the day of Template:Val seconds). This transferred the uncertainty from the gravitational constant to uncertainty in the semi-major axis of the Earth-Sun system, which was no longer exactly one au (the au being defined as depending on the value of the gravitational constant). The astronomical unit thus became a measured quantity rather than a defined, fixed one.<ref name="Herrick65">Template:Cite journal</ref>
In 1976, the IAU reconfirmed the Gaussian constant's status at the XVIth General Assembly in Grenoble,<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> declaring it to be a defining constant, and that
From this definition, the mean distance of Earth from the Sun works out to 1.000 000 03 au, but with perturbations by the other planets, which do not average to zero over time, the average distance is 1.000 000 20 au.<ref name=" Clemence65"/>
AbandonmentEdit
In 2012, the IAU, as part of a new, self-consistent set of units and numerical standards for use in modern dynamical astronomy, redefined the astronomical unit as<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> Template:Quote
and hence abandoned the Gaussian constant as an indirect definition of scale in the Solar System, recommending
The value of k based on the defined value for the astronomical unit would now be subject to the measurement uncertainty of the standard gravitational parameter, <math>k = \sqrt{G M_\odot } \cdot \text{au}^{-1.5} \cdot \text{d} = {1.32712440018(9)}^{0.5} \cdot 1.495978707^{-1.5} \cdot 8.64 \cdot 10^{-2.5} = 0.0172020989484(6).</math>
Units and dimensionsEdit
Template:Mvar is given as a unit-less fraction of the order of 1.7%, but it can be considered equivalent to the square root of the gravitational constant,<ref>Template:Cite book</ref> in which case it has the units of auTemplate:Frac⋅d−1⋅Template:Solar mass−Template:Frac,<ref name="Clemence65"/> where
- au is the distance for which Template:Mvar takes its value as defined by Gauss—the distance of the unperturbed circular orbit of a hypothetical, massless body whose orbital period is Template:Math days,<ref name="Herrick65"/>
- d is the mean solar day (86,400 seconds),
- Template:Solar mass is the mass of the Sun.
Therefore, the dimensions of Template:Mvar are<ref>Template:Cite book</ref>
- lengthTemplate:Frac time−1 mass−Template:Frac or Template:Math.
In spite of this Template:Mvar is known to much greater accuracy than Template:Mvar (or the square root of Template:Mvar). The absolute value of Template:Mvar is known to an accuracy of about 10−4, but the product Template:Math (the gravitational parameter of the Sun) is known to an accuracy better than 10−10.
DerivationEdit
Gauss's originalEdit
Gauss begins his Theoria Motus by presenting without proof several laws concerning the motion of bodies about the Sun.<ref name="Gauss"/> Later in the text, he mentions that Pierre-Simon Laplace treats these in detail in his Mécanique Céleste.<ref>Template:Cite book</ref> Gauss's final two laws are as follows:
- The area swept by a line joining a body and the Sun divided by the time in which it is swept gives a constant quotient. This is Kepler's second law of planetary motion.
- The square of this quotient is proportional to the parameter (that is, the latus rectum) of the orbit and the sum of the mass of the Sun and the body. This is a modified form of Kepler's third law.
He next defines:
- Template:Math as the parameter (i.e., the latus rectum) of a body's orbit,
- Template:Mvar as the mass of the body, where the mass of the Sun = 1,
- Template:Math as the area swept out by a line joining the Sun and the body,
- Template:Mvar as the time in which this area is swept,
and declares that
- <math>\frac{g}{t\sqrt{p}\sqrt{1+\mu}}</math>
is "constant for all heavenly bodies". He continues, "it is of no importance which body we use for determining this number," and hence uses Earth, defining
- unit distance = Earth's mean distance (that is, its semi-major axis) from the Sun,
- unit time = one solar day.
He states that the area swept out by Earth in its orbit "will evidently be" Template:Math, and uses this to simplify his constant to
- <math>\frac{2\pi}{t\sqrt{1+\mu}}.</math>
Here, he names the constant Template:Mvar and plugging in some measured values, Template:Mvar = Template:Val days, Template:Mvar = Template:Sfrac solar masses, achieves the result Template:Math = Template:Val.
In modern termsEdit
Gauss is notorious for leaving out details, and this derivation is no exception. It is here repeated in modern terms, filling out some of the details.
Define without proof
- <math>h=2\frac{dA}{dt},</math>
where<ref name="Smart77-100">Template:Cite book</ref>
- Template:Math is the time rate of sweep of area by a body in its orbit, a constant according to Kepler's second law, and
- Template:Mvar is the specific angular momentum, one of the constants of two-body motion.
Next define
- <math>h^2=\mu p,</math>
where<ref name="Smart77-101">Smart, W. M. (1977). p. 101.</ref>
- Template:Math, a gravitational parameter,<ref group=note>Do not confuse Template:Mvar the gravitational parameter with Gauss's notation for the mass of the body.</ref> where
- Template:Mvar is Newton's gravitational constant,
- Template:Mvar is the mass of the primary body (i.e., the Sun),
- Template:Mvar is the mass of the secondary body (i.e., a planet), and
- Template:Mvar is the semi-parameter (the semi-latus rectum) of the body's orbit.
Note that every variable in the above equations is a constant for two-body motion. Combining these two definitions,
- <math>\left(2\frac{dA}{dt}\right)^2=G(M+m)p,</math>
which is what Gauss was describing with the last of his laws. Taking the square root,
- <math>2\frac{dA}{dt}=\sqrt{G}\sqrt{M+m}\sqrt{p},</math>
and solving for Template:Math,
- <math>\sqrt{G}=\frac{2dA}{dt\sqrt{M+m}\sqrt{p}}.</math>
At this point, define Template:Math.<ref name="Smart53"/> Let Template:Mvar be the entire area swept out by the body as it orbits, hence Template:Math, the area of an ellipse, where Template:Mvar is the semi-major axis and Template:Mvar is the semi-minor axis. Let Template:Math, the time for the body to complete one orbit. Thus,
- <math>k=\frac{2\pi ab}{P\sqrt{M+m}\sqrt{p}}.</math>
Here, Gauss decides to use Earth to solve for Template:Mvar. From the geometry of an ellipse, Template:Math.<ref name="Smart77-99">Smart, W. M. (1977). p. 99.</ref> By setting Earth's semi-major axis, Template:Math, Template:Mvar reduces to Template:Math and Template:Math. Substituting, the area of the ellipse "is evidently" Template:Math, rather than Template:Math. Putting this into the numerator of the equation for Template:Mvar and reducing,
- <math>k=\frac{2\pi}{P\sqrt{M+m}}.</math>
Note that Gauss, by normalizing the size of the orbit, has eliminated it completely from the equation. Normalizing further, set the mass of the Sun to 1,
- <math>k=\frac{2\pi}{P\sqrt{1+m}},</math>
where now Template:Mvar is in solar masses. What is left are two quantities: Template:Mvar, the period of Earth's orbit or the sidereal year, a quantity known precisely by measurement over centuries, and Template:Mvar, the mass of the Earth–Moon system. Again plugging in the measured values as they were known in Gauss's time, Template:Mvar = Template:Val days, Template:Mvar = Template:Sfrac solar masses,Template:Clarify yielding the result Template:Mvar = Template:Val.
Gauss's constant and Kepler's third lawEdit
The Gaussian constant is closely related to Kepler's third law of planetary motion, and one is easily derived from the other. Beginning with the full definition of Gauss's constant,
- <math>k=\frac{2\pi ab}{P\sqrt{M+m}\sqrt{p}},</math>
where
- Template:Mvar is the semi-major axis of the elliptical orbit,
- Template:Mvar is the semi-minor axis of the elliptical orbit,
- Template:Mvar is the orbital period,
- Template:Mvar is the mass of the primary body,
- Template:Mvar is the mass of the secondary body, and
- Template:Mvar is the semi-latus rectum of the elliptical orbit.
From the geometry of an ellipse, the semi-latus rectum, Template:Mvar can be expressed in terms of Template:Mvar and Template:Mvar thus: Template:Mvar.<ref name="Smart77-99"/> Therefore,
- <math>\sqrt{p}=\frac{b}{\sqrt{a}}.</math>
Substituting and reducing, Gauss's constant becomes
- <math>k=\frac{2\pi}{P}\sqrt{\frac{a^3}{M+m}}.</math>
From orbital mechanics, Template:Math is just Template:Mvar, the mean motion of the body in its orbit.<ref name="Smart77-100"/> Hence,
- <math>\begin{align}
k&=n\sqrt{\frac{a^3}{M+m}},\\[8pt] k^2&=\frac{n^2a^3}{M+m},\\[8pt] k^2(M+m)&=n^2a^3, \end{align}</math>
which is the definition of Kepler's third law.<ref name="Smart77-101"/><ref>Template:Cite book </ref> In this form, it is often seen with Template:Mvar, the Newtonian gravitational constant in place of Template:Math.
Setting Template:Math, Template:Math, Template:Math, and Template:Mvar in radians per day results in Template:Math, also in units of radians per day, about which see the relevant section of the mean motion article.
Other definitionsEdit
Template:See The value of Gauss's constant, exactly as he derived it, had been used since Gauss's time because it was held to be a fundamental constant, as described above. The solar mass, mean solar day and sidereal year with which Gauss defined his constant are all slowly changing in value. If modernTemplate:Clarify values were inserted into the defining equation, a value of Template:Val would result.{{ safesubst:#invoke:Unsubst||date=__DATE__ |$B= Template:Fix }}<ref name="Danby88">Template:Cite book</ref>
It is also possible to set the gravitational constant, the mass of the Sun, and the astronomical unit to 1. This defines a unit of time with which the period of the resulting orbit is equal to Template:Math. These are often called canonical units.<ref name="Danby88"/>
See alsoEdit
- Gravitational constant
- Standard gravitational parameter
- Kepler's laws of planetary motion
- Mean motion