Next Article in Journal
Physiochemical Properties and Oxidation Status of Pork from Three Rearing Systems
Previous Article in Journal
Long-Term Deflection Analysis of Large-Span Continuous Prestressed Concrete Rigid-Frame Bridges Based on a Refined Modeling Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Effect of the Back Radiation from Disk to Head after Laser Heating in HAMR

State Key Laboratory of Robotics and System, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(17), 9730; https://doi.org/10.3390/app13179730
Submission received: 9 August 2023 / Revised: 23 August 2023 / Accepted: 24 August 2023 / Published: 28 August 2023

Abstract

:
In a heat-assisted magnetic recording, the thermal effect of the head/disk interface has an important influence on the stability of the recording data. In this paper, we will discuss the thermal radiation from the disk, more specifically, the magnetic recording layer, which is at high temperature after laser heating, to the magnetic head, which is at room temperature. The radiative heat flux can be represented by the Poynting vector. In the near-field band, an effective way to obtain the electromagnetic fields is to solve the Maxwell’s equations combined with fluctuational electrodynamics. The near-field back radiation between specific head and disk material is calculated by the fluctuation-volume-current method. The radiative heat energy will induce the thermal deformation of the magnetic head, which will be discussed by the simulation, laying the foundation for adjusting and controlling the flying status.

1. Introduction

One of the effective methods to increase the storage density in a magnetic recording is to decrease the size of the magnetic particle. When the magnetic particle gets smaller, the coercivity of which gets higher and causes difficulty when writing data [1]. Heat-assisted magnetic recording (HAMR) is a new technique in which the magnetic recording media is heated by a focused laser to reduce the coercivity; hence, a high storage density can be achieved. It has great potential to increase the areal density in HAMR. The main challenge of HAMR is how to quantify the nanoscale near-field heat transfer and increase the reliability of the head/disk interface due to the back heating of the disk [2,3].
As we know, the disk, more specifically the magnetic recording layer, is heated by the focused laser to above the Curie temperature during the writing process. When the writing process is over, the disk is at a high temperature; meanwhile, the magnetic head and the air gap between the head and disk, i.e., the air bearing, are at a lower temperature. Thermal energy will transfer back from the disk to the magnetic head and the air bearing by conduction, convection, and radiation, though mainly by radiation. After absorbing the thermal energy from the hot disk, the slider can generate thermal deformation, which can change the flying height as well as the flying stability of the slider. Hence, it is very important to investigate the back heating of the HAMR disk.
However, current thermal analysis in HAMR mostly focuses on the heating process [4]. Research on the energy transfer in the head/disk interface is relatively rare. There are three layers in the head/disk interface, which are the following: (a) the lubricant layer [5]; (b) the near-field transducer [6]; and (c) the air bearing [7,8,9]. As discussed before, the thermal effects on the magnetic head have a great influence on the air bearing indirectly. Research on the thermal effects of the slider on the head is mainly focused on the results considering multiple factors rather than a single factor [10]. The net temperature rise of the slider caused by the laser absorption and the heat radiation of the locally heated recording magnetic disk are studied in Ref. [9], where the influence of the laser absorption and the heat radiation is considered combined. A new model, together with a finite element (FE) model for an integrated HAMR slider and an optical absorption model for a temperature profile of media hot spot, is then used to solve for the temperature rise and thermal protrusion on the slider body in Ref. [10]. However, because of the limitations of the experimental technology, the accuracy of the theoretical calculation is difficult to verify. As the development of the high-precision micro-electromechanical system (MEMS), tip-enhanced Raman spectroscopy (TERS) temperature measurement and other experimental technology progressed in recent years, the theoretical calculation of the near-field radiative heat transfer is gradually revised and improved so that the particular research on the influence of the near-field radiative heat transfer on the slider is badly needed to be supplemented.
The flying height, i.e., the distance between the head and disk, has been reduced to be around 2 nm, which is much smaller than the characteristic wavelength of the electromagnetic waves. Due to the scale effect, such as the surface wave, evanescent wave, radiation tunneling, and so on, traditional radiation theories, such as Kirchhoff’s law and Planck’s law, cannot be used to predict the radiative heat transfer. Research on the near-field heat transfer is concentrated on the theory of heat transport by phonon tunneling, see Refs. [11,12,13]. In recent years [14], the theoretical and experimental studies on the near-field radiative heat transfer have obviously increased. For different theoretical models, the calculation methods of the near-field radiative heat transfer are different for different geometries, such as heat transfer between two close infinite bulks [15], a small spherical particle and a bulk [16,17], two close infinite bulks with coats on the surface [18,19], several arbitrary shape objects [20,21], nanoparticles [22,23], thin films [24,25], and so on. The main difference between the calculation methods is focused on the method to solve the electric field and magnetic field on the target position, i.e., the solving methods of the stochastic the Maxwell’s equations are different. The most common method to solve the Maxwell’s equations is the discrete approximation method [26], such as the finite-different time-domain (FDTD) method, finite-difference frequency-domain (FDFD) method, boundary finite element (BEM) method, and the thermal dipole discrete approximation (T-DDA) method [20,21]. The limitation of the existing discrete approximation method is that the calculation is of great complexity and time-consuming. Another mainstream method to solve the Maxwell’s equations is to bring in the fluctuation–dissipation theory [27,28], the main idea of which is to regard an object as a collection of many dipoles in constant motion and then to integrate over all sources in the volume of the object to obtain the electromagnetic field on the target position. For simple shape objects, the volume of the object is generally integrated directly, known as the fluctuate-volume-current (FVC) method [25,28,29]. And for complex shape objects, the volume integrals can be converted into surface integrals on the surface of objects to calculate the integral, known as the fluctuate-surface-current (FSC) method [28]. Back to the head/disk model discussed in this paper, the fluctuate-volume-current method can be chosen to calculate the near-field radiative heat transfer.
The paper is constructed as follows: The theoretical analysis of the fluctuate-volume-current method to calculate the near-field radiative heat flux by the Poynting’s theory from disk to head is presented in Section 2. The properties of the imitating material are considered, and then the numerical results and discussion of the near-field radiative heat flux from disk to head in HAMR are given in Section 3. At last, the simulation results and discussion of the thermal deformation of the magnetic head are given in Section 4.

