Method of calculating magnetic interaction of molecules using localized magnetic orbital

Abstract

A magnetic interaction is calculated by expressing a localized magnetic orbital as a linear combination of molecular orbitals and determining a localized magnetic orbital to be calculated from a maximum overlap condition evaluated between a reference orbital localized in anatomic orbital and the localized magnetic orbital. Another magnetic interaction is calculated by expressing a localized magnetic crystal orbital as a linear combination of crystal orbitals and determining a localized magnetic crystal orbital to be calculated from the maximum overlap condition.

Description

BACKGROUND OF THE INVENTION

[0001] 1) Field of the Invention

[0002] The present invention relates to a technology for calculating magnetic interaction of molecules, crystalline compounds, or high polymer compounds each exhibiting ferromagnetism.

[0003] 2) Description of the Related Art

[0004] Conventionally, an evaluation theory of magnetic characteristics has developed by a Heisenberg model which targets at a localized electron system and a Stoner model which targets at an itinerant electron system. However, since the Heisenberg model executes rough approximation and cannot be easily handled, the Heisenberg model is culled out by the Stoner model.

[0005] The Stoner model has passed through a physical examination in a Hubbard model, an Anderson model, or the like which are phenomenological models and, at present, has settled to a band calculation method based on a density functional theory. In the recent band calculation method based on the density functional theory, an effective exchange interaction energy can be obtained at calculation accuracy which is satisfactory to some extent (Sakuma Akimasa, Hitachi Metals Technical Journal, 16, 55 to 60 (2000)).

[0006] However, the conventional evaluation method (band calculation method) targets, for calculation, at an itinerant electron system which uses a Bloch function as a base, and therefore a localized electron system or an intermediate electron system between the itinerant electron system and the localized electron system has been difficult to be handled.

[0007] The band calculation method can be easily applied to an extremely small crystal system such as iron or nickel crystals and oxides thereof. However, even though the band calculation method is applied to a small crystal system, an enormous calculation quantity is required. Therefore, the following problem is left. That is, for example, it is substantially difficult that the band calculation method is applied to an actual huge crystal system such as a neodymium-based rare earth magnetic material (Nd2Fe14B crystal) from the viewpoint of calculation quantity.

[0008] On the other hand, in consideration of the circumstances, a method of easily predicting ferromagnetism by changing molecular orbitals is proposed (Japanese Patent Application Laid-Open No. H11-6825). However, magnetic characteristics detected in an experiment cannot be calculated using this method, and a subject to be calculated is limited to an organic high polymer compound. For this reason, the method poses a problem in a practical use.

SUMMARY OF THE INVENTION

[0009] It is an object of this invention to obtain a method of calculating a magnetic interaction in not only an itinerant electron system but also a localized electron system or an intermediate electron system across both the electron systems.

[0010] The magnetic interaction calculation method according to one aspect of this invention comprises calculating a magnetic interaction of molecules using a localized magnetic orbital ηa which satisfies a maximum overlap condition 1$\begin{array}{c}\sum _a^{}\left|⟨{\eta }_a|{\omega }_a⟩\right|=\mathrm{maximum},\mathrm{wherein}\\ {\eta }_a=\sum _i^{}{T}_{\mathrm{ai}}{\phi }_i\end{array}$

[0011] and ωa is a reference orbital given by 2${\omega }_a=\sum _i^{}{c}_{\mathrm{ai}}^{-1}{\phi }_i$$\mathrm{where}$${\phi }_i=\sum _a{c}_{\mathrm{ai}}{\chi }_a$

[0012] is a molecular orbital expressed as a linear combination of an atomic orbital χa, cai is an expansion coefficient, cai−1 is an inverse matrix of the expansion coefficient cai, and Tai is the expansion coefficient of the a-th orbital ηa for the i-th molecular orbital φi.

