Abstract
We report on the observation of the size effect of thermal magnetization in Au nanoparticles. The thermal deviation of the saturation magnetization departs substantially from that predicted by the Bloch T^{3/2}law, indicating the existence of magnetic anisotropic energy. The results may be understood using the uniaxial anisotropy Heisenberg model, in which the surface atoms give rise to polarized moments while the magnetic anisotropic energy decreases as the size of the Au nanoparticles is reduced. There is a significant maximum magnetic anisotropic energy found for the 6 nm Au nanoparticles, which is associated with the deviation of the lattice constant due to magnetocrystalline anisotropy.
Keywords:
Nanoparticles; Magnetic anisotropy; Magnetic properties; Spin wavesIntroduction
Metal nanoparticles of Pd, Au, and Cu have been extensively studied, because, due to a reduction in dimensionality, their ferromagnetic polarizations are quite different from those observed in transition metals [16]. The most frequent effects of the small size are lattice rearrangement, crystalline imperfections, a higher degree of localization, and narrowed valence band width. It has been reported in previous studies [2,4] that individual Pd and Au nanoparticles may reach their ferromagnetic moment at low temperatures, and that, theoretically, there may be a slight enhancement of the 4d localization, although Pd and Au are both characterized by diamagnetism in the bulk state. Bulk Au metal also demonstrates a typical diamagnetic response of −1.42 × 10^{−6} emu/g [7], when the [Xe]4f^{14}5d^{10}6 s^{1} Au configuration has a closed d shell and a single s electron. Finitesize effects play a dominant role in determining the magnetic properties. A decrease in size can lead to unusual ferromagnetic and diamagnetic properties. The origin of the ferromagnetism observed in filled 4d or 5d electron nanoparticle systems can be explained as due to giant magnetic anisotropy [8] and Fermihole effects [9] that influence the evolution from the surface polarization spins to the diamagnetic bulk state. In this letter, we discuss the effects of surface polarization and weak magnetic anisotropy in Au nanoparticles, which indicate the appearance of ferromagnetic spin polarization and magnetic anisotropic energy at low temperatures. Moreover, the strain induced by the lattice can be used to tune the magnetic anisotropic energy, which is obtained from the quantum spin wave theory and the anisotropic Heisenberg ferromagnetic model.
Experimental Details
The Au nanoparticles used in the present study were fabricated by the thermal evaporation method. Highpurity gold ingots (99.999%) were evaporated in the range of 0.1–2 T. The Ar gas was fed at a rate of ~0.1 Å/s. To avoid contamination by magnetic impurities originating from the stainless steel plate the samples were collected by a rotating silicon substrate maintained at the temperature of liquid nitrogen. The resultant samples consisted of collections of individual Au nanoparticles in the form of dried powder. The morphology and structures of the prepared nanoparticles were then characterized using transmission electron microscopy (TEM, JEM1400 JEOL).
Results and Discussion
Structural Analysis
It is clearly evident in the portion of the TEM images shown in Fig. 1a that the nanoparticles are spherical and well separated. The interconnecting nanoparticles are stuck together in clusters due to electrostatic effects as well as an artifact of the drying of aqueous suspensions. The size and distribution of the nanoparticles can be calculated. An examination of the portion of the TEM image shown in the Fig. 1b clearly shows that size and distribution are quite asymmetric and can be described using a lognormal distribution function. The lognormal distribution is defined as follows: where <d> is the mean value and σ is the standard deviation of the function. The mean diameter and standard deviation obtained from the fits are <d> = 3.7(9) nm and σ = 0.251, respectively. The small standard deviation (σ < 0.5) of the function indicates that the distribution is confined to a limited range. The broadening of the width of the distribution profile is due to crystalline and nanoparticle aggregation effects. The electron diffraction pattern corresponding to a selected area in the 3.7(9) nm Au nanoparticles is shown in Fig. 1c and clearly reveals the crystalline nature of the sample. The pattern of the main spots can easily be indexed as basically cubic in structure with a space group of Fm−3m and a lattice parameter of a = 4.07(4) Å. This is consistent with earlier data for bulk Au [7]. The diameters of the nanocrystals as determined from TEM images of the samples used in this study were approximately 3.7(9), 4.3(6), 5.6(4), 6.0(3), and 7.9(1) nm.
Figure 1. a TEM images of Au nanoparticles; b size distribution obtained from a portion of a TEM image of Au nanoparticles, which can be described using a lognormal distribution function, as indicated by the solid curve; c electron diffraction patterns of a selected area 3.7(9) nm Au nanoparticles, revealing the cubic structure of the Au nanoparticles
Magnetization
Magnetization measurements were performed using the conventional superconducting quantum interference device (Quantum Design, MPMS5) set up with magnetic fields from −5 to 5 T, covering a temperature range from 5 to 300 K. The Au nanoparticle sample was mounted in a sample holder capsule. Figure 2a shows the applied fields and the resultant magnetization of Au nanoparticles (with a mean diameter of 3.7(9) nm) obtained at the eight temperatures. When a lower field is applied the magnetization increases rapidly with the field; the increase follows a curved path, revealing that the magnetization follows a Langevin profile. At high temperatures, magnetization saturation is reached at around H_{a} = 0.4 T which is followed by a highfield linear decrease when the applied field H_{a} reaches 1.2 T. There are no significant differences in the magnetization measurements between the fieldincreasing and the fielddecreasing loops found above 25 K, which are consistent with the superparamagnetic behaviors. Figure 2b shows the representative M(T) curves taken for the 3.7(9) nm Au nanoparticle assemblies at the selected applied magnetic field H_{a} = 0.5 T, revealing a superparamagnetic behavior that can be described by temperature dependence Langevin function. The resultant fitting parameters are shown in Fig. 2b. This suggests the existence of two different magnetic components in the sample—a superparamagnetic and a diamagnetic component. A representative hysteresis loop taken at 5 K is shown in Fig. 3. A distinguishable asymmetric coercivity H_{c} = 175 Oe can be observed in the low H_{a} regime, which signals the existence of ferromagnetic spin in the 3.7(9) nm Au nanoparticles. The value obtained for coercivity is close to other previously published values for similar Au nanoparticles [10]. The asymmetric characteristics are assumed to originate from competition between the unidirectional and uniaxial anisotropy [11,12].
Figure 2. a Effects of the various temperatures plotted in relation to the magnetization. The solid lines represent the fitted results. b Temperature dependence of the M(T) curve taken for the 3.7(9) nm Au nanoparticle assemblies at the selected applied magnetic field H_{a} = 0.5 T
Figure 3. Magnetization loops at 5 K for the 3.7(9) nm Au nanoparticles revealing the appearance of magnetic hysteresis in the low field regime
Consequently, we can describe the superparamagnetic system using a Langevin function in combination with a linear component associated with diamagnetism [13,14]. The resultant total magnetization can be expressed as
Here L(x) = coth(x)−1/x is the Langevin function, M_{s} is the saturation magnetization, k_{B} is the Boltzmann’s constant, and χ_{D} is the diamagnetic susceptibility term. The analysis relevant to Eq. 1 is based on a model which ignores the interparticle interactions and the contributions of the distributions of the magnetic moment due to the lognormal size distribution of the nanoparticle system [15]. It can be seen that the fitted curves (solid line in Fig. 2a) are quite consistent with the experimental data. The mechanism often invoked to explain the occurrence of surfacespin polarization effects in nonmagnetic particles [4] is that the shell of the particle is ordered as a ferromagnetic shell, while the core of each Au nanoparticle still behaves as a diamagnetic single domain. Indeed, there is a discrepancy between the data and the Langevin profile shown in the M(H) curves taken at the low field regime. One possible cause of this difference in the fit is the production of a nonmagnetic surface layer by the chemical interaction between the particle and the oxidation. In light of the results obtained in various studies [16], we believe that this difference has a different origin. An alternate explanation has been made by Berkowitz [17,18], who attributed the reduction in the expected magnetization at low temperature to difficulty in reaching saturation, because of a large surface anisotropy.
Magnetic Anisotropic Energy Gap
The thermal deviation of the saturation magnetization can be used to identify the anisotropic energy gap. Figure 4 shows the dependency of the thermal energy of the thermal deviation ΔM_{s}(T) = [M_{s}(5K)−M_{s}(T)]/M_{th} on the saturation magnetization, as obtained from the fitting of Eq. 1, where M_{th} is the saturated magnetization taken at room temperature. The ΔM_{S}(T) curve follows the Bloch T^{3/2}law (dashed line) expected for ferromagnetic isotropic systems below 10 K [14,19] but departs from the curve in the high thermal energy regime, signaling the onset of magnetic anisotropy [15,2023], presumably due to the high surfacetovolume ratio of the nanoparticles. The discrepancy of fit above 10 K (for low applied fields) may be associated with the effects of uniaxial anisotropy [2022] and with inhomogeneities in the magnetic moments [23], which have been ignored in the above analysis.
Figure 4. Plot of the dependency of the thermal energy on the saturation magnetization M_{s}(T) together with the thermal deviation ΔM_{s}(T) due to the saturation magnetization. The solid lines represent the fitted results
Here, we consider the surface and anisotropic effects in a quantum spin wave model for the Heisenberg ferromagnetic model [24,25]. We can incorporate the spin–spin effects and anisotropy between coupling constants, that are known to be important in a nanosize system, into the anisotropic Hamiltonian, but do not include diamagnetic effects [2628]
where the sum in the first term is the anisotropic ferromagnetic Heisenberg exchange interaction (J_{z} and J_{xy}) between nearestneighbor spins on a nanoparticle; S denotes the spin component of the electrons; and the last part corresponds to the Zeeman energy (mS is the magnetic moment per atom). The theory and the method of calculation have already been described in detail elsewhere [2628], therefore only a few basic steps will be given here. We utilize an external perturbation method and calculate the energy in the ground state of the spin wave with wave vector q and dispersion relation ε_{0}(q). We can now rewrite the equation as
where δ is the nearestneighbor vector. The lattice constants for Au facecentered cubic nanoparticles with Fm−3m symmetry are a = 4.07(4) Å, with a nearestneighbor spacing of The dispersion relation can be rewritten as
This anisotropy in the coupling constants produces an energy gap in the spin wave spectrum of Δ = 12S(J_{z} − J_{xy}). The gap leads to an exponential dependence of the order parameter on the thermal energy k_{B}T:
The solid lines indicate the results from the fit of Eq. 5; the fitting parameters are listed in Table 1; the energy gap obtained from the fit is plotted with the diameter in the Fig. 5a (right panel). In the case of 3.7(9) nm, at higher thermal energy k_{B}T ~ 8 meV, the monotonic change of ΔΜ_{s} is closed to one and will be overcame by the thermal energy. The direction of the magnetization of each Au nanoparticle simply follows the direction of the applied magnetic field. Consequently, the magnetization becomes superparamagnetic and shows paramagnetic properties. The larger the size of the nanoparticle, the higher the magnetic anisotropic energy, which therefore increases with increasing particle size, until reaching the maximum magnetic anisotropic energy: Δ = 6.527 meV in the 6.0(3) nm Au nanoparticles. The results are in good agreement with the molecular field theories, which predict linear or exponential variations for large and small anisotropic energies, depending on whether a classical or quantized system is used for the magnetic moment [29]. In general, magnetic anisotropy means the dependence of the internal energy of a system on the direction of the spontaneous magnetization. Most kinds of magnetic anisotropy are related to the deviations in the lattice constant of the strain, known as magnetocrystalline anisotropy [30]. Figure 5 shows the strain as a function of mean diameter <d>. Shown in Fig. 5b (left panel), the relative strain can be estimated from the change in the aaxis lattice constant of Au nanoparticles
where a and a_{0} (4.076 Å for bulk Au) indicate the lattice constants of the strained and unstrained crystal, respectively. In general, the spin–orbit interaction will induce a small orbital momentum, which couples the magnetic moment to the crystal axes. In a negative strained nanocrystalline system, the wavefunctions between neighboring atoms will overlap and reduce the magnetic anisotropy. A reduction in the size of the nanoparticles (ε < 0) results in unit cell contraction, which increases the stability of the higher symmetry lattice and the coupling strength of wavefunctions, shown in Fig. 5b. In a positive strained ε > 0 nanocrystalline system, shown in Fig. 5c, the lattice expands and decreases the coupling of wavefunctions. In the case of our Au nanoparticles, the maximum positive deviation in strain was observed when the mean size was 6.0(3) nm. The tendency of strain of size effects was similar with the results of anisotropic energy. However, one possible explanation for the higher strain state accompany with higher magnetic anisotropy energy is an indicative of lattice and magneticanisotropy for Au nanoparticles.
Table 1. Summary of the size and fitting results for Au nanoparticles
Figure 5. a Plots of the variation in Δ and strain ε with mean diameter, revealing the increase in magnetic anisotropic energy with increasing particle size. b Schematic plots for negative strain (ε < 0), and c positive strain (ε > 0)
Conclusions
An analysis of the results leads to an interesting conclusion: that nanosized transition metal Au particles exhibit both ferromagnetism and superparamagnetism, which are in contrast to the metallic diamagnetism characteristic of bulk Au. The superparamagnetic component of Au nanoparticles shows an anomalous temperature dependence that can be well explained by the modified Langevin function theory. Weak magnetic anisotropy was observed in the mean deviation magnetization. The energy of the magnetic anisotropic can be determined from the fitting of the anisotropic Heisenberg model and related with the change of strain. One possible explanation for the origin of the observed superparamagnetic component of the magnetization would be the existences of nonlocalized holes and charge transfer which would signify that deviation from stoichiometry would make only a small paramagnetic contribution to the magnetization [31].
Acknowledgments
We appreciate the financial support of this research from the National Science Council of the Republic of China under grant No. NSC972112M259004MY3.
References

