Bose-Einstein Condensation of Confined Magnons in Nanostructures

A number of diverse kinds of magnetic experiments are shown to lead to the observation of a Bose-Einstein condensation of magnons in some nanostructures, especially in Fe and in Co/Pt. The experiments include 1) the measurement of magnetic aftereffect, 2) measurements of Néel’s fluctuation field, including an estimate of activation volume, 3) the measurement of the energy gap which provides a lower limit on the BEC temperature, 4) using Lagrange multipliers in the quantum thermodynamics of magnons, and 5) the observation of a visible anomaly in the Bloch T3/2 law for the temperature dependence of magnetization of nanostructured ferromagnets.


Introduction
The occupation of a single quantum state by a large fraction of bosons at low temperatures [1] was predicted by Bose and Einstein in the 1920s.The quest for Bose-Einstein condensation (BEC) in a dilute atomic gas was achieved [2] in 1995 using laser-cooling to reach ultra-cold temperatures of 20 nK.BEC of dilute atomic gases, now regularly created in a number of laboratories around the world, have led to a wide range of unanticipated applications.Especially exciting is the effort to use BEC for the manipulation of quantum information, entanglement, and topological order.
The study of atomic BEC has yielded rich dividends.A promising extension is to magnons-spin-wave quanta that behave as bosonic quasiparticles-in magnetic nanoparticles.This system has unique characteristics differentiating it from atomic BEC, creating the potential for a whole new variety of interesting behaviors and applications that include high-temperature Bose-Einstein condensation and novel nanomagnetic devices.Bose-Einstein condensation in atoms [2] and photons [3] is now an accepted phenomenon.Recent research [4]- [10] focuses on gaining understanding in different quasi-particles, including photons, phonons, gluons, and magnons.
Magnons are quantized spin waves created when the magnetic moments of an electron spin change due to variation in magnetic field across a crystal lattice.The annihilation and the creation of magnons thus occur when the electron undergoes a spin-flip.The magnetization of the material also depends upon the magnon density [11], which increases or decreases with temperature.The magnon density is also affected by the electron-magnon scattering in the unoccupied region of the transition metals [12].To undergo BEC, it is necessary to either increase the density of bosons or lower the temperature.The consensus on whether these magnons obey Bose-Einstein statistics is still to be reached.Numerous claims were made to prove it [4] [6] [13] [14] but still the true nature of it is yet to be understood.In a quasi-equilibrium state, quasiparticles such as magnons show BEC [15] experimentally.
The BEC of quasi-equilibrium magnons has been achieved at room temperature by exciting more magnons into a gas of quasi-equilibrium magnons with a non-zero chemical potential through microwave pumping.However, some physicists [16] [17] are hesitant to accept the phenomenon observed by Demokritov et al. [6], as BEC and are more willing to perceive it as coherent quantum state or BEC like behavior.Zvyagin [18] observes that "Despite the question of whether Bose-Einstein condensation is the right term to apply to the observed phenomena, coherent quantum states of magnons, either near phase transitions [4] [9] [18] [19], or induced by the paramagnetic pumping [6], are among the most interesting phenomena of modern condensed matter physics".
Magnons and their strong confinement [6] induce gaps in their energy spectra.The lowest lying of the gaps provides a first approximation to the temperature of BEC formation.However, BEC formation is also influenced by the density of the magnons.This modifies the energy gap, which can cause the significant changes in the predicted BEC temperature.The behavior and the magnetic properties of the magnonic crystals vary significantly with the magnon density and temperature [13].For the material to behave as BEC, it must be cooled into its lowest possible energy state.However, if the temperature is increased beyond the gap temperature (the temperature that moves the material from lowest energy state to the next) [20] [21], the behavior changes and it is considered as lower limit on attaining BEC.For a particular nanometer size length of material, the material is in its ferromagnetic state above its gap temperature and is a Bose-Einstein condensate below it.By creating higher magnon density [22], the Bose-Einstein temperature and associated phenomenon could occur at even higher temperatures.
In nanoparticles, finite size and surface effects, anisotropy and influence of geometry, all have significant impact on magnon characteristics, including lifetimes and scattering processes.At nanometers length [6] [23] with different lattice types and particle shapes, the blocking temperature becomes important.If the temperature is increased beyond the gap temperature (the temperature that moves the material from lowest energy state to the next) [24] [25], the behavior changes and it is a lower limit on attaining BEC.For a particular nanometer size length of material, the material is in its ferromagnetic state above its gap temperature and is Bose-Einstein condensate below it.The rough shape of a phase diagram describing the critical behavior of magnons in magnetic nanoparticles, delineating the 1) Bose-Einstein condensate, 2) a superparamagnetic, and 3) ferromagnetic phases in nanostructures has been pictured [26] and shown in Figure 1.
In nanoparticles, finite size and surface effects, anisotropy, and influence of geometry have significant impact on magnon characteristics, including lifetimes and scattering processes.Also, at nanometers length [23] [27] [28] with different lattice types and particle shapes.The material is a Bose-Einstein condensate below the blocking temperature and is in a superparamagnetic state above that.
In a study of magnetic aftereffect we observed non-Arrhenius behavior, in that the decay rates did not increase monotonically with temperature as expected, but instead had a low temperature peak.We analyzed the behavior of the fluctuation fields and again found non-Arrhenius behavior.We attributed the peak to a possible Bose-Einstein condensation.We therefore computed the energy gap in a confined space in order to determine the condensation temperature.Since magnons are bosons, we assumed that they obeyed Bose-Einstein statistics and found that the observed behavior could be explained by a linear variation of the chemical potential with temperature.The temperature intercept of this variation was a positive absolute temperature.Below that temperature the chemical potential is zero.A nonzero chemical potential is possible in quasi-equilibrium, where the number is approximately constant.The slope of the linear portion above T B according to Maxwell's equations is equal to the entropy of the system.The entropy was the same in all the aftereffect experiments, since the samples were all initially DC demagnetized.A further test of this analysis was to measure Bloch's T 3l2 law where we found an anomaly.We again explained this anomaly by Bose-Einstein condensation and got good agreement with experiment.

