Following recent interest in the qualitative analysis of some optimal control and shape optimisation problems, we provide in this article a detailed study of the optimisation of Robin boundary conditions in PDE constrained calculus of variations. Our main model consists of an elliptic PDE of the form ââu ÎČ = f (x, u ÎČ) endowed with the Robin boundary conditions âÎœ u ÎČ +ÎČ(x)u ÎČ = 0. The optimisation variable is the function ÎČ, which is assumed to take values between 0 and 1 and to have a fixed integral. Two types of criteria are under consideration: the first one is non-energetic criteria. In other words, we aim at optimising functionals of the form J (ÎČ) =\int_{⊠or ââŠ} j(u ÎČ). We prove that, depending on the monotonicity of the function j, the optimisers may be of bang-bang type (in other words, the optimisers write 1Î for some measurable subset Î of ââŠ) or, on the contrary, that they may only take values strictly between 0 and 1. This has consequence for a related shape optimisation problem, in which one tries to find where on the boundary Neumann (âÎœ u = 0) and constant Robin conditions (âÎœ u + u = 0) should be placed in order to optimise criteria. The proofs for this first case rely on new fine oscillatory techniques, used in combination with optimality conditions. We then investigate the case of compliance-type functionals. For such energetic functionals, we give an in-depth analysis and even some explicit characterisation of optimal ÎČ * .