Stability of Lagrange Points:James Webb Space Telescope Cacolici,Gianna Nicole;Hanson,Jake;Lejoly,Cassandra;Pearson,Kyle Alexander;Reynolds,Katherine ABSTRACT We continue to present a detailed procedure for analyzing the stability of Lagrange points both analytically and computationally.We verified the stability both analytically in a 2D Earth-Sun system, and computationally with a 3D N-body code.We tested the stability over several periods of time and plotted the path the orbits would take in such case.We also discuss the consequences of these stabilities with respect to the James Webb Space Telescope. I.INTRODUCTION Astronomy has been as a subject studied by humans for a few thousands of years.One of the most fascinating parts of astronomy is the study of the retrograde motion of the planets in the sky. Among the earliest recorded studies of the motion of planets was by Ptolemy around the year 150 CE (Jones).Ptolemy's geocentric model of the universe was accepted for many centuries following his death.This model placed the Earth at the center of the universe with the Sun and planets rotating around the Earth in epicycles.It was not until Copernicus came along (1473-1543)that the idea of a heliocentric system came along.He believed that the size of the planets'orbits was proportional to the distance from the Sun(Rabin 2004).This model,which features the Sun in the center and the planets rotating around it,was not accepted initially,but became the accepted belief of the universe. It was not long until the telescope was invented and people started to truly record the motion of the planets.Tycho Brahe was one of the first to build an observatory to study the planets.He made very detailed observations of the planets and was able to predict their motions with an accuracy never seen before(Van Helden 1995).His results were published after his death and were used by Johannes Kepler,his student,to create the laws of planetary motion,which make use of elliptical orbits.To this day,this is the accepted idea of how planets orbit our Sun.These motions,although predicted by Kepler's Laws,can be predicted by Newton's Laws of Gravity.These explain how orbits can be stable. These orbits are not the only places where objects can exist.Joseph-Louis Lagrange was the first to publish his theories of how three bodies can interact in a system(Lagrange 1772).Although the three body problem is not solvable analytically,Lagrange was able to determine points where all the forces balance out.An object placed at these"Lagrange points"feels no force.These are the points we will be studying in our paper
Stability of Lagrange Points: James Webb Space Telescope Cacolici, Gianna Nicole; Hanson, Jake; Lejoly, Cassandra; Pearson, Kyle Alexander; Reynolds, Katherine ABSTRACT We continue to present a detailed procedure for analyzing the stability of Lagrange points both analytically and computationally. We verified the stability both analytically in a 2D EarthSun system, and computationally with a 3D Nbody code. We tested the stability over several periods of time and plotted the path the orbits would take in such case. We also discuss the consequences of these stabilities with respect to the James Webb Space Telescope. I. INTRODUCTION Astronomy has been as a subject studied by humans for a few thousands of years. One of the most fascinating parts of astronomy is the study of the retrograde motion of the planets in the sky. Among the earliest recorded studies of the motion of planets was by Ptolemy around the year 150 CE (Jones). Ptolemy’s geocentric model of the universe was accepted for many centuries following his death. This model placed the Earth at the center of the universe with the Sun and planets rotating around the Earth in epicycles. It was not until Copernicus came along (14731543) that the idea of a heliocentric system came along. He believed that the size of the planets’ orbits was proportional to the distance from the Sun (Rabin 2004). This model, which features the Sun in the center and the planets rotating around it, was not accepted initially, but became the accepted belief of the universe. It was not long until the telescope was invented and people started to truly record the motion of the planets. Tycho Brahe was one of the first to build an observatory to study the planets. He made very detailed observations of the planets and was able to predict their motions with an accuracy never seen before (Van Helden 1995). His results were published after his death and were used by Johannes Kepler, his student, to create the laws of planetary motion, which make use of elliptical orbits. To this day, this is the accepted idea of how planets orbit our Sun. These motions, although predicted by Kepler’s Laws, can be predicted by Newton’s Laws of Gravity. These explain how orbits can be stable. These orbits are not the only places where objects can exist. JosephLouis Lagrange was the first to publish his theories of how three bodies can interact in a system (Lagrange 1772). Although the three body problem is not solvable analytically, Lagrange was able to determine points where all the forces balance out. An object placed at these “Lagrange points” feels no force. These are the points we will be studying in our paper
14 L3 Sun L1 L2 Earth L5 (http://www.labspaces.net/blog/895/The_Lagrange_points Nature s parking_spots) There are five points in which all the forces balance out between two large bodies.For the purpose of this project,the Lagrange points considered will be those caused by the Sun and Earth.The Lagrange points all occur on the same 2D plane.The first three are in a straight line that passes through the Sun and Earth.The first point,labeled L1,occurs between the two bodies,L2 occurs just beyond the Earth,with both points occurring an equal distance from the Earth.L3 is found outside Earth's orbital path beyond the Sun.The last two points are located in the Earth's path via equilateral triangles. L4 is in line with where the Earth is rotating to and L5 is found trailing behind.All of these points are in constant rotation along with the Earth.The points are able to be viewed as they are in the diagram above by adopting a co-rotating frame.This frame holds the Sun and Earth fixed allowing for an easier visualization of the problem.The positions of the five Lagrange points are as follows: ((uming that=0) Point Location LI (R1-()],) L2 (R[1+()门,0) L3 (R[1-(ě)],) L4 ([*]r) L5 ([*]-r)
(http://www.labspaces.net/blog/895/The_Lagrange_points_____Nature___s_parking_spots) There are five points in which all the forces balance out between two large bodies. For the purpose of this project, the Lagrange points considered will be those caused by the Sun and Earth. The Lagrange points all occur on the same 2D plane. The first three are in a straight line that passes through the Sun and Earth. The first point, labeled L1, occurs between the two bodies, L2 occurs just beyond the Earth, with both points occurring an equal distance from the Earth. L3 is found outside Earth’s orbital path beyond the Sun. The last two points are located in the Earth’s path via equilateral triangles. L4 is in line with where the Earth is rotating to and L5 is found trailing behind. All of these points are in constant rotation along with the Earth. The points are able to be viewed as they are in the diagram above by adopting a corotating frame. This frame holds the Sun and Earth fixed allowing for an easier visualization of the problem. The positions of the five Lagrange points are as follows: (assuming that α = ) M2 M1+M2 Point Location L1 (R[1 − (3 α ) 3 1 ],0) L2 (R[1 + (3 α ) 3 1 ],0) L3 (− R[1 − (12 5α ) 3 1 ], 0) L4 ( R 2 R [M1+M2 M1−M2 ], 2 √3 ) L5 ( − R 2 R [M1+M2 M1−M2 ], 2 √3 )
Our study of Lagrange points can be applied to the James Webb Space Telescope(JWST) since it will be placed at the Earth-Sun L2 point.JWST is set to launch in 2018 and will be the successor to the Hubble Space Telescope and Spitzer Space Telescope.Since 1990,Hubble has been capturing light from the UV,visible,and near-IR,covering wavelengths from 0.1 to 1.7 microns.JWST will aim to improve previous work by acquiring more precise measurements in the IR spanning 0.6 to 29 microns. As light from stars and galaxies travels through space,its wavelength increases.Therefore,by capturing light with longer-wavelengths,JWST will be able to see objects that are farther away from Earth and from farther back in the history of the universe.Since it will be able to capture images from farther back in time,the goal of the JWST will be to investigate the formation of stars and galaxies,as well as what happened after the formation of the universe.Specifically,it has four main science themes. The theme"First Light"will be investigating the formation of the first bright objects in the universe. "Assembly of Galaxies"will explore how galaxies and dark matter formed.The other two science themes are "The Birth of Stars and Protoplanetary Systems,"and"Planetary Systems and the Origins of Life."If the next"Earth 2.0"hasn't been found before JWST is launched,it will surely find it.Having a greater sensitivity at IR wavelengths than Hubble and Spitzer will allow the JWST to measure the transmission and emission spectra of planets more precisely.Measuring the amount of light at specific wavelengths that gets absorbed by a planet's atmosphere as it passes in front of a star in its line of sight can give insight on the chemical composition of the atmosphere.By specifically targeting IR wavelengths the JWST will be able to detect the presence of H20,CH4,CO and CO2 in the atmospheres of planets outside the solar system. L3 5 The James Webb Space Telescope will orbit the L2 Lagrange point.(wst.nasa.gov)
Our study of Lagrange points can be applied to the James Webb Space Telescope (JWST) since it will be placed at the EarthSun L2 point. JWST is set to launch in 2018 and will be the successor to the Hubble Space Telescope and Spitzer Space Telescope. Since 1990, Hubble has been capturing light from the UV, visible, and nearIR, covering wavelengths from 0.1 to 1.7 microns. JWST will aim to improve previous work by acquiring more precise measurements in the IR spanning 0.6 to 29 microns. As light from stars and galaxies travels through space, its wavelength increases. Therefore, by capturing light with longerwavelengths, JWST will be able to see objects that are farther away from Earth and from farther back in the history of the universe. Since it will be able to capture images from farther back in time, the goal of the JWST will be to investigate the formation of stars and galaxies, as well as what happened after the formation of the universe. Specifically, it has four main science themes. The theme “First Light” will be investigating the formation of the first bright objects in the universe. “Assembly of Galaxies” will explore how galaxies and dark matter formed. The other two science themes are “The Birth of Stars and Protoplanetary Systems,” and “Planetary Systems and the Origins of Life.” If the next “Earth 2.0” hasn’t been found before JWST is launched, it will surely find it. Having a greater sensitivity at IR wavelengths than Hubble and Spitzer will allow the JWST to measure the transmission and emission spectra of planets more precisely. Measuring the amount of light at specific wavelengths that gets absorbed by a planet’s atmosphere as it passes in front of a star in its line of sight can give insight on the chemical composition of the atmosphere. By specifically targeting IR wavelengths the JWST will be able to detect the presence of H20, CH4, CO and CO2 in the atmospheres of planets outside the solar system. The James Webb Space Telescope will orbit the L2 Lagrange point. (jwst.nasa.gov)
Currently,the Hubble Space Telescope orbits the Earth at a distance of about 570 kilometers. JWST,however,will orbit the L2 Lagrange point,about 1.5 million kilometers away,as shown in the figure above.As the Earth orbits the Sun,JWST will orbit with it,maintaining a constant position relative to the Sun and Earth.In addition to stability,this position has several advantages for an IR telescope.Its constant relative position means the Earth and Moon will never obscure the telescope's view.It is important for an IR telescope to stay cool(below 50 K on the instrument side of the telescope,as shown in the figure below),so the JWST makes use of a solar shield,essentially a large sheet of reflective material.Since the light from the Sun will come only from one direction,the solar shield can effectively block it. The Two Sides of the Webb Telescope Hot side Cold side 185Fahrenheit 388°Fahrenheit (85°Celsius) (-233°Celsius) Solar panel Science instruments: Communications Detectors filters antenna Mimrors Computer Steering: Reaction wheels&jets light from the sun At the L2 Lagrange point,the solar shield can block the light and heat from the Sun. (jwst.nasa.gov) In the second section,the stability of the Lagrange points will be analyzed analytically.The third section will show the results of the numerical model for the stability,and the timescales on which these are valid.The fourth section will outline the consequences of these results for the James Webb Space Telescope. Il.Analytic Stability of Lagrange Points Having previously solved for the location of the five Lagrange points in the Earth-Sun system, our next task is to evaluate the stability of each.Since force is the negative gradient of a scalar potential, usually one can solve for stability by analyzing the shape of the effective potential at the point of interest. In our case,however,this method of analysis does not work because the coriolis force leads to a velocity dependent potential.Therefore,in order to analyze the stability we must linearize the equations of motion about each Lagrange point,and solve for small perturbations from equilibrium.In other words,we begin with:
Currently, the Hubble Space Telescope orbits the Earth at a distance of about 570 kilometers. JWST, however, will orbit the L2 Lagrange point, about 1.5 million kilometers away, as shown in the figure above. As the Earth orbits the Sun, JWST will orbit with it, maintaining a constant position relative to the Sun and Earth. In addition to stability, this position has several advantages for an IR telescope. Its constant relative position means the Earth and Moon will never obscure the telescope’s view. It is important for an IR telescope to stay cool (below 50 K on the instrument side of the telescope, as shown in the figure below), so the JWST makes use of a solar shield, essentially a large sheet of reflective material. Since the light from the Sun will come only from one direction, the solar shield can effectively block it. At the L2 Lagrange point, the solar shield can block the light and heat from the Sun. (jwst.nasa.gov) In the second section, the stability of the Lagrange points will be analyzed analytically. The third section will show the results of the numerical model for the stability, and the timescales on which these are valid. The fourth section will outline the consequences of these results for the James Webb Space Telescope. II. Analytic Stability of Lagrange Points Having previously solved for the location of the five Lagrange points in the EarthSun system, our next task is to evaluate the stability of each. Since force is the negative gradient of a scalar potential, usually one can solve for stability by analyzing the shape of the effective potential at the point of interest. In our case, however, this method of analysis does not work because the coriolis force leads to a velocity dependent potential. Therefore, in order to analyze the stability we must linearize the equations of motion about each Lagrange point, and solve for small perturbations from equilibrium. In other words, we begin with:
x=xi+6x,Ux=OVx, x=Vi+oy,vy=ovy, Since we are evaluating the stability about each Lagrange point,(x)refers to the coordinates of the Lagrange point whose stability we wish to calculate.For small perturbations,the linearized equation of motion is as follows: 0 0 0 0 0 d dVa dUa dr2 0 251 血的 6 PUa U dy? -20 0 dudr It is important to note the the lower right block is due to the velocity dependent coriolis force, which did not make an appearance in solving for the Lagrange points because at that time we were solving for locations where the velocity was zero.With this in mind,we realize that the x and y component of the original inertial force are no more than the x and y derivative of the effective potential. Thus,we only need to take the x and y derivative of these components to solve the lower left block. With an analytic expression for the linearized evolution matrix,the next step is to evaluate the matrix at each Lagrange point,then use the eigenvalues of the equation to describe the stability. Using the previously derived equation for the inertial force, Fo=02(-(ctaR)R a(-R)R (x+aR)2+y2) e-月 +22(y- c-0+vg月 ayR3 (x+aR2+y2) we take the partial derivatives with respect to x and y for the x and y components to obtain expressions for the lower left block of the linearized evolution matrix 品 Q2(x-(EtaRR a(r-R)R +o2+9G-m24) =22 BR3 38R3(aR+z)2 aR+2P+y)7十(aR+品+n- aR3 3aR3(z-BR)2 2-8P+vn+(-9+n+1
Since we are evaluating the stability about each Lagrange point, (x , )refers to the coordinates i yi of the Lagrange point whose stability we wish to calculate. For small perturbations, the linearized equation of motion is as follows: It is important to note the the lower right block is due to the velocity dependent coriolis force, which did not make an appearance in solving for the Lagrange points because at that time we were solving for locations where the velocity was zero. With this in mind, we realize that the x and y component of the original inertial force are no more than the x and y derivative of the effective potential. Thus, we only need to take the x and y derivative of these components to solve the lower left block. With an analytic expression for the linearized evolution matrix, the next step is to evaluate the matrix at each Lagrange point, then use the eigenvalues of the equation to describe the stability. Using the previously derived equation for the inertial force, we take the partial derivatives with respect to x and y for the x and y components to obtain expressions for the lower left block of the linearized evolution matrix
82Ug 8r8y 22(x- (rtaR)R3 a(z-R)R3 (e+a)2+y2) (e-R)2+y2)8 =22 3Ry B(aR+x】 a(z-BR) (aR+)2+y2)2+ (x-BR)2+y2)572 82Ug 22(y yR ayR3 8y8x (z+aR)2+y2)号 (x-)2+y2)2 =22 3Ry B(aR+E) a(E=BR】 (aR+z)2+y2)572+ (z-8)2+y2)572 82Va n2(y- R ayR 0w2 By (e+aR)2+y2) =22 8 33R3y2 3aRy2 (aR+)2+y2)3+ (aR+x)2+y2572- 2-82+y2+e-97+1 Evaluating these expressions at the symmetrical points LI and L2,we obtain the following results: L1: L2: 82Uua=-922 8r2 82U=922 8x2 82Un=0 Bxoy 82U=0 Bxay 82=0 8yox 82Ug=0 bybx 20=302 82=-322 0y2 Therefore,the eigenvalues of the linearized evolution matrix are: 入1,2=±V1+2W7 Due to the symmetry of LI and L2,it is to be expected that the eigenvalues for each are the same.Since one eigenvalue is real and positive and the other real and negative,we conclude that LI and L2 are saddle points.They are not dynamically stable,meaning that small deviations from the exact equilibrium point will grow over time.Without corrections,an object perturbed slightly from LI or L2 would drift away from the stable equilibrium Next,we analyze the stability at L3
Evaluating these expressions at the symmetrical points L1 and L2, we obtain the following results: Therefore, the eigenvalues of the linearized evolution matrix are: Due to the symmetry of L1 and L2, it is to be expected that the eigenvalues for each are the same. Since one eigenvalue is real and positive and the other real and negative, we conclude that L1 and L2 are saddle points. They are not dynamically stable, meaning that small deviations from the exact equilibrium point will grow over time. Without corrections, an object perturbed slightly from L1 or L2 would drift away from the stable equilibrium. Next, we analyze the stability at L3
L3: 82=-3n2 0x2 828=0 8xOy 82U9=0 By8: 9=胎2 The eigenvalues from the linearized evolution matrix are: 入1,2=0√器 Here a and B refer to the mass ratios of the Earth and Sun,respectively: ME Q=ME+M5 BMMs Since the Earth has a mass far smaller than that of the Sun,the eigenvalues are small numbers, approximately.625106.As for LI and L2,evaluating the evolution matrix for L3 yields one real and positive eigenvalue and one real and negative eigenvalue,indicating a saddle point.It is similarly dynamically unstable.However,since the eigenvalue is much smaller in magnitude than those for LI and L2,the instability is less dramatic.An object displaced slightly from L3 willd drift away from equilibrium much more slowly than it would if displaced the same amount from LI or L2 Analysis of L4 and L5 is slightly more complicated,as none of the partial derivatives cancel to zero.Since L4 and L5 are symmetric,the entries in the matrix are similar. L4: L5: 2=n2 -n2 =3¥(3-a)2 8xOy g=3Y(a-8)n2 8rOy 89=3¥(8-a)22 8yOE 80=3¥(a-9)22 0yox 歌=n2 50=n2 The eigenvalues,then,are:
The eigenvalues from the linearized evolution matrix are: Here α and β refer to the mass ratios of the Earth and Sun, respectively: Since the Earth has a mass far smaller than that of the Sun, the eigenvalues are small numbers, approximately . As for L1 and L2, evaluating the evolution matrix for L3 yields one real and positive eigenvalue and one real and negative eigenvalue, indicating a saddle point. It is similarly dynamically unstable. However, since the eigenvalue is much smaller in magnitude than those for L1 and L2, the instability is less dramatic. An object displaced slightly from L3 willd drift away from equilibrium much more slowly than it would if displaced the same amount from L1 or L2. Analysis of L4 and L5 is slightly more complicated, as none of the partial derivatives cancel to zero. Since L4 and L5 are symmetric, the entries in the matrix are similar. The eigenvalues, then, are:
d1,2=士i号V2-√27B-a2-28 Because the mass of the Sun is much greater than the mass of the Earth,B is much greater than a,and the eigenvalues are all imaginary,with no real components.This result indicates that L4 and L5 are stable. It is interesting to note that evaluation of only the bottom left 2 x 2 portion of the matrix, involving the effective potential,yields two positive eigenvalues.This would indicate an unstable node However,because the system is rotating and not stationary,the coriolis force affects objects near L4 and L5,pulling them toward the point if they are perturbed.Therefore,L4 and L5 are actually stable equilibrium points. III.NUMERICAL SIMULATION TESTING STABILITY In order to verify our analytical solution we created a 3D N-body code to assess the stability regions of our solution.Our code is written using a combination of Python and C/C++.All of the plotting and setting up of initial conditions was done in Python and we outsource the computational work to C++.To conserve energy and accuracy in our code we invoke a fourth order method for updating the position and velocity.The fourth order method is formally known as the fourth order Hermite which is an implicit method that uses a predictor-corrector scheme.Please reference our midterm report for the exact computational formulation we adopt. We tested two different models of the Earth-Sun system to assess the stability of each of the Lagrange points.The first model we tested was one where the initial positions of the bodies placed at the Lagrange points were not perturbed.See Figure 1.Using this model we can compute the stability of each of the points by analyzing the drift timescale or e-folding time.Small departures from equilibrium will grow exponentially with an e-folding time of roughly 1/eigenvalue for the associated Lagrange point. By fitting a function of the form A*exp(B*t)to the distance of a body from its initial point as a function of time we can obtain the drift timescale as 1/B.Below is a table containing the drift times for the 3 unstable Lagrange points.The optimal function was found using the Levenberg-Marquardt minimization algorithm which uses a damped least squares method to solve nonlinear least squares problems
Because the mass of the Sun is much greater than the mass of the Earth, β is much greater than α, and the eigenvalues are all imaginary, with no real components. This result indicates that L4 and L5 are stable. It is interesting to note that evaluation of only the bottom left 2 x 2 portion of the matrix, involving the effective potential, yields two positive eigenvalues. This would indicate an unstable node. However, because the system is rotating and not stationary, the coriolis force affects objects near L4 and L5, pulling them toward the point if they are perturbed. Therefore, L4 and L5 are actually stable equilibrium points. III. NUMERICAL SIMULATION TESTING STABILITY In order to verify our analytical solution we created a 3D Nbody code to assess the stability regions of our solution. Our code is written using a combination of Python and C/C++. All of the plotting and setting up of initial conditions was done in Python and we outsource the computational work to C++. To conserve energy and accuracy in our code we invoke a fourth order method for updating the position and velocity. The fourth order method is formally known as the fourth order Hermite which is an implicit method that uses a predictorcorrector scheme. Please reference our midterm report for the exact computational formulation we adopt. We tested two different models of the EarthSun system to assess the stability of each of the Lagrange points. The first model we tested was one where the initial positions of the bodies placed at the Lagrange points were not perturbed. See Figure 1. Using this model we can compute the stability of each of the points by analyzing the drift timescale or efolding time. Small departures from equilibrium will grow exponentially with an efolding time of roughly 1/eigenvalue for the associated Lagrange point. By fitting a function of the form A*exp(B*t) to the distance of a body from its initial point as a function of time we can obtain the drift timescale as 1/B. Below is a table containing the drift times for the 3 unstable Lagrange points. The optimal function was found using the LevenbergMarquardt minimization algorithm which uses a damped least squares method to solve nonlinear least squares problems
Point Drift Time (Computational) Drift Time (Analytical) LI 24.9 Days -23 Days L2 24.4 Days ~23 Days L3 43.12 Years ~150 Years From this result our computationally derived drift time more or less coincides with our analytically derived answer.You can see that a body placed at LI or L2 will wander off after a few months unless some corrections are made.However if you have a body at L3 it will need far less adjustment compared to LI or L2 to maintain a steady orbit.The reason being is that the L3 point is not influenced by the gravitational potential of the Earth as much as LI or L2. In the second model each body was randomly perturbed using a gaussian distribution centered on the position of the associated Lagrange point and given a variance equal to 0.01%of its position element.See Figure 2 for one of the trial perturbation runs.In this model we tested the boundary of each Lagrange point by perturbing the initial position of the body.Our results are rather hand-wavy since a full Monte Carlo simulation is needed to robustly determine the sensitivity of the boundaries around the Lagrange points.Assessing the energy contour plot is better suited for determining the orbital paths of various bodies.See Figure 3
Point Drift Time (Computational) Drift Time (Analytical) L1 24.9 Days ~23 Days L2 24.4 Days ~23 Days L3 43.12 Years ~150 Years From this result our computationally derived drift time more or less coincides with our analytically derived answer. You can see that a body placed at L1 or L2 will wander off after a few months unless some corrections are made. However if you have a body at L3 it will need far less adjustment compared to L1 or L2 to maintain a steady orbit. The reason being is that the L3 point is not influenced by the gravitational potential of the Earth as much as L1 or L2. In the second model each body was randomly perturbed using a gaussian distribution centered on the position of the associated Lagrange point and given a variance equal to 0.01% of its position element. See Figure 2 for one of the trial perturbation runs. In this model we tested the boundary of each Lagrange point by perturbing the initial position of the body. Our results are rather handwavy since a full Monte Carlo simulation is needed to robustly determine the sensitivity of the boundaries around the Lagrange points. Assessing the energy contour plot is better suited for determining the orbital paths of various bodies. See Figure 3
1.5 1.0 0.5 0.0 -0.51 -1.01 -1.5 -1.0 -0.50.0 0.5 1.0 1.5 Figure 1.Orbits of Lagrange points in a co-rotating frame after 10 years where the initial positions were not perturbed from equilibrium
Figure 1. Orbits of Lagrange points in a corotating frame after 10 years where the initial positions were not perturbed from equilibrium