2. Numerical Analysis

2.1. Theoretical Model

The HAMR disk is consisted of the following several layers: lubricant film (~1 nm), carbon overcoat (~2 nm), magnetic recording layer (~15 nm), interlayer (~30 nm), heat sink (~15 μm), and substrate. A schematic diagram of an HAMR recording system is shown in Figure 1. The thickness of the magnetic recording layer is around 15 nm, much smaller than the horizontal size of the disk so that it can be simplified as a nanoscale film while the head/disk interface system is simplified as a multilayered structure, shown in Figure 2. Note that to analyze the thermal effect of radiation individually, in order to make targeted corrections to the effects of radiation in the future, the conduction and convection heat transfer is ignored here. Meanwhile, due to the neglect of the conduction heat transfer, assume that the lower surface of the magnetic recording layer is thermal insulation material. In Figure 2, there are 4 layers, numbered from 0 to 3 from bottom to top. The half space 0 under the magnetic recording film is the dielectric substrate, which is in thermal equilibrium, at a temperature of T 0 , the permittivity is ε 0 . The film layer 1 represents for the magnetic recording film, which is assumed to be in thermal equilibrium and at a temperature of T 1 . The magnetic recording layer is made of magnetic material, the permittivity and permeability of which is represented as ε 1 and μ 1 , respectively. The thickness of film layer 1 is assumed to be d 1 . The film 2 between the head and disk is vacuum, in a thickness of d 2 . Similarly, the half space 3 represents for the head, which is in thermal equilibrium, and the temperature, which is T 3 , the permittivity and permeability is represented as ε 3 and μ 3 , respectively.

2.2. Spectral Heat Flux

