Models of advanced recording systems: A multi-timescale micromagnetic code for granular thin film magnetic recording systems

Micromagnetic modelling provides the ability to simulate large magnetic systems reliably without the computational cost limitation imposed by atomistic modelling. Through micromagnetic modelling it is possible to simulate systems consisting of thousands of grains over a time range of nanoseconds to years, depending upon the solver used. Here we present the creation and release of an open-source multi-timescale micromagnetic code combining three key solvers: Landau-Lifshitz-Gilbert; Landau-Lifshitz-Bloch; Kinetic Monte Carlo. This code, called MARS (Models of Advanced Recording Systems), is capable of accurately simulating the magnetisation dynamics in large and structurally complex single- and multi-layered granular systems as is shown by comparison to established atomistic simulation results. The short timescale simulations are achieved for systems far from and close to the Curie point via the implemented Landau-Lifshitz-Gilbert and Landau-Lifshitz-Bloch solvers respectively. This enables read/write simulations for general perpendicular magnetic recording and also state of the art heat assisted magnetic recording (HAMR). The long timescale behaviour is simulated via the Kinetic Monte Carlo solver, enabling investigations into signal-to-noise ratio and data longevity. The combination of these solvers opens up the possibility of multi-timescale simulations within a single software package. For example the entire HAMR process from initial data writing and data read back to long term data storage is possible via a single simulation using MARS. The use of atomistic parameterisation for the material input of MARS enables highly accurate material descriptions which provide a bridge between atomistic simulation and real world experimentation. Thus MARS is capable of performing simulations for all aspects of recording media research and development. This ranges from material characterisation and optimisation to system design and implementation. The object orientated nature of MARS is structured to facilitate quick and simple development and easy implementation of user defined custom simulation types which can utilise either timescale or a combination of both timescales. Program summary: Program title: MARS CPC Library link to program files: https://doi.org/10.17632/8mx7cndcdx.1 Developer's repository link: https://bitbucket.org/EwanRannala/mars/ Code Ocean capsule: https://codeocean.com/capsule/2549929 Licensing provisions: MIT Programming language: C++ Supplementary material: MARS testing methodology (PDF), HAMR simulation example video. Nature of problem: A combined model that enables the complete modelling of magnetic recording processes at elevated temperatures covering all time scales from writing (nanoseconds) up to long term data storage (years). The model must also accurately describe the granular nature of the recording media as grain sizes are reduced to a few nanometres. Solution method: Short timescale behaviours are captured via the Landau-Lifshitz-Gilbert and Landau-Lifshitz-Bloch solvers for low and high temperature systems respectively. The long time scale behaviours are captured via a kinetic Monte Carlo solver. To enable complex models which account for mixed timescale behaviours the solvers are implemented as a single class structure which allows for dynamic solver selection. The granular structure is generated via a Laguerre-Voronoi tessellation with a custom implemented packing algorithm to produce highly realistic grain size distributions. Complex thermal dependencies of materials can be incorporated via atomistic parameterisation forming a multi-timescale model of the material.