Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2111.07621 · Physical review. B./Physical review. B

Critical spin dynamics of Heisenberg ferromagnets revisited

Dmytro Tarasevych, Peter Kopietz

Published Jan 4, 2022• 9 citations
arXivHTMLPDF

Abstract

We calculate the dynamic structure factor $S(\mathbit{k},\ensuremath{\omega})$ in the paramagnetic regime of quantum Heisenberg ferromagnets for temperatures $T$ close to the critical temperature ${T}_{c}$ using our recently developed functional renormalization group approach to quantum spin systems. In $d=3$ dimensions, we find that for small momenta $\mathbit{k}$ and frequencies $\ensuremath{\omega}$ the dynamic structure factor assumes the scaling form $S(\mathbit{k},\ensuremath{\omega})=(\ensuremath{\tau}TG(\mathbit{k})/\ensuremath{\pi})\mathrm{\ensuremath{\Phi}}(k\ensuremath{\xi},\ensuremath{\omega}\ensuremath{\tau})$, where $G(\mathbit{k})$ is the static spin-spin correlation function, $\ensuremath{\xi}$ is the correlation length, and the characteristic timescale $\ensuremath{\tau}$ is proportional to ${\ensuremath{\xi}}^{5/2}$. We explicitly calculate the dynamic scaling function $\mathrm{\ensuremath{\Phi}}(x,y)$ and find satisfactory agreement with neutron scattering experiments probing the critical spin dynamics in EuO and EuS. Precisely at the critical point where $\ensuremath{\xi}=\ensuremath{\infty}$, our result for the dynamic structure factor can be written as $S(\mathbit{k},\ensuremath{\omega})={(\ensuremath{\pi}{\ensuremath{\omega}}_{k})}^{\ensuremath{-}1}{T}_{c}G(\mathbit{k}){\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\omega}/{\ensuremath{\omega}}_{k})$, where ${\ensuremath{\omega}}_{k}\ensuremath{\propto}{k}^{5/2}$. We find that ${\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\nu})$ vanishes as ${\ensuremath{\nu}}^{\ensuremath{-}13/5}$ for large $\ensuremath{\nu}$, and as ${\ensuremath{\nu}}^{3/5}$ for small $\ensuremath{\nu}$. While the large-frequency behavior of ${\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\nu})$ is consistent with calculations based on mode-coupling theory and with perturbative renormalization group calculations to second order in $\ensuremath{\epsilon}=6\ensuremath{-}d$, our result for small frequencies disagrees with previous calculations. We argue that until now, neither experiments nor numerical simulations have been sufficiently accurate to determine the low-frequency behavior of ${\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\nu})$. We also calculate the low-temperature behavior of $S(\mathbit{k},\ensuremath{\omega})$ in one- and two-dimensional ferromagnets and find that it satisfies dynamic scaling with exponent $z=2$ and exhibits a pseudogap for small frequencies.

DOI
10.1103/physrevb.105.024403

Related papers

Dissipative spin dynamics in hot quantum paramagnets

Dmytro Tarasevych, Peter Kopietz

We use the functional renormalization approach for quantum spin systems developed by Krieg and Kopietz [Phys. Rev. B 99, 060403(R) (2019)] to calculate the spin-spin correlation function $G(\mathbit{k},\ensuremath{\omega})$ of quantum Heisenberg magnets at infinite temperature. For small wave vectors $\mathbit{k}$ and frequencies $\ensuremath{\omega}$ we find that $G(\mathbit{k},\ensuremath{\omega})$ assumes in dimensions $d>2$ the diffusive form predicted by hydrodynamics. In three dimensions our result for the spin-diffusion coefficient $\mathcal{D}$ is somewhat smaller than previous theoretical predictions based on the extrapolation of the short-time expansion, but is still about $30%$ larger than the measured high-temperature value of $\mathcal{D}$ in the Heisenberg ferromagnet ${\mathrm{Rb}}_{2}{\mathrm{CuBr}}_{4}\ifmmode\cdot\else\textperiodcentered\fi{}2{\mathrm{H}}_{2}\mathrm{O}$. In reduced dimensions $d\ensuremath{\le}2$ we find superdiffusion characterized by a frequency-dependent complex spin-diffusion coefficient $\mathcal{D}(\ensuremath{\omega})$ which diverges logarithmically in $d=2$, and as a power-law $\mathcal{D}(\ensuremath{\omega})\ensuremath{\propto}{\ensuremath{\omega}}^{\ensuremath{-}1/3}$ in $d=1$. Our result in one dimension implies scaling with dynamical exponent $z=3/2$, in agreement with recent calculations for integrable spin chains. Our approach is not restricted to the hydrodynamic regime and allows us to calculate the dynamic structure factor $S(\mathbit{k},\ensuremath{\omega})$ for all wave vectors. We show how the short-wavelength behavior of $S(\mathbit{k},\ensuremath{\omega})$ at high temperatures reflects the relative sign and strength of competing exchange interactions.

