Heavy flavor muons at forward rapidities in $\sqrt{s_{NN}}=200$ GeV p+p, d+A, and A+A collisions

The dynamics of heavy flavor production and suppression in $A$+$A$ collisions play an important role in unraveling the properties of the quark-gluon plasma produced at RHIC. The related observables offer direct insight into strongly-coupled nature of the medium. It has also become increasingly apparent that the production and suppression in the absence of a hot medium need to be studied in detail, via comprehensive measurements in $p$+$p$, $p$($d$)+$A$, and $A$+$A$ collisions of varying size and energy at RHIC and the LHC. I will present the latest heavy flavor invariant yields and nuclear modification factors at forward rapidity in $\sqrt{s_{NN}}=200$ GeV $p$+$p$, $d$+Au, and $A$+$A$ collisions recorded at PHENIX and discuss their implications for our current understanding of heavy flavor production.