France



In recent years, magnetoelectric ΔE-effect sensors have been investigated for the detection of small amplitude and low-frequency magnetic fields that typically occur in many biomedical and diagnostic applications. Such sensors are based on magnetoelectric composite resonators, which consist of mechanically coupled magnetostrictive and piezoelectric materials. They can be processed on a large scale by microelectromechanical system (MEMS) technology with dimensions of a few millimeters down to a few micrometers and are compatible with complementary metal-oxide-semiconductor (CMOS) electronics. Furthermore, they can be operated at room temperature with a broad-bandwidth and a large dynamic range and are robust against microphony and mechanical noise. Such sensors utilize the ΔE-effect of the magnetostrictive material, i.e. the change of the mechanical stiffness tensor with magnetization and hence, with an applied magnetic field.&lt;br/&gt;Since our first publication of this sensor concept in 2011 [1], much progress has been made in understanding the complex interplay of magnetic, mechanical, and electrical properties and their influence on signal and noise of the sensor system. This holds for different designs of magnetoelectric MEMS cantilevers, but also for surface acoustic wave devices [2]. Here, we present recent experimental and theoretical results on the ΔE-effect, the sensitivity, and the noise of magnetoelectric ΔE-effect sensors. The sensors consist of magnetoelectric MEMS resonators with sputter-deposited AlN or Al&lt;sub&gt;73&lt;/sub&gt;Sc&lt;sub&gt;27&lt;/sub&gt;N as a piezoelectric material and soft-magnetic (Fe&lt;sub&gt;90&lt;/sub&gt;Co&lt;sub&gt;10&lt;/sub&gt;)&lt;sub&gt;78&lt;/sub&gt;Si&lt;sub&gt;12&lt;/sub&gt;B&lt;sub&gt;10&lt;/sub&gt;-based single or multilayers.&lt;br/&gt;A magnetoelastic macrospin model of the ΔE-effect in all stiffness tensor components is presented and combined with a finite element electro-mechanical model of the resonator. First and higher order torsion (Fig. 1) and bending modes are compared and analyzed regarding the sensitivity and the ΔE-effect [3]. Measurements and simulations reveal the dependency of the sensitivity on the mode number, caused by the geometry, mode shape, and the ΔE-effect in different stiffness tensor components. Specific rules for the resonator design are derived as a result.&lt;br/&gt;Besides the geometric and magnetic material properties, we find that the resonator loss and the piezoelectric material have a decisive impact on signal and noise. The influence of the quality factor Q [4] and the piezoelectric material [5] on signal and noise are investigated experimentally and explained by combining the previous model with a signal-and-noise equivalent circuit model (Fig. 2). Both quantities, Q and the piezoelectric material show the promising potential of significantly increasing the sensitivity and improving the limit of detection (LOD). Yet, an improvement in the LOD is currently limited by magnetic noise that becomes dominant at large operating voltage amplitudes [6]. Hence, understanding the specific origins of magnetic noise in such magnetostrictive resonant structures is crucial for sensor improvement in the future.&lt;br/&gt;A physical magnetic noise model is presented to identify the origin of stress-amplified, thermal-magnetic noise. It combines large scale finite differences micromagnetism with the mechanical equations of motion to capture magnetostrictive self-energy and frequency effects. Including the piezoelectric constitutive relation and the operating electronics the sensor’s output voltage noise density can be estimated. With the model, we evaluate and discuss experimental results and explain the dependency of magnetic noise on the excitation amplitude and the magnetic bias field. The model provides detailed insights into general noise limits, which arise from using the ΔE-effect as the sensing principle.&lt;br/&gt;Finally, we briefly present first applications in cell imaging [7] and show that small signal measurements with simultaneous localization are enabled by a dual-mode operation scheme with compact self-biased ΔE-effect sensors [8]. This is of particular interest for typical biomedical inverse solution problems where the precise knowledge of the sensor’s position and orientation is essential. Such applications could benefit from sensor arrays that provide many measurements with a high spatial resolution. The noise equivalent model is extended to arrays of parallel sensor elements and compared with measurements. Implications and prospects for sensor developments in the future are discussed in a conclusion.

**References:**

[1] B. Gojdka, R. Jahns, K. Meurisch, et al., APL Vol. 99, p.1-4 (2011); Nature Vol. 480, p.150 (2011)&lt;br/&gt;[2] J. Schmalz, A. Kittmann, P. Durdaut, et al., Sensors, Vol. 20, p.1–17 (2020)&lt;br/&gt;[3] B. Spetzler, C. Kirchhof, E. Quandt, et al., Phys. Rev. Appl., Vol. 12, p.1-13 (2019)&lt;br/&gt;[4] B. Spetzler, C. Kirchhof, J. Reermann, et al, APL., Vol. 114, p.183504 (2019)&lt;br/&gt;[5] B. Spetzler, J. Su, R.-M. Friedrich, et al., APL Materials., invited, submitted (2021)&lt;br/&gt;[6] P. Durdaut, E. Rubiola, J. Friedt, et al., J. Microelectromech. Syst., Vol. 29, p.1347 - 1361 (2020)&lt;br/&gt;[7] N. Lukat, R. M. Friedrich, B. Spetzler, et al., Sensors Actuators A Phys., Vol. 309, 1–8 (2020)&lt;br/&gt;[8] B. Spetzler, C. Bald, P. Durdaut, et al., Sci. Rep., accepted, (2021)

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/1ddee42922239ae2a14790371e63be23.jpg)
Normalized resonance frequencies&lt;sub&gt; &lt;/sub&gt;&lt;i&gt;f&lt;/i&gt;&lt;sub&gt;r&lt;/sub&gt;/&lt;i&gt;f&lt;/i&gt;&lt;sub&gt;r,max&lt;/sub&gt; of the first three torsion modes (TM1-3) as functions of the magnetic flux density &lt;i&gt;µ&lt;/i&gt;&lt;sub&gt;0&lt;/sub&gt;&lt;i&gt;H &lt;/i&gt;applied along the long axis of the cantilever; inset: top photograph of the MEMS cantilever with visible electrodes. A soft magnetic FeCoSiB-Cr multilayer is sputter deposited on the other side.