The numerical problem is to discuss the radiative heat flux flow into half space 3 from film 1, as shown in Figure 2. In fluctuational electrodynamics, the source of thermal radiation is modeled as stochastic electric and magnetic current densities j e and j m , respectively. The electric and magnetic fields at point r due to electric source j e located at r inside a body of volume V are given by the following:
E e ( r , ω ) = i μ v μ s ω V G ¯ ¯ E e ( r , r , ω ) · j e ( r , ω ) d V
H e ( r , ω ) = V G ¯ ¯ H e ( r , r , ω ) · j e ( r , ω ) d V
Similarly, the electric and magnetic fields at point r due to magnetic source j m located at r inside a body of volume V are given by the following:
E m ( r , ω ) = V G ¯ ¯ E m ( r , r , ω ) · j m ( r , ω ) d V
H m ( r , ω ) = i ε v ε s ω V G ¯ ¯ H m ( r , r , ω ) · j m ( r , ω ) d V
where μ v and ε v represent for the permittivity and permeability in vacuum, generally μ v = 4 π × 10 7 NA 2 and ε v = 8.854187817 × 10 12   F / m , respectively. Similarly, μ s and ε s represent for the permittivity and permeability of media s , respectively. Terms G ¯ ¯ E e ( r , r , ω ) , G ¯ ¯ E m ( r , r , ω ) , G ¯ ¯ H e ( r , r , ω ) , and G ¯ ¯ H m ( r , r , ω ) are the electric and magnetic DGFs relating the fields with frequency ω at point r due to a source located at r , respectively. The first E / H in the superscript of DGFs represents for the type of the field at point r , and the second e / m in the superscript of DGF represents for the type of the source at point r . The total electric field at point r due to sources located at r inside a body of volume V is the sum of two types of sources, i.e., the sum of Equations (1) and (3), namely, the total magnetic field.
Generally, the heat flux of the electromagnetic fields can be calculated by the Poynting vector, which is the rate of energy flow per unit area in a plane electromagnetic wave. The electromagnetic fields radiated from a thermal object changes over time and can be represented by the Poynting vector. So that the instantaneous energy flow density, i.e., the instantaneous heat flux, is given by the following:
S = E × H
The monochromatic radiative heat flux is given by the time-averaged Poynting vector as follows:
S ( r , ω ) = 1 2 Re [ E ( r , ω ) × H ( r , ω ) ]
After ignoring the horizontal boundary, there is no energy transfer along the x and y directions in Cartesian coordinates, so that only the z-direction heat flux is considered. Since only positive frequencies are considered in the Fourier decomposition of the time-dependent fields into frequency-dependent quantities, the spectral heat flux is given by the following:
q ( z , ω ) = S z ( r , ω ) = 4 × 1 2 Re [ E x H y E y H x ]
The ensemble average of the spatial correlation functions of the stochastic electric and magnetic current densities are given by the following [18,25]:
j α e ( r , ω ) j β e * ( r , ω ) = ω ε v Im ( ε s ) Θ ( ω , T ) π δ ( r r ) δ ( ω ω ) δ α β
j α m ( r , ω ) j β m * ( r , ω ) = ω μ v Im ( μ s ) Θ ( ω , T ) π δ ( r r ) δ ( ω ω ) δ α β
j α e ( r , ω ) j β m * ( r , ω ) = 0
j α m ( r , ω ) j β e * ( r , ω ) = 0
where the mean energy of a Planck’s oscillator in thermal equilibrium at frequency ω and temperature T is [27,28,30]
Θ ( ω , T ) = ω exp ( ω / k B T ) 1 + 1 2 ω
where ω / 2 represents for the term of vacuum fluctuation, i.e., the zero-point energy, the zero-point energy of two objects can be cancelled out; therefore, it can be ignored during the calculation.
Substitution of the electric and magnetic fields Equations (1)–(4) and the ensemble average of the spatial correlation functions Equations (8)–(11) into heat flux Equation (7), Equation (7) can be reduced into the following:
q ( z , ω ) = 2 k v 2 Θ ( ω , T ) π Re [ i Im ( ε s ) μ s V d V ( G ¯ ¯ x α E e G ¯ ¯ y α H e G ¯ ¯ y α E e G ¯ ¯ x α H e ) i ε s Im ( μ s ) V d V ( G ¯ ¯ x α E m G ¯ ¯ y α H m G ¯ ¯ y α E m G ¯ ¯ x α H m ) ]
where the subscript α refers to all the orthogonal components, i.e., α = x , y , z in Cartesian coordinates. k v represents for the wavevector in vacuum, and there is a relationship k v 2 = ω 2 μ v ε v .
To solve Equation (13), the dyadic correlation term V d V G ¯ ¯ u α E e G ¯ ¯ v α H e and V d V G ¯ ¯ u α E m G ¯ ¯ v α H m need to be determined first. Here, the subscripts u , v = x , y , u v and α = x , y , z have been discussed before. Because of the transition of the Cartesian coordinate to the cylindrical coordinate, the subscripts hereafter stand for u , v = ρ , θ , u v and α = ρ , θ , z . Therefore, the dyadic correlation term in DGFs of fields generated by the electric sources leads to the following:
V d V G ¯ ¯ u α E e G ¯ ¯ v α H e = - d k ρ ( 2 π ) 2 z d z g u α E e g v α H e = 1 2 π 0 k ρ d k ρ z d z g u α E e g v α H e
Similarly, the dyadic correlation term in DGFs of fields generated by the magnetic sources leads to the following:
V d V G ¯ ¯ u α E m G ¯ ¯ v α H m = 1 2 π 0 k ρ d k ρ z d z g u α E m g v α H m
The source object in Figure 2 is the film 1, so that take s = 1 into Equation (13). And then substitution of Equations (14) and (15) into Equation (13), the monochromatic radiative heat flux under cylindrical coordinate is given by the following:
q ( z , ω ) = k v 2 Θ ( ω , T ) π 2 Re [ i Im ( ε 1 ) μ 1 0 k ρ d k ρ z d z ( g ρ ρ E e g θ ρ H e g θ θ E e g ρ θ H e + g ρ z E e g θ z H e ) i ε 1 Im ( μ 1 ) 0 k ρ d k ρ z d z ( g θ ρ E m g ρ ρ H m + g ρ θ E m g θ θ H m g θ z E m g ρ z H m ) ]

2.3. Dyadic Green’s Function