[0013] The magnetic interaction calculation method according to another aspect of this invention comprises calculating a magnetic interaction of crystalline compounds or high polymer compounds using a localized magnetic crystal orbital ηa(k) which satisfies a maximum overlap condition 3$\begin{array}{c}\sum _a\left|⟨{\eta }_a|{\omega }_a⟩\right|=\mathrm{maximum},\mathrm{wherein}\\ {\eta }_a\left(k\right)=\sum _i{T}_{\mathrm{ai}}\left(k\right){\phi }_i\left(k\right)\end{array}$

[0014] and ωa(k) is a reference orbital given by 4${\omega }_a\left(k\right)=\sum _i{c_{\mathrm{ai}}\left(k\right)}^{-1}{\phi }_i\left(k\right)$$\mathrm{where}$${\phi }_i\left(k\right)=\sum _a{c}_{\mathrm{ai}}\left(k\right)\sum _R\exp \left(i\mathrm{kR}\right){\chi }_a\left(R\right)$

[0015] is a crystal orbital expressed as a linear combination of an atomic orbital χa, cai (k) is an expansion coefficient, cai (k)−1 is an inverse matrix of the expansion coefficient cai (k), Tai(k) is the expansion coefficient of the a-th localized magnetic orbital ηa(k) for the i-th crystal orbital φi(k), k is a wave vector ranging from −π/c to π/c where c is primitive vectors, and R is a vector representing crystal periods.

[0016] The magnetic characteristic value calculation method accordingto still another aspect of this invention comprises calculating a magnetic characteristic value based on each of the magnetic interaction calculation methods.

[0017] The program according to still another aspect of this invention causes, when executed on a computer, the computer to execute each of the magnetic interaction calculation methods.

[0018] The program according to still another aspect of this invention causes, when executed on a computer, the computer to execute each of the magnetic characteristic value calculation methods. Each of the magnetic characteristic value calculation methods comprises calculating a magnetic characteristic value based on each of the magnetic interaction calculation methods.

[0019] These and other objects, features and advantages of the present invention are specifically set forth in or will become apparent from the following detailed descriptions of the invention when read in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] FIG. 1 is a diagram which shows a unit cell structure in an Nd2Fe14B crystal used in calculation when the magnetic interaction calculation method of the present invention is evaluated,

[0021] FIG. 2 is a comparative graph which shows Co-dopant quantity dependence of a Curie temperature in Nd2(Fe1−xCox)14B,

[0022] FIG. 3 is a comparative graph which shows Co-dopant quantity dependence of a residual magnetic flux density in Nd2(Fe1−xCox)14B, and

[0023] FIG. 4 is a comparative graph which shows Dy-dopant quantity dependence of residual magnetic flux density in an (Nd1−xDyx)2Fe14B crystal.

DETAILED DESCRIPTION

[0024] Embodiments of the magnetic interaction calculation method or the magnetic characteristic value calculation method according to the present invention will be described below with reference to the accompanying drawings. It should be noted that the present invention is not limited to this embodiment.

[0025] Molecular orbitals φi can be expressed as a linear combination of atomic orbitals χa by the following equation: 5${\phi }_i=\sum _a{c}_{\mathrm{ai}}{\chi }_a$

[0026] Where reference symbol cai denotes an expansion coefficient, and satisfies the following standardizing condition: 6$\begin{array}{cc}\sum _a\sum _b{c}_{\mathrm{ai}}^T{s}_{\mathrm{ab}}{c}_{\mathrm{bi}}=1& \left(1\right)\end{array}$

[0027] where Sab is an overlap integral between χa and χb.

[0028] Because the molecular orbitals φi are canonical orbitals to establish a standard orthogonal system (perfect system), localized magnetic orbitals ηa can be developed as expressed in the following equation (2) by using φi as a base: 7$\begin{array}{cc}{\eta }_a=\sum _i{T}_{\mathrm{ai}}{\phi }_i& \left(2\right)\end{array}$