Non-Arrhenius Behavior in Magnetic Aftereffect
When a ferromagnetic material is magnetized, a reversed field can more or less demagnetize it.Maintaining the field and temperature constant, the magnetization then decays slowly (magnetic aftereffect) towards the anhysteretic ground state.It is expected, and it is found in many magnetic materials, that the decay of magnetization, measured by thermal magnetic aftereffect experiments, becomes smaller when the temperature is cooled.This Arrhenius behavior is a consequence of Maxwell-Boltzmann statistics.In some nanostructured ferromagnets, a peculiar temperature dependence of the thermal magnetic aftereffect has been reported [29], namely an increase in decay rate as the temperature is lowered from room temperature, peaking at some low temperature.This non-Arrhenius behavior could be explained by assuming that the magnons responsible for the decay obeyed Bose-Einstein (BE) statistics with a non-zero chemical potential, and that the peak is connected with the Bose-Einstein condensation temperature, T B .In a classic paper on magnetic aftereffect [30], the Preisach-Arrhenius (PA) model is an approach built on fundamental thermodynamic principles and predicts the decay rate at any holding field very accurately.The PA model uses the Arrhenius-Néel law [31] to determine the aftereffect.

Decay Rates
Magnetic aftereffect, the slow decay of the magnetization toward the anhysteretic state, is a non-equilibrium process that is controlled by 1) overcoming a distribution of energy barriers by thermal activation, and by 2) the chemical potential, ζ, that determines how many thermal quanta are available.The former is predominant in the remanence state and the latter in the demagnetized state.Arrhenius behavior implies a monotonic increase in the decay rate with increasing temperature.Normally, this is just what is seen in most materials.However, when the magnon energy distribution, which obeys quantum statistics, has a finite Bose-Einstein condensation at a BEC temperature, T B , a peak is observed in the maximum decay rate vs. temperature at T B .Such a peak has been observed in some nano-magnetic materials.For example, decay rate studies of the temperature dependence of the magnetic aftereffect of Fe nano-particles [32] shows a large peak at ∼30 K, as seen in Figure 2 for nano-particle powder.Assuming [33] that the magnons responsible for the decay obeyed Bose-Einstein (BE) statistics with a linear temperature-dependent chemical potential, ζ, obeying [34] , where H, M, S, v, p, and T are the magnetic field, the magnetization, the entropy, the volume, the pressure, and the temperature, respectively.Thus, under constant pressure and magnetization ( ) for 0 for Unless S ≡ 0, ζ cannot be zero for all temperatures, and it is well known that the entropy is not zero, at least in the demagnetized state.For any given magnetic state, the entropy, S, is a constant.Assuming a chemical potential variation from Equation ( 2) with two adjustable parameters, S and T B , for the Fe nanoparticles yields Figure 3, which shows the chemical potential needed to provide the fit shown in Figure 2. Similar behavior in metal particle recording tape is shown in Figure 4 and the corresponding chemical potential is shown in Figure 7.
To understand fully this non-Arrhenius behavior, it is useful to elaborate their thermodynamic framework.For simplicity, assume that the basic particles are all identical and perfectly ordered.If the aftereffect is measured at each temperature, by starting from the same magnetic state, then the distribution will have the same entropy, and thus the same chemical potential.Hence, the slope of the chemical potential with temperature should be constant, as found experimentally.
Similar results were found in nano-particulate (metal-particle tape) media [35].In this case, as in the nanoparticles, we found that the maximum decay rate, R, had a low-temperature peak in contrast to the expected monotonic Arrhenius behavior, as seen in Figure 4. Again, except at the lowest temperatures, the chemical potential is linear with temperature (Figure 7).Note that [36], in neither the nanograin iron, nor the nano-particulate media, is there any other explanation for the peak, i.e., there are no phase changes, no changes in the easy direction of magnetization, no movement of interstitial atoms, or other significant crystallographic or magnetic changes.
The magnetic aftereffect can be modeled as due to interactions between magnons (spin waves) and hysterons (magnetic entities acting as a single magnetic moment).If the energy of the magnon is less than that of the barrier energy for a reversal of the magnetization, then the collision is purely elastic.On the other hand, if the magnon has greater energy than the barrier energy, then there is a finite probability that the hysteron will reverse.The frequency of particle reversal, f, due to magnons is the product of the attempt frequency f A and the probability that the magnon has sufficient energy to overcome the energy barrier to reversal, p(E > E B ).If the energy distribution obeys Maxwell-Boltzmann (MB) statistics and the density of states in energy space is uniform, then the rate equation is given by the Arrhenius law.On the contrary, if the energy distribution obeys quantum statistics or the density of states in energy space is not uniform, non-Arrhenius behavior would be observed.In general, the number of particles that in a given collision have energy greater than E B is obtained by integrating over all possible energies, that is where n(E) is the expectation value of the number of particles in a state that has energy E, and δ(E) is the density of states in energy space.

