ELUCID. IX. Recovering the Substructures and History of the Coma Cluster

We employ constrained hydrodynamic simulations of the Coma galaxy cluster from the ELUCID project to study its substructures and assembly history. Our simulations accurately reproduce the global properties of Coma, including its position, virial mass, radius, and surrounding large-scale filaments. Using a combined HBT+SKID method, we obtain a total intracluster light (ICL) fraction of 12.2%-23.5%, consistent with recent observations. The mass-weighted ICL profiles of velocity dispersion, stellar age, and iron abundance generally agree with MaNGA measurements, and its clumpy east-west elongation is closely linked to the merger history of the two brightest cluster galaxies (BCGs). The simulation also successfully reproduces the north and west intracluster filaments (ICFs) detected by weak lensing, and the surrounding galaxy groups are distributed in alignment with the directions toward the nearby A2199 and A1367 clusters. The simulation further predict a complex cluster assembly history, with two major mergers at z=0.74 and 0.45, and a pericentric passage of the two BCGs occurring about 0.64 Gyr before its current state. Our results demonstrate that constrained simulations are a powerful tool for connecting observed structures to the unobservable assembly histories of individual galaxy clusters.