In this talk, we introduce methods from convex optimization to solve the multi-marginal transport type problems arise in the context of density functional theory. Convex relaxations are used to provide outer approximation to the set of N-representable 2-marginals and 3-marginals, which in turn provide lower bounds to the energy. We further propose rounding schemes to obtain upper bound to the energy.