Abstract
Phase transition in d=3 ferromagnetic Ising models with random fields is analyzed directly at the zero-temperature critical point. Critical behavior is extracted from correlation functions averaged over an ensemble of exact ground states obtained with a new integer optimization algorithm. For Gaussian distribution of random fields finite-size scaling demonstrates a continuous phase transition with effective disconnected susceptibility exponent η̃0.9, correlation-length exponent ν1.0, and magnetization exponent β0.05.