Teranishi T, Kiyokawa I, Miyake M:
Adv. Mater.. 1998, 10:596.
COI number [1:CAS:528:DyaK1cXjsl2nurc%3D]
Publisher Full Text 
Hori H, Teranishi T, Nakae Y, Seino Y, Miyake M, Yamada S:
Phys. Lett. A. 1999, 263:406.
COI number [1:CAS:528:DC%2BD3cXnvFer]; Bibcode number [1999PhLA..263..406H]
Publisher Full Text 
Li WH, Wu SY, Yang CC, Lai SK, Lee KC, Huang HL, Yang HD:
Phys. Rev. Lett.. 2002, 89:135504.
Bibcode number [2002PhRvL..89m5504L]
PubMed Abstract  Publisher Full Text 
Yamamoto Y, Miura T, Suzuki M, Kawamura N, Miyagawa H, Nakamura T, Kobayashi K, Teranishi T, Hori H:
Phys. Rev. Lett.. 2004, 93:116801.
COI number [1:STN:280:DC%2BD2cvotVGitQ%3D%3D]; Bibcode number [2004PhRvL..93k6801Y]
PubMed Abstract  Publisher Full Text 
Hori H, Yamamoto Y, Iwamoto T, Miura T, Teranishi T, Miyake M:
Phys. Rev. B. 2004, 69:174411.
Bibcode number [2004PhRvB..69q4411H]
Publisher Full Text 
Crespo P, Litrán R, Rojas TC, Multigner M, de la Fuente JM, SánchezLópez JC, García MA, Hernando A, Penadés S, Fernández A:
Phys. Rev. Lett.. 2004, 93:087204.
COI number [1:STN:280:DC%2BD2cvotVKqsw%3D%3D]; Bibcode number [2004PhRvL..93h7204C]
PubMed Abstract  Publisher Full Text 
Phys. Rev.. 1931, 38:2051.
COI number [1:CAS:528:DyaA38XhtV2kuw%3D%3D]; Bibcode number [1931PhRv...38.2051S]
Publisher Full Text 
Hernando A, Crespo P, García MA:
Phys. Rev. Lett.. 2006, 96:057206.
COI number [1:STN:280:DC%2BD287gtF2mtQ%3D%3D]; Bibcode number [2006PhRvL..96e7206H]
PubMed Abstract  Publisher Full Text 
Phys. Rev. B. 1988, 37:745.
Bibcode number [1988PhRvB..37..745H]
Publisher Full Text 
Garitaonandia JS, Insausti M, Goikolea E, Suzuki M, Cashino JD, Kawamura N, Ohsawa H, de Muro IG, Suzuki K, Plazaola F, Rojo T:
Nano Lett.. 2008, 8:661.
COI number [1:CAS:528:DC%2BD1cXptl2jsg%3D%3D]; Bibcode number [2008NanoL...8..661G]
PubMed Abstract  Publisher Full Text 
J. Phys. Condens. Matter. 2007, 19:236225.
Bibcode number [2007JPCM...19w6225C]
Publisher Full Text 
J. Appl. Phys.. 2001, 89:3403.
COI number [1:CAS:528:DC%2BD3MXhsFGgtLs%3D]; Bibcode number [2001JAP....89.3403X]
Publisher Full Text 
J. Appl. Phys.. 1959, 30:S120.
Bibcode number [1959JAP....30S.120B]
Publisher Full Text 
Kittel C: Introduction to solid state physics. Wiley, New York; 1996.