Fluctuation Fields
A different type of decay was observed in measurements from 4 K to room temperature [31] on a film that was nominally a nanostructured (3 nm Co/2 nm Pt) [15] magnetic multilayer film grown on a sapphire substrate at 573 K. Analysis of this film showed that it actually had a nanopillar structure [37] of about 30 nm height and a similar diameter.The magnetization is perpendicular to the film and arises from these aligned nano-pillars.The magnetization decay rate in the cobalt-platinum is much faster than those of the aforementioned nano-particulate systems.In this case, it is possible to view the entire sigmoidal shaped decay curve.Such decay curves, in log time, are shown in Figure 5 for various fields near the coercivity.This permits an alternate method of analysis: characterizing the aftereffect in terms of the fluctuation field introduced by Néel [33].This experiment provides an especially robust observation of the BE condensation of the magnons.The fluctuation field, H f , as expanded by Street and Brown [38], can be viewed as the driving force in a magnetic aftereffect.It is normally given by 0 where V is the activation volume of the individual magnetic entity ("particle") and M s is the saturation magnetization.Many investigators incorporate the chemical potential term into the attempt frequency, hence they find the activation volume is a significant function of temperature, which is contrary to having a given physical volume [39].
The technique to obtain the fluctuation field from the decay curves involves scaling and translating them.If there is a negligible reversible magnetization and negligible inter-particle interactions, then it can be shown [40] that all the aftereffect curves and the major loop are similar in a geometric sense.The similarity of the major loop with the aftereffect curve is a direct method of confirming the Wohlfarth hypothesis since the field that maximizes the slope of the major loop, which is the susceptibility, also maximizes the slope of the aftereffect decay.
Since this Co/Pt satisfied the criteria in the previous paragraph, a new technique was developed [41] to measure the fluctuation field that consisted of finding two parameters: a horizontal translation and a scaling to account for the difference in anhysteretic magnetization at different holding fields.
The technique is performed by mapping each of the curves onto the major hysteresis loop, in order to normalize them so that they all extend from positive saturation to negative saturation.Since the horizontal translation is directly proportional to the ratio of the holding field to the fluctuation field, it is possible to determine the fluctuation field from the slope of the translation versus the holding field, as obtained from a linear regression of the data in Figure 6.
The fluctuation field at each temperature is obtained from the slope of the required translation as a function of the holding field.As seen in Figure 7, the relaxation does not slow down as the temperature is lowered but, instead increases and reaches a peak at 14 K, which we identify as a Bose-Einstein condensation temperature, T B .The fluctuation field is theoretically (under some simplifying assumptions) directly proportional to the absolute temperature below the BE condensation temperature, T B , of the magnons (see Figure 9).Below T B the fluctuation field mimics the Arrhenius law.The field has a sharp peak at the T B , then decreases until it approaches a constant value asymptotically from below for higher temperatures.The chemical potential, ζ, obtained from the fit shown in Figure 8 using the data in Figure 9.
Note that ζ = 0 below T B (which is a requirement for having a BEC), and linear above.As already noted, these peaks in the magnetic aftereffect present particularly robust indications of a Bose-Einstein phase transition.
For a given material, there is a critical condensation temperature below which there is no peak in the fluctuation field.This is illustrated in Figure 10, where this temperature is 10 K.When the condensation temperature is above this value, the fluctuation field increases linearly, i.e., following the Arrhenius law, until that temperature and above it, the fluctuation field decreases as the temperature continues to increase.This calculation assumes that the energy barrier is temperature independent.For a temperature-dependent energy barrier, the curve will eventually be monotonically increasing.

