Magnetic avalanches in granular ferromagnets: Thermal activated collective behavior Gia-Wei Chern1 1Department of Physics, University of Virginia, Charlottesville, VA 22904, USA We present a numerical study on the thermal activated avalanche dynamics in granular materials composed of ferromagnetic clusters embedded in a non-magnetic matrix. A microscopic dynami- cal simulation based on the reaction-diffusion process is developed to modeling the magnetization process of such systems. The large-scale simulations presented here explicitly demonstrate inter- granular collective behavior induced by thermal activation of spin tunneling. In particular, we observe an intriguing criticality controlled by the rate of energy dissipation. We show that ther- mal activated avalanches can be understood in the framework of continuum percolation and the emergent dissipation induced criticality is in the universality class of three-dimensional percolation transition. Implications of these results to the phase-separated states of colossal magnetoresistance materials and other artificial granular magnetic systems are also discussed. 1. Introduction tion, the spin dynamics is dominated by single-site spin tunneling assisted by thermal activation. Consequently the avalanche behavior is mostly controlled by the diffu- Avalanche phenomena are ubiquitous in Nature [1]. sion and dissipation of thermal energy. Fundamentally, When systems are driven by external forces, their re- magnetic deflagration can be understood in the frame- sponses are often in the form of intermittent bursts of work of reaction-diffusion process [10{12] which also de- activities or avalanches. One well-studied example is the scribes avalanche phenomena, ranging from phase tran- Barkhausen noise in magnetic system which corresponds sition dynamics to pattern formation and propagation of to avalanches of domain movements. As avalanche be- epidemic waves. haviors often extend over a huge range of sizes, It has A key component of thermal activated avalanches is been suggested that avalanche phenomena can be clas- the self-sustainability of excess thermal energies that can sified into universality classes that are independent of be used to overcome the energy barriers. When a system the microscopic details. A canonical example is the ran- is trapped in a meta-stable state, a small perturbation dom field Ising model (RFIM) which describes a system or a quick variation of the external conditions can induce of interacting Ising spins subject to site-dependent ran- relaxation of a portion of the system through thermal dom fields [2{4]. When this Ising system is driven by a activated tunneling. The energy released from this initial magnetic field, it responds in the form of sporadic large- relaxation can be reabsorbed by the rest of the system scale spin reconfigurations. These spin avalanches de- and used to overcome the barrier again. An avalanche pend strongly on the level of disorder: while the magne- is then ignited when the relaxation process becomes self- tization dynamics is controlled by large system-spanning sustained. In this scenario, the rate of energy dissipation events at low disorder, small avalanche clusters dominate naturally plays an important factor as thermal activation the magnetic behavior when disorder is increased. Re- becomes negligible if a significant amount of the released markably, the system exhibits a power-law distribution energy is removed from the system. of avalanche sizes at a critical disorder. This nonequilib- Thermal coupling has also been shown to affect the rium criticality induced by quenched disorder is charac- magnetization dynamics of an important class of materi- terized by a unique set of critical exponents as in equilib- als known to exhibit a colossal magentoresistence (CMR), rium phase transitions [3, 4]. The RFIM has since been including perovskite manganites and cobaltites [13{16]. applied to a number of other systems beyond spins. Temperature-dependent ultrasharp magnetization steps Conventional wisdom also suggests that the avalanche have been observed in several doped Pr0:5Ca0:5MnO3 behaviors are unaffected by thermal fluctuations and are manganites. Moreover, magnetic deflagration induced arXiv:1611.01098v1 [cond-mat.dis-nn] 3 Nov 2016 governed by a deterministic dynamics that depends on by surface acoustic waves has been reported in La-based the static, quenched disorder in the system [1, 2]. In this manganites [17]. Recent experiments have also observed regard, avalanche is an intrinsic many-body phenomenon unusual sweep-rate dependence of the avalanche behav- in which spin-spin interactions play an indispensable role. ior in CMR materials [18, 19]. These observations all In many systems, thermal effects are indeed negligible as point to the important role of thermal activations in the the energy barriers of spin-flips are often too large for avalanche dynamics of these perovskite oxides. Moti- thermal activation to be effective. On the other hand, vated by these experiments, the goal of this paper is to recent experimental and theoretical studies on magnetic present a theoretical investigation of thermal activated avalanches in molecular magnets have shown that ther- collective behaviors in CMR-related compounds. mal activation is the dominant factor [5{9]. In partic- A crucial realization in the study of CMR materials is ular, in the process known as the magnetic deflagra- the key role of heterogeneity in understanding the CMR 2 phenomena [20]. For example, doped cobaltites such as (a) La1 xSrxCoO3 have been shown to exhibit a particu- larly− clear form of this inhomogeneity by electron mi- croscopy [21] and nuclear magnetic resonance [22]. Rele- vant to our study here is the small doping regime, where the ground state can be viewed as consisting of nano- scale ferromagnetic metallic clusters embedded in a non- magnetic insulating matrix. Such phase-separated state thus leads to a natural formation of a granular ferro- magnet, analogous to the artificial granular systems that exhibits a giant magneto-resistance (GMR) behavior [23{ 26]. These GMR meta-materials are composed of ferro- magentic clusters deposited in a non-magnetic metallic or insulating matrix. Indeed, lightly doped perovskite cobaltites exhibit a hysteretic magnetoresistance behav- -1 ior with temperature and field dependence similar to the 10 artificial GMR materials [27]. -3 s/s In the present work we perform microscopic dynamical 10 e− m simulations to investigate the intrinsic effects of thermal activation in the avalanche dynamics of granular ferro- 10-5 magnets. The model of magnetic clusters used in this P (s) work is specifically designed such that inter-granular co- 10-7 operative phenomena are absent in the zero-temperature athermal dynamics. We develop a numerical scheme 10-9 based on a microscopic implementation of the reaction- (b) diffusion process modified for the granular systems. The 200 400 600 large-scale simulations on lattice of N ∼ 108 sites clearly s demonstrate thermal activated collective phenomena in FIG. 1: (a) Schematic diagram of magnetic clusters embed- the magnetization dynamics of such materials. In ad- ded in a non-magnetic matrix. The gray cells denote mag- dition, we observe an intriguing non-monotonic depen- netic ions, while white cells represent non-magnetic sites. dence of the avalanche behaviors on the dissipation co- The arrows correspond to the Ising variable σ = ±1 of the efficient, indicating the existence of emergent criticality magnetic ions. (b) Distribution of magnetic clusters gener- controlled by the rate of energy dissipations in granular ated through a site-percolation process with probability p = ferromagnets. A continuum percolation model is then de- 0.05. The distribution at large s follows a exponential decay veloped to understand the numerical results. In partic- P (s) ∼ exp(−s=sm) with sm ≈ 36. ular, the emergent out-of-equilibrium critical avalanche corresponds to the percolation transition in this picture. shell. This LS state is energetically close to two other spin-state configurations [28, 29]: an intermediate-spin 2. Model of magnetic clusters (IS) state with partially filled eg and t2g orbitals and spin degree of freedom S = 1, and a high-spin (HS) state The granular ferromagnets considered in this work with active t2g orbitals and S = 2. Spin-state crossover are composed of magnetic clusters embedded in a non- can be induced either by enhanced temperatures or the magnetic matrix. In order to perform microscopic simu- application of a magnetic field. lations, we define the granular model on a simple cubic lattice. A schematic diagram of the granular magnet is A large magnetoresistance observed in hole doped shown in Fig. 1(a). To be more specific, here we will use cobaltites La1 xSrxCoO3 has reinvigorated the research − the case of La1 xSrxCoO3 as the example [21], although of this perovskite system [31]. The Co ions in the parent − our analysis and conclusion to be discussed in the next compound LaCoO3 occupy the B-sites of the perovskite section does not depend on details of how the clusters structure and form a cubic lattice. The La atoms sit at are generated. The parent compound LaCoO3 itself has the A-sites that correspond to the body centers of indi- attracted considerable attention because of an intriguing vidual Co cubes. Substitution of Sr2+ for La3+ intro- thermal induced insulator-metal crossover that is accom- duces holes into the CoO3 array and locally distorts the panied by a spin-state transition of Co ions [28]. At low lattice such that the neighboring eight Co ions are ele- temperatures and small fields, the trivalent Co ions are vated to the IS state with magnetic moment S = 1 [30]. in the non-magnetic low-spin (LS) state with a filled t2g As the holes are confined to the IS clusters, the itinerant 3 1 electrons interact with the localized spins of IS ions via (a) (b) the Hund's coupling. The metallic IS clusters can thus 20 be thought of as nanoscale ferromagnets stabilized by the 0.5 double exchange mechanism [32, 33]. 15 In our simulations, the metallic clusters are generated 0 through a site percolation process [34] as follows.