Abstract This paper presents the muon momentum calibration and performance studies for the ATLAS detector based on the pp collisions data sample produced at $$\sqrt{s}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msqrt> <mml:mi>s</mml:mi> </mml:msqrt> </mml:math> = 13 TeV at the LHC during Run 2 and corresponding to an integrated luminosity of 139 $${\textrm{fb}}^{-1}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msup> <mml:mrow> <mml:mtext>fb</mml:mtext> </mml:mrow> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> </mml:math> . An innovative approach is used to correct for potential charge-dependent momentum biases related to the knowledge of the detector geometry, using the $$Z\rightarrow \mu ^{+}\mu ^{-}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>Z</mml:mi> <mml:mo>→</mml:mo> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>-</mml:mo> </mml:msup> </mml:mrow> </mml:math> resonance. The muon momentum scale and resolution are measured using samples of $$J/\psi \rightarrow \mu ^{+}\mu ^{-}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>J</mml:mi> <mml:mo>/</mml:mo> <mml:mi>ψ</mml:mi> <mml:mo>→</mml:mo> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>-</mml:mo> </mml:msup> </mml:mrow> </mml:math> and $$Z\rightarrow \mu ^{+}\mu ^{-}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>Z</mml:mi> <mml:mo>→</mml:mo> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>-</mml:mo> </mml:msup> </mml:mrow> </mml:math> events. A calibration procedure is defined and applied to simulated data to match the performance measured in real data. The calibration is validated using an independent sample of $$\Upsilon \rightarrow \mu ^{+}\mu ^{-}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>Υ</mml:mi> <mml:mo>→</mml:mo> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:msup> <mml:mi>μ</mml:mi> <mml:mo>-</mml:mo> </mml:msup> </mml:mrow> </mml:math> events. At the Z $$(J/\psi )$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mo>(</mml:mo> <mml:mi>J</mml:mi> <mml:mo>/</mml:mo> <mml:mi>ψ</mml:mi> <mml:mo>)</mml:mo> </mml:mrow> </mml:math> peak, the momentum scale is measured with an uncertainty at the 0.05% (0.1%) level, and the resolution is measured with an uncertainty at the 1.5% (2%) level. The charge-dependent bias is removed with a dedicated in situ correction for momenta up to 450 GeV with a precision better than 0.03 $${\textrm{TeV}}^{-1}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msup> <mml:mrow> <mml:mtext>TeV</mml:mtext> </mml:mrow> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> </mml:math> .