Energy Gap in a Confined Space
There have been several discussions related [42]- [44] to the energy gap in a finite size material.The gap tem-     perature, T gap , is the lower bound for the BEC temperature.This means that in Figure 3, the line separating the BEC from the ferromagnetic region can be substantially higher.
The energy E i of a magnon in state i, for a simple body-centered or face centered cubic lattice, is given by [45] where J is the exchange constant, k i is the wave vector associated with state i, and a is the lattice spacing, that is, the distance between corner atoms in the unit cell.The wave vector can be expressed in terms of the de Broglie wavelength, λ i , as The smallest non-zero wave vector is obtained by setting λ i equal to the largest dimension of the particle.The BEC temperature is estimated by computing the gap temperature, T gap .This temperature is a lower bound on T B , which could be appreciably larger.It is obtained by setting the lowest non-zero magnon energy equal to the thermal energy.Using the equipartition of energy, this energy is 1/2k B T, where k B is Boltzmann's constant and T is the temperature in K, the thermal de Broglie wavelength for a magnon in a spherical particle 2 2π and the gap temperature, T gap , is found by setting λ equal to the particle's diameter, d, and solving for T we obtain T gap , which is a lower bound T B .Thus, For iron, using the lattice constant a of 0.4 nm, and the accepted [46] J value of 2.16 × 10 -21 J, this would become 12,270 (a/d) 2 .These parameters, J and a, are approximately characteristic of the iron and cobalt nanostructures.

