A search for neutral Higgs bosons of the minimal supersymmetric standard model (MSSM) and for a heavneutral Zâ² boson is performed using a data sample corresponding to an integrated luminosity of 3.2Â fb- 1 from protonâproton collisions at s=13 TeV recorded by the ATLAS detector at the LHC. The heavy resonance is assumed to decay to a Ï+Ï- pair with at least one Ï lepton decaying to final states with hadrons and a neutrino. The search is performed in the mass range of 0.2â1.2 TeV for the MSSM neutral Higgs bosons and 0.5â2.5 TeV for the heavy neutral Zâ² boson. The data are in good agreement with the background predicted by the Standard Model. The results are interpreted in MSSM and Zâ² benchmark scenarios. The most stringent constraints on the MSSM mAâtan Î² space exclude at 95Â % confidence level (CL) tan Î²> 7.6 for mA= 200 GeV in the mhmod+ MSSM scenario. For the Sequential Standard Model, a ZSSMâ² mass up to 1.90 TeV is excluded at 95Â % CL and masses up to 1.82â2.17 TeV are excluded for a ZSFMâ² of the strong flavour model.