The result of a search for the pair production of the lightest supersymmetric partner of the bottom quark b ˜ 1 using 139 fb −1 of proton-proton data collected at √s = 13 TeV by the ATLAS detector is reported. In the supersymmetric scenarios considered both of the bottom-squarks decay into a b-quark and the second-lightest neutralino, b ˜ 1 → b + χ ˜ 2 0. Each χ ˜ 2 0 is assumed to subsequently decay with 100% branching ratio into a Higgs boson (h) like the one in the Standard Model and the lightest neutralino: χ ˜ 2 0 → h + χ ˜ 1 0. The χ ˜ 1 0 is assumed to be the lightest supersymmetric particle (LSP) and is stable. Two signal mass configurations are targeted: the first has a constant LSP mass of 60 GeV; and the second has a constant mass difference between the χ ˜ 2 0 and χ ˜ 1 0 of 130 GeV. The final states considered contain no charged leptons, three or more b-jets, and large missing transverse momentum. No significant excess of events over the Standard Model background expectation is observed in any of the signal regions considered. Limits at the 95% confidence level are placed in the supersymmetric models considered, and bottom-squarks with mass up to 1.5 TeV are excluded.