Muon relaxation functions Stephen J. Blundell Clarendon Laboratory, - - PowerPoint PPT Presentation

muon relaxation functions
SMART_READER_LITE
LIVE PREVIEW

Muon relaxation functions Stephen J. Blundell Clarendon Laboratory, - - PowerPoint PPT Presentation

Muon relaxation functions Stephen J. Blundell Clarendon Laboratory, Department of Physics, University of Oxford, UK Muon training course - 2018 (Thanks to Francis Pratt for a few of the later slides on muonium-like states.) Muon training course


slide-1
SLIDE 1

Stephen J. Blundell

Clarendon Laboratory, Department of Physics, University of Oxford, UK

Muon training course - 2018

Muon relaxation functions

(Thanks to Francis Pratt for a few of the later slides on muonium-like states.)

Muon training course - 2018

slide-2
SLIDE 2

see later (data analysis sessions) Wimda

Other packages are available

slide-3
SLIDE 3

Lecture plan

  • Static distributions - what is a Kubo-Toyabe?
  • Gaussian or Lorentzian?
  • Dynamic relaxation functions - what happens when

the muons get a bit jumpy?

  • Stretched exponentials - dangerous evil or answer

to all problems?

  • When quantum mechanics shines on the

experiment!

  • Where is your muon?
slide-4
SLIDE 4

Muon spin precession

Response to a Static Field

slide-5
SLIDE 5

Response to a Static Field

slide-6
SLIDE 6

Response to a Static Field

slide-7
SLIDE 7

Distribution of Static Fields

Assume that the Bx, By and Bz components are each distributed according to a Gaussian distribution, e.g. The overall field distribution peaks near 21/2Δ/γµ

slide-8
SLIDE 8

Distribution of Static Fields

slide-9
SLIDE 9

Ryogo Kubo (1920-1995)

Static Kubo-Toyabe Function

slide-10
SLIDE 10

Minimum at t=√3/Δ 1/3 tail

Static Kubo-Toyabe Function

slide-11
SLIDE 11

Kubo-Toyabe in Field

slide-12
SLIDE 12

Kubo-Toyabe in Longitudinal Field

slide-13
SLIDE 13

2

LF Kubo-Toyabe Function

slide-14
SLIDE 14

Gaussian or Lorentzian Field Distribution?

(Dense spins)

slide-15
SLIDE 15

Gaussian or Lorentzian Field Distribution?

Dilute spins:

slide-16
SLIDE 16

Lorentzian Kubo-Toyabe (LKT)

slide-17
SLIDE 17

The ‘Dilute’ Spin Condition

e.g. metallic random alloy spin glasses e.g. complex molecular magnets

A broad range of couplings from the muon to the nearest spin is the key here

slide-18
SLIDE 18

All these functions available in Wimda

F.L. Pratt, Physica B 289-290, 710 (2000) http://shadow.nd.rl.ac.uk/wimda/

Other packages are available

slide-19
SLIDE 19

One can interpolate between statics and dynamics using a dynamical Kubo-Toyabe function

Introduce Dynamics

slide-20
SLIDE 20

Introduce Dynamics

slide-21
SLIDE 21

Introduce Dynamics

slide-22
SLIDE 22

Introduce Dynamics

slide-23
SLIDE 23

Introduce Dynamics

slide-24
SLIDE 24

Dynamical Kubo-Toyabe (DKT)

slide-25
SLIDE 25

Dynamical Kubo-Toyabe (DKT)

Since the Gz integral depends only on Gz at earlier times and the known static function gz , it can be built up sequentially in time

slide-26
SLIDE 26

Slow Hopping

slide-27
SLIDE 27

Fast Hopping

slide-28
SLIDE 28

Static Dynamic Effect of Longitudinal Field

slide-29
SLIDE 29

Zero-field Longitudinal-field Effect of Dynamics

slide-30
SLIDE 30

The Keren Function

