Using hydrodynamical simulations to understand the inflow of gas into the inner ~200 pc of the Milky Way due to a satellite galaxy merger.

Link to paper

CODE: Periodic Orbits in Non-Axisymmetric Potentials

C++ code used in Gallego & Cuadra 2017 to estimate period orbits at 200 pc of the Galactic Center, where the bar potential dominates. These orbits were used as initial conditions for a hydrodynamical simulation with Gadget-2. The method to identify periodic orbits is called "shooting technique". This starts with an assumption for the initial velocity in the x−y phase-space, in which the x-velocity component is zero on the semimajor axis of the orbit. The equations of motion are then integrated with a fourth-order Runge–Kutta method on the reference system of the rotating bar. The initial velocity assumption is then optimized until obtaining a periodic orbit with a desired fractional precision. Download the code here: find_periodic_orbits.cpp