[0029] It is assumed that, as a sum in the equation (2), φi is taken about occupied orbitals when the localized magnetic orbitals ηa are occupied orbitals (ηaocc), and φi is taken about unoccupied orbitals when the localized magnetic orbitals ηa is unoccupied orbitals ηavac.

[0030] A reference orbital ωa which is localized on the atomic orbitals χa is defined by using an inverse matrix cai−1 of the expansion coefficient cai as follows: 8${\omega }_a=\sum _i{c}_{\mathrm{ai}}^{-1}{\phi }_i$

[0031] In this equation, the reference orbitals ωa are an asymptotic solution. The inverse matrix cai−1 can be calculated from the equation (1).

[0032] The localized magnetic orbitals ηa can be obtained from the reference orbitals ωa and the next maximum overlap condition: 9$\begin{array}{cc}\sum _a\left|⟨{\eta }_a|{\omega }_a⟩\right|=\mathrm{maximum}& \left(3\right)\end{array}$

[0033] In addition, the equation (3) is resolved to a maximum problem of Ω(T) expressed by the following equation (4): 10$\begin{array}{cc}\begin{array}{c}\Omega \left(T\right)=\sum _a|⟨{\eta }_a|{\omega }_a⟩{|}^2\\ =\sum _a|\sum _i{T}_{\mathrm{ai}}*⟨{\phi }_i|{\omega }_a⟩{|}^2\\ =\sum _a|\sum _i\sum _j{T}_{\mathrm{ai}}*c_{\mathrm{aj}}^{-1}⟨{\phi }_i|{\phi }_j⟩{|}^2\\ =\sum _a|\sum _i{T}_{\mathrm{ai}}*c_{\mathrm{ai}}^{-1}{|}^2\end{array}& \left(4\right)\end{array}$

[0034] As is apparent from the equation (4), the localized magnetic orbitals ηa can be expressed as a subspace of the inverse matrix cai−1. Therefore, when a subspace matrix of the inverse matrix cai−1 defined by the equation (4) is represented by dia, according to the equation (2), the localized magnetic orbitals ηa can be expressed with respect to occupied orbitals ηaocc and unoccupied orbitals ηavac given by the following equations: 11$\begin{array}{c}{\eta }_a^{\mathrm{occ}}={\left(\sum _i^{\mathrm{occ}}{d}_{\mathrm{ia}}^2\right)}^{1/2}\sum _i^{\mathrm{occ}}{d}_{\mathrm{ia}}{\phi }_i\\ {\eta }_a^{\mathrm{vac}}={\left(\sum _i^{\mathrm{vac}}{d}_{\mathrm{ia}}^2\right)}^{-1/2}\sum _i^{\mathrm{vac}}{d}_{\mathrm{ia}}{\phi }_i\end{array}$

[0035] The occupied orbitals ηaocc and the unoccupied orbitals ηavac have one-to-one correspondence (isomorphism mapping) with respect to the reference orbitals ωa localized in the atomic orbitals χa.

[0036] An electron count Laocc and a hole count Lavac held by the occupied orbitals ηaocc and the unoccupied orbitals ηavac are given by the following equations, respectively: 12$\begin{array}{c}{L}_a^{\mathrm{occ}}=\sum _i^{\mathrm{occ}}{d}_{\mathrm{ia}}^2\\ L_a^{\mathrm{vac}}=\sum _i^{\mathrm{vac}}{d}_{\mathrm{ia}}^2\end{array}$

[0037] When the orbital energies of the occupied orbitals ηaocc and the unoccupied orbitals ηavac are represented by λaocc and λavac, respectively, the orbital energies λaocc and λavac are given by the following equations: 13$\begin{array}{c}{\lambda }_a^{\mathrm{occ}}={\left(\sum _i^{\mathrm{occ}}{d}_{\mathrm{ia}}^2\right)}^{-1/2}\sum _i^{\mathrm{occ}}{d}_{\mathrm{ia}}{\epsilon }_i\\ {\lambda }_a^{\mathrm{vac}}={\left(\sum _i^{\mathrm{vac}}{d}_{\mathrm{ia}}^2\right)}^{-1/2}\sum _i^{\mathrm{vac}}{d}_{\mathrm{ia}}{\epsilon }_i\end{array}$