The sources in the film generate two kinds of electromagnetic waves, emitting to the forward and backward directions. Meanwhile, resulting from multiple reflections between the upper and lower surfaces of the film, the electromagnetic waves travel both forward and backward in the end. So that the electromagnetic waves consist of four parts, i.e., the forward travelling waves arising from a source emitting in the forward direction, the backward travelling waves arising from a source emitting in the forward direction, the forward travelling waves arising from a source emitting in the backward direction, and the backward travelling waves arising from a source emitting in the backward direction, as shown in Figure 3.
DGF between source r located in the magnetic recording layer film 1 and the target position r located in magnetic head half space 3 is given by the following:
G ¯ ¯ E e ( r , r , ω ) = i 8 π 2 d k ρ k z 1 e i k ρ ( R R ) [ ( t 13 s + + s ^ s ^ + t 13 p + + p ^ 3 + p ^ 1 + ) e i k z 3 ( z z 3 ) i k z 1 z + ( t 13 s + s ^ s ^ + t 13 p + p ^ 3 p ^ 1 + ) e i k z 3 ( z z 3 ) i k z 1 z + ( t 13 s + s ^ s ^ + t 13 p + p ^ 3 + p ^ 1 ) e i k z 3 ( z z 3 ) + i k z 1 z + ( t 13 s s ^ s ^ + t 13 p p ^ 3 p ^ 1 ) e i k z 3 ( z z 3 ) + i k z 1 z ]
where the term t s l p + + p ^ l + p ^ s + stands for the direction changes of the polarized waves, for example, the term t 13 p + + p ^ 3 + p ^ 1 + means the direction of p-polarized waves in film 1 is p ^ 1 + , and it changes to p ^ 3 + in bulk 3, meanwhile the amplitude has a coefficient t 13 p + + , the superscript refers to the type of the electromagnetic, i.e., p + + stands for the forward travelling waves arising from a source emitting in the forward direction. To make it clear, the amplitude coefficient t 12 s , p + + , t 12 s , p + , t 12 s , p + , t 12 s , p can be represented as A s , p , B s , p , C s , p , D s , p . And substitution of the unit vectors of s-polarized and p-polarized into Equation (17), DGFs of the electric field due to electric sources are given by the following:
G ¯ ¯ E e ( r , r , ω ) = i 8 π 2 d k ρ k z 1 e i k ρ ( R R ) × [ k z 1 k z 3 k 1 k 3 ( A p e i k z 3 ( z z 3 ) i k z 1 z B p e i k z 3 ( z z 3 ) i k z 1 z C p e i k z 3 ( z z 3 ) + i k z 1 z + D p e i k z 3 ( z z 3 ) + i k z 1 z ) ρ ^ ρ ^ + k ρ k z 3 k 1 k 3 ( A p e i k z 3 ( z z 3 ) i k z 1 z + B p e i k z 3 ( z z 3 ) i k z 1 z C p e i k z 3 ( z z 3 ) + i k z 1 z + D p e i k z 3 ( z z 3 ) + i k z 1 z ) ρ ^ z ^ + ( A s e i k z 3 ( z z 3 ) i k z 1 z + B s e i k z 3 ( z z 3 ) i k z 1 z + C s e i k z 3 ( z z 3 ) + i k z 1 z + D s e i k z 3 ( z z 3 ) + i k z 1 z ) θ ^ θ ^ + k ρ k z 1 k 1 k 3 ( A p e i k z 3 ( z z 3 ) i k z 1 z B p e i k z 3 ( z z 3 ) i k z 1 z + C p e i k z 3 ( z z 3 ) + i k z 1 z + D p e i k z 3 ( z z 3 ) + i k z 1 z ) z ^ ρ ^ + k ρ 2 k 1 k 3 ( A p e i k z 3 ( z z 3 ) i k z 1 z + B p e i k z 3 ( z z 3 ) i k z 1 z + C p e i k z 3 ( z z 3 ) + i k z 1 z + D p e i k z 3 ( z z 3 ) + i k z 1 z ) z ^ z ^ ]
It is easy to deduce the Wyle components respectively by Equation (18). And the Wyle components of other DGFs can be obtained similarly. Substituting the DGFs into the heat flux Equation (16) and noting that the integral interval of the magnetic recording layer is z = [ z 1 , z 2 ] = [ 0 , d ] . Finally, the spectral heat flux is given by the following:
q ( z , ω ) = k v 2 Θ ( ω , T ) π 2 Re [ i Im ( ε 1 ) μ 1 0 k ρ d k ρ 0 d d z ( g ρ ρ E e g θ ρ H e g θ θ E e g ρ θ H e + g ρ z E e g θ z H e ) i ε 1 Im ( μ 1 ) 0 k ρ d k ρ 0 d d z ( g θ ρ E m g ρ ρ H m + g ρ θ E m g θ θ H m g θ z E m g ρ z H m ) ]

3. Results and Discussion

3.1. Materials in HAMR

The magnetic recording layer in HAMR is made out of magnetic materials, which is in progress along with the development of the magnetic recording technology. Until now, FePt has become the next generation of high-density magnetic recording media material. FePt has significant advantages in magnetic recording because it has several excellent characteristics, such as high anisotropy, similar saturation magnetization, small superparamagnetic critical size, and the ability to keep high thermal stability when the size of the particle is reduced to 3–4 nm [31,32,33,34]. To meet the practical requirement of a magnetic recording medium, the alloy needs to be modified by doping. But, in theoretical calculations, the modifying of the material can be ignored [32,35]. Several studies on the properties of alloys show that the dielectric function of thin alloy films of different alloy ratios basically fluctuates in a similar trend and within a small range [36,37]. To simplify the calculation, the parameters of a single metal material can be used instead of a complex alloy material.
ε = 1 ω p 2 ω 2 + i γ e ω
where ω p and γ e stand for the equivalent plasma frequency and damping factor of the material, respectively. The two parameters of metal Fe are ω p = 0.62204 × 10 16   rad / s and γ e = 0.27709 × 2 π × 10 14   rad / s . In the radiation band, the susceptibility of the metal is very small so that the permeability can be regarded as constant by ignoring the influence on the radiative heat flux contributed by the magnetism.
The permittivity of the material of the magnetic head for calculation, Al 2 O 3 , is approximated by the following:
ε = ε j Ω j L O 2 ω 2 + i γ j L O ω Ω j T O 2 ω 2 + i γ j T O ω
where Ω j L O and Ω j T O are longitudinal and transverse resonance frequency respectively, γ j L O and γ j T O are longitudinal and transverse damping factor respectively. These parameters are related to the temperature, as discussed in Ref. [38].

3.2. Spectral Heat Flux of the Back Radiation in HAMR

