Orbits of the Galilean Moons: Difference between revisions

From PC5214 wiki
Jump to navigation Jump to search
Matthew (talk | contribs)
Added image of Jupiter live view
Matthew (talk | contribs)
moved image of moon and jupiter to conclusion
 
(114 intermediate revisions by the same user not shown)
Line 1: Line 1:
[[File:MWGalileanSatellites.jpg|right|thumb|300px|Jupiter and the four Galilean moons]]
Galileo Galilei’s discovery of celestial bodies that orbit something other than the Earth marked the beginning of the end of the geocentric model of the universe. In this project, we will perform the same observations on those moons as Galileo did 400 years ago.
Galileo Galilei’s discovery of celestial bodies that orbit something other than the Earth marked the beginning of the end of the geocentric model of the universe. In this project, we will perform the same observations on those moons as Galileo did 400 years ago.


==Team Members==
==Team Members==
* A0168128Y [[User:Matthew|Matthew Wee]]
* A0168128Y [[User:Matthew|Matthew Wee]]
==Objectives==
* To observe the planet [https://en.wikipedia.org/wiki/Jupiter Jupiter] and its four largest moons—[https://en.wikipedia.org/wiki/Io_(moon) Io], [https://en.wikipedia.org/wiki/Europa_(moon) Europa], [https://en.wikipedia.org/wiki/Ganymede_(moon) Ganymede], and [https://en.wikipedia.org/wiki/Callisto_(moon) Callisto]
* To measure the orbital parameters of those four moons, such as orbital [https://en.wikipedia.org/wiki/Orbital_period period] and [https://en.wikipedia.org/wiki/Semi-major_axis radius]
* To explore the [https://en.wikipedia.org/wiki/Orbital_resonance orbital resonances] of the inner three of those moons—Io, Europa, and Ganymede
* To verify [https://en.wikipedia.org/wiki/Kepler%27s_laws_of_planetary_motion#Third_law Kepler’s Third Law]
* To assess the viability of consumer equipment in performing precise scientific measurements
==History==
[[File:MWSidereusNunciusAnimation.gif|right|thumb|Animation of Galileo's observations as drawn in his book [https://en.wikipedia.org/wiki/Sidereus_Nuncius ''Sidereus Nuncius''] in 1610]]
'''Jupiter''' is a planet that was known to the ancients as it is the fourth brightest object in the sky behind the Sun, the Moon, and Venus, and is visible to the naked eye. In 1609, astronomer [https://en.wikipedia.org/wiki/Galileo_Galilei Galileo Galilei] made improvements to his telescope, which was the best in the world at the time. With the telescope, over the course of several weeks, he observed four “stars” moving in a line around Jupiter and was persuaded after just four days of observations that they were not stationary stars but instead objects orbiting the planet Jupiter. German astronomer [https://en.wikipedia.org/wiki/Simon_Marius Simon Marius] also made this discovery independently and at the same time, but he is lesser known because he published his observations only after Galileo did. He had the privilege, however, of having his names of these moons (Io, Europa, Ganymede, Callisto) be widely accepted, much to Galileo's dismay.
The discovery that there exist objects that do not orbit the Earth, along with several other breakthroughs, was the cornerstone of the [https://en.wikipedia.org/wiki/Copernican_Revolution Copernican Revolution], which was the transition from the geocentric to the heliocentric model, with the Sun instead of the Earth at the center of the universe.
Since then, astronomers have discovered a total of [https://en.wikipedia.org/wiki/Moons_of_Jupiter 80 moons that orbit Jupiter], 23 of which have no official names.
===Background===
The four Galilean moons, in increasing order of orbital radius, are '''Io''', '''Europa''', '''Ganymede''', and '''Callisto''', of which, Io, Callisto, and Ganymede are larger than the Earth's [https://en.wikipedia.org/wiki/Moon Moon] in both mass and diameter. Ganymede is even larger than the planet [https://en.wikipedia.org/wiki/Mercury_(planet) Mercury]. Because of their size, they can be observed from Earth with a small telescope or even a pair of powerful binoculars.
The main challenge faced when trying to observe these moons is the sheer brightness of Jupiter relative to the dim moons, and their small angular separation from Jupiter. This is especially true for the innermost moons Io and Europa.


==Apparatus==
==Apparatus==
[[File:JupiterScreenshot.jpg|right|thumb|300px|Screenshot of a camera live view of Jupiter]]
[[File:MWCameraSetup.jpg|right|thumb|300px|Setup used to image Jupiter and its satellites]]
No specialized equipment was used in this experiment to make observations of the planet and its moons.
 
* Canon EOS 60D DSLR<ref name="canoncamera">Canon Camera Museum (2010) [https://global.canon/en/c-museum/product/dslr805.html ''Interchangeable Lens Digital Cameras - Digital SLR Camera: EOS 60D Specifications'']</ref>
** APS-C sensor (1.6× crop)
** Pixel size 4.31µm
* Sigma 150–600mm f/5–6.3<ref name="sigmalens">Sigma Corporation (2018) [https://www.sigma-global.com/en/lenses/c015_150_600_5_63/ 150-600mm F5-6.3 DG OS HSM Contemporary Lens Specifications]</ref>
** Lens diameter 95mm
** Field of view: 16.4º–4.1º (effective 10.2º–2.6º)
* Low-end consumer tripod
** Spirit level
* Remote shutter
 
===Exposure Settings===
* Manual
* Focal length: 600mm (effective 960mm)
* 1/15s exposure
* f/6.3 aperture
* ISO 3200 sensitivity
* Manual focus
* Image stabilization off
 
The exposure time cannot be set for too long as motion blur due to the Earth's rotation will be visible in the image.
 
==Data==
[[File:MWJupiterScreenshot.jpg|right|thumb|300px|Screenshot of a camera live view of Jupiter]]
[[File:MWSampleExposure.jpg|right|thumb|300px|An actual image captured on 1 Feb (click to zoom in)]]
===Procedure===
A series of observations were performed with the telephoto lens mounted on the camera. Weather permitting, exposures were taken every day within an hour of sunset for two and a half calendar weeks from 31 January 2022 to 17 February 2022.
 
At the time of these observations, Earth's motion around the Sun (and to a lesser extent, Jupiter's as well) caused the apparent position of Jupiter in the sky to drift closer to the sun every day, forcing observations to occur closer and closer to sunset. The last observation on 17 February was taken at 7:59 PM, when the sun had barely set and the sky was still relatively bright (compared to the dim moons). It was decided that no more observations could be taken thereafter.
 
The moons in the captured images were identified using [https://stellarium-web.org Stellarium] or [https://icandiapps.com/night-sky-4/ Night Sky].
 
By inspecting the moons and their positions relative to Jupiter in the captured images, their apparent angular separation can be determined for every observation. The data are presented in the table below.
 
===Raw Data===
 
{| class="wikitable" style="text-align: center
|-
! scope="col" style="width: 20%" rowspan="2"| Observation
! scope="col" colspan="4"| Angular distance between moon and center of Jupiter (arbitrary units)
|-
! scope="col" style="width: 20%"| Io
! scope="col" style="width: 20%"| Europa
! scope="col" style="width: 20%"| Ganymede
! scope="col" style="width: 20%"| Callisto
|-
! scope="row"| 2022-01-31
| 68.9
| style="font-style: italic; opacity: 0.5;"| occultation
| 102.7
| -254.2
|-
! scope="row"| 2022-02-01
| -32.2
| 155.4
| 245.7
| -371.0
|-
! scope="row"| 2022-02-02
| style="font-style: italic; opacity: 0.5;" colspan="4" | no data collected
|-
! scope="row"| 2022-02-03
| 50.6
| -127.9
| 23.8
| -437.9
|-
! scope="row"| 2022-02-04
| -76.3
| 113.8
| -175.9
| -392.3
|-
! scope="row"| 2022-02-05
| 96.4
| 83.2
| -251.8
| -284.4
|-
! scope="row"| 2022-02-06
| style="font-style: italic; opacity: 0.5;" colspan="4" | no data collected
|-
! scope="row"| 2022-02-07
| 86.2
| -31.0
| 62.2
| 25.5
|-
! scope="row"| 2022-02-08
| -55.7
| 154.5
| 224.8
| 183.0
|-
! scope="row"| 2022-02-09
| style="font-style: italic; opacity: 0.5;" colspan="4" | no data collected
|-
! scope="row"| 2022-02-10
| style="font-style: italic; opacity: 0.5;"| occultation
| -145.8
| 68.0
| 405.6
|-
! scope="row"| 2022-02-11
| -52.4
| 87.2
| -140.1
| 437.5
|-
! scope="row"| 2022-02-12
| 81.0
| 108.1
| -246.5
| 408.5
|-
! scope="row"| 2022-02-13
| style="font-style: italic; opacity: 0.5;" colspan="4" | no data collected
|-
! scope="row"| 2022-02-14
| style="font-style: italic; opacity: 0.5;" colspan="4" | no data collected
|-
! scope="row"| 2022-02-15
| -75.2
| 145.8
| 195.9
| 40.0
|-
! scope="row"| 2022-02-16
| style="font-style: italic; opacity: 0.5;" colspan="4" | no data collected
|-
! scope="row"| 2022-02-17
| style="font-style: italic; opacity: 0.5;"| transit
| -147.2
| 101.8
| -252.2
|-
|}
 
When a moon is [https://en.wikipedia.org/wiki/Occultation '''occulted'''] by Jupiter, it is not visible and its position cannot be determined with precision. Similarly, during a [https://en.wikipedia.org/wiki/Transit_(astronomy) '''transit'''], Jupiter outshines the moon and its precise position also cannot be determined. Either way, that day's observations are not included in that moon's dataset.
 
Inclement weather and other circumstances have also prevented data collection on certain days as indicated.
 
While angular distances do not have negative values, the signs represent whether that moon has a greater or lower [https://en.wikipedia.org/wiki/Horizontal_coordinate_system elevation] in the sky compared to Jupiter; this helps differentiate in which part of their orbits the moons are for the purposes of modeling.
 
==Modeling==
[[File:MWSineProjection.gif|right|thumb|450px|'''Left''': actual position of Jupiter and one of its moons as viewed from the north<br>'''Center''': apparent positions as viewed by a distant observer from Earth<br>'''Right''': the apparent position plotted over time]]
 
The equation used to model the motion of the moons is
 
<math display="block">d\left(t\right)=a\sin\left(\frac{2\pi}T\cdot t+\phi\right)</math>
 
where <math>d\left(t\right)</math> is the angular displacement between the moon and Jupiter as a function of time <math>t</math>, <math>a</math> is the amplitude of the moon's oscillation about Jupiter, proportional to its [https://en.wikipedia.org/wiki/Semi-major_and_semi-minor_axes semi-major axis], <math>T</math> is the period of the orbit, and <math>\phi</math> is the phase of the orbit.
 
A simple sinusoid is used to the model the moons' motions. By projecting uniform circular motion unto the [https://en.wikipedia.org/wiki/Celestial_sphere plane of the sky], a pure sinusoidal wave would be obtained. In order to use this model, the following assumptions must be made:
* The [https://en.wikipedia.org/wiki/Orbital_inclination inclination] of the moons' orbits are 90º.
* The orbital paths of the moons are circular, and as a result
** the moons have constant angular velocity.
** the moons maintain a constant distance from Jupiter.
* The distance between the Earth and the Jovian system remains constant.
 
The validity of these assumptions will be discussed [[#Validity of Model|later]].
 
Amplitude <math>a</math> is expressed in units of pixels, and is approximately proportional to angular distance. Time <math>t</math> is expressed in units of calendar days from the midnight starting 31 January 2022, Singapore civil time. For example, if an observation is made on 31 January 6:00 PM, <math>t=0.75\mathrm d</math>. Period <math>T</math> is also expressed in units of days (d). The phase <math>\phi</math> serves only to fit the sinusoidal function to the data and is not physically meaningful in this analysis.
 
==Results==
The following analyses were done in Wolfram Mathematica and Microsoft Excel.
 
===Fitted Parameters===
 
Using the Mathematica [https://reference.wolfram.com/language/ref/FindFit.html <code>FindFit</code>] function, the [[#Raw Data|collected data]] were fitted to the model, obtaining the following parameters, rounded to four significant figures:
 
{| class="wikitable" style="text-align: center
|-
! scope="col" style="width: 25%" | Satellite
! scope="col" style="width: 25%" | Amplitude<br/><math display="inline">a</math> (a.u.)
! scope="col" style="width: 25%" | Period<br/><math>T</math> (d)
! scope="col" style="width: 25%" | Phase<br/><math>\phi</math> (rad)
|-
! scope="row"| Io
| 98.38
| 1.770
| 5.701
|-
! scope="row"| Europa
| 154.7
| 3.555
| 4.807
|-
! scope="row"| Ganymede
| 249.3
| 7.196
| 5.975
|-
! scope="row"| Callisto
| 441.7
| 16.92
| 3.433
|}
 
Angular distance is expressed in arbitrary units (a.u.), not [https://en.wikipedia.org/wiki/Astronomical_unit astronomical units] (AU).
 
===Plots===
<gallery widths=400px heights=300px caption="Plots of Observations and Fitted Models>
File:MWIoPlot.png
File:MWEuropaPlot.png
File:MWGanymedePlot.png
File:MWCallistoPlot.png
</gallery>
All dates and times are expressed in Singapore standard time (UTC+08:00).
 
The red dots plot the observations made while the grey lines denote the fitted model. Note the differences in vertical axis scaling.
 
===Comparisons with Literature Values===
Visually, from the plots, the model of the moons' motions seem to fit the observed data well.
 
====Orbital Period====
{| class="wikitable" style="text-align: center
|-
! scope="col" style="width: 25%" | Satellite
! scope="col" style="width: 25%" | Observed Period<br>(d)
! scope="col" style="width: 25%" | Actual Period<br>(d)<ref name="nasamoonsofjupiter">NASA, [https://en.wikipedia.org/wiki/Jet_Propulsion_Laboratory JPL] (July 2013) [https://www.nasa.gov/sites/default/files/files/Moons_of_Jupiter_Lithograph.pdf ''Galilean Moons of Jupiter'']</ref>
! scope="col" style="width: 25%" | Discrepancy
|-
! scope="row"| Io
| 1.770
| 1.769
| 0.063%
|-
! scope="row"| Europa
| 3.555
| 3.551
| 0.104%
|-
! scope="row"| Ganymede
| 7.196
| 7.155
| 0.570%
|-
! scope="row"| Callisto
| 16.92
| 16.69
| 1.407%
|}
 
The discrepancies between the empirical data and the literature values are small but increase with increasing period.
 
====Orbital Radii====
Because the distance to Jupiter is not known in this exercise, the actual orbital radius (in kilometers) of each moon cannot be inferred from the observations alone. Instead the ratios of the orbital radii of the three outer moons to Io are taken and compared to the literature values.
 
{| class="wikitable" style="text-align: center
|-
! scope="col" style="width: 17%" | Satellite
! scope="col" style="width: 17%" | Observed Amplitude<br>(a.u.)
! scope="col" style="width: 16%" | Observed Ratio<br>(relative to Io)
! scope="col" style="width: 17%" | Actual Semi-major Axis<br>(km)<ref name="nasamoonsofjupiter"/>
! scope="col" style="width: 16%" | Actual Ratio<br>(relative to Io)
! scope="col" style="width: 17%" | Discrepancy
|-
! scope="row"| Io
| 98.38
| 1
| 422 000
| 1
| -
|-
! scope="row"| Europa
| 154.7
| 1.573
| 671 000
| 1.590
| 1.079%
|-
! scope="row"| Ganymede
| 249.3
| 2.534
| 1 070 000
| 2.536
| 0.075%
|-
! scope="row"| Callisto
| 441.7
| 4.490
| 1 883 000
| 4.462
| 0.625%
|}
 
Again, the discrepancies between the empirical data and the literature values are small.
 
====Orbital Resonance====
By expressing the orbital periods in multiples of the period of Io, the [https://en.wikipedia.org/wiki/Orbital_resonance orbital resonance phenomenon] between Io, Europa, and Ganymede can be identified.
 
{| class="wikitable" style="text-align: center
|-
! scope="col" style="width: 25%" | Satellite
! scope="col" style="width: 25%" | Observed Period<br>(relative to Io)
! scope="col" style="width: 25%" | Actual Period<br>(relative to Io)
! scope="col" style="width: 25%" | Discrepancy
|-
! scope="row"| Io
| 1
| 1
| -
|-
! scope="row"| Europa
| 2.008
| 2.007
| 0.041%
|-
! scope="row"| Ganymede
| 4.065
| 4.045
| 0.506%
|-
! scope="row"| Callisto
| 9.561
| 9.434
| 1.343%
|}
 
The relative periods between the three innermost moons are small integer multiples of each other—characteristic of bodies in orbital resonance. The oscillating gravitational influence they experience between themselves is a stable interaction that ensures their relative positions maintain over time. Callisto does not seem to participate in this arrangement.
 
====Kepler's Third Law====
[https://en.wikipedia.org/wiki/Kepler%27s_laws_of_planetary_motion#Third_law Kepler's third law] of planetary motion, in mathematical terms, reads that for all objects orbiting the same primary massive object,
 
<math display="block">T^2\propto a^3</math>
 
where <math>T</math> is the period of the orbit and <math>a</math> is the semi-major axis of the orbit.
 
Kepler published these laws with the Sun and planets in mind (hence "Kepler's laws of planetary motion"), but they also apply to moons orbiting planets as well.
 
{| class="wikitable" style="text-align: center
|-
! scope="col" style="width: 20%" | Satellite
! scope="col" style="width: 20%" | Observed Amplitude<br><math>a</math> (a.u.)
! scope="col" style="width: 20%" | Observed Period<br><math>T</math> (d)
! scope="col" style="width: 20%" | Constant of proportionality<br><math>\frac{a^3}{T^2}</math>
! scope="col" style="width: 20%" | Discrepancy from Mean
|-
! scope="row"| Io
| 98.38
| 1.770
| 303 900
| 1.539%
|-
! scope="row"| Europa
| 154.7
| 3.555
| 293 300
| 2.022%
|-
! scope="row"| Ganymede
| 249.3
| 7.193
| 299 100
| 0.065%
|-
! scope="row"| Callisto
| 441.7
| 16.92
| 301 000
| 0.547%
|}
 
=====Derivation=====
By equating [https://en.wikipedia.org/wiki/Newton's_law_of_universal_gravitation Newton's Law of universal gravitation] with [https://en.wikipedia.org/wiki/Centripetal_force centripetal acceleration],
 
<math display="block">G\frac{Mm}{r^2}=mr\omega^2</math>
 
where <math>G</math> is the [https://en.wikipedia.org/wiki/Gravitational_constant gravitational constant], <math>M</math> is the mass of Jupiter, <math>m</math> is the mass of the moon, <math>r</math> is the radius of the moon's orbit, and <math display="inline">\omega=\frac{2\pi}T</math> is the angular velocity of the moon.
 
[[File:MWKeplerPlot.png|right|thumb|400px|Log-log plot of observed period against orbital radius; the grey line of best fit connecting the four moons illustrates Kepler's third law]]
 
By expressing angular velocity <math>\omega</math> in terms of orbital period <math>T</math> and rearranging,
 
<math display="block>\begin{align}G\frac{Mm}{r^2}&=mr\left(\frac{2\pi}T\right)^2\\T^2&=\left(\frac{4\pi^2}{GM}\right)r^3\\T^2&\propto r^3\end{align}</math>
 
So proving Kepler's third law. Additionally, by linearizing the equation,
 
<math display="block">\begin{align}\lg\left[T^2\right]&=\lg\left[\left(\frac{4\pi^2}{GM}\right)r^3\right]\\2\lg T&=3\lg r+\lg\left(\frac{4\pi^2}{GM}\right)\\\lg T&=\frac32\lg r+\lg\left(\frac{2\pi}{\sqrt{GM}}\right)\end{align}</math>
 
By plotting <math>\lg T</math> against <math>\lg r</math>, a straight-line graph with gradient <math display="inline">\frac32</math> would be obtained with vertical intercept <math display="inline">\lg\left(\frac{2\pi}{\sqrt{GM}}\right)</math>. If the actual orbital radii of the moons were known (they cannot be calculated from these experimental observations alone <sup>[[#Orbital Radii|[see why]]]</sup>), the mass of Jupiter <math>M</math> could be calculated from the plot.
 
The plot on the right has a line of best fit with gradient 1.494, within 0.414% of the theoretical prediction of 1.5.
 
==Discussion==
 
===Validity of Model===
 
The [[#Modeling|model used]] is a very simplistic approximation of the complex motions traced by the Galilean moons. Yet, the fitted parameters agree well with the literature values and the resulting equation has reasonably good predictive power and seems to fit the observed data well.
 
However, on closer inspection, several problems with the model can be identified.
 
====Inclination====
 
[[File:MWInclinationAnimation.gif|right|thumb|450px|
'''Top''': profile view of the orbital plane as viewed from the side with Earth on the left; the orange dot is Jupiter while the red dot is one of its moons<br>'''Middle''': apparent positions as viewed by a distant observer from Earth<br>'''Bottom''': the apparent angular separation in arbitrary units and its plot over time]]
 
The model assumes that the [https://en.wikipedia.org/wiki/Orbital_inclination inclination] of the moons' orbits relative to the line of sight with the Earth is zero. In other words, it assumes the moon should pass directly over the center of Jupiter when in actuality it will not, leading to inaccurate estimations of the orbital phase, especially so when the moon appears near to the center of Jupiter.
 
However, the inclination of these moons are small (between 0.04º and 0.47º<ref name="nasasmallworlds">NASA, [https://en.wikipedia.org/wiki/Goddard_Space_Flight_Center GSFC], [https://en.wikipedia.org/wiki/NASA_Space_Science_Data_Coordinated_Archive NSSDC] (April 2016) [https://nssdc.gsfc.nasa.gov/planetary/factsheet/galileanfact_table.html ''Solar System Small Worlds Fact Sheet'']</ref>) and this effect can be safely ignored.
 
====Motion of the Earth====
The model ignores the fact that the Earth had been moving away from Jupiter during the duration of observation and so the apparent orbital radii would appear to decrease over time.
 
However, since the duration of observation is not very long (about two weeks), and since the Earth–Jupiter distance was about to reach its maximum (stationary point when Jupiter is on the opposite side of the Sun from the Earth), this effect not very noticeable.
 
===Sampling Frequency===
Because the sampling frequency (once per day at best) is not much greater than the orbital frequencies of the moons (once per 1.7 days for Io), the [https://en.wikipedia.org/wiki/Nyquist–Shannon_sampling_theorem Nyquist–Shannon sampling theorem] implies that the precise orbital period <math>T</math> cannot be accurately and uniquely determined for Io and perhaps even Europa (once per 3.55 days) as well.
 
The sampling rate is adequately high compared to the orbital frequencies of Ganymede (once per 7.16 days) and Callisto (once per 16.69 days) such that the computer can determine their orbital parameters autonomously. However, for Io and Europa, a good starting guess must be provided before the computer could find a good fit.
 
To improve the computer's estimates and ability to guess, more samples can be taken over a longer period of time or, perhaps more effectively, over the same period of time but with greater frequency. This was not possible in this experiment , however, as the geometry of the Sun–Earth–Jupiter system at the time of the observations meant that Jupiter and its moons were only visible for an hour or so per day. In future, this experiment should be run during a part of the calendar year where Jupiter is in [https://en.wikipedia.org/wiki/Opposition_(astronomy) opposition] (directly overhead at solar midnight), allowing for multiple hours of observations per night. Observing Jupiter in opposition also has the added advantage that the Earth–Jupiter distance would be the smallest, thus Jupiter appears biggest in the sky (0.0136º<ref>Wolfram Research (2014) [https://reference.wolfram.com/language/ref/PlanetData.html <code>PlanetData</code>], Wolfram Language function</ref> at opposition on 26 September 2022<ref>Space.com (December 2021) [https://www.space.com/39240-when-to-see-planets-in-the-sky.html#section-jupiter ''When, where and how to see the planets in the 2022 night sky'']</ref> vs 0.00897º at [https://en.wikipedia.org/wiki/Conjunction_(astronomy)#Superior_and_inferior solar conjunction] on 5 March 2022).
 
===Diffraction Limits===
The physical capabilities of the lens–camera system are assessed in this section.
 
The effective diameter of the [https://en.wikipedia.org/wiki/Aperture aperture]<math>D</math> can be found using the [https://en.wikipedia.org/wiki/F-number f-number] formula,
 
<math display="block">\begin{align}D&=\frac fN\\&=\frac{600}{6.3}\\&=95\;\mathrm{mm}\end{align}</math>
 
where <math>f</math> is the [https://en.wikipedia.org/wiki/Focal_length focal length] of the lens system and <math>N</math> is the f-number. 95mm is also the diameter of the frontmost glass element on the lens.
 
An optics system is considered [https://en.wikipedia.org/wiki/Diffraction-limited_system diffraction limited] when the resolution performance is limited by physics of diffraction. An [https://en.wikipedia.org/wiki/Airy_disk Airy disk] is the smallest sized disk to which a point light source can be focused on the imaging plane. The angular half-width of an Airy disk <math>\theta</math> can be calculated with
 
<math display="block">\begin{align}\theta&=\arcsin\left(1.22\frac\lambda D\right)\\
&=\arcsin\left(1.22\frac{600\;\mathrm{nm}}{95\;\mathrm{mm}}\right)\\
&=7.7\times10^{-6}\;\mathrm{rad}\end{align}</math>
 
where <math>\lambda=600\;\mathrm{nm}</math> is the wavelength of light, in this case taken to be orange as Jupiter is somewhat orange<ref>NASA, [https://en.wikipedia.org/wiki/European_Space_Agency ESA], [https://en.wikipedia.org/wiki/Goddard_Space_Flight_Center GSFC] (2014) [https://en.wikipedia.org/wiki/Hubble_Space_Telescope Hubble Space Telescope] image: [https://esahubble.org/images/heic1410a/ ''Jupiter and its shrunken Great Red Spot'']</ref>. From there, using geometry, the diameter <math>d_A</math> of the Airy disk on the sensor would be
 
<math display="block">\begin{align}d_A&=2f\tan\theta\\
&=2\left(600\;\mathrm{mm}\right)\tan\left(7.7\times10^{-6}\;\mathrm{rad}\right)\\
&=9.2\;\mathrm{\mu m}\\
&>4.31\;\mathrm{\mu m}\text{ (pixel size)}\end{align}</math>
 
Since the Airy disk more than twice the size of two pixels on the sensor, the optical system '''is diffraction limited'''.
 
While this does not affect the precision to which the positions of the moons can be determined, it is intriguing to note that consumer-grade, unspecialized hardware is limited not by the quality of optics or sensor technology, but the fundamental laws of physics.
 
By enlarging the diameter of the aperture, a higher resolution can be achieved. Additionally, more light can be collected, increasing the signal-to-noise ratio and improving image quality. However, this would lead to increased bulk in the glass lens and, of course, substantially increased cost. This is why almost all modern telescopes are [https://en.wikipedia.org/wiki/Reflecting_telescope reflecting] (using large mirrors to focus light) instead of refracting (using large pieces of glass to focus light). This also has the added benefit of eliminating many drawbacks of large glass such as [https://en.wikipedia.org/wiki/Chromatic_aberration chromatic aberration].
 
===Angular Field of View===
A few methods can be used to determine the horizontal field of view of the setup. Being able to measure angles precisely is an important part of astronomy and can be used to infer physical distances using trigonometry assuming the absolute distances to the objects of interest are known.
 
====Manufacturer's Claim====
According to the manufacturer, the lens has a field of view (after accounting for the [https://en.wikipedia.org/wiki/APS-C#Multiplier_factors APS-C crop factor]) of <math display="inline">\phi=\frac{4.1}{1.6}=2.56^\circ</math><ref name="sigmalens"/>.
 
====Focal Length====
 
Using geometry, the angular resolution <math>\rho</math> of the setup can be determined.
 
<math display="block">\begin{align}\rho&=\arctan\left(\frac{\text{pixel size}}f\right)\\
&=\arctan\left(\frac{4.31\;\mathrm{\mu m}}{600\;\mathrm{mm}}\right)\\
&=0.000412^\circ\\
&=0^\circ\,0'\,1.48''\end{align}</math>
 
This is the angle corresponding to a single pixel in a captured image. Multiplying that with the number of pixels horizontally <math>n=5184\;\mathrm{px}</math><ref name="canoncamera"/>, the field of view would be <math>\phi=0.000412^\circ\times5184=2.13^\circ</math>.
 
====Imaging====
 
[[File:MWMoonCapture.jpg|right|thumb|300px|Image of the full moon used to calculate the angular field of view of the setup]]
 
A photo of the full moon was taken on 17 Feb. Knowing the distance to the moon <math>r_\mathrm{Moon}=389\,110\;\mathrm{km}</math><ref>Wolfram Research (2014) [https://reference.wolfram.com/language/ref/PlanetaryMoonData.html <code>PlanetaryMoonData</code>], Wolfram Language function</ref> and the diameter of the moon <math>d_\mathrm{Moon}=3476.2\;\mathrm{km}</math><ref>NASA, [https://en.wikipedia.org/wiki/Goddard_Space_Flight_Center GSFC], [https://en.wikipedia.org/wiki/NASA_Space_Science_Data_Coordinated_Archive NSSDC] (December 2021) [https://nssdc.gsfc.nasa.gov/planetary/factsheet/galileanfact_table.html ''Moon Fact Sheet'']</ref>, the angular diameter of the moon <math>\theta_\mathrm{Moon}</math> can be calculated by
 
<math display="block">\begin{align}\theta_\mathrm{Moon}&=\arctan\left(\frac{3476.2}{389\,110}\right)\\
&=0.5119^\circ\end{align}</math>
 
Since the moon occupies <math>n=1124\;\mathrm{px}</math> of this image horizontally, the field of view would be <math display="inline">\phi=\frac{5184}{1124}\times0.5119^\circ=2.36^\circ</math>.
 
====Discussion====
 
These methods assume that the image contains no [https://en.wikipedia.org/wiki/Distortion_(optics) lens distortion], especially at the fringes of the image, far from the [https://en.wikipedia.org/wiki/Optical_axis optical axis], which may not be true and lead to the discrepancies between the derived values above.
 
==Conclusion==
 
[[File:MWMoonJupiterConjunction.jpg|right|thumb|300px|Conjunction between the moon and the Jovian system, showing their relative sizes; the four Galilean moons are also visible in the full-sized image]]
 
Through the course of this project, the four largest moons of Jupiter were surveyed over the course of two weeks. Therefrom, the following results were obtained:
 
* The orbital periods of the four Galilean moons, Io, Europa, Ganymede, and Callisto were found to be, respectively, <math>T_\mathrm{Io}=1.770\;\mathrm{d}</math>, <math>T_\mathrm{Europa}=3.555\;\mathrm{d}</math>, <math>T_\mathrm{Ganymede}=7.196\;\mathrm{d}</math>, and <math>T_\mathrm{Callisto}=16.92\;\mathrm{d}</math>, within 1.407% of literature values.
 
* Their relative orbital radii, in multiples of the orbital radius of Io were found to be, respectively, <math>r_\mathrm{Io}=1\;r_\mathrm{Io}</math>, <math>r_\mathrm{Europa}=1.573\;r_\mathrm{Io}</math>, <math>r_\mathrm{Ganymede}=2.534\;r_\mathrm{Io}</math>, <math>r_\mathrm{Callisto}=4.490\;r_\mathrm{Io}</math>, within 1.079% of literature values.
 
* Their relative orbital periods, in multiples of the orbital period of Io were found to be, respectively, <math>T_\mathrm{Io}=1\;T_\mathrm{Io}</math>, <math>T_\mathrm{Europa}=2.008\;T_\mathrm{Io}</math>, <math>T_\mathrm{Ganymede}=4.065\;T_\mathrm{Io}</math>, <math>T_\mathrm{Callisto}=9.561\;T_\mathrm{Io}</math>, within 1.343% of literature values, in so doing demonstrating the orbital resonance phenomenon.
 
* Kepler's third law was validated to within 2.022%.
 
The fact that these observations were performed with consumer equipment not specifically designed for astronomical observations, that yet still managed to produce results reasonably close to the established literature values, shows that astronomy is among the most accessible in all sub-fields of physics. It requires no specialized equipment available only to scientists in laboratories, only access to the night sky, the [https://en.wikipedia.org/wiki/Common_heritage_of_humanity common heritage of humanity].
 
==Citations==

Latest revision as of 04:39, 30 April 2022

Jupiter and the four Galilean moons

Galileo Galilei’s discovery of celestial bodies that orbit something other than the Earth marked the beginning of the end of the geocentric model of the universe. In this project, we will perform the same observations on those moons as Galileo did 400 years ago.

Team Members

Objectives

History

Animation of Galileo's observations as drawn in his book Sidereus Nuncius in 1610

Jupiter is a planet that was known to the ancients as it is the fourth brightest object in the sky behind the Sun, the Moon, and Venus, and is visible to the naked eye. In 1609, astronomer Galileo Galilei made improvements to his telescope, which was the best in the world at the time. With the telescope, over the course of several weeks, he observed four “stars” moving in a line around Jupiter and was persuaded after just four days of observations that they were not stationary stars but instead objects orbiting the planet Jupiter. German astronomer Simon Marius also made this discovery independently and at the same time, but he is lesser known because he published his observations only after Galileo did. He had the privilege, however, of having his names of these moons (Io, Europa, Ganymede, Callisto) be widely accepted, much to Galileo's dismay.

The discovery that there exist objects that do not orbit the Earth, along with several other breakthroughs, was the cornerstone of the Copernican Revolution, which was the transition from the geocentric to the heliocentric model, with the Sun instead of the Earth at the center of the universe.

Since then, astronomers have discovered a total of 80 moons that orbit Jupiter, 23 of which have no official names.

Background

The four Galilean moons, in increasing order of orbital radius, are Io, Europa, Ganymede, and Callisto, of which, Io, Callisto, and Ganymede are larger than the Earth's Moon in both mass and diameter. Ganymede is even larger than the planet Mercury. Because of their size, they can be observed from Earth with a small telescope or even a pair of powerful binoculars.

The main challenge faced when trying to observe these moons is the sheer brightness of Jupiter relative to the dim moons, and their small angular separation from Jupiter. This is especially true for the innermost moons Io and Europa.

Apparatus

Setup used to image Jupiter and its satellites

No specialized equipment was used in this experiment to make observations of the planet and its moons.

  • Canon EOS 60D DSLR[1]
    • APS-C sensor (1.6× crop)
    • Pixel size 4.31µm
  • Sigma 150–600mm f/5–6.3[2]
    • Lens diameter 95mm
    • Field of view: 16.4º–4.1º (effective 10.2º–2.6º)
  • Low-end consumer tripod
    • Spirit level
  • Remote shutter

Exposure Settings

  • Manual
  • Focal length: 600mm (effective 960mm)
  • 1/15s exposure
  • f/6.3 aperture
  • ISO 3200 sensitivity
  • Manual focus
  • Image stabilization off

The exposure time cannot be set for too long as motion blur due to the Earth's rotation will be visible in the image.

Data

Screenshot of a camera live view of Jupiter
An actual image captured on 1 Feb (click to zoom in)

Procedure

A series of observations were performed with the telephoto lens mounted on the camera. Weather permitting, exposures were taken every day within an hour of sunset for two and a half calendar weeks from 31 January 2022 to 17 February 2022.

At the time of these observations, Earth's motion around the Sun (and to a lesser extent, Jupiter's as well) caused the apparent position of Jupiter in the sky to drift closer to the sun every day, forcing observations to occur closer and closer to sunset. The last observation on 17 February was taken at 7:59 PM, when the sun had barely set and the sky was still relatively bright (compared to the dim moons). It was decided that no more observations could be taken thereafter.

The moons in the captured images were identified using Stellarium or Night Sky.

By inspecting the moons and their positions relative to Jupiter in the captured images, their apparent angular separation can be determined for every observation. The data are presented in the table below.

Raw Data

Observation Angular distance between moon and center of Jupiter (arbitrary units)
Io Europa Ganymede Callisto
2022-01-31 68.9 occultation 102.7 -254.2
2022-02-01 -32.2 155.4 245.7 -371.0
2022-02-02 no data collected
2022-02-03 50.6 -127.9 23.8 -437.9
2022-02-04 -76.3 113.8 -175.9 -392.3
2022-02-05 96.4 83.2 -251.8 -284.4
2022-02-06 no data collected
2022-02-07 86.2 -31.0 62.2 25.5
2022-02-08 -55.7 154.5 224.8 183.0
2022-02-09 no data collected
2022-02-10 occultation -145.8 68.0 405.6
2022-02-11 -52.4 87.2 -140.1 437.5
2022-02-12 81.0 108.1 -246.5 408.5
2022-02-13 no data collected
2022-02-14 no data collected
2022-02-15 -75.2 145.8 195.9 40.0
2022-02-16 no data collected
2022-02-17 transit -147.2 101.8 -252.2

When a moon is occulted by Jupiter, it is not visible and its position cannot be determined with precision. Similarly, during a transit, Jupiter outshines the moon and its precise position also cannot be determined. Either way, that day's observations are not included in that moon's dataset.

Inclement weather and other circumstances have also prevented data collection on certain days as indicated.

While angular distances do not have negative values, the signs represent whether that moon has a greater or lower elevation in the sky compared to Jupiter; this helps differentiate in which part of their orbits the moons are for the purposes of modeling.

Modeling

Left: actual position of Jupiter and one of its moons as viewed from the north
Center: apparent positions as viewed by a distant observer from Earth
Right: the apparent position plotted over time

The equation used to model the motion of the moons is

where is the angular displacement between the moon and Jupiter as a function of time , is the amplitude of the moon's oscillation about Jupiter, proportional to its semi-major axis, is the period of the orbit, and is the phase of the orbit.

A simple sinusoid is used to the model the moons' motions. By projecting uniform circular motion unto the plane of the sky, a pure sinusoidal wave would be obtained. In order to use this model, the following assumptions must be made:

  • The inclination of the moons' orbits are 90º.
  • The orbital paths of the moons are circular, and as a result
    • the moons have constant angular velocity.
    • the moons maintain a constant distance from Jupiter.
  • The distance between the Earth and the Jovian system remains constant.

The validity of these assumptions will be discussed later.

Amplitude is expressed in units of pixels, and is approximately proportional to angular distance. Time is expressed in units of calendar days from the midnight starting 31 January 2022, Singapore civil time. For example, if an observation is made on 31 January 6:00 PM, . Period is also expressed in units of days (d). The phase serves only to fit the sinusoidal function to the data and is not physically meaningful in this analysis.

Results

The following analyses were done in Wolfram Mathematica and Microsoft Excel.

Fitted Parameters

Using the Mathematica FindFit function, the collected data were fitted to the model, obtaining the following parameters, rounded to four significant figures:

Satellite Amplitude
(a.u.)
Period
(d)
Phase
(rad)
Io 98.38 1.770 5.701
Europa 154.7 3.555 4.807
Ganymede 249.3 7.196 5.975
Callisto 441.7 16.92 3.433

Angular distance is expressed in arbitrary units (a.u.), not astronomical units (AU).

Plots

All dates and times are expressed in Singapore standard time (UTC+08:00).

The red dots plot the observations made while the grey lines denote the fitted model. Note the differences in vertical axis scaling.

Comparisons with Literature Values

Visually, from the plots, the model of the moons' motions seem to fit the observed data well.

Orbital Period

Satellite Observed Period
(d)
Actual Period
(d)[3]
Discrepancy
Io 1.770 1.769 0.063%
Europa 3.555 3.551 0.104%
Ganymede 7.196 7.155 0.570%
Callisto 16.92 16.69 1.407%

The discrepancies between the empirical data and the literature values are small but increase with increasing period.

Orbital Radii

Because the distance to Jupiter is not known in this exercise, the actual orbital radius (in kilometers) of each moon cannot be inferred from the observations alone. Instead the ratios of the orbital radii of the three outer moons to Io are taken and compared to the literature values.

Satellite Observed Amplitude
(a.u.)
Observed Ratio
(relative to Io)
Actual Semi-major Axis
(km)[3]
Actual Ratio
(relative to Io)
Discrepancy
Io 98.38 1 422 000 1 -
Europa 154.7 1.573 671 000 1.590 1.079%
Ganymede 249.3 2.534 1 070 000 2.536 0.075%
Callisto 441.7 4.490 1 883 000 4.462 0.625%

Again, the discrepancies between the empirical data and the literature values are small.

Orbital Resonance

By expressing the orbital periods in multiples of the period of Io, the orbital resonance phenomenon between Io, Europa, and Ganymede can be identified.

Satellite Observed Period
(relative to Io)
Actual Period
(relative to Io)
Discrepancy
Io 1 1 -
Europa 2.008 2.007 0.041%
Ganymede 4.065 4.045 0.506%
Callisto 9.561 9.434 1.343%

The relative periods between the three innermost moons are small integer multiples of each other—characteristic of bodies in orbital resonance. The oscillating gravitational influence they experience between themselves is a stable interaction that ensures their relative positions maintain over time. Callisto does not seem to participate in this arrangement.

Kepler's Third Law

Kepler's third law of planetary motion, in mathematical terms, reads that for all objects orbiting the same primary massive object,

where is the period of the orbit and is the semi-major axis of the orbit.

Kepler published these laws with the Sun and planets in mind (hence "Kepler's laws of planetary motion"), but they also apply to moons orbiting planets as well.

Satellite Observed Amplitude
(a.u.)
Observed Period
(d)
Constant of proportionality
Discrepancy from Mean
Io 98.38 1.770 303 900 1.539%
Europa 154.7 3.555 293 300 2.022%
Ganymede 249.3 7.193 299 100 0.065%
Callisto 441.7 16.92 301 000 0.547%
Derivation

By equating Newton's Law of universal gravitation with centripetal acceleration,

where is the gravitational constant, is the mass of Jupiter, is the mass of the moon, is the radius of the moon's orbit, and is the angular velocity of the moon.

Log-log plot of observed period against orbital radius; the grey line of best fit connecting the four moons illustrates Kepler's third law

By expressing angular velocity in terms of orbital period and rearranging,

So proving Kepler's third law. Additionally, by linearizing the equation,

By plotting against , a straight-line graph with gradient would be obtained with vertical intercept . If the actual orbital radii of the moons were known (they cannot be calculated from these experimental observations alone [see why]), the mass of Jupiter could be calculated from the plot.

The plot on the right has a line of best fit with gradient 1.494, within 0.414% of the theoretical prediction of 1.5.

Discussion

Validity of Model

The model used is a very simplistic approximation of the complex motions traced by the Galilean moons. Yet, the fitted parameters agree well with the literature values and the resulting equation has reasonably good predictive power and seems to fit the observed data well.

However, on closer inspection, several problems with the model can be identified.

Inclination

Top: profile view of the orbital plane as viewed from the side with Earth on the left; the orange dot is Jupiter while the red dot is one of its moons
Middle: apparent positions as viewed by a distant observer from Earth
Bottom: the apparent angular separation in arbitrary units and its plot over time

The model assumes that the inclination of the moons' orbits relative to the line of sight with the Earth is zero. In other words, it assumes the moon should pass directly over the center of Jupiter when in actuality it will not, leading to inaccurate estimations of the orbital phase, especially so when the moon appears near to the center of Jupiter.

However, the inclination of these moons are small (between 0.04º and 0.47º[4]) and this effect can be safely ignored.

Motion of the Earth

The model ignores the fact that the Earth had been moving away from Jupiter during the duration of observation and so the apparent orbital radii would appear to decrease over time.

However, since the duration of observation is not very long (about two weeks), and since the Earth–Jupiter distance was about to reach its maximum (stationary point when Jupiter is on the opposite side of the Sun from the Earth), this effect not very noticeable.

Sampling Frequency

Because the sampling frequency (once per day at best) is not much greater than the orbital frequencies of the moons (once per 1.7 days for Io), the Nyquist–Shannon sampling theorem implies that the precise orbital period cannot be accurately and uniquely determined for Io and perhaps even Europa (once per 3.55 days) as well.

The sampling rate is adequately high compared to the orbital frequencies of Ganymede (once per 7.16 days) and Callisto (once per 16.69 days) such that the computer can determine their orbital parameters autonomously. However, for Io and Europa, a good starting guess must be provided before the computer could find a good fit.

To improve the computer's estimates and ability to guess, more samples can be taken over a longer period of time or, perhaps more effectively, over the same period of time but with greater frequency. This was not possible in this experiment , however, as the geometry of the Sun–Earth–Jupiter system at the time of the observations meant that Jupiter and its moons were only visible for an hour or so per day. In future, this experiment should be run during a part of the calendar year where Jupiter is in opposition (directly overhead at solar midnight), allowing for multiple hours of observations per night. Observing Jupiter in opposition also has the added advantage that the Earth–Jupiter distance would be the smallest, thus Jupiter appears biggest in the sky (0.0136º[5] at opposition on 26 September 2022[6] vs 0.00897º at solar conjunction on 5 March 2022).

Diffraction Limits

The physical capabilities of the lens–camera system are assessed in this section.

The effective diameter of the aperture can be found using the f-number formula,

where is the focal length of the lens system and is the f-number. 95mm is also the diameter of the frontmost glass element on the lens.

An optics system is considered diffraction limited when the resolution performance is limited by physics of diffraction. An Airy disk is the smallest sized disk to which a point light source can be focused on the imaging plane. The angular half-width of an Airy disk can be calculated with

where is the wavelength of light, in this case taken to be orange as Jupiter is somewhat orange[7]. From there, using geometry, the diameter of the Airy disk on the sensor would be

Since the Airy disk more than twice the size of two pixels on the sensor, the optical system is diffraction limited.

While this does not affect the precision to which the positions of the moons can be determined, it is intriguing to note that consumer-grade, unspecialized hardware is limited not by the quality of optics or sensor technology, but the fundamental laws of physics.

By enlarging the diameter of the aperture, a higher resolution can be achieved. Additionally, more light can be collected, increasing the signal-to-noise ratio and improving image quality. However, this would lead to increased bulk in the glass lens and, of course, substantially increased cost. This is why almost all modern telescopes are reflecting (using large mirrors to focus light) instead of refracting (using large pieces of glass to focus light). This also has the added benefit of eliminating many drawbacks of large glass such as chromatic aberration.

Angular Field of View

A few methods can be used to determine the horizontal field of view of the setup. Being able to measure angles precisely is an important part of astronomy and can be used to infer physical distances using trigonometry assuming the absolute distances to the objects of interest are known.

Manufacturer's Claim

According to the manufacturer, the lens has a field of view (after accounting for the APS-C crop factor) of [2].

Focal Length

Using geometry, the angular resolution of the setup can be determined.

This is the angle corresponding to a single pixel in a captured image. Multiplying that with the number of pixels horizontally [1], the field of view would be .

Imaging

Image of the full moon used to calculate the angular field of view of the setup

A photo of the full moon was taken on 17 Feb. Knowing the distance to the moon [8] and the diameter of the moon [9], the angular diameter of the moon can be calculated by

Since the moon occupies of this image horizontally, the field of view would be .

Discussion

These methods assume that the image contains no lens distortion, especially at the fringes of the image, far from the optical axis, which may not be true and lead to the discrepancies between the derived values above.

Conclusion

Conjunction between the moon and the Jovian system, showing their relative sizes; the four Galilean moons are also visible in the full-sized image

Through the course of this project, the four largest moons of Jupiter were surveyed over the course of two weeks. Therefrom, the following results were obtained:

  • The orbital periods of the four Galilean moons, Io, Europa, Ganymede, and Callisto were found to be, respectively, , , , and , within 1.407% of literature values.
  • Their relative orbital radii, in multiples of the orbital radius of Io were found to be, respectively, , , , , within 1.079% of literature values.
  • Their relative orbital periods, in multiples of the orbital period of Io were found to be, respectively, , , , , within 1.343% of literature values, in so doing demonstrating the orbital resonance phenomenon.
  • Kepler's third law was validated to within 2.022%.

The fact that these observations were performed with consumer equipment not specifically designed for astronomical observations, that yet still managed to produce results reasonably close to the established literature values, shows that astronomy is among the most accessible in all sub-fields of physics. It requires no specialized equipment available only to scientists in laboratories, only access to the night sky, the common heritage of humanity.

Citations

  1. 1.0 1.1 Canon Camera Museum (2010) Interchangeable Lens Digital Cameras - Digital SLR Camera: EOS 60D Specifications
  2. 2.0 2.1 Sigma Corporation (2018) 150-600mm F5-6.3 DG OS HSM Contemporary Lens Specifications
  3. 3.0 3.1 NASA, JPL (July 2013) Galilean Moons of Jupiter
  4. NASA, GSFC, NSSDC (April 2016) Solar System Small Worlds Fact Sheet
  5. Wolfram Research (2014) PlanetData, Wolfram Language function
  6. Space.com (December 2021) When, where and how to see the planets in the 2022 night sky
  7. NASA, ESA, GSFC (2014) Hubble Space Telescope image: Jupiter and its shrunken Great Red Spot
  8. Wolfram Research (2014) PlanetaryMoonData, Wolfram Language function
  9. NASA, GSFC, NSSDC (December 2021) Moon Fact Sheet