HMI Data Driven MHD Model Predicted Active Region Photospheric Heating Rates: Their Scale Invariant Power Law Distributions, and Their Plausible Temporal Correlation With Flares