Perturbation expansion for Pz(t) gives an analytical result valid for ν>Δ Amit Keren ZF limit: (LF Abragam) Fast ν limit: (ν>ωL)

PRB 50,10039 (94)

slide-31
SLIDE 31

Stretched Exponential Functions

slide-32
SLIDE 32

Stretched Exponential Functions Gaussian β=2 Lorentzian β=1 Stretched β<1

Stretched exponentials generally arise from: 1) Distribution of relaxation times 2) Distribution of couplings

lineshape parameter

slide-33
SLIDE 33

Distribution of Relaxation Time

slide-34
SLIDE 34

Distribution of Relaxation Times: Diffusion

Risch-Kehr function

PRB 46, 5246 (1992)

1D:

Klaus Kehr (1934-2000)

slide-35
SLIDE 35

B-dependent Relaxation and Spectral Density Φ(t) =

Correlation function for field fluctuations: Fourier transform of Φ(t) gives the spectral density S(ω) S(ω) = ν / (ν2 + ω2) λ is proportional to S(ωL) Complex relaxation processes such as those based on diffusion typically give power laws for Φ(t) and S(ω)

slide-36
SLIDE 36

B-dependent Relaxation and Spectral Density

Fe19

Blundell et al (2003)

Polyaniline

Stretched relaxation is due to the coupling distribution Stretched relaxation is due to the 1D diffusion process

Single molecule magnet

slide-37
SLIDE 37

Muons that stop closer to magnetic ions “see” a wider local field distribution (which extends to higher fields) than muons which stop at a greater distance Y.J. Uemura et al, PRB 31, 546 (1985)

Distribution of Couplings Spin Glasses

slide-38
SLIDE 38

The correct relaxation function must therefore be an average over distribution widths Δ. This leads to a root-exponential relaxation function: G(t) = G(0) exp(-(λt)1/2) where the relaxation rate λ is inversely proportional to the fluctuation rate ν.

Distribution of Couplings

slide-39
SLIDE 39

Non-magnetic host

Magnetic Non-magnetic

Distribution of Couplings

slide-40
SLIDE 40

Spin glass

Muon stops close to magnetic ion

Magnetic Non-magnetic

Distribution of Couplings

slide-41
SLIDE 41

Spin glass

Muon stops well away from magnetic ion

Magnetic Non-magnetic

Distribution of Couplings

slide-42
SLIDE 42

Muons that stop closer to magnetic ions “see” a wider local field distribution (which extends to higher fields) than muons which stop at a greater distance

Y.J. Uemura et al, PRB 31, 546 (1985)

Distribution of Couplings Spin Glasses

slide-43
SLIDE 43

Distribution of Couplings

slide-44
SLIDE 44

……but is the dogma correct?

Distribution of Couplings

slide-45
SLIDE 45

Monte-Carlo calculation

  • f distribution of Δ
  • S. J. Blundell, T. Lancaster, F. L. Pratt, C. A. Steer, M. L. Brooks and J. F. Letard, J. Phys. IV France 114, 601 (2004)

Distribution of Couplings

slide-46
SLIDE 46
  • S. J. Blundell, T. Lancaster, F. L. Pratt, C. A. Steer, M. L. Brooks and J. F. Letard, J. Phys. IV France 114, 601 (2004)

Distribution of Couplings

slide-47
SLIDE 47

1. 2. 3.

slide-48
SLIDE 48

Hints of Quantum Coherence

slide-49
SLIDE 49

F- µ+

The F-µ-F State

slide-50
SLIDE 50

The F-µ-F State

slide-51
SLIDE 51

Pz(t) = 1 6 3 + cos √ 3ωt + ✓ 1 − 1 √ 3 ◆ cos " 3 − √ 3 2 ωt # + ✓ 1 + 1 √ 3 ◆ cos " 3 + √ 3 2 ωt #!

coherent

  • scillations arising

from the magnetic dipolar interaction

F- F-

µ+

The F-µ-F State

slide-52
SLIDE 52

Fluorine: small, high nuclear moment abundant species very sensitive to r1/r2 and α

