This Letter presents a search for heavy charged long-lived particles produced in proton-proton collisions at $\sqrt{s} = 13$ TeV at the LHC using a data sample corresponding to an integrated luminosity of 36.1 fb$^{-1}$ collected by the ATLAS experiment in 2015 and 2016. These particles are expected to travel with a velocity significantly below the speed of light, and therefore have a specific ionisation higher than any high-momentum Standard Model particle of unit charge. The pixel subsystem of the ATLAS detector is used in this search to measure the ionisation energy loss of all reconstructed charged particles which traverse the pixel detector. Results are interpreted assuming the pair production of $R$-hadrons as composite colourless states of a long-lived gluino and Standard Model partons. No significant deviation from Standard Model background expectations is observed, and lifetime-dependent upper limits on $R$-hadron production cross-sections and gluino masses are set, assuming the gluino always decays in two quarks and a stable neutralino. $R$-hadrons with lifetimes above 1.0 ns are excluded at the 95% confidence level, with lower limits on the gluino mass ranging between 1290 GeV and 2060 GeV. In the case of stable $R$-hadrons, the lower limit on the gluino mass at the 95% confidence level is 1890 GeV.
A search for new phenomena in final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair, jets, and large missing transverse momentum is presented. This analysis makes use of proton--proton collision data with an integrated luminosity of $36.1 \; \mathrm{fb}^{-1}$, collected during 2015 and 2016 at a centre-of-mass energy $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. The search targets the pair production of supersymmetric coloured particles (squarks or gluinos) and their decays into final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair and the lightest neutralino ($\tilde{\chi}_1^0$) via one of two next-to-lightest neutralino ($\tilde{\chi}_2^0$) decay mechanisms: $\tilde{\chi}_2^0 \rightarrow Z \tilde{\chi}_1^0$, where the $Z$ boson decays leptonically leading to a peak in the dilepton invariant mass distribution around the $Z$ boson mass; and $\tilde{\chi}_2^0 \rightarrow \ell^+\ell^- \tilde{\chi}_1^0$ with no intermediate $\ell^+\ell^-$ resonance, yielding a kinematic endpoint in the dilepton invariant mass spectrum. The data are found to be consistent with the Standard Model expectation. Results are interpreted using simplified models, and exclude gluinos and squarks with masses as large as 1.85 TeV and 1.3 TeV at 95% confidence level, respectively.