In Figure 2, assume that the thickness of film 1 is 15 nm, the flying height, i.e., the thickness of air film 2 between head and disk is 2 nm, the radiative heat flux from film 1 to bulk 3 is calculated at full frequency after ignoring the frequencies where the heat flux is several orders of magnitude below the maximum value. The frequency–flux graph at an effective frequency band, from ω = 1 × 10 11   rad / s to ω = 1 × 10 15   rad / s , is plotted shown in Figure 4. Overall, the contribution of the near-field radiative heat flux (the red line) is several orders of magnitude larger than that of the far-field radiative heat flux (the blue line), related to the frequency band; thus, the contribution of far-field for radiative heat flux can be ignored here. Moreover, the near-field heat transfer is enhanced at the resonant frequency, which is related to the properties of the material. Notice the near-field part, there is a remarkable fluctuation of the heat flux during the frequencies from ω 1 = 1.10 × 10 13   rad / s to ω 2 = 2.04 × 10 13   rad / s , which is because of the electromagnetic properties of the materials. The abrupt change in the heat flux at frequencies ω 1 = 1.10 × 10 13   rad / s (Figure 5a) and ω 2 = 2.04 × 10 13   rad / s (Figure 5b) is because the sign of the imaginary part of the electrical permittivity changes at these two frequencies, i.e., when ω < ω 1 = 1.10 × 10 13   rad / s , the imaginary part of the electrical permittivity is negative, and it changes to positive when ω 1 = 1.10 × 10 13   rad / s < ω < ω 2 = 2.04 × 10 13   rad / s , and changes again to negative when ω > ω 2 = 2.04 × 10 13 rad / s , as shown in Figure 5. Noticed that there is a deviation between the peaks of the heat flux and the peaks of the electrical permittivity of Al 2 O 3 during the frequency band from ω 1 = 1.10 × 10 13   rad / s to ω 2 = 2.04 × 10 13   rad / s , which is because of the influence of the imaginary part of the electrical permittivity of Fe, as shown in Figure 6. Generally, the electromagnetic properties of the materials, particularly the electrical permittivity of the materials, have great influence on the spectral heat flux. In addition, the surface phonon–polariton resonance is visible during the frequencies ω 2 = 2.04 × 10 13   rad / s to ω 4 = 4.40 × 10 13   rad / s , and the large value of the near-field heat flux is at the frequency ω 3 = 3.42 × 10 13   rad / s .
The spectral heat flux at an effective frequency band under varying flying heights from 1 nm to 10 nm is plotted in Figure 7 to discuss the influence of the flying height on the spectral heat flux. As shown in Figure 7, the surface phonon–polariton resonance frequency band of the heat flux in the extreme near-field distance, i.e., d = 1 nm and 2 nm, is narrower than that in the far-field band, i.e., d = 5 nm and 10 nm. In addition, there is a remarkable end of the resonance frequency for the extreme near-field radiative heat transfer.

3.3. The Influence of the Flying Height on the Back Radiation in HAMR

As shown in Figure 7, obviously, the radiative heat transfer is decreased when the flying height is increased, and the peaks appear in the same frequencies, but the peaks become smoother. By integrating over the whole frequencies, the total near-field back radiation heat flux can be calculated. Taking the varying flying height from 1 nm to 10 nm, the graph of heat flux varying with flying height is plotted in Figure 8. Obviously, when the flying height is small enough, the contribution of the far-field heat flux can be ignored, and the contribution of the near-field heat flux plays a major role in the total flux. As shown in Figure 8, the near-field heat flux (the red line) is almost overlapped with the total heat flux (the black line). Even under the extremely small flying height, when the flying height is d = 1   nm , the heat flux is almost two orders of magnitude larger than that when the flying height is d = 10   nm . So, because the near-field radiative heat transfer is very sensitive to the changes in the flying height, slight fluctuations in the flying height can lead to huge near-field radiative heat transfer changes, resulting in thermal deformation of the head, the failure of the near-field transducer, the instability of the flying status, and so on.

4. Simulation Results and Discussion of the Thermal Deformation of Head

The back radiative heat flux from disk to head, which can be regarded as the thermal source on the disk heating the slider on the head, will cause a thermal deformation of the slider. The slider is modelled to find out the maximal deformation of the head. It makes no sense to consider the whole size of the magnetic head because, in comparison with the heated area, the magnetic head is too big; therefore, a calculation diameter of 4   μ m for the cylinder head is taken. The heated area is taken to a 200   nm × 200   nm square in the center of the head. The thickness of the cylinder head is much larger than the size of the heated area too, so a thickness of 10 times the size of the heated area, 2   μ m , is taken. After heating the square area by the radiative heat flux calculated in Section 2, the flying height of 1 nm was picked, under which distance the radiative heat flux is large enough to obviously observe the deformation, the distribution of temperature is plotted in Figure 9. The high temperature appears inside the heating square. The temperature is higher inside and lower outside the square. The heating has almost no influence on the other area in the head model.
The deformation caused by the radiative heat flux is shown in Figure 10. As expected, the largest deformation is in the center of the head model (Figure 10b), and the deformation in the square area is larger than the other area in the head model (Figure 10a). It shows that when the flying height is 1 nm, the maximum deformation is about 0.058 nm. It is going to make a big difference on the air bearing, i.e., the 1 nm-thick air film.
The thermal deformation and deformation rate, i.e., the rate of the deformation under different flying heights from 1 nm to 10 nm, is plotted in Figure 11. As shown in Figure 11, the deformation rate is increased when the flying height gets smaller, especially when the flying height is smaller than 5 nm, the deformation rate is too big to be ignored. It is because when the flying height is increased, the radiative heat flux is decreased, as discussed in Section 2. Ref. [10] shows that radiation at very small gap spacing could have a more significant effect on the thermal protrusions (thermal deformation) by comparing the thermal deformation with and without back radiation considered. The simulation results in Figure 11 basically coincide with the results in Ref. [10]. Now that the flying height reaching 2 nm, the discussion of the thermal deformation is meaningful.

