Shell model Monte Carlo techniques are used to calculate two-neutrino double beta decay matrix elements. We test the approach against direct diagonalization for 48 Ca in the complete pf shell using the KB3 interaction. The method is then applied to the decay of 76 Ge in the ͑0f 5͞2 , 1p, 0g 9͞2 ͒ model space using a newly calculated realistic interaction. Our result for the matrix element is 0.12 6 0.05 MeV 21 , in reasonable agreement with the experimental value.PACS numbers: 21.60. Cs, 21.60.Ka, 27.40.+z, The double beta (bb) decay of a nucleus is a rare second order weak process [1,2]. The as-yet unobserved neutrinoless mode is of fundamental interest, as it would signal a neutrino mass, lepton number nonconservation, or admixtures of right handed weak currents. In contrast, the existence of the 2n mode has been firmly established (see the review in Ref. [2]). The ability to accurately describe this latter process is an important element in the interpretation of limits on neutrinoless decays. Unfortunately, it seems that 2n matrix elements are highly suppressed and so depend sensitively on small, poorly determined parts of the nuclear wave functions.Most recent calculations of 2nbb matrix elements for nuclei heavier than 48 Ca rely on the quasiparticle random phase approximation (QRPA) [2]. While this approach is computationally simple and includes many features of the two-body interaction known to be relevant for bb decay, the calculated matrix elements are uncertain because of their great sensitivity to the J 1 1 , T 0 particleparticle interaction [3]. The interacting shell model offers a more microscopic approach to the problem. Complete 0hv shell model calculations [4] not only recover more quenching of Gamow-Teller (GT) transitions than QRPA calculations, but also are in agreement with observations (after the universal renormalization of g A to 1.0). However, computational limitations have restricted shell model calculations of the 2nbb decay matrix element to 48 Ca, the lightest of all bb candidates.In this Letter we show how shell model Monte Carlo (SMMC) methods can be used to calculate 2nbb decay matrix elements. We first calculate the decay of 48 Ca in the complete pf shell and validate our method against direct diagonalization. We then present results for the decay of 76 Ge, one of the few nuclei where the 2nbb decay has been precisely measured and where the best limits on the 0n decay have been obtained [5]. Our calculation in the complete ͑0f 5͞2 , 1p, 0g 9͞2 ͒ orbitals, which is impractical using traditional shell model methods, is the first of 2nbb decay in such a large model space.The 2nbb matrix elements between the 0 1 ground states of the initial and final even-even nuclei is given by [6]Here, j0 ͘ and E an average energy denominator. As there is no prescription for choosing the average energy denominator (and even the closure matrix element is usually calculated in a severely truncated basis), the uncertainty in this approximation is difficult to estimate.To calculate the exact 2nbb...