Supplementary Materials for

Similar documents
Chapter 13: Multiple Regression

DUE: WEDS FEB 21ST 2018

arxiv:cs.cv/ Jun 2000

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

RECEIVED. Negative Transverse Impedance

Uncertainty in measurements of power and energy on power networks

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

Kernel Methods and SVMs Extension

Comparison of Regression Lines

Numerical Transient Heat Conduction Experiment

Temperature. Chapter Heat Engine

PES 1120 Spring 2014, Spendier Lecture 6/Page 1

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Quantum spin system with on-site exchange in a magnetic field

MAGNETISM MAGNETIC DIPOLES

Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields

Relaxation in water /spin ice models

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Indeterminate pin-jointed frames (trusses)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

More metrics on cartesian products

CONDUCTORS AND INSULATORS

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

NP-Completeness : Proofs

Experiment 1 Mass, volume and density

V.C The Niemeijer van Leeuwen Cumulant Approximation

Monte Carlo simulation study on magnetic hysteresis loop of Co nanowires

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Negative Binomial Regression

SUPPLEMENTARY INFORMATION

Appendix B. The Finite Difference Scheme

Constitutive Modelling of Superplastic AA-5083

Pulse Coded Modulation

Oguchi Approximation of a mixed spin-2 and spin-5/2 Blume-Capel Ising ferrimagnetic system

Week 9 Chapter 10 Section 1-5

Computer Simulations of Parallel-to-Series Conversion in Solid State Frame Transfer Image Sensors. J. Bisschop

Electron-Impact Double Ionization of the H 2

Susceptibility and Inverted Hysteresis Loop of Prussian Blue Analogs with Orthorhombic Structure

Electricity and Magnetism - Physics 121 Lecture 10 - Sources of Magnetic Fields (Currents) Y&F Chapter 28, Sec. 1-7

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Enhancing spontaneous emission rates of molecules using nanopatterned multilayer hyperbolic metamaterials

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Supplementary Figures

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

Inductance Calculation for Conductors of Arbitrary Shape

Analysis of Unsteady Aerodynamics of a Car Model with Radiator in Dynamic Pitching Motion using LS-DYNA

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

CHAPTER 14 GENERAL PERTURBATION THEORY

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

Lecture 14: Forces and Stresses

Fundamental loop-current method using virtual voltage sources technique for special cases

Supplemental document

Modeling acoustic transducer surface waves by Transmission Line Matrix method

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

SUPPLEMENTARY INFORMATION

Pressure Measurements Laboratory

NUMERICAL DIFFERENTIATION

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Neutral-Current Neutrino-Nucleus Inelastic Reactions for Core Collapse Supernovae

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

5.62 Physical Chemistry II Spring 2008

Multiple Sound Source Location in 3D Space with a Synchronized Neural System

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

RG treatment of an Ising chain with long range interactions

8. Superfluid to Mott-insulator transition

The Feynman path integral

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Spin waves in two and three dimensional magnetic materials. H. S. Wijesinhe, K. A. I. L. Wijewardena Gamalath

PORE STRUCTURE AND THERMAL CONDUCTIVITY OF BURNT CLAY BRICKS INTRODUCTION

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

Lecture 4. Macrostates and Microstates (Ch. 2 )

Dr. Fritz Wilhelm, Physics 230 E:\Excel files\230 lecture\ch26 capacitance.docx 1 of 13 Last saved: 12/27/2008; 8:40 PM. Homework: See website.

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Mathematical Preparations

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

High resolution entropy stable scheme for shallow water equations

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

Grover s Algorithm + Quantum Zeno Effect + Vaidman

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Semiclassical Model of Electron Dynamics. 8.1 Description of the Semiclassical Model

Quantitative Discrimination of Effective Porosity Using Digital Image Analysis - Implications for Porosity-Permeability Transforms

THE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME

THE CURRENT BALANCE Physics 258/259

One Dimensional Axial Deformations

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE.

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Supporting information

Resource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis

Note on the Electron EDM

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Errors for Linear Systems