![](https://assets.underline.io/uploads/markdown_image/1/image/9c28120fbc2c2f322ee01c1b390d8acb.jpg)
Measured and simulated limit of detection (LOD) of two sensors of identical design that differs only in the piezoelectric material, which is AlN and AlScN, respectively. At large excitation voltage amplitudes &lt;i&gt; û&lt;/i&gt;&lt;sub&gt;ex&lt;/sub&gt; magnetic noise imparis the LOD.

INTERMAG 2021

 ΔE-Effect Magnetic Field Sensors

INTERMAG 2021 Welcome Video



	
<figure class="video_container">
 <iframe src="https://screencast-o-matic.com/player/crflqFVnUpZ?v=6&ff=1&title=0&controls=1" width=640 height=360 frameborder="0" allowfullscreen="true"> </iframe>
</figure>

Want to buy a ticket? 

Register for INTERMAG 2021 here

[![](https://assets.underline.io/uploads/markdown_image/1/image/17ce76e7960db633258024eef8d6da7f.png)](https://intermag.org/register)

Got your ticket, but still can’t access the event?

Contact us and our team will look into it.

In recent years, magnetoelectric ΔE-effect sensors have been investigated for the detection of small amplitude and low-frequency magnetic fields that typically occur in many biomedical and diagnostic applications. Such sensors are based on magnetoelectric composite resonators, which consist of mechanically coupled magnetostrictive and piezoelectric materials. They can be processed on a large scale by microelectromechanical system (MEMS) technology with dimensions of a few millimeters down to a few micrometers and are compatible with complementary metal-oxide-semiconductor (CMOS) electronics. Furthermore, they can be operated at room temperature with a broad-bandwidth and a large dynamic range and are robust against microphony and mechanical noise. Such sensors utilize the ΔE-effect of the magnetostrictive material, i.e. the change of the mechanical stiffness tensor with magnetization and hence, with an applied magnetic field. Since our first publication of this sensor concept in 2011 [1], much progress has been made in understanding the complex interplay of magnetic, mechanical, and electrical properties and their influence on signal and noise of the sensor system. This holds for different designs of magnetoelectric MEMS cantilevers, but also for surface acoustic wave devices [2]. Here, we present recent experimental and theoretical results on the ΔE-effect, the sensitivity, and the noise of magnetoelectric ΔE-effect sensors. The sensors consist of magnetoelectric MEMS resonators with sputter-deposited AlN or Al73Sc27N as a piezoelectric material and soft-magnetic (Fe90Co10)78Si12B10-based single or multilayers. A magnetoelastic macrospin model of the ΔE-effect in all stiffness tensor components is presented and combined with a finite element electro-mechanical model of the resonator. First and higher order torsion (Fig. 1) and bending modes are compared and analyzed regarding the sensitivity and the ΔE-effect [3]. Measurements and simulations reveal the dependency of the sensitivity on the mode number, caused by the geometry, mode shape, and the ΔE-effect in different stiffness tensor components. Specific rules for the resonator design are derived as a result. Besides the geometric and magnetic material properties, we find that the resonator loss and the piezoelectric material have a decisive impact on signal and noise. The influence of the quality factor Q [4] and the piezoelectric material [5] on signal and noise are investigated experimentally and explained by combining the previous model with a signal-and-noise equivalent circuit model (Fig. 2). Both quantities, Q and the piezoelectric material show the promising potential of significantly increasing the sensitivity and improving the limit of detection (LOD). Yet, an improvement in the LOD is currently limited by magnetic noise that becomes dominant at large operating voltage amplitudes [6]. Hence, understanding the specific origins of magnetic noise in such magnetostrictive resonant structures is crucial for sensor improvement in the future. A physical magnetic noise model is presented to identify the origin of stress-amplified, thermal-magnetic noise. It combines large scale finite differences micromagnetism with the mechanical equations of motion to capture magnetostrictive self-energy and frequency effects. Including the piezoelectric constitutive relation and the operating electronics the sensor’s output voltage noise density can be estimated. With the model, we evaluate and discuss experimental results and explain the dependency of magnetic noise on the excitation amplitude and the magnetic bias field. The model provides detailed insights into general noise limits, which arise from using the ΔE-effect as the sensing principle. Finally, we briefly present first applications in cell imaging [7] and show that small signal measurements with simultaneous localization are enabled by a dual-mode operation scheme with compact self-biased ΔE-effect sensors [8]. This is of particular interest for typical biomedical inverse solution problems where the precise knowledge of the sensor’s position and orientation is essential. Such applications could benefit from sensor arrays that provide many measurements with a high spatial resolution. The noise equivalent model is extended to arrays of parallel sensor elements and compared with measurements. Implications and prospects for sensor developments in the future are discussed in a conclusion.

**References:**

[1] B. Gojdka, R. Jahns, K. Meurisch, et al., APL Vol. 99, p.1-4 (2011); Nature Vol. 480, p.150 (2011) [2] J. Schmalz, A. Kittmann, P. Durdaut, et al., Sensors, Vol. 20, p.1–17 (2020) [3] B. Spetzler, C. Kirchhof, E. Quandt, et al., Phys. Rev. Appl., Vol. 12, p.1-13 (2019) [4] B. Spetzler, C. Kirchhof, J. Reermann, et al, APL., Vol. 114, p.183504 (2019) [5] B. Spetzler, J. Su, R.-M. Friedrich, et al., APL Materials., invited, submitted (2021) [6] P. Durdaut, E. Rubiola, J. Friedt, et al., J. Microelectromech. Syst., Vol. 29, p.1347 - 1361 (2020) [7] N. Lukat, R. M. Friedrich, B. Spetzler, et al., Sensors Actuators A Phys., Vol. 309, 1–8 (2020) [8] B. Spetzler, C. Bald, P. Durdaut, et al., Sci. Rep., accepted, (2021)

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/1ddee42922239ae2a14790371e63be23.jpg)
Normalized resonance frequencies fr/fr,max of the first three torsion modes (TM1-3) as functions of the magnetic flux density µ0H applied along the long axis of the cantilever; inset: top photograph of the MEMS cantilever with visible electrodes. A soft magnetic FeCoSiB-Cr multilayer is sputter deposited on the other side.

![](https://assets.underline.io/uploads/markdown_image/1/image/9c28120fbc2c2f322ee01c1b390d8acb.jpg)
Measured and simulated limit of detection (LOD) of two sensors of identical design that differs only in the piezoelectric material, which is AlN and AlScN, respectively. At large excitation voltage amplitudes ûex magnetic noise imparis the LOD.

technical paper

Recently, the problem of information storage becomes more serious due to the information explosion. It is widely known that future advances in the areal density of hard disk drives (HDDs) are essentially limited by the magnetic recording trilemma arising from the conflicting requirements between signal-to-noise ratio, writeability, and stability[1]. Energy-assisted magnetic recording (EAMR) such as heat-assisted magnetic recording[2] and microwave-assisted magnetic recording[3] is thought to solve the trilemma by improving writability and enhancing areal recording density. Two-dimensional magnetic recording with shingled writing[4] has also been proposed for increasing track density, which is enabled by partial overlapping of recorded tracks. EAMR combined with SMR is expected to be achieved for the next generation recording schemes of HDDs. On the other hand, servo signal writing is one of serious problems for high areal recording density in hard disks. The self-servo track writing technology, which regenerated from previous written servo patterns, has been adopted during head-disk assembly. However, since the servo signal writing in the present hard disks takes more than a few days per one drives on manufacturing process, and the supply of HDDs can not keep up with the social demand. Magnetic printing with perpendicularly magnetized patterns is a strong candidate for servo track writing with extremely high speed and low cost[5]. However, the recording field strength by utilizing the conventional magnetic printing master medium is not enough to write servo-signals on HDDs with large coercivity for EAMR, and it becomes much weaker as the pattern size or the track pitch decreases. The decrease in the recording field or the field gradient deteriorates the signal-to-noise ratio of servo-signals, and this problem becomes more serious in HDDs which has narrower track width. In order to accurately write servo-signals on HDDs for EAMR combined with SMR, it is essential to enhance the recording field. I propose a new master structure, herein called double magnet master (DMM) medium[6], for magnetic printing technique to improve recording characteristics. Figure 1 shows the schematic illustration of double magnet master medium and recording layer. The DMM medium consists of two magnetic materials with different coercivity. Both materials have perpendicular magnetic anisotropy but exchange interaction between both materials does not act. One material has a coercivity higher than the coercivity of the recording layer, while the other has a lower coercivity. During application of printing field Hp, the magnetization of low coercivity parts (magnet 1) orients to the direction of printing field while high coercivity parts (magnet 2) is not subject to magnetization reversal. As a result, the recording field can be enhanced. The recording field distributions of L/S patterns with pattern length of 10 nm, which corresponds to bit length, were calculated by micromagnetic simulation for conventional perpendicular master (CPM)[5] and DMM media. Tha magnetic spacing is 2 nm in this case. The printing field Hp was varied from 5 kOe to 50 kOe. In CPM, the recording field Hr consists of the magnetostatic field generated from only the magnet 1, while in DMM, the recording field Hr consists of the magnetostatic fields generated from both magnets 1 and 2. The recording field difference ΔHr of about 11 kOe for DMM is twice as that for CPM in various printing field. It was clearly found that the recording field can be enhanced by DMM. Since the magnetostatic field generated from double magnets is effective when the magnetization in the harder magnet does not reverse, the printing field is less than the coercivity of the harder magnet, which is about 73kOe by using L10-FePt patterms as magnet 2. In order to evaluate printing characteristics of CPM and DMM, the printed magnetization distributions of ECC-type EAMR media [7] with the coercivities of 15kOe and 22kOe were calculated. Figure 2 shows printing performance [8] as a function of printing field Hp. The printing performance of 1 means the perfect printing. The optimum printing fields are slightly less than those of the coercivity of recording layer. The printing performances of DMM were dramatically improved comparing with CPM printing characteristics. It was clearly confirmed that the printing performances of DMM are superior to that of CPM. Therefore, the magnetic printing by utilizing DMM is a promising way to realized high-speed servo-track writing for EAMR combined with SMR.

**References:**

[1] H. J. Richter, A. Y. Dobin, J. Magn. Magn. Mater. 287 (2005) 4150. [2] M. H. Kryder, E. C. Gage, T. W. McDaniel, W. A. Challener, R. E. Rottmayer, G. Ju, Y. T. Hsia, M. F. Erden, Proc. IEEE 96 (2008) 1810-1835. [3] J.-G. Zhu, X. Zhu, Y. Tang, IEEE Trans. Magn. 44 Iss. 1 (2008) 125-131. [4] A. Amer, J. Holliday, D. D. E. Long, E. L. Miller, J.-F. Paris, T. Schwarz, IEEE Trans. Magn. 47, Iss. 10 (2011) 3691-3697. [5] N. Sheeda, M. Nakazawa, H. Konishi, T. Komine, and R. Sugita, IEEE Trans. Magn. 45 (2009) 3676-3678. [6] T. Komine, submitted to J. Magn. Magn. Mater. 2021. [7] T. Tanaka, S. Kashiwagi, Y. Kanai, K. Matsuyama, J. Magn. Magn. Mater. 416 (2016) 188-193. [8] T. Komine, T. Murata, Y. Sakaguchi, R. Sugita, IEEE Trans. Magn. 44 Iss. 11 (2008) 3416-3418.

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/d8f8ad00cc0a9b7cd3f35f27da53b2ef.jpg)
Fig.1 Schematic illustration of magnetic printing process by utilizing double magnet master medium.

![](https://assets.underline.io/uploads/markdown_image/1/image/dda58c7ac541a3940c784a3e79047b17.jpg)
Fig.2 Printed magnetization distribution of recording layer and printing performance as a function of printing field.

 Double magnet master media for magnetic printing onto energy-assisted magnetic recording media

The ability to engineer magnetic phase distributions over nanometer length scales is invaluable in the development of spintronic and data storage technologies with enhanced efficiency. Equiatomic FeRh with B2 crystallography is a useful material for local manipulation of magnetic phases, owing to its unusual first-order metamagnetic phase transition (MPT) from antiferromagnetic (AF) to ferromagnetic (FM) ordering, which occurs at ~100 °C [1]. It has been shown that in addition to the thermally driven MPT, there is a disorder-induced MPT from AF to FM ordering, and that with further disordering a second MPT to PM behaviour can be induced [2]. Bennett et al. [3] have demonstrated that by irradiating FeRh with noble gas ions, the disorder-induced MPTs can be exploited to produce a rich variation in the magnetism as a function of depth. In the present work [4], we combine several non-destructive experimental techniques with extensive simulations to identify a direct correspondence between the nature and concentration of induced defects and the resulting modification of the magnetic depth profile in 40 nm FeRh thin films as the fluence of 25 keV Ne+ ions increases. Polarised neutron reflectometry (PNR) was used to determine the magnetisation depth profiles for a range of sample temperatures under an applied field of 1 kOe. The data were fitted to dynamical simulations of the nuclear and magnetic scattering length density (n/mSLD) profiles (Figure 1), which are proportional to the density and in-plane magnetisation, Ms, respectively. The data demonstrate that it is possible to produce a uniform and thermally stable magnetisation with a single-step irradiation, as in Figure 1(b), or to produce a simple bilayer phase structures, as in Figure 1(c) and (d). To elucidate the interplay between the magnetisation depth distributions, the coexistence of magnetic phases, and the ion-induced damage profiles, conversion electron Mössbauer spectroscopy (CEMS) and Doppler broadening variable energy positron annihilation spectroscopy (DB-VEPAS) were performed on each sample. DB-VEPAS allows relative changes of defect concentration to be determined via the line broadening of gamma photons produced from the annihilation of positrons with varying energy Ep. By fitting the DB-VEPAS data, the mean positron diffusion length as a function of increasing fluence is obtained, which is inversely proportional to the concentration of open volume defects. Thus, we demonstrate that the defect concentration monotonically increases with the irradiating fluence. To complement the DB-VEPAS measurements, variable energy positron annihilation lifetime spectroscopy (VEPALS) measurements were also performed on each sample. Here, variations in positron lifetime with increasing fluence signify changes to the size of open volumes; thus, changes in the concentration and the type of defects are determined. By comparing the measured data to theoretical lifetime values calculated using density functional theory (DFT), we identify the predominant defect type for each irradiation condition with depth sensitivity. The PNR, CEMS, and positron data allow the effect of increasing fluence to be explained as shown schematically in Figure 2, wherein the onset of FM and PM phases are associated with the emergence of static disorder and monovacancy defects respectively. In the as-grown case, there is a substantial concentration of monovacancy defects. At the lowest value of the fluence (7.5 × 1013 ions cm-2), the defect concentration increases and sub-unit cell atomic displacements begin to dominate. The generation of static disorder continues as the fluence is increased to 3 × 1014 ions cm-2, with no change to the predominant defect type. On increasing the fluence to 6 × 1014 ions cm-2, atomic displacements give way to monovacancies, with some static disorder remaining at the substrate interface. For all samples, our data are well described assuming AF, FM, and PM phases only, and do not present evidence for states such as spin glasses or magnetic dead layers. This feature, coupled to the fact that the mSLD profiles of the irradiated samples comprise relatively few sublayers with discrete magnetisations which are stable with respect to increasing temperature, highlights the effectiveness of sensitive and controlled disordering of ordered alloy precursors for producing technologically useful magnetic configurations. Our results provide the basis for the patterning of single layer or multilayer structures that incorporate FeRh to create devices on sub-100 nm length scales with AF and FM order in close physical proximity.

**References:**

[1] M. Fallot, Ann. Phys., vol. 10, pp. 291–332 (1938) [2] A. Iwase et al., Nucl. Instrum. Methods Phys. Res., vol. 256, pp. 429–433 (2007) [3] S. P. Bennett et al., Mater. Res. Lett., vol. 6, pp. 106–112 (2018) [4] W. Griggs et al., APL Mater., vol. 8, p. 121103 (2020)

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/69319fe9359c4803c0c33b75650a7a95.jpg)
Fig. 1. (a)-(d) Magnetic (coloured) and nuclear (grey) scattering length density (mSLD and nSLD) profiles corresponding to the fitted PNR data for each sample across the range of measured temperatures.

![](https://assets.underline.io/uploads/markdown_image/1/image/85121693078388560eeaf904d351ce49.png)
Fig. 2. The effect of increasing the irradiating fluence on the depth-dependent magnetisation, magnetic phase fractions, and concentration and type of structural defects is illustrated.

 Depth Selective Magnetic Phase Coexistence in FeRh Thin Films

As the current recording media for the heat assisted magnetic recording (HAMR), the granular L10-FePt thin film needs to have uniform grain size with small pitch distance and well-defined grain boundaries. To achieve the targeted microstructure with proper L10 ordering using sputtering techniques, FePt along with grain boundary materials (segregants) are co-deposited on MgO underlayers at elevated temperature. The columnar growth of FePt is desired to have sufficient grain volume without grain size expansion. The fabrication of grains at small grain pitch with well-defined grain boundaries has been proven to be challenging.[1] Over the past decade, segregant materials, such as carbon, have attracted extensive attention for producing isolated FePt grains with well-defined grain boundaries but has experienced difficulty at reducing grain pitch distance and making grains with large aspect ratio. Though the columnar growth of FePt can be achieved by using some oxide segregant materials, such as TaOx[2] and TiOx[3], the resulting films always show poor L10 ordering and poor magnetic properties. Proper ordering of FePt can be achieved with SiOx[4], but the resulting microstructure tends to be maze-like. Therefore, exploiting multilayer deposition is a promising strategy to achieve the desired microstructure with proper ordering and magnetic properties. In this study, we utilize a dual-layer strategy of FePt-BN/FePt-SiOx on polycrystalline MgO underlayer in order to obtain the desired media microstructure with good L10 ordering and magnetic properties. The film stack of Ta(5)|Cr(40)|MgO(6)|Media(6.5) was deposited on Lotus NXT glass substrate with the numbers in parentheses in nm. Lotus NXT glass is one of the promising glass substrates for high-temperature deposition of HAMR media. The Ta and Cr layers were DC sputtered at room temperature and 280°C, respectively. The samples were subsequently annealed at 650°C for 1 hour in order to promote the texture and grain size of Cr buffer layer. The MgO underlayer was RF deposited on Cr buffer layer after the sample cooled down to room temperature. There are two magnetic layers in the media layer, where FePt was co-sputtered with BN and SiOx, respectively, for the first (M1) and the second (M2) magnetic layer. The deposition temperatures of the M1 layer and the M2 layer were fixed at 700°C and 550°C, respectively. From the microstructural point of view, the main objective of the M1 layer is to obtain well-isolated FePt grains to set the grain center-to-center pitch distance, whereas the M2 layer aims to suppress the lateral growth and to promote the columnar growth of the FePt grains. Fig. 1 (a) shows the in-plane HAADF image of the film stack Lotus NXT glass |Ta(5)|Cr(40)|MgO(6)|Media(6.5) with its (b) grain size distribution profile and (c) grain center-to-center pitch distance distribution profile. The figure demonstrates well-isolated FePt grains with mean grain size of 6.3 nm and mean pitch distance of 8.3 nm. As shown in the XRD pattern of this media in Fig. 1 (d), there are only (001) and (002) peaks from FePt, indicating the c-axis of FePt grains are predominantly perpendicular to the film plane. The broad peak at around 23.5° corresponds to the amorphous peak from Lotus NXT glass substrate. The order parameter (S) considering the relative integrated intensity ratio between (001) and (002), I001/I002, with the absorption factor and Lorentz factor is 0.7, suggesting good chemical ordering in FePt grains. Fig. 1 (e) shows the perpendicular and in-plane hysteresis loops of this media. The perpendicular coercivity (Hc⊥) is 26 kOe. The in-plane loop shows a small opening with the coercivity (Hc//) of 2.9 kOe and the normalized remanence magnetization (Mr//) of 0.08. The remanent ratio, which is Mr// /Mr⊥, is 0.09, suggesting a minor in-plane component of FePt grains in the sample. By exploiting the dual-layer structure of FePt-BN/FePt-SiOx, we have demonstrated well-isolated columnar FePt grains with small mean pitch distance of 8.3 nm. The XRD pattern shows that the grains are well-textured and highly ordered. As a result, the perpendicular coercivity reaches 26 kOe. The utilization of multilayer deposition strategy with segregants BN and SiOx provides a promising technique to achieve desired HAMR media microstructure without sacrificing FePt chemical ordering and magnetic properties.

**References:**

[1] Y. Kubota et al., IEEE Trans. Magn., vol. 54, no. 11, p. 3201206, 2018 [2] B. Lim et al., J. Appl. Phys., vol. 105, no. 7, p. 07A730, 2009 [3] B. S. D. C. S. Varaprasad et al., IEEE Trans. Magn., vol. 49, no. 2, p. 718–722, 2013 [4] H. Ho et al., Appl. Phys. Lett., vol. 102, no. 11, p. 112411, 2013

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/44d762af9dd6d7af7c27e30c471fccc1.jpg)
Fig. 1 (a) In-plane HAADF image of the film stack Lotus NXT glass |Ta(5)|Cr(40)|MgO(6)|FePt-38vol.% BN(2)|FePt-35vol.%SiOx with its corresponding (b) grain size distribution and (c) grain center-to-center pitch distance. Its XRD pattern and magnetic properties are in (d) and (e) respectively.

 Fabrication of FePt-BN/FePt-SiOx dual-layer structure for HAMR media on Corning Lotus™ NXT glass substrate

I. INTRODUCTION In recent years, spin-torque transfer magnetic random access memory (STT-MRAM) has received tremendous attention due to its low switching energy, fast read/write speed, and high scalability [1]. However, thermal fluctuation and process variation have a detrimental effect on the data recovery of STT-MRAM, leading to both read errors and write errors [2]. Error correction codes (ECCs), such as the Hamming codes and Bose-Chaudhuri-Hoquenghem (BCH) codes have been employed to improve the reliability of STT-MRAM [2, 3]. Meanwhile, the channel quantizer that quantizes the signal read back from the STT-MRAM cell is critical to support the above ECCs, since high-precision analog-to-digital converters (ADCs) are not applicable for STT-MRAM. Various criteria have been proposed for designing the quantizer for STT-MRAM, including the Maximizing-Mutual-Information (MMI) criterion, the Maximizing-Cutoff-Rate (MCR) criterion, and the Optimizing- Polyanskiy-Poor-Verdu -Bound (PPVB) criterion [4]. However, the final choice of the quantization scheme is still relying on the decoding error rate performance of ECCs. It is typically evaluated by using computer simulations, which are too slow to reach the target error rate level of 10-12 or below of data storage systems. In this work, we propose a union bound analysis which can accurately predict the word error rates (WERs) of ECCs with maximum-likelihood (ML) decoding over the quantized STT-MRAM channel. The derived bound provides a theoretical tool for comparing different quantization schemes at the targeted low error rate levels of data storage systems. Moreover, we also present a new criterion to design the quantizer by minimizing the WERs of ECC decoding that are obtained from the union bound analysis. Numerical results show that the union bound optimized quantizer can achieve better WER performance than the state-of-art quantizers for STT-MRAM. II. UNION BOUND ANALYSIS AND CHANNEL QUANTIZATION The union bound analysis is a theoretical technique that is used for bounding the error rate of ML decoding of ECCs. It was mainly developed for codes over symmetric channels [5]. In [6] and [7], more sophisticated union bound techniques were proposed for simple asymmetric channels, i.e. the Z-channel and the general binary asymmetric channel (BAC). In this work, we consider the quantized STT-MRAM channel with both write errors and read errors, which can be modelled as a concatenation of a BAC and a binary-input multiple-output discrete memoryless channel. Our proposed union bound analysis first utilizes a function of a predefined multi-dimensional distance between two codewords to calculate the pairwise error probability (PEP) [7] of the quantized STT-MRAM channel. Note that this is different from the cases of the symmetric channels where the PEP is uniquely determined by the Hamming distance of the code. Based on the PEP, we can derive a union bound of the WER of ML decoding as a function of the weight distributions of ECC. Next, we further propose to use the WER obtained from the above union bound analysis as the criterion to design the channel quantizer of STT-MRAM. We also develop effective differential evolution algorithm to determine the optimum quantization boundaries for the multi-bit quantizers. Since our proposed union bound optimized quantizer takes into consideration the weight spectrum of ECC, it can achieve better error rate performance than the prior-art quantizers. III. NUMERICAL AND SIMULATION RESULTS The system parameters of STT-MRAM are taken from [2]. In the quantized STT-MRAM channel model, we assume the write error rate of 10-5, and vary the mean normalized resistance spreads σ0/μ0=σ1/μ1 of the low and high resistances to account for the influence of different process variations. The (72,64) extended Hamming code [3] is adopted to validate our analysis, and the simulated ML decoding performance is obtained using the Bahl, Cocke, Jelinek and Raviv (BCJR) decoder [5]. In Fig.1, we compare the decoding WERs provided by the union bound with those achieved by the BCJR decoder through simulations, with the proposed union bound optimized quantizer with different number of quantization bits q. Observe that there is a noticeable gap between the union bound and simulated WERs for the case of q=1. However, with two or more quantization bits, the performance gap becomes negligible, and when the WER is less than 10-6, the union bound is almost identical to the simulated WER. In Fig. 2, we fix the quantization bit to be q=2, and compare the WER performance with different types of quantizers. A nice agreement is observed again between the WERs predicted by the union bound and those obtained from simulations, for the various types of quantizers. Moreover, our proposed union bound optimized quantizer outperforms all the other quantizers. There is a larger performance improvement over the MMI quantizer, which is most widely adopted in the literature. This demonstrates the potential of our proposed union bound analysis in predicting the decoding error rate performance and in guiding the design of channel quantizer for STT-MRAM. Acknowledgement: This work is supported by RIE2020 Advanced Manufacturing and Engineering (AME) programmatic grant A18A6b0057 and Singapore MOE Tier 2 fund MOE2019-T2-2-123.

**References:**

[1] M. Marinella, “The future of memory,” in Proc. Aerospace Conf, Big Sky, MN, USA, Mar. 2012, pp. 1–11. [2] K. Cai, and K. A. S. Immink, “Cascaded channel modeling, analysis, and hybrid decoding for spin-torque transfer magnetic random access memory (STT-MRAM),” IEEE Trans. Magn., vol. 53, no. 11, pp. 1–11, 2017. [3] K. Cai, Z. Qin, and B. Chen, “Channel capacity and soft-decision decoding of LDPC codes for spin-torque transfer magnetic random access memory (STT-MRAM),” in Proc. Int. Conf. Comput., Netw., Commun. (ICNC), San Diego, CA, USA, Jan. 2013, pp. 550–554. [4] Z. Mei, K. Cai and L. Shi, “Information Theoretic Bounds Based Channel Quantization Design for Emerging Memories,” in Proc. IEEE Information Theory Workshop (ITW), Guang Zhou, China, Nov. 2018, pp. 1–5. [5]W. E. Ryan and S. Lin, Channel Codes: Classical and Modern, Cambridge, Cambridge University Press, 2009. [6] A. Barbero, P. Ellingsen, S. Spinsante, and O. Ytrehus, “Maximum likelihood decoding of codes on the Z-channel,” in Proc. IEEE International Conference on Communications (ICC), Istanbul, Turkey, Jun. 2006, pp. 1–5. [7] G. Song, K. Cai and J. Cheng, “A union bound analysis for codes over binary asymmetric channels,” in Proc. IEEE International Conference on Communications (ICC), Paris, France, May. 2017, pp. 1–5.

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/7687adc035569fbf9f2d0102c1f4556f.jpg)
Fig. 1. Comparison between analytical and simulation WERs, with union bound optimized quantizer with q= 1,2,3.

![](https://assets.underline.io/uploads/markdown_image/1/image/ff97a646ffd6d77b07f4e57f875be682.jpg)
Fig. 2. Comparison between analytical and simulation WERs with different types of quantizers, for q=2.

 Union Bound Analysis for Spin-Torque Transfer Magnetic Random Access Memory with Channel Quantization

Two-dimensional magnetic recording (TDMR) with laterally-displaced multiple read heads has been developed to alleviate the impact of inter-track interference (ITI) and track edge effects. However, the cross-track (XT) position of the displaced 2nd head varies significantly with the skew angle of head assembly from ID to OD across disk radius [1]. As a result, each disk platter is divided into hundreds of groups of sectors, or zones, to accommodate the skew angle variations. The 2-dimentional equalization technique used in conventional data detection channels requires different channel parameters at different skew angles/zones for channel optimization, and the number of different channel parameter sets across the disk platter are in hundreds. Furthermore, resulting performance gain is further diminished by random fluctuations such as non-repeatable disk runout. In this situation, if we employ machine learning technique on TDMR data detection, one could imagine the possibility of achieving optimal performance in resolving ITI and track edge effects across the entire disk surface while using “learning” to replace relatively complex parameterization of physical modeling. In our previous work, we have developed deep neural network (DNN) for disk drive data detection channels with performance (detection error-rate vs. signal-to-noise ratio) slightly more superior than the state-of-the-art disk drive channels [2,3]. In the work presented here, we train a convolution neural network (CNN) in TDMR system to detect data under a wide range of ITI from different head skew angles. We demonstrate that in the medium noise dominant situation, CNN can learn all ITI in the training dataset and completely recover ITI-free BER without any physical modeling. Our eventual goal is to explore if we can train a single CNN detector to learn and optimize for all head skew angles/zones automatically, thereby lifting the complexity of designing different channel parameter sets at each head skew angle. The recorded tracks and CNN detection scheme are shown in Fig. 1. Two adjacent data tracks are written with random data and arbitrary phase. The width of the track is denoted TW, and no erase-band nor transition curvature is included. Two XT overlapped heads of width greater than TW are used. Here we simplify the modeling and assume each head only captures signals under its XT coverage. Since the relative positions of head-1 and head-2 are locked, we characterize ITI of the reader array with a single variable: the ratio between head-1 coverage of Track-II (Δx1) versus the track width (TW): δ1=Δx1/TW. Here we assume the ITI stays fixed within each zone, and different zones will have different ITI. In this work, we only investigate the medium noise dominated situation: the jitter noise is modeled as the uncertainty of the transition locations following a Gaussian distribution. Since both heads “see” the same magnetic patterns from the two tracks, the jitter noise is correlated noise. A CNN is used as the data detector in this work. It receives signals from both heads as input, and we only train it to detect data for Track-I. The rest of the parameters are listed in Fig. 1 caption. To simulate reading data at different skew angles, we evaluate the CNN detection BER with signals of 11 different ITI from δ1=0.2 to 0.8 as in Fig. 2. The black dashed line benchmarks the ITI-free BER. For the blue curve, the CNN detector is trained using signals from one skew angle/zone with a fixed ITI, δ1=0.5 . In this case, CNN can only eliminate the trained ITI, and the detection accuracy degrades significantly whenever the testing ITI is different. In contrast, the red curve plots the BER when we train the CNN using signals with mixed ITI same as the testing dataset. This time CNN can successfully eliminate ITI for all testing head skew angles/zones. Our results suggest that CNN has the capacity to learn very complicated correlations and is solely data-driven. With the learning process, the CNN-based ML channel presented here is capable of resolving the interference from ITI infested signals, automatically rejecting them, and arriving at correct detection. This approach is fundamentally different from the conventional channel design that relies on sophisticated physical modeling and channel parametrization. In addition, the fact that a single CNN can learn a wide range of ITI states that it can optimize the channel performance for head skew angles using a single CNN eliminating the need to use multiple channel parameter sets. Not only this machine learning channel can eliminate the complex physical characterization process required for conventional approach, more importantly, the TDMR performance gain can be fully realized without compromising. The potential of fully eliminating ITI also presents the possibility to have reader width wider than track width which could be critically important as track pitches in HDDs continue to decrease.

**References:**

[1] Suzuto R, Nakamura Y, Osawa H, et al. Effect of reader sensitivity rotation in TDMR with head skew. IEEE Transactions on Magnetics, 2016, 52(7): 1-4. [2] Qin Y, Zhu J G. Deep Neural Network: Data Detection Channel for Hard Disk Drives by Learning. IEEE Transactions on Magnetics, 2020, 56(2): 1-8. [3] Qin Y, Zhu J G. Automatically Resolving Inter-Track Interference with Convolution Neural Network Detection Channel in TDMR. IEEE Transactions on Magnetics, 2020.

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/1636e831501326557b33b2a856085d21.jpg)
CNN detector in TDMR. The heads have DT sensitivity of PW50/bit-length = 2.5. The green(yellow) signals are from head-1(head-2), and they are sampled 1 bit/period. The jitter width is assumed 1.4nm in this work. The CNN channel has 3 hidden layers (blue) with 50 filters each. The kernel sizes are 3-3-3-1. ReLU is used as the non-linear activation function, and cross-entropy is the loss function.

![](https://assets.underline.io/uploads/markdown_image/1/image/3a69ad272e492b195d37c70a85e1048b.jpg)
Detection BER at different zone/ITI. Black dashed curve: ITI-free detection BER. Blue: CNN is only trained with signals of fixed ITI δ1=0.5 . Red: CNN is trained with signals with 11 different ITI (same as the testing dataset).

 Enable TDMR Gain with Convolution Neural Network Based Machine Learning Algorithm

I. Introduction Track misregistration (TMR) and two-dimensional (2-D) interference that consists of an inter-track interference (ITI) and inter-symbol interference (ISI) are the significant problems in an ultra-high bit-patterned media recording system (BPMR) which easily degrade the bit-error-rate (BER) performance of the system, because of the extremely narrow distance between bit-islands. It very challenges to control the reader to fly above the center of the desired track. Normally the TMR effect can be controlled by servomechanism [1] which required extra area to store the redundancy bits for tracking the data track. Instead, several works use only readback signals to estimate and correct the TMR effect [2-4]. For example, we proposed a TMR prediction method by using the energy ratio of readback signals, and then the TMR effect can be dealt with utilizing the appropriate 2-D equalizers together with the use of soft-information adjustment (SIA) technique [5]. II. Proposed Methods To improve the TMR estimation accuracy and BER performances of BPMR system; therefore, we propose to use the neural network models called a multi-layer perceptron (MLP) in both TMR level prediction and TMR effect correction processes. We consider the dual-reader two-track reading under BPMR media at an areal density (AD) of 3.0 Tbits/in2, which track-pitch, Tz and bit-period, Tx are set to be 14.5 nm. In the reading process, 2 readers simultaneously read the medium of upper- and lower-tracks, while the sidetracks located above and below two main tracks will be also read, which denoted as r-1,k, r0,k, r1,k, and r2,k. Then, 4 readback signals are equalized by 1-D (one-dimensional) equalizers which are designed according to signal with no TMR effect. The equalized data sequences will then be multiplied with their absolute value for each data sample to extend their amplitude before sending them to the TMR level prediction and TMR effect correction processes, respectively. In TMR level prediction, four modified readback signals, u-1,k, u0,k, u1,k, and u2,k with the size of 4×32,760 bits will then be grouped into 4×3 bits in each group before sending to the TMR level prediction process. That means each considered group consists of 12 data bits and so there are 10,920 groups in one set of readback sequence, the number group represented by j ∈ {1, 2, 3, …, 10920}. Fig. 1(a) shows the structure of the MLP TMR estimator which consists of 1 input layer and 12 input nodes, 5 hidden layers consist of 100 hidden nodes, and 1 output layer with a single node. This MLP TMR estimator receives a group of readback signals at a time for predicting the estimated TMR level, TMR'j of the jth group. Then, the estimated TMR levels of all groups were averaged to obtain the estimated TMR level, TMR' of the entire sequence. Here, the accuracy percentage of TMR' is measured as accuracy (%) = 100 – (|TMR'-TMR| / TMR × 100) [5]. we also propose to use the MLP data detector which may be employed instead of the conventional detectors-based PRML system. Fig. 1(b) shows the structure of the MLP data detector, which is identical to MLP TMR estimator except for the number of input and output nodes. The MLP data detector receives a group of modified readbacks similar to the MLP TMR estimator, but it also receives an additional estimated TMR level, TMR' that was produced from MLP TMR estimator. Here, it is fixed as one of the input data of the MLP data detector. Finally, we can obtain the estimated user bits of both upper- and lower-tracks, i.e., a'0,k and a'1,k, in the size of 2×3 bits for all 10,920 groups. In this study, the leaky rectified linear unit (ReLU) is used as an activation function for the hidden nodes and linear function in the output nodes for both MLP models. The linear regressions are then adopted for calculating their outputs. The Adam algorithm is used to optimize the weights and bias of both MLP models. III. Results and Discussion In this study, we start with the consideration of an accuracy percentage of TMR level estimation in various TMR and SNR levels where the position jitter is given to be 0%. The TMR level is varied from 0% to 15% of a track pitch as mentioned above. The signal-to-noise ratio (SNR) that is defined as 10log10(1/σ2), where σ is a standard deviation of electronics noise, is first considered at 10, 15, 20, and 25 dBs. We found that the proposed MLP TMR estimator can provide an accuracy percentage higher than 80% for all 4 considered SNRs and all considered TMR levels. Especially, at the higher SNR, e.g., 15, 20, and 25 dBs, we can get more than 90% accuracy for all TMR levels. Moreover, we also investigate the BER of our proposed MLP detector, which is compared with the convention partial response maximum likelihood (PRML) system that performs based on 2-D equalizers and 2-D Viterbi detectors together with the SIA technique [5] where it is represented with “conventional”. As shown in Fig. 2, it seems that at higher SNR, our proposed system can provide better performances for all TMR levels. For example, at BER = 10-4 and TMR level is 0%, our proposed system can provide more than 5 dB over the conventional system. The system performance implies that it is possible to use the MLP neural networks for signal processing in BPMR system.

**References:**

[1] Y. B. Chang, D. K. Park, and N. C. Park, IEEE Trans. Magn., Vol. 38, no. 2, p. 1441 (2002) [2] C. Warisarn, IEEE Trans. on Magn., Vol. 55, no. 7, p. 1 (2019) [3] W. Busyatras, C. Warisarn, and Y. Okamoto, AIP Adv., Vol. 7, p. 056501 (2016) [4] S. Han, G. Kong, and S. Choi, IEEE Transactions on Magnetics, Vol. 55, no. 7, p. 1 (2019) [5] K. Kankhunthod, W. Busyatras, and C. Warisarn, 17th ECTI-CON., p. 735 (2020)

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/09de0a5a9c8f4396146c47ebd77277d5.jpg)
Fig. 1: MLP structures of the TMR level estimator and TMR detector.

![](https://assets.underline.io/uploads/markdown_image/1/image/fd1a13ccd5cbab14c903baba723d02af.jpg)
Fig. 2: BER performance comparison between the proposed MLP based TMR correction method versus the conventional method at the AD of 3.0 Tbits/in2, here the position jitter is fixed to be 0%.

 Multilayer Perceptron Based Method for Track Misregistration Correcting in Dual-reader Two-track Reading BPMR Systems

I. INTRODUCTION In the head-disk interface (HDI) of heat-assisted magnetic recording (HAMR), a laser is employed over the read/write area to locally heat the recording media to its Curie temperature to facilitate data writing. During this process, the lubricant on the disk undergoes harsh cyclic thermal processes so that the lubricant is partially depleted by the laser [1]. Thus, the lubricant thickness is locally reduced during the writing, but afterwards it recovers gradually due to lubricant reflow [2]. In this study, the effect of the lubricant thickness on the nanoscale heat transfer across the HDI is studied using a thermal fly-height control (TFC) heater and an embedded contact sensor (ECS). The dR/dP, the rate of change in the ECS resistance with respect to the TFC power, was measured to indicate the HDI heat transfer coefficient [3]. It is demonstrated that the ECS is a good thermal sensor for the lubricant thickness through the contact cooling effect. II. EXPERIMENT SETUP Experiments were performed on a component-level spin-stand stage. When flying over the HAMR disks with the lubricant thickness from 8.0 Å to 12.5 Å, the TFC heater was energized until head-disk touchdown occurred. During the touchdown process, the ECS, biased with a constant DC current, was used to monitor the head temperature. An acoustic emission (AE) sensor was used to detect the head-disk contact. The ECS resistance versus the TFC power relation (R-P) was measured. Then the R-P curve’s derivative, dR/dP, was obtained to indicate the head-disk nanoscale heat transfer. III. RESULTS AND DISCUSSION Fig. 1 shows the results of dR/dP versus backoff TFC power using different ECS bias currents, where zero backoff power is the TFC power at the head-disk contact onset detected by the AE sensor. The figure presents the contact cooling effect around +1 mW backoff TFC power, which appears as the local minimum. It is seen from the figure that the cooling effect is stronger with a higher ECS bias current, since the hotter ECS undergoes a larger temperature drop when touching the rotating disk surface. Thus, the ECS biased with a higher current delivers better sensitivity for the experiments. Fig. 2 shows the dR/dP curves for the disks with the lubricant thickness varying from 8.0 Å to 12.5 Å. As the lubricant thickness increases, the dR/dP local minimum rises and hence the contact cooling effect weakens. Therefore, the lubricant behaves as a thermal barrier to the head-disk heat transfer and its effect can be quantified by the ECS. IV. CONCLUSIONS In this digest, touchdown experiments are introduced to quantify the effect of the lubricant on the nanoscale heat transfer at the head-disk interface. The experimental results show that a higher ECS bias current ensures a better sensitivity for the experiments and that the lubricant layer weakens the head-disk heat transfer. Taking advantage of this effect, the ECS can be further used to measure the local lubricant thickness depletion and reflow during HAMR operations.

**References:**

[1] M. Kryder, E. Gage, T. McDaniel, Proc. IEEE, 96(11), 1810-1835 (2008). [2] H. Wu, A. Mendez, S. Xiong, J. Appl. Phys., 117(17), 17E310 (2015). [3] N. Zuckerman, J. Kiely, M. Anaya-Dufresne, IEEE Trans. Magn., 55(3), 1-5 (2018).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/90768b5c5efe6453870c0a23b326801e.jpg)
Fig. 1. The dRdP versus backoff TFC power using different ECS bias currents.

![](https://assets.underline.io/uploads/markdown_image/1/image/3d1ddabf5d889ce0eb44baa841bc9121.jpg)
Fig. 2. The dRdP versus backoff TFC power using different lubricant thicknesses (ECS bias: 2.50 mA).

 Effect of Lubricant Thickness on the Nanoscale Heat Transfer at the Head-Disk Interface

The current exponential growth in data in combination with the slow-down in areal density and cost scaling of hard disk drives (HDDs) is fueling demand for cost effective storage technology. Magnetic tape storage has a very low total cost of ownership and is suited for the storage of less frequently accessed data where tape's higher access latency is an acceptable trade-off for the reduced cost. Today's best of breed tape drives provide a 20 TB native cartridge capacity at an areal density of 11.7 Gbit/in2. This is about two orders of magnitude lower than the areal density of state-of-the-art HDDs and implies that tape has significant potential for continued scaling. This is in line with the projections of the Information Storage Industry Consortium (INSIC) Tape Roadmap that projects a 34% annual growth in tape areal density leading to a predicted 278 Gb/in2 in 2029 [1]. In this work we describe a tape areal density demonstration using a prototype perpendicularly oriented strontium ferrite (SrFe) tape that shows the potential of recording at 317 Gbit/in2. Key properties of the media include: a particle volume of 950 nm3, a particle coercivity of Hc=208 kA/m, a perpendicular squareness ratio of 0.87, a surface roughness Ra=1.1 nm and Rz=17 nm measured by AFM, and a tape thickness of 4.3 μm. We investigated the recording potential of the SrFe media by writing repeating 255-bit pseudo-random binary sequences (PRBS) at linear densities ranging from 500 to 768 kbpi. The data was read back using a 29 nm side shielded tunneling magneto resistive (TMR) reader from an HDD and the captured readback signal was processed offline using a software read channel. Fig. 1(a) shows the average SNR of the captured data measured at the input of an EPR4 detector and plotted versus linear density. Fig. 1(b) plots the byte-error rate (BER) versus linear density, obtained by processing the data using three detection algorithms: 1) 8-state EPR4 detection, 2) adaptive 16-state noise-predictive maximum-likelihood (NPML) detection and 3) extended 64-state data-dependent NPML detection that tracks the first and second order statistics of the data-dependent noise (D3-NPML). Each data point results from processing about five million bits. To set a target threshold for the BER at the output of the detector, we analyzed the performance of the Reed-Solomon product code used in state-of-the-art tape drives combined with the iterative decoding scheme described in [2]. This analysis indicates a post-error-correction-coding byte-error rate (UBER) of less than 1e-20, equivalent to that of an enterprise class tape drive, can be achieved with two full decoding iterations if the BER at the output of the detector is &lt; 4.5e-2. Hence, we take a BER &lt; 4.5e-2 as the target threshold. The D3-NPML detector provides significantly better performance than the NPML detector and crosses the BER threshold of 4.5e-2 at a linear density of 750 kbpi. At a linear density of 702 kbpi the D3-NPML detector provides a BER of 2.8e-2 and the NPML detector also has a BER performance below the target threshold. The average SNR at 750 kbpi and 702 kpbi are 9.9 dB and 10.6 dB, respectively. Thus, the D3-NPML detector achieves our target 1e-20 UBER at an operating point of 702 kbpi with an SNR margin of about 0.7 dB. In the area of track following, we made a set of advances to enable the rapid track-density scaling projections of the INSIC roadmap. First, we developed a new timing-based servo pattern that is shown in Fig. 2(a). Second, we developed a field programmable gate array prototyping board that implements four servo channels and a quad channel averaging scheme where the tape head position is calculated by combining position estimates from the readback signals of four servo patterns, see Fig. 2(b). Third, we enhanced the tape transport described in [2] with two 20 mm diameter air bearing tape guides and a prototype track-following actuator that enables a lower spacing between the tape guides adjacent to the head. This decreased tape guide spacing reduces compressional wave disturbances that can degrade the track-following performance. Fourth, we developed a low friction tape head that reduces friction by about 75% compared to the heads used in our previous work [2],[3] and further reduces compressional wave disturbances. Finally, we used the H-infinity framework to design a set of tape velocity optimized track-following controllers. Combining these technologies, we achieved a position error signal (PES) characterized by its standard deviation of σPES ≤ 3.18 nm over the tape velocity range of 1.2 to 4.1 m/s as shown in Fig. 2(c). We estimate the minimum reliable track width supported by this tracking accuracy with the INSIC model used in our previous work [2], [3], i.e.: track width = 2 x (2)0.5 x 3σPES + reader width. Taking the 29nm width of the reader used for the recording experiments and a worst case σPES = 3.2 nm results in a 56.2 nm track width and a track density of 451.9 ktpi. Combined with the linear density of 702 kbpi demonstrated with the 29 nm reader enables a potential areal recording density of 317 Gb/in2.