Chantrell RW, Popplewell J, Charles SW:
IEEE Trans. Magn.. 1978, MAG14:975.
COI number [1:CAS:528:DyaE1cXlvFSjt7c%3D]; Bibcode number [1978ITM....14..975C]
Publisher Full Text 
IEEE Trans. Magn.. 1981, MAG17:2923.
COI number [1:CAS:528:DyaL38XitV2gsbs%3D]; Bibcode number [1981ITM....17.2923P]
Publisher Full Text 
Du YW, Lu HX, Wang YQ, Wang TX:
J. Magn. Magn. Mater.. 1983, 31:896.
Bibcode number [1983JMMM...31..896D]
Publisher Full Text 
Berkowitz AE, Lahut JA, Jacobs IS, Levinson LM, Forester DW:
Phys. Rev. Lett.. 1975, 34:594.
COI number [1:CAS:528:DyaE2MXhsVyku70%3D]; Bibcode number [1975PhRvL..34..594B]
Publisher Full Text 
Europhys. Lett.. 2006, 75:618.
COI number [1:CAS:528:DC%2BD28XhtVWgs7jK]; Bibcode number [2006EL.....75..618M]
Publisher Full Text 
ElHilo M, Chantrell RW, O’Grady K:
J. Appl. Phys.. 1998, 84:5114.
COI number [1:CAS:528:DyaK1cXms1Cisr8%3D]; Bibcode number [1998JAP....84.5114E]
Publisher Full Text 
Wiekhorst F, Shevchenko E, Weller H, Kötzler J:
Phys. Rev. B. 2003, 67:224416.
Bibcode number [2003PhRvB..67v4416W]
Publisher Full Text 
J. Phys. Confer. Ser.. 2005, 17:1.
Bibcode number [2005JPhCS..17....1H]
Publisher Full Text 
Silva NJO, Amaral VS, Carlos LD:
Phys. Rev. B. 2005, 71:184408.
Bibcode number [2005PhRvB..71r4408S]
Publisher Full Text 
Phys. Rev.. 1956, 102:1217.
Bibcode number [1956PhRv..102.1217D]
Publisher Full Text 
Phys. Rev.. 1962, 127:1508.
Bibcode number [1962PhRv..127.1508K]
Publisher Full Text 
Hendriksen PV, Linderoth S, LindgÅrd PA:
Phys. Rev. B. 1993, 48:7259.
COI number [1:CAS:528:DyaK3sXmvVCmurg%3D]; Bibcode number [1993PhRvB..48.7259H]
Publisher Full Text 
Hendriksen PV, Linderoth S, LindgÅrd PA:
J. Phys. Condens. Matter. 1993, 5:5675.
COI number [1:CAS:528:DyaK3sXms1Wjtbc%3D]; Bibcode number [1993JPCM....5.5675H]
Publisher Full Text 
Plischke M, Bergersen B: Equilibrium statistical physics. World Scientific, Singapore; 1994:393.

Taylor PL: A Quantum approach to the solid state. Prentice Hall, New Jersey; 1970.

Yin LF, Wei DH, Lei N, Zhou LH, Tian CS, Dong GS, Jin XF, Guo LP, Jia QJ, Wu RQ:
Phys. Rev. Lett.. 2006, 97:067203.
COI number [1:STN:280:DC%2BD28rpsVantw%3D%3D]; Bibcode number [2006PhRvL..97f7203Y]
PubMed Abstract  Publisher Full Text 
Stevens DK, Cleland JW, Crawford JH, Schweinler HC:
Phys. Rev.. 1955, 100:1084.
COI number [1:CAS:528:DyaG28XhvVGntQ%3D%3D]; Bibcode number [1955PhRv..100.1084S]
Publisher Full Text