Abstract A search for pair-produced scalar or vector leptoquarks decaying into a b -quark and a $$\tau $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>τ</mml:mi> </mml:math> -lepton is presented using the full LHC Run 2 (2015–2018) data sample of 139 fb $$^{-1}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msup> <mml:mrow /> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> </mml:math> collected with the ATLAS detector in proton–proton collisions at a centre-of-mass energy of $$\sqrt{s} =13$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msqrt> <mml:mi>s</mml:mi> </mml:msqrt> <mml:mo>=</mml:mo> <mml:mn>13</mml:mn> </mml:mrow> </mml:math> TeV. Events in which at least one $$\tau $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>τ</mml:mi> </mml:math> -lepton decays hadronically are considered, and multivariate discriminants are used to extract the signals. No significant deviations from the Standard Model expectation are observed and 95% confidence-level upper limits on the production cross-section are derived as a function of leptoquark mass and branching ratio $$\mathcal {B}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>B</mml:mi> </mml:math> into a $$\tau $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>τ</mml:mi> </mml:math> -lepton and b -quark. For scalar leptoquarks, masses below 1460 GeV are excluded assuming $$\mathcal {B}=100$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>B</mml:mi> <mml:mo>=</mml:mo> <mml:mn>100</mml:mn> </mml:mrow> </mml:math> %, while for vector leptoquarks the corresponding limit is 1650 GeV (1910 GeV) in the minimal-coupling (Yang–Mills) scenario.