A search optimized for new heavy particles decaying to two b-quarks and produced in association with additional b-quarks is reported. The sensitivity is improved by b-tagging at least one lower-pT jet in addition to the two highest-pT jets. The data used in this search correspond to an integrated luminosity of 103 fb−1 collected with a dedicated trijet trigger during the 2017 and 2018 s=13 TeV proton-proton collision runs with the ATLAS detector at the LHC. The search looks for resonant peaks in the b-tagged dijet invariant mass spectrum over a smoothly falling background. The background is estimated with an innovative data-driven method based on orthonormal functions. The observed b-tagged dijet invariant mass spectrum is compatible with the background-only hypothesis. Upper limits at 95% confidence level on a heavy vector-boson production cross section times branching ratio to a pair of b-quarks are derived.Received 23 August 2021Accepted 11 November 2021DOI:https://doi.org/10.1103/PhysRevD.105.012001Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP3.© 2022 CERN, for the ATLAS CollaborationPhysics Subject Headings (PhySH)Research AreasExtensions of Higgs sectorExtensions of gauge sectorParticles & Fields