2026ApJ..1000..100D

The nascent methodology of applying the principles of causal discovery to astrophysical data has produced affirming results about deeply held theories concerning the causal nature behind the observed coevolution of supermassive black holes (SMBHs) with their host galaxies. The key results from observations have demonstrated an apparent causal reversal across different galaxy morphologies-SMBHs causally influence the evolution of the physical parameters of their spiral galaxy hosts, whereas SMBHs in elliptical galaxies are passive companions that grow in near lockstep with their hosts. To further explore and ascertain insights, it is necessary to utilize galaxy simulations to track the time evolution of the observed causal relations to learn more about the temporal nature of the changing SMBH-galaxy evolutionary directions. We conducted experiments with the NIHAO suite of cosmological zoom-in hydrodynamical simulations to follow the evolution of individual galaxies along with their central SMBH masses (M_•) and properties including central stellar velocity dispersion (σ₀). We reproduce the causal results from real galaxies, but add clarity by observing the SMBH-galaxy causal directions are noticeably inverted between the epochs before and after the peak of star formation. The implications for causal reversal of the M_•-σ₀ relation portend larger concerns about the reliability of SMBH masses estimated at high redshifts and presumptions of overmassive black holes at early epochs. Toward this problem, we apply updated causally informed scaling relations that predict high-z black hole masses that are approximately 2 orders of magnitude less massive, and thus not overmassive with respect to local z = 0 SMBH-galaxy mass ratios.