entanglement

The F-µ-F State

slide-53
SLIDE 53

state found in many ionic fluorides, and also teflon (PTFE)

The F-µ-F State

slide-54
SLIDE 54

Cs2AgF4

  • T. Lancaster, S. J. Blundell, et al.

Tc=13.95(3) K

ωd = 2π x 0.211(1) MHz F-µ separation 1.19(1) Å.

  • Phys. Rev. B 75, R220408

(2007)

The F-µ-F State

slide-55
SLIDE 55

CuF2(H2O)2(pyz)

  • Phys. Rev. Lett. 99, 267601 (2007)
  • 1. Interaction

with a single fluorine ion

The F-µ-F State

slide-56
SLIDE 56

[CuNO3(pyz)2]PF6

  • 2. Crooked FµF

bond close to a PF6 ion

  • Phys. Rev. Lett. 99, 267601 (2007)

The F-µ-F State

slide-57
SLIDE 57

[Cu(HF2)(pyz)2]X

  • 3. Interaction

with a HF2

  • ion
  • Phys. Rev. Lett. 99, 267601 (2007)

The F-µ-F State

slide-58
SLIDE 58

Analysing Asymmetry : Magnets

Polycrystalline samples

b = B / B0

LF decoupling below Tc

  • F. L. Pratt
slide-59
SLIDE 59

Analysing Asymmetry : Muonium-like States

LF decoupling

  • r ‘repolarisation’

Hyperfine tensor (A,D1,D2)

F . L. Pratt, Phil. Mag. Lett. 75, 371 (1997)

Nuclear couplings AN

  • Z. Phys. B 22, 109 (1975)

D1 D2 A (500,20,1) MHz A AN

slide-60
SLIDE 60

Analysing Asymmetry : Muonium-like States

Avoided level-crossing resonances:

Δ1 muon flip Δ0 muon-nuclear flip-flop

  • F. L. Pratt
slide-61
SLIDE 61

All these functions available in Wimda

F.L. Pratt, Physica B 289-290, 710 (2000) http://shadow.nd.rl.ac.uk/wimda/

Other packages are available

slide-62
SLIDE 62

DFT+µ = (density functional theory +µ )

  • numerically solve (lattice) structures
  • determine muon site
  • quantify perturbations

DFT+µ began with two papers (Oxford + Parma groups) studying fluorides: J.S. Mö̈ller et al., Phys. Rev. B 87, 121108(R) (2013).

  • F. Bernadini, et al., Phys. Rev. B 87, 115148 (2013).

This work has been extended to many other systems, see e.g. S.J. Blundell et al, Phys. Rev. B 88, 064423 (2013). J.S. Mö̈ller et al., Phys. Scr. 88, 068510 (2013).

  • F. Xiao et al., Phys. Rev. B 91, 144417 (2015).
  • P. Bonfà et al., J. Phys. Chem. C 119, 4278 (2015).
  • F. Lang et al., Phys. Rev. B 94, 020407(R) (2016).
  • P. Bonfà et al., J. Phys. Soc. Jpn. 85, 091014 (2016)

DFT+µ

slide-63
SLIDE 63

DFT+µ = (density functional theory +µ )

  • numerically solve (lattice) structures
  • determine muon site
  • quantify perturbations

DFT+µ can not only assess the muon site, but also any muon- induced distortion. A worst-case scenario is where magnetism arises from a non-Kramers ground state. This leads to our study

  • f quantum spin ice.
  • F. R. Foronda et al., Phys. Rev. Lett. 114, 017602 (2015).

Challenges for pyrochlores:

  • 88 atoms per unit cell
  • 4f valence electrons
  • ~102-104 cpu hours per calculation

Results:

  • typical O-H like bond with length 1 Å
  • 4f electrons influence negligible
  • r4f ≈ 2×r5s ≈ 5×(Pr-µ distance)

DFT+µ

slide-64
SLIDE 64

…it’s now time for some relaxation! After all those relaxation functions…..