arXiv:2104.04270•Physical review. B./Physical review. B•Jul 16, 2021•17 citations

Spin-dynamics study of the dynamic critical behavior of the three-dimensional classical Heisenberg ferromagnet

Kun Chen, David Paul Landau

Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L body-centered-cubic classical Heisenberg ferromagnet with L\ensuremath{\le}40 in the vicinity of the critical point ${\mathit{T}}_{\mathit{c}}$. The temporal evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method with initial spin configurations generated by Monte Carlo methods. The space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms were calculated to determine the neutron-scattering functions. We developed a dynamic finite-size scaling theory for the neutron-scattering function at ${\mathit{T}}_{\mathit{c}}$ and used it to extract the dynamic critical exponent z. Within our resolution limit, the value of z was estimated to be 2.478(28), in excellent agreement with the dynamic scaling prediction, and dynamic scaling was found to be valid for momentum transfer q up to 0.4\ensuremath{\pi} in the (1,0,0) direction.

Physical review. B, Condensed matter•Feb 1, 1994•88 citations

Spin-dynamics simulations of the magnetic dynamics of RbMnF3 and direct comparison with experiment

Shan-Ho Tsai, Alex Bunker, David Paul Landau

Spin-dynamics techniques have been used to perform large-scale simulations of the dynamic behavior of the classical Heisenberg antiferromagnet in simple cubic lattices with linear sizes $L<~60.$ This system is widely recognized as an appropriate model for the magnetic properties of ${\mathrm{RbMnF}}_{3}.$ Time evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using an algorithm implemented by Krech et al., which is based on fourth-order Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical temperature. The dynamic critical exponent was estimated to be $z=(1.43\ifmmode\pm\else\textpm\fi{}0.03),$ which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Direct, quantitative comparisons of both the dispersion curve and the line shapes obtained from our simulations with very recent experimental results for ${\mathrm{RbMnF}}_{3}$ are presented.

arXiv:cond-mat/9910429•Physical review. B, Condensed matter•Jan 1, 2000•28 citations

Critical dynamics in the two-dimensional classicalXYmodel: A spin-dynamics study

Hans Gerd Evertz, David Paul Landau

Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the classical three component XY model (i.e., the anisotropic limit of an easy-plane Heisenberg ferromagnet), on square lattices of size up to ${192}^{2}$, for several temperatures below, at, and above ${\mathit{T}}_{\mathrm{KT}}$. The temporal evolution of spin configurations was determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by a hybrid Monte Carlo algorithm. The neutron-scattering function S(q,\ensuremath{\omega}) was calculated from the resultant space-time displaced spin-spin correlation function. Pronounced spin-wave peaks were found both in the in-plane and the out-of-plane scattering function over a wide range of temperatures. The in-plane scattering function ${\mathit{S}}^{\mathit{xx}}$ also has a large number of clear but weak additional peaks, which we interpret to come from two-spin-wave scattering. In addition, we observed a small central peak in ${\mathit{S}}^{\mathit{xx}}$, even at temperatures well below the phase transition. We used dynamic finite-size scaling theory to extract the dynamic critical exponent z. We find z=1.00(4) for all T\ensuremath{\le}${\mathit{T}}_{\mathrm{KT}}$, in excellent agreement with theoretical predictions, although the shape of S(q,\ensuremath{\omega}) is not well described by current theory. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9607063•Physical review. B, Condensed matter•Nov 1, 1996•90 citations

Critical dynamics of the body-centered-cubic classical Heisenberg antiferromagnet

Alex Bunker, Kun Chen, David Paul Landau

