Make a file of new moons over the past 524 years

Start in January 1502 and end in December 2025. That works out to be 16 of the Mayans’ cycles of 405 synodic months.
In[]:=
newmoons={MoonPhaseDate[DateObject[{1502,1,1}],"NewMoon"]}
Out[]=

Sat 18 Jan 1502 18:53:37GMT-6

In[]:=
n=405*16
Out[]=
6480
In[]:=
For[i=0,i<n,i++,AppendTo[newmoons,MoonPhaseDate[DatePlus[Last[newmoons],1],"NewMoon"]]]
In[]:=
​​Last[newmoons]
Out[]=
Fri 19 Dec 2025 19:43:20GMT-6
Save the list of new moons to a file so later calculations aren’t delayed by making this list.
In[]:=
Save[FileNameJoin[NotebookDirectory[],"newmoons.wl"],newmoons]