Quantum Thermodynamics of Magnons
Magnons are quasi-particles that interact with the spin of electrons to change the magnetization of the material.In order to conserve angular momentum, since the electrons have spin 1 2 , the difference between two electron states must be an integer.Therefore, magnons are bosons.When the magnons are in equilibrium, they obey Bose-Einstein statistics with zero chemical potential, but when they are in quasi-equilibrium, but not in true equilibrium.That is when the number of magnons is essentially constant, they obey Bose-Einstein statistics with nonzero chemical potential, ζ.The number of magnons can be kept essentially constant if new magnons are created at the same rate at which they are destroyed.This can be described by noting that the magnons are in quasi-equilibrium.
The chemical potential has to be adjusted so that the sum of ( 1) over all states is equal to the total number of particles.In order to find the quasi-equilibrium state, we have to find the state with the maximum thermodynamic probability, that i the entropy, subject to the constraints that the total energy and the number of particles is constant.This is generally solved using Lagrange multipliers so that the problem can be solved as an unconstrained maximization.The Lagrange multiplier that keeps the number of magnons constant is the chemical potential.
Thus, the number of magnons with energy E i , is given by Then the chemical potential is chosen by letting where N is the total number of particles.Under quasi-equilibrium, ζ is negative and its magnitude grows larger with increasing temperature.In that case, since E i is positive, as the temperature increases, the statistics approach classical statistics as in the correspondence principle, since the exponential term is so large that the −1 can be neglected.Another way of looking at it is, if the exponential term in the denominator is very large, the magnon wavelength is small and so there is negligible overlap of the wave functions, and it is possible to distinguish the particles without violating the exclusion principle.
Since the chemical potential is a negative number, it has the effect of reducing the number of magnons in all energy states.This is necessary in order to explain the extremely slow decay of aftereffect at room temperature in the presence of an attempt frequency of the order of 10 9 .
Kittel has shown [47]- [50] that the number of magnons at any temperature is directly proportional to the difference between the saturation magnetization occurring at that temperature to that at 0 K.For ferromagnetic materials such as iron, this is only a few percent of the saturation magnetization at room temperature.Furthermore, since the change in magnetization with time in these materials is extremely slow, the number of magnons is essentially a constant.While T B can be made larger by increasing the magnon density, this review only ex-amines the theoretical dependence of T B on the nanoparticle size, which is easily accessible experimentally.
Magnons (spin-wave quanta) behave like a gas of weakly interacting particles obeying Bose-Einstein statistics.When confined in a small nanoparticle (e.g., an iron sphere with diameter on the order of 1 -10 nm), the allowed magnon wave vectors k i are quantized, giving a discrete spectrum E i with an energy gap, E gap = E 1 − E 0 , between the lowest and next to lowest allowed states.If the magnons are in equilibrium at temperature T, then the number in a state i, is given by the Bose-Einstein distribution function.

Quasi-Equilibrium-Number of Magnons
In measurements of thermal magnetic aftereffect, the magnetization decay may be so small, that for some measurements the decay is only a small fraction of a per cent.Thus, when measuring the magnetization as a function of temperature, at any given temperature, the magnetization is essentially constant.Such a system is said to be in quasi-equilibrium.Since the fractional change in magnetization is equal to the number of magnons divided by the number of electrons, the number of magnons in these cases is also constant.Thus, when measuring the saturation magnetization as a function of temperature, as in the case of the Co/Pt nanostructure, the magnetization is not a measurable function of time and thus, the number of magnons is constant.In the case of the decay measurements, if the observation window is short enough, the magnetization changes almost imperceptibly, and again the number of magnons is constant.When maximizing the thermodynamic probability, the method of Lagrange multipliers is used to permit only variations, if the number of particles is constant, and ζ is that multiplier.
In an open system, or in one that bosons can be created or destroyed, the chemical potential must be zero.However, if the system is far from equilibrium, but is changing sufficiently slowly so that it reaches a quasiequilibrium, a non-zero chemical potential is permitted.Can quasiparticles, such as magnons, whose numbers are not conserved but have finite lifetimes undergo Bose-Einstein condensation?The consensus on this issue [16], reached about a decade ago, is that if the lifetime of the particles is much longer than the time they need to scatter with each other, condensation is possible.It is thus appropriate to say that if the observations do not include all possible tests, the system is either BEC or BEC-like.

Anomaly in the Bloch T 3/2 Law in Nanostructures
Bloch's [48] T 3/2 law of variation of saturation moment with temperature [7] is regarded as the fundamental law of ferromagnetism.The decrease of the saturation moment with increasing temperature in a bulk ferromagnetic material is determined by the "stiffness" of the exchange coupling which tends to align the spins of all the atoms parallel [8].The law has been tested in a careful set of experiments [9] in single crystals of ferromagnetic metals.
The law assumed that the chemical potential was zero.A correction to the chemical potential variation, characteristic of BE condensation, suggested above showed that the magnetization obeyed the Bloch's T 3l2 law up to the condensation temperature, but above that temperature it is asymptotic to another T 3l2 law, with different coefficients.
The fractional decrease in saturation magnetization, M, is 3 2 which is the same as the traditionally derived case discussed by Kittel [47].If S = 0 then, ζ = 0, as in perfect single crystals, this becomes ( ) Measurements in nanostructures [4] showed that this correction, illustrated in Figure 11, is accurate virtually to room temperature, further demonstrating that these confined magnons form a Bose-Einstein condensate.

