<p>I was able to see this using skygazer 4.5M, quite frankly, I am new to astronomy and inspired by other's (mainly Nilesh Oak) work I started doing this.</p>
<p>He had already proved this and I wanted to see if everything that he said is true. writing code to simulate this in astropy is just to much work and I am not that familiar with the package so, for now I cannot do anything about it(but I have not quit on this yet)</p>
<p>For this to work we need to consider the change in pole star, from Polaris to Vega which happens every 13000 years(half of full precession of equinoxes).</p>
<p>Now again, since we are looking at the sky with naked eye and it is from earth, we can only say how it might have looked from earth at specific latitude and longitude, so if we consider all this and movement of stars and the movement of earth itself, it will look like the star Alcor is crossing the meridian first before Mizar, really fascinating observation.</p>
<p>If you are interested, here is the link to the video where Nilesh Oak, talks about it:
<div class="youtube-embed"><div>
                <iframe width="640px" height="395px" src="https://www.youtube.com/embed/RedV48OCEFg?start=0"></iframe>
            </div></div></p>
<p>Mizar is known as Vasishtha and Alcor as Arundhati in the video.</p>
<p>Both of the answers has points that helped me I will give up votes for both and accept mine as answer to close the question.</p>

CC BY-SA 4.0
2021-10-24T11:00:20.567Z
2021-10-24T16:20:11.153Z
2021-10-24T16:20:11.153Z
Prashanth Kumar B
<p>The projected separation of Mizar and Alcor is 12 arcminutes. If they are both at a distance of 25 pc (suggested by the Hipparcos satellite), then this equates to a physical separation of at least 18,000 au. Using Kepler's third law we can then estimate an orbital period of <em>at least</em> 800,000 years (assuming a total system mass of around <span class="math-container">$9M_\odot$</span>. Thus over a period of time of 13,000 years or so I think you can safely neglect the acceleration caused by the interaction of Mizar and Alcor.</p>
<p>Similarly over the course of 13,000 years, the distance to the stars will only change by about 0.0013 pc per km/s of radial velocity. The radial velocities listed in SIMBAD are <span class="math-container">$-6.3 \pm 0.4$</span> and <span class="math-container">$-8.9 \pm 0.9$</span> km/s for Mizar and Alcor respectively, and proper motions in RA of <span class="math-container">$119.0 \pm 1.5$</span> and <span class="math-container">$120.21 \pm 0.12$</span> mas/year and in Declination of <span class="math-container">$-26.0 \pm 1.7$</span> and <span class="math-container">$-16.04 \pm 0.14$</span> mas/year. Assuming the velocity tangential to the Sun remains constant then these low radial velocities will not change the measured proper motions very significantly over that time and certainly not by more than the current uncertainties in the relative distance to the two stars or their proper motions.</p>
<p>Thus all you need to do then to work out (relative) positions on the sky is extrapolate the current positions and proper motions backwards for the requisite amount of time assuming that they they do not change. How accurate these predictions are will depend on how accurate the proper motions are.</p>
<p>Roughly speaking (and you should do this properly using the astropy apply_space_motion routines) the change in Alcor's position <em>relative</em> to Mizar is about <span class="math-container">$1.2 \pm 1.5$</span> mas/year in RA and <span class="math-container">$10.0 \pm 1.7$</span> mas/year in Dec. One simply multiplies these numbers (and their uncertainties) by -10 to get the change of position of Alcor relative to Mizar 10,000 years ago in arcseconds. Given that the current separation is 12 arc<em>minutes</em>, then there wouldn't really be much change in their <em>relative</em> positions visible to the naked eye (I think they would get about 1.6 arcminutes closer together in Declination).</p>

CC BY-SA 4.0
2021-10-20T16:59:23.010Z
2021-10-20T17:22:39.160Z
2021-10-20T17:22:39.160Z
ProfRob
<p>The answer probably depends on what you want, what is 'long time' for you. Astropy offers a method to handle coordinates and proper motion: <a href="https://docs.astropy.org/en/stable/coordinates/apply_space_motion.html" rel="nofollow noreferrer">https://docs.astropy.org/en/stable/coordinates/apply_space_motion.html</a>.</p>
<p>However, unless you have actual positional information for times you are interested in, you have to forward or backward integrate from what you know and basically have to do what you always have to do in such case: obtain the 3D positional information and 3D velocity information and then integrate forward or backward in time.</p>
<p>In practice you probably will get the proper motion of stars in right ascension and declination and the 3rd velocity coordinate along the line of sight is hard to come by and might not be available at all; similarily and additionally the distance might be the least well-constrained part of the stellar coordinates which adds uncertainty to the conversion of proper motion to 3D velocity. You might in many cases get acceptable results, if you ignore the radial component in both spatial and velocity coordinate and keep your integration time small compared to the galactic year (revolution time of the sun around the galactic center, ~220 Myr).</p>
<p>In order to simulate the view of the night sky from Earth 100s or 1000s of years ago or in the future you will also want to simulate the precession of the Earth's rotational axis which will change the <em>apparent</em> position of the stars in the sky by changing the point on the celestial sphere the sky seems to revolve around.</p>

