We report results on the proton mass decomposition and also on related quark and glue momentum fractions. The results are based on overlap valence fermions on four ensembles of N f = 2 + 1 DWF configurations with three lattice spacings and three volumes, and several pion masses including the physical pion mass. With fully non-perturbative renormalization (and universal normalization on both quark and gluon), we find that the quark energy and glue field energy contribute 33(4)(4)% and 37(5)(4)% respectively in the M S scheme at µ = 2 GeV. A quarter of the trace anomaly gives a 23 (1)(1)% contribution to the proton mass based on the sum rule, given 9(2)(1)% contribution from the u, d, and s quark scalar condensates. The u, d, s and glue momentum fractions in the M S scheme are in good agreement with global analyses at µ = 2 GeV.Introduction: In the standard model, Higgs boson provides the origin of quark masses. But how it is related to the proton mass and thus the masses of nuclei and atoms is another question. The masses of the valence quarks in the proton are just ∼3 MeV per quark which is directly related to the Higgs boson, while the total proton mass is 938 MeV. The percentage of the quark and gluon contributions to the proton mass can only be provided by solving QCD non-perturbatively, and/or with information from experiment. With phenomenological input, the first decomposition was carried out by Ji [1]. As in Refs. [1,2], the Hamiltonian of QCD can be decomposed asin the rest frame of the hadron state where M is the hadron mass, T µν is the energy momentum tensor of QCD with T 44 as its expectation value in the hadron, and the trace anomaly gives