Midterm/Final Team Projects

Physics 141/241

Winter 2020

 

Dear students:

The link to galactics2020Final.tar.gz has a compiled working Galactics package under the Sierra Mac OS X platform and I also tested it on the HighSierra Mac OS X platform with the same binaries. There are pre-prepared default galaxies and their bulge-disk-halo galaxy components in the appropriate directories. You should be able to prepare new non-default galaxies with "make galaxy" in non-lab enviroment in any of the four galaxy source directories if working online outside the lab. This would provide galaxy preparation outside the lab environment. Running under Linux might require some adjustments in the Makefile. Mac laptops or iMacs should work without modifications.

The link to treecode2020Final.tar.gz has a compiled treecode package under the HighSierra Mac platform. The command "treecode" will run a default test on a Plummer model with 4096 particles. The code is very fast and you can run the final projects successfully with it outside the lab environment. In case you prefer to work online, not using the lab environment, we will not require the GyrfalcON version which only runs directly on the lab iMacs. The link to the treecode manual treeguide_c.pdf describes how to run "treecode" with your custom prepared input data files of the galaxies and also further describes the input-output part in detail. Running under Linux might require some adjustments in the Makefile. Mac laptops or iMacs should work without modifications.

 

Midterm/Final:

The teams work on the Midterm projects within the framework of the Toomre&Toomre model of the restricted three-body problem. The final projects are the best full simulations of the teams using some bulge-disk-halo galaxy collison model.

 

Midterms will be presented on March 10.

Final projects will be presented on March 20 between 2:00 pm - 5:00 pm.

 

Useful links

 

midterm examples:

 

Mice powerpoint summary

Mice view1

Mice view2

Antennae report sample

Antennae research note

Antennae view1

Antennae view2

 

Andromeda report sample

Cartwheel report sample

Cartwheel massive head-on

Cartwheel massive on-center

Cartwheel massless head-on

Cartwheel massless on-center

 

PhD thesis web site

pdf file of the PhD thesis


===============================================================================

 

The Mice Project (NGC 4676)

Mice Teams

Mice A team: Aguilar, George, Chen

 

Mice B team: Clinton, Hah, Ishkhanyan

 

Link to Mice Galaxy collision

Mice(Hubble)

For convenience, Toomre&Toomre chose the time unit to be MATH, the pericenter distance $R_{\min }=25kpc$ and the heavier mass in each encounter, $10^{11}M_{solar}$, chosen as the mass unit. The collision orbits (elliptic, or parabolic) are identified by choosing (1) $R_{\min }=r_{p}$ and (2) the initial separation $R_{init}$, or the equivalent time $t_{init}$, with $t=0$ chosen at pericenter. Initial velocities are dependent on the masses of the galaxies. Disk orientations and other input parameters are determined from the reference material of the projects as requested in the project descriptions. For the Midterm project Toomre&Toomre disk galaxy models are prepared with test particles to run the restricted 3-body integrator. In the Final project realistic Galaxies are crafted at the origin of the coordinate system using the code of Konrad Kuijken and John Dubinski to generate disk-bulge-halo galaxies. The code is discussed in the paper by Kuijken&Dubinski and the tar file of the code is packaged galactics.tar.gz.

Midterm: Mice Project

Figure 22 on page 657 in Toomre&Toomre summarizes the initial setup of their most successive Mice encounter. Reproduce their result of Fig.22 using 297 disk particles in each disk as described on page 20. They fill 11 concentric rings of radii MATH. Increase the filling of the rings by increments of three test particles and starting with $12$ particles on the innermost ring. The outermost ring has $42$ particles. Run your restricted 3-body code where the test particles are tracking the two elliptic orbits of the core galaxies. The test particles are interacting with the cores only. Work in the center-of-mass coordinate system where the relative coordinate vector $\overrightarrow{r}$ traces the elliptic Kepler orbit. For the elliptic orbits, use the appropriate parametrization as it was outlined in Lecture 5:


Figure

The origin of the CM coordinate system cn be chosen to coincide with the Sun location of the figure. Two equal disks of radius $0.7R_{\min }$ experience an elliptic encounter with $e=0.6$, having begun flat and circular at the time $t=-16.4$ of the last apocenter. As viewed from either disk, the adopted node to peri angles MATH with inclinations MATH. The angles are discussed in Fig. 6a,b (see also Galaxy Rotations note and collision orbit note). The resulting composite object is shown at $t=6.086$ in Fig. 22.