CC BY-SA 4.0
2021-10-20T14:35:57.587Z
What I am trying to do is simulating the night sky view from earth, yes, I understand the earth's postion also becomes a factor here, this is what I wanted to know, how to introduce this, or does astropy handles this....not sure....!

I have update the initial question as well, added further details, about what I am trying to do.
Prashanth Kumar B
CC BY-SA 4.0
2021-10-20T13:45:54.007Z
2021-10-20T13:45:54.007Z
planetmaker
<p>Is there any particular methods that help us calculate the position of a star at a certain point in time? This <a href="https://astronomy.stackexchange.com/questions/38823/how-can-i-calculate-the-apparent-position-of-stars-as-seen-from-a-location-on-ea">question</a> is what I am trying to achieve, but for me, I do not want to use the data because the data that we can access to calculate the position is limited to only a few years. I want a method which uses the current position, speed, velocity of stars and calculate its position for many years back (like in 11000 BCE to 13000 BCE). I am trying to use astropy and I am really new to this, so I just need a push in the right direction; helpful resources will serve a great deal as well, thank you.</p>
<p>Edited:
There are a few stars' positions that are recorded and mentioned in one of our history books. SO what I am trying to do is collect all the recorded stars' positions and try to map them in real-time at a particular year and see if the recorded data exists or not.</p>
<p>Now out of this data, the one that stands out as puzzling as is was the movements of Alcor and Mizar in the sky.</p>
<p>If we observe these 2 stars from earth, for a naked eye it would seem like Mizar is always in front and Alcor is always in following Mizar.</p>
<p>But in one of our history books, there is a recorded event where Mizar follows Alcor (it looks like it is following, I know they just revolve around each other), so what I am trying, is to see if there was any particular period where this would be possible.</p>
<p>This is the starting problem there are many observations like this so it would be really helpful if I could get any info on achieving this.</p>

CC BY-SA 4.0
2021-10-20T13:08:45.543Z
Roughly how many years is "many years back"? 10? 100? 1 000? 10 000?
notovny
CC BY-SA 4.0
2021-10-20T14:24:26.897Z
Should have been clear, many years in the sense back to 11000 BCE to 13000 BCE.
Prashanth Kumar B
CC BY-SA 4.0
2021-10-20T16:41:06.970Z
Can you be clear what you mean by "positions".
ProfRob
<star><positional-astronomy><astropy><skyfield>
Calculating the position of stars at a given point in time
Prashanth Kumar B
# Calculating the position of stars at a given point in time

(Asked by: username_0 on 2021-10-20)