We have used spin dynamics techniques to perform large-scale simulations of the dynamic behavior of the L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L body-centered-cubic classical Heisenberg antiferromagnet with L\ensuremath{\le}48 at a range of temperatures above and below as well as at the critical point ${\mathit{T}}_{\mathit{c}}$. The temporal evolutions of the spin configurations were determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by Monte Carlo simulations. The neutron scattering function S(q,\ensuremath{\omega}) was calculated from the space- and time-displaced spin-spin correlation function. We used a previously developed dynamic finite-size scaling theory to extract the dynamic critical exponent z from S(q,\ensuremath{\omega}) at ${\mathit{T}}_{\mathit{c}}$. Our results are in agreement with the theoretical prediction of z=1.5 and with experimental results; however, we find that the asymptotic regime was only entered at L\ensuremath{\approxeq}30. In the analysis of the form of the transverse and longitudinal components of S(q,\ensuremath{\omega}) we found that a central diffusion peak appears below ${\mathit{T}}_{\mathit{c}}$ predominantly in the longitudinal component and remains present through and above ${\mathit{T}}_{\mathit{c}}$. The transverse component of the spin-wave peak is Lorentzian below ${\mathit{T}}_{\mathit{c}}$ but for T\ensuremath{\geqslant}${\mathit{T}}_{\mathit{c}}$ is described best by a more complex functional form. Below ${\mathit{T}}_{\mathit{c}}$ we see evidence of multiple spin-wave peaks in the longitudinal component. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Oct 1, 1996•31 citations

SPIN DYNAMICS SIMULATIONS OF CLASSICAL, THREE-DIMENSIONAL HEISENBERG MAGNETS

D. P. Landau, Alex Bunker, KUN CHEN

Spin dynamics simulations have been used to study dynamic critical behavior of classical Heisenberg magnets. The temporal evolutions of the spin configurations were determined numerically from coupled equations of motion by a fourth-order predictor corrector method, with initial spin configurations generated by Monte-Carlo simulations. The neutron scattering function, S (q, ω), was calculated from the space and time displaced spin-spin correlation function and the dynamic critical exponent z was extracted using dynamic finite size scaling theory. For both ferromagnetic and antiferromagnetic models we find good agreement with theoretical predictions and experimental results.