**References:**

[1] Information Storage Industry Consortium (INSIC) 2019-2029 Tape Roadmap, http://www.insic.org/wp-content/uploads/2019/07/INSIC-Technology-Roadmap-2019.pdf [2] S. Furrer et al., "85.9 Gb/in2 recording areal density on barium ferrite tape ", IEEE Trans. Magn., vol. 51, no. 4, Apr. 2015, Art. ID 3100207. [3] M. A. Lantz et al., "123 Gb/in2 recording areal density on barium ferrite tape ", IEEE Trans. Magn., vol. 51, no. 11, Nov. 2015, Art. ID 3101304.

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/ea49499647aa7c96c27439679edc4096.jpg)
Fig 1. (a) SNR vs. linear density, (b) BER vs. linear density for three detectors. The dashed line in the BER plot indicates the 4.5e-2 target byte-error rate threshold.

![](https://assets.underline.io/uploads/markdown_image/1/image/1705d36ba0e3724ce22ce7522c845fd5.jpg)
Fig 2. (a) Geometry of the new servo pattern, (b) Read back signals from four servo tracks formatted with the new servo pattern, at 2.15m/s tape velocity, (c) Standard deviation of the PES versus tape velocity.

 317 Gb/in2 Recording Areal Density on Strontium Ferrite Tape INVITED

Introduction and Background 60 Zettabytes (ZB or 10^21 Bytes) will be created in 2021 and grow to 160 ZB by 2025.(1) To store this data, Hard Disk Drives (HDD) areal density must continue to increase from current products near 1.1Tbpsi (2) to the demonstrated AD of 2.77Tbpsi (3) and beyond. To support this growth, reader technology must also advance to higher AD with manufacturable solutions. Remarkably, HDD reader technology still utilizes the Spin Valve reader concept that was first introduced in 1997 at ~ 1 Gbpsi. This 1000x increase in reader AD capability has been achieved by scaling of critical dimensions, the introduction of TMR, advances in stack materials/ processing and utilization of multiple read sensors. The extendibility of the SV TMR design and supporting innovations has been remarkable - but likely not without limits. The HDD industry’s Advanced Storage Research Consortium (ASRC, http://asrc.idema.org/) has sought to understand and address read sensor requirements for HDD Areal Densities (AD) up to 4Tbpsi and higher. In this paper we release the ASRC reader roadmap. We identify critical challenges and key gaps for the SV TMR type reader to support 4Tbpsi. We then outline what may be required for &gt; 4Tbpsi and explore how to proceed. Approach and Methods We begin by seeking to understand the potential extendibility of the Reader Spin Valve design with TMR stack. We partition the problem between Resolution and SNR. Multiple reader design geometries are considered. (Figure 1.) The reader design resolution capabilities were determined using a new 2D micromagnetic reciprocity approach. The Bottom SV and derivative Recessed AFM design appear to be sufficient for 4Tbpsi. To understand the reader SNR requirements, one must begin with the recorded data that one needs to readback. A top-down system level approach and bottom-up micromagnetic modeling approach are combined and iterated to identify the suitable reader requirements. A RSS model is used to estimate reader requirements between 2 and 6 Tbpsi.(4) A combined finite difference micromagnetic model and analytical reader signal-noise model is used to perform DOE’s to determine the closest fit solutions. Results and Discussion The results of this study show that reader widths must shrink to near 10nm and shield-to-shield spacing approach 15nm to provide cross-track and down-track resolution to support track densities near 1130 kTPI and linear densities &gt; 3600 kFCI for ~ 4Tbpsi. Meeting these targets will demand improved process capabilities and new innovations in device designs. SNR requirements are used to project TMR reader stack magnetoresistance (MR) and resistance-area (RA) product targets. Within reasonable assumptions for sensor bias, stability and material properties, we show stack requirements are well in excess of 100% MR at 0.10 uOhm-cm RA. (Figure 2) This requirement is significantly more challenging than previously projected and exceeds best reports to date for TMR readers (5). In addition to seeking high signal, it is also vital to understand and manage the noise components. We show thermal magnetic noise is a growing concern but we further show that other contributions such as preamp noise could dominate as reader dimensions shrink and resistance grows. Roadmap Challenges and Gaps An overall reader spec is provided with proposed geometric, magnetic, and electronic targets. We also explore potential future designs to meet requirements &gt; 4Tbpsi. A summary of key issues, possible solutions and important research topics is reviewed. We hope this study can illuminate and motivate critical research into these key reader technologies. Meeting these challenges will be critical to continued HDD AD growth to 4Tbpsi and beyond.

**References:**

 1) Data Age 2025, Dec 2018 https://www.seagate.com/our-story/data-age-2025/ 2) Forbes, July 2020 https://www.forbes.com/sites/tomcoughlin/2020/07/08/western-digital-data-center-16-20-tb-hdds/?sh=4f56c3f87ba9 3) TMRC 2020, Paper A7, S. Grantz et al, "Heat Assisted Magnetic Recording Dependence on Reader for Conventional and Shingled Magnetic Recording" see https://blog.seagate.com/enterprises/energy-assisted-magnetic-recording-will-solve-the-need-for-capacity/ 4) S. Hernandez, et. al, "Geometrical Scaling Limits of Heat-Assisted Magnetic Recording", IEEE Trans. Magn., 2020 (accepted). 5) Nakatani et al, “Read sensor technology for ultrahigh density magnetic recording” MRS Bulletin 2018.

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/f9d9c340bac329ce8b71e72a67697970.jpg)
A range of reader designs were considered to satisfy the down-track and cross-track resolution requirements. The Bottom SV and derivative Recessed AFM design may be sufficient near 4Tbpsi.

