A general mapping between the energy of pertinent magnetic solutions and the diagonal terms of the spin Hamiltonian in a local representation provides the first general framework to extract accurate values for the many body terms of extended spin Hamiltonians from periodic first-principle calculations. Estimates of these terms for La 2 CuO 4 , the paradigm of high-T c superconductor parent compounds, and for the SrCu 2 O 3 ladder compound are reported. For La 2 CuO 4 , present results support experimental evidence by Toader et al. [Phys. Rev. Lett. 94, 197202 (2005)]. For SrCu 2 O 3 even larger four-body spin amplitudes are found together with J l =J r 1 and non-negligible ferromagnetic interladder exchange. DOI: 10.1103/PhysRevLett.97.087003 PACS numbers: 74.20.Mn, 71.15.ÿm, 74.25.Jb, 75.30.Et The strong antiferromagnetic interactions observed in lamellar cuprates are fundamental ingredients of the high-T c (HTC) superconductivity microscopic mechanism [1,2]. Magnetic interactions arise from the particular crystal and electronic structure of these cuprates with Cu 2 ions arranged in edge sharing Cu 4 O 4 plaquettes. The electronic ground state involves a single d x2ÿy2 -type hole in each Cu3d shell leading to a network of effective spin S 1=2 particles. Nevertheless, these systems are strongly correlated in nature, making standard band theory techniques unable to accurately describe either their valence or low energy spectrum [3,4]. The low energy spectrum and collective properties of these compounds are assumed to be governed by a Heisenberg Hamiltonian as in the first term of Eq. (1) accounting for the magnetic coupling J ij between nearestneighbor (NN) centers i and j only. This is in agreement with the widely accepted general picture for HTC superconductivity involving a ''Heisenberg sea'' where holes are introduced by doping the perfect structures. However, to fully understand the magnetic excitations and the infrared and neutron scattering spectra of 2D [5][6][7][8][9][10][11] and spin ladder cuprates [12,13], it has been necessary to extend the spin Hamiltonian as in Eq. (1), [18] provided definitive evidence of its existence in La 2 CuO 4 . They suggest J ring 0:5J, comparable to the pairing energies, and propose that the resulting circulating currents could have an important role in the mechanism of superconductivity. Notwithstanding, previous estimates of J ring for 2D and spin ladder cuprates, obtained from either indirect measurements or numerical simulations with an extended Heisenberg model, propose substantially smaller amplitudes with J ring 0:30J [7,9,10,12,19] Clearly, a bottom-up accurate and independent determination of all important terms in the spin Hamiltonian in Eq. (1) for La 2 CuO 4 -or any other similar system-is PRL 97,