International Journal of Modern Physics C•Jun 1, 1996•3 citations
We calculate the dynamic structure factor S(\mathbitk,\ensuremathω)S(\mathbit{k},\ensuremath{\omega})S(\mathbitk,\ensuremathω) in the paramagnetic regime of quantum Heisenberg ferromagnets for temperatures TTT close to the critical temperature Tc{T}_{c}Tc​ using our recently developed functional renormalization group approach to quantum spin systems. In d=3d=3d=3 dimensions, we find that for small momenta \mathbitk\mathbit{k}\mathbitk and frequencies \ensuremathω\ensuremath{\omega}\ensuremathω the dynamic structure factor assumes the scaling form S(\mathbitk,\ensuremathω)=(\ensuremathτTG(\mathbitk)/\ensuremathπ)\ensuremathΦ(k\ensuremathξ,\ensuremathω\ensuremathτ)S(\mathbit{k},\ensuremath{\omega})=(\ensuremath{\tau}TG(\mathbit{k})/\ensuremath{\pi})\mathrm{\ensuremath{\Phi}}(k\ensuremath{\xi},\ensuremath{\omega}\ensuremath{\tau})S(\mathbitk,\ensuremathω)=(\ensuremathτTG(\mathbitk)/\ensuremathπ)\ensuremathΦ(k\ensuremathξ,\ensuremathω\ensuremathτ), where G(\mathbitk)G(\mathbit{k})G(\mathbitk) is the static spin-spin correlation function, \ensuremathξ\ensuremath{\xi}\ensuremathξ is the correlation length, and the characteristic timescale \ensuremathτ\ensuremath{\tau}\ensuremathτ is proportional to \ensuremathξ5/2{\ensuremath{\xi}}^{5/2}\ensuremathξ5/2. We explicitly calculate the dynamic scaling function \ensuremathΦ(x,y)\mathrm{\ensuremath{\Phi}}(x,y)\ensuremathΦ(x,y) and find satisfactory agreement with neutron scattering experiments probing the critical spin dynamics in EuO and EuS. Precisely at the critical point where \ensuremathξ=\ensuremath∞\ensuremath{\xi}=\ensuremath{\infty}\ensuremathξ=\ensuremath∞, our result for the dynamic structure factor can be written as S(\mathbitk,\ensuremathω)=(\ensuremathπ\ensuremathωk)\ensuremath−1TcG(\mathbitk)\ensuremathΨc(\ensuremathω/\ensuremathωk)S(\mathbit{k},\ensuremath{\omega})={(\ensuremath{\pi}{\ensuremath{\omega}}_{k})}^{\ensuremath{-}1}{T}_{c}G(\mathbit{k}){\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\omega}/{\ensuremath{\omega}}_{k})S(\mathbitk,\ensuremathω)=(\ensuremathπ\ensuremathωk​)\ensuremath−1Tc​G(\mathbitk)\ensuremathΨc​(\ensuremathω/\ensuremathωk​), where \ensuremathωk\ensuremath∝k5/2{\ensuremath{\omega}}_{k}\ensuremath{\propto}{k}^{5/2}\ensuremathωk​\ensuremath∝k5/2. We find that \ensuremathΨc(\ensuremathν){\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\nu})\ensuremathΨc​(\ensuremathν) vanishes as \ensuremathν\ensuremath−13/5{\ensuremath{\nu}}^{\ensuremath{-}13/5}\ensuremathν\ensuremath−13/5 for large \ensuremathν\ensuremath{\nu}\ensuremathν, and as \ensuremathν3/5{\ensuremath{\nu}}^{3/5}\ensuremathν3/5 for small \ensuremathν\ensuremath{\nu}\ensuremathν. While the large-frequency behavior of \ensuremathΨc(\ensuremathν){\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\nu})\ensuremathΨc​(\ensuremathν) is consistent with calculations based on mode-coupling theory and with perturbative renormalization group calculations to second order in \ensuremathϵ=6\ensuremath−d\ensuremath{\epsilon}=6\ensuremath{-}d\ensuremathϵ=6\ensuremath−d, our result for small frequencies disagrees with previous calculations. We argue that until now, neither experiments nor numerical simulations have been sufficiently accurate to determine the low-frequency behavior of \ensuremathΨc(\ensuremathν){\mathrm{\ensuremath{\Psi}}}_{c}(\ensuremath{\nu})\ensuremathΨc​(\ensuremathν). We also calculate the low-temperature behavior of S(\mathbitk,\ensuremathω)S(\mathbit{k},\ensuremath{\omega})S(\mathbitk,\ensuremathω) in one- and two-dimensional ferromagnets and find that it satisfies dynamic scaling with exponent z=2z=2z=2 and exhibits a pseudogap for small frequencies.
We use the functional renormalization approach for quantum spin systems developed by Krieg and Kopietz [Phys. Rev. B 99, 060403(R) (2019)] to calculate the spin-spin correlation function $G(\mathbit{k},\ensuremath{\omega})$ of quantum Heisenberg magnets at infinite temperature. For small wave vectors $\mathbit{k}$ and frequencies $\ensuremath{\omega}$ we find that $G(\mathbit{k},\ensuremath{\omega})$ assumes in dimensions d>2d>2d>2 the diffusive form predicted by hydrodynamics. In three dimensions our result for the spin-diffusion coefficient D\mathcal{D}D is somewhat smaller than previous theoretical predictions based on the extrapolation of the short-time expansion, but is still about 3030%30 larger than the measured high-temperature value of D\mathcal{D}D in the Heisenberg ferromagnet ${\mathrm{Rb}}_{2}{\mathrm{CuBr}}_{4}\ifmmode\cdot\else\textperiodcentered\fi{}2{\mathrm{H}}_{2}\mathrm{O}$. In reduced dimensions $d\ensuremath{\le}2$ we find superdiffusion characterized by a frequency-dependent complex spin-diffusion coefficient $\mathcal{D}(\ensuremath{\omega})$ which diverges logarithmically in d=2d=2d=2, and as a power-law $\mathcal{D}(\ensuremath{\omega})\ensuremath{\propto}{\ensuremath{\omega}}^{\ensuremath{-}1/3}$ in d=1d=1d=1. Our result in one dimension implies scaling with dynamical exponent z=3/2z=3/2z=3/2, in agreement with recent calculations for integrable spin chains. Our approach is not restricted to the hydrodynamic regime and allows us to calculate the dynamic structure factor $S(\mathbit{k},\ensuremath{\omega})$ for all wave vectors. We show how the short-wavelength behavior of $S(\mathbit{k},\ensuremath{\omega})$ at high temperatures reflects the relative sign and strength of competing exchange interactions.
Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L body-centered-cubic classical Heisenberg ferromagnet with L\ensuremath{\le}40 in the vicinity of the critical point Tc{\mathit{T}}_{\mathit{c}}Tc​. The temporal evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method with initial spin configurations generated by Monte Carlo methods. The space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms were calculated to determine the neutron-scattering functions. We developed a dynamic finite-size scaling theory for the neutron-scattering function at Tc{\mathit{T}}_{\mathit{c}}Tc​ and used it to extract the dynamic critical exponent z. Within our resolution limit, the value of z was estimated to be 2.478(28), in excellent agreement with the dynamic scaling prediction, and dynamic scaling was found to be valid for momentum transfer q up to 0.4\ensuremath{\pi} in the (1,0,0) direction.
Spin-dynamics techniques have been used to perform large-scale simulations of the dynamic behavior of the classical Heisenberg antiferromagnet in simple cubic lattices with linear sizes L< 60.L<~60.L< 60. This system is widely recognized as an appropriate model for the magnetic properties of RbMnF3.{\mathrm{RbMnF}}_{3}.RbMnF3​. Time evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using an algorithm implemented by Krech et al., which is based on fourth-order Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical temperature. The dynamic critical exponent was estimated to be $z=(1.43\ifmmode\pm\else\textpm\fi{}0.03),$ which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Direct, quantitative comparisons of both the dispersion curve and the line shapes obtained from our simulations with very recent experimental results for RbMnF3{\mathrm{RbMnF}}_{3}RbMnF3​ are presented.
Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the classical three component XY model (i.e., the anisotropic limit of an easy-plane Heisenberg ferromagnet), on square lattices of size up to 1922{192}^{2}1922, for several temperatures below, at, and above TKT{\mathit{T}}_{\mathrm{KT}}TKT​. The temporal evolution of spin configurations was determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by a hybrid Monte Carlo algorithm. The neutron-scattering function S(q,\ensuremath{\omega}) was calculated from the resultant space-time displaced spin-spin correlation function. Pronounced spin-wave peaks were found both in the in-plane and the out-of-plane scattering function over a wide range of temperatures. The in-plane scattering function Sxx{\mathit{S}}^{\mathit{xx}}Sxx also has a large number of clear but weak additional peaks, which we interpret to come from two-spin-wave scattering. In addition, we observed a small central peak in Sxx{\mathit{S}}^{\mathit{xx}}Sxx, even at temperatures well below the phase transition. We used dynamic finite-size scaling theory to extract the dynamic critical exponent z. We find z=1.00(4) for all T\ensuremath{\le}TKT{\mathit{T}}_{\mathrm{KT}}TKT​, in excellent agreement with theoretical predictions, although the shape of S(q,\ensuremath{\omega}) is not well described by current theory. \textcopyright{} 1996 The American Physical Society.
We have used spin dynamics techniques to perform large-scale simulations of the dynamic behavior of the L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L body-centered-cubic classical Heisenberg antiferromagnet with L\ensuremath{\le}48 at a range of temperatures above and below as well as at the critical point Tc{\mathit{T}}_{\mathit{c}}Tc​. The temporal evolutions of the spin configurations were determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by Monte Carlo simulations. The neutron scattering function S(q,\ensuremath{\omega}) was calculated from the space- and time-displaced spin-spin correlation function. We used a previously developed dynamic finite-size scaling theory to extract the dynamic critical exponent z from S(q,\ensuremath{\omega}) at Tc{\mathit{T}}_{\mathit{c}}Tc​. Our results are in agreement with the theoretical prediction of z=1.5 and with experimental results; however, we find that the asymptotic regime was only entered at L\ensuremath{\approxeq}30. In the analysis of the form of the transverse and longitudinal components of S(q,\ensuremath{\omega}) we found that a central diffusion peak appears below Tc{\mathit{T}}_{\mathit{c}}Tc​ predominantly in the longitudinal component and remains present through and above Tc{\mathit{T}}_{\mathit{c}}Tc​. The transverse component of the spin-wave peak is Lorentzian below Tc{\mathit{T}}_{\mathit{c}}Tc​ but for T\ensuremath{\geqslant}Tc{\mathit{T}}_{\mathit{c}}Tc​ is described best by a more complex functional form. Below Tc{\mathit{T}}_{\mathit{c}}Tc​ we see evidence of multiple spin-wave peaks in the longitudinal component. \textcopyright{} 1996 The American Physical Society.
Spin dynamics simulations have been used to study dynamic critical behavior of classical Heisenberg magnets. The temporal evolutions of the spin configurations were determined numerically from coupled equations of motion by a fourth-order predictor corrector method, with initial spin configurations generated by Monte-Carlo simulations. The neutron scattering function, S (q, ω), was calculated from the space and time displaced spin-spin correlation function and the dynamic critical exponent z was extracted using dynamic finite size scaling theory. For both ferromagnetic and antiferromagnetic models we find good agreement with theoretical predictions and experimental results.