Results
Bose-Einstein condensation of magnons confined in nanostructures has been established by a variety of different experiments all giving the same results for a particular specimen.These results are obtained in quasi-equilibrium.For a typical magnon, magnon-phonon processes are slower than magnon-magnon processes.This means that spin-spin equilibrium is established more rapidly than spin-lattice equilibrium [52].
Of special note is the identification of the peak seen in thermal magnetic aftereffect versus temperature.As the temperature is lowered, "The elegant and wide-ranging Arrhenius equation tells us how much slower" [53] Atzmony, et al. [12], noticed the peaks and said that they were non-Arrhenius, but had no explanation.A number of other reports of similar peaks were never explained or even commented on.The Arrhenius equation is based on classical thermodynamics.Since the magnons are bosons, Bose-Einstein thermodynamics provides excellent agreement between experiment and theory [4]

Discussion
The study of atomic BEC has yielded rich dividends.A promising extension is to magnons-spin-wave quanta that behave as bosonic quasiparticles-in magnetic nanoparticles.This system has unique characteristics differrentiating it from atomic BEC, creating the potential for a whole new variety of interesting behaviors and applications that include novel nanomagnetic devices.
Further work is needed to determine the relationship between nanostructural morphology (size, shape, etc.) and the Bose-Einstein condensation temperature.
A most important point is that a magnon propagates spatially all over the magnet.By this propagation, quantum coherence is established between spatially separated points.Therefore, by exciting a macroscopic number of magnons, one can easily construct states with huge entanglement [54].Entanglement is a basis for quantum computation.
The BEC has important significance to magnetic resonance: Below T B , the entropy is zero, and is essentially constant above that.If the Gilbert damping is related to the entropy, then very narrow line widths are possible below T B .This has implications for low loss or sharp filters for microwave devices.
A condensate is a quantum system.It follows, therefore, that its properties depend in large part on the measurements you can perform on it.Condensates are open systems; some dynamical processes cause members of the condensate to be lost while other processes add particles to it.This does not imply a zero chemical potential, e.g., introductory chemistry students routinely calculate the chemical potential for the reaction: hydrogen oxygen water + ⇔ [55].Since BE condensates are not made in pure states, at the most fundamental level, there should be a highly entangled state of the particles forming the condensate plus the environment.

Conclusion
Bose-Einstein (or Bose-Einstein-like) condensation in metallic nanostructures has been observed using a variety of magnetic signatures [4] [11] [15] [26] [29] [30] [32] [33] [35]- [37] [39]- [41].All of the methods gave comparable results.Especially illuminating was the measured graph of chemical potential versus temperature, which in some alloys was exactly as statistical mechanics text books described the ideal case of BEC, i.e., zero chemical potential up to the BEC transition point, and then negative and linear above that temperature.

Figure 1 .
Figure 1.Approximate shape of a phase diagram delineating the boundaries of the regions.

Figure 2 .
Figure 2. Maximum decay rate, R, vs. temperature for an Fe nano-particle powder, showing a peak at ∼30 K.The solid line is a fit to the theory.

Figure 3 .Figure 2 .
Figure 3.Chemical potential required to obtain the fit of R shown in Figure 2. The curvature near ζ = 0 is attributed to size dispersion.

Figure 4 .
Figure 4.The maximum decay coefficient as a function of temperature for the particulate iron in a recording tape media.The dotted line is a plot of the decay assuming classical statistics and using measured values of the parameters.The inset shows a typical measurement (M vs. T) of the decay coefficient.

Figure 6 .
Figure 6.The translation, in log time, versus the holding field, in mT, from the CoPt nanopillar data of Figure 5.The solid line is a linear regression of the data.

Figure 7 .
Figure 7. Variation of ζ/k, in units of 10 3 K, for the recording tape.

Figure 8 .
Figure 8.The translation in log time, versus the holding field, in mT, from the Co/Pt nanopillar data of Figure 6.The solid line is a linear regression of the data.

Figure 9 .
Figure 9. Temperature variation of the fluctuation field, H f , in a Co/Pt nanostructure.

Figure 10 .
Figure 10.The effect of TBE on the fluctuation field, when H f and E B are constant.