5. Conclusions

In this paper, the fluctuate-volume-current (FVC) method is used to calculate the near-field back radiative heat transfer from the disk to the magnetic head in heat-assisted magnetic recording (HAMR). And the thermal deformation of the slider on the magnetic head caused by the back radiative heat from the hot magnetic recording layer is then analyzed.
The formulation for calculating the near-field radiative heat transfer by the fluctuate-surface-current method is based on the fluctuation–dissipation theory, which was brought in to solve Maxwell’s equations to obtain the electromagnetic fields on the target position in the magnetic head, and then the heat flux using the Poynting’s theory is obtained.
Materials Fe and Al 2 O 3 are chosen to represent the magnetic recording film and the slider on the magnetic head, respectively. The radiative heat flux from the hot magnetic recording layer to the cold slider is calculated. It is shown in the spectral heat flux graph that the fluctuations of the heat flux are caused by the electromagnetic properties of the materials, and the frequencies reached the surface phonon–polariton resonance.
It is shown in the graph of heat flux varying with flying height that, under the small flying height band, the contribution of the far-field heat flux can be ignored. And even under the extremely small flying height from 1 nm to 10 nm, the heat flux has almost two orders of magnitude changes.
The thermal analysis and the structural analysis of the slider on the magnetic head are used to discuss the influence of the back radiation on the structure of the slider. It is shown that when the distance between the head and disk gets larger, the deformation rate decreases and becomes small enough to be ignored when the distance is larger than 5 nm. Now that the flying height is reaching 2 nm, the thermal deformation may cause an obvious decrease in the air bearing (i.e., the gap between slider and disk), leading to the fluctuation of the head flying. More research on eliminating the negative effects of the back radiation thermal effects should be conducted in the future, such as how to reduce the near-field radiative heat transfer, how to achieve ultra-fast heat dissipation, and how to maintain the dynamic balance of the flying height while the slider is deformed by the back heating from the disk. The analysis of the thermal effect of the back radiation in this paper will be a preparation for these studies.

6. Patents

This section is not mandatory but may be added if there are patents resulting from the work reported in this manuscript.

Author Contributions