[0038] where åi is an orbital energy of the molecular orbitals φi.

[0039] According to the above equations, an ionization potential IPa and an electron affinity EAa in the localized magnetic orbitals can be expressed from the orbital energy λaocc and the orbital energy λavac as follows: 14$\begin{array}{c}{\mathrm{IP}}_a=-{\lambda }_a^{\mathrm{occ}}\\ {\mathrm{EA}}_a=-{\lambda }_a^{\mathrm{vac}}\end{array}$

[0040] In the above description, the molecular models are explained. However, since the present invention can be similarly extended to corresponding crystal orbitals in a crystalline compound or a high polymer compound, a description thereof will be omitted.

[0041] A magnetic interaction based on localized magnetic orbitals will be described below. An effective exchange interaction Jabeff representing a magnetic interaction is approximated by the sum of a direct exchange interaction JabEX and a super exchange interaction JabSE. 15$J_{\mathrm{ab}}^{\mathrm{eff}}=J_{\mathrm{ab}}^{\mathrm{EX}}+J_{\mathrm{ab}}^{\mathrm{SE}}$

[0042] The direct (potential) exchange interaction JabEX based on the localized magnetic orbitals is expressed by the following equation: 16$J_{\mathrm{ab}}^{\mathrm{EX}}=2K_{\mathrm{ab}}{L}_a^{\mathrm{occ}}{L}_b^{\mathrm{occ}}$

[0043] In this equation, an exchange integral Kab is given by the following equation: 17$K_{\mathrm{ab}}=\int \int {\eta }_a^{*}\left(1\right){\eta }_b^{*}\left(2\right)\frac{1}{r_{12}}{\eta }_a\left(2\right){\eta }_b\left(1\right)d{\tau }_{1}d{\tau }_2$

[0044] and Laocc and Lbocc denote the numbers of electrons of the occupied orbitals ηaocc and ηbocc in the localized magnetic orbitals.

[0045] A super exchange interaction JabSE based on the localized magnetic orbitals is expressed by the following equation: 18$J_{\mathrm{ab}}^{\mathrm{SE}}=K_{\mathrm{ab}}^{\mathrm{SE}}{L}_a^{\mathrm{occ}}{L}_b^{\mathrm{vac}}$

[0046] In this equation, by using a transfer integral tab and an overlap integral sab, a kinematic integral KabSE is given by the following equation: 19$K_{\mathrm{ab}}^{\mathrm{SE}}=2t_{\mathrm{ab}}{S}_{\mathrm{ab}}$

[0047] According to the above description, a model which calculates a magnetic interaction can be achieved. Therefore, a magnetic characteristic value is then calculated by using the model, and is compared with an experimental value to be verified. Here, as the magnetic characteristic value, a Curie temperature or a residual magnetic flux density are evaluated.

[0048] In a first embodiment, in order to examine the reliability and versatility of the present invention, the physicality value of a neodymium-based rare earth magnet was evaluated. The neodymium-based rare earth magnet is an intermediate electron system having itinerant electrons and localized electrons. A unit cell in the crystal of the electron system consists of 68 atoms. This electron system is a huge system which cannot be easily calculated by a conventional band calculation method.

[0049] FIG. 1 is a diagram which shows a unit cell structure in an Nd2Fe14B crystal used in calculation. Here, the theoretical value of a Curie temperature obtained when Fe is displaced with Co is compared with a calculated value based on the present invention. The Curie temperature is a target to be evaluated because the Curie temperature is an important measure with which an effective exchange interaction energy obtained from localized magnetic crystal orbitals is compared with the experimental value. Since Fe can be completely replaced with Co, the calculation result and the experiment result can be effectively compared with each other by a replacement ratio (dopant ratio) of Co. The effect of the dopant of Co is considered by using imaginary crystal approximation which is a kind of mean field approximation.

[0050] More specifically, the Curie temperature in an Nd2(Fe1−xCox)14B crystal was evaluated. FIG. 2 is a comparative graph which shows Co-dopant quantity dependence of the Curie temperature in Nd2(Fe1−xCox)14B. The experiment result provided here is quoted from a reference “Y. Matsuura et al. Appl. Phys. Lett., 46, 308 (1985)”. As is apparent from FIG. 2, the Curie temperature with respect to the Co-dopant quantity exhibits the same trend in the experiment and the calculation.

[0051] In particular, in an area in which the Co-dopant quantity is 25% or less, the experimental value and the calculated value of the Curie temperature are almost equal to each other. Therefore, it was confirmed that the Curie temperature calculated from the present invention could be “quantitatively” compared with the experiment result. More specifically, it was confirmed that an effective exchange interaction energy based on the present invention was reliable.

[0052] Since it was understood by the present invention that the Curie temperature was quantitatively approximate to an actually measured value, simulation for Co-dopant dependence of a residual magnetic flux density was performed in a similar system. FIG. 3 is a comparative graph which shows Co-dopant dependence of a residual magnetic flux density in Nd2(Fe1−xCox)14B. The experiment result provided here was measured at 4.2 K, and is quoted from a reference “S. Hirosawa et al. J. Appl. Phys., 59,873 (1986)”. Only the experiment result related to Nd2Fe14B is described because the residual magnetic flux density of the Nd2(Fe1−xCox)14B crystal is not reported.

[0053] As shown in FIG. 3, the experimental value and the calculated value of a residual magnetic flux density in Nd2Fe14B were 18.50 kG and 16.60 kG, respectively. It was confirmed that these values were almost equal to each other. As shown in FIG. 3, the residual magnetic flux density decreases as the Co-dopant quantity increases. This downward trend is the same as that of a result related to a saturated magnetic flux density. Therefore, a calculated value (physicality value) taken along the experimental value can be extremely effectively obtained based on the method of calculating a magnetic interaction according to the present invention.

[0054] Although the calculation based on the present invention was performed by an engineering workstation, the computation time was 50 minutes. Since it was estimated that the band calculation method based on the conventional density functional theory required at least a computation time of about several weeks, it was confirmed that magnetic characteristics could be easily calculated by the calculation method according to the present invention.

[0055] In a second embodiment, as in the first embodiment, comparative evaluation of a residual magnetic flux density of a Dy2Fe14B crystal was performed. In addition, a calculation was also performed to an (Nd1−xDyx)2Fe14B crystal. FIG. 4 is a comparative graph which shows a Dy-dopant dependence of the residual magnetic flux density in the (Nd1−xDyx)2Fe14B crystal. The experiment result provided here is quoted from the reference “S. Hirosawa et al. J. Appl. Phys., 59,873 (1986)”.

[0056] As shown in FIG. 4, the experimental value and the calculated value of a residual magnetic flux density in Dy2Fe14B were 5.73 kG and 8.48 kG, respectively. When this calculated value is compared with the calculated value in Nd2Fe14B, it is understood that the calculated value can be quantitatively compared with the experimental value. The computation time required for this calculation was55 minutes. Since it was estimated that the band calculation method based on the conventional density functional theory required at least a computation time of about several weeks, it was confirmed that magnetic characteristics could be easily calculated by the calculation method according to the present invention.

[0057] The embodiments according to the present invention have been described above. The concrete example of configuration is not limited to the embodiments, and the design of the configuration may be changed without departing from the spirit and scope of the invention.

[0058] For example, a computer executable program which causes a computer to execute the magnetic interaction calculation method or the magnetic value calculation method, and this program is recorded on a computer readable recording medium. The program recorded on the recording medium is loaded on the computer to be executed, so that the respective functions may be realized.