(1) Compare Toomre&Toomre to your simulation.

(2) Animate and make an Mpeg movie.

(3) Give a simple explanation of the tidal tails and explain whether the tidal tails are in the plane of the rotating disks, or the collision plane, or neither.

(4) Write a project report and prepare a 20 minutes presentation.




Final: Project Mice

Barnes in his paper argued that a realistic full simulation of the encounter using realistic bulge-disk-halo galaxy model calls for parabolic encounter. He gave arguments that the bulge-disk-halo galaxy structure requires to replace the elliptic encounter of Toomre&Toomre with parabolic encounter. For the initial condition, consider to separate galaxies on parabolic collision course. Galaxy A and galaxy B do not overlap initially, and the center of masses of the individual galaxies would continue on a parabolic orbit, if the interactions of the stars would not change the picture when the galaxies begin to overlap.


mice_top.gif

First, we would like to set up the appropriate initial condition for the collision as depicted in the figure above. Notes on the setup are given in collision orbit note. Notes on the disk angles are given in Galaxy Rotations note. Barnes chose in the paper the length unit to be $35.7kpc$, the velocity unit as MATH, and the mass unit as MATH. One time unit then works out to be $194Myr$. Barnes chose for the galaxy components: MATH Parabolic orbit was chosen with $r_{p}=0.25,$ $r_{disk}=1/12,$ $a_{b}=0.04168,$ $z_{disk}=0.005,$ $a_{halo}=0.1$. The viewing time of the encounter is $t=0.875$ after pericenter. The disk angles are MATH.

(1) Use the Galactics package to fabricate the initial galaxies. Choose your resolution for the project.

(2) Run Gyrfalcon to reproduce the results of Barnes in the article.

(3) Animate and make Mpeg movies of the project.

(4) Give again simple explanation of the consistency of tidal tails with Toomre&Toomre and explain again whether the tidal tails are in the plane of the rotating disks, or the collision plane, or neither.

(5) Write a project report and prepare a 20 minutes presentation.


===============================================================================

 

The Antennae Project (NGC 4038/9)

Antennae Teams

 

Antennae A team: Jiang, Johnson

Antennae B team: Stanfield, Suharnoko, Zeng

 

Link to Antennae Galaxy collision

Antennae(Hubble)

 

 

For convenience, Toomre&Toomre chose the time unit to be MATH, the pericenter distance $R_{\min }=25kpc$ and the heavier mass in each encounter, $10^{11}M_{solar},$ chosen as the unit mass. The collision orbits (elliptic, or parabolic) are identified by choosing (1) $R_{\min }=r_{p}$ and (2) the initial separation $R_{init}$, or the equivalent time $t_{init}$, with $t=0$ chosen at pericenter. Initial velocities are dependent on the masses of the galaxies. Disk orientations and other input parameters are determined from the reference material of the projects as requested in the project descriptions. For the Midterm project Toomre&Toomre disk galaxy models are prepared with test particles to run the restricted 3-body integrator. In the Final project realistic Galaxies are crafted at the origin of the coordinate system using the code of Konrad Kuijken and John Dubinski to generate disk-bulge-halo galaxies. The code is discussed in the paper by Kuijken&Dubinski and the tar file of the code is packaged Galactics.tar.gz.

Midterm: Antennae Project

Figure 23 on page 37 in Toomre&Toomre summarizes the initial setup of their symmetric Antennae encounter. Reproduce their result of Fig.23 using 345 disk particles in each disk as described on page 20. They fill 12 concentric rings of radii MATH. Increase the filling of the rings by increments of three test particles added to the next ring and starting with $12$ particles on the innermost ring. The outermost ring has $45$ particles. Run your restricted 3-body code where the test particles are tracking the two elliptic orbits of the core galaxies with $e=0.5$. The test particles are interacting with the cores only and the interaction is softened by MATH at the closest range to mimic the distributed mass of the real object. Work in the center-of-mass coordinate system where the relative coordinate vector $\overrightarrow{r}$ traces the elliptic Kepler orbit. For the elliptic orbits, use the appropriate parametrization as it was outlined in Lecture 5:


Figure

The origin of the CM coordinate system can be chosen to coincide with the Sun location of the figure. Two equal disks of radius $0.75R_{\min }$ experience an elliptic encounter, having begun flat and circular at the time $t=-16.4$ of the last apocenter. As viewed from either disk, the adopted node to peri angles MATH with inclinations MATH. The angles are discussed in Fig. 6a,b (see also Galaxy Rotations note and collision orbit note). The viewing time of the resulting composite object is $t=15$ in Fig. 23.

(1) Compare Toomre&Toomre to your simulation.

(2) Animate and make an Mpeg movie.

(3) Give a simple explanation of the tidal tails and explain whether the tidal tails are in the plane of the rotating disks, or the collision plane, or neither.

(4) Write a project report and prepare a 20 minutes presentation.




Final: Antennae Project

Barnes in his paper on the Antennae designed a realistic full simulation of the encounter using realistic bulge-disk-halo galaxy model with elliptic encounter. For the initial condition, consider to separate galaxies on elliptic collision course with $e=0.5$. Galaxy A and galaxy B do not overlap initially, and the center of masses of the individual galaxies would continue on an elliptic orbit, if the interactions of the stars would not change the picture when the galaxies begin to overlap.


mice_top.gif

First, we would like to set up the appropriate initial condition for the collision as depicted in the figure above, but the parbolic orbits replaced by elliptic ones. Notes on the setup for parbolic orbits are given in collision orbit note and have to be replaced by the properties of elliptic orbits. Notes on the disk angles are given in Galaxy Rotations note. Barnes chose in the paper the length unit to be $40kpc$, and the mass unit as MATH. One time unit is $250$ $Myr$. Barnes chose for the galaxy components: MATH with a total mass MATH Elliptic orbits were chosen and started with $r_{p}=0.5,$ $r_{disk}=1/12,$ node to peri angles MATH with inclinations MATH and time $t=1$ to epicenter to produce a slow and symmetrically prograde encounter with the two disks inclined so as to sling tidal tails high above the orbital plane where they will eventually be seen in projection as crossing each other. The viewing time of the encounter is close to the nex apocenter.

(1) Use the Galactics package to fabricate the initial galaxies. Choose your resolution for the project.

(2) Run Gyrfalcon to reproduce the results of Barnes in the article.

(3) Animate and make Mpeg movies of the project.

(4) Give again simple explanation of the consistency of tidal tails with Toomre&Toomre and explain again whether the tidal tails are in the plane of the rotating disks, or the collision plane, or neither.

(5) Write a project report and prepare a 20 minutes presentation.


================================================================================

 

 

Cartwheel Project

Cartwheel Teams

Cartwheel A team:Guiang, Fajardo, Arora

 

Cartwheel B team: McElwee, Wheeler, Nguyen

 

Link to Cartwheel Galaxy collision

Cartwheel(Hubble)

 

 

 

Midterm: Cartwheel

In the spirit of the Toomre&Toomre modeling of the Mice, Antennae, and Milky Way mergers, study the paper by Lynds&Toomre to set up and reproduce their model.

(1) Obtain Lynds&Toomre models mimicking the mergers of Figs. 5 and 6 in the paper.

(2) Animate and make Mpeg movies of your mergers.

(3) Give a simple explanation of the Cartwheel.

(4) Write a project report and prepare a 20 minutes presentation.




Final: Cartwheel

In the spirit of realistic modeling of galaxy mergers, like the Mice, Antennae, and Milky Way mergers, study the paper by Athanassoula et al. to set up and reproduce their model of the Cartwheel using bulge-disk-halo galaxies.

(1) Use the Galactics package to fabricate the initial galaxies. Choose your resolution for the project.

(2) Run Gyrfalcon to get close to the results of Athanassoula et al. as reported in their paper. Consider also the Aarseth code, if it suits better the project, in case of demand on accuracy.

(3) Animate and make Mpeg movies of the project.

(4) Give again simple explanation of how the Cartwheel merger happened.

(5) Write a project report and prepare a 20 minutes presentation.

=============================================================================

 

 

Andromeda Project

Andromeda Team

Andromeda A team: Kennedy, Kambalur

 

 

Links to Andromeda Galaxy collision

 

 

 

Andromeda (Hubble)

 

 

Journal publication on Andromeda collision

 

Andromeda collision: Master Thesis I.

 

 

Andromeda collision: Master Thesis II.

 

 

 

=============================================================================