  We investigate the double parton distributions (DPDs) for a positronium-like
bound state using light-front QED. We incorporate the higher Fock three
particle component of the state, that includes a photon. We obtain the overlap
representation of the DPDs in terms of the three-particle light-front wave
functions (LFWFs). Our calculation explores the correlations between the
momentum fractions of the particles probed and the transverse distance between
them, without any assumption of factorization between them. We also investigate
the behavior of the DPDs near the kinematical boundary when the sum of the
momentum fractions is close to one.
