Historically, the even-order PN equations have been considered a less accurate approximation to the transport equation than the odd-order PN-1 equations. This perception has stemmed from two apparent conceptual difficulties imposed by the even-order PN methods— the difficulty in prescribing rigorous boundary conditions for even-order PN equations that contain the odd number of angular flux moments and the discontinuous character of the even-order PN solutions at material interfaces. With the first one of the mentioned even-order PN conceptual problems, a presentation is made of a straightforward and physically-motivated variational procedure based on a new functional that leads from a multigroup planar geometry transport problem to a multigroup P2 problem with clearly and rigorously defined multigroup boundary conditions. These boundary conditions are new and allow neutron transfer between energy groups at the boundary. These boundary conditions are tested by comparing P2, P1, and SN calculations. Our results show that in the test problems considered, the multigroup P2 equations with variational boundary conditions are always more accurate than the P1 multigroup equations with Federighi-Pomraning or Marshak boundary conditions applied to each energy group.