that would change the norm from L2 towards L-infinity, which shouldn’t change the shape of the response peaks (so not what you want), but it should change the math fundamentally.
if you just wanted to apply an extra squaring on the response array, that wouldn’t affect its monotonicity.
what you need is to find local maxima. non-maximum suppression can be constructed from a grayscale dilation and equality comparison.