Is there any particular methods that help us calculate the position of a star at a certain point in time? This [question](https://astronomy.stackexchange.com/questions/38823/how-can-i-calculate-the-apparent-position-of-stars-as-seen-from-a-location-on-ea) is what I am trying to achieve, but for me, I do not want to use the data because the data that we can access to calculate the position is limited to only a few years. I want a method which uses the current position, speed, velocity of stars and calculate its position for many years back (like in 11000 BCE to 13000 BCE). I am trying to use astropy and I am really new to this, so I just need a push in the right direction; helpful resources will serve a great deal as well, thank you.

Edited:
There are a few stars' positions that are recorded and mentioned in one of our history books. SO what I am trying to do is collect all the recorded stars' positions and try to map them in real-time at a particular year and see if the recorded data exists or not.

Now out of this data, the one that stands out as puzzling as is was the movements of Alcor and Mizar in the sky.

If we observe these 2 stars from earth, for a naked eye it would seem like Mizar is always in front and Alcor is always in following Mizar.

But in one of our history books, there is a recorded event where Mizar follows Alcor (it looks like it is following, I know they just revolve around each other), so what I am trying, is to see if there was any particular period where this would be possible.

This is the starting problem there are many observations like this so it would be really helpful if I could get any info on achieving this.

**Tags:** star, positional-astronomy, astropy, skyfield
> **Comments:**
> **username_1**: Roughly how many years is "many years back"? 10? 100? 1 000? 10 000?  
> **username_0**: Should have been clear, many years in the sense back to 11000 BCE to 13000 BCE.  
> **username_2**: Can you be clear what you mean by "positions".  

---

## 3 answers

### Answer by username_0 on 2021-10-24. Score: 0 | Accepted answer

I was able to see this using skygazer 4.5M, quite frankly, I am new to astronomy and inspired by other's (mainly Nilesh Oak) work I started doing this.

He had already proved this and I wanted to see if everything that he said is true. writing code to simulate this in astropy is just to much work and I am not that familiar with the package so, for now I cannot do anything about it(but I have not quit on this yet)

For this to work we need to consider the change in pole star, from Polaris to Vega which happens every 13000 years(half of full precession of equinoxes).

Now again, since we are looking at the sky with naked eye and it is from earth, we can only say how it might have looked from earth at specific latitude and longitude, so if we consider all this and movement of stars and the movement of earth itself, it will look like the star Alcor is crossing the meridian first before Mizar, really fascinating observation.

If you are interested, here is the link to the video where Nilesh Oak, talks about it:

Mizar is known as Vasishtha and Alcor as Arundhati in the video.

Both of the answers has points that helped me I will give up votes for both and accept mine as answer to close the question.


---

### Answer by username_2 on 2021-10-20. Score: 3

The projected separation of Mizar and Alcor is 12 arcminutes. If they are both at a distance of 25 pc (suggested by the Hipparcos satellite), then this equates to a physical separation of at least 18,000 au. Using Kepler's third law we can then estimate an orbital period of *at least* 800,000 years (assuming a total system mass of around $9M\_\odot$. Thus over a period of time of 13,000 years or so I think you can safely neglect the acceleration caused by the interaction of Mizar and Alcor.

Similarly over the course of 13,000 years, the distance to the stars will only change by about 0.0013 pc per km/s of radial velocity. The radial velocities listed in SIMBAD are $-6.3 \pm 0.4$ and $-8.9 \pm 0.9$ km/s for Mizar and Alcor respectively, and proper motions in RA of $119.0 \pm 1.5$ and $120.21 \pm 0.12$ mas/year and in Declination of $-26.0 \pm 1.7$ and $-16.04 \pm 0.14$ mas/year. Assuming the velocity tangential to the Sun remains constant then these low radial velocities will not change the measured proper motions very significantly over that time and certainly not by more than the current uncertainties in the relative distance to the two stars or their proper motions.

Thus all you need to do then to work out (relative) positions on the sky is extrapolate the current positions and proper motions backwards for the requisite amount of time assuming that they they do not change. How accurate these predictions are will depend on how accurate the proper motions are.

Roughly speaking (and you should do this properly using the astropy apply\_space\_motion routines) the change in Alcor's position *relative* to Mizar is about $1.2 \pm 1.5$ mas/year in RA and $10.0 \pm 1.7$ mas/year in Dec. One simply multiplies these numbers (and their uncertainties) by -10 to get the change of position of Alcor relative to Mizar 10,000 years ago in arcseconds. Given that the current separation is 12 arc*minutes*, then there wouldn't really be much change in their *relative* positions visible to the naked eye (I think they would get about 1.6 arcminutes closer together in Declination).


---

### Answer by username_3 on 2021-10-20. Score: 2

The answer probably depends on what you want, what is 'long time' for you. Astropy offers a method to handle coordinates and proper motion: <https://docs.astropy.org/en/stable/coordinates/apply_space_motion.html>.

However, unless you have actual positional information for times you are interested in, you have to forward or backward integrate from what you know and basically have to do what you always have to do in such case: obtain the 3D positional information and 3D velocity information and then integrate forward or backward in time.

In practice you probably will get the proper motion of stars in right ascension and declination and the 3rd velocity coordinate along the line of sight is hard to come by and might not be available at all; similarily and additionally the distance might be the least well-constrained part of the stellar coordinates which adds uncertainty to the conversion of proper motion to 3D velocity. You might in many cases get acceptable results, if you ignore the radial component in both spatial and velocity coordinate and keep your integration time small compared to the galactic year (revolution time of the sun around the galactic center, ~220 Myr).

In order to simulate the view of the night sky from Earth 100s or 1000s of years ago or in the future you will also want to simulate the precession of the Earth's rotational axis which will change the *apparent* position of the stars in the sky by changing the point on the celestial sphere the sky seems to revolve around.

> **Comments:**
> **username_0**: What I am trying to do is simulating the night sky view from earth, yes, I understand the earth's postion also becomes a factor here, this is what I wanted to know, how to introduce this, or does astropy handles this....not sure....!  I have update the initial question as well, added further details, about what I am trying to do.  

---