Open string operator quantization

Effect of Losses in a Layered Structure Containing DPS and DNG Media

Chapter 8. Potential Energy and Conservation of Energy

Lecture 7: Boltzmann distribution & Thermodynamics of mixing

Transcription:

advances.scencemag.org/cg/content/full/2/7/e1600304/dc1 Supplementary Materals for Interface-drven topologcal Hall effect n SrRuO3-SrIrO3 blayer Jobu Matsuno, Naok Ogawa, Kenj Yasuda, Fumtaka Kagawa, Wataru Koshbae, Naoto Nagaosa, Yoshnor Tokura, Masash Kawasak Publshed 8 July 2016, Sc. Adv. 2, e1600304 (2016) DOI: 10.1126/scadv.1600304 Ths PDF fle ncludes: secton SI. Comparson between the Kerr rotaton angle and the magnetzaton secton SII. Detals of the calculaton of the skyrmon stablty n the multlayer secton SIII. Prelmnary magnetc force mcroscopy mages fg. S1. Kerr rotaton, magnetzaton, and topologcal Hall resstvty. fg. S2. Multlayer model and calculated skyrmon radus. fg. S3. Magnetc force mcroscopy mages for m = 5 measured at 50 K and varous magnetc feld.

secton SI. Comparson between the Kerr rotaton angle and the magnetzaton We performed the magneto-optc Kerr effect (MOKE) measurements wth a laser at 690 nm (= 1.80 ev), correspondng to hgher energy nterband transtons; an essental energy scale of spnorbt nteracton responsble for MOKE can be regarded as perturbaton and therefore the Kerr rotaton angle θk s theoretcally justfed to be proportonal to macroscopc magnetzaton M [P. N. Argyres, Phys. Rev. 97, 334 (1955).] Snce AHE s usually defned to be proportonal to magnetzaton M as a functon of magnetc feld at a fxed temperature, the AHE should be always proportonal to θk, rrespectve of the detaled magnetc structure. We also expermentally confrmed ths; n fg. S1A, we plot the Kerr rotaton angle θk measured by MOKE and the magnetzaton M measured by SQUID magnetometry for the m = 5 sample at 80 K, whch s exactly the same condton wth those n Fg. 2C. We can clearly see the good concdence between θk and M. Subsequently, the extracted THE from MOKE agrees well wth that from SQUID as ndcated n fg. S1B. We can hereby justfy usng θk nstead of M also n our case. The reason why we used the Kerr rotaton angle s that magnetzaton measured for thn flms nevtably ncludes a farly large contrbuton of the substrates, makng t dffcult to measure flm magnetzaton n ultrathn lmt of m = 4. Please note that the M-H measurement s much more dffcult than M-T one shown n Fg. 1A. Whle the latter s measured at 0.05 T, we need 3 T or more for the former; ths ndcates 60 tmes larger background sgnals derved from damagnetc substrates, whch made the flm contrbuton hardly dscernble. Therefore, we used Kerr rotaton angle as M for the sake of consstency through the data wth dfferent thckness (m). fg. S1. Kerr rotaton, magnetzaton, and topologcal Hall resstvty. (A) Magnetc feld dependence of Kerr rotaton angle (θk, orange) and magnetzaton (M, green) of m = 5 at 80 K. (B) Topologcal Hall resstvty deduced from MOKE (orange) and SQUID (green), respectvely.

