Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Park Park Park Park Park Park Park Park
Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Park Park Park Park Park Park Park Park
Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Kyungwha Park Park Park Park Park Park Park Park
Kyungwha Park
Department of Physics, Virginia Tech
Supported by Jeffress Memorial Trust Fund NCSA clusters, Virginia Tech System X, VT clusters
Outline
What are single-molecule magnets (SMMs)? Motivation Methodology: density-functional theory (DFT) Review of properties of isolated SMM Mn12 SMM Mn12 deposited on a gold surface
- Method and model - Changes in electronic structure (orbital broadening, charge transfer) - Changes in magnetic properties
Mn12
1.5nm S=10
Mn4
Mn3+ Mn4+ O
Cl N C
S=9/2
Examples: Mn12 , Mn4 , Fe4, Fe8, Co4, Ni4, cyanide-bridged molecules, Mn84 D. Hendrickson J. Am. Chem. Soc. 114, 2455 (1992) torus, etc.
S=6
Fe2Ni2
Quantum mechanical
Energy = Dms2 ; D > 0 for B=0
Spin projection ms
Barrier
Energy
S=10
"up"
In ground state
Ms
down up
tobs=500 s
Steps: tunneling on Plateaus: tunneling off Classical ferromagnets at 300 K HL (T) Off resonance down up
B. Barbara et al., J. Phys. Soc. Jpn. (2000); J. Friedman et al., PRL (1996); L. Thomas et al., Nature (1996)
= 0 exp(U/ kBT),
C. Paulsen et al., J. Mag. Mag. Mat. (1995) A. Caneschi et al., Nature (1993)
Magnetic measurement:
Properties of Mn12 monolayers differ from those of bulk Mn12
Naitabdi et al., Adv. Mater. 17. 1612 (2005) Salman et al., Nano Lett. (2007)
STM image
del Pennino et al., Surf. Sci. (2006); Voss et al., PRB (2007)
DFT
Kohn & Pople (1998, Nobel prize in chemistry)
N non-interacting single-electron problems: Hi=ii where i=i (x,y,z) r 2 H= / 2 + Veff (r ) where Veff (r) is an effective potential Cost: NM3 for M x M x M grids (e.g.: 10 x103 for N=10, M=10 )
r r r 2 (r ) = Ck f k exp[ k (r Rk ) ]
Solve N nonnon-interacting one-electron equations under effective potential self-consistently: Veff (r)= Vnucl+VH[(r)]+Vxc[(r)]
NO
Energy is converged?
r [ / 2 + Veff (r )]i = i i
2
YES
Calculate forces, (r), geometry
NO YES Geometry
Force is small? is optimized!
Mn12: [Mn12O12(O2CCH3)16(H2O)4]
T.Lis Acta Cryst. B36, 2042 (1980)
1.5nm
y x
Mn4+ Mn3+
Experimental data
Single molecule: S4 symmetry 4 Mn4+ (3d3 ,S=3/2) ions in cube 8 Mn3+ (3d4 ,S=2) ions in outer crown Easy axis: z axis
M.R. Pederson & K.A. Jackson, PRB (1990); ibid, PRB (1991); K.A. Jackson & M.R. Pederson, PRB (1990); D. Porezag & M.R. Pederson, PRB (1996)
How does a metal surface affect electronic structure and magnetic anisotropy of Mn12?
gold surface
Experiment
Monolayer of Mn12 on gold surface
L. Zobbi et al., Chem. Comm. 1640 (2005).
~20 nm
S. Barraza-Lopez et al., Phys. Rev. B (2007).
Theory
Consider slab calculation: Mn12-S2gold surface
short link
- All-electron wavefunctions are available - To take into account spin-orbit coupling - Accuracy is improved in magnetic materials
Au: s1 d10 Mn: 3p6 4s2 3d5 O: s2 p4 C: s2 p2 S: s2 p4
Au Fermi level
LUMO HOMO
Different scale
Orbital broadening
Mn12 molecular orbitals broaden upon adsorption Orbital broadening is much less than charging energy Weak coupling
Isolated
Isolated
AB C
Mn atoms in region C
Mn(I): inner, -2.588 B Mn(II): outer, 3.524 B Mn(III): outer, 3.543 B monoclinic unit cell 17.7 x 17.7 x 34.0 3
2.0
-e
v
Anisotropic
Mn
r r r r r B = v E = v ( )
: Coulomb potential due to atomic nuclei and electrons
Although total magnetic moment differs for Mn12 and Geo.2 (S-terminated Mn12), the MAB of ordinary Mn12 is the same as that of Geo. 2 Magnetic anisotropy barrier for Geo. 3 gets reduced by about 9%
Geometry 3
Mn12 66.7
Geo 2 66.9
Geo 3 60.7
S H
The spatial charge and magnetic moment distributes in a similar way for whole structure and Geo. 3. So [magneitc ansiotropy barrier (MAB) of Geo 3] = [MAB of whole structure] The magnetic anisotropy barrier for the whole structure is reduced by 9% of that for the isolated Mn12.
2.0
2.0 B
Mn12 66.7
V. I. Anisimov et al., PRB 44, 44 943 (1991) V. I. Anisimov et al., J. Phys.: Condens. Matter 9, 767 (1997)
U=6*
PBE GGA
Mn12 HOMO Mn12 LUMO gap Geo 2 HOMO Geo 2 LUMO gap
1.35
0.24
-6.37 -6.66 -5.37 -6.07 -6.05 -5.51 -5.22 -4.84 0.32 0.15
Summary
Modeled a monolayer of Mn12 adsorbed on a Au surface via a thiol group using DFT Electronic structure: broadening of molecular orbitals was small, weak coupling between Mn12 and a Au surface even with the short link Charge transfer from Au to Mn12 Total moment for whole structure back to 20 B Magnetic anisotropy barrier for whole structure is reduced by 9% compared to that for isolated Mn12 molecule Effect of on-site U on orbitals for isolated molecules
Barraza-Lopez, Avery and Park: PRB 76 224413 (2007), JAP 103 07B907 (2008)
Check convergence with # of k points and energy cutoff for plane waves With calculated lattice constant, we vary # of Au layers and # of vertical vacuum layers Au surface energy saturates with 8 Au monolayers and 7 vacuum layers To cover Mn12, at least 36 surface gold atoms per monolayer are needed Use six gold monolayers (36 x 6 = 216 gold atoms) Relax gold slab (w/o molecule) until max forces are less than 0.01 eV/
vacuum layer
3d
3d
3d
Broadening of orbitals is much less than Coulomb charging energy: Weak coupling
2p
HOMO
+2.0