The non-linear, non-isothermal, magnetohydrodynamic (MHD) laminar convection flows of Buongiorno's nanofluid past a vertical surface with Darcy-Forchheimer model is mathematically investigated in the present article. Keller's Box implicit finite difference technique is utilized to solve the dimensionless conservation equations. Graphical and tabulated results are analyzed to study the behavior of primary and secondary velocity, temperature, nanoparticle volume fraction, shear stress rate, heat transfer rate and mass transfer rate for various emerging thermos-physical parameters. The Hall current and ion-slip current effects are also considered. Validations of earlier solutions in the literature is also included. The study finds applications in nanomaterial fabrication processing, biomedical, polymer processing, chemical engineering, crude oil purifying, etc. 1 Nomenclature: σ, electric conductivity of the fluid; η, non-dimensional radial coordinate; μ, dynamic viscosity; ξ, non-dimensional tangential coordinate; ψ, non-dimensional stream function; ν, Kinematic viscosity; ϕ, dimensionless concentration; θ, dimensionless temperature; τ, ratio of effective heat capacity of nanoparticle to the heat capacity of the fluid; (ρc) m , effective heat capacity; β e , Hall current parameter; ω e , electron frequency; τ e , electorn collision time; ρ f , density of the fluid; β i , ion-slip parameter; ρ p , density of the fluid; ∞, free stream condition; a, constant; B, magnetic field; b, Forchheimer inertial drag parameter; B 0 , imposed magnetic field;