  Within a simulated Milky Way-like galaxy, we identify and analyse analogues
of the Gaia-Enceladus (GE), Kraken and Sequoia mergers that each matches
remarkably well observational results, including in velocity and chemical
abundance space, and their distributions in the $j_{z}$-Energy plane. The
Kraken analogue is the earliest merger and has the highest total mass ratio.
Consistent with previous studies, it is chemically indistinguishable from old
in-situ stars at the time of its accretion. The GE and Sequoia analogue events
accrete at similar times in our simulation, both along filaments but from
opposite sides of the main galaxy. The mean stellar ages of the GE and Sequoia
analogues are both similar and, from our simulation results, we see that they
can be separate entities and still naturally reproduce the observed properties
of their stellar remnants at the present day, including the significant
retrograde velocities of the Sequoia analogue remnant stars and the difference
in the tracks of the two galaxies through chemical abundance space. Our results
provide supporting information about the properties of these three merger
events, and show for the first time that they can all be reproduced with a
fully cosmological simulation, providing a possible self consistent
evolutionary pathway for the Milky Way's formation.

Astrophysics of Galaxies