Methodology, investigation, and writing, Y.Z.; supervision and funding acquisition, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kryder, M.H.; Gage, E.C.; McDaniel, T.W.; Challener, W.A.; Rottmayer, R.E.; Ganping, J.; Yiao-Tee, H.; Erden, M.F. Heat Assisted Magnetic Recording. Proc. IEEE 2008, 96, 1810–1835. [Google Scholar] [CrossRef]
  2. Rottmayer, R.E.; Batra, S.; Buechel, D.; Challener, W.A.; Hohlfeld, J.; Kubota, Y.; Li, L.; Lu, B.; Mihalcea, C.; Mountfield, K.; et al. Heat-Assisted Magnetic Recording. IEEE Trans. Magn. 2006, 42, 2417–2421. [Google Scholar] [CrossRef]
  3. Pan, L.; Bogy, D.B. Data storage: Heat-assisted magnetic recording. Nat. Photonics 2009, 3, 189–190. [Google Scholar] [CrossRef]
  4. Hu, Y.; Wu, H.; Meng, Y.; Bogy, D.B. Nanoscale thermal analysis for heat-assisted magnetic recording. J. Appl. Phys. 2017, 122, 134303. [Google Scholar] [CrossRef]
  5. Zeng, Y.; Zhou, W.; Huang, X.; Yu, S. Numerical study on thermal-induced lubricant depletion in laser heat-assisted magnetic recording systems. Int. J. Heat Mass Transf. 2012, 55, 886–896. [Google Scholar] [CrossRef]
  6. Trinh, T.D.; Rajauria, S.; Smith, R.; Schreck, E.; Dai, Q.; Talke, F.E. Temperature-Induced Near-Field Transducer Failure in Heat-Assisted Magnetic Recording. IEEE Trans. Magn. 2020, 56, 3300304. [Google Scholar] [CrossRef]
  7. Zheng, H.; Li, H.; Talke, F.E. Numerical simulation of thermal flying height control sliders in heat-assisted magnetic recording. Microsyst. Technol. 2012, 18, 1731–1739. [Google Scholar] [CrossRef]
  8. Xu, B.X.; Yuan, H.X.; Zhang, J.; Yang, J.P.; Ji, R.; Chong, T.C. Thermal effect on slider flight height in heat assisted magnetic recording. J. Appl. Phys. 2008, 103, 07F525. [Google Scholar] [CrossRef]
  9. Xu, B.X.; Hu, S.B.; Yuan, H.X.; Zhang, J.; Chen, Y.J.; Ji, R.; Miao, X.S.; Chen, J.S.; Chong, T.C. Thermal effects of heated magnetic disk on the slider in heat-assisted magnetic recording. J. Appl. Phys. 2006, 99, 08N102. [Google Scholar] [CrossRef]
  10. Myo, K.S.; Zhou, W.; Huang, X.; Yu, P. Back heating effect of media hot spot at nanoscale head–disk interface. Microsyst. Technol. 2016, 23, 2587–2597. [Google Scholar] [CrossRef]
  11. Budaev, B.V.; Bogy, D.B. Mechanisms of heat transport across a nano-scale gap in heat assisted magnetic recording. J. Appl. Phys. 2012, 111, 124508. [Google Scholar] [CrossRef]
  12. Budaev, B.V.; Bogy, D.B. Heat transport by phonon tunneling across layered structures used in heat assisted magnetic recording. J. Appl. Phys. 2015, 117, 104512. [Google Scholar] [CrossRef]
  13. Wu, H.; Bogy, D. Use of an embedded contact sensor to study nanoscale heat transfer in heat assisted magnetic recording. Appl. Phys. Lett. 2017, 110, 033104. [Google Scholar] [CrossRef]
  14. Basu, S.; Zhang, Z.M.; Fu, C.J. Review of near-field thermal radiation and its application to energy conversion. Int. J. Energy Res. 2009, 33, 1203–1232. [Google Scholar] [CrossRef]
  15. Hu, L.; Narayanaswamy, A.; Chen, X.; Chen, G. Near-field thermal radiation between two closely spaced glass plates exceeding Planck’s blackbody radiation law. Appl. Phys. Lett. 2008, 92, 133106. [Google Scholar] [CrossRef]
  16. Mulet, J.-P.; Joulain, K.; Carminati, R.; Greffet, J.-J. Nanoscale radiative heat transfer between a small particle and a plane surface. Appl. Phys. Lett. 2001, 78, 2931–2933. [Google Scholar] [CrossRef]
  17. Shen, S.; Narayanaswamy, A.; Chen, G. Surface phonon polaritons mediated energy transfer between nanoscale gaps. Nano Lett. 2009, 9, 2909–2913. [Google Scholar] [CrossRef]
  18. Biehs, S.A. Thermal heat radiation, near-field energy density and near-field radiative heat transfer of coated materials. Eur. Phys. J. B 2007, 58, 423–431. [Google Scholar] [CrossRef]
  19. Fu, C.J.; Tan, W.C. Near-field radiative heat transfer between two plane surfaces with one having a dielectric coating. J. Quant. Spectrosc. Radiat. Transf. 2009, 110, 1027–1036. [Google Scholar] [CrossRef]
  20. Edalatpour, S.; Francoeur, M. Near-field radiative heat transfer between arbitrarily shaped objects and a surface. Phys. Rev. B 2016, 94, 045406. [Google Scholar] [CrossRef]
  21. Luo, M.; Dong, J.; Zhao, J.; Liu, L.; Antezza, M. Radiative heat transfer between metallic nanoparticle clusters in both near field and far field. Phys. Rev. B 2019, 99, 134207. [Google Scholar] [CrossRef]
  22. Chapuis, P.-O.; Laroche, M.; Volz, S.; Greffet, J.-J. Near-field induction heating of metallic nanoparticles due to infrared magnetic dipole contribution. Phys. Rev. B 2008, 77, 125402. [Google Scholar] [CrossRef]
  23. Chapuis, P.-O.; Laroche, M.; Volz, S.; Greffet, J.-J. Radiative heat transfer between metallic nanoparticles. Appl. Phys. Lett. 2008, 92, 201906. [Google Scholar] [CrossRef]
  24. Francoeur, M.; Mengüç, M.P.; Vaillon, R. Near-field radiative heat transfer enhancement via surface phonon polaritons coupling in thin films. Appl. Phys. Lett. 2008, 93, 043109. [Google Scholar] [CrossRef]
  25. Francoeur, M.; Pinar Mengüç, M.; Vaillon, R. Solution of near-field thermal radiation in one-dimensional layered media using dyadic Green’s functions and the scattering matrix method. J. Quant. Spectrosc. Radiat. Transf. 2009, 110, 2002–2018. [Google Scholar] [CrossRef]
  26. Yurkin, M.A.; Hoekstra, A.G. The discrete dipole approximation: An overview and recent developments. J. Quant. Spectrosc. Radiat. Transf. 2007, 106, 558–589. [Google Scholar] [CrossRef]
  27. Rytov, S.M.; Kravtsov, Y.A.; Tatarskii, V.I.; Today, K.J.P. Principles of Statistical Radiophysics; Springer: Berlin/Heidelberg, Germany, 1989; Volume 42, p. 65. [Google Scholar]
  28. Joulain, K.; Mulet, J.-P.; Marquier, F.; Carminati, R.; Greffet, J.-J. Surface electromagnetic waves thermally excited: Radiative heat transfer, coherence properties and Casimir forces revisited in the near field. Surf. Sci. Rep. 2005, 57, 59–112. [Google Scholar] [CrossRef]
  29. Zheng, Z.; Xuan, Y. Theory of near-field radiative heat transfer for stratified magnetic media. Int. J. Heat Mass Transf. 2011, 54, 1101–1110. [Google Scholar] [CrossRef]
  30. Zheng, Z.; Xuan, Y. Near-field radiative heat transfer between general materials and metamaterials. Chin. Sci. Bull. 2011, 56, 2312–2319. [Google Scholar] [CrossRef]
  31. Fernandez, R.; Teweldebrhan, D.; Zhang, C.; Balandin, A.; Khizroev, S. A comparative analysis of Ag and Cu heat sink layers in L10-FePt films for heat-assisted magnetic recording. J. Appl. Phys. 2011, 109, 07B763. [Google Scholar] [CrossRef]
  32. Cen, Z.H.; Xu, B.X.; Hu, J.F.; Li, J.M.; Cher, K.M.; Toh, Y.T.; Ye, K.D.; Zhang, J. Optical property study of FePt-C nanocomposite thin film for heat-assisted magnetic recording. Opt. Express 2013, 21, 9906–9914. [Google Scholar] [CrossRef] [PubMed]
  33. Chernyshov, A.; Treves, D.; Le, T.; Zong, F.; Ajan, A.; Acharya, R. Measurement of FePt thermal properties relevant to heat-assisted magnetic recording. J. Appl. Phys. 2014, 115, 17B735. [Google Scholar] [CrossRef]
  34. Willig, L.; Reppert, A.v.; Deb, M.; Ganss, F.; Hellwig, O.; Bargheer, M. Finite-size effects in ultrafast remagnetization dynamics of FePt. Phys. Rev. B 2019, 100, 224408. [Google Scholar] [CrossRef]
  35. Chen, S.; Shu, X.; Xie, Q.; Zhou, C.; Zhou, J.; Deng, J.; Guo, R.; Peng, Y.G.; Ju, G.; Chen, J. Structure, magnetic and thermal properties of FePt–C–BN granular films for heat assisted magnetic recording. J. Phys. D Appl. Phys. 2020, 53, 135002. [Google Scholar] [CrossRef]
  36. Rioux, D.; Vallières, S.; Besner, S.; Muñoz, P.; Mazur, E.; Meunier, M. An Analytic Model for the Dielectric Function of Au, Ag, and their Alloys. Adv. Opt. Mater. 2014, 2, 176–182. [Google Scholar] [CrossRef]
  37. Tikuišis, K.K.; Beran, L.; Cejpek, P.; Uhlířová, K.; Hamrle, J.; Vaňatka, M.; Urbánek, M.; Veis, M. Optical and magneto-optical properties of permalloy thin films in 0.7–6.4 eV photon energy range. Mater. Des. 2017, 114, 31–39. [Google Scholar] [CrossRef]
  38. Gervais, F.; Piriou, B. Anharmonicity in several-polar-mode crystals: Adjusting phonon self-energy of LO and TO modes in Al2O3 and TiO2 to fit infrared reflectivity. J. Phys. C Solid State Phys. 1974, 7, 2374–2386. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of an HAMR recording system.
