Reply to "Comment on 'Nature and entropy content of the ordering transitions in RCo2'"

In the previous Comment, Forker and co-workers claim that perturbed angular correlation (PAC) data leave no alternative to the conclusion that the spontaneous magnetization of PrCo2 and NdCo2 undergoes a discontinuous, first-order phase transition at TC. We show here that their claim is in clear contradiction with a wealth of experimental evidence, including our own. Finally, we propose a possible origin for the disagreement between their interpretation of the PAC results and the literature on this subject.

Evidence of the anomalous charge state 57Fe4+ in the nuclear decay of 57Co3+

The first observation of the elusive Fe4+ charge state coming from the nuclear decay of 57Co3+ has been found in the Mössbauer emission spectra of 57Co:La2Li0.5Co0.5O4. A Ti-doped sample was prepared in order to show that the Fe4+ fraction can be conveniently monitored. Both results were predicted on the basis of the electronic energy-band scheme of these oxides.

Towards and ab initio description of magnetism in ionic solids

The physical contributions to the KNiF3 magnetic exchange coupling integral have been obtained from specially designed ab initio cluster model calculations. Three important mechanisms have been identified. These are the delocalization of the magnetic orbitals into the anion p band, the variational contribution of the second-order interactions, and the many-body terms hidden in the two-body operator and the Heisenberg Hamiltonian.

First-principles periodic calculation of four-body spin terms in high-Tc cuprate superconductors

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 La2CuO4, the paradigm of high-Tc superconductor parent compounds, and for the SrCu2O3 ladder compound are reported. For La2CuO4, present results support experimental evidence by Toader et al. [Phys. Rev. Lett. 94, 197202 (2005)]. For SrCu2O3 even larger four-body spin amplitudes are found together with Jl/Jr=1 and non-negligible ferromagnetic interladder exchange.

First-principles periodic calculation of four-body spin terms in high-Tc Cuprate superconductors

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 La2CuO4, the paradigm of high-Tc superconductor parent compounds, and for the SrCu2O3 ladder compound are reported. For La2CuO4, present results support experimental evidence by Toader et al. [Phys. Rev. Lett. 94, 197202 (2005)]. For SrCu2O3 even larger four-body spin amplitudes are found together with Jl/Jr=1 and non-negligible ferromagnetic interladder exchange.

Towards and ab initio description of magnetism in ionic solids

The physical contributions to the KNiF3 magnetic exchange coupling integral have been obtained from specially designed ab initio cluster model calculations. Three important mechanisms have been identified. These are the delocalization of the magnetic orbitals into the anion p band, the variational contribution of the second-order interactions, and the many-body terms hidden in the two-body operator and the Heisenberg Hamiltonian.

Ground and excited states of KNiF3: An ab initio cluster-model approach

Finite cluster models and a variety of ab initio wave functions have been used to study the electronic structure of bulk KNiF3. Several electronic states, including the ground state and some charge-transfer excited states, have been considered. The study of the cluster-model wave functions has permitted an understanding of the nature of the chemical bond in the electronic ground state. This is found to be highly ionic and the different ionic and covalent contributions to the bonding have been identified and quantified. Finally, we have studied the charge-transfer excited states leading to the optical gap and have found that calculated and experimental values are in good agreement. The wave functions corresponding to these excited states have also been analyzed and show that although KNiF3 may be described as a ligand-to-metal charge-transfer insulator there is a strong configuration mixing with the metal-to-metal charge-transfer states.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Absence of collective effects in Heisenberg systems with localized magnetic moments

Existence of collective effects in magnetic coupling in ionic solids is studied by mapping spin eigenstates of the Heisenberg and exact nonrelativistic Hamiltonians on cluster models representing KNiF3, K2NiF4, NiO, and La2CuO4. Ab initio techniques are used to estimate the Heisenberg constant J. For clusters with two magnetic centers, the values obtained are about the same for models having more magnetic centers. The absence of collective effects in J strongly suggests that magnetic interactions in this kind of ionic solids are genuinely local and entangle only the two magnetic centers involved.

Magnetic coupling in the weac ferromagnet CuF2

CuF2 is known to be an antiferromagnetic compound with a weak ferromagnetism due to the anisotropy of its monoclinic unit cell (Dzialoshinsky-Moriya mechanism). We investigate the magnetic ordering of this compound by means of ab initio periodic unrestricted Hartree-Fock calculations and by cluster calculations which employ state-of-the-art configuration interaction expansions and modern density functional theory techniques. The combined use of periodic and cluster models permits us to firmly establish that the antiferromagnetic order arises from the coupling of one-dimensional subunits which themselves exhibit a very small ferromagnetic coupling between Cu neighbor cations. This magnetic order could be anticipated from the close correspondence between CuF2 and rutile crystal structures.

Ab Initio study of magnetic interactions in the KCuF3 and K2CuF4 low-dimensional Systems

The ab initio cluster model approach has been used to study the electronic structure and magnetic coupling of KCuF3 and K2CuF4 in their various ordered polytype crystal forms. Due to a cooperative Jahn-Teller distortion these systems exhibit strong anisotropies. In particular, the magnetic properties strongly differ from those of isomorphic compounds. Hence, KCuF3 is a quasi-one-dimensional (1D) nearest neighbor Heisenberg antiferromagnet whereas K2CuF4 is the only ferromagnet among the K2MF4 series of compounds (M=Mn, Fe, Co, Ni, and Cu) behaving all as quasi-2D nearest neighbor Heisenberg systems. Different ab initio techniques are used to explore the magnetic coupling in these systems. All methods, including unrestricted Hartree-Fock, are able to explain the magnetic ordering. However, quantitative agreement with experiment is reached only when using a state-of-the-art configuration interaction approach. Finally, an analysis of the dependence of the magnetic coupling constant with respect to distortion parameters is presented.

The role of dynamical polarization of the ligand-to-metal charge transfer excitations in the ab initio determination of effective exchange parameters

The role of the bridging ligand on the effective Heisenberg coupling parameters is analyzed in detail. This analysis strongly suggests that the ligand-to-metal charge transfer excitations are responsible for a large part of the final value of the magnetic coupling constant. This permits us to suggest a variant of the difference dedicated configuration interaction (DDCI) method, presently one of the most accurate and reliable for the evaluation of magnetic effective interactions. This method treats the bridging ligand orbitals mediating the interaction at the same level than the magnetic orbitals and preserves the high quality of the DDCI results while being much less computationally demanding. The numerical accuracy of the new approach is illustrated on various systems with one or two magnetic electrons per magnetic center. The fact that accurate results can be obtained using a rather reduced configuration interaction space opens the possibility to study more complex systems with many magnetic centers and/or many electrons per center.