We study multi-Higgs final states in vector boson fusion (VBF) processes at the LHC and at future proton-proton colliders, focusing on the prospects for measurements at 27 TeV and at 100 TeV. We use an effective Lagrangian which includes higher-dimensional operators in the mass eigenstates which are relevant to VBF processes and relate this to specific param-At the LHC, the dominant SM process of multi-Higgs production is gluon-gluon-fusion (ggF) via a heavy top-quark loop. Current LHC data constrain the triple Higgs self-coupling only very weakly [3], and not at all the quartic Higgs self-coupling. Multiple groups have evaluated the potential for a first meaningful measurement of the triple Higgs self-coupling at future highluminosity runs of the LHC [4-10]. The considered decay channels of the Higgs pair include and bbbb [18, 22, 23]. It is expected that the triple self-coupling can be constrained within 40% accuracy after collecting 3 ab −1 of data at the 14 TeV LHC [24]. There will be no access to the quartic self-coupling. Beyond the LHC, at a future 100 TeV hadron collider, the Higgs pair-production rate is enhanced significantly [9,13,[25][26][27][28][29][30], allowing for a more accurate determination of the Higgs potential.A subdominant process of multi-Higgs production in hadron collisions is vector-bosonfusion (VBF) [31], W + W − → hh, where the W bosons are effectively radiated from incoming quarks. In addition to its dependence on the Higgs self-interaction, this process also gives access to the hhW W coupling, which is another SM interaction that has not been accessible by other means. In fact, the hVV (V = W ± , Z ) couplings can be derived from the precise measurement of the decay branching fractions h → W W * and h → Z Z * at the LHC, up to a common normalization factor. The current LHC data on the decay branching fractions for these channels are consistent with the SM predictions [32, 33]. By contrast, there are no significant constraints from data on hhV V direct couplings.The VBF mode of double-Higgs production at hadron colliders has been studied in Refs. [34][35][36][37][38]. Beyond tree level, the NLO QCD correction enhances the cross section by ∼ 7% [39,40]. In the high-luminosity mode of the LHC (HL-LHC) with 3000 fb −1 at 14 TeV, the hhV V interaction can be constrained to 20%. A 100 TeV hadron collider has the potential to reduce the uncertainty down to 1% [37]. The hhW W coupling is also accessible in the W ± W ± h final state. In Ref. [41] it is found that this particular final state can constrain this coupling to O(100%) at the HL-LHC, and to 20% at a 100 TeV collider.The other modes of double Higgs production are tt hh or V hh production. The production rates of these mode are smaller than ggF and VBF mode. The discussion of these modes is beyond the scope of current work. For interested readers, we refer to [39,40,[42][43][44][45][46][47][48].Measuring the quartic Higgs self-coupling in the SM is a much more challenging task. At the 14 TeV LHC, the cross section of g g → hhh is onl...