![](https://assets.underline.io/uploads/markdown_image/1/image/0ed76210a04117fadf5171ea35378555.jpg)
Reader stack requirements at 4Tbpsi for BAR 2.4 (left) and 3.2 (right). MR requirement varies with free layer damping constant alpha as illustrated. The best fit is shown by design point circle on graph .

 HDD Reader Technology Roadmap to an Areal Density of 4 Tbpsi and Beyond INVITED

The development of high-density magnetic recording media is limited by the superparamagnetism in very small ferromagnetic crystals. Hard magnetic materials with strong perpendicular anisotropy offer stability and high recording density. To overcome the difficulty of writing media with a large coercivity, heat assisted magnetic recording (HAMR) has been developed, rapidly heating the media to the Curie temperature Tc before writing, followed by rapid cooling. Requirements are a suitable Tc, coupled with anisotropic thermal conductivity and hard magnetic properties. Here we introduce Rh2CoSb as a new hard magnet with potential for thin film magnetic recording. A magnetocrystalline anisotropy of 3.6 MJm-3 is combined with a saturation magnetization of μ0Ms = 0.52 T at 2 K (2.2 MJm-3 and 0.44 T at room-temperature). The magnetic hardness parameter of 3.7 at room temperature is the highest observed for any rare-earth free hard magnet. The anisotropy is related to an unquenched orbital moment of 0.42 µB on Co, which is hybridized with neighbouring Rh atoms with a large spin-orbit interaction. Moreover, the pronounced temperature-dependence of the anisotropy that follows from its Tc of 450 K, together with a high thermal conductivity of 20 Wm-1K-1, makes Rh2CoSb a candidate for development for heat assisted writing with a recording density in excess of 10 Tb/in2.

**References:**

Adv. Mater. 2020, 32, 2004331

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/e4f78a1c92343bbb531bc8b98bb97109.png)
Magnetic properties of Rh2CoSb.

![](https://assets.underline.io/uploads/markdown_image/1/image/bd747b83e3690c6f8fce758097e85b48.png)
Comparison of the magnetic hardness parameter κ with other hard magnets.

Downloads

Next from INTERMAG 2021

Double magnet master media for magnetic printing onto energy-assisted magnetic recording media

Similar lecture

Impedance of Planar Structures with Radial Current Distribution

.css-70qvj9{display:-webkit-box;display:-webkit-flex;display:-ms-flexbox;display:flex;-webkit-align-items:center;-webkit-box-align:center;-ms-flex-align:center;align-items:center;}Downloads

Next from INTERMAG 2021

Double magnet master media for magnetic printing onto energy-assisted magnetic recording media

Similar lecture

Impedance of Planar Structures with Radial Current Distribution

Downloads