[0059] As described above, according to the present invention, the magnetic interaction in not only the itinerant electron system but also the localized electron system or the intermediate electron system across both the electron systems can be advantageously calculated by using the localized magnetic orbitals or the localized magnetic crystal orbitals. Therefore, the magnetic characteristic value in the localized electron system or the intermediate electron system across the both electron systems can be advantageously calculated. Furthermore, the magnetic characteristic values of large molecules, a crystalline compound, and a high polymer compound can be easily calculated and predicted advantageously.

[0060] Although the invention has been described with respect to a specific embodiment for a complete and clear disclosure, the appended claims are not to be thus limited but are to be construed as embodying all modifications and alternative constructions that may occur to one skilled in the art which fairly fall within the basic teaching herein set forth.

Claims

1. A magnetic interaction calculation method comprising: calculating a magnetic interaction of molecules using a localized magnetic orbital ηa which satisfies a maximum overlap condition 20∑a ⁢ ⁢|⟨ηa|ωa⟩|=maximum,whereinηa=∑i⁢ ⁢Tai⁢φi and ωa is a reference orbital given by 21ωa=∑i⁢cai-1⁢φi⁢ ⁢where⁢ ⁢φi=cai⁢χa is a molecular orbital expressed as a linear combination of an atomic orbital χa, cai is an expansion coefficient, cai−1 is an inverse matrix of the expansion coefficient cai, and Tai is an expansion coefficient of the a-th orbital ηa for the i-th molecular orbital φi.

2. The method according to claim 1, wherein as the magnetic interaction, an orthogonal exchange interaction 2KabLaoccLbocc obtained by a potential exchange integral Kab in localized magnetic orbitals ηa and ηb and electron counts Laocc and Lbocc of occupied orbitals, is employed.

3. The method according to claim 1, wherein as the magnetic interaction, a super exchange interaction KSEabLaoccLbvac obtained by a kinematic exchange integral KSEab in localized magnetic orbitals ηa and ηb, an electron count Laocc of occupied orbitals, and a hole count Lbvac of unoccupied orbitals, is employed.

4. A magnetic interaction calculation method comprising: calculating a magnetic interaction of crystalline compounds or high polymer compounds using a localized magnetic crystal orbital ηa(k) which satisfies a maximum overlap condition 22∑a⁢&LeftBracketingBar;⟨ηa⁢&LeftBracketingBar;ωa⟩&RightBracketingBar;=maximum,wherein⁢ ⁢ ⁢ηa⁡(k)=∑i⁢Ta⁢ ⁢i⁡(k)⁢φi⁡(k) and ωa(k) is a reference orbital given by 23ωa⁡(k)=∑i⁢cai⁡(k)-1⁢φi⁡(k)⁢ ⁢whereφi⁡(k)=∑a⁢cai⁡(k)⁢∑R⁢exp⁡(ⅈ⁢ ⁢k⁢ ⁢R)⁢χa⁡(R) is a crystal orbital expressed as a linear combination of an atomic orbital χa, cai (k) is an expansion coefficient, cai (k)−1 is an inverse matrix of the expansion coefficient cai (k), Tai(k) is an expansion coefficient of the a-th localized magnetic orbital ηa(k) for the i-th crystal orbital φi(k), k is a wave vector ranging from −π/c to π/c where c is primitive vectors, and R is a vector representing crystal periods.

5. The method according to claim 4, wherein as the magnetic interaction, an orthogonal exchange interaction 2KabLaoccLbocc obtained by a potential exchange integral Kab in localized magnetic crystal orbitals ηa(k) and ηb(k) and electron counts Laocc and Lbocc of occupied orbitals, is employed.

6. The method according to claim 4, wherein as the magnetic interaction, a super exchange interaction KSEabLaoccLbvac obtained by a kinematic exchange integral KSEab in localized magnetic crystal orbitals ηa(k) and ηb(k), an electron count Laocc of occupied orbitals, and a hole count Lbvac of unoccupied orbitals, is employed.