The direct detection of gravitational waves opens new perspectives for measuring properties of gravitationally bound compact objects. It is then important to investigate black holes and neutron stars in alternative theories of gravity, since they can have features that make them observationally distinguishable from their General Relativity (GR) counterparts. In this work, we examine a special case of vector Galileons, a vector-tensor theory of gravity with interesting cosmological properties, which consists of a one parameter modification of the Einstein-Maxwell action. Within this theory, we study configurations describing asymptotically flat, spherically symmetric black holes and neutron stars. The set of black hole solutions in this theory is surprisingly rich, generalising results found in GR or in related scalar-tensor theories. We investigate the properties and conserved charges of black holes, using both analytical and numerical techniques, highlighting configurations that are more compact than in GR. We then study properties of neutron stars, showing how the vector profile can influence the star internal structure. Depending on properties of matter and fields inside the star, neutron stars can be more massive than in GR, and they can be even more compact than Schwarzschild black holes, making these objects observationally interesting. We also comment on possible extensions of our configurations to magnetically charged or rotating configurations.

In the context of a cubic Galileon model in which the Vainshtein mechanism suppresses the scalar field interactions with matter, we study low-density stars with slow rotation and static relativistic stars. We develop an expansion scheme to find approximated solutions inside the Vainshtein radius, and show that deviations from General Relativity (GR), while considering rotation, are also suppressed by the Vainshtein mechanism. In a quadratic coupling model, in which the scalarisation effect can significantly enhance deviations from GR in normal scalar tensor gravity, the Galileon term successfully suppresses the large deviations away from GR. Moreover, using a realistic equation of state, we construct solutions for a relativistic star, and show that deviations from GR are more suppressed for higher density objects. However, we found that the scalar field solution ceases to exist above a critical density, which roughly corresponds to the maximum mass of a neutron star. This indicates that, for a compact object described by a polytropic equation of state, the configuration that would collapse into a black hole cannot support a non-trivial scalar field.