First-principles calculations, in combination with the four-state energy
mapping method, are performed to extract the magnetic interaction parameters of
multiferroic BiFeO3. Such parameters include the symmetric exchange (SE)
couplings and the Dzyaloshinskii-Moriya (DM) interactions up to second nearest
neighbors, as well as the single ion anisotropy (SIA). All magnetic parameters
are obtained not only for the R3c structural ground state, but also for the
R3m and R3ˉc phases in order to determine the effects of
ferroelectricity and antiferrodistortion distortions, respectively, on these
magnetic parameters. In particular, two different second-nearest neighbor
couplings are identified and their origins are discussed in details. Moreover,
Monte-Carlo (MC) simulations using a magnetic Hamiltonian incorporating these
first-principles-derived interaction parameters are further performed. They
result (i) not only in the accurate prediction of the spin-canted G-type
antiferromagnetic structure and of the known magnetic cycloid propagating along
a direction, as well as their unusual characteristics (such
as a weak magnetization and spin-density-waves, respectively); (ii) but also in
the finding of another cycloidal state of low-energy and that awaits to be
experimentally confirmed. Turning on and off the different magnetic interaction
parameters in the MC simulations also reveal the precise role of each of them
on magnetism