secton SII. Detals of the calculaton of the skyrmon stablty n the multlayer We numercally examned a sngle skyrmon stablty n the model defned by the followng Hamltonan H J m n l l Dˆ ˆ ˆ, n ', ' y n1, n1, x x n1, n1, yˆ h ll' ' l1 n z l, The model s schematcally shown n fg. S2A, demonstratng that the DM nteracton D exsts only on the frst layer (l = 1). The drecton s fxed by the broken nverson symmetry whle we do not know f D s postve or negatve. The normalzed magnetc moment at the ste on the layer l s denoted as nl,. The unt vectors xˆ and ŷ defne the two-dmensonal square lattce on a layer. The calculaton was performed for a fnte sze system where the layer has 300 300 lattces wth perodc boundary condton. The last term represents the Zeeman energy wth external magnetc feld h perpendcular to the multlayer. Fgure S2B shows the layer-dependent radus of the stablzed sngle skyrmon wth a parameter set of J = 1.0, D = 0.3, h = 0.01. The skyrmon radus at layer l among the total m layers, rm(l), s defned by 2 r ( ) m l 1 z n l, 0 One can see that the layer l dependence of rm(l) s not strong; the skyrmon s nearly cylndrcal. Wthn the whole parameter range nvestgated, the skyrmon s always runnng through all the layers once t s stablzed and hence the monopole structure never appears. In such a stuaton, nl, s are almost parallel to each other between the layers,.e., nl, = n1,. Therefore, the orgnal Hamltonan turns nto the effectve one H eff mj n1, n1, ' ( D / m) yˆ n n xˆ xˆ 1, 1, n1, n1, yˆ ' h n z 1, We can mmedately derve Deff = D / m from ths two-dmensonal Hamltonan when one consders the Hamltonan per one layer,.e., Heff / m.

fg. S2. Multlayer model and calculated skyrmon radus. (A) Schematcs of the multlayer model used n the calculaton of the skyrmon stablty. (B) Estmated skyrmon radus at layer l among the total m (m = 6) layers, rm=6(l).

secton SIII. Prelmnary magnetc force mcroscopy (MFM) mages We measured the MFM mages for the m = 5 sample at 50 K at whch sgnfcant THE sgnal s observed. The result s shown wth the correspondng topologcal Hall resstvty and the Kerr rotaton angle n fg. S3. We observed the magnetzaton reversal process between 0.08 T and 0.50 T, correspondng to the θk( M )-H curve. The strkng feature s the sngle bubble doman of ~30 nm, for example observed at 0.20 T. Consderng the small nature of the doman, we speculate that t does not contan any Bloch lnes and hence can be a skyrmon; the observed mages are compatble wth the skyrmon pcture. Nevertheless, we would lke to pont out that ths cannot be a defnte proof of the skyrmon for the followng reasons: 1) Snce the MFM mages are lackng n the nformaton on n-plane magnetzaton, we cannot expermentally tell whether the bubble has a topologcal charge or not. 2) The spatal resoluton (= 20 nm) s nsuffcent to detect sngle skyrmon features wth the expected sze of ~10 nm. Measurements MFM was conducted under hgh vacuum wth a commercally avalable scannng probe mcroscope (attocube AFM/MFM I). We used SSS-QMFMR cantlever from Nanosensors. The measurement was performed at 50 K wth the cantlever exctaton ampltude of 20 nm and the lft heght of ~30 nm. The resonance frequency was 69.795 khz and the Q-factor was 3.3 10 4 under the measurement condton. We observed magnetzaton reversal of the cantlever at ~0.05 T and therefore the magnetzaton drecton was pontng upward durng the measurement ( 0.08 T). In order to extract ntrnsc magnetc contrast, we have subtracted the background data. For the mages measured at 0.14. 0.20, and 0.26 T, the background s the 0.08-T mage, whle we used the 0.50-T mage for 0.26, 0.32, 0.38, and 0.44 T. The reason why we used the two background data s that there s a magnetc-feld dependent offset n the MFM mages. At 0.26 T, we show the two mages from whch we subtracted each background data. We observed that they are n prncple consstent, despte the small dfference n the appearance.

fg. S3. Magnetc force mcroscopy mages for m = 5 measured at 50 K and varous magnetc feld. Magnetc feld was controlled from 2 T to 0.5 T. Raw mage measured at 0.08 T has been subtracted as a background from upper three mages, whle the background s the 0.50-T mage for lower four mages. At 0.20 T, the sze of the bubble doman enclosed by the dashed square s roughly 30 nm. The magnetc-feld dependence of the correspondng topologcal Hall resstvty and the Kerr rotaton angle s also shown.