In this dissertation, we discuss the prospect of discovering flavor changing neutral Higgs (FCNH) interactions with quarks and leptons at current and future hadron colliders. Particularly, we have looked for t→ch0 and ϕ0→τ±μ∓, where ϕ0 could be a CP-even scalar \bigg[h0 (lighter), H0 (heavier) \bigg ] or a CP-odd pseudo-scalar (A0). A general two Higgs doublet model (gTHDM) is used to simulate t→ch0 and $\phi^0\to \tau^{\pm} \mu^{\mp} $ decays. The LHC measurements of the light Higgs boson (h0) favor the decoupling limit of gTHDM, in which the couplings of h0 approach Standard Model values. In this limit, FCNH couplings of the light Higgs boson h0 are naturally suppressed by a small mixing parameter cos⁡(β−α), while the FCNH couplings of heavier neutral Higgs bosons H0,andA0 are sustained by sin⁡(β−α)∼1. Promising results are found for the LHC collision energies s=13 TeV and 14 TeV. In addition, we study the discovery potential of future pp colliders, with $\sqrt{s} = $ 27 TeV and 100 TeV. For ϕ0→τ±μ∓,
we evaluate the production rate of physics background from dominant processes (τ+τ−,WW,ZZ,Wq,Wg,tt¯) with realistic acceptance cuts and tagging efficiencies. For t→ch0, where top is coming from top pair production, we consider h0→WW∗→ℓ+ℓ−+\slashedET and h0→τ+τ−→ℓ+ℓ−+\slashedET and another top decaying hadronically to a b quark and two light jets. For this report we have studied h0→WW∗ and h0→τ+τ− separately. Our analysis suggests a reach of 5σ or better, with integrated luminosity L = 3 ab−1 and s = 13, 14 and 27 TeV for λtc≤0.064, under the current ATLAS limits for both light Higgs decay mode separately. For h0→WW∗ we have also presented the discovery potential at 100 TeV.