Figure 1. Schematic diagram of an HAMR recording system.
Applsci 13 09730 g001
Figure 2. Simplified model of head/disk structure.
Figure 2. Simplified model of head/disk structure.
Applsci 13 09730 g002
Figure 3. Schematic of the four parts of electromagnetic waves.
Figure 3. Schematic of the four parts of electromagnetic waves.
Applsci 13 09730 g003
Figure 4. Spectral heat flux for a flying height of 2 nm.
Figure 4. Spectral heat flux for a flying height of 2 nm.
Applsci 13 09730 g004
Figure 5. The sign of the imaginary part of the electrical permittivity changes at the frequency.
Figure 5. The sign of the imaginary part of the electrical permittivity changes at the frequency.
Applsci 13 09730 g005
Figure 6. The imaginary part of electrical permittivity for Al 2 O 3 (red line) and Fe (blue line).
Figure 6. The imaginary part of electrical permittivity for Al 2 O 3 (red line) and Fe (blue line).
Applsci 13 09730 g006
Figure 7. Spectral heat flux for varying flying height and frequency.
Figure 7. Spectral heat flux for varying flying height and frequency.
Applsci 13 09730 g007
Figure 8. Graph of heat flux varying with flying height.
Figure 8. Graph of heat flux varying with flying height.
Applsci 13 09730 g008
Figure 9. The distribution of temperature after radiative heating.
Figure 9. The distribution of temperature after radiative heating.
Applsci 13 09730 g009
Figure 10. The deformation caused by the radiative heat flux.
Figure 10. The deformation caused by the radiative heat flux.
Applsci 13 09730 g010
Figure 11. (a) The deformation and (b) the deformation rate under varying flying height.
Figure 11. (a) The deformation and (b) the deformation rate under varying flying height.
Applsci 13 09730 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, Y.; Zhang, G. Thermal Effect of the Back Radiation from Disk to Head after Laser Heating in HAMR. Appl. Sci. 2023, 13, 9730. https://doi.org/10.3390/app13179730

AMA Style

Zhao Y, Zhang G. Thermal Effect of the Back Radiation from Disk to Head after Laser Heating in HAMR. Applied Sciences. 2023; 13(17):9730. https://doi.org/10.3390/app13179730

Chicago/Turabian Style

Zhao, Yu, and Guangyu Zhang. 2023. "Thermal Effect of the Back Radiation from Disk to Head after Laser Heating in HAMR" Applied Sciences 13, no. 17: 9730. https://doi.org/10.3390/app13179730

APA Style

Zhao, Y., & Zhang, G. (2023). Thermal Effect of the Back Radiation from Disk to Head after Laser Heating in HAMR. Applied Sciences, 13(17), 9730. https://doi.org/10.3390/app13179730

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop