Midterm and Final Team Projects



Physics 141/241

Winter 2017


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 the bulge-disk-halo Dubinski model for their galaxy encounter projects.

 

Progress on first phase of projects will be presented on March 13 and March 15.

Final projects will be presented on March 24 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:Alolayan,Bosey,Li,Ye

 

Mice B team:Aldas,An,Lund,Lewis

 

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:Whiteaker,Wei,Mendoza,Delight

 

Antennae B team:Kwong,Guerro,Park,Shakouri

 

 

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:Au,Chen,Howland,Moody

 

Cartwheel B team:Pappas,Wu,Wylie-Modro,Zhai

 

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 Teams

Andromeda A team:Chandrasekar-Iyer,Kumar,Nayak,Singh

 

Andromeda B team:

 

Links to Andromeda Galaxy collision

 

 

 

Andromeda (Hubble)

 

 

Final paper on Andromeda collision

 

Journal publication on Andromeda collision

 

Andromeda collision: Master Thesis

 

 

 

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

 



CUDA software/hardware support:

 

ACCESS to the igpu6-210.ucsd.edu server with 8 Titan cards:

The Titan box is now up and running as hostname 'igpu6-210.ucsd.edu',
and students in ph141/ph241 may login with their course credentials. It can be accessed through ssh.

Students should have access to their class storage directories.

Fast local file storage is in the form of an 80gb /scratch/ filesystem
located on the local SSD drives. This data is not backed up.

While igpu6-210 is configured as 64-bit, it mounts some 32-bit software
under /usr/local, including an updated version of gcc. Before compiling
any of the nvidia examples you'll need to switch back to the system gcc:

[agt@igpu6-210]:~:500$ which gcc
/software/common/gcc-4.8.1/bin/gcc
[agt@igpu6-210]:~:501$ module unload gcc
[agt@igpu6-210]:~:502$ which gcc
/usr/bin/gcc
[agt@igpu6-210]:~:503$ which g++
/usr/bin/g++

ACS ran into trouble getting the virtualized access to the Titan cards
functioning reliably, so for now, we've just got one timeshared Linux
system hosting all 8 cards, running on the bare hardware.

Down the road we anticipate we'll want to install "cuda_wrapper" which
allows cards to be assigned to specific users and/or groups.

CUDA nbody code on the lab iMac:

Path to nbody source code directory: /Developer/NVIDIA/CUDA-8.0/samples/5_Simulations/nbody

use the command make to compile the code

binary (executable) runs from: /Developer/NVIDIA/CUDA-8.0/samples/bin/x86_64/darwin/release/nbody

nbody CUDA paper from NVIDIA

CUDA TreeCodes:

BONSAI CUDA code

BONSAI CUDA paper

Bonsai student project (courtesy of Michael Folkert)

Bonsai student project (zip file)

 

OctGrav CUDA code

OctGrav CUDA paper

CUDA exact Nbody Code:

Sapporo CUDA library code

Sapporo CUDA paper