Higgs boson decays in flexible brane world models with stable, massive gravi-vectors are considered. Such vectors couple bilinearly to the standard model fields through either the standard model energymomentum tensor, the weak hypercharge field strength, or the Higgs scalar. The role of the coupling involving the extrinsic curvature is highlighted. It is found that within the presently allowed parameter space, the decay rate of the Higgs into two gravi-vectors (which would appear as an invisible Higgs decay) can be comparable to the rate for any of the standard model decay modes.Various theoretical extensions of theories of gravity include vector particles. In particular, such gravi-vectors [1] appear in flexible brane world models in which a fourdimensional space-time is embedded in a higherdimensional space-time thus breaking the extradimensional spatial translation symmetries [2]. When these symmetries are made local thereby including higherdimensional gravitational interactions, the erstwhile Nambu-Goldstone scalar degrees of freedom [3] associated with the higher-dimensional spatial translation symmetry breakdown become the longitudinal components of the now massive vector particles, X i . For N ! 2 additional compactified isotropic spatial dimensions, these vectors, which are completely neutral under the standard model gauge group, carry an additional SOðNÞ quantum number, labeled by i ¼ 1; . . . ; N, which reflects the isometry of the codimensional space when the four-dimensional brane is embedded in the larger-dimensional space-time. On the other hand, all standard model particles are SOðNÞ singlets. Consequently, SOðNÞ invariant interactions of these vector Proca fields require them to appear in pairs and they are thus massive, stable physical degrees of freedom. For a single codimension, N ¼ 1, the vector is stable provided there is a unbroken parity with respect to the extra dimension under which the brane vector is odd. Being stable, the vectors are also candidates for the dark matter of the Universe and are thus subject to the appropriate constraints [4].The method of nonlinear realizations is a powerful tool which can be used to extract the generic form of vector interactions in a model independent given only the structure of the particular symmetry breakdown. In previous work [1], we applied this coset construction to the case of local space-time symmetries in brane world models. The resulting interaction terms of the resultant vector field will then contain various model-dependent parameters such as masses and dimensionless couplings which are not fixed by the method of nonlinear realizations but must be determined either using explicit models or by experiment. In the present case, the strength of the vector interactions is governed by the ratiowhere M X is the vector mass whose nonzero value is model dependent and F 4 X is the brane tension. In the flexible brane limit where the brane tension is much smaller than the D-dimensional Planck scale, the Kaluza-Klein modes of higher-dimensional gravity...